USGS Science⌕ Search

SEARCH · USGS Science

Results for “The Earth Observer”

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 883 records · Page 49Linked to original sources
USGS2018

The Polar WRF downscaled historical and projected twenty-first century climate for the coast and foothills of Arctic Alaska

Climate change is most pronounced in the northern high latitude region. Yet, climate observations are unable to fully capture regional-scale dynamics due to the sparse weather station coverage, which limits our ability to make reliable climate-based assessments. A set of simulated data products was therefore developed for the North Slope of Alaska through a dynamical downscaling approach. The polar-optimized Weather Research and Forecast (Polar WRF) model was forced by three sources: The ERA-interim reanalysis data (for 1979–2014), the Community Earth System Model 1.0 (CESM1.0) historical simulation (for 1950–2005), and the CESM1.0 projected (for 2006–2100) simulations in two Representative Concentration Pathways (RCP4.5 and RCP8.5) scenarios. Climatic variables were produced in a 10-km grid spacing and a 3-h interval. The ERA-interim forced WRF (ERA-WRF) proves the value of dynamical downscaling, which yields more realistic topographical-induced precipitation and air temperature, as well as corrects underestimations in observed precipitation. In summary, dry and cold biases to the north of the Brooks Range are presented in ERA-WRF, while CESM forced WRF (CESM-WRF) holds wet and warm biases in its historical period. A linear scaling method allowed for an adjustment of the biases, while keeping the majority of the variability and extreme values of modeled precipitation and air temperature. CESM-WRF under RCP 4.5 scenario projects smaller increase in precipitation and air temperature than observed in the historical CESM-WRF product, while the CESM-WRF under RCP 8.5 scenario shows larger changes. The fine spatial and temporal resolution, long temporal coverage, and multi-scenario projections jointly make the dataset appropriate to address a myriad of physical and biological changes occurring on the North Slope of Alaska.

Alaska↗
USGS2014

Residual shear strength variability as a primary control on movement of landslides reactivated by earthquake-induced ground motion: Implications for coastal Oregon, U.S.

Most large seismogenic landslides are reactivations of preexisting landslides with basal shear zones in the residual strength condition. Residual shear strength often varies during rapid displacement, but the response of residual shear zones to seismic loading is largely unknown. We used a ring shear apparatus to perform simulated seismic loading tests, constant displacement rate tests, and tests during which shear stress was gradually varied on specimens from two landslides to improve understanding of coseismic landslide reactivation and to identify shear strength models valid for slow gravitational failure through rapid coseismic failure. The landslides we studied represent many along the Oregon, U.S., coast. Seismic loading tests resulted in (1) catastrophic failure involving unbounded displacement when stresses represented those for the existing landslides and (2) limited to unbounded displacement when stresses represented those for hypothetical dormant landslides, suggesting that coseismic landslide reactivation may be significant during future great earthquakes occurring near the Oregon Coast. Constant displacement rate tests indicated that shear strength decreased exponentially during the first few decimeters of displacement but increased logarithmically with increasing displacement rate when sheared at 0.001 cm s −1 or greater. Dynamic shear resistance estimated from shear strength models correlated well with stresses observed during seismic loading tests, indicating that displacement rate and amount primarily controlled failure characteristics. We developed a stress-based approach to estimate coseismic landslide displacement that utilizes the variable shear strength model. The approach produced results that compared favorably to observations made during seismic loading tests, indicating its utility for application to landslides.

Oregon↗
USGS1990

Textural development of clayey and quartzofeldspathic fault gouges relative to their sliding behavior

