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,729 records · Page 96Linked to original sources
USGS1996

Crustal structure of the Colorado Plateau, Arizona: Application of new long-offset seismic data analysis techniques

The Colorado Plateau is a large crustal block in the southwestern United States that has been raised intact nearly 2 km above sea level since Cretaceous marine sediments were deposited on its surface. Controversy exists concerning the thickness of the plateau crust and the source of its buoyancy. Interpretations of seismic data collected on the plateau vary as to whether the crust is closer to 40 or 50 km thick. A thick crust could support the observed topography of the Colorado Plateau isostatically, while a thinner crust would indicate the presence of an underlying low-density mantle. This paper reports results on long-offset seismic data collected during the 1989 segment of the U.S. Geological Survey Pacific to Arizona Crustal Experiment that extended from the Transition Zone into the Colorado Plateau in northwest Arizona. We apply two new methods to analyze long-offset data that employ finite difference travel time calculations: (1) a first-arrival time inverter to find upper crustal velocity structure and (2) a forward-modeling technique that allows the direct use of the inverted upper crustal solution in modeling secondary reflected arrivals. We find that the crustal thickness increases from 30 km beneath the metamorphic core complexes in the southern Basin and Range province to about 42 km beneath the northern Transition Zone and southern Colorado Plateau margin. We observe some crustal thinning (to ∼37 km thick) and slightly higher lower crustal velocities farther inboard; beneath the Kaibab uplift on the north rim of the Grand Canyon the crust thickens to a maximum of 48 km. We observe a nonuniform crustal thickness beneath the Colorado Plateau that varies by ∼15% and corresponds approximately to variations in topography with the thickest crust underlying the highest elevations. Crustal compositions (as inferred from seismic velocities) appear to be the same beneath the Colorado Plateau as those in the Basin and Range province to the southwest, implying that the plateau crust represents an unextended version of the Basin and Range. Some of the variability in crustal structure appears to correspond to preserved lithospheric discontinuities that date back to the Proterozoic Era.

Journal of Geophysical Research B: Solid Earth↗
USGS1972

Subduction zones: Not relevant to present-day problems of waste disposal

SUBDUCTION zones are considered to be sites of disposal for vast areas of the Earth's surface 1 , while new surface is generated simultaneously at rise crests 2 . Bostrom and Sherif 3 suggest that the world's industrial and domestic waste be dumped into subduction zones at deep sea trenches to allow nature to complete the recycling process at geologically rapid rates of 5 to 10 cm/yr. They also point out that trenches are often sites of rapid rates of deposition and suggest that the dumped wastes would, speaking geologically, soon be buried. Francis 4 suggests that canisters of toxic chemical and radioactive wastes could be dumped onto trench sediments and be expected to sink at rates of 20 m/yr, assuming that the mass of turbidites in the trench fill often spontaneously liquefies on shaking by earthquakes. The assumption is based on the supposed lack of evidence for deformed sediment in trenches. I will argue that the suggestion of Bostrom and Sherif 3 is not useful for the next few dozen generations of human populations and will point out observational evidence to show that Francis's 4 assumption is incorrectly founded.

Nature↗
USGS2016

An empirical approach for estimating stress-coupling lengths for marine-terminating glaciers

Variability in the dynamic behavior of marine-terminating glaciers is poorly understood, despite an increase in the abundance and resolution of observations. When paired with ice thicknesses, surface velocities can be used to quantify the dynamic redistribution of stresses in response to environmental perturbations through computation of the glacier force balance. However, because the force balance is not purely local, force balance calculations must be performed at the spatial scale over which stresses are transferred within glacier ice, or the stress-coupling length (SCL). Here we present a new empirical method to estimate the SCL for marine-terminating glaciers using high-resolution observations. We use the empirically-determined periodicity in resistive stress oscillations as a proxy for the SCL. Application of our empirical method to two well-studied tidewater glaciers (Helheim Glacier, SE Greenland, and Columbia Glacier, Alaska, USA) demonstrates that SCL estimates obtained using this approach are consistent with theory (i.e., can be parameterized as a function of the ice thickness) and with prior, independent SCL estimates. In order to accurately resolve stress variations, we suggest that similar empirical stress-coupling parameterizations be employed in future analyses of glacier dynamics.

