Geology ReportsSearch

SEARCH · Geology Reports

Results for “Climate of the Past”

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

A comparative analysis of OpenET for evaluating evapotranspiration in California almond orchards

The almond industry in California faces water management challenges that are being exacerbated by droughts, climate change, and groundwater sustainability legislation. The Tree-crop Remote sensing of Evapotranspiration eXperiment (T-REX) aims to explore opportunities to improve precision irrigation management for woody perennial cropping systems. Almond orchards in the California Central Valley were equipped with eddy covariance flux measurements to evaluate satellite remote sensing-based evapotranspiration (RSET) models. OpenET provides high-resolution (30-m spatial and daily temporal) RSET data, synthesizing decades of research for practical water management. This study provides an evaluation of OpenET performance at six almond sites covering a large range in soils, age, and variety. It also compares OpenET ensemble evapotranspiration (ET) data with applied irrigation and precipitation records over an additional 148 almond orchards located in the Central Valley of California. Results show OpenET models, including the ensemble ET value, produced reasonable and actionable ET values, with overall coefficient of determination (R 2 ) and mean absolute error values of 0.73- and 0.95-mm d −1 at the daily time step, respectively. However, given the temporal sampling of Landsat (8-day revisit) and the interpolation methods used, the assessed ET models had difficulty in capturing short-term variability in almond ET; for example, the rapid decline in measured ET observed as a response to lack of irrigation preceding and during almond harvest. The study also drew attention to the spatial complexity in scenarios where irrigated orchards are surrounded by hot/dry areas, causing discrepancies between measured and modeled ET values. In comparison with irrigation records, OpenET ensemble ET was capable of quantifying water input (applied irrigation + precipitation) in almond orchards to within 13 % when evaluating monthly data. Initial results presented here reinforce the idea that RSET models, such as in OpenET, are powerful tools, yet their application requires nuanced understanding and careful consideration of local conditions.

California

Simulation of groundwater flow in the Long Island, New York regional aquifer system for pumping and recharge conditions from 1900 to 2019

The U.S. Geological Survey has developed a transient, groundwater-flow model that simulates hydrologic conditions in the Long Island aquifer system as part of an ongoing (since 2016) multiyear, cooperative investigation with the New York State Department of Environmental Conservation. The goals of this investigation are to assist stakeholders and resource managers to evaluate the response of the hydrologic system to changes in future hydraulic stresses. Responses in the hydrologic system include changes in water levels in the hydrogeologic units; discharge to streams, coastal waters, and subsurface infrastructure; and the extent of saline groundwater in the aquifers. Hydraulic stresses include future water-supply management and changes in land use and infrastructure. The numerical model synthesizes a diverse set of physiographic, geologic, climatic, land-use, and historical population, water use, and infrastructure data to physically represent the Long Island aquifer system from land surface to bedrock and to simulate annual hydrologic conditions between 1900 and 2019. A three-dimensional hydrogeologic framework was developed from existing and recently collected borehole geologic and geophysical data collected as part of a companion drilling program. Water-transmitting properties of the principal aquifer sediments were defined in three dimensions from new and existing lithologic logs. The distribution of recharge from precipitation was estimated from landscape characteristics and climate data. Anthropogenic recharge from wastewater, leaky infrastructure, and storm runoff were estimated from population, infrastructure, and pumping data. Water-use data, including well locations, depths, and pumping rates, were obtained from historical sources and records and used to estimate pumping stresses continuously in time and space, at an annual average time scale. The data were incorporated into a three-dimensional numerical model using the U.S. Geological Survey finite difference modeling code MODFLOW 6; the model encompassed all of Long Island and surrounding surface waters and simulated historical hydrologic conditions from 1900 to 2019. The calibration process involved trial and error adjustments using prior knowledge to improve general fit to observations followed by an inverse calibration to update and optimize input parameters, using an iterative ensemble smoother algorithm implemented in PEST++ version 5.0. This resulted in a model that generally was in good agreement with observed, dynamically varying hydrologic conditions from 1900 to 2019. The calibrated model was used to develop two base-case models for scenario testing of future, hypothetical conditions where one represented average-annual conditions, and one represented average-seasonal conditions from 2010 to 2019. The model representing average-annual conditions was modified further to represent an alternate sea-level position of 6 feet above the North American Vertical Datum of 1988, and the model representing average-seasonal conditions was modified to represent the average seasonal effects of a 5-year drought imposed upon current hydrologic conditions. Recharge is the sole source of water to the aquifer system; groundwater discharges to coastal water and streams and is withdrawn by pumped wells. Model-estimated annual recharge ranged from about 11 inches in 1965 to 41 inches in 1983. On average, from 2010 to 2019, about 23 percent of water was pumped from wells, and about 47 and 27 percent discharged to coastal waters and streams, respectively; the remaining 4 percent was water that moved into storage in the aquifer matrix. Water levels on Long Island vary naturally during time in response to changes in recharge; the amount of variation is largest in the interior of the island, in areas with highest water table altitudes near groundwater divides and lowest near streams and the coastal waters. The total range of water table altitudes on Long Island between 1900 and 2019 ranged from near 0 to more than 70 feet in western parts of Long Island. The largest range in altitudes is in New York City and is associated with areas of large historical withdrawals between the 1920s and the late 1980s. Water table altitudes generally varied by less than 10 feet in eastern Suffolk County, where the aquifer is under more natural conditions. Saltwater intrusion is of great concern on Long Island, particularly in western Long Island where both the unconfined and confined parts of the aquifer system have been intruded in response to large-scale groundwater withdrawals; however, the volume of freshwater in the islandwide aquifer system only has changed by about 5 percent between 1900 and 2019. The decadal change in the freshwater volume was largest during the early and mid-20th century, corresponding to the largest historical pumping, but that volume change did not exceed 1 percent. The negligible change in freshwater volume suggests that saltwater intrusion as of 2019 was limited at an islandwide scale but continues to occur in local areas of Queens and Nassau Counties, adversely affecting current water supplies and limiting future water supplies for affected communities. The regional groundwater model developed for this investigation is a tool that can be used to help determine the viability of current and future water supplies at a regional scale and can be used to support development of additional models at finer scale to support more focused assessments of groundwater sustainability.

