USGS ScienceSearch

SEARCH · USGS Science

Results for “Frontiers in Ecology and Environment”

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,673 records · Page 32Linked to original sources

Prioritizing resource protection and understanding potential susceptibility of springs to surficial changes in a low-temperature geothermal system

Geothermal systems are vulnerable to changes in water budget and composition, requiring science-based management. This study uses a dataset of spring water temperatures, time series of groundwater residence time tracers (tritium and carbon-14), and stable isotopes of water to understand geothermal flow in a low-temperature geothermal system in north west Colorado, United States (Steamboat Springs). The geothermal system is bisected by the Yampa River, necessitating a stream mass balance approach to quantify total discharge. Time series analysis of water temperature data provides a ranked list of features more susceptible to surficial changes, which is corroborated using time series of tritium which indicate spatially distinct patterns of mixing between modern and pre-modern groundwater. All springs contain a portion of pre-modern groundwater that is thousands to tens of thousands of years old, a period coinciding with melting of extensive Pleistocene glaciers that was likely one of the recharge sources to the geothermal system. Stream mass balance indicates that greater than 80% of the total geothermal discharge is derived from diffuse or small springs, highlighting the extensive nature of the geothermal outflow zone and the association with local geologic structures. This study provides baseline data to support management of the Steamboat Springs geothermal system and indicates the utility of these approaches in developing science-based geothermal management.

Colorado

Recent scientific contributions by the U.S. Geological Survey in the San Francisco Bay and Sacramento–San Joaquin Delta Estuary

Introduction The San Francisco Bay and Sacramento–San Joaquin Delta Estuary (hereafter, Bay-Delta) is the largest estuary on the West Coast of the United States. The Bay-Delta covers more than 1,600 square miles and drains a watershed of more than 75,000 square miles, which is greater than 40 percent of California. The region surrounding the Bay-Delta is home to about 10 million people, and its habitats (fig. 1) support more than 800 plant and animal species. The waterways of the Bay-Delta are the central hub of California’s extensive freshwater delivery system, supplying water to more than 27 million Californians and 4 million acres of farmland in the Central Valley. The U.S. Geological Survey (USGS) is the primary science agency of the U.S. Department of the Interior, providing a broad range of Earth, water, biological, and mapping data and expertise to inform natural resource management across the country. This fact sheet focuses on research conducted by the USGS in the Bay-Delta region, mostly within the past 5 years. The fact sheet is organized across five major themes: water flow and water quality, fish and wildlife, wetland restoration, invasive species, and hazards. In each of these areas, the USGS works closely with Federal, State, and local agencies; academia; and non-governmental organizations to provide objective, evidence-based science. The data and knowledge gained from USGS research and monitoring are publicly available and are routinely used by agencies, including the U.S. Bureau of Reclamation, U.S. Army Corps of Engineers, U.S. Fish and Wildlife Service, California Department of Water Resources, California Department of Fish and Wildlife, California State Water Resources Control Board, Delta Stewardship Council, and many additional organizations and stakeholders to inform management of the Bay-Delta region and its resources.

California

Quaternary-active faults and the role of inherited structures in the Sacramento-San Joaquin Delta, western Central Valley, northern California

Seismic sources and their associated hazards within the Sacramento-San Joaquin Delta region of north-central California are relatively poorly characterized as compared to other, more heavily studied regions of northern California, such as the San Francisco Bay Area. Here we present a synthesis of subsurface, bedrock geology, and geodetic datasets from the Delta and from the Coast Ranges and Diablo Range to the northwest and southwest, respectively. We integrate these data and our own surface geologic and geomorphic observations to present a comprehensive review of faults in the Delta that exhibit Quaternary activity. Structural geologic data from the surrounding region highlight the significant influence that Late Cretaceous-to-Paleogene forearc structures exert on the geometry and kinematics of major Quaternary-active structures within the Delta. These inherited structures — including the Pittsburg-Kirby Hills Fault, Midland Fault, and Great Valley Fault System — exhibit a range of geometries and kinematics. Analysis of geomorphology along these structures suggests that these structures combine to accommodate Quaternary strain across the Delta region. A clearer understanding of subsurface geometries and structural relationships, built upon the regional tectonic history, provides insight into modern deformation accommodated on older structures and helps inform interpretations of seismic hazard within the Delta.