Frontiers in Earth Science↗
USGS2014

The geology and geochemistry of Isla Floreana, Galápagos: A different type of late-stage ocean island volcanism

Isla Floreana, the southernmost volcano in the Galápagos Archipelago, has erupted a diverse suite of alkaline basalts continually since 1.5 Ma. Because these basalts have different compositions than xenoliths and older lavas from the deep submarine sector of the volcano, Floreana is interpreted as being in a rejuvenescent or late-stage phase of volcanism. Most lavas contain xenoliths, or their disaggregated remains. The xenolithic debris and large ranges in composition, including during single eruptions, indicate that the magmas do not reside in crustal magma chambers, unlike magmas in the western Galápagos. Floreana lavas have distinctive trace element compositions that are rich in fluid-immobile elements (e.g., Ta, Nb, Th, Zr) and even richer in fluid-mobile elements (e.g., Ba, Sr, Pb). Rare earth element (REE) patterns are light REE-enriched and distinctively concave-up. Neodymium isotopic ratios are comparable to those from Fernandina, at the core of the Galápagos plume, but Floreana has the most radiogenic Sr and Pb isotopic ratios in the archipelago. These trace element patterns and isotopic ratios are attributed to a mixed source originating within the Galápagos plume, which includes depleted upper mantle, plume material rich in TITAN elements (Ti, Ta, Nb), and recycled oceanic crust that has undergone partial dehydration in an ancient subduction zone. Because Floreana lies at the periphery of the Galápagos plume, melting occurs mostly in the spinel zone, and enriched components dominate; the Floreana recycled mantle component influence is detectable in volcanoes along the entire southern periphery of the archipelago as well. Floreana is the only Galápagos volcano known to have undergone late-stage volcanism. Here, however, the secondary stage activity is more compositionally enriched than the shield-building phase, in contrast to what is observed in Hawai‘i, suggesting that the mechanism driving late-stage volcanism may vary among ocean island provinces.

Galápagos, Isla Floreana↗
USGS2016

The role of crystallization-driven exsolution on the sulfur mass balance in volcanic arc magmas

The release of large amounts of sulfur to the stratosphere during explosive eruptions affects the radiative balance in the atmosphere and consequentially impacts climate for up to several years after the event. Quantitative estimations of the processes that control the mass balance of sulfur between melt, crystals, and vapor bubbles is needed to better understand the potential sulfur yield of individual eruption events and the conditions that favor large sulfur outputs to the atmosphere. The processes that control sulfur partitioning in magmas are (1) exsolution of volatiles (dominantly H 2 O) during decompression (first boiling) and during isobaric crystallization (second boiling), (2) the crystallization and breakdown of sulfide or sulfate phases in the magma, and (3) the transport of sulfur-rich vapor (gas influx) from deeper unerupted regions of the magma reservoir. Vapor exsolution and the formation/breakdown of sulfur-rich phases can all be considered as closed-system processes where mass balance arguments are generally easier to constrain, whereas the contribution of sulfur by vapor transport (open system process) is more difficult to quantify. The ubiquitous “excess sulfur” problem, which refers to the much higher sulfur mass released during eruptions than what can be accounted for by amount of sulfur originally dissolved in erupted melt, as estimated from melt inclusion sulfur concentrations (the “petrologic estimate”), reflects the challenges in closing the sulfur mass balance between crystals, melt, and vapor before and during a volcanic eruption. In this work, we try to quantify the relative importance of closed- and open-system processes for silicic arc volcanoes using kinetic models of sulfur partitioning during exsolution. Our calculations show that crystallization-induced exsolution (second boiling) can generate a significant fraction of the excess sulfur observed in crystal-rich arc magmas. This result does not negate the important role of vapor migration in sulfur mass balance but rather points out that second boiling (in situ exsolution) can provide the necessary yield to drive the excess sulfur to the levels observed for crystal-rich systems. In contrast, in crystal-poor systems, magma recharge that releases sulfur-rich bubbles is necessary and most likely the primary contributor to sulfur mass balance. Finally, we apply our model to account for the effect of sulfur partitioning during second boiling and its impact on sulfur released during the Cerro Galan supereruption in Argentina (2.08 Ma) and show the potential importance of second boiling in releasing a large amount of sulfur to the atmosphere during the eruption of large crystal-rich ignimbrites.

