USGS ScienceSearch

SEARCH · USGS Science

Results for “One Ecosystem”

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.

1,216 records · Page 39Linked to original sources

New methodology for assessing underground natural gas storage resources – Example from Michigan Basin, United States

Energy consumption in the United States (U.S.) and across the world is shifting away from traditional fossil fuels like coal and oil, and towards natural gas and renewable sources, including hydrogen. Because gas demand is typically greatest during cold seasons and renewable sources sometimes produce variable supplies, it is important to store energy for use when demand exceeds supply. Whereas batteries and tanks typically store energy above ground, geologic (underground) storage may be able to retain much greater quantities of energy over much longer durations (e.g., Matos et al. 2019; Buursink et al. 2023). Consequently, the U.S. Geological Survey (USGS) is developing a methodology to locate new geologic energy storage sites with an initial focus on depleted hydrocarbon reservoirs followed by estimates of associated pore space or gas storage capacity. Depleted hydrocarbon fields in Michigan Basin in the U.S. already host large quantities of injected natural gas (about 681 billion cubic feet of seasonal gas capacity in 42 gas storage facilities; U.S. Energy Information Administration 2024). Therefore, we used this basin as an example to assess potential new storage resources alongside existing storage capacity and carbon dioxide sequestration resources (e.g., Katz and Coats 1968; Haagsma et al. 2020). In Michigan Basin gas storage is taking place in reservoirs within several key hydrocarbon producing formations with different lithologies, including the informal Stray sandstone of the Michigan Formation and the pinnacle reefs of the Niagara Group (Gautier et al. 1995; Swezey et al. 2015). The assessment methodology introduced here consists of two steps relying on both well information and reservoir data.

Michigan

Hydrogeologic atlas of aquifers in Indiana

Aquifers in 12 water-management basins of Indiana are identified in a series of 104 hydrogeologic sections and 12 maps that show the thickness and configuration of aquifers. The vertical distribution of water-bearing units and a generalized potentiometric profile are shown along 3,500 miles of section lines that were constructed from drillers' logs of more than 4,200 wells. The horizontal scale of the sections is 1:125,000. Maps of aquifers showing the areal distribution of each aquifer type were drawn at a scale of 1:500,000. Unconsolidated aquifers are the most widely used aquifers in Indiana and include surficial, buried, and discontinuous layers of sand and gravel. Most of the surficial sand and gravel is in large outwash plains in northern Indiana and along the major rivers. Buried sand and gravel aquifers are interbedded with till deposits in much of the northern two-thirds of Indiana. Discontinuous sand and gravel deposits are present as isolated lenses, primarily in glaciated areas. The bedrock aquifers generally have lower yields than most of the sand and gravel aquifers; however, bedrock aquifers are areally widespread and are an important source of water. Bedrock aquifer types consist of carbonates; sandstones; complexly interbedded sandstones, siltstones, shales, limestones, and coals; and an upper weathered zone in low permeability rock. Carbonate aquifers underlie about one-half of Indiana and are the most productive of the bedrock aquifers. The other principal bedrock aquifer type, sandstone, underlies large areas in the southwestern one-fifth of Indiana. No aquifer is known to be present in the southeastern corner of Indiana.

Water-Resources Investigations Report

Evaluation of models for estimating hydraulic conductivity in glacial aquifers from NMR logging

Nuclear magnetic resonance (NMR) logging is a promising method for estimating hydraulic conductivity ( K ). During the past ∼60 years, NMR logging has been used for petroleum applications, and different models have been developed for deriving estimates of permeability. These models involve calibration parameters whose values were determined through decades of research on sandstones and carbonates. We assessed the use of five models to derive estimates of K in glacial aquifers from NMR logging data acquired in two wells at each of two field sites in central Wisconsin, USA. Measurements of K , obtained with a direct push permeameter (DPP), K DPP , were used to obtain the calibration parameters in the Schlumberger-Doll Research, Seevers, Timur-Coates, Kozeny-Godefroy, and sum-of-echoes (SOE) models so as to predict K from the NMR data; and were also used to assess the ability of the models to predict K DPP . We obtained four well-scale calibration parameter values for each model using the NMR and DPP measurements in each well; and one study-scale parameter value for each model by using all data. The SOE model achieved an agreement with K DPP that matched or exceeded that of the other models. The Timur-Coates estimates of K were found to be substantially different from K DPP . Although the well-scale parameter values for the Schlumberger-Doll, Seevers, and SOE models were found to vary by less than a factor of 2, more research is needed to confirm their general applicability so that site-specific calibration is not required to obtain accurate estimates of K from NMR logging data.

Wisconsin

Streams, springs, and volcanic lakes for volcano monitoring