California

Petrogenesis and mineralization potential of spinifex komatiitic basalts in the Bradley Peak greenstone terrane, Wyoming Province

Komatiitic volcanic rocks are important hosts of Ni sulfide mineralization and record early Earth evolution; however, those in the well-studied Archean Wyoming Province have received little attention. Here, we elucidate the timing and petrogenesis of the Bradley Peak komatiitic volcanic rocks using field and textural observations, geochronology, and geochemistry. Detrital and igneous zircon U-Pb ages for two samples from previously undated units support published age determinations, placing the eruption age at 2.72 Ga. Stratigraphy of the volcanic flows was mapped and 36 samples including cumulates, greenschists, and spinifex-textured rocks were collected. Whole-rock geochemistry was used to classify the spinifex-textured samples as Al-undepleted komatiitic basalts (11–17 wt% MgO). Platinum-group element concentrations (n = 25) are like those in global Al-undepleted komatiitic basalts, and PGE/Ti ratios do not indicate the volcanic flows likely host sulfide mineralization. Initial ε Nd values of −0.5 to +4.7 (n = 16), indicate that these lavas were derived from a depleted mantle source and have negligible evolved crust contamination. The primary magma to the komatiitic basalt flows is estimated to have had 19 wt% MgO and be derived from ∼15 to 25 % mantle partial melting at 3–4 GPa. Trace element chemistry and thermodynamic modeling suggest the primary melt assimilated local banded iron formation. Although the Bradley Peak komatiitic basalts do not contain positive evidence of magmatic sulfide deposits, depleted Au in the flows suggests they could be source rocks for nearby orogenic gold deposits.

Wyoming

Interactive effects of salinity and hydrology on radial growth of bald cypress (Taxodium distichum (L.) Rich.) in coastal Louisiana, USA

Tidal freshwater forests are usually located at or above the level of mean high water. Some Louisiana coastal forests are below mean high water, especially bald cypress ( Taxodium distichum (L.) Rich.) forests because flooding has increased due to the combined effects of global sea level rise and local subsidence. In addition, constructed channels from the coast inland act as conduits for saltwater. As a result, saltwater intrusion affects the productivity of Louisiana’s coastal bald cypress forests. To study the long-term effects of hydrology and salinity on the health of these systems, we fitted dendrometer bands on selected trees to record basal area increment as a measure of growth in permanent forest productivity plots established within six bald cypress stands. Three stands were in freshwater sites with low salinity rooting zone groundwater (0.1–1.3 ppt), while the other three had higher salinity rooting zone groundwater (0.2–4.9 ppt). Water level was logged continuously, and salinity was measured monthly to quarterly on the surface and in groundwater wells. Higher groundwater salinity levels were related to decreased bald cypress radial growth, while higher freshwater flooding increased radial growth. With these data, coastal managers can model rates of bald cypress forest change as a function of salinity and flooding.

Louisiana

Identifying precursors and tracking pulses of magma ascent in multidisciplinary data during the 2018–2023 phreatomagmatic eruption at Semisopochnoi Island, Alaska

