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 1,441 records · Page 80Linked to original sources
USGS2005

Source mechanism of Vulcanian degassing at Popocatépetl Volcano, Mexico, determined from waveform inversions of very long period signals

The source mechanism of very long period (VLP) signals accompanying volcanic degassing bursts at Popocatépetl is analyzed in the 15–70 s band by minimizing the residual error between data and synthetics calculated for a point source embedded in a homogeneous medium. The waveforms of two eruptions (23 April and 23 May 2000) representative of mild Vulcanian activity are well reproduced by our inversion, which takes into account volcano topography. The source centroid is positioned 1500 m below the western perimeter of the summit crater, and the modeled source is composed of a shallow dipping crack (sill with easterly dip of 10°) intersecting a steeply dipping crack (northeast striking dike dipping 83° northwest), whose surface extension bisects the vent. Both cracks undergo a similar sequence of inflation, deflation, and reinflation, reflecting a cycle of pressurization, depressurization, and repressurization within a time interval of 3–5 min. The largest moment release occurs in the sill, showing a maximum volume change of 500–1000 m 3 , pressure drop of 3–5 MPa, and amplitude of recovered pressure equal to 1.2 times the amplitude of the pressure drop. In contrast, the maximum volume change in the dike is less (200–300 m 3 ), with a corresponding pressure drop of 1–2 MPa and pressure recovery equal to the pressure drop. Accompanying these volumetric sources are single-force components with magnitudes of 10 8 N, consistent with melt advection in response to pressure transients. The source time histories of the volumetric components of the source indicate that significant mass movement starts within the sill and triggers a mass movement response in the dike within a few seconds. Such source behavior is consistent with the opening of a pathway for escape of pent-up gases from slow pressurization of the sill driven by magma crystallization. The opening of this pathway and associated rapid evacuation of volcanic gases induces the pressure drop. Pressure recovery in the magma filling the sill is driven by diffusion of gases from the resulting supersaturated melt into bubbles. Assuming a penny-shaped crack at ambient pressure of 40 MPa, the observed pressure and volume variations can be modeled with the following attributes: crack radius (100 m), crack aperture (5 m), bubble number density (10 10 –10 12 m −3 ), initial bubble radius (10 −6 m), final bubble radius (∼10 −5 m), and net decrease of gas concentration in the melt (0.01 wt %).

Popocatépetl Volcano↗
USGS2008

Postearthquake relaxation after the 2004 M6 Parkfield, California, earthquake and rate-and-state friction

An unusually complete set of measurements (including rapid rate GPS over the first 10 days) of postseismic deformation is available at 12 continuous GPS stations located close to the epicenter of the 2004 M6.0 Parkfield earthquake. The principal component modes for the relaxation of the ensemble of those 12 GPS stations were determined. The first mode alone furnishes an adequate approximation to the data. Thus, the relaxation at all stations can be represented by the product of a common temporal function and distinct amplitudes for each component (north or east) of relaxation at each station. The distribution in space of the amplitudes indicates that the relaxation is dominantly strike slip. The temporal function, which spans times from about 5 min to 900 days postearthquake, can be fit by a superposition of three creep terms, each of the form ??l loge(1 + t/??l), with characteristic times ??, = 4.06, 0.11, and 0.0001 days. It seems likely that what is actually involved is a broad spectrum of characteristic times, the individual components of which arise from afterslip on different fault patches. Perfettini and Avouac (2004) have shown that an individual creep term can be explained by the spring-slider model with rate-dependent (no state variable) friction. The observed temporal function can also be explained using a single spring-slider model (i.e., single fault patch) that includes rate-and-state-dependent friction, a single-state variable, and either of the two commonly used (aging and slip) state evolution laws. In the latter fits, the rate-and-state friction parameter b is negative.

Journal of Geophysical Research B: Solid Earth↗
USGS1999

Overview of the Mars Pathfinder Mission: Launch through landing, surface operations, data sets, and science results