Introduction Volcanic unrest can trigger appreciable change to surface waters such as streams, springs, and volcanic lakes. Magma degassing produces gases and soluble salts that are absorbed into groundwater that feeds streams and lakes. As magma ascends, the amount of heat and degassing will increase, and so will any related geochemical and thermal signal. Subsurface magma movement can cause pressurization that alters hydrostatic head and may induce groundwater discharge. Fluid-pressure changes have been linked to distal volcano-tectonic earthquakes (White and McCausland, 2016; Coulon and others, 2017) and phreatic eruptions (for example, Yamaoka and others, 2016). Clearly, changes in groundwater and surface waters are both indicators of unrest and clues to how and where magma is rising toward the surface. Where possible, it is prudent to incorporate real-time hydrologic data into multiparameter monitoring of restless volcanoes. Hydrologic dynamics can also be tracked by changes in groundwater levels that are commonly measured in shallow boreholes (see of Flinders, A.F., Lowenstern, J.B., Coombs, M.L., and Poland, M.P., eds., Recommended capabilities and instrumentation for volcano monitoring in the United States: U.S. Geological Survey Scientific Investigations Report 2024–5062–K, 5 p., https://doi.org/10.3133/sir20245062k. ">chapter K , this volume, on boreholes; Hurwitz and Lowenstern, 2024). Although inferred to be common, relatively few volcano-hydrology anomalies are well documented, and many are essentially anecdotal (Newhall and others, 2001), reflecting the fact that high-resolution time series remain rare. Extreme examples include the 2008 eruption of Nevado del Huila, Colombia, where relatively minor phreatomagmatic eruptions were accompanied by expulsion of as much as 300 million cubic meters of groundwater from fissures high on the volcano (Worni and others, 2011), generating large lahars. Substantial decreases in flow rate from springs about 8 kilometers from the summit of Mayon Volcano, Philippines, have been noted before most eruptions in the 20th century (Newhall and others, 2001). Stream monitoring at Redoubt Volcano in 2009 allowed Werner and others (2012) to recognize that groundwater was unable to absorb (or scrub) the high flux of volcanic gas and that a high CO 2 /SO 2 precursor signal had been evident for 5 months prior to the eruption. A key to better interpreting hydrologic anomalies—or even identifying them—is therefore obtaining adequate baseline data. Most hydrologic monitoring at U.S. volcanoes has been accomplished by intermittent sampling surveys with annual or less frequent sampling (for example, https://hotspringchem.wr.usgs.gov/index.php ). More frequent sampling, however, generally is needed to establish reliable baselines. A recent hydrologic and hydrothermal monitoring experiment at 25 sites and 10 of the 12 level 4 (very high threat) volcanoes in the U.S. portion of the Cascade Range demonstrated that there is sufficient temporal variability in hydrothermal fluxes, even during quiescent periods, that one-time measurements will commonly have limited interpretive value (Crankshaw and others, 2018). Thus, surveys are best augmented with data from streamgages (for example, Evans and others, 2004; Bergfeld and others, 2008). Streamflow (water discharge) data allow measured temperature and specific conductance to be converted to heat and solute mass fluxes, which could be insightful parameters for detecting anomalous activity (McCleskey and others, 2012). At the Yellowstone Caldera, long-term monitoring of river solutes has allowed calculation of the chloride flux, a proxy for heat discharge (Hurwitz and others, 2007; McCleskey and others, 2016) from the subsurface magma. This is readily accomplished because data from streamgages are continuously recorded and archived by the U.S. Geological Survey (USGS) National Water Information System (NWIS) (USGS, 2024). Similar studies on stratovolcanoes or shield volcanoes would be scientifically useful, and yet are logistically challenging, requiring streamgages on numerous radial drainages complemented by either frequent manual sampling or numerous deployments of equipment to measure water temperature and specific conductance as a proxy for water chemistry. Another challenge is that some volcanic areas, especially shield volcanoes, are characterized by near-surface porous rocks and soils, such that surface streams are rare and replaced by distant, dilute large-volume springs with only a trace of any original volcanically sourced water (Manga, 2001; Hurwitz and others, 2021). Volcanic lakes are worthy of special attention for monitoring efforts, as their temperature and composition can provide evidence of increased flux of volatile-rich fluids from below. Quantifying changes in volatile and heat release from magma can be simpler in lakes than for volcanoes with radial drainages and no major lakes. Moreover, volcanic lakes pose a range of hazards themselves, including phreatomagmatic eruptions, debris flows, flank collapse, tsunamis, and toxic gas release (Mastin and Witter, 2000; Delmelle and others, 2015; Manville, 2015; Rouwet and others, 2015)—hazards that have historically been responsible for substantial loss of life at many volcanoes worldwide (Manville, 2015). Catastrophic CO 2 release at Lake Nyos, Cameroon, in 1986 suffocated about 1,750 people and about 3,500 livestock and was probably triggered by a large landslide into the gas-saturated lake (Kling and others, 1987; Evans and others, 1993). Gas-charged springs in Soda Bay within Clear Lake (California) have caused almost a dozen deaths to bathers in the past hundred years (ABC News, 2000). A 2005 example of lake overturn and abundant gas release was documented at Mount Chiginagak in Alaska (Schaefer and others, 2008) but did not result in any human casualties. Although thermally stratified lakes, which promote trapping of exsolved magmatic gas, tend to develop in tropical regions, the phenomenon can also arise where salinity creates meromixis (a condition in which a lake does not mix completely), as occurs in Mono Lake, California (Jellison and Melack, 1993; Jellison and others, 1998). If magma erupts or flows into a lake, the interaction between hot magma and cold water can be explosive (Mastin and others, 2004; Zimanowski and others, 2015) and substantially expand the area affected by the eruption. Another hazard is the breaching of crater rims by landslides triggered by volcanic and (or) seismic activity. Under some circumstances, substantial volumes of water can be displaced, leading to large floods and lahars. Late Holocene lake flooding from Aniakchak Crater in the Alaska Peninsula (Waythomas, 2022) and from Paulina Lake in Newberry Crater, Oregon (Chitwood and Jensen, 2000), caused by the failure of outlet sills, testify to the substantial hazards at lake-filled calderas. Several volcanic systems in the United States host lakes known to receive heat and gas from underlying magma. These lakes vary widely in area, depth, and chemical composition. Lakes are present at level 4 volcanoes, including Crater Lake and Newberry Volcano in Oregon; Yellowstone Caldera in Wyoming; Long Valley Caldera, Clear Lake volcanic field, Medicine Lake, and Salton Buttes in California; and Aniakchak Crater, Mount Katmai, Fisher Caldera, Mount Okmok, and Kaguyak Crater, among others, in Alaska. A water lake was present in Halemaʻumaʻu, the crater of Kīlauea, Hawai‘i (fig. F1), from October 2019 to December 2020. Level 3 volcanoes with lakes include Mono Lake volcanic field (Calif.), Mount Bachelor (Ore.), Ukinrek Maars and Mount Chiginagak (Alaska), and Soda Lake (Nevada). In addition, there are lakes at many levels 1 and 2 volcanoes. In the United States, there are no strongly acidic lakes that receive abundant input of magmatic gas, such as those found at Mount Ruapehu (New Zealand), Ijen and Kelud (Indonesia), and Poás (Costa Rica). Nevertheless, many contain fluids that provide clues to magmatic processes below. Since publication of a previous report on recommended instrumentation for volcano monitoring (Moran and others, 2008), continuous hydrologic monitoring has become increasingly feasible. However, changes in water pressure, temperature, and chemistry remain, in general, poorly studied phenomena at volcanoes (Sparks, 2003; National Academies of Sciences, Engineering, and Medicine, 2017). Recent efforts by the USGS have included the temporary study of Cascade Range volcanoes, which included frequent (15 minute to hourly) temporal sampling of temperature, depth, and conductivity (Crankshaw and others, 2018; Ingebritsen and Evans, 2019). At Yellowstone Caldera, many streamgages have now added thermistors and specific conductance sensors, allowing estimation of time-dependent chloride flux as a proxy for variations in subsurface heat flux (McCleskey and others, 2012, 2016). Efforts to better understand lakes have also accelerated, with bathymetric mapping and sampling carried out at several locations in the United States. Especially thorough work was done at Yellowstone Lake thanks to the Hydrothermal Dynamics of Yellowstone Lake (HD-YLAKE, https://hdylake.org ) project, funded primarily by the National Science Foundation. In addition to geophysical surveys and recovery of cores and other samples, HD-YLAKE investigations included remotely operated vehicle (ROV) investigations of hydrothermal vents on the lake floor (fig. F2). Data collected by the ROV provided a better understanding of the thermal and chemical influx from lake-bottom hydrothermal systems (Sohn and others, 2017). In this chapter, we focus on detecting changes in the chemistry, temperature, discharge, or water levels of streams, springs, and lakes that can be caused by seismicity, volumetric strains, or increases in gas flux associated with ascending magma. There is unavoidable overlap with other chapters of this report. Samples of water and gas can also be obtained in boreholes ( of Flinders, A.F., Lowenstern, J.B., Coombs, M.L., and Poland, M.P., eds., Recommended capabilities and instrumentation for volcano monitoring in the United States: U.S. Geological Survey Scientific Investigations Report 2024–5062–K, 5 p., https://doi.org/10.3133/sir20245062k. ">chapter K , this volume; Hurwitz and Lowenstern, 2024), both shallow and deep. Gas monitoring ( of Flinders, A.F., Lowenstern, J.B., Coombs, M.L., and Poland, M.P., eds., Recommended capabilities and instrumentation for volcano monitoring in the United States: U.S. Geological Survey Scientific Investigations Report 2024–5062–E, 11 p., https://doi.org/10.3133/sir20245062e.">chapter E , this volume; Lewicki and others, 2024) relies in part on samples from springs and wells, particularly where measurable gas plumes are absent. Water acts as a trigger and lubricant for landslides and sediment-rich floods, and so hydrology has obvious relevance for lahar monitoring, as discussed in of Flinders, A.F., Lowenstern, J.B., Coombs, M.L., and Poland, M.P., eds., Recommended capabilities and instrumentation for volcano monitoring in the United States: U.S. Geological Survey Scientific Investigations Report 2024–5062–H, 6 p., https://doi.org/10.3133/sir20245062h. ">chapter H (this volume; Thelen and others, 2024). Shared situational awareness among scientists engaged in geophysical, gas, and hydrologic monitoring will improve overall understanding of the volcanic hazard.

Scientific Investigations Report

Ductile and brittle Rio Grande Rift deformation in Oligocene granite records a two-stage rift history in southern Colorado

The timing and nature of early deformation in the Rio Grande Rift remains poorly constrained. We present evidence for the earliest structural signature of rift extension in the Sangre de Cristo Range, southern Colorado, based on new geologic mapping, structural analysis, rock magnetic data, and thermochronology. These analyses focus on the ~30.0 Ma granite of Chokecherry Canyon, which hosts discrete low-angle mylonitic shear zones and a distributed, gently SW-dipping protomylonitic fabric. Incremental stretching axes, stretching lineations, and Kmax magnetic lineations plunge gently WSW. Quartz microstructures and crystallographic orientations indicate dominantly coaxial strain in the protomylonite and general shear in the discrete shear zones. Quartz c-axis opening-angle thermometry suggests deformation at ~420–540°C. Thermal modeling of ⁴⁰Ar/³⁹Ar K-feldspar data indicates rapid post magmatic cooling below the brittle–plastic transition, supporting shear-zone formation immediately after emplacement. Slow cooling from ~20–13 Ma was followed by renewed rapid cooling at ~13 Ma, interpreted as the onset of extensional exhumation along the Sangre de Cristo Fault System. These results show that extension in the northern Rio Grande Rift was active by ~30 Ma, earlier than previously recognized. We propose a two-stage model for northern Rio Grande Rift evolution: Stage I (30–23 Ma) records ENE–WSW extension localized in low-angle mylonitic shear zones associated with mid-crustal intrusions; Stage II (≤18 Ma) reflects brittle high-angle normal faulting, focused exhumation, and rift narrowing. Stage I magmatism and deformation along the western range front likely established crustal weaknesses that guided later fault development.

Colorado

A decade of shaking in the Garden City: The dynamics of preparedness, perceptions, and beliefs in Canterbury, New Zealand, and implications for earthquake information

This study explored earthquake preparedness over time - before, during, and 10 years after the Canterbury Earthquake Sequence (CES) in Aotearoa New Zealand (NZ). Surveys of Canterbury residents were conducted in 2009, 2013 and again in 2021, using variables derived from Community Engagement Theory (CET). The surveys measured earthquake perceptions and beliefs, participation and engagement, and preparedness actions. Results were compared across the three samples. Findings indicate that perceptions and beliefs (e.g. risk perception, outcome expectancy beliefs), and types of preparedness actions taken (e.g. collection of survival items, structural preparedness, community and agency relationships), differed over time, depending on people’s experiences before, during, and after the CES. For example, during and after the CES people were more likely to believe that preparing provided a benefit to daily life, but less likely to think it could reduce property damage, perhaps due to people’s experiences of disruption and damage during the earthquakes. An understanding of such dynamics can assist with the provision and timing of risk and preparedness information. This study highlights the importance of providing applicable and actionable preparedness information, that is relevant to people’s experiences, throughout an earthquake sequence. Such information might evolve and change in focus over time depending on risks and needs. Focus could also be given to information that builds peoples beliefs and capacities to undertake preparedness in evolving situations. Understanding preparedness in the context of different experiences and timeframes is useful in helping update models such as the CET, where the dynamics of time might be better incorporated.

Canterbury

Earthquake ground-motion model adjustments for the San Francisco Bay area

We develop adjustments to ergodic ground‐motion models (GMMs) to improve their performance in the San Francisco Bay Area (SFBA). GMMs are widely used in hazard assessments to estimate characteristics of ground shaking based on known properties of the source, path, and site. Such models are often developed using datasets containing records from various regions, resulting in models that represent median ground‐motion behavior, which may not adequately represent ground motions within subregions. This is true for the SFBA, where ground motions attenuate more rapidly with distance than in many other parts of California that dominate GMM databases. To support improved seismic hazard estimates in the SFBA, we calculate regional constants and anelastic attenuation coefficient adjustments relative to two commonly used ergodic GMMs: BSSA14 ( Boore et al. , 2014 ) and ASK14 ( Abrahamson et al. , 2014 ). These adjustments are obtained for a suite of ground‐motion intensity measures (peak ground acceleration, peak ground velocity, and 5%‐damped pseudospectral acceleration at oscillator periods ranging from 0.075 to 10 s) using mixed‐effects regression. Use of the regionally adjusted models reduces the overall bias by up to 0.5 natural log units for BSSA14 and up to 0.6 natural log units for ASK14. We demonstrate one application of our attenuation adjustments and their implications in an earthquake early warning case study of the 2014 M 6.0 South Napa earthquake. The predicted extent of shaking using the adjusted models better matches observed shaking at large source‐to‐site distances, especially for lower shaking intensities, thus potentially reducing overalerting. We encourage the use of our model adjustments when ergodic models are considered for seismic hazard studies in the SFBA.

California

Localization of spatiotemporally heterogeneous subsurface flows using autoencoder-based deep learning framework for time-lapse self-potential tomography

Self-potential (SP) monitoring has emerged as a valuable method for characterizing subsurface hydrogeological features and processes due to its sensitivity to fluid-induced electrokinetic effects. Despite advancements in SP inversion, challenges remain in imaging groundwater dynamics from SP activities due to complex hydrological settings and transient noise. In this study, a deep learning autoencoder (AE)-based framework is proposed for the spatiotemporal localization of subsurface fluid movement from time-lapse SP tomography. Temporal segments of time-lapse numerical inversions were first derived from long-term SP monitoring conducted from a floodplain site in Oak Ridge, Tennessee, known for active hyporheic exchange. Subsequently, AE models based on vision transformer (ViT), convolutional long short-term memory (ConvLSTM), convolutional neural network, and temporal convolutional network were individually trained and compared on the SP tomography segments for reconstruction performance. Finally, the reconstruction error over time serves as an anomaly score to identify moments of active SP variation, whereas spatial distributions of errors within these moments are analyzed to image and localize regions associated with anomalous subsurface fluid movement. The results demonstrate that ConvLSTM- and ViT-AE are most capable for the localization task with contrasting error distributions and consistent delineation of anomalies. Applying the method to both SP arrays parallel and perpendicular to the stream produced consistent anomaly zones near a fault or karst feature, validating the robustness and generalization of the approach. These results demonstrate the potential of the proposed framework as a scalable and interpretable tool for spatiotemporal analysis of subsurface flow dynamics in complex hydrogeological systems.

Tennessee

U-Pb scheelite ages of tungsten and antimony mineralization in the Stibnite-Yellow Pine district, central Idaho

The Stibnite-Yellow Pine district contains the largest antimony resource in the United States, as well as significant gold, and is a historic producer of tungsten. Application of in situ laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS) direct dating of scheelite from two Au-Sb-W ore deposits, Yellow Pine and Hangar Flats, yielded an older group of U-Pb ages in the range of 60.0 ± 2.8 to 57.0 ± 1.1 Ma and a younger U-Pb age for scheelite intergrown with stibnite of 47.4 ± 1.1 Ma. These in situ analyses were calibrated by isotope dilution-thermal ionization mass spectrometry (ID-TIMS) U-Pb lower intercept ages of two coarsely crystalline scheelite samples that yielded ages of 57.52 ± 0.22 and 56.62 ± 0.16 Ma. Scheelite of the latter age is of sufficient quality to serve as a primary reference material for LA-ICP-MS scheelite U-Pb geochronology. The group of older U-Pb scheelite ages agrees with 40 Ar/ 39 Ar ages of 56.9 ± 1.2 to 56.38 ± 0.54 Ma on adularia from Yellow Pine and Hangar Flats, whereas the younger U-Pb scheelite age is similar to an 40 Ar/ 39 Ar age of 46.00 ± 0.40 Ma on adularia from an epithermal gold-silver deposit in the adjacent Thunder Mountain caldera. Our results indicate that the main stage of tungsten mineralization occurred at ca. 57 Ma, whereas the main stage of antimony mineralization occurred at ca. 47 Ma—thereby providing first-time age constraints for antimony and tungsten mineralization in the Stibnite-Yellow Pine district.