Many of the secondary fault structures developed during triaxial friction experiments have been generally correlated with the structures of natural fault zones. Therefore, any physical differences that can be found between laboratory samples that slide stably and those that show stick-slip motion may help to identify the cause of earthquakes. We have examined petrographically the run products of many triaxial friction experiments using clayey and quartzofeldspathic gouges, which comprise the principal types of natural fault gouge material. The examined samples were tested under a wide range of temperature, confining and fluid pressure, and velocity conditions. The clayey and quartzofeldspathic gouges show some textural differences, owing to their different mineral contents and grain sizes and shapes. In the clayey gouges, for example, a clay mineral fabric and kink band sets are commonly developed, whereas in the quartzofeldspathic gouges fracturing and crushing of the predominately quartz and feldspar grains are important processes. For both types of gouge, however, and whatever the pressure-temperature-velocity conditions of the experiments, the transition from stable sliding to stick-slip motion is correlated with: (i) a change from pervasive deformation of the gouge layer to localized slip in subsidiary shears; and (ii) an increase in the angle betweem the shears that crosscut the gouge layer (Riedel shears) and ones that form along the gouge-rock cylinder boundaries (boundary shears). This suggests that the localization of shear within a fault zone combined with relatively high Riedel-shear angles are somehow connected with earthquakes. Secondary fracture sets similar to Riedel shears have been identified at various scales in major strike-slip faults such as the San Andreas of the western United States (Wallace, 1973) and the Luhuo and Fuyun earthquake faults of China (Deng and Zhang, 1984; Deng et al. , 1986). The San Andreas also contains locked and creeping sections that correspond to the stick-slip and stably sliding experimental samples, respectively. We plan to study the physical structure of the San Andreas fault, to see if the experimentally observed differences related to sliding behavior can also be distinguished in the field.

Physics and Chemistry of the Earth↗
USGS2011

Process-based modeling of tsunami inundation and sediment transport

The infrequent and unpredictable nature of tsunamis precludes the use of field experiments to measure the hydrodynamic and sediment transport processes that occur. Instead, these processes are often approximated from laboratory, numerical, and theoretical studies or inferred from observations of the resultant sediment deposits. Here Delft3D, a three-dimensional numerical model, is used to simulate the inundation and sediment transport of a tsunami similar in magnitude to the 26 December 2004 Indian Ocean tsunami over one measured and three idealized morphologies. The model is first shown to match well the observations taken at Kuala Meurisi, Sumatra, and then used to examine in detail the processes that occur during the tsunami. The model predicts that at a given cross-shore location the onshore flow accelerates rapidly to a maximum as the wavefront passes, and then gradually decelerates before reversing direction and flowing offshore. The onshore flow does not tend to zero everywhere at maximum inundation, but instead flow reversal occurs near the shoreline even as the wavefront continues to inundate landward. While some sediment is eroded by the passing wavefront, the suspension of sandy sediment is dominated by the long-duration, high-velocity backwash that occurs along the beach face and offshore of the shoreline. Some of the sediment suspended during backwash is advected shoreward by the subsequent wave, creating large spatial gradients in the suspended sediment concentrations, which may not be in equilibrium with the local hydrodynamics. The inundation and transport of sediment during a tsunami can be affected by complexities in the morphological profile and interactions between multiple waves, and many of the hydrodynamic and sediment transport processes predicted here are similar to analogous processes previously observed in the swash zone. Copyright 2011 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗
USGS1998

Poroelastic rebound along the Landers 1992 earthquake surface rupture

