USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysical Research - Planets”

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.

At least 271 records · Page 15Linked to original sources
USGS2017

Quantifying the heterogeneity of the tectonic stress field using borehole data

The heterogeneity of the tectonic stress field is a fundamental property which influences earthquake dynamics and subsurface engineering. Self-similar scaling of stress heterogeneities is frequently assumed to explain characteristics of earthquakes such as the magnitude-frequency relation. However, observational evidence for such scaling of the stress field heterogeneity is scarce. We analyze the local stress orientations using image logs of two closely spaced boreholes in the Coso Geothermal Field with sub-vertical and deviated trajectories, respectively, each spanning about 2 km in depth. Both the mean and the standard deviation of stress orientation indicators (borehole breakouts, drilling-induced fractures and petal-centerline fractures) determined from each borehole agree to the limit of the resolution of our method although measurements at specific depths may not. We find that the standard deviation in these boreholes strongly depends on the interval length analyzed, generally increasing up to a wellbore log length of about 600 m and constant for longer intervals. We find the same behavior in global data from the World Stress Map. This suggests that the standard deviation of stress indicators characterizes the heterogeneity of the tectonic stress field rather than the quality of the stress measurement. A large standard deviation of a stress measurement might be an expression of strong crustal heterogeneity rather than of an unreliable stress determination. Robust characterization of stress heterogeneity requires logs that sample stress indicators along a representative sample volume of at least 1 km.

Journal of Geophysical Research B: Solid Earth↗
USGS2017

A decade of induced slip on the causative fault of the 2015 Mw 4.0 Venus earthquake, northeast Johnson County, Texas

On 7 May 2015, a M w 4.0 earthquake occurred near Venus, northeast Johnson County, Texas, in an area of the Bend Arch-Fort Worth Basin that reports long-term, high-volume wastewater disposal and that has hosted felt earthquakes since 2009. In the weeks following the M w 4.0 earthquake, we deployed a local seismic network and purchased nearby active-source seismic reflection data to capture additional events, characterize the causative fault, and explore potential links between ongoing industry activity and seismicity. Hypocenter relocations of the resulting local earthquake catalog span ~4–6 km depth and indicate a fault striking ~230 ° , dipping to the west, consistent with a nodal plane of the M w 4.0 regional moment tensor. Fault plane solutions indicate normal faulting, with B axes striking parallel to maximum horizontal compressive stress. Seismic reflection data image the reactivated basement fault penetrating the Ordovician disposal layer and Mississippian production layer, but not displacing post-Lower Pennsylvanian units. Template matching at regional seismic stations indicates that low-magnitude earthquakes with similar waveforms began in April 2008, with increasing magnitude over time. Pressure data from five saltwater disposal wells within 5 km of the active fault indicate a disposal formation that is 0.9–4.8 MPa above hydrostatic. We suggest that the injection of 28,000,000 m 3 of wastewater between 2006 and 2015 at these wells led to an increase in subsurface pore fluid pressure that contributed to inducing this long-lived earthquake sequence. The 2015 M w 4.0 event represents the largest event in the continuing evolution of slip on the causative fault.

Texas↗
USGS2017

Magnetotelluric imaging of lower crustal melt and lithospheric hydration in the Rocky Mountain Front transition zone, Colorado, USA

We present an electrical resistivity model of the crust and upper mantle from two‐dimensional (2‐D) anisotropic inversion of magnetotelluric data collected along a 450 km transect of the Rio Grande rift, southern Rocky Mountains, and High Plains in Colorado, USA. Our model provides a window into the modern‐day lithosphere beneath the Rocky Mountain Front to depths in excess of 150 km. Two key features of the 2‐D resistivity model are (1) a broad zone (~200 km wide) of enhanced electrical conductivity (<20 Ωm) in the midcrust to lower crust that is centered beneath the highest elevations of the southern Rocky Mountains and (2) hydrated lithospheric mantle beneath the Great Plains with water content in excess of 100 ppm. We interpret the high conductivity region of the lower crust as a zone of partially molten basalt and associated deep‐crustal fluids that is the result of recent (less than 10 Ma) tectonic activity in the region. The recent supply of volatiles and/or heat to the base of the crust in the late Cenozoic implies that modern‐day tectonic activity in the western United States extends to at least the western margin of the Great Plains. The transition from conductive to resistive upper mantle is caused by a gradient in lithospheric modification, likely including hydration of nominally anhydrous minerals, with maximum hydration occurring beneath the Rocky Mountain Front. This lithospheric “hydration front” has implications for the tectonic evolution of the continental interior and the mechanisms by which water infiltrates the lithosphere.