The 2018–2023 phreatomagmatic eruptions at Semisopochnoi Island, Alaska produced abundant long-period (LP) seismicity, harmonic and broadband tremor, and explosion signals over several well-monitored periods of eruption and quiescence. The corresponding dataset provides an excellent opportunity to investigate precursory and syn -eruptive geophysical signals of long-lived phreatomagmatic eruptions using multiparameter observations. We generated explosion and LP event catalogs through novel implementations of the REDPy ( Hotovec-Ellis, 2024 ) repeating event detector in mid-2021 following a network upgrade and the onset of a new phase of the eruption. The hundreds of detected explosions show a high degree of infrasound waveform similarity over more than a year, indicating a repeating source mechanism likely associated with explosive magma-water interaction. The seismic LP catalog shows that events began over a month prior to renewed explosive activity at the beginning of August 2021, and that lower frequency index (FI) LPs were generated in the week prior to the onset of explosions. We applied a recently developed machine learning tool (VOISS-Net, Tan et al., 2024 ) to catalog abundant broadband and harmonic seismic tremor recorded before and during the renewed explosive activity, along with LPs and explosions. The tremor catalogs complement the LP and explosion catalogs by filling out the seismic sequence with the dominant signal types. Together, these catalogs reveal a seismic sequence of renewed unrest that started with several weeks of LP events, followed by LPs with lower FI values and harmonic tremor in the days prior to explosive activity, and finally the onset of discrete explosions and broadband eruption tremor. We interpret this sequence as the ascent of a new pulse of magma that first interacted with the hydrothermal/groundwater system to produce LPs, followed by harmonic tremor, and that ultimately drove explosive magma-water interactions and periods of continuous ash emissions. The 2021 seismic sequence, in combination with long-term records of satellite SO 2 emissions, deformation from interferometric synthetic aperture radar (InSAR) analysis, ash sample analysis, infrasound, and volcano tectonic seismicity, allows us to interpret the entire 9-year period of unrest and eruption that began with an intrusion and earthquake swarm in 2014.

Alaska

Factors affecting the distribution of water-bearing fractures in the bedrock aquifers of West Virginia

Bedrock aquifers cover 23,601 square miles within the State of West Virginia and comprise 97.4 percent of the surficial area within the State; the remaining 2.6 percent (621 square miles) consists of alluvial sand-and-gravel and glacial outwash aquifers bordering the State’s major rivers. While West Virginia’s alluvial aquifers have been studied extensively, bedrock aquifers have only been characterized for studies completed in a few areas in Jefferson, McDowell, and Monroe Counties. Bedrock aquifers are water supplies for public supply, agriculture, industry, and residential homeowner use. In this study, the U.S. Geological Survey, in cooperation with the West Virginia Department of Environmental Protection Division of Water and Waste Management, provides a statewide assessment of the occurrence and distribution of fractures within bedrock aquifers of the State and the various topographic, physiographic, and lithologic influences controlling the occurrence and distribution of bedrock fractures. The results of this study provide an increased understanding of the distribution of fractures in bedrock aquifers in West Virginia and help to verify trends that have been suspected for many years but were never well documented or verified by data. The results confirmed that the density of fractures and those that were determined to be water bearing decrease significantly with depth. A statistically significant difference in the density of fractures was observed at a depth of 215 feet for wells in the Appalachian Plateaus Physiographic Province’s and in the Valley and Ridge Physiographic Province’s aquifers; a higher density of fractures and water-bearing fractures were above a depth of 215 feet than below that depth. This is an important consideration when drilling wells for residential, commercial, industrial, or agricultural water supply. Abandoned underground coal mines are commonly believed to form large pools of water in the interconnected mine entries in abandoned room and pillar coal mines. Such pools of water can and do exist in abandoned underground coal mines, but many mines lack open entries and are held up by overburden strata and pillars that can collapse and form aquifers comprised of vast interconnected rubble zones (gob), especially in older mines. Data assessed for this study showed that shale-corrected values of effective porosity for limestone aquifers in West Virginia had a median value of 2 percent and an average value of 4 percent and generally are mineralized with low effective porosity. Argillaceous or sandy limestone has a median shale-corrected porosity of 4 percent and an average shale-corrected porosity of 5 percent. The median and average shale-corrected porosity of sandstone aquifers was estimated to be 14 percent, but the median shale-corrected porosity for argillaceous or calcareous sandstone was 5 percent and the average shale-corrected porosity for argillaceous or calcareous sandstone was 6 percent. Even though shale has a relatively high total sonic porosity compared to other lithologies, shale and siltstone had relatively low shale-corrected porosity, ranging from 0 to 2 percent. Well yields were previously documented to be highest in valley settings, lowest on hilltops, and intermediate on hillsides. Transmissivity data provided by this study confirm this general pattern within the Appalachian Plateaus Province; however, the Valley and Ridge Province does not follow this pattern. While still lowest on hilltop settings, the highest well yields were in hillside settings. The trend for the Valley and Ridge Province was likely skewed because of 9 high-yield wells specifically targeting deeper thin limestone units, such as the Tonoloway and Helderberg Limestones, at depths with transmissivity in excess of 2,000 feet squared per day in Mineral County, West Virginia, or targeting karst aquifers in Berkeley, Jefferson, or Greenbrier Counties, West Virginia. Finally, water-bearing fractures have been hypothesized to comprise a small number of all fractures within a typical bedrock well in West Virginia. Data collected for this study support this theory. A total of 3,403 fractures were identified during this study; 3,151 (92.6 percent) of those fractures are low-transmissive, and only 252 (7.4 percent) fractures are water-bearing. Even though a well may contain many fractures, less than 8 percent are considered water-bearing fractures.

