SEARCH · USGS Science

Results for “Journal of Geophysical Research - Atmospheres”

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

SourceResearch connectionFrom yearThrough year

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

At least 145 records · Page 8Linked to original sources

USGS2014

Radiocarbon age-offsets in an arctic lake reveal the long-term response of permafrost carbon to climate change

Continued warming of the Arctic may cause permafrost to thaw and speed the decomposition of large stores of soil organic carbon (OC), thereby accentuating global warming. However, it is unclear if recent warming has raised the current rates of permafrost OC release to anomalous levels or to what extent soil carbon release is sensitive to climate forcing. Here we use a time series of radiocarbon age-offsets ( 14 C) between the bulk lake sediment and plant macrofossils deposited in an arctic lake as an archive for soil and permafrost OC release over the last 14,500 years. The lake traps and archives OC imported from the watershed and allows us to test whether prior warming events stimulated old carbon release and heightened age-offsets. Today, the age-offset (2 ka; thousand of calibrated years before A.D. 1950) and the depositional rate of ancient OC from the watershed into the lake are relatively low and similar to those during the Younger Dryas cold interval (occurring 12.9–11.7 ka). In contrast, age-offsets were higher (3.0–5.0 ka) when summer air temperatures were warmer than present during the Holocene Thermal Maximum (11.7–9.0 ka) and Bølling-Allerød periods (14.5–12.9 ka). During these warm times, permafrost thaw contributed to ancient OC depositional rates that were ~10 times greater than today. Although permafrost OC was vulnerable to climate warming in the past, we suggest surface soil organic horizons and peat are presently limiting summer thaw and carbon release. As a result, the temperature threshold to trigger widespread permafrost OC release is higher than during previous warming events.

Alaska↗

USGS2014

Using dissolved organic matter age and composition to detect permafrost thaw in boreal watersheds of interior Alaska

Recent warming at high latitudes has accelerated permafrost thaw, which can modify soil carbon dynamics and watershed hydrology. The flux and composition of dissolved organic matter (DOM) from soils to rivers are sensitive to permafrost configuration and its impact on subsurface hydrology and groundwater discharge. Here, we evaluate the utility of DOM composition and age as a tool for detecting permafrost thaw in three rivers (Beaver, Birch, and Hess Creeks) within the discontinuous permafrost zone of interior Alaska. We observed strong temporal controls on Δ 14 C content of hydrophobic acid isolates (Δ 14 C-HPOA) across all rivers, with the most enriched values occurring during spring snowmelt (75 ± 8&permil;) and most depleted during winter flow (−21 ± 8&permil;). Radiocarbon ages of winter flow samples ranged from 35 to 445 yr BP, closely tracking estimated median base flow travel times for this region (335 years). During spring snowmelt, young DOM was composed of highly aromatic, high molecular-weight compounds, whereas older DOM of winter flow had lower aromaticity and molecular weight. We observed a significant correlation between Δ 14 C-HPOA and UV absorbance coefficient at 254 nm ( α 254 ) across all study rivers. Using α 254 as an optical indicator for Δ 14 C-HPOA, we also observed a long-term decline in α 254 during maximum annual thaw depth over the last decade at the Hess Creek study site. These findings suggest a shift in watershed hydrology associated with increasing active layer thickness. Further development of DOM optical indicators may serve as a novel and inexpensive tool for detecting permafrost degradation in northern watersheds.

Alaska↗

USGS2014

Uncertainty in the fate of soil organic carbon: A comparison of three conceptually different soil decomposition models