Colorado↗
USGS2017

Decaying lava extrusion rate at El Reventador Volcano, Ecuador measured using high-resolution satellite radar

Lava extrusion at erupting volcanoes causes rapid changes in topography and morphology on the order of tens or even hundreds of meters. Satellite radar provides a method for measuring changes in topographic height over a given time period to an accuracy of meters, either by measuring the width of radar shadow cast by steep sided features, or by measuring the difference in radar phase between two sensors separated in space. We measure height changes, and hence estimate extruded lava volume flux, at El Reventador, Ecuador, between 2011 and 2016, using data from the RADARSAT-2 and TanDEM-X satellite missions. We find that 39 new lava flows were extruded between 9 February 2012 and 24 August 2016, with a cumulative volume of 44.8M m 3 dense rock equivalent and a gradually decreasing eruption rate. The average dense rock rate of lava extrusion during this time is 0.31 ± 0.02 m 3 s −1 , which is similar to the long-term average from 1972 to 2016. Apart from a volumetrically small dyke opening event between 9 March and 10 June 2012, lava extrusion at El Reventador is not accompanied by any significant magmatic ground deformation. We use a simple physics-based model to estimate that the volume of the magma reservoir under El Reventador is greater than 3 km 3 . Our lava extrusion data can be equally well fit by models representing a closed reservoir depressurising during the eruption with no magma recharge, or an open reservoir with a time-constant magma recharge rate of up to 0.35 ± 0.01 m 3 s −1 .

El Reventador Volcano↗
USGS2018

Temporal stress changes caused by earthquakes: A review

Earthquakes can change the stress field in the Earth’s lithosphere as they relieve and redistribute stress. Earthquake-induced stress changes have been observed as temporal rotations of the principal stress axes following major earthquakes in a variety of tectonic settings. The stress changes due to the 2011 Mw9.0 Tohoku-Oki, Japan, earthquake were particularly well documented. Earthquake stress rotations can inform our understanding of earthquake physics, most notably addressing the long-standing problem of whether the Earth’s crust at plate boundaries is “strong” or “weak.” Many of the observed stress rotations, including that due to the Tohoku-Oki earthquake, indicate near-complete stress drop in the mainshock. This implies low background differential stress, on the order of earthquake stress drop, supporting the weak crust model. Earthquake stress rotations can also be used to address other important geophysical questions, such as the level of crustal stress heterogeneity and the mechanisms of postseismic stress reloading. The quantitative interpretation of stress rotations is evolving from those based on simple analytical methods to those based on more sophisticated numerical modeling that can capture the spatial-temporal complexity of the earthquake stress changes.

Journal of Geophysical Research B: Solid Earth↗
USGS2017

Source spectral properties of small-to-moderate earthquakes in southern Kansas

The source spectral properties of injection-induced earthquakes give insight into their nucleation, rupture processes, and influence on ground motion. Here we apply a spectral decomposition approach to analyze P -wave spectra and estimate Brune-type stress drop for more than 2000 M L 1.5–5.2 earthquakes occurring in southern Kansas from 2014 to 2016. We find that these earthquakes are characterized by low stress drop values (median ∼0.4MPa) compared to natural seismicity in California. We observe a significant increase in stress drop as a function of depth, but the shallow depth distribution of these events is not by itself sufficient to explain their lower stress drop. Stress drop increases with magnitude from M1.5–M3.5, but this scaling trend may weaken above M4 and also depends on the assumed source model. Although we observe a nonstationary, sequence-specific temporal evolution in stress drop, we find no clear systematic relation with the activity of nearby injection wells.

Kansas↗
USGS2018

Deformation of the Pacific/North America plate boundary at Queen Charlotte Fault: The possible role of rheology