Journal of Geophysical Research B: Solid Earth↗
USGS1999

Nature and origin of upper crustal seismic velocity fluctuations and associated scaling properties: Combined stochastic analyses of KTB velocity and lithology logs

The main borehole of the German Continental Deep Drilling Program (KTB) extends over 9000 m into a crystalline upper crust consisting primarily of interlayered gneiss and metabasite. We present a joint analysis of the velocity and lithology logs in an effort to extract the lithology component of the velocity log. Covariance analysis of lithology log, approximated as a binary series, indicates that it may originate from the superposition of two Brownian stochastic processes (fractal dimension 1.5) with characteristic scales of ∼2800 m and ∼150 m, respectively. Covariance analysis of the velocity fluctuations provides evidence for the superposition of four stochastic process with distinct characteristic scales. The largest two scales are identical to those derived from the lithology, confirming that these scales of velocity heterogeneity are caused by lithology variations. The third characteristic scale, ∼20 m, also a Brownian process, is probably related to fracturing based on correlation with the resistivity log. The superposition of these three Brownian processes closely mimics the commonly observed 1/ k decay (fractal dimension 2.0) of the velocity power spectrum. The smallest scale process (characteristic scale ∼1.7 m) requires a low fractal dimension, ∼1.0, and accounts for ∼60% of the total rms velocity variation. A comparison of successive logs from 6900–7140 m depth indicates that such variations are not repeatable and thus probably do not represent true velocity variations in the crust. The results of this study resolve disparity between the differing published estimates of seismic heterogeneity based on the KTB sonic logs, and bridge the gap between estimates of crustal heterogeneity from geologic maps and borehole logs.

Journal of Geophysical Research B: Solid Earth↗
USGS2008

Variability of passive gas emissions, seismicity, and deformation during crater lake growth at White Island Volcano, New Zealand, 2002-2006

We report on 4 years of airborne measurements of CO2, SO2, and H2S emission rates during a quiescent period at White Island volcano, New Zealand, beginning in 2003. During this time a significant crater lake emerged, allowing scrubbig processes to be investigated. CO2 emissions varied from a baseline of 250 to >2000 t d-1 and demonstrated clear annual cycling that was consistent with numbers of earthquake detections and annual changes in sea level. The annual variability was found to be most likely related to increases in the strain on the volcano during sea level highs, temporarily causing fractures to reduce in size in the upper conduit. SO2 emissions varied from 0 to >400 t d-1 and were clearly affected by scrubbing processes within the first year of take development. Scrubbing caused increases of SO42- and Cl- in lake waters, and the ratio of carbon to total sulphur suggested that elemental sulphur deposition was also significant in the lake during the first year. Careful measurements of the lake level and chemistry allowed estimates of the rate of H2O(g) and HCl(g) input into the lake and suggested that the molar abundances of major gas species (H2O, CO2, SO2, and HCl) during this quiescent phase were similar to fumarolic ratios observed between earlier eruptive periods. The volume of magma estimated from CO2 emissions (0.0 15-0.04 km3) was validated by Cl- increases in the lake, suggesting that the gas and magma are transported from deep to shallow depths as a closed system and likely become open in the upper conduit region. The absence of surface deformation further leads to a necessity of magma convection to supply and remove magma from the degassing depths. Two models of convection configurations are discussed. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2023

The role of giant impacts in planet formation