Mars Pathfinder successfully landed at Ares Vallis on July 4, 1997, deployed and navigated a small rover about 100 m clockwise around the lander, and collected data from three science instruments and ten technology experiments. The mission operated for three months and returned 2.3 Gbits of data, including over 16,500 lander and 550 rover images, 16 chemical analyses of rocks and soil, and 8.5 million individual temperature, pressure and wind measurements. Path‐finder is the best known location on Mars, having been clearly identified with respect to other features on the surface by correlating five prominent horizon features and two small craters in lander images with those in high‐resolution orbiter images and in inertial space from two‐way ranging and Doppler tracking. Tracking of the lander has fixed the spin pole of Mars, determined the precession rate since Viking 20 years ago, and indicates a polar moment of inertia, which constrains a central metallic core to be between 1300 and ∼2000 km in radius. Dark rocks appear to be high in silica and geochemically similar to anorogenic andesites; lighter rocks are richer in sulfur and lower in silica, consistent with being coated with various amounts of dust. Rover and lander images show rocks with a variety of morphologies, fabrics and textures, suggesting a variety of rock types are present. Rounded pebbles and cobbles on the surface as well as rounded bumps and pits on some rocks indicate these rocks may be conglomerates (although other explanations are also possible), which almost definitely require liquid water to form and a warmer and wetter past. Air‐borne dust is composed of composite silicate particles with a small fraction of a highly magnetic mineral, interpreted to be most likely maghemite; explanations suggest iron was dissolved from crustal materials during an active hydrologic cycle with maghemite freeze dried onto silicate dust grains. Remote sensing data at a scale of a kilometer or greater and an Earth analog correctly predicted a rocky plain safe for landing and roving with a variety of rocks deposited by catstrophic floods, which are relatively dust free. The surface appears to have changed little since it formed billions of years ago, with the exception that eolian activity may have deflated the surface by ∼3–7 cm, sculpted wind tails, collected sand into dunes, and eroded ventifacts (fluted and grooved rocks). Pathfinder found a dusty lower atmosphere, early morning water ice clouds, and morning near‐surface air temperatures that changed abruptly with time and height. Small scale vortices, interpreted to be dust devils, were observed repeatedly in the afternoon by the meteorology instruments and have been imaged.

Journal of Geophysical Research E: Planets↗
USGS1996

Data base for a national mineral-resource assessment of undiscovered deposits of gold, silver, copper, lead, and zinc in the conterminous United States

For this assessment, the conterminous United States was divided into 12 regions Adirondack Mountains, Central and Southern Rocky Mountains, Colorado Plateau, East Central, Great Basin, Great Plains, Lake Superior, Northern Appalachians, Northern Rocky Mountains, Pacific Coast, Southern Appalachians, and Southern Basin and Range. The assessment, which was conducted by regional assessment teams of scientists from the USGS, was based on the concepts of permissive tracts and deposit models. Permissive tracts are discrete areas of the United States for which estimates of numbers of undiscovered deposits of a particular deposit type were made. A permissive tract is defined by its geographic boundaries such that the probability of deposits of the type delineated occurring outside the boundary is neglible. Deposit models, which are based on a compilation of worldwide literature and on observation, are sets of data in a convenient form that describe a group of deposits which have similar characteristics and that contain information on the common geologic attributes of the deposits and the environments in which they are found. Within each region, the assessment teams delineated permissive tracts for those deposit models that were judged to be appropriate and, when the amount of information warranted, estimated the number of undiscovered deposits. A total of 46 deposit models were used to assess 236 separate permissive tracts. Estimates of undiscovered deposits were limited to a depth of 1 km beneath the surface of the Earth. The estimates of the number of undiscovered deposits of gold, silver, copper, lead, and zinc were expressed in the form of a probability distribution. Commonly, the number of undiscovered deposits was estimated at the 90th, 50th, and 10th percentiles. A Monte Carlo simulation computer program was used to combine the probability distribution of the number of undiscovered deposits with the grade and tonnage data sets associated with each deposit model to obtain the probability distribution for undiscovered metal.

conterminous United States↗
USGS1915

Erosion and sedimentation in Chesapeake Bay around the mouth of Choptank River

With the unfolding of geologic knowledge during the last century the processes of denudation, transportation of sediments, and sedimentation have become better understood, and to some extent their relative effects in bringing about the present configuration of the earth's surface have been determined. The nature of these processes has been studied in many parts of the globe, but owing to the large size of the units affected and the slowness of the processes there has been little opportunity to collect quantitative data. Indeed, the data available are very largely conjectural, their degree of accuracy being only that of good guessing. However, with the advent of accurate topographic and hydrographic surveys, the first steps toward actual measurements of some of the many interesting surficial changes that are in progress have been taken, and it only remains for time and additional observations to afford opportunities for comparison. Resurveys, particularly by the United States Coast and Geodetic Survey, of parts of the eastern coast of the United States at different intervals have shown in many places marked changes of both shore line and sea bottom, as on Cape Cod, on Nantucket Island, and in Delaware Bay. Comparisons of this kind have so far been incidental to work in engineering and coast surveying, so that there has been little opportunity for the selection of localities particularly adapted to studies of erosion and sedimentation.