West Virginia

Earthquake scenario development in conjunction with the 2023 USGS National Seismic Hazard Model

We present earthquake scenarios developed to accompany the release of the 2023 update to the US Geological Survey National Seismic Hazard Model (NSHM). Scenarios can serve a range of local and regional needs, from developing proactive-targeted mitigation strategies for minimizing impending risk to aiding emergency management planning. These deterministic scenarios can also be used to communicate seismic hazard and risk to audiences who are not well versed in methods, such as probabilistic seismic hazard analyses. Specifically, we discuss the scenarios developed, challenges, and lessons learned in the development process, and how this work aided the development of the 2023 NSHM itself. In total, 28 scenarios were developed for Hawaii, Utah, Alaska, and Virginia considering the 2023 NSHM science, past scenario efforts, and input from local experts and stakeholders. Finally, we investigate how NSHM modeling decisions can change estimated impacts to Utah and Hawaii in more detail showing, for example, that a shallower dip of the Wasatch fault under Salt Lake City can increase predicted ground-motion intensities and therefore estimated losses and deaths.

Earthquake Spectra

Latitudinal gradients of snow contamination in the Rocky Mountains associated with anthropogenic sources

Seasonal snow is an important source of drinking water and recreation, and for agriculture in the Rocky Mountain region. Monitoring snow-water quality can inform on the effects to the albedo and energy balance of the snowpack, and the sources of natural and anthropogenic aerosol and gases. This study analyzed metals in the seasonal snowpack from water year (WY) 2018 for 49 sites. Calcium, lanthanum, and cerium concentrations support the importance of mineral dust to the southern Rocky Mountains. Mercury (Hg), zinc (Zn), and cadmium (Cd) concentrations showed a similar spatial pattern to mineral dust, whereas antimony (Sb) concentrations were highest in the northern Rocky Mountains. To assess the relative contributions from dust versus anthropogenic contaminant sources, enrichment factors (EF) were calculated, with values above 10 indicating anthropogenic contamination. For Cd, Hg, Sb, and Zn, EF values exceeded 10 at northern sites. These observations were compared to spatial trends of EF values of Hg from WY2009 to WY2018, regional monitoring networks, and back trajectory analyses. The agreement between these datasets revealed temporally consistent contaminant sources and/or transport processes to the northern Rocky Mountains snowpack. Sources include current and historical mining and smelting in the region. Strategies to limit the emissions of these metals to the Northern Rockies could benefit from focusing on remediation of contaminated sites, and continued monitoring and mitigation of active mining and smelting.

Colorado , Idaho, Montana, New Mexico, Utah, Wyomi

Reduced injection rates and shallower depths mitigated induced seismicity in Oklahoma