New York

Assessment of the interconnection between Tampa Bay and the Floridan aquifer system: Historical groundwater data compilation and analysis, 1976–2022

The U.S. Geological Survey used existing data collected after the last major navigational channel modification in the mid-1980s to investigate groundwater levels and chloride concentrations in wells in the Floridan aquifer system and other aquifers beneath and near Tampa Bay. Tampa Bay is located on the west-central coast of Florida and provides access for commercial shipping. In 2021, the U.S. Army Corps of Engineers began to investigate alternatives to improve the efficiency of the deep-draft navigation channels within Tampa Bay, reduce costs, improve safety, and analyze the implications of modifying navigational channels. The Floridan aquifer system underlies Tampa Bay and is the primary source of public water supply in the region. Modifications to the channels have the potential to expose the Floridan aquifer system to the saltwater in Tampa Bay, with the potential to increase salinity in wells in the region. Other factors affecting the groundwater levels and the location of the freshwater/saltwater interface were also examined, including changes in sea level, groundwater extraction, and variations in climate. Groundwater levels and well-construction reports were used to identify whether different aquifer units are well-connected. Twelve wells had available data before and after the last major channel modifications, which took place in the 1980s, with six datasets of chloride concentration available in areas along the northern and eastern coastline of Tampa Bay, which is nearest to historical dredging activities. Of these six, Kendall’s τ and p -values indicated increasing trends in chloride concentration for three datasets (TR 11-2, TR 10-2, 51), no trend in chloride concentration for two datasets (TR 9-3, 50), and a decreasing trend in chloride concentration for one dataset (TR 9-1). The upward trends in chloride concentration observed for TR 10-2 and 51 are likely the result of changes in local groundwater withdrawals. Well TR 11-2 had a gradual increasing trend in chloride concentration, fresh groundwater throughout the period of record, and a 3- to 4-foot increase in hydraulic head during the period of record, possibly caused by the construction and control of the Tampa Bypass Canal, resulting in changes to the regional potentiometric surface.

Florida

James Buttle review: The characteristics of baseflow resilience across diverse ecohydrological terrains

The dynamic storage of aquifers is the portion of groundwater that can potentially drain to any given point along a stream to create baseflow. Baseflow typically occurs year-round in perennial streams, though the characteristics and stability of dynamic storage are often most important to instream processes during extended dry periods (without precipitation and snowmelt) when runoff and quickflows are minimised. The term ‘baseflow resilience’ is defined for this review as the tendency of baseflow in streams to maintain a consistent volume and water quality year to year while under stress from climate variability and extremes, along with anthropogenic stressors such as water withdrawals, land use change, and water quality degradation. ‘Baseflow resilience’ has, in part, a user-defined meaning spanning water supply and water quality variables of primary interest. Watershed characteristics that directly impact resilience can often produce non-intuitive feedbacks that enhance some attributes of baseflow while simultaneously impairing others. For example, permeable stream corridor geology creates strong stream-groundwater hydrologic connectivity, yet fast groundwater drainage via preferential high-permeability flowpaths can lead to streamflow not being sustained during extended dry periods. Also, shallow groundwater sources are generally more immediately vulnerable to extreme events, warming, salinization, transpiration, and precipitation drought, compared to deeper groundwater. Yet baseflow drought in streams influenced by deeper groundwater can lag precipitation drought by years, and contaminant legacies may propagate through deep groundwater flowpaths to receiving waters for decades to centuries. Finally, irrigation withdrawals can intercept groundwater that would have drained to streams, and the application of irrigation may leach contaminants from the soil zone by unnaturally raising water tables, yet irrigation return flows can sustain baseflow and groundwater-dependent habitats in semiarid areas. This review covers the concept of hydrologic resilience in the context of stream baseflow processes and summarises the common hydrogeological controls on, and multiscale stressors of, dynamic groundwater storage. Further, we present several quantitative metrics to assess a range of water supply to water quality baseflow characteristics using both broadly available and boutique data types, a subset of which are demonstrated using data from the Delaware River Basin, USA.

Delaware River Basin

Using stable oxygen isotope dual-inlet isotope-ratio mass spectrometry to elucidate uranium transport and mixed 230Th/U calcite formation ages at the seminal Devils Hole, Nevada, natural laboratory

Rationale Vein calcite in Devils Hole has been precipitating continuously in oxygen-isotope equilibrium at a constant temperature for over 500 000 years, providing an unmatched δ 18 O paleoclimate time series. A substantial issue is that coeval calcite (based on matching δ 18 O values) has uranium-series ages differing by 12 000 years. Methods An unparalleled high-accuracy δ 18 O chronology series from continuously submerged calcite was used to correct the published uranium-series ages of non-continuously formed calcite in two cores, cyclically exposed by water-table decline during glacial–interglacial transitions. This method relies on the premise that the δ 18 O values of coevally precipitated calcite are identical, allowing matching calcite δ 18 O values to establish formation ages. Results Exposed calcite can have apparent ages that are 12 000 years too young due to unrecognized uranium mobility and resulting mixed ages identified in over 50 mixed uranium-series ages from previous studies. Secondary uranium in fluids, sourced from the formation or dissolution of porous carbonate deposits (folia) with high uranium-238 ( 238 U) concentrations, has migrated up to 10 mm into vein calcite. Conclusions The continuously submerged Devils Hole δ 18 O chronology is not explained by orbital forcing. Rather, this chronology represents a regional climate record in the southern Great Basin of sea-surface-temperature (SST) variations off California, variations that preceded the last and penultimate deglaciations by 5000 to approximately 10 000 years. Temporal discrepancies between the continuously submerged Devils Hole chronology and other regional δ 18 O records (e.g., the Leviathan chronology) can be explained by unrecognized cryptic, pernicious uranium mobility, leading to model estimations that may be thousands of years younger than actual ages. Consequently, paleo-moisture availability, water-table, and groundwater recharge models based on these mixed uranium-series ages are too young by as much as 12 000 years. The potential for post-formation uranium addition in subaerial cores and speleothems underscores the need for caution in uranium-series dating, highlighting δ 18 O time-series comparisons as a method for identifying mixed ages.