Conventional Q10 soil organic matter decomposition models and more complex microbial models are available for making projections of future soil carbon dynamics. However, it is unclear (1) how well the conceptually different approaches can simulate observed decomposition and (2) to what extent the trajectories of long-term simulations differ when using the different approaches. In this study, we compared three structurally different soil carbon (C) decomposition models (one Q10 and two microbial models of different complexity), each with a one- and two-horizon version. The models were calibrated and validated using 4 years of measurements of heterotrophic soil CO 2 efflux from trenched plots in a Dahurian larch ( Larix gmelinii Rupr.) plantation. All models reproduced the observed heterotrophic component of soil CO 2 efflux, but the trajectories of soil carbon dynamics differed substantially in 100 year simulations with and without warming and increased litterfall input, with microbial models that produced better agreement with observed changes in soil organic C in long-term warming experiments. Our results also suggest that both constant and varying carbon use efficiency are plausible when modeling future decomposition dynamics and that the use of a short-term (e.g., a few years) period of measurement is insufficient to adequately constrain model parameters that represent long-term responses of microbial thermal adaption. These results highlight the need to reframe the representation of decomposition models and to constrain parameters with long-term observations and multiple data streams. We urge caution in interpreting future soil carbon responses derived from existing decomposition models because both conceptual and parameter uncertainties are substantial.

Journal of Geophysical Research G: Biogeosciences↗

USGS2014

Nitrogen cycling processes and microbial community composition in bed sediments in the Yukon River at Pilot Station

Information on the contribution of nitrogen (N)-cycling processes in bed sediments to river nutrient fluxes in large northern latitude river systems is limited. This study examined the relationship between N-cycling processes in bed sediments and N speciation and loading in the Yukon River near its mouth at the Bering Sea. We conducted laboratory bioassays to measure N-cycling processes in sediment samples collected over distinct water cycle seasons. In conjunction, the microbial community composition in the bed sediments using genes involved in N-cycling ( narG , napA , nosZ , and amoA ) and 16S rRNA gene pyrosequences was examined. Temporal variation was observed in net N mineralization, nitrate uptake, and denitrification rate potentials and correlated strongly with sediment carbon (C) and extractable N content and microbial community composition rather than with river water nutrient concentrations. The C content of the bed sediment was notably impacted by the spring flood, ranging from 1.1% in the midst of an ice-jam to 0.1% immediately after ice-out, suggesting a buildup of organic material (OM) prior to scouring of the bed sediments during ice break up. The dominant members of the microbial community that explained differences in N-processing rates belonged to the genera Crenothrix , Flavobacterium , and the family of Comamonadaceae. Our results suggest that biogeochemical processing rates in the bed sediments appear to be more coupled to hydrology, nutrient availability in the sediments, and microbial community composition rather than river nutrient concentrations at Pilot Station.

Alaska↗

USGS2015

Complex terrain alters temperature and moisture limitations of forest soil respiration across a semiarid to subalpine gradient

Forest soil respiration is a major carbon (C) flux that is characterized by significant variability in space and time. We quantified growing season soil respiration during both a drought year and a nondrought year across a complex landscape to identify how landscape and climate interact to control soil respiration. We asked the following questions: (1) How does soil respiration vary across the catchments due to terrain-induced variability in moisture availability and temperature? (2) Does the relative importance of moisture versus temperature limitation of respiration vary across space and time? And (3) what terrain elements are important for dictating the pattern of soil respiration and its controls? Moisture superseded temperature in explaining watershed respiration patterns, with wetter yet cooler areas higher up and on north facing slopes yielding greater soil respiration than lower and south facing areas. Wetter subalpine forests had reduced moisture limitation in favor of greater seasonal temperature limitation, and the reverse was true for low-elevation semiarid forests. Coincident climate poorly predicted soil respiration in the montane transition zone; however, antecedent precipitation from the prior 10 days provided additional explanatory power. A seasonal trend in respiration remained after accounting for microclimate effects, suggesting that local climate alone may not adequately predict seasonal variability in soil respiration in montane forests. Soil respiration climate controls were more strongly related to topography during the drought year highlighting the importance of landscape complexity in ecosystem response to drought.

Journal of Geophysical Research: Biogeosciences↗

USGS2015

Soil greenhouse gas emissions and carbon budgeting in a short-hydroperiod floodplain wetland

