SEARCH · USGS Science

Results for “Geophysical Journal”

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 559 records · Page 31Linked to original sources

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

Source limitation of carbon gas emissions in high-elevation mountain streams and lakes

Inland waters are an important component of the global carbon cycle through transport, storage, and direct emissions of CO 2 and CH 4 to the atmosphere. Despite predictions of high physical gas exchange rates due to turbulent flows and ubiquitous supersaturation of CO 2 —and perhaps also CH 4 —patterns of gas emissions are essentially undocumented for high mountain ecosystems. Much like other headwater networks around the globe, we found that high-elevation streams in Rocky Mountain National Park, USA, were supersaturated with CO 2 during the growing season and were net sources to the atmosphere. CO 2 concentrations in lakes, on the other hand, tended to be less than atmospheric equilibrium during the open water season. CO 2 and CH 4 emissions from the aquatic conduit were relatively small compared to many parts of the globe. Irrespective of the physical template for high gas exchange (high k ), we found evidence of CO 2 source limitation to mountain streams during the growing season, which limits overall CO 2 emissions. Our results suggest a reduced importance of aquatic ecosystems for carbon cycling in high-elevation landscapes having limited soil development and high CO 2 consumption via mineral weathering.

Journal of Geophysical Research G: Biogeosciences↗

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↗

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↗

USGS2015

Low resistivity and permeability in actively deforming shear zones on the San Andreas Fault at SAFOD

The San Andreas Fault Observatory at Depth (SAFOD) scientific drillhole near Parkfield, California crosses the San Andreas Fault at a depth of 2.7 km. Downhole measurements and analysis of core retrieved from Phase 3 drilling reveal two narrow, actively deforming zones of smectite-clay gouge within a roughly 200 m-wide fault damage zone of sandstones, siltstones and mudstones. Here we report electrical resistivity and permeability measurements on core samples from all of these structural units at effective confining pressures up to 120 MPa. Electrical resistivity (~10 ohm-m) and permeability (10-21 to 10-22 m2) in the actively deforming zones were one to two orders of magnitude lower than the surrounding damage zone material, consistent with broader-scale observations from the downhole resistivity and seismic velocity logs. The higher porosity of the clay gouge, 2 to 8 times greater than that in the damage zone rocks, along with surface conduction were the principal factors contributing to the observed low resistivities. The high percentage of fine-grained clay in the deforming zones also greatly reduced permeability to values low enough to create a barrier to fluid flow across the fault. Together, resistivity and permeability data can be used to assess the hydrogeologic characteristics of the fault, key to understanding fault structure and strength. The low resistivities and strength measurements of the SAFOD core are consistent with observations of low resistivity clays that are often found in the principal slip zones of other active faults making resistivity logs a valuable tool for identifying these zones.

Journal of Geophysical Research↗

USGS2015

A rapid estimation of near field tsunami run-up

Many efforts have been made to quickly estimate the maximum run-up height of tsunamis associated with large earthquakes. This is a difficult task, because of the time it takes to construct a tsunami model using real time data from the source. It is possible to construct a database of potential seismic sources and their corresponding tsunami a priori.However, such models are generally based on uniform slip distributions and thus oversimplify the knowledge of the earthquake source. Here, we show how to predict tsunami run-up from any seismic source model using an analytic solution, that was specifically designed for subduction zones with a well defined geometry, i.e., Chile, Japan, Nicaragua, Alaska. The main idea of this work is to provide a tool for emergency response, trading off accuracy for speed. The solutions we present for large earthquakes appear promising. Here, run-up models are computed for: The 1992 Mw 7.7 Nicaragua Earthquake, the 2001 Mw 8.4 Perú Earthquake, the 2003Mw 8.3 Hokkaido Earthquake, the 2007 Mw 8.1 Perú Earthquake, the 2010 Mw 8.8 Maule Earthquake, the 2011 Mw 9.0 Tohoku Earthquake and the recent 2014 Mw 8.2 Iquique Earthquake. The maximum run-up estimations are consistent with measurements made inland after each event, with a peak of 9 m for Nicaragua, 8 m for Perú (2001), 32 m for Maule, 41 m for Tohoku, and 4.1 m for Iquique. Considering recent advances made in the analysis of real time GPS data and the ability to rapidly resolve the finiteness of a large earthquake close to existing GPS networks, it will be possible in the near future to perform these calculations within the first minutes after the occurrence of similar events. Thus, such calculations will provide faster run-up information than is available from existing uniform-slip seismic source databases or past events of pre-modeled seismic sources.