The Pacific/North America (PA/NA) plate boundary between Vancouver Island and Alaska is similar to the PA/NA boundary in California in its kinematic history and the rate and azimuth of current relative motion, yet their deformation styles are distinct. The California plate boundary shows a broad zone of parallel strike slip and thrust faults and folds, whereas the 49‐mm/yr PA/NA relative plate motion in Canada and Alaska is centered on a single, narrow, continuous ~900‐km‐long fault, the Queen Charlotte Fault (QCF). Using gravity analysis, we propose that this plate boundary is centered on the continent/ocean boundary (COB), an unusual location for continental transform faults because plate boundaries typically localize within the continental lithosphere, which is weaker. Because the COB is a boundary between materials of contrasting elastic properties, once a fault is established there, it will probably remain stable. We propose that deformation progressively shifted to the COB in the wake of Yakutat terrane's northward motion along the margin. Minor convergence across the plate boundary is probably accommodated by fault reactivation on Pacific crust and by an eastward dipping QCF. Underthrusting of Pacific slab under Haida Gwaii occurs at convergence angles >14°–15° and may have been responsible for the emergence of the archipelago. The calculated slab entry dip (5°–8°) suggests that the slab probably does not extend into the asthenosphere. The PA/NA plate boundary at the QCF can serve as a structurally simple site to investigate the impact of rheology and composition on crustal deformation and the initiation of slab underthrusting.

Journal of Geophysical Research B: Solid Earth↗
USGS2018

Mass addition at Mount St. Helens, Washington, inferred from repeated gravity surveys

Measurements of subtle changes in the Earth’s gravityfield can provideinformation on the addition/loss of mass (e.g., magma or aqueousfluids) beneath a volcano. In this study, wemeasured gravity at Mount St. Helens from 2010 to 2016 to investigate possible mass changes followingthe 2004–2008 dome-forming eruption. The raw gravity measurements were corrected for changes in themass and shape of Crater Glacier and for deformation of the volcanic edifice to obtain residual gravity values.We found positive residual gravity changes that we interpret as evidence for partial recharge of the magmareservoir that fed the 2004–2008 eruption and/or for accumulation of groundwater in one or more shallowaquifers. Most of the gravity signal can be explained by groundwater accumulation; magma recharge is notrequired. However, there is additional evidence for at least partial recharge from seismic, deformation, andgeochemical observations

Washington↗
USGS2018

Constraints on friction, dilatancy, diffusivity, and effective stress from low-frequency earthquake rates on the deep San Andreas Fault

Families of recurring low-frequency earthquakes (LFEs) within nonvolcanic tremor on the San Andreas Fault in central California are sensitive to tidal stresses. LFEs occur at all levels of the tides, are strongly correlated and in phase with the ~200 Pa shear stresses, and weakly and not systematically correlated with the ~2 kPa tidal normal stresses. We assume that LFEs are small sources that repeatedly fail during shear within a much larger scale, aseismically slipping fault zone and consider two different models of the fault slip: (1) modulation of the fault slip rate by the tidal stresses or (2) episodic slip, triggered by the tides. LFEs are strongly clustered with duration much shorter than the semidiurnal tide; they cannot be significantly modulated on that time scale. The recurrence times of clusters, however, are many times longer than the semidiurnal, leading to an appearance of tidal triggering. In this context we examine the predictions of laboratory-observed triggered frictional (dilatant) fault slip. The undrained end-member model produces no sensitivity to the tidal normal stress, and slip onsets are in phase with the tidal shear stress. The tidal correlation constrains the diffusivity to be less than ~1 × 10 −6 /s and the product of the friction and dilatancy coefficients to be at most 5 × 10 −7 , orders of magnitude smaller than observed at room temperature. In the absence of dilatancy the effective normal stress at failure would be about ~55 kPa. For this model the observations require intrinsic weakness, low dilatancy, and lithostatic pore fluid.

Journal of Geophysical Research↗
USGS2018

Combining multiphase groundwater flow and slope stability models to assess stratovolcano flank collapse in the Cascade Range