Understanding the controls on floodplain carbon (C) cycling is important for assessing greenhouse gas emissions and the potential for C sequestration in river-floodplain ecosystems. We hypothesized that greater hydrologic connectivity would increase C inputs to floodplains that would not only stimulate soil C gas emissions but also sequester more C in soils. In an urban Piedmont river (151 km 2 watershed) with a floodplain that is dry most of the year, we quantified soil CO 2 , CH 4 , and N 2 O net emissions along gradients of floodplain hydrologic connectivity, identified controls on soil aerobic and anaerobic respiration, and developed a floodplain soil C budget. Sites were chosen along a longitudinal river gradient and across lateral floodplain geomorphic units (levee, backswamp, and toe slope). CO 2 emissions decreased downstream in backswamps and toe slopes and were high on the levees. CH 4 and N 2 O fluxes were near zero; however, CH 4 emissions were highest in the backswamp. Annual CO 2 emissions correlated negatively with soil water-filled pore space and positively with variables related to drier, coarser soil. Conversely, annual CH 4 emissions had the opposite pattern of CO 2 . Spatial variation in aerobic and anaerobic respiration was thus controlled by oxygen availability but was not related to C inputs from sedimentation or vegetation. The annual mean soil CO 2 emission rate was 1091 g C m −2  yr −1 , the net sedimentation rate was 111 g C m −2  yr −1 , and the vegetation production rate was 240 g C m −2  yr −1 , with a soil C balance (loss) of −338 g C m −2  yr −1 . This floodplain is losing C likely due to long-term drying from watershed urbanization.

Virginia↗

USGS2015

Cascadia subducting plate fluids channelled to fore-arc mantle corner: ETS and silica deposition

In this study we first summarize the constraints that on the Cascadia subduction thrust, there is a 70 km gap downdip between the megathrust seismogenic zone and the Episodic Tremor and Slip (ETS) that lies further landward; there is not a continuous transition from unstable to conditionally stable sliding. Seismic rupture occurs mainly offshore for this hot subduction zone. ETS lies onshore. We then suggest what does control the downdip position of ETS. We conclude that fluids from dehydration of the downgoing plate, focused to rise above the fore-arc mantle corner, are responsible for ETS. There is a remarkable correspondence between the position of ETS and this corner along the whole margin. Hydrated mineral assemblages in the subducting oceanic crust and uppermost mantle are dehydrated with downdip increasing temperature, and seismic tomography data indicate that these fluids have strongly serpentinized the overlying fore-arc mantle. Laboratory data show that such fore-arc mantle serpentinite has low permeability and likely blocks vertical expulsion and restricts flow updip within the underlying permeable oceanic crust and subduction shear zone. At the fore-arc mantle corner these fluids are released upward into the more permeable overlying fore-arc crust. An indication of this fluid flux comes from low Poisson's Ratios (and V p / V s ) found above the corner that may be explained by a concentration of quartz which has exceptionally low Poisson's Ratio. The rising fluids should be silica saturated and precipitate quartz with decreasing temperature and pressure as they rise above the corner.

Journal of Geophysical Research↗

USGS2015

The Palos Verdes Fault offshore southern California: late Pleistocene to present tectonic geomorphology, seascape evolution and slip rate estimate based on AUV and ROV surveys

