Search USGSSearch

SEARCH · Search USGS

Results for “Journal of Geophysical Research: Solid Earth”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

825 records · Page 46Linked to original sources

Seismic velocity anisotropy and heterogeneity beneath the Mantle Electromagnetic and Tomography Experiment (MELT) region of the East Pacific Rise from analysis of P and S body waves

We use teleseismic P and S delay times and shear wave splitting measurements to constrain isotropic and anisotropic heterogeneity in the mantle beneath the southern East Pacific Rise (SEPR). The data comprise 462 P and S delay times and 18 shear wave splitting observations recorded during the Mantle Electromagnetic and Tomography (MELT) Experiment. We estimate the mantle melt content (F) and temperature (T) variation from the isotropic velocity variation. Our results indicate that the maximum variation in F beneath our array is between zero and ???1.2%, and maximum variation in T is between zero and ???100 K. We favor an explanation having partial contributions from both T and F. We approximate the seismic anisotropy of the upper mantle with hexagonal symmetry, consistent with the assumption of two dimensionality of mantle flow. Our new tomographic technique uses a nonlinear inversion of P and slow S polarization delay times to simultaneously solve for coupled VP and VS heterogeneity throughout the model and for the magnitude of anisotropy within discrete domains. The domain dimensions and the dip of the anisotropy are fixed for each inversion but are varied in a grid search, obtaining the misfit of the models to the body wave delay data and to split times of vertically propagating S waves. The data misfit and the isotropic heterogeneity are sensitive to domain dimensions and dip of anisotropy. In a region centered beneath the SEPR the best average dip of the hexagonal symmetry axis is horizontal or dipping shallowly (<30??) west. Given the resolution of our data, a subaxial region characterized by vertically aligned symmetry axes may exist but is limited to be <80 km deep. We infer that the mantle flow beneath the SEPR is consistent with shallow asthenospheric return flow from the direction of the South Pacific superswell.

Journal of Geophysical Research B: Solid Earth

Slip distribution of the 1952 Tokachi-Oki earthquake (M 8.1) along the Kuril Trench deduced from tsunami waveform inversion

We inverted 13 tsunami waveforms recorded in Japan to estimate the slip distribution of the 1952 Tokachi-Oki earthquake (M 8.1), which occurred southeast off Hokkaido along the southern Kuril subduction zone. The previously estimated source area determined from tsunami travel times [Hatori, 1973] did not coincide with the observed aftershock distribution. Our results show that a large amount of slip occurred in the aftershock area east of Hatori's tsunami source area, suggesting that a portion of the interplate thrust near the trench was ruptured by the main shock. We also found more than 5 m of slip along the deeper part of the seismogenic interface, just below the central part of Hatori's tsunami source area. This region, which also has the largest stress drop during the main shock, had few aftershocks. Large tsunami heights on the eastern Hokkaido coast are better explained by the heterogeneous slip model than previous uniform-slip fault models. The total seismic moment is estimated to be 1.87 ?? 1021 N m, giving a moment magnitude of Mw = 8.1. The revised tsunami source area is estimated to be 25.2 ?? 103 km2, ???3 times larger than the previous tsunami source area. Out of four large earthquakes with M ??? 7 that subsequently occurred in and around the rupture area of the 1952 event, three were at the edges of regions with relatively small amount of slip. We also found that a subducted seamount near the edge of the rupture area possibly impeded slip along the plate interface.

Journal of Geophysical Research B: Solid Earth

Shallow soil CO2 flow along the San Andreas and Calaveras Faults, California

We evaluate a comprehensive soil CO2 survey along the San Andreas fault (SAF) in Parkfield, and the Calaveras fault (CF) in Hollister, California, in the context of spatial and temporal variability, origin, and transport of CO2 in fractured terrain. CO2 efflux was measured within grids with portable instrumentation and continously with meteorological parameters at a fixed station, in both faulted and unfaulted areas. Spatial and temporal variability of surface CO2 effluxes was observed to be higher at faulted SAF and CF sites, relative to comparable background areas. However, ??13C (-23.3 to - 16.4???) and ??14C (75.5 to 94.4???) values of soil CO2 in both faulted and unfaulted areas are indicative of biogenic CO2, even though CO2 effluxes in faulted areas reached values as high as 428 g m-2 d-1. Profiles of soil CO2 concentration as a function of depth were measured at multiple sites within SAF and CF grids and repeatedly at two locations at the SAF grid. Many of these profiles suggest a surprisingly high component of advective CO2 flow. Spectral and correlation analysis of SAF CO2 efflux and meteorological parameter time series indicates that effects of wind speed variations on atmospheric air flow though fractures modulate surface efflux of biogenic CO2. The resulting areal patterns in CO2 effluxes could be erroneously attributed to a deep gas source in the absence of isotopic data, a problem that must be addressed in fault zone soil gas studies.