Nevada

Comparative assessment of STIC sensors, streamflow and rain gauges for quantifying river connectivity in intermittent systems

In intermittent stream systems, including those occurring in Texas, USA, the severity of low-flow conditions, duration of seasonal disconnection, and frequency of no-flow events have been amplified by drought. Documentation of these no-flow events is necessary to evaluate ecosystem health. However, many intermittent reaches remain un-gauged given that perennial river sec-tions are often prioritized for gauge placement. Our objectives were to 1) document stream flow using Stream Temperature, Intermittency, and Conductivity (STIC) loggers to determine the frequency and duration of no-flow events in intermittent tributaries of the Colorado River, Texas and 2) compare logger data to publicly available data from streamflow discharge and precipitation gauge networks to understand differences among these data types for drying event characterization. We use these comparisons to summarize benefits and limitations of the application of in-stream data loggers. STIC loggers were deployed at 19 sites, one in each pool and riffle habitat of a stream reach. STIC loggers recorded a measurement of relative conductance every six hours from June 2022 to March 2024, which was used to determine the presence or absence of flow connectivity in a reach. No-flow duration among intermittent reaches varied between 37 and 270 days across tributaries during an ongoing drought in the study area. Overall, logger data was more precise than discharge data for characterizing no-flow events or precipitation data when documenting presence of water in the stream channel due to runoff. Lack of discharge gauges in intermittent tributaries left large sections of stream reaches undocumented and resulted in mischaracterization of flow patterns. Drought severity across the tributaries did not follow longitudinal patterns that would be expected by the climatic precipitation gradient of the study area. More research is needed to determine if factors such as population size affect severity. Likewise, precipitation data did not correlate well with logger water presence data, lacking consideration for groundwater recharge, soil hydrophobicity, and surface compaction. This study shows that to monitor no-flow events, detailed spatial datasets are necessary and that STIC loggers are useful tools that provide data to fill spatial information gaps and facilitate more accurate flow characterization and water presence data in intermittent systems.

Texas

Connecting flood-related fluvial erosion and deposition with vulnerable downstream road-stream crossings

Fluvial erosion is increasingly responsible for infrastructure and building damages associated with floods as the intensity of extreme rainfalls hit rural and urban rivers in a variety of climate settings across the United States. Extreme floods in 2016 and 2018 caused widespread culvert blockages and road failures, including extensive damage along steep tributaries and ravines in the Marengo River, Wisconsin, watershed during 2016 and 2018. A study conducted by the U.S. Geological Survey (USGS), Wisconsin Wetlands Association (WWA), Ashland County, and the Northwest Wisconsin Regional Planning Commission (NWRPC) investigated the special concern of fluvial erosion hazards (FEHs) associated with gullying, streamside landslides, and the loss of wetland storage in headwaters. In 2019, a pilot study was begun to map and classify ephemeral and perennial streams and wetlands in terms of their sensitivity to FEHs. This study combined data from field-based rapid geomorphic assessments (RGAs) coupled with a stream network-wide geographic information system (GIS) approach for mapping stream segments, referred to as fluvial process zones (FPZ), sensitive to erosion, deposition, and channel change. The GIS approach used nationally available 10-meter (m) resolution topology and an extended stream network to map FPZs based on Strahler stream order, stream power, channel slope, presence of adjacent steep valley sides and headwater flats, and adjacent landform setting. Bankfull channel widths derived from RGA-based hydraulic geometry curves combined with drainage areas, an estimate of bankfull flow, and channel slope were used to calculate specific stream power for the FPZs. Lastly, the FPZs were characterized by their location within three major landform settings that affect erosion potential. The resulting vulnerability maps provided a screening framework to identify FPZs that are sensitive to incision, gullying and mass wasting along steep headwater ephemeral channels, as well as downstream perennial channels that have the potential for valley-side landslides, coarse sediment deposition, and channel change. Lastly, each FPZ was characterized in terms of hydrologic alteration associated with ditching. The vulnerability mapping products and rankings of sensitivity of FPZs will ultimately be used by Ashland County and their collaborators to prioritize natural flood management projects that mitigate FEHs, restore hydrology, and reconnect channels with adjacent wetlands and floodplains.

Wisconsin

Deep-ocean macrofaunal assemblages on ferromanganese and phosphorite-rich substrates in the Southern California Borderland

Mineral-rich hardgrounds, such as ferromanganese (FeMn) crusts and phosphorites, occur on seamounts and continental margins, gaining attention for their resource potential due to their enrichment in valuable metals in some regions. This study focuses on the Southern California Borderland (SCB), an area characterized by uneven and heterogeneous topography featuring FeMn crusts, phosphorites, basalt, and sedimentary rocks that occur at varying depths and are exposed to a range of oxygen concentrations. Due to its heterogeneity, this region serves as an optimal setting for investigating the relationship between mineral-rich hardgrounds and benthic fauna. This study characterizes the density, diversity, and community composition of macrofauna (>300 μm) on hardgrounds as a function of substrate type and environment (depth and oxygen ranges). Rocks and their macrofauna were sampled quantitatively using remotely operated vehicles (ROVs) during expeditions in 2020 and 2021 at depths above, within, and below the oxygen minimum zone (OMZ). A total of 3,555 macrofauna individuals were counted and 416 different morphospecies (excluding encrusting bryozoans and hydrozoans) were identified from 82 rocks at depths between 231 and 2,688 m. Average density for SCB macrofauna was 11.08 ± 0.87 ind. 200 cm −2 and mean Shannon-Wiener diversity per rock (H′ [loge] ) was 2.22 ± 0.07. A relationship was found between substrate type and macrofaunal communities. Phosphorite rocks had the highest H′ of the four substrates compared on a per-rock basis. However, when samples were pooled by substrate, FeMn crusts had the highest H′ and rarefaction diversity. Of all the environmental variables examined, water depth explained the largest variance in macrofaunal community composition. Macrofaunal density and diversity values were similar at sites within and outside the OMZ. This study is the first to analyze the macrofaunal communities of mineral-rich hardgrounds in the SCB, which support deep-ocean biodiversity by acting as specialized substrates for macrofaunal communities. Understanding the intricate relationships between macrofaunal assemblages and mineral-rich substrates may inform effects from environmental disruptions associated with deep-seabed mining or climate change. The findings contribute baseline information useful for effective conservation and management of the SCB and will support scientists in monitoring changes in these communities due to environmental disturbance or human impact in the future.