The Palos Verdes Fault (PVF) is one of few active faults in Southern California that crosses the shoreline and can be studied using both terrestrial and subaqueous methodologies. To characterize the near-seafloor fault morphology, tectonic influences on continental slope sedimentary processes and late Pleistocene to present slip rate, a grid of high-resolution multibeam bathymetric data, and chirp subbottom profiles were acquired with an autonomous underwater vehicle (AUV) along the main trace of PVF in water depths between 250 and 600 m. Radiocarbon dates were obtained from vibracores collected using a remotely operated vehicle (ROV) and ship-based gravity cores. The PVF is expressed as a well-defined seafloor lineation marked by subtle along-strike bends. Right-stepping transtensional bends exert first-order control on sediment flow dynamics and the spatial distribution of Holocene depocenters; deformed strata within a small pull-apart basin record punctuated growth faulting associated with at least three Holocene surface ruptures. An upper (shallower) landslide scarp, a buried sedimentary mound, and a deeper scarp have been right-laterally offset across the PVF by 55 ± 5, 52 ± 4 , and 39 ± 8 m, respectively. The ages of the upper scarp and buried mound are approximately 31 ka; the age of the deeper scarp is bracketed to 17–24 ka. These three piercing points bracket the late Pleistocene to present slip rate to 1.3–2.8 mm/yr and provide a best estimate of 1.6–1.9 mm/yr. The deformation observed along the PVF is characteristic of strike-slip faulting and accounts for 20–30% of the total right-lateral slip budget accommodated offshore Southern California.

California↗

USGS2015

Monitoring changes in seismic velocity related to an ongoing rapid inflation event at Okmok volcano, Alaska

Okmok is one of the most active volcanoes in the Aleutian Arc. In an effort to improve our ability to detect precursory activity leading to eruption at Okmok, we monitor a recent, and possibly ongoing, GPS-inferred rapid inflation event at the volcano using ambient noise interferometry (ANI). Applying this method, we identify changes in seismic velocity outside of Okmok’s caldera, which are related to the hydrologic cycle. Within the caldera, we observe decreases in seismic velocity that are associated with the GPS-inferred rapid inflation event. We also determine temporal changes in waveform decorrelation and show a continual increase in decorrelation rate over the time associated with the rapid inflation event. Themagnitude of relative velocity decreases and decorrelation rate increases are comparable to previous studies at Piton de la Fournaise that associate such changes with increased production of volatiles and/ormagmatic intrusion within the magma reservoir and associated opening of fractures and/or fissures. Notably, the largest decrease in relative velocity occurs along the intrastation path passing nearest to the center of the caldera. This observation, along with equal amplitude relative velocity decreases revealed via analysis of intracaldera autocorrelations, suggests that the inflation sourcemay be located approximately within the center of the caldera and represent recharge of shallow magma storage in this location. Importantly, there is a relative absence of seismicity associated with this and previous rapid inflation events at Okmok. Thus, these ANI results are the first seismic evidence of such rapid inflation at the volcano.

Alaska↗

USGS2015

Electrical properties of methane hydrate + sediment mixtures

Knowledge of the electrical properties of multicomponent systems with gas hydrate, sediments, and pore water is needed to help relate electromagnetic (EM) measurements to specific gas hydrate concentration and distribution patterns in nature. Toward this goal, we built a pressure cell capable of measuring in situ electrical properties of multicomponent systems such that the effects of individual components and mixing relations can be assessed. We first established the temperature-dependent electrical conductivity ( σ ) of pure, single-phase methane hydrate to be ~5 orders of magnitude lower than seawater, a substantial contrast that can help differentiate hydrate deposits from significantly more conductive water-saturated sediments in EM field surveys. Here we report σ measurements of two-component systems in which methane hydrate is mixed with variable amounts of quartz sand or glass beads. Sand by itself has low σ but is found to increase the overall σ of mixtures with well-connected methane hydrate. Alternatively, the overall σ decreases when sand concentrations are high enough to cause gas hydrate to be poorly connected, indicating that hydrate grains provide the primary conduction path. Our measurements suggest that impurities from sand induce chemical interactions and/or doping effects that result in higher electrical conductivity with lower temperature dependence. These results can be used in the modeling of massive or two-phase gas-hydrate-bearing systems devoid of conductive pore water. Further experiments that include a free water phase are the necessary next steps toward developing complex models relevant to most natural systems.

Journal of Geophysical Research↗

USGS2015

Holocene geologic slip rate for the Banning strand of the southern San Andreas Fault, southern California