Journal of Geophysical Research B: Solid Earth

Late Pleistocene to Holocene slip rates for the Gurvan Bulag thrust fault (Gobi-Altay, Mongolia) estimated with 10Be dates

We surveyed morphotectonic markers along the central part of the Gurvan Bulag thrust, a fault that ruptured with the Bogd fault during the Gobi-Altay earthquake (1957, M 8.3), to document climatic and tectonic processes along the fault for the late Pleistocene- Holocene period. The markers were dated using 10Be produced in situ. Two major periods of alluviation ended at 131 ?? 20 and 16 ?? 4.8 ka. These appear to be contemporaneous with global climatic changes at the terminations of marine isotope stages (MIS) 6 and 2. The vertical slip rates, determined from offset measurements and surfaces ages, are 0.14 ?? 0.03 mm/yr over the late Pleistocene-Holocene and between 0.44 ?? 0.11 and 1.05 ?? 0.25 mm/yr since the end of the late Pleistocene. The higher of these slip rates for the last ???16 kyr is consistent with paleoseismic investigations along the fault [Prentice et al., 2002], and suggests that, at the end of late Pleistocene, the fault evolved from quiescence to having recurrence intervals of 4.0 ?? 1.2 kyr for surface ruptures with ???4 m vertical offset (similar to that of 1957). The inferred recurrence interval is comparable to that of the Bogd fault (3.7 ?? 1.3 kyr) suggesting that the two faults may have ruptured together also earlier during the last ???16 kyr.

Journal of Geophysical Research B: Solid Earth

Tsunamis generated by subaerial mass flows

Tsunamis generated in lakes and reservoirs by subaerial mass flows pose distinctive problems for hazards assessment because the domain of interest is commonly the "near field," beyond the zone of complex splashing but close enough to the source that wave propagation effects are not predominant. Scaling analysis of the equations governing water wave propagation shows that near-field wave amplitude and wavelength should depend on certain measures of mass flow dynamics and volume. The scaling analysis motivates a successful collapse (in dimensionless space) of data from two distinct sets of experiments with solid block "wave makers." To first order, wave amplitude/water depth is a simple function of the ratio of dimensionless wave maker travel time to dimensionless wave maker volume per unit width. Wave amplitude data from previous laboratory investigations with both rigid and deformable wave makers follow the same trend in dimensionless parameter space as our own data. The characteristic wavelength/water depth for all our experiments is simply proportional to dimensionless wave maker travel time, which is itself given approximately by a simple function of wave maker length/water depth. Wave maker shape and rigidity do not otherwise influence wave features. Application of the amplitude scaling relation to several historical events yields "predicted" near-field wave amplitudes in reasonable agreement with measurements and observations. Together, the scaling relations for near-field amplitude, wavelength, and submerged travel time provide key inputs necessary for computational wave propagation and hazards assessment.

Journal of Geophysical Research B: Solid Earth

Topographically driven groundwater flow and the San Andreas heat flow paradox revisited

Evidence for a weak San Andreas Fault includes (1) borehole heat flow measurements that show no evidence for a frictionally generated heat flow anomaly and (2) the inferred orientation of ??1 nearly perpendicular to the fault trace. Interpretations of the stress orientation data remain controversial, at least in close proximity to the fault, leading some researchers to hypothesize that the San Andreas Fault is, in fact, strong and that its thermal signature may be removed or redistributed by topographically driven groundwater flow in areas of rugged topography, such as typify the San Andreas Fault system. To evaluate this scenario, we use a steady state, two-dimensional model of coupled heat and fluid flow within cross sections oriented perpendicular to the fault and to the primary regional topography. Our results show that existing heat flow data near Parkfield, California, do not readily discriminate between the expected thermal signature of a strong fault and that of a weak fault. In contrast, for a wide range of groundwater flow scenarios in the Mojave Desert, models that include frictional heat generation along a strong fault are inconsistent with existing heat flow data, suggesting that the San Andreas Fault at this location is indeed weak. In both areas, comparison of modeling results and heat flow data suggest that advective redistribution of heat is minimal. The robust results for the Mojave region demonstrate that topographically driven groundwater flow, at least in two dimensions, is inadequate to obscure the frictionally generated heat flow anomaly from a strong fault. However, our results do not preclude the possibility of transient advective heat transport associated with earthquakes.