Idaho

Sensitivity of Alabama freshwater gastropod species to nickel exposure

Snails are effective bioindicators due to their prolific distribution, high level of endemism, and capacity to accumulate contaminants. Freshwater snails have unique ecological niches which are imperiled by land‐use change and the introduction of hazardous chemicals. To assess how environmental alterations affect gastropods, lab‐based studies are needed to characterize the toxicity of specific stressors. This can help guide policy decisions and remediation efforts. The aim of this research was to characterize acute toxicity of nickel (Ni) on endemic snails ( Somatogyrus georgianus [Walker, 1904], Elimia cahawbensis [Lea, 1861], and Elimia spp.) and measure the accumulation of Ni and mineral elements including calcium (Ca), magnesium, potassium, and sodium (Na). Snails were exposed to six concentrations (25–800 µg/L) of Ni for 96 h. Among the studied snail species, E. cahawbensis was the most sensitive to Ni, with the lowest lethal concentration where 50% of the organisms died (LC50) at 88.88 µg/L Ni after 96 h. The LC50 at 96 h for S. georgianus was 167.78 µg/L Ni, and 393.13 μg/L Ni for Elimia spp. Except for Elimia spp., mortality of the other two snail species corresponded to the whole‐body uptake of Ni. Nickel exposure also influenced Ca and Na uptake for Elimia spp. All three endemic species are potential candidate species for evaluating localized effects of human activities, and the present study provides a first step in characterizing how snails would be affected by environmental alterations. More research could further characterize potential effects of other human stressors on these endemic snail species. Future research into subindividual responses and routes of exposure can further elucidate variations in species sensitivity.