Maryland↗
USGS2011

The dynamical core, physical parameterizations, and basic simulation characteristics of the atmospheric component AM3 of the GFDL global coupled model CM3

The Geophysical Fluid Dynamics Laboratory (GFDL) has developed a coupled general circulation model (CM3) for the atmosphere, oceans, land, and sea ice. The goal of CM3 is to address emerging issues in climate change, including aerosol–cloud interactions, chemistry–climate interactions, and coupling between the troposphere and stratosphere. The model is also designed to serve as the physical system component of earth system models and models for decadal prediction in the near-term future—for example, through improved simulations in tropical land precipitation relative to earlier-generation GFDL models. This paper describes the dynamical core, physical parameterizations, and basic simulation characteristics of the atmospheric component (AM3) of this model. Relative to GFDL AM2, AM3 includes new treatments of deep and shallow cumulus convection, cloud droplet activation by aerosols, subgrid variability of stratiform vertical velocities for droplet activation, and atmospheric chemistry driven by emissions with advective, convective, and turbulent transport. AM3 employs a cubed-sphere implementation of a finite-volume dynamical core and is coupled to LM3, a new land model with ecosystem dynamics and hydrology. Its horizontal resolution is approximately 200 km, and its vertical resolution ranges approximately from 70 m near the earth’s surface to 1 to 1.5 km near the tropopause and 3 to 4 km in much of the stratosphere. Most basic circulation features in AM3 are simulated as realistically, or more so, as in AM2. In particular, dry biases have been reduced over South America. In coupled mode, the simulation of Arctic sea ice concentration has improved. AM3 aerosol optical depths, scattering properties, and surface clear-sky downward shortwave radiation are more realistic than in AM2. The simulation of marine stratocumulus decks remains problematic, as in AM2. The most intense 0.2% of precipitation rates occur less frequently in AM3 than observed. The last two decades of the twentieth century warm in CM3 by 0.32°C relative to 1881–1920. The Climate Research Unit (CRU) and Goddard Institute for Space Studies analyses of observations show warming of 0.56° and 0.52°C, respectively, over this period. CM3 includes anthropogenic cooling by aerosol–cloud interactions, and its warming by the late twentieth century is somewhat less realistic than in CM2.1, which warmed 0.66°C but did not include aerosol–cloud interactions. The improved simulation of the direct aerosol effect (apparent in surface clear-sky downward radiation) in CM3 evidently acts in concert with its simulation of cloud–aerosol interactions to limit greenhouse gas warming.

Journal of Climate↗
USGS2011

Selected Images of the Effects of the October 15, 2006, Kiholo Bay-Mahukona, Hawai'i, Earthquakes and Recovery Efforts

Early on the morning of October 15, 2006, two moderate earthquakes—the largest in decades—struck the Island of Hawai‘i. The first of these, which occurred at 7:07 a.m., HST (1707 UTC), was a magnitude ( M ) 6.7 earthquake, centered beneath Kīholo Bay on the northwestern coast of the island (19.878°N, 155.935°W), at a depth of 39 km. The second earthquake, which struck 6 minutes, 24 seconds later, at 7:14 a.m., HST (1714 UTC), was located 28 km to the north-northwest of Kīholo Bay (20.129°N, 155.983°W), centered at a depth of 19 km. This M 6.0 earthquake has since been referred to as the Māhukona earthquake. Losses from the combined effects of these earthquakes are estimated to be $200 million—the most costly events, by far, in Hawai‘i’s earthquake history. Although the vast majority of earthquakes in the State of Hawaii are closely related to the active volcanism associated with the southeastern part of the Island of Hawai‘i, the October 2006 Kīholo Bay and Māhukona earthquakes clearly suggest the devastating potential of deeper lithospheric earthquakes. Large earthquakes thought to be nearly M 7 have struck near the islands of Lāna‘i (1871) and Maui (1938). It is thought that these, like the 2006 earthquakes, were deep lithospheric flexure earthquakes (Wyss and Koyanagi, 1992; Klein and others, 2001). Thus, it is important to recognize the potential seismic hazard posed by such earthquakes beneath the older Hawaiian Islands. The data and observations afforded by the 2006 earthquakes promise to improve probabilistic seismic hazards modeling in Hawai‘i. The effects of the October 15, 2006, Kīholo Bay-Māhukona earthquakes are shown in images taken from the coastal route along the northern half of the Island of Hawai‘i, where damage was the most concentrated. The direction of presentation is counter-clockwise, from Pa‘auilo on the eastern or windward (Hāmākua) side to Kealakekua Bay on the western or leeward (Kona) side. A list of sites, their locations, coordinates, and distance from the epicenter at Kīholo Bay are given in table 1. A Google Earth map (fig. 7) and a topographic map (fig. 8) pinpoint the 36 sites where damage was documented and digital images were compiled for this collection.