Journal of Geophysical Research B: Solid Earth

Evidence for a Battle Mountain-Eureka crustal fault zone, north-central Nevada, and its relation to Neoproterozoic-Early Paleozoic continental breakup

Combined evidence from gravity, radiogenic isotope, and magnetotelluric (MT) data indicates a crustal fault zone that coincides with the northwest-trending Battle Mountain-Eureka (BME) mineral trend in north-central Nevada, USA. The BME crustal fault zone likely originated during Neoproterozoic-Early Paleozoic rifting of the continent and had a large influence on subsequent tectonic events, such as emplacement of allochthons and episodic deformation, magmatism, and mineralization throughout the Phanerozoic. MT models show the fault zone is about 10 km wide, 130-km long, and extends from 1 to 5 km below the surface to deep crustal levels. Isotope data and gravity models imply the fault zone separates crust of fundamentally different character. Geophysical evidence for such a long-lived structure, likely inherited from continental breakup, defies conventional wisdom that structures this old have been destroyed by Cenozoic extensional processes. Moreover, the coincidence with the alignment of mineral deposits supports the assertion by many economic geologists that these alignments are indicators of buried regional structures.

Journal of Geophysical Research B: Solid Earth

Slip rate and earthquake recurrence along the central Septentrional fault, North American-Caribbean plate boundary, Dominican Republic

The Septentrional fault zone (SFZ) is the major North American-Caribbean, strike-slip, plate boundary fault at the longitude of eastern Hispaniola. The SFZ traverses the densely populated Cibao Valley of the Dominican Republic, forming a prominent scarp in alluvium. Our studies at four sites along the central SFZ are aimed at quantifying the late Quaternary behavior of this structure to better understand the seismic hazard it represents for the northeastern Caribbean. Our investigations of excavations at sites near Rio Cenovi show that the most recent ground-rupturing earthquake along this fault in the north central Dominican Republic occurred between A.D. 1040 and A.D. 1230, and involved a minimum of ???4 m of left-lateral slip and 2.3 m of normal dip slip at that site. Our studies of offset stream terraces at two locations, Rio Juan Lopez and Rio Licey, provide late Holocene slip rate estimates of 6-9 mm/yr and a maximum of 11-12 mm/yr, respectively, across the Septentrional fault. Combining these results gives a best estimate of 6-12 mm/yr for the slip rate across the SFZ. Three excavations, two near Tenares and one at the Rio Licey site, yielded evidence for the occurrence of earlier prehistoric earthquakes. Dates of strata associated with the penultimate event suggest that it occurred post-A.D. 30, giving a recurrence interval of 800-1200 years. These studies indicate that the SFZ has likely accumulated elastic strain sufficient to generate a major earthquake during the more than 800 years since it last slipped and should be considered likely to produce a destructive future earthquake.

Journal of Geophysical Research B: Solid Earth

Three-dimensional crustal structure of Long Valley caldera, California, and evidence for the migration of CO2 under Mammoth Mountain

A temporary network of 69 three-component seismic stations captured a major seismic sequence in Long Valley caldera in 1997. We performed a tomographic inversion for crustal structure beneath a 28 km ?? 16 km area encompassing part of the resurgent dome, the south moat, and Mammoth Mountain. Resolution of crustal structure beneath the center of the study volume was good down to ???3 km below sea level (???5 km below the surface). Relatively high wave speeds are associated with the Bishop Tuff and lower wave speeds characterize debris in the surrounding moat. A low-Vp/Vs anomaly extending from near the surface to ???1 km below sea level beneath Mammoth Mountain may represent a CO2 reservoir that is supplying CO2-rich springs, venting at the surface, and killing trees. We investigated temporal variations in structure beneath Mammoth Mountain by differencing our results with tomographic images obtained using data from 1989/1990. Significant changes in both Vp and Vs were consistent with the migration of CO2 into the upper 2 km or so beneath Mammoth Mountain and its depletion in peripheral volumes that correlate with surface venting areas. Repeat tomography is capable of detecting the migration of gas beneath active silicic volcanoes and may thus provide a useful volcano monitoring tool.

Journal of Geophysical Research B: Solid Earth