Alabama

Nanometer-scale relationships between sedimentary organic matter molecular composition, fluorescence, cathodoluminescence, and reflectance: The importance of oxygen content at low thermal maturities

Molecular characterization of sedimentary organic matter (SOM), termed macerals, is a common goal when seeking to understand petroleum generation as well as other geologic processes in deep time. However, unambiguous measurement of discrete macerals is challenging due to the small size of organic particles in sedimentary rocks, the proximity of different organic matter types to one another, mineral-organic matter interactions, and maceral mixing that occurs during SOM isolation prior to ex situ analysis. The recent advent of infrared spectrometers capable of nanometer-scale resolution and the application of these technologies to geologic samples has enabled advances in rapid, in situ molecular characterization of SOM allowing for insights into paleoenvironmental processes, such as organic matter productivity and preservation, among others. Here we employ one such technology, optical photothermal infrared (OPTIR) spectroscopy, to map SOM functional group distributions at 500-nm resolution in a sample from the Lower Cretaceous Sunniland Limestone of the South Florida Basin. Examined fields of view include occurrences of amorphous organic matter (AOM), inertinite, micrinite, solid bitumen, telalginite, and vitrinite. OPTIR data from these macerals are compared against traditional organic petrographic data from the same organic grains including fluorescence intensity and white light reflectance as well as against cathodoluminescence response, an emerging organic petrographic approach. Maceral oxygen content (using carbonyl functional group abundance as a proxy) is observed to vary widely between maceral types but correlates strongly with fluorescence and cathodoluminescence intensity as well as against reflectance. These findings highlight the important role that oxygen content plays in determining the optical properties of SOM and further demonstrate the ability of OPTIR to discriminate subtle molecular differences between SOM types.