Northwest directed slip from the southern San Andreas Fault is transferred to the Mission Creek, Banning, and Garnet Hill fault strands in the northwestern Coachella Valley. How slip is partitioned between these three faults is critical to southern California seismic hazard estimates but is poorly understood. In this paper, we report the first slip rate measured for the Banning fault strand. We constrain the depositional age of an alluvial fan offset 25 ± 5 m from its source by the Banning strand to between 5.1 ± 0.4 ka (95% confidence interval (CI)) and 6.4 + 3.7/−2.1 ka (95% CI) using U-series dating of pedogenic carbonate clast coatings and 10 Be cosmogenic nuclide exposure dating of surface clasts. We calculate a Holocene geologic slip rate for the Banning strand of 3.9 + 2.3/−1.6 mm/yr (median, 95% CI) to 4.9 + 1.0/−0.9 mm/yr (median, 95% CI). This rate represents only 25–35% of the total slip accommodated by this section of the southern San Andreas Fault, suggesting a model in which slip is less concentrated on the Banning strand than previously thought. In rejecting the possibility that the Banning strand is the dominant structure, our results highlight an even greater need for slip rate and paleoseismic measurements along faults in the northwestern Coachella Valley in order to test the validity of current earthquake hazard models. In addition, our comparison of ages measured with U-series and 10 Be exposure dating demonstrates the importance of using multiple geochronometers when estimating the depositional age of alluvial landforms.

Journal of Geophysical Research B: Solid Earth↗

USGS2015

Crustal deformation in the New Madrid seismic zone and the role of postseismic processes

Global Navigation Satellite System data across the New Madrid seismic zone (NMSZ) in the central United States over the period from 2000 through 2014 are analyzed and modeled with several deformation mechanisms including the following: (1) creep on subsurface dislocations, (2) postseismic frictional afterslip and viscoelastic relaxation from the 1811–1812 and 1450 earthquakes in the NMSZ, and (3) regional strain. In agreement with previous studies, a dislocation creeping at about 4 mm/yr between 12 and 20 km depth along the downdip extension of the Reelfoot fault reproduces the observations well. We find that a dynamic model of postseismic frictional afterslip from the 1450 and February 1812 Reelfoot fault events can explain this creep. Kinematic and dynamic models involving the Cottonwood Grove fault provide minimal predictive power. This is likely due to the smaller size of the December 1811 event on the Cottonwood Grove fault and a distribution of stations better suited to constrain localized strain across the Reelfoot fault. Regional compressive strain across the NMSZ is found to be less than 3 × 10 −9 /yr. If much of the present-day surface deformation results from afterslip, it is likely that many of the earthquakes we see today in the NMSZ are aftershocks from the 1811–1812 New Madrid earthquakes. Despite this conclusion, our results are consistent with observations and models of intraplate earthquake clustering. Given this and the recent paleoseismic history of the region, we suggest that seismic hazard is likely to remain significant.

Journal of Geophysical Research B: Solid Earth↗

USGS2015

Ionospheric current source modeling and global geomagnetic induction using ground geomagnetic observatory data

Long-period global-scale electromagnetic induction studies of deep Earth conductivity are based almost exclusively on magnetovariational methods and require accurate models of external source spatial structure. We describe approaches to inverting for both the external sources and three-dimensional (3-D) conductivity variations and apply these methods to long-period ( T ≥1.2 days) geomagnetic observatory data. Our scheme involves three steps: (1) Observatory data from 60 years (only partly overlapping and with many large gaps) are reduced and merged into dominant spatial modes using a scheme based on frequency domain principal components. (2) Resulting modes are inverted for corresponding external source spatial structure, using a simplified conductivity model with radial variations overlain by a two-dimensional thin sheet. The source inversion is regularized using a physically based source covariance, generated through superposition of correlated tilted zonal (quasi-dipole) current loops, representing ionospheric source complexity smoothed by Earth rotation. Free parameters in the source covariance model are tuned by a leave-one-out cross-validation scheme. (3) The estimated data modes are inverted for 3-D Earth conductivity, assuming the source excitation estimated in step 2. Together, these developments constitute key components in a practical scheme for simultaneous inversion of the catalogue of historical and modern observatory data for external source spatial structure and 3-D Earth conductivity.