Hydrothermal alteration can create low‐permeability zones, potentially resulting in elevated pore‐fluid pressures, within a volcanic edifice. Strength reduction by rock alteration and high pore‐fluid pressures have been suggested as a mechanism for edifice flank instability. Here we combine numerical models of multiphase heat transport and groundwater flow with a slope‐stability code that incorporates three‐dimensional distributions of strength and pore‐water pressure to address the following questions: (1) What permeability distributions and contrasts produce elevated pore‐fluid pressures in a stratovolcano? (2) What are the effects of these elevated pressures on flank stability? (3) Finally, what are the effects of magma intrusion on potential flank failure in an edifice? Simulation results show that under a range of plausible parameters, water tables in a stratovolcano can be elevated or perched. These elevated water tables result in universally lower stability (lower factor of safety) compared with equivalent dry edifices, indicating a higher likelihood of flank collapse. Low‐permeability (<1 × 10 −17 m 2 ) layers such as altered pyroclastic deposits or breccias can result in locally saturated regions (perched water) and lower factors of safety near the ground surface but may actually reduce liquid water saturation and pore pressures in the core of the edifice and thus may favor small, shallow collapses over larger, deeper collapses. Magma intrusion into the base of the edifice increases pore‐fluid pressures and decreases the factor of safety. However, the shear strength of edifice rocks also exerts a significant control on stability, so both mechanical properties and pore‐fluid pressures are important for stability assessments.

Washington↗
USGS2018

Seismicity induced by wastewater injection in Washington County, Ohio: Influence of preexisting structure, regional stress regime, and well operations

Recent seismicity in Washington County, Ohio, has been suggested to be induced by wastewater disposal operations despite injection ~2 km above the Precambrian basement. We investigated the relationships between disposal well locations and operational histories, spatiotemporal patterns of seismicity enhanced by waveform correlation, and mapped subsurface structures. We also analyzed proxies for in situ stress conditions to evaluate existing principal stress orientations and magnitudes. Double‐difference relocations of small ( M L ≤ 2.1) seismicity revealed linear clusters defining a NE‐SW orientation that parallels the trend of small‐amplitude folds in Upper Paleozoic rocks within 1–2 km of the injection well. Hypocenter depths of ~3.1 to 4.3 km indicate seismicity occurred along faults below the Silurian injection interval in lower Paleozoic strata and Precambrian crystalline basement. Similar folds/faults imaged in nearby seismic reflection lines were associated with basement‐involved fault systems that cut the injection interval targeted by the injection well, providing permeability pathways for fluid pressure increases that initiated slip. Interpreted faults were also optimally oriented in the regional stress field for reactivation due to injection‐related decreases in effective stress. Well completion records suggest that during the recorded seismicity, the Long Run‐1 well was injecting disposal fluids close to or above the fracture pressure of the targeted Silurian injection interval. A comparison of observed seismicity rates and monthly injection records from the disposal well also indicate a close correspondence between earthquake frequency and injection rate. Taken together, these results provide a physical framework for understanding why wastewater disposal operations induced the Washington County seismic sequence.

Ohio↗
USGS2017

Interactions of estuarine shoreline infrastructure with multiscale sea level variability

Sea level rise increases the risk of storms and other short‐term water‐rise events, because it sets a higher water level such that coastal surges become more likely to overtop protections and cause floods. To protect coastal communities, it is necessary to understand the interaction among multiday and tidal sea level variabilities, coastal infrastructure, and sea level rise. We performed a series of numerical simulations for San Francisco Bay to examine two shoreline scenarios and a series of short‐term and long‐term sea level variations. The two shoreline configurations include the existing topography and a coherent full‐bay containment that follows the existing land boundary with an impermeable wall. The sea level variability consists of a half‐meter perturbation, with duration ranging from 2 days to permanent (i.e., sea level rise). The extent of coastal flooding was found to increase with the duration of the high‐water‐level event. The nonlinear interaction between these intermediate scale events and astronomical tidal forcing only contributes ∼1% of the tidal heights; at the same time, the tides are found to be a dominant factor in establishing the evolution and diffusion of multiday high water events. Establishing containment at existing shorelines can change the tidal height spectrum up to 5%, and the impact of this shoreline structure appears stronger in the low‐frequency range. To interpret the spatial and temporal variability at a wide range of frequencies, Optimal Dynamic Mode Decomposition is introduced to analyze the coastal processes and an inverse method is applied to determine the coefficients of a 1‐D diffusion wave model that quantify the impact of bottom roughness, tidal basin geometry, and shoreline configuration on the high water events