Data Series↗
USGS2021

Solar-system-wide significance of Mars polar science

1. The North Polar Layered Deposits contain thousands of ice layers that record accumulation and climate history for at least several million years, making the most accessible and most complete climate record aside from the Earth’s – and the only one to record the impact of large obliquity shifts. 2. Mars Polar Science is a diverse and integrated system spanning much of the planet, above and below the surface: geographic poles, mid- and high-latitudes, numerous active volatile processes that interact with the near-surface, and the lower and middle atmosphere circulations, even influencing the upper atmosphere and atmospheric escape. 3. Mars polar processes serves as an analogue and scientific link to numerous solar system bodies. Mars has two volatile cycles that influence surface processes, atmospheric processes, and climate cycles. Mars, with many data sets and observations, serves as a proxy for outer planets and moons that undergo the same physics but with other volatile species. Additionally, technology designed for Mars polar exploration benefits exploration elsewhere. 4. Mars Polar and Atmospheric Science results are critical for future human exploration and in situ resources and the study of planetary habitability.

Bulletin of the AAS↗
USGS2020

North Carolina State climate report

Our scientific understanding of the climate system strongly supports the conclusion that North Carolina’s climate has changed in recent decades and the expectation that large changes—much larger than at any time in the state’s history—will occur if current trends in greenhouse gas concentrations continue. Even under a scenario where emissions peak around 2050 and decline thereafter, North Carolina will experience substantial changes in climate. The projected changes with the highest level of scientific confidence include increases in temperature, increases in summer absolute humidity, increases in sea level, and increases in extreme precipitation. It is also likely that there will be increases in the intensity of the strongest hurricanes. A full appreciation for past and future changes in North Carolina’s climate requires a global perspective. Earth’s climate has warmed substantially since the late 19th century, with most of that warming occurring in the last 50 years. This warming trend is clear from global temperature records and many other indicators, including rising global sea levels and rapid decreases in arctic sea ice cover. Scientists have very high confidence that this warming is largely due to human activities that have significantly increased atmospheric concentrations of carbon dioxide (CO2) and other greenhouse gases. Exhaustive research has examined other potential causes of this warming, and the increase in greenhouse gas concentrations is the only plausible cause that is consistent with the observed data and the physics that governs the climate system.

North Carolina↗
USGS2002

Waveform inversion of oscillatory signatures in long-period events beneath volcanoes

The source mechanism of long-period (LP) events is examined using synthetic waveforms generated by the acoustic resonance of a fluid-filled crack. We perform a series of numerical tests in which the oscillatory signatures of synthetic LP waveforms are used to determine the source time functions of the six moment tensor components from waveform inversions assuming a point source. The results indicate that the moment tensor representation is valid for the odd modes of crack resonance with wavelengths 2 L / n , 2 W / n , n = 3, 5, 7, …, where L and W are the crack length and width, respectively. For the even modes with wavelengths 2 L / n , 2 W / n , n = 2, 4, 6, …, a generalized source representation using higher-order tensors is required, although the efficiency of seismic waves radiated by the even modes is expected to be small. We apply the moment tensor inversion to the oscillatory signatures of an LP event observed at Kusatsu-Shirane Volcano, central Japan. Our results point to the resonance of a subhorizontal crack located a few hundred meters beneath the summit crater lakes. The present approach may be useful to quantify the source location, geometry, and force system of LP events, and opens the way for moment tensor inversions of tremor.

Journal of Geophysical Research B: Solid Earth↗
USGS2021