Planets are expected to conclude their growth through a series of giant impacts: energetic, global events that significantly alter planetary composition and evolution. Computer models and theory have elucidated the diverse outcomes of giant impacts in detail, improving our ability to interpret collision conditions from observations of their remnants. However, many open questions remain, as even the formation of the Moon—a widely suspected giant-impact product for which we have the most information—is still debated. We review giant-impact theory, the diverse nature of giant-impact outcomes, and the governing physical processes. We discuss the importance of computer simulations, informed by experiments, for accurately modeling the impact process. Finally, we outline how the application of probability theory and computational advancements can assist in inferring collision histories from observations, and we identify promising opportunities for advancing giant-impact theory in the future. Giant impacts exhibit diverse possible outcomes leading to changes in planetary mass, composition, and thermal history depending on the conditions. Improvements to computer simulation methodologies and new laboratory experiments provide critical insights into the detailed outcomes of giant impacts. When colliding planets are similar in size, they can merge or escape one another with roughly equal probability, but with different effects on their resulting masses, densities, and orbits. Different sequences of giant impacts can produce similar planets, encouraging the use of probability theory to evaluate distinct formation hypothesis.

Annual Reviews of Earth and Planetary Science↗
USGS2009

The 2006-2007 Kuril Islands great earthquake sequence

[1] The southwestern half of a ∼500 km long seismic gap in the central Kuril Island arc subduction zone experienced two great earthquakes with extensive preshock and aftershock sequences in late 2006 to early 2007. The nature of seismic coupling in the gap had been uncertain due to the limited historical record of prior large events and the presence of distinctive upper plate, trench and outer rise structures relative to adjacent regions along the arc that have experienced repeated great interplate earthquakes in the last few centuries. The intraplate region seaward of the seismic gap had several shallow compressional events during the preceding decades (notably an M S 7.2 event on 16 March 1963), leading to speculation that the interplate fault was seismically coupled. This issue was partly resolved by failure of the shallow portion of the interplate megathrust in an M W = 8.3 thrust event on 15 November 2006. This event ruptured ∼250 km along the seismic gap, just northeast of the great 1963 Kuril Island (M w = 8.5) earthquake rupture zone. Within minutes of the thrust event, intense earthquake activity commenced beneath the outer wall of the trench seaward of the interplate rupture, with the larger events having normal-faulting mechanisms. An unusual double band of interplate and intraplate aftershocks developed. On 13 January 2007, an M W = 8.1 extensional earthquake ruptured within the Pacific plate beneath the seaward edge of the Kuril trench. This event is the third largest normal-faulting earthquake seaward of a subduction zone on record, and its rupture zone extended to at least 33 km depth and paralleled most of the length of the 2006 rupture. The 13 January 2007 event produced stronger shaking in Japan than the larger thrust event, as a consequence of higher short-period energy radiation from the source. The great event aftershock sequences were dominated by the expected faulting geometries; thrust faulting for the 2006 rupture zone, and normal faulting for the 2007 rupture zone. A large intraplate compressional event occurred on 15 January 2009 (M w = 7.4) near 45 km depth, below the rupture zone of the 2007 event and in the vicinity of the 16 March 1963 compressional event. The fault geometry, rupture process and slip distributions of the two great events are estimated using very broadband teleseismic body and surface wave observations. The occurrence of the thrust event in the shallowest portion of the interplate fault in a region with a paucity of large thrust events at greater depths suggests that the event removed most of the slip deficit on this portion of the interplate fault. This great earthquake doublet demonstrates the heightened seismic hazard posed by induced intraplate faulting following large interplate thrust events. Future seismic failure of the remainder of the seismic gap appears viable, with the northeastern region that has also experienced compressional activity seaward of the megathrust warranting particular attention.

Journal of Geophysical Research B: Solid Earth↗
USGS2008

Wind-driven particle mobility on Mars: Insights from Mars Exploration Rover observations at "El Dorado" and surroundings at Gusev Crater