Journal of Geophysical Research: Oceans↗
USGS2017

Persistent shoreline shape induced from offshore geologic framework: Effects of shoreface connected ridges

Mechanisms relating offshore geologic framework to shoreline evolution are determined through geologic investigations, oceanographic deployments, and numerical modeling. Analysis of shoreline positions from the past 50 years along Fire Island, New York, a 50 km long barrier island, demonstrates a persistent undulating shape along the western half of the island. The shelf offshore of these persistent undulations is characterized with shoreface-connected sand ridges (SFCR) of a similar alongshore length scale, leading to a hypothesis that the ridges control the shoreline shape through the modification of flow. To evaluate this, a hydrodynamic model was configured to start with the US East Coast and scale down to resolve the Fire Island nearshore. The model was validated using observations along western Fire Island and buoy data, and used to compute waves, currents and sediment fluxes. To isolate the influence of the SFCR on the generation of the persistent shoreline shape, simulations were performed with a linearized nearshore bathymetry to remove alongshore transport gradients associated with shoreline shape. The model accurately predicts the scale and variation of the alongshore transport that would generate the persistent shoreline undulations. In one location, however, the ridge crest connects to the nearshore and leads to an offshore-directed transport that produces a difference in the shoreline shape. This qualitatively supports the hypothesized effect of cross-shore fluxes on coastal evolution. Alongshore flows in the nearshore during a representative storm are driven by wave breaking, vortex force, advection and pressure gradient, all of which are affected by the SFCR.

New York↗
USGS2017

A Bayesian-based system to assess wave-driven flooding hazards on coral reef-lined coasts

Many low-elevation, coral reef-lined, tropical coasts are vulnerable to the effects of climate change, sea level rise, and wave-induced flooding. The considerable morphological diversity of these coasts and the variability of the hydrodynamic forcing that they are exposed to make predicting wave-induced flooding a challenge. A process-based wave-resolving hydrodynamic model (XBeach Non-Hydrostatic, “XBNH”) was used to create a large synthetic database for use in a “Bayesian Estimator for Wave Attack in Reef Environments” (BEWARE), relating incident hydrodynamics and coral reef geomorphology to coastal flooding hazards on reef-lined coasts. Building on previous work, BEWARE improves system understanding of reef hydrodynamics by examining the intrinsic reef and extrinsic forcing factors controlling runup and flooding on reef-lined coasts. The Bayesian estimator has high predictive skill for the XBNH model outputs that are flooding indicators, and was validated for a number of available field cases. It was found that, in order to accurately predict flooding hazards, water depth over the reef flat, incident wave conditions, and reef flat width are the most essential factors, whereas other factors such as beach slope and bed friction due to the presence or absence of corals are less important. BEWARE is a potentially powerful tool for use in early warning systems or risk assessment studies, and can be used to make projections about how wave-induced flooding on coral reef-lined coasts may change due to climate change.

Journal of Geophysical Research C: Oceans↗
USGS2018

Direct measurements of mean Reynolds stress and ripple roughness in the presence of energetic forcing by surface waves

Direct covariance observations of the mean flow Reynolds stress and sonar images of the seafloor collected on a wave‐exposed inner continental shelf demonstrate that the drag exerted by the seabed on the overlying flow is consistent with boundary layer models for wave‐current interaction, provided that the orientation and anisotropy of the bed roughness are appropriately quantified. Large spatial and temporal variations in drag result from nonequilibrium ripple dynamics, ripple anisotropy, and the orientation of the ripples relative to the current. At a location in coarse sand characterized by large two‐dimensional orbital ripples, the observed drag shows a strong dependence on the relative orientation of the mean current to the ripple crests. At a contrasting location in fine sand, where more isotropic sub‐orbital ripples are observed, the sensitivity of the current to the orientation of the ripples is reduced. Further, at the coarse site under conditions when the currents are parallel to the ripple crests and the wave orbital diameter is smaller than the wavelength of the relic orbital ripples, the flow becomes hydraulically smooth. This transition is not observed at the fine site, where the observed wave orbital diameter is always greater than the wavelength of the observed sub‐orbital ripples. Paradoxically, the dominant along‐shelf flows often experience lower drag at the coarse site than at the fine site, despite the larger ripples, highlighting the complex dynamics controlling drag in wave‐exposed environments with heterogeneous roughness.