California

StreamStats—A quarter century of delivering web-based geospatial and hydrologic information to the public, and lessons learned

StreamStats is a U.S. Geological Survey (USGS) web application that provides streamflow statistics, such as the 1-percent annual exceedance probability peak flow, the mean flow, and the 7-day, 10-year low flow, to the public through a map-based user interface. These statistics are used in many ways, such as in the design of roads, bridges, and other structures; in delineation of floodplains for land-use zoning and setting of insurance rates; for regulatory purposes, such as the permitting of wastewater discharges; and for hydrologic and climate change studies. StreamStats was first developed for Massachusetts and released in 2001. The application provided users with the ability to obtain streamflow statistics computed from data collected at USGS streamgages and to obtain estimates of streamflow statistics for user-selected ungaged sites. Massachusetts StreamStats used geographic information system software and digital mapping to compute drainage-basin characteristics, which were then used in statistical models to estimate streamflow statistics for the user-selected sites. The statistical models were in the form of equations that were developed through a process known as regression analysis. StreamStats was the first known web application with the ability to do interactive geoprocessing. The utility of Massachusetts StreamStats was instantly apparent, leading the USGS to develop a version of StreamStats that could be implemented nationally. USGS State offices normally were required to develop custom regression equations and prepare local digital mapping data needed for implementing StreamStats for their States. Funding needed to complete this work usually was provided through cooperative agreements between the USGS and State agencies. In 2004, Idaho became the first to be released in the national version of StreamStats. By 2023, 44 States were fully implemented and six were undergoing implementation. StreamStats has undergone many modifications over the years to keep up with changes to the underlying software and to add functionality. Customized functionality and separate linked StreamStats applications were developed for several States. Meeting the high demand for additions and improvements to StreamStats while also adhering to budgetary constraints has, at times, been challenging. The StreamStats development team has identified numerous additional improvements that could be made to provide better performance and more functionality. The lessons learned from the experience of building and operating StreamStats for nearly 25 years could be relevant to others interested in pursuing efforts of a similar scale.

Circular

A time-space model of graphite mineral systems

Increasing demand for graphite in energy storage systems warrants review of graphite ore genesis in a mineral systems framework. Orogenic graphite encompasses the metamorphic and orogenic mineral systems that produce flake graphite and hydrothermal vein (lump and chip) graphite deposits, respectively. A common feature of orogenic graphite deposits is an association with upper amphibolite- to granulite-facies metasedimentary rocks in continent-continent or continent-island arc collisional orogens. Orogenic flake graphite deposits form primarily through graphitization of organic carbon during regional metamorphism, but strain localization and partial melting of pelitic protoliths are likely important processes for graphite grade and quality enrichment. Orogenic vein graphite deposits precipitate from hydrothermal fluids, possibly derived from metamorphism or anatexis at depth. Decarbonation reactions in mixed calcareous-carbonaceous metasedimentary sequences are the most likely carbon sources for the veins. In contrast, intrusion-related graphite includes magmatic-hydrothermal and metamorphic mineral systems that form primarily in continental arc settings via the interaction of magmas with carbonaceous sedimentary packages. Magmatic-hydrothermal flake graphite deposits are hosted in plutonic and volcanic rocks, and result from the exsolution of CO 2 -CH 4 -rich fluids from melts contaminated by such packages. Contact metamorphism of carbonaceous sedimentary rocks by plutons produces some microcrystalline (amorphous) graphite deposits, including many in China. Compilation of geologic data from known graphite deposits globally suggests that pulses of carbon deposition in the Paleoproterozoic, Mesoproterozoic, and Neoproterozoic provided source carbon material. Subsequent supercontinent orogenesis at ca. 2,100 to 1,700 Ma (Columbia), ca. 1,300 to 1,000 Ma (Rodinia), and ca. 650 to 500 Ma (Gondwana) resulted in the genesis of orogenic flake and vein graphite deposits, where favorable geologic components overlapped with organic ± carbonate carbon-rich strata. Cryogenian deposition of graphite protoliths and Cryogenian – Cambrian metamorphic mineralization account for nearly 75% of all known resources globally and coincide with profound carbon isotope excursions and climate variability, implying a link with the global carbon budget. Comparatively few graphite deposits are associated with Pangea-forming orogens, attributed to less exhumation and/or denudation. High-temperature metasedimentary belts containing organic carbon-rich protoliths are most favorable for hosting orogenic flake graphite deposits, whereas sequences that also contain carbonate protoliths are favorable for orogenic graphite veins. Continent-scale orogenic belts may host both deposit types along with vanadium deposits. Use of a time-space mineral systems framework for graphite deposits can improve exploration models needed to ensure future supply of this critical mineral and provide insights into Earth’s long-term carbon cycle.

Mineralium Deposita

Assessment and characterization of ephemeral stream channel stability and mechanisms affecting erosion in Grand Valley, western Colorado, 2018–21