Links between atmosphere, ocean, and cryosphere from two decades of microseism observations on the Antarctic Peninsula

The lack of landmasses, climatological low pressure, and strong circumpolar westerly winds between the latitudes of 50°S to 65°S produce exceptional storm‐driven wave conditions in the Southern Ocean. This combination makes the Antarctic Peninsula one of Earth's most notable regions of high‐amplitude wave activity and thus, ocean‐swell‐driven microseism noise in both the primary (direct wave‐coastal region interactions) and secondary (direct ocean floor forcing due to interacting wave trains) period bands. Microseism observations are examined across 23 years (1993–2015) from Palmer Station (PMSA), on the west coast of the Antarctic Peninsula, and from East Falkland Island (EFI). These records provide a spatially integrative measure of both Southern Ocean wave amplitudes and the interactions between ocean waves and the solid Earth in the presence of sea ice, which can reduce wave coupling with the continental shelf. We utilize a spatiotemporal correlation‐based approach to illuminate how the distribution of sea ice influences seasonal microseism power. We characterize primary and secondary microseism power due to variations in sea ice and find that primary microseism energy is both more sensitive to sea ice and more capable of propagating across ocean basins than secondary microseism energy. During positive phases of the Southern Annular Mode, sea ice is reduced in the Bellingshausen Sea and overall storm activity in the Drake Passage increases, thus strongly increasing microseism power levels.

Antarctic Penninsula

Chemical bonding in the outer core: high-pressure electronic structures of oxygen and sulfur in metallic iron

From its density the outer core is believed to be an alloy of iron and a light element such as sulfur or oxygen. The nature of the light element in the core is an important constraint for theories of the Earth's formation. In this paper the electronic structures of oxygen and sulfur impurities in metallic iron are investigated to determine if pressure, temperature, and composition-induced changes in bonding might affect phase equilibria along the Fe-FeS and Fe-FeO binaries. The electronic structure of sulfur in metallic iron is consistent with the miscibility between Fe and FeS liquids. Volume compression strengthens the Fe-S bond, and it is expected that at sufficiently high pressure, sulfur can substitute for Fe and give solid solution behavior between Fe and FeS. In contrast, the electronic structure of oxygen in metallic iron shows that oxygen cannot act as a substitutional impurity (replacing Fe). This explains the observed miscibility gap on the Fe-FeO binary at 1 atm pressure. Volume compression does not greatly change the electronic structure if oxygen substitutes for iron in bcc and fcc iron. Iron-oxygen bonding does occur, however, if oxygen occupies interstitial sites. Insofar as the molar volume of FeO incorporated as interstitial oxygen in metallic iron is smaller than that of pure FeO, the incorporation of oxygen into metallic iron may be favored under the pressures of the Earth's core.

Journal of Geophysical Research

Hydrologically induced deformation in Long Valley Caldera and adjacent Sierra Nevada

Vertical and horizontal components of GNSS displacements in the Long Valley Caldera and adjacent Sierra Nevada range show a clear correlation with hydrological trends at both multiyear and seasonal time scales. We observe a clear vertical and horizontal seasonal deformation pattern primarily attributable to the solid earth response to hydrological surface loading at large-to-regional (Sierra Nevada range) scales. Several GNSS sites, located at the eastern edge of the Sierra Nevada along the southwestern rim of Long Valley Caldera, also show significant horizontal deformation that cannot be explained by elastic deformation from surface loading. Due to the location of these sites and the strong correlation between their horizontal displacements and spring discharge, we hypothesize that this deformation reflects poroelastic processes related to snowmelt runoff water infiltrating into the Sierra Nevada slopes around Long Valley Caldera. Interestingly, this is also an area where water infiltrates to feed the local hydrothermal system, and where snowmelt-induced earthquake swarms have been recently detected.

California, Nevada

Impact crater outflows on Venus: Morphology and emplacement mechanisms