Maps of surface displacement following the 1992 Landers, California, earthquake, generated by interferometric processing of ERS-1 synthetic aperture radar (SAR) images, reveal effects of various postseismic deformation processes along the 1992 surface rupture. The large-scale pattern of the postseismic displacement field includes large lobes, mostly visible on the west side of the fault, comparable in shape with the lobes observed in the coseismic displacement field. This pattern and the steep displacement gradient observed near the Emerson-Camp Rock fault cannot be simply explained by afterslip on deep sections of the 1992 rupture. Models show that horizontal slip occurring on a buried dislocation in a Poisson's material produces a characteristic quadripole pattern in the surface displacement field with several centimeters of vertical motion at distances of 10–20 km from the fault, yet this pattern is not observed in the postseismic interferograms. As previously proposed to explain local strain in the fault step overs [ Peltzer et al ., 1996b], we argue that poroelastic rebound caused by pore fluid flow may also occur over greater distances from the fault, compensating the vertical ground shift produced by fault afterslip. Such a rebound is explained by the gradual change of the crustal rocks' Poisson's ratio value from undrained (coseismic) to drained (postseismic) conditions as pore pressure gradients produced by the earthquake dissipate. Using the Poisson's ratio values of 0.27 and 0.31 for the drained and undrained crustal rocks, respectively, elastic dislocation models show that the combined contributions of afterslip on deep sections of the fault and poroelastic rebound can account for the range change observed in the SAR data and the horizontal displacement measured at Global Positioning System (GPS) sites along a 60-km-long transect across the Emerson fault [ Savage and Svarc , 1997]. Using a detailed surface slip distribution on the Homestead Valley, Kickapoo, and Johnson Valley faults, we modeled the poroelastic rebound in the Homestead Valley pull apart. A Poisson's ratio value of 0.35 for the undrained gouge rocks in the fault zone is required to account for the observed surface uplift in the 3.5 years following the earthquake. This large value implies a seismic velocity ratio V p /V s of 2.1, consistent with the observed low V s values of fault zone guided waves at shallow depth [ Li et al ., 1997]. The SAR data also reveal postseismic creep along shallow patches of the Eureka Peak and Burnt Mountain faults with a characteristic decay time of 0.8 years. Coseismic, dilatant hardening (locking process) followed by post-seismic, pore pressure controlled fault creep provide a plausible mechanism to account for the decay time of the observed slip rate along this section of the fault.

Journal of Geophysical Research B: Solid Earth↗
USGS2023

Steady-state forms of channel profiles shaped by debris flow and fluvial processes

Debris flows regularly traverse bedrock channels that dissect steep landscapes, but our understanding of bedrock erosion by debris flows and their impact on steepland morphology is still rudimentary. Quantitative models of steep bedrock channel networks are based on geomorphic transport laws designed to represent erosion by water-dominated flows. To quantify the impact of debris flow erosion on steep channel network form, it is first necessary to develop methods to estimate spatial variations in bulk debris flow properties (e.g., flow depth, velocity) throughout the channel network that can be integrated into landscape evolution models. Here, we propose and evaluate two methods to estimate spatial variations in bulk debris flow properties along the length of a channel profile. We incorporate both methods into a model designed to simulate the evolution of longitudinal channel profiles that evolve in response to debris flow and fluvial processes. To explore this model framework, we propose a general family of debris flow erosion laws where erosion rate is a function of debris flow depth and channel slope. Model results indicate that erosion by debris flows can explain the occurrence of a scaling break in the slope–area curve at low-drainage areas and that upper-network channel morphology may be useful for inferring catchment-averaged erosion rates in quasi-steady landscapes. Validating specific forms of a debris flow incision law, however, would require more detailed model–data comparisons in specific landscapes where input parameters and channel morphometry can be better constrained. Results improve our ability to interpret topographic signals within steep channel networks and identify observational targets critical for constraining a debris flow incision law.

Earth Surface Dynamics↗
USGS2015

Crustal stress and structure at Kīlauea Volcano inferred from seismic anisotropy

Seismic anisotropy, measured through shear wave splitting (SWS) analysis, can be indicative of the state of stress in Earth's crust. Changes in SWS at Kīlauea Volcano, Hawai‘i, associated with the onset of summit eruptive activity in 2008 hint at the potential of the technique for tracking volcanic activity. To use SWS observations as a monitoring tool, however, it is important to understand the cause of seismic anisotropy at the volcano throughout the eruptive cycle. To address this need, we analyzed SWS results from across Kīlauea in combination with macroscopic surface structures (mapped fractures, faults, and fissures) and stress orientations inferred from fault plane solutions. Seismic anisotropy seems to be due to pervasive aligned structures in most regions of the volcano. The upper East and Southwest Rift Zones, however, show a bimodality in stress and SWS, suggesting a stress discontinuity with depth, perhaps related to magma conduits that trend obliquely to the dominant structure. Other areas in and around Kīlauea Caldera display principal stresses of similar magnitudes, indicating that small stress perturbations can rotate the maximum horizontal compressive stress direction by up to 90 ° . In these locations, static structures generally control SWS, but dynamic conditions due to magmatic activity can override the structural control. Monitoring of SWS may therefore provide important signs of impending volcanism.