Journal of Geophysical Research↗

USGS2016

Vertical deformation associated with normal fault systems evolved over coseismic, postseismic, and multiseismic periods

Vertical deformation of extensional provinces varies significantly and in seemingly contradictory ways. Sparse but robust geodetic, seismic, and geologic observations in the Basin and Range province of the western United States indicate that immediately after an earthquake, vertical change primarily occurs as subsidence of the normal fault hanging wall. A few decades later, a ±100 km wide zone is symmetrically uplifted. The preserved topography of long-term rifting shows bent and tilted footwall flanks rising high above deep basins. We develop finite element models subjected to extensional and gravitational forces to study time-varying deformation associated with normal faulting. We replicate observations with a model that has a weak upper mantle overlain by a stronger lower crust and a breakable elastic upper crust. A 60° dipping normal fault cuts through the upper crust and extends through the lower crust to simulate an underlying shear zone. Stretching the model under gravity demonstrates that asymmetric slip via collapse of the hanging wall is a natural consequence of coseismic deformation. Focused flow in the upper mantle imposed by deformation of the lower crust localizes uplift under the footwall; the breakable upper crust is a necessary model feature to replicate footwall bending over the observed width ( < 10 km), which is predicted to take place within 1-2 decades after each large earthquake. Thus the best-preserved topographic signature of rifting is expected to occur early in the postseismic period. The relatively stronger lower crust in our models is necessary to replicate broader postseismic uplift that is observed geodetically in subsequent decades.

Journal of Geophysical Research↗

USGS2015

Identifying block structure in the Pacific Northwest, USA

We have identified block structure in the Pacific Northwest (west of 116°W between 38°N and 49°N) by clustering GPS stations so that the same Euler vector approximates the velocity of each station in a cluster. Given the total number k of clusters desired, the clustering procedure finds the best assignment of stations to clusters. Clustering is calculated for k= 2 to 14. In geographic space, cluster boundaries that remain relatively stable as k is increased are tentatively identified as block boundaries. That identification is reinforced if the cluster boundary coincides with a geologic feature. Boundaries identified in northern California and Nevada are the Central Nevada Seismic Belt, the west side of the Northern Walker Lane Belt, and the Bartlett Springs Fault. Three blocks cover all of Oregon and Washington. The principal block boundary there extends west-northwest along the Brothers Fault Zone, then north and northwest along the eastern boundary of Siletzia, the accreted oceanic basement of the forearc. East of this boundary is the Intermountain block, its eastern boundary undefined. A cluster boundary at Cape Blanco subdivides the forearc along the faulted southern margin of Siletzia. South of Cape Blanco the Klamath Mountains-Basin and Range block extends east to the Central Nevada Seismic Belt and south to the Sierra Nevada-Great Valley block. The Siletzia block north of Cape Blanco coincides almost exactly with the accreted Siletz terrane. The cluster boundary in the eastern Olympic Peninsula may mark permanent shortening of Siletzia against the Intermountain block.

California, Idaho, Nevada, Oregon, Washington↗

USGS2015

Borehole strainmeter measurements spanning the 2014, M w 6.0 South Napa Earthquake, California: The effect from instrument calibration

The 24 August 2014 M w 6.0 South Napa, California earthquake produced significant offsets on 12 borehole strainmeters in the San Francisco Bay area. These strainmeters are located between 24 and 80 km from the source and the observed offsets ranged up to 400 parts-per-billion (ppb), which exceeds their nominal precision by a factor of 100. However, the observed offsets of tidally calibrated strains differ by up to 130 ppb from predictions based on a moment tensor derived from seismic data. The large misfit can be attributed to a combination of poor instrument calibration and better modeling of the strain fit from the earthquake. Borehole strainmeters require in-situ calibration, which historically has been accomplished by comparing their measurements of Earth tides with the strain-tides predicted by a model. Although the borehole strainmeter accurately measure the deformation within the borehole, the long-wavelength strain signals from tides or other tectonic processes recorded in the borehole are modified by the presence of the borehole and the elastic properties of the grout and the instrument. Previous analyses of surface-mounted, strainmeter data and their relationship with the predicted tides suggest that tidal models could be in error by 30%. The poor fit of the borehole strainmeter data from this earthquake can be improved by simultaneously varying the components of the model tides up to 30% and making small adjustments to the point-source model of the earthquake, which reduces the RMS misfit from 130 ppb to 18 ppb. This suggests that relying on tidal models to calibrate borehole strainmeters significantly reduces their accuracy.

California↗