The Grand Valley in western Colorado is in the semiarid Southwest United States. The north side of the Grand Valley has many ungaged ephemeral streams, which are of particular interest because (1) the underlying bedrock geology, Late Cretaceous Mancos Shale, is a sedimentary rock deposit identified as a major salinity contributor to the Colorado River and (2) despite infrequent streamflows of short duration, monsoon-derived floods in these ephemeral streams can carry substantial amounts of sediment downstream, affecting upstream and downstream banks and channel cross sections. The study area is of interest, because salinity, or the total dissolved solids concentration, in the Colorado River causes an estimated $300 million to $400 million per year in economic damages in the United States, and it is estimated 62 percent of the Upper Colorado River Basin’s total dissolved solid loads originate from geologic sources. In an effort to minimize salt contributions to the Colorado River from public lands administered by the Bureau of Land Management, a comprehensive salinity control approach is typically used to reduce nonpoint sources of salinity through land management techniques and practices. In 2018, the U.S. Geological Survey, in cooperation with the Bureau of Land Management, began an assessment of ephemeral streams located on the north side of the Grand Valley, western Colorado, to characterize stream channel stability and identify mechanisms affecting erosion. The U.S. Geological Survey developed a method for automatically extracting channel cross-section geometry from existing remotely sensed terrain models. Based on estimated flood stage and surrogate streamflows, hydraulic characteristics were calculated. Furthermore, the channel geometries and hydraulic characteristics were used to estimate channel stability using a statistical model. Cross-section stabilities were determined from a stream channel stability assessment for a subset of 1,406 visited (field observed) locations out of 13,415 cross sections, which were delineated from remotely sensed terrain models. The application of Manning’s resistance equation in combination with multiple logistic regression models demonstrated channel stability can be estimated with a 0.845 goodness of fit for a validation dataset when using a combination of drainage area, width-to-depth ratio, sinuosity, and shear stress as the explanatory variables. Stream channel stability was extrapolated for 13,415 unvisited (not field observed) cross sections using the multiple logistic regression model and defined explanatory variables. Mapping of the ephemeral streams and their associated stabilities may be used by the Bureau of Land Management to prioritize areas for remediation or changes in management strategies to reduce sediment and salinity loading to the Colorado River. The study found channel stability within the ephemeral streams to be spatially variable, longitudinally discontinuous, and dictated by changes in channel bed slope. The stable ephemeral streams were relatively wide and shallow and often had smaller drainage areas with less potential for producing shear stresses capable of overcoming channel adhesion. A change in channel bed slope can provide the means necessary to generate shear stresses appropriate to initiate erosion and a subsequent stability transition to incising channels. Channel widening happens when either or both banks of an incising channel reach a critical height for mass wasting, or when channel curvature causes higher sidewall stress. Regardless, widening channels can promote increases in sinuosity and subsequently reduce steep channel bed slopes. Consequently, stable and widening channels can have comparable bed slopes, making channel bed slope a poor explanatory variable to predict channel stability overall, despite its function to initiate channel instability. The results were based on a surrogate 0.10 annual exceedance probability (AEP; return period equal to the 10-year flood) interval streamflow, although it was recognized fluctuations in streamflow would also affect channel stability. Past and current changes within the study area affect streamflow; therefore, mechanisms affecting erosion include land use disturbances, soil compaction, loss of vegetation cover, drought, less frequent and more extreme precipitation, and fires—which all intensify the potential runoff and erosion within the study area.

Colorado

Editorial: From cold seeps to hydrothermal vents: Geology, chemistry, microbiology, and ecology in marine and coastal environments

This Research Topic compiles contemporary studies on cold seeps, hydrothermal vents, mud volcanoes, and related seafloor features that are associated with focused fluid emissions and the transfer of carbon, other chemical species, and sometimes heat from the geosphere to the ocean. Because these features sometimes tap fluids and gas originating kilometers below the seafloor, they provide an important window into deep processes that are otherwise inaccessible to scientists. At the shallow portion of their journey, migrating fluids nearing the seafloor contribute to a range of unique biological, physical, and chemical processes within the sediments themselves and at the sediment-water interface. Seafloor fluid emissions play a critical role in global biogeochemical cycles, ocean chemistry, and possibly even climate change. Seafloor leakage points often emit hydrocarbon gases (especially methane and CO 2 ) and are sometimes the loci for deposition of seafloor minerals that have economic value. A burgeoning area of research focuses on natural products generated at these features, seeking compounds with potential pharmaceutical or other applications. Multidisciplinary studies have become routine for characterization of seafloor fluid emission sites, attesting to the inseparability of geologic, physical, chemical, and biological processes in these settings. It is increasingly common for researchers to combine in a single research cruise: subbottom imaging and seafloor mapping; porewater and water column geochemistry and gas sampling; sediment retrieval for lithologic, biostratigraphic, and solid phase analyses; and studies of benthic and subseafloor communities at the microbial to macrofaunal scales. This multidisciplinary approach has the advantage of ensuring the spatial and temporal coincidence of surveys and samples, an important factor at highly dynamic seafloor fluid emission sites. In addition, researchers often use remotely operated vehicles (ROVs), autonomous underwater vehicles (AUVs), or human-occupied vehicles (HOVs) to record video of the seafloor, compile photomosaics, collect targeted samples, and survey with high-resolution geophysical near-seafloor systems, providing a degree of detail about seafloor fluid emission sites that is unprecedented compared to most areas of the deep ocean. While rarer, long-term cabled observatories or shorter-term deployments of portable observatories are also used at some loci for seafloor fluid flux and are particularly helpful for capturing temporal variations at these dynamic features. Here we summarize the Research Topic’s contribution to multidisciplinary seafloor emission studies in the categories of cold seeps, mud volcanoes, and hydrothermal vents. Figure 1 shows the geographic distribution of the studies in this Research Topic and key features referred to in this Introduction.

Frontiers in Earth Science

Changes in soil erosion caused by wildfire: A conceptual biogeographic model