Hawaii↗
USGS2026

Modeling the seasonality of wind-driven hydrocarbon waves in Titan’s polar lakes

Titan, the only body in the solar system aside from Earth with standing liquids on its surface, has polar hydrocarbon lakes and seas. As Titan’s atmosphere generates light winds, there should be waves on the surface of these lakes and seas, yet, direct wave observations are scant. We introduce and use PlanetWaves , an open source 4D spectral wave model, to study Titan’s waves and create seasonal maps of wave shape and propagation on Ontario Lacus and Ligeia Mare. Titan’s modeled waves grow up to 30 times larger than terrestrial waves for the same wind speed, are seasonally present and are largest in the spring and summer when winds are strongest. Average daily winds almost never exceed the wave generation threshold of 0.5–0.7 m/s. Average storm winds (∼1.5 m/s) generate waves 15–48 cm in height with a period ranging 6–10.5 s while maximum storm winds (∼4 m/s) generate waves 2.7–3.2 m in height with a period up to 32 s. Titan’s waves become fetch-independent at ∼40 km for average storm winds occurring ∼1% of a Titan year and ∼100 kilometers for maximum storm winds occurring 2-3 times per Titan decade. On Ontario Lacus, storm winds blow nearly parallel to the eastern shore, potentially driving wave modification of the smooth eastern shoreline. On Ligeia Mare, waves rarely propagate toward a hypothesized wave modified shoreline suggesting that another process, such as tectonics, may contribute to a straight shoreline morphology.

JGR Planets↗
USGS1997

Correlated errors in geodetic time series: Implications for time-dependent deformation

Analysis of frequent trilateration observations from the two-color electronic distance measuring networks in California demonstrate that the noise power spectra are dominated by white noise at higher frequencies and power law behavior at lower frequencies. In contrast, Earth scientists typically have assumed that only white noise is present in a geodetic time series, since a combination of infrequent measurements and low precision usually preclude identifying the time-correlated signature in such data. After removing a linear trend from the two-color data, it becomes evident that there are primarily two recognizable types of time-correlated noise present in the residuals. The first type is a seasonal variation in displacement which is probably a result of measuring to shallow surface monuments installed in clayey soil which responds to seasonally occurring rainfall; this noise is significant only for a small fraction of the sites analyzed. The second type of correlated noise becomes evident only after spectral analysis of line length changes and shows a functional relation at long periods between power and frequency of 1/ƒ α , where ƒ is frequency and α≈2. With α=2, this type of correlated noise is termed random-walk noise, and its source is mainly thought to be small random motions of geodetic monuments with respect to the Earth's crust, though other sources are possible. Because the line length changes in the two-color networks are measured at irregular intervals, power spectral techniques cannot reliably estimate the level of 1/ƒ α noise. Rather, we also use here a maximum likelihood estimation technique which assumes that there are only two sources of noise in the residual time series (white noise and random-walk noise) and estimates the amount of each. From this analysis we find that the random-walk noise level averages about 1.3 mm/√yr and that our estimates of the white noise component confirm theoretical limitations of the measurement technique. In addition, the seasonal noise can be as large as 3 mm in amplitude but typically is less than 0.5 mm. Because of the presence of random-walk noise in these time series, modeling and interpretation of the geodetic data must account for this source of error. By way of example we show that estimating the time-varying strain tensor (a form of spatial averaging) from geodetic data having both random-walk and white noise error components results in seemingly significant variations in the rate of strain accumulation; spatial averaging does reduce the size of both noise components but not their relative influence on the resulting strain accumulation model.

Journal of Geophysical Research B: Solid Earth↗
USGS2007

Analysis of continuous GPS measurements from southern Victoria Land, Antarctica