Journal of Geophysical Research↗

USGS2015

Imaging the magmatic system of Mono Basin, California with magnetotellurics in three--dimensions

A three–dimensional (3D) electrical resistivity model of Mono Basin in eastern California unveils a complex subsurface filled with zones of partial melt, fluid–filled fracture networks, cold plutons, and regional faults. In 2013, 62 broadband magnetotelluric (MT) stations were collected in an array around southeastern Mono Basin from which a 3D electrical resistivity model was created with a resolvable depth of 35 km. Multiple robust electrical resistivity features were found that correlate with existing geophysical observations. The most robust features are two 300 ± 50 km3 near-vertical conductive bodies (3–10 Ω·m) that underlie the southeast and north-eastern margin of Mono Craters below 10 km depth. These features are interpreted as magmatic crystal–melt mush zones of 15 ± 5% interstitial melt surrounded by hydrothermal fluids and are likely sources for Holocene eruptions. Two conductive east–dipping structures appear to connect each magma source region to the surface. A conductive arc–like structure (< 0.9 Ω·m) links the northernmost mush column at 10 km depth to just below vents near Panum Crater, where the high conductivity suggests the presence of hydrothermal fluids. The connection from the southernmost mush column at 10 km depth to below South Coulée is less obvious with higher resistivity (200 Ω·m) suggestive of a cooled connection. A third, less constrained conductive feature (4–10 Ω·m) 15 km deep extending to 35 km is located west of Mono Craters near the eastern front of the Sierra Nevada escarpment, and is coincident with a zone of sporadic, long–period earthquakes that are characteristic of a fluid-filled (magmatic or metamorphic) fracture network. A resistive feature (103–105 Ω·m) located under Aeolian Buttes contains a deep root down to 25 km. The eastern edge of this resistor appears to structurally control the arcuate shape of Mono Craters. These observations have been combined to form a new conceptual model of the magmatic system beneath Mono Craters to a depth of 30 km.

California↗

USGS2015

Changes in seismic velocity during the first 14 months of the 2004–2008 eruption of Mount St. Helens, Washington

Mount St. Helens began erupting in late 2004 following an 18 year quiescence. Swarms of repeating earthquakes accompanied the extrusion of a mostly solid dacite dome over the next 4 years. In some cases the waveforms from these earthquakes evolved slowly, likely reflecting changes in the properties of the volcano that affect seismic wave propagation. We use coda-wave interferometry to quantify small changes in seismic velocity structure (usually <1%) between two similar earthquakes and employed waveforms from several hundred families of repeating earthquakes together to create a continuous function of velocity change observed at permanent stations operated within 20 km of the volcano. The high rate of earthquakes allowed tracking of velocity changes on an hourly time scale. Changes in velocity were largest near the newly extruding dome and likely related to shallow deformation as magma first worked its way to the surface. We found strong correlation between velocity changes and the inverse of real-time seismic amplitude measurements during the first 3 weeks of activity, suggesting that fluctuations of pressure in the shallow subsurface may have driven both seismicity and velocity changes. Velocity changes during the remainder of the eruption likely result from a complex interplay of multiple effects and are not well explained by any single factor alone, highlighting the need for complementary geophysical data when interpreting velocity changes.

Journal of Geophysical Research B: Solid Earth↗

USGS2016

Effective stress, friction and deep crustal faulting