Earthquake magnitude distributions on northern Caribbean faults from combinatorial optimization models

On-fault earthquake magnitude distributions are calculated for northern Caribbean faults using estimates of fault slip and regional seismicity parameters. Integer programming, a combinatorial optimization method, is used to determine the optimal spatial arrangement of earthquakes sampled from a truncated Gutenberg-Richter distribution that minimizes the global misfit in slip rates on a complex fault system. Slip rates and their uncertainty on major faults are derived from a previously published GPS block model for the region, with fault traces determined from offshore geophysical mapping and previously published onshore studies. The optimal spatial arrangement of the sampled earthquakes is compared with the 500-year history of earthquake observations. Rupture segmentation of the subduction interface along the Hispaniola-Puerto Rico Trench (PRT) fault and seismic coupling on the PRT fault appear to exert the primary control over this spatial arrangement. Introducing a rupture barrier for the Hispaniola-PRT fault northwest of Mona Passage, based on geophysical and seismicity observations, and assigning a low slip rate of 2 mm/yr on the PRT fault are most consistent with historical earthquakes in the region. The addition of low slip-rate secondary faults as well as segmentation of the Hispaniola and Septentrional strike-slip fault improves the consistency with historical seismicity. An important observation from the modeling is that varying the slip rate on the PRT fault and different segmentation scenarios result in significant changes to the optimal magnitude distribution on faults farther away. In general, optimal on-fault magnitude distributions are more complex and inter-dependent than is typically assumed in probabilistic seismic hazard analysis and probabilistic tsunami hazard analysis.

JGR Solid Earth↗
USGS2005

Ca cycling and isotopic fluxes in forested ecosystems in Hawaii

Biogeochemical processes fractionate Ca isotopes in plants and soils along a 4 million year developmental sequence in the Hawaiian Islands. We observed that plants preferentially take up 40Ca relative to 44Ca, and that biological fractionation and changes in the relative contributions from volcanic and marine sources produce a significant increase in 44Ca in soil exchangeable pools. Our results imply moderate fluxes enriched in 44Ca from strongly nutrient-depleted old soils, in contrast with high 40Ca fluxes in young and little weathered environments. In addition, biological fractionation controls divergent geochemical pathways of Ca and Sr in the plant-soil system. While Ca depletes progressively with increasing soil age, Sr/Ca ratios increase systematically. Sr isotope ratios provide a valuable tracer for provenance studies of alkaline earth elements in forested ecosystems, but its usefulness is limited when deciphering biogeochemical processes involved in the terrestrial Ca cycle. Ca isotopes in combination with Sr/ Ca ratios reveal more complex processes involved in the biogeochemistry of Ca and Sr. Copyright 2005 by the American Geophysical Union.

Geophysical Research Letters↗
USGS1996

Earthquake mechanisms from linear-programming inversion of seismic-wave amplitude ratios

The amplitudes of radiated seismic waves contain far more information about earthquake source mechanisms than do first-motion polarities, but amplitudes are severely distorted by the effects of heterogeneity in the Earth. This distortion can be reduced greatly by using the ratios of amplitudes of appropriately chosen seismic phases, rather than simple amplitudes, but existing methods for inverting amplitude ratios are severely nonlinear and require computationally intensive searching methods to ensure that solutions are globally optimal. Searching methods are particularly costly if general (moment tensor) mechanisms are allowed. Efficient linear-programming methods, which do not suffer from these problems, have previously been applied to inverting polarities and wave amplitudes. We extend these methods to amplitude ratios, in which formulation on inequality constraint for an amplitude ratio takes the same mathematical form as a polarity observation. Three-component digital data for an earthquake at the Hengill-Grensdalur geothermal area in southwestern Iceland illustrate the power of the method. Polarities of P , SH , and SV waves, unusually well distributed on the focal sphere, cannot distinguish between diverse mechanisms, including a double couple. Amplitude ratios, on the other hand, clearly rule out the double-couple solution and require a large explosive isotropic component.

Bulletin of the Seismological Society of America↗
USGS1978

Water-level monitoring in the area of the Palmdale Uplift, Southern California