Several years of continuous data have been collected at remote bedrock Global Positioning System (GPS) sites in southern Victoria Land, Antarctica. Annual to sub-annual variations are observed in the position time-series. An atmospheric pressure loading (APL) effect is calculated from pressure field anomalies supplied by the European Centre for Medium-Range Weather Forecasts (ECMWF) model loading an elastic Earth model. The predicted APL signal has a moderate correlation with the vertical position time-series at McMurdo, Ross Island (International Global Navigation Satellite System Service (IGS) station MCM4), produced using a global solution. In contrast, a local solution in which MCM4 is the fiducial site generates a vertical time series for a remote site in Victoria Land (Cape Roberts, ROB4) which exhibits a low, inverse correlation with the predicted atmospheric pressure loading signal. If, in the future, known and well modeled geophysical loads can be separated from the time-series, then local hydrological loading, of interest for glaciological and climate applications, can potentially be extracted from the GPS time-series.

Open-File Report↗
USGS2016

Toward more realistic projections of soil carbon dynamics by Earth system models

Soil carbon (C) is a critical component of Earth system models (ESMs), and its diverse representations are a major source of the large spread across models in the terrestrial C sink from the third to fifth assessment reports of the Intergovernmental Panel on Climate Change (IPCC). Improving soil C projections is of a high priority for Earth system modeling in the future IPCC and other assessments. To achieve this goal, we suggest that (1) model structures should reflect real-world processes, (2) parameters should be calibrated to match model outputs with observations, and (3) external forcing variables should accurately prescribe the environmental conditions that soils experience. First, most soil C cycle models simulate C input from litter production and C release through decomposition. The latter process has traditionally been represented by first-order decay functions, regulated primarily by temperature, moisture, litter quality, and soil texture. While this formulation well captures macroscopic soil organic C (SOC) dynamics, better understanding is needed of their underlying mechanisms as related to microbial processes, depth-dependent environmental controls, and other processes that strongly affect soil C dynamics. Second, incomplete use of observations in model parameterization is a major cause of bias in soil C projections from ESMs. Optimal parameter calibration with both pool- and flux-based data sets through data assimilation is among the highest priorities for near-term research to reduce biases among ESMs. Third, external variables are represented inconsistently among ESMs, leading to differences in modeled soil C dynamics. We recommend the implementation of traceability analyses to identify how external variables and model parameterizations influence SOC dynamics in different ESMs. Overall, projections of the terrestrial C sink can be substantially improved when reliable data sets are available to select the most representative model structure, constrain parameters, and prescribe forcing fields.

Global Biogeochemical Cycles↗
USGS2012

Differential melt scaling for oblique impacts on terrestrial planets

Analytical estimates of melt volumes produced by a given projectile and contained in a given impact crater are derived as a function of impact velocity, impact angle, planetary gravity, target and projectile densities, and specific internal energy of melting. Applications to impact events and impact craters on the Earth, Moon, and Mars are demonstrated and discussed. The most probable oblique impact (45°) produces ∼1.6 times less melt volume than a vertical impact, and ∼1.6 and 3.7 times more melt volume than impacts with 30° and 15° trajectories, respectively. The melt volume for a particular crater diameter increases with planetary gravity, so a crater on Earth should have more melt than similar-size craters on Mars and the Moon. The melt volume for a particular projectile diameter does not depend on gravity, but has a strong dependence on impact velocity, so the melt generated by a given projectile on the Moon is significantly larger than on Mars. Higher surface temperatures and geothermal gradients increase melt production, as do lower energies of melting. Collectively, the results imply thinner central melt sheets and a smaller proportion of melt particles in impact breccias on the Moon and Mars than on Earth. These effects are illustrated in a comparison of the Chicxulub crater on Earth, linked to the Cretaceous–Tertiary mass extinction, Gusev crater on Mars, where the Mars Exploration Rover Spirit landed, and Tsiolkovsky crater on the Moon. The results are comparable to those obtained from field and spacecraft observations, other analytical expressions, and hydrocode simulations.

Icarus↗
USGS2000

Acoustic properties of a crack containing magmatic or hydrothermal fluids