Studies of crustal faulting and rock friction invariably assume the effective normal stress that determines fault shear resistance during frictional sliding is the applied normal stress minus the pore pressure. Here we propose an expression for the effective stress coefficient α f at temperatures and stresses near the brittle-ductile transition (BDT) that depends on the percentage of solid-solid contact area across the fault. α f varies with depth and is only near 1 when the yield strength of asperity contacts greatly exceeds the applied normal stress. For a vertical strike-slip quartz fault zone at hydrostatic pore pressure and assuming 1 mm and 1 km shear zone widths for friction and ductile shear, respectively, the BDT is at ~13 km. α f near 1 is restricted to depths where the shear zone is narrow. Below the BDT α f = 0 is due to a dramatically decreased strain rate. Under these circumstances friction cannot be reactivated below the BDT by increasing the pore pressure alone and requires localization. If pore pressure increases and the fault localizes back to 1 mm, then brittle behavior can occur to a depth of around 35 km. The interdependencies among effective stress, contact-scale strain rate, and pore pressure allow estimates of the conditions necessary for deep low-frequency seismicity seen on the San Andreas near Parkfield and in some subduction zones. Among the implications are that shear in the region separating shallow earthquakes and deep low-frequency seismicity is distributed and that the deeper zone involves both elevated pore fluid pressure and localization.

Journal of Geophysical Research B: Solid Earth↗

USGS2015

Resolving bathymetry from airborne gravity along Greenland fjords

Recent glacier mass loss in Greenland has been attributed to encroaching warming waters, but knowledge of fjord bathymetry is required to investigate this mechanism. The bathymetry in many Greenland fjords is unmapped and difficult to measure. From 2010 to 2012, National Aeronautics and Space Administration's Operation IceBridge collected a unique set of airborne gravity, magnetic, radar, and lidar data along the major outlet glaciers and fjords in Greenland. We applied a consistent technique using the IceBridge gravity data to create 90 bathymetric profiles along 54 Greenland fjords. We also used this technique to recover subice topography where warm or crevassed ice prevents the radar system from imaging the bed. Here we discuss our methodology, basic assumptions and error analysis. We present the new bathymetry data and discuss observations in six major regions of Greenland covered by IceBridge. The gravity models provide a total of 1950 line kilometers of bathymetry, 875 line kilometers of subice topography, and 12 new grounding line depths.

Journal of Geophysical Research B: Solid Earth↗

USGS2015

Pore-pressure sensitivities to dynamic strains: observations in active tectonic regions

Triggered seismicity arising from dynamic stresses is often explained by the Mohr-Coulomb failure criterion, where elevated pore pressures reduce the effective strength of faults in fluid-saturated rock. The seismic response of a fluid-rock system naturally depends on its hydro-mechanical properties, but accurately assessing how pore-fluid pressure responds to applied stress over large scales in situ remains a challenging task; hence, spatial variations in response are not well understood, especially around active faults. Here I analyze previously unutilized records of dynamic strain and pore-pressure from regional and teleseismic earthquakes at Plate Boundary Observatory (PBO) stations from 2006 through 2012 to investigate variations in response along the Pacific/North American tectonic plate boundary. I find robust scaling-response coefficients between excess pore pressure and dynamic strain at each station that are spatially correlated: around the San Andreas and San Jacinto fault systems, the response is lowest in regions of the crust undergoing the highest rates of secular shear strain. PBO stations in the Parkfield instrument cluster are at comparable distances to the San Andreas fault (SAF), and spatial variations there follow patterns in dextral creep rates along the fault, with the highest response in the actively creeping section, which is consistent with a narrowing zone of strain accumulation seen in geodetic velocity profiles. At stations in the San Juan Bautista (SJB) and Anza instrument clusters, the response depends non-linearly on the inverse fault-perpendicular distance, with the response decreasing towards the fault; the SJB cluster is at the northern transition from creeping-to-locked behavior along the SAF, where creep rates are at moderate to low levels, and the Anza cluster is around the San Jacinto fault, where to date there have been no statistically significant creep rates observed at the surface. These results suggest that the strength of the pore pressure response in fluid-saturated rock near active faults is controlled by shear strain accumulation associated with tectonic loading, which implies a strong feedback between fault strength and permeability: dynamic triggering susceptibilities may vary in space and also in time.

Journal of Geophysical Research B: Solid Earth↗