Organic Geochemistry

Waning greenhouse gas emissions from U.S. Federal lease coal production by the mid-21st century

This study presents estimates of future years (2024–2051) United States Federal lease coal production and the resulting greenhouse gas (GHG) emissions from the combustion, transport, and mining of that fuel. Results from the coal production estimate indicate a decline in production from Federal leases; with known production of 240 million short tons (mtn) in 2023 and a projected decline to 34.0 mtn by 2051, which represents a reduction to 14.2% of the 2023 value. In parallel with this projection, total GHG emissions are estimated to decrease from 402.2 million metric tons of carbon dioxide equivalent (MMT CO 2 eq.) in 2024 to 55.0 MMT CO 2 eq. in 2051, a decline to 13.7% of 2024 emissions estimates. The reductions in coal production and emissions are mainly the result of planned coal combustion power plant closures, with major projected closures in 2037 and 2048. However, GHG emissions estimates for future years can be uncertain as they rely heavily on coal production estimates from operators' public business plans and other publicly available resources. Forward looking plans of this type are subject to significant changes if economic and political factors deviate from current information. Results suggest that average GHG emissions over the time series breakout to 95% end point combustion, 3.7% transportation combustion emissions, and 1.3% fugitive emissions, although there is uncertainty associated with these figures. Uncertainty stemming from production projections, sector distributions, and emissions factors on the future emissions estimates increases with time, ranging from −28% to +48% within the 2024–2051 timeframe.