Soil erosion rates after wildfire are strongly controlled by intrinsic properties such as topography, weather, climate, soil, and vegetation. These landscape and hydroclimatic properties are important in determining post-fire erosion rates; however, their influence on post-fire erosion and their interaction with the intensity of a wildfire remains uncertain. A key limitation in resolving this uncertainty is the lack of conceptual models and frameworks for organising data related to the geomorphic sensitivity of landscapes to wildfire. Our aim is to develop a framework for consolidating understanding of post-fire erosion in the context of hydroclimatic conditions which contribute to system states, for example soil and vegetation properties, and wildfire regime. The framework is developed around a simple conceptual model where the change in erosion due to wildfire is a product of change in runoff generation and sediment supply, which is strongly related to landscape net primary productivity (NPP). We hypothesised that geomorphic sensitivity to wildfire should vary as a unimodal humped relationship across a gradient of NPP, peaking at an intermediate level. To develop this framework and to test the hypothesis, we first review intrinsic soil and vegetation properties related to the supply and transport of sediment from burned and unburned hillslopes. Net primary productivity is systematically related to these intrinsic properties because it integrates many processes involved in soil and vegetation development. Empirical data indicate a trend in the change in surface runoff generation with NPP after wildfire, peaking at an NPP of approximately 15 Mg C ha −1 y −1 . A simple model of fuel availability and soil heating are correlated with a similar “humped” trend in sediment supply. These results are consistent with our conceptual model, which indicates that sediment supply and runoff contribute towards a distinct peak in wildfire effects on erosion at an intermediate level of NPP. We propose that landscapes of intermediate NPP typically have the highest quantity of fuel available to burn, which cause large changes to the soil surface properties. Landscapes at intermediate NPP also tend to produce intrinsic soil and vegetation properties that promote erosion after wildfire. The interplay between these short and long-term landscape characteristics is strongest at intermediate levels of NPP. Our proposed biogeographic model of geomorphic sensitivity to wildfire was supported by erosion data from burned hillslope and zero-order catchments studies from a range fire-prone landscapes in Australia and North America. Our proposed conceptual model will help identify areas most vulnerable to post-fire erosion changes.

Geomorphology

Computation of regional groundwater budgets for the Virginia Coastal Plain aquifer system

Computation of detailed groundwater flow budgets for subdivisions of Virginia’s Coastal Plain aquifer system has enabled quantification and more thorough understanding of groundwater flow within this important water resource. A zone budget analysis conducted on previously published groundwater models of the Virginia Coastal Plain and Virginia Eastern Shore shows that groundwater conditions vary substantially throughout the Coastal Plain aquifer system due to local variations in hydrogeology and historical and ongoing variations in groundwater use and management. Decades of substantial groundwater withdrawal from the Coastal Plain aquifer system have fundamentally altered groundwater flow from pre-development conditions. Rates of sustainable withdrawal are limited because the downward groundwater flow rate into confined aquifers supplying groundwater is a relatively small portion of the total groundwater water budget for the aquifer system. Analyses of groundwater budgets from the Virginia Coastal Plain model show that groundwater flow is generally outward from the surficial aquifer to rivers and coastal water bodies and downward through a series of underlying aquifers and confining units to the Potomac aquifer, which is the deepest aquifer and the source of most groundwater withdrawals. Downward flow into the Potomac aquifer currently is estimated to be only 7 percent of total net precipitation-derived net recharge at the land surface but makes up about 66 percent of inflow to the aquifer in Virginia, with much of the remaining inflow occurring laterally from areas outside of defined groundwater budget regions in Virginia. For several decades prior to 2010, high rates of withdrawal from the Potomac aquifer resulted in substantial decline in groundwater storage in the aquifer and in most overlying aquifers and confining units. From 2010 to 2025, rates of withdrawal substantially lower than the historical maximum have resulted in small net increases in groundwater storage in the confined aquifer system for most regions of the Virginia Coastal Plain. Nevertheless, for the same period, groundwater storage for the entire model domain continues to incrementally decline, indicating that storage recovery in Virginia is offset by a continued decrease in storage in areas beneath the Chesapeake Bay or in adjacent areas of Maryland and North Carolina. Withdrawals from the Potomac aquifer have induced substantial downward flow which is a large part of groundwater budgets for confined aquifers such as the Potomac. Downward groundwater flow continues under current conditions, but because vertical flow rates are a function of the difference between water pressure in the upper surficial systems and lower confined units, those rates are lower than those in earlier decades as the confined water levels partially recover from larger groundwater withdrawals in the past. Geographically, groundwater flow is generally inward from perimeter regions of the Virginia Coastal Plain toward central regions with the largest withdrawal rates. Estimated groundwater inflow from coastal regions could be contributing to saltwater intrusion, though that was not measured directly in this study. Analyses of groundwater budgets from the Virginia Eastern Shore peninsula, a geographic region of the Virginia Coastal Plain, show that groundwater flow for that isolated aquifer system is generally outward from the surficial aquifer to coastal water bodies and downward into the confined Yorktown-Eastover aquifer system, which is the source of most withdrawals. Downward groundwater flow into the confined Yorktown-Eastover aquifer system is estimated to be less than 2 percent of total recharge and less than 9 percent of net recharge at the water table but makes up over 93 percent of all inflow to the confined aquifer system. Decades of substantial but relatively consistent groundwater withdrawals have induced greater downward flow rates into the confined aquifer system but also have resulted in loss of groundwater from storage. Currently, estimated storage loss accounts for slightly under 7 percent of withdrawals from the confined aquifer system. The current withdrawal rate from the confined Yorktown-Eastover system is near the highest reported rate for the Eastern Shore, which means that the storage depletion is expected to continue, even though groundwater levels appear to be relatively stable. Estimated groundwater flow rates upward from the confining unit underlying the Yorktown-Eastover system and small rates of inflow from coastal water bodies underscore ongoing concerns about up-coning and lateral intrusion of salty groundwater.

EarthArXiv

Assessment of coastal and fluvial morphodynamic changes using Structure-for-Motion: A case study of the Sfȃntu Gheorghe Mouth (Danube Delta, Romania)