The ripple field known as 'El Dorado' was a unique stop on Spirit's traverse where dust-raising, active mafic sand ripples and larger inactive coarse-grained ripples interact, illuminating several long-standing issues of Martian dust mobility, sand mobility, and the origin of transverse aeolian ridges. Strong regional wind events endured by Spirit caused perceptible migration of ripple crests in deposits SSE of El Dorado, erasure of tracks in sandy areas, and changes to dust mantling the site. Localized thermal vortices swept across El Dorado, leaving paths of reduced dust but without perceptibly damaging nearly cohesionless sandy ripple crests. From orbit, winds responsible for frequently raising clay-sized dust into the atmosphere do not seem to significantly affect dunes composed of (more easily entrained) sand-sized particles, a long-standing paradox. This disparity between dust mobilization and sand mobilization on Mars is due largely to two factors: (1) dust occurs on the surface as fragile, low-density, sand-sized aggregates that are easily entrained and disrupted, compared with clay-sized air fall particles; and (2) induration of regolith is pervasive. Light-toned bed forms investigated at Gusev are coarse-grained ripples, an interpretation we propose for many of the smallest linear, light-toned bed forms of uncertain origin seen in high-resolution orbital images across Mars. On Earth, wind can organize bimodal or poorly sorted loose sediment into coarse-grained ripples. Coarse-grained ripples could be relatively common on Mars because development of durable, well-sorted sediments analogous to terrestrial aeolian quartz sand deposits is restricted by the lack of free quartz and limited hydraulic sediment processing. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗
USGS2024

Remote sensing large-wood storage downstream of reservoirs during and after dam removal: Elwha River, Washington, USA

Large wood is an integral part of many rivers, often defining river-corridor morphology and habitat, but its occurrence, magnitude, and evolution in a river system are much less well understood than the sedimentary and hydraulic components, and due to methodological limitations, have seldom previously been mapped in substantial detail. We present a new method for this, representing a substantial advance in automated deep-learning-based image segmentation. From these maps, we measured large wood and sediment deposits from high-resolution orthoimages to explore the dynamics of large wood in two reaches of the Elwha River, Washington, USA, between 2012 and 2017 as it adjusted to upstream dam removals. The data set consists of a time series of orthoimages (12.5-cm resolution) constructed using Structure-from-Motion photogrammetry on imagery from 14 aerial surveys. Model training was optimized to yield maximum accuracy for estimated wood areas, compared to manually digitized wood, therefore model development and intended application were coupled. These fully reproducible methods and model resulted in a maximum of 15% error between observed and estimated total wood areas and wood deposit size-distributions over the full spatio-temporal extent of the data. Areal extent of wood in the channel margin approximately doubled in the years following dam removal, with greatest increases in large wood in wider, lower-gradient sections. Large-wood deposition increased between the start of dam removal (2011) and winter 2013, then plateaued. Sediment bars continued to grow up until 2016/17, assisted by a partially static wood framework deposited predominantly during the period up to winter 2013.

Washington↗
USGS2025

Relationship of atmospheric nitrogen deposition to soil nitrogen cycling along an elevation gradient in the Colorado Front Range

Microbial processing of atmospheric nitrogen (N) deposition regulates the retention and mobilization of N in soils, with important implications for water quality. Understanding the links between N deposition, microbial communities, N transformations, and water quality is critical as N deposition shifts toward reduced N and remains persistently high in many regions. Here, we investigated these connections along an elevation transect in the Colorado Front Range. Although rates of N deposition and pools of extractable N increased down the elevation transect, soil microbial communities and N transformation rates did not follow clear elevational patterns. The subalpine microbial community was distinct, corresponding to a high C:N ratio and low pH, while the microbial communities at the lower elevation sites were all very similar. Net nitrification, mineralization, and nitrification potential rates were highest at the Plains (1,700 m) and Montane (2,527 m) sites, suggesting that these ecosystems mobilize N. In contrast, the net immobilization of N observed at the Foothills (1,978 m) and Subalpine (3,015 m) sites suggests that these ecosystems retain N deposition. The contrast in N transformation rates between the plains and foothills, both of which receive elevated N deposition, may be due to spatial heterogeneity not captured in this study and warrants further investigation. Stream N concentrations from the subalpine to the foothills were consistently low, indicating that these soils are currently able to process and retain N deposition, but this may be disrupted if drought, wildfire, or land-use change alter the ability of the soils to retain N.

Colorado↗
USGS2012

Three-dimensional P and S velocity structure in the Coalinga Region, California