Alabama, Colorado, Montana, North Dakota, Utah, Wy

Numerical modeling of Late Pleistocene to Holocene earthquake-induced progressive rock slope damage in the paraglacial Serpentine valley, Prince William Sound, Alaska

Landslides in deglaciating fjords pose a potential tsunami threat to nearby communities; however, processes contributing to long-term progressive rock damage and landslide conditioning remain poorly constrained in paraglacial settings. Here, we analyse the role of earthquake-induced rock mass damage over late Pleistocene to Holocene time scales, as a conditioning factor for modern landslides, using distinct element numerical modelling to assess spatial and temporal patterns of fracture propagation influenced by varying glacier thickness. Conceptualised numerical models were parameterised by in situ rock mass, glacial and topographic conditions in Serpentine valley, located in Prince William Sound, Alaska, where several large landslides are actively developing along the western valley wall. Results show that although rigid glacier buttressing reduces co-seismic rock mass damage, it does not suppress it completely, and damage occurs both above and below the glacier surface elevation. We demonstrate that topography, and especially steep slopes with topographic convexity, as well as preexisting damage of joint networks and faults inherited from tectonic and exhumation induced stresses, exert primary control on the location of new co-seismic damage. Simulations representing a simplified deglaciation sequence over the past 25 ka, with a series of 10 evenly spaced earthquakes, generated rock mass damage patterns that qualitatively match in situ landslide structural and kinematic observations at one instability in Serpentine valley. Our conceptual study helps clarify the role of repeated seismicity over glacial timescales as a long-term conditioning process for paraglacial rock slope failure, highlighting spatial and temporal patterns of progressive damage accumulation, with outcomes relevant for modern landslide hazard assessment.