Abnormal behavior of water levels in wells has been observed prior to a number of earthquakes. For instance, water-level minima have been noted in the Cienega Winery well before earthquakes on the San Andreas fault. Abnormal water-level fluctuations were used in conjunctions with other precursors to predict the February 4, 1975, Haicheng earthquake in northeastern China. That such changes should occur prior to earthquakes is not surprising. Ground water that occupies the void spaces in porous rocks or alluvium can be expected to rise in wells when an aquifer is squeezed and fall when it is distended. COnfined aquifers, in particualr, have been found to be highly sensitive to Earth strain changes.

California↗
USGS2012

Experimental deformation of topaz crystals: Possible embrittlement by intracrystalline water

Crystallographically oriented single‐crystal prisms of gem quality topaz (composition AlSiO (OH− F) where x = 0.04 ± 0.01) were deformed at a confining pressure of 1.50 GPa, a temperature of 800°C, and a strain rate of 2×10 s. Under nearly identical conditions, all crystals of anhydrous rock‐forming minerals that have been tested to date, such as olivine, quartz, feldspars, pyroxenes, and refractory oxides, deform plastically; in contrast, our topaz crystals failed by brittle fracture regardless of the orientation of the compression direction. No optical evidence for plastic deformation was detected. Another suite of experiments with compression perpendicular to the (001) cleavage at = 100°–950°C and a strain rate of 2×10 s displayed two regimes of behavior: (1) at >400°C, fracture strength was independent of temperature, and fracture occurred on one or two surfaces parallel to {103}; (2) at <400°C, the fracture strength increased rapidly with decreasing temperature, no macroscopic stress drop was observed, and many closely spaced conjugate fractures formed on (103) and (103). The anomalous brittleness of topaz compared to anhydrous silicate and oxide crystals indicates that intracrystalline “water” plays a role in the embrittlement. We suggest that water within the topaz crystals promotes fracture in ways similar to the mechanisms of slow crack growth aided by environmental moisture.

Journal of Geophysical Research B: Solid Earth↗
USGS2005

Questa baseline and pre-mining ground-water quality investigation. 5. Well installation, water-level data, and surface- and ground-water geochemistry in the Straight Creek drainage basin, Red River Valley, New Mexico, 2001-03

The U.S. Geological Survey, in cooperation with the New Mexico Environment Department, is investigating the pre-mining ground-water chemistry at the Molycorp molybdenum mine in the Red River Valley, northern New Mexico. The primary approach is to determine the processes controlling ground-water chemistry at an unmined, off-site, proximal analog. The Straight Creek drainage basin, chosen for this purpose, consists of the same quartz-sericite-pyrite altered andesitic and rhyolitic volcanic rock of Tertiary age as the mine site. The weathered and rugged volcanic bedrock surface is overlain by heterogeneous debris-flow deposits that interfinger with alluvial deposits near the confluence of Straight Creek and the Red River. Pyritized rock in the upper part of the drainage basin is the source of acid rock drainage (pH 2.8-3.3) that infiltrates debris-flow deposits containing acidic ground water (pH 3.0-4.0) and bedrock containing water of circumneutral pH values (5.6-7.7). Eleven observation wells were installed in the Straight Creek drainage basin. The wells were completed in debris-flow deposits, bedrock, and interfingering debris-flow and Red River alluvial deposits. Chemical analyses of ground water from these wells, combined with chemical analyses of surface water, water-level data, and lithologic and geophysical logs, provided information used to develop an understanding of the processes contributing to the chemistry of ground water in the Straight Creek drainage basin. Surface- and ground-water samples were routinely collected for determination of total major cations and selected trace metals; dissolved major cations, selected trace metals, and rare-earth elements; anions and alkalinity; and dissolved-iron species. Rare-earth elements were determined on selected samples only. Samples were collected for determination of dissolved organic carbon, mercury, sulfur isotopic composition (34S and 18O of sulfate), and water isotopic composition (2H and 18O) during selected samplings. One set of ground-water samples was collected for helium-3/tritium and chlorofluorocarbon (CFC) age dating. Several lines of evidence indicate that surface water is the primary input to the Straight Creek ground-water system. Straight Creek streamflow and water levels in wells closest to the apex of the Straight Creek debris fan and closest to Straight Creek itself appear to respond to the same seasonal inputs. Oxygen and hydrogen isotopic compositions in Straight Creek surface water and ground water are similar, and concentrations of most dissolved constituents in most Straight Creek surface-water and shallow (debris-flow and alluvial) aquifer ground-water samples correlate strongly with sulfate (concentrations decrease linearly with sulfate in a downgradient direction). After infiltration of surface water, dilution along the flow path is the dominant mechanism controlling ground-water chemistry. However, concentrations of some constituents can be higher in ground water than can be accounted for by concentrations in Straight Creek surface water, and additional sources of these constituents must therefore be inferred. Constituents for which concentrations in ground water can be high relative to surface water include calcium, magnesium, strontium, silica, sodium, and potassium in ground water from debris-flow and alluvial aquifers and manganese, calcium, magnesium, strontium, sodium, and potassium in ground water from the bedrock aquifer. All ground water is a calcium sulfate type, often at or near gypsum saturation because of abundant gypsum in the aquifer material developed from co-existing calcite and pyrite mineralization. Calcite dissolution, the major buffering mechanism for bedrock aquifer ground water, also contributes to relatively higher calcium concentrations in some ground water. The main source of the second most abundant cation, magnesium, is probably dissolution of magnesium-rich carbonates or silicates.