The Coalinga earthquake sequence of 1983 provided a unique opportunity to perform a three-dimensional velocity and hypocenter inversion in an area of complex three-dimensional structure dominated by folding and blind thrusts. Additionally, since other varied geological and geophysical studies have been completed in this area, the three-dimensional inversion solution could be compared to other interpretations. Inversion of 7696 P and 1511 S first arrivals from earthquakes and 696 P first arrivals from refraction shots produced a three-dimensional velocity model with grid spacing of 1–2 km in the hypocentral area. The overall shape and location of velocity features correspond well to the mapped surface geology. The three-dimensional inversion yields details of folds where the resolution is good and the general shape where resolution is lower. The amounts of structural relief inferred for the local folds are similar to values inferred from geologic data. The three-dimensional velocity solution has several distinctive features. There is a linear high-velocity body (6.1–6.5 km/s), about 25 km long, from 6 to 8 km depth, that may represent a fragment of Coast Range ophiolite. A shallow low-velocity zone (LVZ), which extends for 20 km along the fold axis at about 6 km depth and correlates with LVZs observed in both the refraction and the reflection data, may indicate high pore pressure caused by lateral compressive strain. Deeper LVZs occur within inferred Franciscan material and are characterized by horizontal or southwest dipping zones of varied thickness, 4–8 km wide and 5–10 km long. These LVZs may represent multiple thrust faults. Their locations and geometry are consistent with thrust faults inferred with seismic reflection data. The three-dimensional velocity solution compares well to prior two-dimensional seismic reflection and refraction models and observed gravity. The shape of the inferred sedimentary section agrees well with the reflections from the Cenozoic strata. Compared to the refraction model, the three-dimensional solution has similar velocities and similar locations of velocity features but is more detailed in the hypocentral zone where it uses more data. The gravity computed from the three-dimensional velocities is similar to the observed gravity in both shape and amplitude. Both a simple one-dimensional initial model and a complex initial model derived from the refraction interpretation were tried. A simple starting model gave the best results. The S velocity solution has different resolution than the P velocity solution because it uses a different set of stations, and it has lower resolution because it uses fewer arrival times. While the general patterns of velocity variation are similar for both V p and V s , the V s solution tends to have more smearing of velocity features and can have somewhat different locations of velocity features.

Journal of Geophysical Research Solid Earth↗
USGS2000

Stable isotope systematics of sulfate minerals

Stable isotope studies of sulfate minerals are especially useful for unraveling the geochemical history of geological systems. All sulfate minerals can yield sulfur and oxygen isotope data. Hydrous sulfate minerals, such as gypsum, also yield oxygen and hydrogen isotope data for the water of hydration, and more complex sulfate minerals, such as alunite and jarosite also yield oxygen and hydrogen isotope data from hydroxyl sites. Applications of stable isotope data can be divided into two broad categories: geothermometry and tracer studies. The equilibrium partitioning of stable isotopes between two substances, such as the isotopes of sulfur between barite and pyrite, is a function of temperature. Studies can also use stable isotopes as a tracer to fingerprint various sources of hydrogen, oxygen, and sulfur, and to identify physical and chemical processes such as evaporation of water, mixing of waters, and reduction of sulfate to sulfide. Studies of sulfate minerals range from low-temperature surficial processes associated with the evaporation of seawater to form evaporite deposits to high-temperature magmatic-hydrothermal processes associated with the formation of base-and precious-metal deposits. Studies have been conducted on scales from submicroscopic chemical processes associated with the weathering of pyrite to global processes affecting the sulfur budget of the oceans. Sulfate isotope studies provide important information to investigations of energy and mineral resources, environmental geochemistry, paleoclimates, oceanography (past and present), sedimentary, igneous, and metamorphic processes, Earth systems, geomicrobiology, and hydrology. One of the most important aspects of understanding and interpreting the stable isotope characteristics of sulfate minerals is the complex interplay between equilibrium and kinetic chemical and isotopic processes. With few exceptions, sulfate minerals are precipitated from water or have extensively interacted with water at some time in their history. Because of this nearly ubiquitous association with water, the kinetics of isotopic exchange reactions among dissolved species and solids are fundamental in dictating the isotopic composition of sulfate minerals. In general, the heavier isotope of sulfur is enriched in the higher oxidation state, such that under equilibrium conditions, sulfate minerals (e.g. barite, anhydrite) are expected to be enriched in the heavy isotope relative to disulfide minerals (e.g. pyrite, marcasite), which in turn are expected to be enriched relative to monosulfide minerals (e.g. pyrrhotite, sphalerite, galena) ( Sakai 1968 , Bachinski 1969 ). The kinetics of isotopic exchange among minerals with sulfur at the same oxidation state, such as sphalerite, and galena, are such that equilibrium is commonly observed. In contrast, isotopic equilibrium for exchange reactions between minerals of different oxidation states depends on factors such as the pH, time and temperature of reaction, the direction of reaction, fluid composition, and the presence or absence of catalysts ( Ohmoto and Lasaga 1982 ). The kinetics of oxygen isotope exchange between dissolved sulfate and water are extremely sluggish. Extrapolation of the high-temperature (100 to 300°C) isotopic exchange kinetic data of Chiba and Sakai (1985) to ambient temperatures suggests that it would take several billions of years for dissolved sulfate and seawater to reach oxygen isotopic equilibrium. In contrast, the residence time of sulfate in the oceans is only 7.9 million years ( Holland 1978 ). However, at higher temperatures (>200°C), oxygen isotopic exchange is sufficiently rapid to permit application of sulfate isotope geothermometry to geothermal systems and hydrothermal mineral deposits. In general, equilibrium prevails at low pH and high temperatures, whereas kinetic factors preclude equilibrium at low temperatures even at low pH. Thus, the sluggish kinetics of sulfur and oxygen isotope exchange reaction at low temperatures impair the use of these isotopes to understand the conditions of formation of sulfate minerals in these environments. However, because of these slow kinetics, the oxygen and sulfur isotopic compositions of sulfate minerals may preserve a record of the sources and processes that initially produced the dissolved sulfate, because the isotope ratios may not re-equilibrate during fluid transport and mineral precipitation. The first part of this chapter is designed to provide the reader with a basic understanding of the principles that form the foundations of stable isotope geochemistry. Next, an overview of analytical methods used to determine the stable isotope composition of sulfate minerals is presented. This overview is followed by a discussion of geochemical processes that determine the stable isotope characteristics of sulfate minerals and related compounds. The chapter then concludes with an examination of the stable isotope systematics of sulfate minerals in a variety of geochemical environments.