The ability to accurately map erosion, flooding, and habitat loss in coastal environments is crucial for formulating national strategies aimed at preventing and mitigating the impacts of natural disasters. A fundamental component of this process is the implementation of coastal morphodynamics monitoring through Structure-from-Motion (SfM) techniques, utilizing high-resolution 2D/3D data obtained from aerial photogrammetry. To assess morphodynamic changes over a three-year period (2022 – 2024), several SfM-based photogrammetric studies were conducted, each year, in the Romanian sector of the Danube-Black Sea coastal zone, specifically at the mouth of one of the Danube River distributaries (Sf Gheorghe branch) into the Black Sea, and along the left bank, near Sf Gheorghe locality, located within the Danube Delta Biosphere Reserve (DDBR). The essential equipment for aerial photogrammetry comprises Unmanned Aerial Vehicles (UAVs) and Global Navigation Satellite Systems (GNSS). In this study, the UAV used was a DJI Mavic 3T (Enterprise/Thermal) drone, complemented by two Trimble R12i and R4 GNSS systems, as well as approximately 10 Ground Control Points (GCPs). Data acquisition and processing were carried out using specialized photogrammetric software (Agisoft Metashape) along with various GIS tools (e.g., Blue Marble Geographics Global Mapper and ESRI ArcMap). The photogrammetric products generated for the study, as detailed in this paper, include Digital Elevation Models (DEMs), Digital Terrain Models (DTMs), orthomosaics (orthophotos), and others. At Sfântu Gheorghe beach, a comparison between 2023 and 2024 photogrammetric surveys revealed that the left bank of the Sf. Gheorghe Arm, at the river mouth into the Black Sea, suffered from a twist (erosion) of up to 64 metres. Additionally, on the selected perimetre (total area of 31,910 square meters ) from the beach and dune zone of Sf. Gheorghe, an area of up to 16,202 square meters was eroded between 2023 and 2024. This contrasts with the period between 2022 and 2023, during which deposition predominated. Erosion at the Danube mouths and the adjacent Black Sea coastline is driven by a complex interaction of natural and anthropogenic factors. Natural processes, including subsidence, sea-level rise, and episodic extreme storm events, contribute significantly to coastal dynamics. Meanwhile, human-induced factors, such as upstream hydrotechnical works that limits sediment transport, cutting of navigation canals, as well as the exacerbating effects of climate change, further accelerate erosion. The recent Structure-from-Motion (SfM) surveys provide essential quantitative data, enabling a detailed analysis of both short-term and long-term morphodynamic changes influenced by seasonal variations and extreme hydrometeorological events in this highly dynamic coastal system.

Sfȃntu Gheorghe Mouth (Danube Delta)

Computation of regional groundwater budgets for the Virginia Coastal Plain aquifer system

Computation of detailed groundwater flow budgets for subdivisions of the Virginia Coastal Plain aquifer system has enabled quantification and more thorough understanding of groundwater flow within this important water resource. A zone budget analysis based on previously published groundwater models of the Virginia Coastal Plain and Virginia Eastern Shore indicates that groundwater conditions vary substantially throughout the Coastal Plain aquifer system because of local variations in hydrogeology and historical and ongoing variations in groundwater use and management. Decades of substantial groundwater withdrawal from the Coastal Plain aquifer system have altered groundwater flow from predevelopment conditions. Rates of sustainable withdrawal are limited because the downward groundwater flow rate into confined aquifers is a relatively small part of the total groundwater budget for the aquifer system compared to the rate of recharge at the land surface. Analyses of groundwater budgets from the Virginia Coastal Plain model indicate that groundwater flow is generally outward from the surficial aquifer to rivers and coastal waterbodies and downward through a series of underlying aquifers and confining units to the Potomac aquifer, which is the deepest aquifer and the source of most groundwater withdrawals. Downward flow into the Potomac aquifer is estimated to be only 7 percent of total net precipitation-derived net recharge at the land surface but makes up about 66 percent of inflow to the aquifer in Virginia, with much of the remaining inflow occurring laterally from outside of defined groundwater budget regions in Virginia. For several decades prior to 2010, high rates of withdrawal from the Potomac aquifer resulted in substantial decline in groundwater storage in the aquifer and in most overlying aquifers and confining units. From 2010 to 2023, rates of withdrawal substantially lower than the historical maximum resulted in small net increases in groundwater storage in the confined aquifer system for most regions of the Virginia Coastal Plain. Nevertheless, for the same period, groundwater storage for the entire model domain continues to incrementally decline, indicating that storage recovery in Virginia is offset by a continued decrease in storage in areas beneath the Chesapeake Bay or adjacent areas of Maryland and North Carolina. Withdrawals from the Potomac aquifer have induced substantial downward flow which is a large part of groundwater budgets for confined aquifers such as the Potomac. For the most recent simulated conditions (2023) downward groundwater flow continues, but because vertical flow rates are a function of the difference between water pressure in the upper surficial systems and lower confined units, rates of downward flow are lower than those in earlier decades as the confined water levels partially recover from larger groundwater withdrawals in the past. Geographically, groundwater flow is generally inward from perimeter regions of the Virginia Coastal Plain toward central regions with the largest withdrawal rates. Groundwater inflow from coastal regions could be contributing to saltwater intrusion, even though that was not measured in this study. Analyses of groundwater budgets from the Virginia Eastern Shore peninsula, a geographic region of the Virginia Coastal Plain, indicate that groundwater flow for that isolated aquifer system is generally outward from the surficial aquifer to coastal water bodies and downward into the confined Yorktown-Eastover aquifer system, which is the source of most withdrawals. Downward groundwater flow into the confined Yorktown-Eastover aquifer system is estimated to be less than 2 percent of total recharge and less than 9 percent of net recharge at the water table but makes up more than 93 percent of all inflow to the confined aquifer system. Decades of substantial but relatively consistent groundwater withdrawals have induced greater downward flow rates into the confined aquifer system but also have resulted in loss of groundwater from storage. For the most recent simulated period (2023), estimated storage loss accounts for slightly under 7 percent of withdrawals from the confined aquifer system. The reported withdrawal rate for this period from the confined Yorktown-Eastover system is near the highest reported rate for the Virginia Eastern Shore, which means that the storage depletion is expected to continue, even though groundwater levels appear to be relatively stable. Estimated groundwater flow rates upward from the confining unit underlying the Yorktown-Eastover system and low rates of inflow from coastal water bodies underscore ongoing concerns about up-coning and lateral intrusion of salty groundwater.