Alaska

Comparisons of shoreline positions from satellite-derived and traditional field- and remote-sensing techniques

Satellite-derived shorelines (SDS) have the potential to help researchers answer critical coastal science questions and support work to predict coastal change by filling in the spatial and temporal gaps present in current field-based and remote-sensing data collection methods. The U.S. Geological Survey conducted comparison analyses of traditionally sourced shorelines and SDS in diverse coastal landscapes to determine how SDS could be used in ongoing and future work across varied coastal environments and provided some initial findings that could be used for implementation. Using CoastSeg, a browser-based program for SDS detection and mapping, SDS for the period 1984–2023 for multiple locations across the United States were compared to shoreline positions from traditionally sourced shoreline data. In this report, the authors present these comparisons alongside lessons learned and challenges encountered when building SDS workflows in different coastal locations. Results show that individual SDS have larger uncertainty and yet produced similar linear trends to sparser, traditionally sourced shoreline data; because SDS methods provide orders of magnitude more data than traditional shoreline-detection methods, they can be used to evaluate shoreline behaviors. Refining average scalar slopes used in tidal corrections did not result in substantial decreases in uncertainty. Using lessons from this work to outline needs for regional implementation, initial setup time would be considerable, being on the order of weeks. However, once complete, shoreline detections and analyses are fast (on the order of minutes to hours) and achievable using a desktop computer.

Alaska, Florida, Massachusetts, Washington

Design, revisions, and considerations for continued use of a ground-water-flow model of the coastal plain aquifer system in Virginia

A digital numerical model of the ground-water-flow system in the Virginia Coastal Plain was constructed as part of the Regional Aquifer System Analysis program of the U.S. Geological Survey and subsequently revised based on reinterpretations of the hydrogeologic framework. The revised model has been incorporated by the Virginia Department of Environmental Quality as a means to evaluate the effects of existing and proposed ground-water withdrawals. The revised model, which is documented here, includes many design features of the original model. Principal modifications are (1) reassignment of some model layers to represent different aquifers, (2) replacement of applied recharge and the uppermost surface-water layer by constant-head cells to represent the unconfined aquifer, and (3) revision of the lateral extents and transmissivity and vertical leakance values assigned to model layers. Ground-water levels simulated by the revised model under steady-state conditions, using withdrawal rates representing the period 1978-80, generally are similar to the water levels simulated by the original model under transient conditions using the same withdrawal rates. Differences probably result largely from changes in transmissivity and vertical leakance values. T he revised model enables the continued evaluation of development of the large and complex aquifer system. Improved information on parts of some aquifers has been incorporated. Only large, regional-scale trends in ground-water water levels, flow directions, and flow rates are simulated, however, and not local trends and (or) short-term changes in water level and flow. Changing hydraulic stresses, knowledge of hydrogeologic conditions, and resource-management needs require periodic revision of the model to maintain its usefulness for the management of ground-water resources in the Virginia Coastal Plain.

Virginia

Gas emissions from the Sulphur Bank Mercury Mine hydrothermal system, Clear Lake volcanic field, California