We estimate the acoustic properties of a crack containing magmatic or hydrothermal fluids to quantify the source properties of long-period (LP) events observed in volcanic areas assuming that a crack-like structure is the source of LP events. The tails of synthetic waveforms obtained from a model of a fluid-driven crack are analyzed by the Sompi method to determine the complex frequencies of one of the modes of crack resonance over a wide range of the model parameters α/ a and ρ f /ρ s , where αis the P wave velocity of the rock matrix, a is the sound speed of the fluid, and ρ f and p s are the densities of the fluid and rock matrix, respectively. The quality factor due to radiation loss ( Q r ) for the selected mode almost monotonically increases with increasing α/ a , while the dimensionless frequency ( v ) of the mode decreases with increasing α/ a and ρ f /ρ s . These results are used to estimate Q and v for a crack containing various types of fluids (gas-gas mixtures, liquid-gas mixtures, and dusty and misty gases) for values of a , ρ f , and quality factor due to intrinsic losses ( Q i ) appropriate for these types of fluids, in which Q is given by Q −l = Q r −l + Q i −1 . For a crack containing such fluids, we obtain Q ranging from almost unity to several hundred, which consistently explains the wide variety of quality factors measured in LP events observed at various volcanoes. We underscore the importance of dusty and misty gases containing small-size particles with radii around 1 μm to explain long-lasting oscillations with Q significantly larger than 100. Our results may provide a basis for the interpretation of spatial and temporal variations in the observed complex frequencies of LP events in terms of fluid compositions beneath volcanoes.

Journal of Geophysical Research B: Solid Earth↗
USGS2021

Improving the usability of Galileo and Voyager images of Jupiter’s moon, Europa

NASA's Voyager 1, Voyager 2, and Galileo spacecraft acquired hundreds of images of Jupiter's moon Europa. These images provide the only moderate- to high-resolution views of the moon's surface and are therefore a critical resource for scientific analysis and future mission planning. Unfortunately, uncertain knowledge of the spacecraft's position and pointing during image acquisition resulted in significant errors in the location of the images on the surface. The result is that adjacent images are poorly aligned, with some images displaced by more than 100 km from their correct location. These errors severely degrade the usability of the Voyager and Galileo imaging data sets. To improve the usability of these data sets, we used the U.S. Geological Survey Integrated Software for Imagers and Spectrometers to build a nearly global image tie-point network with more than 50,000 tie points and 135,000 image measurements on 481 Galileo and 221 Voyager images. A global least-squares bundle adjustment of our final Europa tie-point network calculated latitude, longitude, and radius values for each point by minimizing residuals globally, and resulted in root mean square (RMS) uncertainties of 246.6 m, 307.0 m, and 70.5 m in latitude, longitude, and radius, respectively. The total RMS uncertainty was 0.32 pixels. This work enables direct use of nearly the entire Galileo and Voyager image data sets for Europa. We are providing the community with updated NASA Navigation and Ancillary Information Facility Spacecraft, Planet, Instrument, C-matrix (pointing), and Events kernels, mosaics of Galileo images acquired during each observation sequence, and individual processed and projected level 2 images.

Earth and Space Science↗
USGS2012

Deformation from 1973 to 1987 in the epicentral area of the 1959 Hebgen Lake, Montana, earthquake (Ms = 7.5)

A 40‐km aperture trilateration network centered on the 1959 Hebgen Lake earthquake epicenter has been surveyed in 1973, 1974, 1976, 1978, 1981, 1984, and 1987. The deformation inferred from those surveys is described roughly by a uniaxial, 0.266 ± 0.014 μstrain/yr, N15°E ± 1°extension that is uniform in both time and space. That extension is orthogonal to the strike (N78°W ± 5°) of the 1959 rupture plane inferred from focal mechanism solutions. The absence of strain accumulation in the N75°W direction suggests that the source of deformation must be greatly elongated in that direction. Looked at in greater detail, the deformation is found to be concentrated in the northern part of the network. The zone of concentrated deformation coincides with a part of a well‐defined trend in seismicity that extends at least 100 km N75°W from the Sour Creek resurgent dome in the Yellowstone caldera. Extension perpendicular to that trend is indicated by focal mechanism solutions. Thus, the zone of rifting identified within the Hebgen Lake network apparently extends east southeast to the Yellowstone caldera. About 8 mm/yr extension is observed across a 20‐km width within the zone.