Journal of Geophysical Research C: Oceans↗
USGS2017

Vulnerability of coral reefs to bioerosion from land-based sources of pollution

Ocean acidification (OA), the gradual decline in ocean pH and [ ] caused by rising levels of atmospheric CO 2 , poses a significant threat to coral reef ecosystems, depressing rates of calcium carbonate (CaCO 3 ) production, and enhancing rates of bioerosion and dissolution. As ocean pH and [ ] decline globally, there is increasing emphasis on managing local stressors that can exacerbate the vulnerability of coral reefs to the effects of OA. We show that sustained, nutrient rich, lower pH submarine groundwater discharging onto nearshore coral reefs off west Maui lowers the pH of seawater and exposes corals to nitrate concentrations 50 times higher than ambient. Rates of coral calcification are substantially decreased, and rates of bioerosion are orders of magnitude higher than those observed in coral cores collected in the Pacific under equivalent low pH conditions but living in oligotrophic waters. Heavier coral nitrogen isotope (δ 15 N) values pinpoint not only site-specific eutrophication, but also a sewage nitrogen source enriched in 15 N. Our results show that eutrophication of reef seawater by land-based sources of pollution can magnify the effects of OA through nutrient driven-bioerosion. These conditions could contribute to the collapse of coastal coral reef ecosystems sooner than current projections predict based only on ocean acidification.

Hawaii↗
USGS2018

Investigating the origin of continual radio frequency impulses during explosive volcanic eruptions

Volcanic lightning studies have revealed that there is a relatively long‐lasting source of very high frequency radiation associated with the onset of explosive volcanic eruptions that is distinct from radiation produced by lightning. This very high frequency signal is referred to as “continual radio frequency (CRF)” due to its long‐lasting nature. The discharge mechanism producing this signal was previously hypothesized to be caused by numerous, small (10–100 m) leader‐forming discharges near the vent of the volcano. To test this hypothesis, a multiparametric data set of electrical and volcanic activity occurring during explosive eruptions of Sakurajima Volcano in Japan was collected from May to June 2015. Our observations show that a single CRF impulse has a duration on the order of 160 ns (giving an upper limit on discharge length of 10 m) and is distinct from near‐vent lightning discharges that were on the order of 30 m in length. CRF impulses did not produce discernible electric field changes and occurred in the absence of a net static electric field. Lightning mapping data and infrared video observations of the eruption column showed that CRF impulses originated from the gas thrust region of the column. These observations indicate that CRF impulses are not produced by small, leader‐forming discharges but rather are more similar to a streamer discharge, likely on the order of a few meters in length.

Sakurajima Volcano↗
USGS2017

Presentation and analysis of a worldwide database of earthquake-induced landslide inventories

Earthquake-induced landslide (EQIL) inventories are essential tools to extend our knowledge of the relationship between earthquakes and the landslides they can trigger. Regrettably, such inventories are difficult to generate and therefore scarce, and the available ones differ in terms of their quality and level of completeness. Moreover, access to existing EQIL inventories is currently difficult because there is no centralized database. To address these issues, we compiled EQIL inventories from around the globe based on an extensive literature study. The database contains information on 363 landslide-triggering earthquakes and includes 66 digital landslide inventories. To make these data openly available, we created a repository to host the digital inventories that we have permission to redistribute through the U.S. Geological Survey ScienceBase platform. It can grow over time as more authors contribute their inventories. We analyze the distribution of EQIL events by time period and location, more specifically breaking down the distribution by continent, country, and mountain region. Additionally, we analyze frequency distributions of EQIL characteristics, such as the approximate area affected by landslides, total number of landslides, maximum distance from fault rupture zone, and distance from epicenter when the fault plane location is unknown. For the available digital EQIL inventories, we examine the underlying characteristics of landslide size, topographic slope, roughness, local relief, distance to streams, peak ground acceleration, peak ground velocity, and Modified Mercalli Intensity. Also, we present an evaluation system to help users assess the suitability of the available inventories for different types of EQIL studies and model development.

Journal of Geophysical Research F: Earth Surface↗