The proximity of wastewater disposal to the Precambrian basement is a critical factor influencing induced earthquake rates in the Central United States, but the impact of reducing injection depths has not been widely demonstrated. Beginning in 2015, state regulatory efforts in Oklahoma and Kansas mandated that wells injecting into the lower Arbuckle Group, a basal sedimentary unit, be backfilled with cement (i.e. “plugged back”) so that they inject into shallower formations. This plug back activity gives us a unique opportunity to investigate the relationship between injection depth and induced seismicity rate. To evaluate the impact that decreased injection rates and plug backs had on the seismicity rates, we create a suite of rate-state earthquake models. Observed seismicity rates are best fit when only lower Arbuckle volumes are considered, suggesting the lower Arbuckle injectors were primarily responsible for the seismicity and that plug backs were effective at isolating the injected volumes to shallower formations. Our models demonstrate that if these wells had not been plugged back, seismicity rates would be multiple times larger than they are today. We find that the combination of well plug backs and injection volume decreases can be effective strategies for reducing induced seismicity rates.

Oklahoma

‘The fish that stop’: Drivers of historical decline for Pacific cod and implications for modern management in an era of rapidly changing climate

n the Gulf of Alaska, a series of marine heat waves depleted Pacific cod ( Gadus macrocephalus ) biomass to the lowest abundance ever recorded and led to the fishery’s closure in 2020. Although the fishery has been productive for decades, this collapse may have historical precedents. Traditional knowledge holders refer to cod as ‘the fish that stop’, and there is a suggested period of decline in the 1930s. Here we conduct a catch reconstruction of the early commercial fishery (1864–1950), confirming a rapid catch decline in the 1920s and 1930s. Next, we evaluate evidence for possible drivers. We document changes to demand and technology that contributed to declining catch. However, we also find both qualitative and quantitative evidence of depletion, suggesting catch declines were not driven entirely by social factors. Overfishing may have contributed to localized catch declines as evidenced by declining catch rates in heavily fished localities. We also find evidence for climate as a driver of regional decline, with the period of catch decline characterized by up to 2°C higher temperatures as compared to the earlier period of high fisheries production. Our analysis underscores the importance of understanding long-term drivers of fisheries productivity and the value of linking fisheries and climate histories.

Alaska

Methodology for inclusion of produced and stored carbon dioxide in the U.S. Geological Survey Federal lands greenhouse gas inventory

The U.S. Geological Survey (USGS) has developed two new carbon dioxide (CO2) emissions and sequestration accounting methods for use in future reports. The first method is a Federal lease-produced CO2 emissions calculation for an update of the report, “Federal Lands Greenhouse Gas Emissions and Sequestration in the United States.” The methodology to incorporate Federal lease CO2 production emissions into the updated report relies on CO2 sales royalty data from the Office of Natural Resources Revenue (ONRR). The end usage points for the gas include enhanced oil recovery with CO2 (CO2-EOR), food and beverage, and chemical production. CO2-EOR is the main end point for natural CO2 production in the United States; it accounted for 94% of usage in 2022 [1]. Federal lands emissions from this sector are estimated at 460 metric tons of CO2 in 2022, a very small amount relative to most other Federal lands emissions sector estimates. The second new method, planned for a separate report, is a calculation of the geologic storage of CO2 on Federal lands. The second method estimates the CO2 stored under Federal surface lands and documents Federal climate change mitigation efforts. Currently, there is no storage of CO2 at an industrial level on Federal lands, however multiple proposals and projects are planned. This method was developed on non-Federal lands datasets in an effort to prepare for when these activities on Federal lands will require accounting. National estimates for CO2 geologic storage using this method, but without a Federal lands filtering step, totaled 8.0 million metric tons (Mt) in 2022. The two methods described here are new benchmark methods in a collection of accounting procedures to document the current state of greenhouse gas emissions and their storage on Federal lands. These benchmarks can then be used to measure any subsequent changes in emissions from or carbon storage beneath Federal lands. While the magnitude of the values is currently non-existent to small, emissions mitigation goals established by decision makers indicate that these values will grow, and their documentation will take on greater value and use.

continental United States

Simulation of the impacts of spring diversions on streamflow in the Strawberry Creek watershed, San Bernardino County, California, using an integrated hydrological model