Reviews in Mineralogy and Geochemistry↗
USGS2012

Advances in the simulation and automated measurement of well-sorted granular material: 2. Direct measures of particle properties

1. In this, the second of a pair of papers on the structure of well-sorted natural granular material (sediment), new methods are described for automated measurements from images of sediment, of: 1) particle-size standard deviation (arithmetic sorting) with and without apparent void fraction; and 2) mean particle size in material with void fraction. A variety of simulations of granular material are used for testing purposes, in addition to images of natural sediment. Simulations are also used to establish that the effects on automated particle sizing of grains visible through the interstices of the grains at the very surface of a granular material continue to a depth of approximately 4 grain diameters and that this is independent of mean particle size. Ensemble root-mean squared error between observed and estimated arithmetic sorting coefficients for 262 images of natural silts, sands and gravels (drawn from 8 populations) is 31%, which reduces to 27% if adjusted for bias (slope correction between observed and estimated values). These methods allow non-intrusive and fully automated measurements of surfaces of unconsolidated granular material. With no tunable parameters or empirically derived coefficients, they should be broadly universal in appropriate applications. However, empirical corrections may need to be applied for the most accurate results. Finally, analytical formulas are derived for the one-step pore-particle transition probability matrix, estimated from the image's autocorrelogram, from which void fraction of a section of granular material can be estimated directly. This model gives excellent predictions of bulk void fraction yet imperfect predictions of pore-particle transitions.

Journal of Geophysical Research F: Earth Surface↗
USGS2020

Shallow basin structure and attenuation are key to predicting long shaking duration in Los Angeles Basin