Virginia

Chemostratigraphy of the Cretaceous Hue Shale in Arctic Alaska: Exploring paleoceanographic controls on trace element enrichment, organic matter accumulation, and source-rock evolution

We document chemostratigraphy in an outcrop of late Albian to early Campanian (∼103–82 Ma) marine source rocks to better understand paleoenvironmental controls on trace element (TE) enrichment and organic matter accumulation in the distal Colville foreland basin of Arctic Alaska and how those drivers are linked to arc volcanism and successions of Cretaceous oceanographic and climatic biogeochemical events. This unique, 113-m-thick section of Cretaceous Hue Shale deposited during a series of previously undocumented Arctic Cretaceous oceanic anoxic events (Lease et al., 2024) is the only known exposure of thermally immature (0.48–0.52% R o , random vitrinite reflectance) Hue Shale in Arctic Alaska. Strata comprise mainly clay-rich mudstone with elevated total organic carbon (TOC) and hydrogen index values reaching 26.3 wt% (mean = 7.5 wt%) and 689 mg hydrocarbon (HC)/g TOC (mean = 385 mg HC/g TOC), respectively. Maceral composition consists predominantly of fluorescent amorphous organic matter, with abundant brightly fluorescent alginite, including Tasmanites , acritarchs, and Leiosphaeridia . Discrete layers of volcanic ash (preserved as bentonite) are present throughout the section and provide quantitative age control based on U–Pb dates. Chemostratigraphic trends are investigated to advance our understanding of local oceanographic conditions and controls on stratigraphic and temporal heterogeneity of Brookian source rocks. Concurrent sedimentary enrichment in Mo, U, V, Pb, and Cu across the Albian–Cenomanian boundary of the exposed basal gamma-ray zone, may reflect anoxic to euxinic benthic redox conditions favoring organic matter accumulation and preservation. Fluctuating degrees of anoxia-euxinia are inferred throughout the overlying Hue Shale succession, reflected by varying patterns of TE enrichment and TE–TOC covariation. Elevated C org /P molar ratios (>250) across most of the section, with several values exceeding 690, signify that enhanced biological productivity is sustained throughout deposition. Enhanced productivity, recorded by both C org /P and excess Ba, also parallels increases in source rock richness (elevated TOC and S 2 values) during the late Albian–early Cenomanian and late Cenomanian – Turonian. Enhanced productivity and variations in oceanic circulation/stratification likely both drove changes in benthic redox conditions that favored organic carbon accumulation and preservation. Increased continental arc volcanism (e.g., Okhotsk–Chukotka volcanic belt) and High Arctic Large Igneous Province magmatic eruptions throughout the Cretaceous, inferred to have influenced nutrient cycling and local aqueous nutrient availability, also have been invoked as potential drivers of organic carbon burial and source-rock development across the sedimentary sequence. Results presented here document the organic-rich and oil-prone source-rock quality of the Hue Shale in the distal part of the Colville foreland basin and bolster the potential for a Cretaceous petroleum system beneath the eastern North Slope.

Alaska

Precipitation-based flood-inundation maps for the East Fork Little Blue River and tributaries at Lee’s Summit, Missouri, 2024

The U.S. Geological Survey, in cooperation with the City of Lee’s Summit, Missouri, assessed flooding of the East Fork Little Blue River and tributaries for varying precipitation magnitudes and durations, varying antecedent runoff conditions, and projected climate-change conditions. The precipitation scenarios were used to develop a library of flood-inundation maps for a 2.95-mile reach of the East Fork Little Blue River and tributaries within the city. A two-dimensional U.S. Army Corps of Engineers Hydrologic Engineering Center–River Analysis System (HEC–RAS; ver. 6.5) rain-on-grid model was calibrated to selected runoff events representing a range of antecedent runoff conditions and hydrologic responses. Lowest adjacent grades for structures within the nearby study area were incorporated into the terrain, and depth grids and water-surface elevation grids were developed for the study area. Simulated velocities at selected bridge locations were also developed from the model. The model was calibrated using water-surface elevation data collected from water-level loggers (pressure transducers) and streamflow measurements and water-surface elevation measurements made at a reference point during runoff events. The calibrated HEC–RAS model was used to simulate streamflows from design rainfall events of 15-minute to 24-hour durations and ranging from a 100- to 0.1-percent annual exceedance probability (1-year to 1,000-year recurrence intervals). Flood-inundation maps were produced for depths at a reference location of 3 to 16 feet, or a depth exceeding the 0.1-percent annual exceedance probability interval precipitation. The results of each precipitation duration-frequency value were represented by a 1-foot-increment inundation map based on the generated peak streamflow from that rainfall event and the corresponding water-surface elevation at the East Fork Little Blue River reference location. Within the HEC–RAS model, 240 scenarios were developed from the design rainfall events with each of 3 antecedent conditions. Additional scenarios were created to simulate the effects of projected precipitation scenarios on the 100-year recurrence interval, 24-hour storm and the 100-year recurrence interval, 6-hour storm. All simulation results were assigned to a flood-inundation map condition based on the generated peak flow and corresponding water-surface elevation at the East Fork Little Blue River reference location. The flood-inundation maps are shown on a web mapping application made available to the public through the City of Lee’s Summit (hyperlink will be added when available). The flood-inundation maps are tied to real-time precipitation data obtained from the Automated Surface Observing System weather station at the Lee’s Summit Municipal Airport, accessible at https://mesonet.agron.iastate.edu/request/download.phtml?network=MO_ASOS . The availability of these maps, along with information regarding observed rainfall, could help provide emergency management personnel and residents with information that is critical for flood-response activities, such as evacuations and road closures, and for postflood recovery efforts.

Missouri