The Strawberry Creek watershed, in the San Bernardino Mountains of southern California, features a group of natural springs known as Arrowhead Springs that have been augmented with diversions in the form of sub-horizontal borings and tunnels. Understanding the effects of these structures on streamflow through groundwater capture is crucial for managing surface-water resources in this watershed. In this study, we constructed the Strawberry Creek integrated hydrological model to increase this understanding. The Strawberry Creek integrated hydrological model is an integrated surface runoff and groundwater model, GSFLOW to simulate surface runoff and infiltration and groundwater flow. The model has 263 rows and 176 columns in each of its three layers. The model area includes the Strawberry Creek and four adjacent watersheds. The precipitation-runoff modeling system was calibrated using two streamflow gaging stations, and the groundwater and surface-water flow model was calibrated to the reported spring diversion discharge and a sparse number of groundwater-level measurements. The Strawberry Creek integrated hydrological model was run with and without diversions active, and simulated streamflow was compared, indicating that in the headwaters of Strawberry Creek, about 35 percent of the diversion flow was captured from base flow.

California

Evidence of mineral alteration in a salt marsh subterranean estuary: Implications for carbon and trace element cycling

Subterranean estuaries (STE) in salt marshes are biogeochemically active zones where interactions between terrestrial groundwater and seawater drive complex cycling of carbon and trace elements, influenced by mineral dissolution. These systems, characterized by fine-grained organic-rich peat overlying permeable coastal aquifers, play a crucial role as a blue carbon sink, yet their geochemical dynamics remain poorly understood. We investigated dissolved trace elements, carbon, silica, and radium isotopes in a salt marsh STE (Sage Lot Pond, Waquoit Bay, MA) over seasonal and annual cycles. Our results reveal that groundwater and estuarine water circulation through marsh peat and aquifer sediments leads to enrichments of dissolved organic and inorganic carbon (DOC and DIC), Si, Ba, and Mn, with variable source/sink behavior of Fe and net removal of U. Submarine groundwater discharge dominated Ba fluxes, whereas pore water drainage from marsh peat acted as the main sink for U and source of Si. Fe cycling was variable, with terrestrial Fe largely removed as groundwater passed through the STE, consistent with Fe-sulfide and amorphous phase formation. Radium isotope ratios identified two distinct subsurface flow pathways, influenced by metal-oxide cycling and organic matter breakdown. Si production was decoupled from DIC, suggesting Si originates from mineral alteration, whereas DIC results from both mineral weathering and microbial respiration. Silicate mineral alteration, coupled with marsh pore water drainage, accounts for up to 16% of annual DIC exports (66 g C m −2 y −1 ), highlighting the importance of STEs in coastal carbon and trace element cycling, especially as marshes face environmental change.

Massachusetts

Evaluation of best management practices at an edge-of-field site in the Eagle Creek watershed, Ohio, 2012–20

In 2010, the U.S. Geological Survey worked in partnership with the Great Lakes Restoration Initiative and the Natural Resources Conservation Service to identify farm fields in priority watersheds—watersheds critically important to Great Lake health. Intensive best management practices were implemented in five States to test the efficacy of best management practices, also referred to as agricultural conservation practices, to reduce sediment and nutrient runoff on agricultural fields. These fields were chosen as representatives of the priority watersheds because their farming practices and geographical conditions were common among farms in those watersheds. A farm site located near Findlay, Ohio, in the Eagle Creek watershed used two conservation practices, cover crop (Natural Resources Conservation Service practice standard 340) and variable rate technology (VRT; Natural Resources Conservation Service practice standard 590), which represent common practices in the Maumee priority watershed. This study monitored surface runoff and subsurface tile runoff at the Eagle Creek watershed site. The effects of cover crop and VRT nutrient application on sediment and nutrient runoff were assessed from October 2012 to September 2020. Cover crops were applied in fall of 2016, 2017, and 2018, coinciding with VRT applications. The following parameters were analyzed as part of the study: runoff discharge; peak discharge; and concentrations and loads of total phosphorus, particulate phosphorus, dissolved reactive phosphorus, and suspended sediment. High peak runoff events disproportionally affected nutrient concentrations and loads in runoff. A threshold of the 85th percentile of the peak discharge was chosen to assess the effect of cover crop and VRT nutrient application on nutrient mitigation at each gage. Runoff events below the 85th percentile were considered small, and runoff events over the 85th percentile were considered large. Concentrations from collected water-quality samples were used for the analysis; estimated concentrations were excluded from the statistical evaluation. All surface runoff parameters during small runoff events, except surface runoff volume, were significantly lower in the period with cover crops and VRT implementation. Contrarily, subsurface runoff parameters during small events were not statistically different between the periods with and without cover crops and VRT implementation. Surface and subsurface parameters during large events were also not statistically different between the periods with and without cover crops and VRT implementation. These results indicate cover crops and VRT may improve water quality during small runoff events, but additional best management practices that mitigate large runoff events may lead to greater water-quality improvement given the contribution of large events to overall losses.