Montana, Wyoming↗
USGS2000

Land-subsidence and ground-water storage monitoring in the Tucson Active Management Area, Arizona

The Tucson Active Management Area (TAMA) comprises two basins--Tucson Basin and Avra Valley. The TAMA has been directed by Arizona ground-water law to attain an annual balance between groundwater withdrawals and recharge by the year 2025. This balance is defined by the statute as "safe yield." Current ground-water withdrawals exceed recharge, resulting in conditions of ground-water overdraft, which causes removal of water from ground-water storage and subsidence of the land surface. Depletion of storage and associated land subsidence will not be halted until all discharge from the system, both natural and human induced, is balanced by recharge. The amount of the ground-water overdraft has been difficult to estimate until recently because it could not be directly measured. Overdraft has been estimated using indirect water-budget methods that rely on uncertain estimates of recharge. As a result, the status of the ground-water budget could not be known with great certainty. Gravity methods offer a means to directly measure ground-water overdraft through measurement of changes in the gravitational field of the Earth that are caused by changes in the amount of water stored in the subsurface. Changes in vertical position also affect the measured gravity value and thus subsidence also must be monitored. The combination of periodic observations of gravity and vertical positions provide direct measures of changes in stored ground water and land subsidence.

Arizona↗
USGS2012

Short-period strain (0.1–105 s): Near-source strain field for an earthquake (ML 3.2) near San Juan Bautista, California

Measurements of dilational earth strain in the frequency band 25–10 −5 Hz have been made on a deep borehole strainmeter installed near the San Andreas fault. These data are used to determine seismic radiation fields during nuclear explosions, teleseisms, local earthquakes, and ground noise during seismically quiet times. Strains of less than 10 −10 on these instruments can be clearly resolved at short periods (< 10 s) and are recorded with wide dynamic range digital recorders. This permits measurement of the static and dynamic strain variations in the near field of local earthquakes. Noise spectra for earth strain referenced to 1 (strain) 2 /Hz show that strain resolution decreases at about 10 dB per decade of frequency from −150 dB at 10 −4 Hz to −223 dB at 10 Hz. Exact expressions are derived to relate the volumetric strain and displacement field for a homogeneous P wave in a general viscoelastic solid as observed on colocated dilatometers and seismometers. A rare near-field recording of strain and seismic velocity was obtained on May 26, 1984, from an earthquake ( M L 3.2) at a hypocentral distance of 3.2 km near the San Andreas fault at San Juan Bautista, California. While the data indicate no precursory strain release at the 5 × 10 −11 strain level, a coseismic strain release of 1.86 nanostrain was observed. This change in strain is consistent with that calculated from a simple dislocation model of the event. Ground displacement spectra, determined from the downhole strain data and instrument-corrected surface seismic data, suggest that source parameters estimated from surface recordings may be contaminated by amplification effects in near-surface low-velocity materials.

California↗
USGS2007

The role of fault zone drilling

The objective of fault-zone drilling projects is to directly study the physical and chemical processes that control deformation and earthquake generation within active fault zones. An enormous amount of field, laboratory, and theoretical work has been directed toward the mechanical and hydrological behavior of faults over the past several decades. Nonetheless, it is currently impossible to differentiate between – or even adequately constrain – the numerous conceptual models of active faults proposed over the years. For this reason, the Earth science community is left in the untenable position of having no generally accepted paradigm for the mechanical behavior of faults at depth. One of the primary causes for this dilemma is the difficulty of either directly observing or inferring physical properties and deformation mechanisms along faults at depth, as well as the need to observe directly key parameters such as the state of stress acting on faults at depth, pore fluid pressure (and its possible variation in space and time), and processes associated with earthquake nucleation and rupture. Today, we know very little about the composition of active faults at depth, their constitutive properties, the state of in situ stress or pore pressure within fault zones, the origin of fault-zone pore fluids, or the nature and significance of time-dependent fault-zone processes.

Book chapter↗