SEARCH · Search USGS
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.
Semidiurnal water‐level fluctuations of tidal period have been observed in wells completed in the Chickamauga limestone of Middle and Late Ordovician age in east Tennessee. The periodic oscillations of the water level are similar to fluctuations in artesian pressure that in other localities have been recognized and correlated with earth tides.
Coseismic stress changes have been the primary physical principle used to explain aftershocks and triggered earthquakes. However, this method does not adequately forecast earthquake rates and diverse rupture populations when subjected to formal testing. We show that earthquake forecasts can be impaired by assumptions made in physics-based models, such as the existence of hypothetical optimal faults, and regional scale invariability of the stress field. We compare calculations made under these assumptions along with different realizations of a new conceptual triggering model that features a complete assay of all possible ruptures. In this concept, there always exists a set of theoretical planes that has positive failure stress conditions under a combination of background and coseismic static stress change. In the Earth, all of these theoretical planes may not exist, and if they do, they may not be ready to fail. Thus the actual aftershock plane may not correspond to the plane with the maximum stress change value. This is consistent with observations that mainshocks commonly activate faults with exotic orientations and rakes. Our testing ground is the M=7.2, 2010 El Mayor-Cucapah earthquake sequence that activated multiple diverse fault populations across the USA-Mexico border in California and Baja California. We carry out a retrospective test involving 748 M≥3.0 triggered earthquakes that occurred during a 3-yr period after the mainshock. We find that a probabilistic expression of possible aftershock planes constrained by pre-mainshock rupture patterns is strongly favoured (89% of aftershocks consistent with static stress triggering) versus an optimal fault implementation (35% consistent). Results show that coseismic stress change magnitudes do not necessarily control earthquake triggering, instead we find that the summed background stress and coseismic stress change promotes diverse ruptures. Our model can thus explain earthquake triggering in regions where optimal plane mapping shows coseismic stress reduction.
The current (March 2008 to February 2009) summit eruptive activity at Kilauea Volcano is characterized by explosive degassing bursts accompanied by very long period (VLP) seismic signals. We model the source mechanisms of VLP signals in the 10–50 s band using data recorded for 15 bursts with a 10‐station broadband network deployed in the summit caldera. To determine the source centroid location and source mechanism, we minimize the residual error between data and synthetics calculated by the finite difference method for a point source embedded in a homogeneous medium that takes topography into account. The VLP signals associated with the bursts originate in a source region ∼1 km below the eastern perimeter of Halemaumau pit crater. The observed waveforms are well explained by the combination of a volumetric component and a vertical single force component. For the volumetric component, several source geometries are obtained which equally explain the observed waveforms. These geometries include (1) a pipe dipping 64° to the northeast; (2) two intersecting cracks including an east striking crack (dike) dipping 80° to the north, intersecting a north striking crack (another dike) dipping 65° to the east; (3) a pipe dipping 58° to the northeast, intersecting a crack dipping 48° to the west–southwest; and (4) a pipe dipping 57° to the northeast, intersecting a pipe dipping 58° to the west–southwest. Using the dual‐crack model as reference, the largest volume change obtained among the 15 bursts is ∼24,400 m 3 , and the maximum amplitude (peak to peak) of the force is ∼20 GN. Each burst is marked by a similar sequence of deflation and inflation, trailed by decaying oscillations of the volumetric source. The vertical force is initially upward, synchronous with source deflation, then downward, synchronous with source reinflation, followed by oscillations with polarity opposite to the volumetric oscillations. This combination of force and volume change is attributed to pressure and momentum changes induced during a fluid dynamic source mechanism involving the ascent, expansion, and burst of a large slug of gas within the upper ∼150 m of the magma conduit. As the slug expands upon approach to the surface and more liquid becomes wall supported by viscous shear forces, the pressure below the slug decreases, inducing conduit deflation and an upward force on the Earth. The final rapid slug expansion and burst stimulate VLP and LP oscillations of the conduit system, which slowly decay due to viscous dissipation and elastic radiation. Consideration of the fluid dynamic arguments leads us to prefer the dual‐crack VLP source model as it is the only candidate model capable of producing plausible values of length scales and pressure changes. The magnitudes of the vertical forces observed in the 15 bursts appear consistent with slug masses of 10 4 to 10 6 kg.
The algorithm used to transform velocity signals from beam coordinates to earth coordinates in an acoustic Doppler current profiler (ADCP) relies on the assumption that the currents are uniform over the horizontal distance separating the beams. This condition may be violated by (nonlinear) internal waves, which can have wavelengths as small as 100-200 m. In this case, the standard algorithm combines velocities measured at different phases of a wave and produces horizontal velocities that increasingly differ from true velocities with distance from the ADCP. Observations made in Massachusetts Bay show that currents measured with a bottom-mounted upward-looking ADCP during periods when short-wavelength internal waves are present differ significantly from currents measured by point current meters, except very close to the instrument. These periods are flagged with high error velocities by the standard ADCP algorithm. In this paper measurements from the four spatially diverging beams and the backscatter intensity signal are used to calculate the propagation direction and celerity of the internal waves. Once this information is known, a modified beam-to-earth transformation that combines appropriately lagged beam measurements can be used to obtain current estimates in earth coordinates that compare well with pointwise measurements. ?? 2005 American Meteorological Society.
The initial coda of teleseismic P-waves contains considerable information about the crust and upper mantle structure directly beneath a receiver. When this information can be recovered for a dense network of seismographs much can be learned about the structure of the earth. Data from the high quality broadband seismic stations of the SKIPPY and KIMBA projects along with permanent stations are used to investigate the upper crustal structure of Australia. A dataset of 65 shear-velocity models derived from receiver functions has enabled the sedimentary and upper crustal structure of Australia to be summarised. Regions of thick soft sediment show good agreement with topographical lows. A simple relation between upper-crustal velocity and magnetisation, as has been suggested by other investigators, has not been observed, but this may be due to the magnetic signal being muted by overlying sediments. A prominent mid-crustal discontinuity is apparent in the Tasman and New England mega-elements. This may represent a mid-crustal decollement that had structural control during accretion.
Over 3 500 measurements of surface heat-flux have been catalogued and analyzed to study the large-scale variations of terrestrial heat-flow. It was found that heat-flow values are correlated with major geologic provinces: higher averages and scattered values in active tectonic regions, and lower averages and more uniform values in stable areas. Analyzing the data in the light of new global tectonics shows that the variations of heat-flow are consistent with the hypotheses of sea-floor spreading and plate tectonics. The observed heat-flow across the mid-oceanic ridges can be accounted for by a simple model of a spreading sea floor.
Global-scale studies suggest that dryland ecosystems dominate an increasing trend in the magnitude and interannual variability of the land CO 2 sink. However, such model-based analyses are poorly constrained by measured CO 2 exchange in open shrublands, which is the most common global land cover type, covering ∼14% of Earth’s surface. Here we evaluate how the amount and seasonal timing of water availability regulate CO 2 exchange between shrublands and the atmosphere. We use eddy covariance data from six US sites across the three warm deserts of North America with observed ranges in annual precipitation of ∼100–400mm, annual temperatures of 13–18°C, and records of 2–8 years (33 site-years in total). The Chihuahuan, Sonoran and Mojave Deserts present gradients in both mean annual precipitation and its seasonal distribution between the wet-winter Mojave Desert and the wet-summer Chihuahuan Desert. We found that due to hydrologic losses during the wettest summers in the Sonoran and Chihuahuan Deserts, evapotranspiration (ET) was a better metric than precipitation of water available to drive dryland CO 2 exchange. In contrast with recent synthesis studies across diverse dryland biomes, we found that NEP could not be directly predicted from ET due to wintertime decoupling of the relationship between ecosystem respiration (R eco ) and gross ecosystem productivity (GEP). Ecosystem water use efficiency (WUE=GEP/ET) did not differ between winter and summer. Carbon use efficiency (CUE=NEP/GEP), however, was greater in winter because R eco returned a smaller fraction of carbon to the atmosphere (23% of GEP) than in summer (77%). Combining the water-carbon relations found here with historical precipitation since 1980, we estimate that lower average winter precipitation during the 21st century reduced the net carbon sink of the three deserts by an average of 6.8TgC yr 1 . Our results highlight that winter precipitation is critical to the annual carbon balance of these warm desert shrublands.
Detailed monitoring at Mount St. Helens since 1980 has enabled prediction of the intermittent eruptive activity (mostly dome growth) with unprecedented success. During 1981 and 1982, accelerating deformation of the crater floor around the vent (including radial cracks, thrust faults, and ground tilt) was the earliest indicator of impending activity. Numerical experiments using the finite element method to model the mechanics of the crater floor show that all deformation features can be explained by a uniform shear-stress boundary condition along the conduit wall. The magnitude of the shear stress required to match observed displacements (1–7 MPa) is inversely proportional to the conduit diameter (estimated to be 25–100 m). The most probable source of this shear stress is the flow of viscous magma up the conduit and into the lava dome. We propose a model in which the accelerating deformation, beginning as much as 4 weeks before extrusions, is caused by the increasing velocity of ascending magma in the conduit. This model is examined by using deformation data of the dome before four extrusions in 1981 and 1982 to estimate the volumetric flow rate through the conduit. This flow rate and an estimate of the effective viscosity of the magma enable calculation of an ascent velocity and an applied shear stress that, again, depend on the conduit diameter. The results of these calculations are consistent with the finite element experiments and show that the proposed model is feasible. In light of this model, events observed just before or near the time extrusions began, such as reversals of ground tilt direction from outward to inward and the sudden decrease in the number of shallow earthquakes, may indicate an abrupt decrease of shear stress in the conduit. This could be explained by a decrease in either the ascent velocity, or the effective viscosity of the magma ascending through the shallow conduit, or both, near the time of extrusion. Precursory deformation like that measured at Mount St. Helens should be observable at similar volcanoes elsewhere because it is caused by the fundamental process of magma ascent.
A combined aeromagnetic and radio echo ice-sounding survey made in 1978 in Antarctica over the Dufek layered mafic intrusion suggests a minimum area of the intrusion of about 50,000 square kilometers, making it comparable in size with the Bushveld Complex of Africa. Comparisons of the magnetic and subglacial topographic profiles illustrate the usefulness of this combination of methods in studying bedrock geology beneath ice-covered areas. Magnetic anomalies range in peak-to-trough amplitude from about 50 nanoteslas over the lowermost exposed portion of the section in the Dufek Massif to about 3600 nanoteslas over the uppermost part of the section in the Forrestal Range. Theoretical magnetic anomalies, computed from a model based on the subice topography fitted to the highest amplitude observed magnetic anomalies, required normal and reversed magnetizations ranging from 10 -3 to 10 -2 electromagnetic units per cubic centimeter. This result is interpreted as indicating that the Dufek intrusion cooled through the Curie isotherm during one or more reversals of the earth's magnetic field.
Inter-calibration of satellite instruments is critical for detection and quantification of changes in the Earth’s environment, weather forecasting, understanding climate processes, and monitoring climate and land cover change. These applications use data from many satellites; for the data to be inter-operable, the instruments must be cross-calibrated. To meet the stringent needs of such applications requires that instruments provide reliable, accurate, and consistent measurements over time. Robust techniques are required to ensure that observations from different instruments can be normalized to a common scale that the community agrees on. The long-term reliability of this process needs to be sustained in accordance with established reference standards and best practices. Furthermore, establishing physical meaning to the information through robust Système International d'unités (SI) traceable Calibration and Validation (Cal/Val) is essential to fully understand the parameters under observation. The processes of calibration, correction, stability monitoring, and quality assurance need to be underpinned and evidenced by comparison with “peer instruments” and, ideally, highly calibrated in-orbit reference instruments. Inter-calibration between instruments is a central pillar of the Cal/Val strategies of many national and international satellite remote sensing organizations. Inter-calibration techniques as outlined in this paper not only provide a practical means of identifying and correcting relative biases in radiometric calibration between instruments but also enable potential data gaps between measurement records in a critical time series to be bridged. Use of a robust set of internationally agreed upon and coordinated inter-calibration techniques will lead to significant improvement in the consistency between satellite instruments and facilitate accurate monitoring of the Earth’s climate at uncertainty levels needed to detect and attribute the mechanisms of change. This paper summarizes the state-of-the-art of post-launch radiometric calibration of remote sensing satellite instruments, through inter-calibration.
Humans now play a major role in altering Earth and its biota. Finding ways to ameliorate human impacts on biodiversity and to sustain and restore the ecosystem services on which we depend is a grand scientific and societal challenge. Conservation paleobiology is an emerging discipline that uses geohistorical data to meet these challenges by developing and testing models of how biota respond to environmental stressors. Here we ( a ) describe how the discipline has already provided insights about biotic responses to key environmental stressors, ( b ) outline research aimed at disentangling the effects of multiple stressors, ( c ) provide examples of deliverables for managers and policy makers, and ( d ) identify methodological advances in geohistorical analysis that will foster the next major breakthroughs in conservation outcomes. We highlight cases for which exclusive reliance on observations of living biota may lead researchers to erroneous conclusions about the nature and magnitude of biotic change, vulnerability, and resilience.
With the recent successful accounting of basin depth ground-motion adjustments in seismic hazard analyses for select areas of the western United States, we move toward implementing similar adjustments in the Atlantic and Gulf Coastal Plains by constructing a sediment thickness model and evaluating multiple relevant site amplification models for central and eastern United States seismic hazard analyses. We digitize and combine existing sediment thickness data sets into a composite surface that delineates the base of Cretaceous sediments under the Atlantic Coastal Plain and the base of Mesozoic sediments under the Gulf Coastal Plain. Amplification models dependent on sediment thickness, site natural period, and source-to-site path length are compared with data sets of observed ground motions to evaluate the ability of the new models to improve ground motion estimates. We find that the amplification models can account for observed trends in sediment-thickness and period-dependent residuals, but some tuning is required. For example, the model of Chapman and Guo requires a reference V S 30 , the time-averaged shear-wave velocity within 30 m of the Earth’s surface, for non-Coastal Plain sites, which we estimate to be between about 1 and 2 km/s. Along with our sediment thickness model, we estimate a velocity profile for application to the Harmon et al. site-natural-period-based model in order to best match the Chapman and Guo period dependence for a broad range of sediment thicknesses. The Next Generation of Attenuation models for the eastern United States Gulf Coast path-based adjustment models can also account for seismic attenuation in the Coastal Plain sediments and reduce the standard deviation of total residuals. If enacted in the U.S. Geological Survey National Seismic Hazard Model, these amplification models will reduce predicted short-period (<1 s) and increase predicted long-period (>1 s) ground motions in the Coastal Plains appreciably.
Anthropogenic climate change is driving sea level rise, leading to numerous impacts on the coastal zone, such as increased coastal flooding, beach erosion, cliff failure, saltwater intrusion in aquifers, and groundwater inundation. Many beaches around the world are currently experiencing chronic erosion as a result of gradual, present-day rates of sea level rise (about 3 mm/year) and human-driven restrictions in sand supply (e.g., harbor dredging and river damming). Accelerated sea level rise threatens to worsen coastal erosion and challenge the very existence of natural beaches throughout the world. Understanding and predicting the rates of sea level rise and coastal erosion depends on integrating data on natural systems with computer simulations. Although many computer modeling approaches are available to simulate shoreline change, few are capable of making reliable long-term predictions needed for full adaption or to enhance resilience. Recent advancements have allowed convincing decadal to centennial-scale predictions of shoreline evolution. For example, along 500 km of the Southern California coast, a new model featuring data assimilation predicts that up to 67% of beaches may completely erode by 2100 without large-scale human interventions. In spite of recent advancements, coastal evolution models must continue to improve in their theoretical framework, quantification of accuracy and uncertainty, computational efficiency, predictive capability, and integration with observed data, in order to meet the scientific and engineering challenges produced by a changing climate.
Results from computational morphodynamics modeling of coupled flow-bed-sediment systems are described for ten applications as a review of recent advances in the field. Each of these applications is drawn from solvers included in the public-domain International River Interface Cooperative (iRIC) software package. For mesoscale river features such as bars, predictions of alternate and higher mode river bars are shown for flows with equilibrium sediment supply and for a single case of oversupplied sediment. For microscale bed features such as bedforms, computational results are shown for the development and evolution of two-dimensional bedforms using a simple closure-based two-dimensional model, for two- and three-dimensional ripples and dunes using a three-dimensional large-eddy simulation flow model coupled to a physics-based particle transport model, and for the development of bed streaks using a three-dimensional unsteady Reynolds-averaged Navier-Stokes solver with a simple sediment-transport treatment. Finally, macroscale or channel evolution treatments are used to examine the temporal development of meandering channels, a failure model for cantilevered banks, the effect of bank vegetation on channel width, the development of channel networks in tidal systems, and the evolution of bedrock channels. In all examples, computational morphodynamics results from iRIC solvers are shown to compare well to observations of natural bed morphology. For each of the three scales investigated here, brief suggestions for future work and potential research directions are offered.
Obstacle marks are instream bedforms, typically composed of an upstream frontal scour hole and a downstream sediment accumulation in the vicinity of an obstacle. Local scouring at infrastructure (e.g. bridge piers) is a well‐studied phenomenon in hydraulic engineering, while less attention is given to the time‐dependent evolution of frontal scour holes at instream boulders and their geometric relations (depth to width, and length ratio). Furthermore, a comparison between laboratory studies and field observations is rare. Therefore, the morphodynamic importance of such scour features to fluvial sediment transport and morphological change is largely unknown. In this study, obstacle marks at boulder‐like obstructions were physically modelled in 30 unscaled process‐focused flume experiments (runtime per experiment ≥ 5760 min) at a range of flows (subcritical, clear‐water conditions, emergent and submerged water levels) and boundary conditions designed to represent the field setting (i.e. obstacle tilting, and limited thickness of the alluvial layer). Additionally, geometries of scour holes at 90 in‐situ boulders (diameter ≥ 1 m) located in a 50‐km segment of the Colorado River in Marble Canyon (AZ) were measured from a 1 m‐resolution digital elevation model. Flume experiments reveal similar evolution of local scouring, irrespective of hydraulic conditions, controlled by the scour incision, whereas the thickness of the alluvial layer and obstacle tilting into the evolving frontal scour hole limit incision. Three temporal evolution phases—(1) rapid incision, (2) decreasing incision, and (3) scour widening—are identified based on statistical analysis of spatiotemporal bed elevation time series. A quantitative model is presented that mechanistically predicts enlargement in local scour length and width based on (1) scour depth, (2) the inclination of scour slopes, and (3) the planform area of the frontal scour hole bottom. The comparison of field observations and laboratory results demonstrates scale invariance of geometry, which implies similitude of processes and form rather than equifinality.
Over the past several decades, the Yellowstone caldera has experienced frequent earthquake swarms and repeated cycles of uplift and subsidence, reflecting dynamic volcanic and tectonic processes. Here, we examine the detailed spatial-temporal evolution of the 2010 Madison Plateau swarm, which occurred near the northwest boundary of the Yellowstone caldera. To fully explore the evolution of the swarm, we integrated procedures for seismic waveform-based earthquake detection with precise double-difference relative relocation. Using cross-correlation of continuous seismic data and waveform templates constructed from cataloged events, we detected and precisely located 8710 earthquakes during the three-week swarm, nearly four times the number of events included in the standard catalog. This high-resolution analysis reveals distinct migration of earthquake activity over the course of the swarm. The swarm initiated abruptly on January 17, 2010 at about 10 km depth and expanded dramatically outward (both shallower and deeper) over time, primarily along a NNW-striking, ~55º ENE-dipping structure. To explain these characteristics, we hypothesize that the swarm was triggered by the rupture of a zone of confined high-pressure aqueous fluids into a pre-existing crustal fault system, prompting release of accumulated stress. The high-pressure fluid injection may have been accommodated by hybrid shear and dilatational failure, as is commonly observed in exhumed hydrothermally affected fault zones. This process has likely occurred repeatedly in Yellowstone as aqueous fluids exsolved from magma migrate into the brittle crust, and it may be a key element in the observed cycles of caldera uplift and subsidence.
Under favorable circumstances the time-dependent aseismic deformation resulting from the loading of the lithosphere by the stress drop of large dip slip earthquakes can be used to determine both the effective elastic plate thickness and the asthenospheric viscosity. The deformation has several similarities with the deflection of the lithosphere by surface loads and with movements due to postglacial rebound. Level changes obtained in the 80 years since the M = 7.5, 1896 Riku-u earthquake, an intraplate thrust event in northern Honshu, provide convincing evidence that asthenospheric readjustments are responsible for the observed movements. Leveling surveys crossing the zone of surface faulting have been repeated five times since 1900 and delineate a localized depression that has subsided at a continually decreasing rate. The depression is centered close to the 1896 faulting, and its shape and width, about 75 km, are matched by our model using a plate thickness of 30 km. The decaying subsidence rate constrains the viscosity of the uppermost asthenosphere to be 1×10 20 P. A linear viscous rheology matches the observed decay quite well, although measurements are sparse during the several decades following the earthquake.
Many subsurface processes involve transitions in granular material states, from arrested to creeping to flowing. Experiments and frameworks for idealized systems reveal that granular fabrics develop during shearing, co-evolve with applied stress, and govern such transitions. We use microtomography to test whether fabrics at two San Andreas fault sites reflect slip history and whether idealized frameworks extend to nature. Near-surface sediments within the fault zone transition between deformation patterns over the seismic cycle, including bulk/localized grain re-arrangements, individual grain fracturing, and localized zones of fracturing. Aseismic and co-seismic shearing produce distinct preferred grain orientations. Co-seismic fabrics can be preserved after centuries of aseismic strain, aseismic fabrics may be overprinted, and grain size and coordination number influence the fabrics. Idealized frameworks, namely anisotropic critical state theory, frictional jamming, and material memory, can explain our observations, and fault zone sediments likely undergo cycles of memory creation and erasure that influence rigidity spatiotemporally.