Ohio

Geologic framework, hydrostratigraphy, and ichnology of the Blanco, Payton, and Rough Hollow 7.5-minute quadrangles, Blanco, Comal, Hays, and Kendall Counties, Texas

During 2023, the U.S. Geological Survey, in cooperation with the Edwards Aquifer Authority, revised an initial characterization completed during 2015–16 of the geologic framework, hydrostratigraphy, and ichnology of the Edwards and Trinity aquifers in the Blanco, Payton, and Rough Hollow 7.5-minute quadrangles in Blanco, Comal, Hays, and Kendall Counties, Texas. The purpose of this report is to present the updated geologic framework, hydrostratigraphy, and ichnology of the Trinity and Edwards Groups in those quadrangles. Rocks exposed in the study area are of the Lower Cretaceous Trinity Group and lower part of the Fort Terrett Formation of the Lower Cretaceous Edwards Group. The faulting and fracturing in the study area are part of the Balcones Fault Zone, an extensional system of faults active during the Paleocene to middle Eocene that generally trends southwest to northeast in south-central Texas. Hydrostratigraphically, the rocks exposed in the study area contain a section of the Edwards (Balcones Fault Zone) aquifer, upper zone of the Trinity aquifer, and middle zone of the Trinity aquifer. In the study area, the only hydrostratigraphic units of the Edwards (Balcones Fault Zone) aquifer remaining are (from top to bottom) VII and VIII, which cap several hills. The mapped hydrostratigraphic units of the upper zone of the Trinity aquifer are (from top to bottom) the cavernous, Camp Bullis, upper evaporite, fossiliferous, and lower evaporite. The mapped hydrostratigraphic units of the middle zone of the Trinity aquifer are (from top to bottom) the Bulverde, Little Blanco, Twin Sisters, Doeppenschmidt, Herff Falls (where present), Rust, Honey Creek, Hensell, and Cow Creek hydrostratigraphic units.

Texas

Simulation of the impacts of spring diversions on streamflow in the Strawberry Creek watershed, San Bernardino County, California, using an integrated hydrological model

The Strawberry Creek watershed, situated in the San Bernardino Mountains of southern California, features a group of natural springs known as Arrowhead Springs that have been augmented with diversions in the form of sub-horizontal borings and tunnels. Understanding the impact of these structures on streamflow through groundwater capture is crucial for managing surface-water resources in this watershed. In this study we constructed the Strawberry Creek integrated hydrological model (SCIHM) to increase this understanding. The SCIHM is an integrated surface runoff and groundwater model that uses the coupled groundwater and surface-water flow model (GSFLOW), which is based on the integration of the precipitation-runoff modeling system (PRMS) and the modular groundwater flow model commonly called MODFLOW, version MODFLOW-2005 software to simulate surface runoff and infiltration and groundwater flow. The model has three layers, 263 rows, and 176 columns. The model area includes the Strawberry Creek and four adjacent watersheds. The PRMS model was calibrated using two streamflow gaging stations and the GSFLOW model was calibrated to reported spring diversion discharge and a sparse number of groundwater-level measurements. The SCIHM was run with and without diversions active and simulated streamflow was compared, finding that in the headwaters of Strawberry Creek about 35 percent of the diversion flow was captured from base flow.

EarthArXiv

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