Many of the 932 impact craters discovered by the Magellan spacecraft at Venus are associated with lobate flows that originate at or near the crater rim. They extend for several to several hundred kilometers from the crater, and they commonly have a strong radar backscatter. A morphologic study of all identifiable crater outflows on Venus has revealed that many individual flows each consist of two areas, defined by distinct morphologic features. These two areas appear to represent two stages of deposition for each flow. The part of the flow that is generally deposited closest to the crater tends to be on the downrange side of the crater, flows in the downrange direction, and it is interpreted to be a late-stage ejecta. In many cases, this proximal part of the flow is too thin to completely bury the large blocks in subjacent ejecta deposits. Dendritic channels, present in many proximal flows, appear to have drained liquid from the proximal part in the downhill direction, and they debouch to feed the outer part of the flows. This distal part flows downhill, fills small grabens, and is ponded by ridges, behavior that mimics that of volcanic lava flows. The meandering and dendritic channels and the relation of the distal flows to topography strongly suggest that the distal portion is the result of coalescence and slow drainage of impact melt from the proximal portion. Impact melt forms a lining to the transient crater and mixes turbulently with solid clasts, and part of this mixture may be ejected to form the proximal part of the flow during the excavation stage of crater development. A statistical study of the Venusian craters has revealed that, in general, large craters produced by impacts with relatively low incidence angles to the surface are more likely to produce flows than small craters produced by higher-angle impacts. The greater flow production and downrange focusing of the proximal flows with decreasing incidence angle indicate a strong control of the flows by the impactor flight direction, and a high downrange velocity imparted to the proximal flow material in lower angle impacts. On the Moon, small flows interpreted to be composed of impact melt are observed atop the ejecta of large, fresh craters; on Earth, melt-rich suevite deposits form the uppermost layer of ejecta of some fresh craters. These features, albeit much smaller, may be analogous to the flows on Venus. Numerical models have predicted that larger volumes of impact melt would be produced on Venus than on the cooler terrestrial bodies due to high atmospheric and target temperatures, perhaps 3 times the volume produced on the Moon for a given crater diameter.

Journal of Geophysical Research E: Planets

Compositional variations on the Moon: Recalibration of Galileo solid‐state imaging data for the Orientale region and farside

Updated radiometric calibration and systematic processing procedures for Galileo solid‐state imaging (SSI) data from the first (1990) Earth‐Moon encounter are presented. These procedures were applied to a whole‐disk imaging sequence of the Moon centered near Mare Orientale, called Lunmap 14 (L14). Processing of L14 data included radiometric calibration, subpixel coregistration, scattered light removal, geometric control and reprojection, photometric normalization, and calibration to Earth‐based spectra. Coregistration and scattered‐light removal procedures are improvements over the initial calibration of the SSI mosaics. The effects of scattered light correction are best seen using a whole‐Moon view such as L14; resolution of the debated amounts of light scattering from within or outside the camera field of view is not necessary. Scattered light removal particularly affects the 1‐μm spectral region and has implications for interpretation of mafic mineral signatures in mare deposits of the lunar limb and farside. Recalibrated spectra indicate that mare ponds of the limb show moderately deep 1‐μm absorptions, and thus mafic mineral contents, comparable to those of other nearside basalts. Mafic mineral contents of Schiller‐Shickard cryptomaria are higher than previously thought and are similar to some low‐Ti nearside basalts. Many of the recalibrated spectra from South Pole/Aitken are similar to those of Schiller‐Schickard cryptomaria, suggesting that many of these soils represent a mixed mare/highland lithology. The hypothesis that there is an olivine enrichment in southern South Pole/Aitken basin is not supported strongly by spectra shown here.

Journal of Geophysical Research: Planets

Advances in the simulation and automated measurement of well-sorted granular material: 1. Simulation

1. In this, the first of a pair of papers which address the simulation and automated measurement of well-sorted natural granular material, a method is presented for simulation of two-phase (solid, void) assemblages of discrete non-cohesive particles. The purpose is to have a flexible, yet computationally and theoretically simple, suite of tools with well constrained and well known statistical properties, in order to simulate realistic granular material as a discrete element model with realistic size and shape distributions, for a variety of purposes. The stochastic modeling framework is based on three-dimensional tessellations with variable degrees of order in particle-packing arrangement. Examples of sediments with a variety of particle size distributions and spatial variability in grain size are presented. The relationship between particle shape and porosity conforms to published data. The immediate application is testing new algorithms for automated measurements of particle properties (mean and standard deviation of particle sizes, and apparent porosity) from images of natural sediment, as detailed in the second of this pair of papers. The model could also prove useful for simulating specific depositional structures found in natural sediments, the result of physical alterations to packing and grain fabric, using discrete particle flow models. While the principal focus here is on naturally occurring sediment and sedimentary rock, the methods presented might also be useful for simulations of similar granular or cellular material encountered in engineering, industrial and life sciences.

Journal of Geophysical Research F: Earth Surface