New Mexico↗
USGS2021

Limited mantle hydration by bending faults at the Middle America Trench

Seismic anisotropy measurements show that upper mantle hydration at the Middle America Trench (MAT) is limited to serpentinization and/or water in fault zones, rather than distributed uniformly. Subduction of hydrated oceanic lithosphere recycles water back into the deep mantle, drives arc volcanism, and affects seismicity at subduction zones. Constraining the extent of upper mantle hydration is an important part of understanding many fundamental processes on Earth. Substantially reduced seismic velocities in tomography suggest that outer rise plate‐bending faults provide a pathway for seawater to rehydrate the slab mantle just prior to subduction. Estimates of outer‐rise hydration based on tomograms vary significantly, with some large enough to imply that, globally, subduction has consumed more than two oceans worth of water during the Phanerozoic. We found that, while the mean upper mantle wavespeed is reduced at the MAT outer rise, the amplitude and orientation of inherited anisotropy are preserved at depths >1 km below the Moho. At shallower depths, relict anisotropy is replaced by slowing in the fault‐normal direction. These observations are incompatible with pervasive hydration but consistent with models of wave propagation through serpentinized fault zones that thin to <100‐m in width at depths >1 km below Moho. Confining hydration to fault zones reduces water storage estimates for the MAT upper mantle from ∼3.5 wt% to <0.9 wt% H 2 0. Since the intermediate thermal structure in the ∼24 Myr‐old MAT slab favors serpentinization, limited hydration suggests that fault mechanics are the limiting factor, not temperatures. Subducting mantle may be similarly dry globally.

Journal of Geophysical Research↗
USGS2023

Quantification of geodetic strain rate uncertainties and implications for seismic hazard estimates

Geodetic velocity data provide first-order constraints on crustal surface strain rates, which in turn are linked to seismic hazard. Estimating the 2-D surface strain tensor everywhere requires knowledge of the surface velocity field everywhere, while geodetic data such as Global Navigation Satellite System (GNSS) only have spatially scattered measurements on the surface of the Earth. To use these data to estimate strain rates, some type of interpolation is required. In this study, we review methodologies for strain rate estimation and compare a suite of methods, including a new implementation based on the geostatistical method of kriging, to compare variation between methods with uncertainty based on one method. We estimate the velocity field and calculate strain rates in southern California using a GNSS velocity field and five different interpolation methods to understand the sources of variability in inferred strain rates. Uncertainty related to data noise and station spacing (aleatoric uncertainty) is minimal where station spacing is dense and maximum far from observations. Differences between methods, related to epistemic uncertainty, are usually highest in areas of high strain rate due to differences in how gradients in the velocity field are handled by different interpolation methods. Parameter choices, unsurprisingly, have a strong influence on strain rate field, and we propose the traditional L -curve approach as one method for quantifying the inherent trade-off between fit to the data and models that are reflective of tectonic strain rates. Doing so, we find total variability between five representative strain rate models to be roughly 40 per cent, a much lower value than roughly 100 per cent that was found in previous studies (Hearn et al .). Using multiple methods to tune parameters and calculate strain rates provides a better understanding of the range of acceptable models for a given velocity field. Finally, we present an open-source Python package (Materna et al .) for calculating strain rates, Strain_2D, which allows for the same data and model grid to be used in multiple strain rate methods, can be extended with other methods from the community, and provides an interface for comparing strain rate models, calculating statistics and estimating strain rate uncertainty for a given GNSS data set.

Geophysical Journal International↗