The Sulphur Bank Mercury Mine (SBMM) hydrothermal system offers insights into active degassing processes in the Clear Lake volcanic field (CLVF), a high-threat region based on its record of Holocene eruptions and proximity to populated areas. Here we present chemical and isotopic analyses of gas samples collected between 2015 and 2023, along with the first comprehensive CO 2 flux survey of the SBMM area conducted in 2023. Sampled gases are CO 2 - and CH 4 -rich (≥84 and 6 mol% in dry gas, respectively) with high mantle-derived helium contributions ( 3 He/ 4 He = 6.54–7.86 R C /R A ). Carbon isotopic compositions of CO 2 (δ 13 C = −10.0 to −9.5 ‰) and CH 4 (δ 13 C = −35.8 ‰) indicate mixed sources, with significant contributions from metamorphism of organic-rich Franciscan Complex rocks hosting the hydrothermal system. Modeling of gas compositions shows that scrubbing by interaction with air-saturated groundwater strongly influences observed compositional variability. From our CO₂ flux measurements, we estimate the deeply derived CO 2 emission rate from the SBMM hydrothermal area (0.2 km 2 ) at 240 t d −1 , comparable to many quiescently degassing volcanoes worldwide. We also provide a first-order estimate of CH 4 emissions at approximately 0.5 t d −1 . Our findings establish crucial baseline data for future volcanic monitoring efforts, enhancing detection capabilities for potential changes in this active hydrothermal system. This work contributes to the broader understanding of volatile contributions from volcanic and metamorphic sources to the global carbon budget, while highlighting the strong influence of bedrock geology on gas compositions in the CLVF.

California

Techniques for simulating flood hydrographs and estimating flood volumes for ungaged basins in east and west Tennessee

A dimensionless hydrograph developed for a variety of basin conditions in Georgia was tested for its applicability to streams in East and West Tennessee by comparing it to a similar dimensionless hydrograph developed for streams in East and West Tennessee. Hydrographs of observed discharge at 83 streams in East Tennessee and 38 in West Tennessee were used in the study. Statistical analyses were performed by comparing simulated (or computed) hydrographs, derived by application of the Georgia dimensionless hydrograph, and dimensionless hydrographs developed from Tennessee data, with the observed hydrographs at 50 and 75% of their peak-flow widths. Results of the tests indicate that the Georgia dimensionless hydrography is virtually the same as the one developed for streams in East Tennessee, but that it is different from the dimensionless hydrograph developed for streams in West Tennessee. Because of the extensive testing of the Georgia dimensionless hydrograph, it was determined to be applicable for East Tennessee, whereas the dimensionless hydrograph developed from data on streams in West Tennessee was determined to be applicable in West Tennessee. As part of the dimensionless hydrograph development, an average lagtime in hours for each study basin, and the volume in inches of flood runoff for each flood event were computed. By use of multiple-regression analysis, equations were developed that relate basin lagtime to drainage area size, basin length, and percent impervious area. Similarly, flood volumes were related to drainage area size, peak discharge, and basin lagtime. These equations, along with the appropriate dimensionless hydrograph, can be used to estimate a typical (average) flood hydrograph and volume for recurrence-intervals up to 100 years at any ungaged site draining less than 50 sq mi in East and West Tennessee.

Tennessee

Critical review of mercury methylation and methylmercury demethylation rate constants in aquatic sediments for biogeochemical modeling

Mercury is a toxin that causes neurological impairments in adults, is particularly harmful for fetuses and children, and is deadly in severe cases, making it a worldwide health concern. Methylmercury (MeHg) is the environmentally relevant form of mercury (Hg) because it biomagnifies along the food chain. Methylmercury is mainly produced in aquatic sediments via methylation of inorganic Hg (Hg(II)) and transformed back via demethylation. Because transformation rates determine MeHg concentrations, quantification of methylation and demethylation rates is needed to inform management of MeHg. Published rate constants for Hg(II) methylation ( 𝑘 𝑚 ) and MeHg demethylation ( 𝑘 𝑑 ) vary greatly, stemming partly from differences in experimental methods. We conducted a comprehensive review of rate laws, evaluated published rate constants, and performed biogeochemical simulations to assess variability in reported 𝑘 𝑚 and 𝑘 𝑑 . Based on selected studies employing the same pseudo-first-order rate law and similar experimental methods, we found that 𝑘 𝑚 = 0.04 ± 0.03 d −1 is a reasonable range for wetland sediments. Over a number of environments, maximum 𝑘 𝑑 was smaller at sites without Hg source ( 𝑘 𝑑 = 0.5 d −1 ) than at sites with identified Hg source ( 𝑘 𝑑 = 1.8 d −1 ). Larger variability and higher uncertainty in 𝑘 𝑑 compared to 𝑘 𝑚 highlight the need for more research on MeHg demethylation rates. This critical review: (a) aids the design of future experimental studies of 𝑘 𝑚 and 𝑘 𝑑 ; (b) provides guidance for comparing rate constants from different studies; (c) presents a biogeochemical reaction model to assess rate constants; and (d) informs selection of 𝑘 𝑚 and 𝑘 𝑑 values from the literature for use in model simulations.

Critical Reviews in Environmental Science and Tech