Ground motions in the Los Angeles Basin during large earthquakes are modulated by earthquake ruptures, path effects into the basin, basin effects, and local site response. We analyzed the direct effect of shallow basin structures on shaking duration at a period of 2–10 s in the Los Angeles region through modeling small magnitude, shallow, and deep earthquake pairs. The source depth modulates the basin response, particularly the shaking duration, and these features are a function of path effect and not site condition. Three-dimensional simulations using the CVM-S4.26.M01 velocity model show good fitting to the initial portion of the waveforms at periods of 5 s and longer but fail to predict the long shaking duration during shallow events, especially at periods less than 5 s. Simulations using CVM-H do not match the timing of the initial arrivals as well as CVM-S4.26.M01, and the strong late arrivals in the CVM-H simulation travel with an apparent velocity slower than observed. A higher-quality factor than traditionally assumed may produce synthetics with longer durations but is unable to accurately match the amplitude and phase. Beamforming analysis using dense array data further reveals the long duration surface waves have the same back azimuth as the direct arrivals and are generated at the basin edges, while the later coda waves are scattered from off-azimuth directions, potentially due to strong, sharp boundaries offshore. Improving the description of these shallow basin structures and attenuation model will enhance our capability to predict long-period ground motions in basins.

California↗
USGS2025

Timescales of surface faulting preservation in low-strain intraplate regions from landscape evolution modeling and the geomorphic and historical record

Large surface-rupturing intraplate earthquakes in stable continental regions (SCRs) are uncommon globally and have recurrence intervals of thousands to hundreds of thousands of years based on the paleoseismic and geomorphic record, challenging accurate active fault identification in these regions. To constrain the timescales of preservation for scarps created by surface ruptures from dip-slip earthquakes, we use a two-dimensional scarp diffusion model for typical intraplate settings and explore which parameters influence fault scarp preservation. These parameters include the coseismic vertical surface offset, the recurrence interval of similar magnitude earthquakes, diffusivity (as a proxy for mean annual precipitation rate), and the erodibility of the surficial material. We constrain parameter ranges from a compilation of historical surface ruptures in intraplate settings in a variety of climates, including the Central and Eastern United States, Australia, Europe, Central Asia (Mongolia, China), India, and West Africa. The timescales of scarp preservation from landscape evolution modeling agree well with observations of scarp preservation in low-strain SCR and intraplate tectonic settings, with some notable exceptions for Australian scarps. We find that the erodibility of the surficial material and earthquake recurrence interval have a stronger effect on the timescales of scarp preservation than diffusivity or coseismic vertical surface offset. Our model results may aid in identifying and characterizing subtle, slow-moving active faults in low-strain SCR and intraplate tectonic settings for different tectonic, geomorphic, and climatic characteristics. Accurate fault locations and characterization from the landscape record has implications for both probabilistic seismic and fault displacement hazard analyses.

Journal of Geophysical Research Solid Earth↗
USGS2023

Linking seed size and number to trait syndromes in trees

Aim Our understanding of the mechanisms that maintain forest diversity under changing climate can benefit from knowledge about traits that are closely linked to fitness. We tested whether the link between traits and seed number and seed size is consistent with two hypotheses, termed the leaf economics spectrum and the plant size syndrome, or whether reproduction represents an independent dimension related to a seed size–seed number trade-off. Location Most of the data come from Europe, North and Central America and East Asia. A minority of the data come from South America, Africa and Australia. Time period 1960–2022. Major taxa studied Trees. Methods We gathered 12 million observations of the number of seeds produced in 784 tree species. We estimated the number of seeds produced by individual trees and scaled it up to the species level. Next, we used principal components analysis and generalized joint attribute modelling (GJAM) to map seed number and size on the tree traits spectrum. Results Incorporating seed size and number into trait analysis while controlling for environment and phylogeny with GJAM exposes relationships in trees that might otherwise remain hidden. Production of the large total biomass of seeds [product of seed number and seed size; hereafter, species seed productivity (SSP)] is associated with high leaf area, low foliar nitrogen, low specific leaf area (SLA) and dense wood. Production of high seed numbers is associated with small seeds produced by nutrient-demanding species with softwood, small leaves and high SLA. Trait covariation is consistent with opposing strategies: one fast-growing, early successional, with high dispersal, and the other slow-growing, stress-tolerant, that recruit in shaded conditions. Main conclusions Earth system models currently assume that reproductive allocation is indifferent among plant functional types. Easily measurable seed size is a strong predictor of the seed number and species seed productivity. The connection of SSP with the functional traits can form the first basis of improved fecundity prediction across global forests.

Global Ecology and Biogeography↗