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 46Linked to original sources

Diel temperature signals track seasonal shifts in localized groundwater contributions to headwater streamflow generation at network scale

Groundwater contributions to streamflow sustain aquatic ecosystem resilience; streams without significant groundwater inputs often have well-coupled air and water temperatures that degrade cold-water habitat during warm low flow periods. Widespread uncertainty in stream-groundwater connectivity across space and time has created disparate predictions of energy and nutrient fluxes across headwater networks, hindering predictions of cold-water habitat resilience under climate change scenarios. Recently, annual paired air and water temperature signals have been harnessed to indicate stream water thermal sensitivity and the dominance of deep versus shallow groundwater influence, although the utility of diel air–water temperature signal metrics for hydrologic inference has remained unexplored. Here we analyzed two consecutive years of locally paired, air–water temperature data from 47 headwater stream sites in the Catskill Mountains, New York, USA, and discovered characteristic seasonal patterns in diel temperature signal sinusoid metrics (amplitude ratio, phase lag, and mean ratio) driven by shifts in streamflow generation mechanisms and stream network position. Hydrologic interpretations of observed patterns were supported by stream heat budget model scenarios and additional analysis of paired air–water temperature data from two streams in Shenandoah National Park, Virginia, USA, with well characterized stream-groundwater connectivity. We found that within smaller tributaries, streamflow generation transitions from runoff to groundwater dominance were driven by hillslope drying during seasonal periods of lower precipitation. This was evidenced by significant correlations (p < 0.01) between daily water:air temperature signal amplitudes (non-linear decreases of ∼ 50 %) and derived base-flow index at 22 of the 28 sites, indicating enhanced local groundwater influence on streamflow promotes decoupling of diel air–water temperature signals. Additionally, ratios between daily water:air temperature signal means were lower in tributaries (∼0.68) when compared to main-stem (∼0.8) sites, increasing linearly throughout the observational period. In conceptual stream heat budget models, groundwater inflow had minimal effects on daily phase lags (∼0.2 hr), but increases in fractional groundwater discharge (0–50 %) depressed daily amplitude (∼20 % to 50 %) and mean ratios (∼15 %), supporting the sensitivity of daily metrics to interpreted changes in seasonal groundwater contributions to streamflow. During observational periods (i.e., April through October 2021 and 2022), significant differences (p < 0.01) between tributary and main-stem air–water metrics occurred when base-flow contributions were highest (∼0.93 vs. ∼ 0.68), as sites lower in the network had daily temperature metrics dominated by stream channel thermal inertia, rather than local groundwater connectivity, showing enhanced air–water diel signal coupling during warmer, drier periods. Divergent air temperature coupling across the network was interpreted as being driven by distance from local groundwater source zones, additional lateral groundwater inflows do not contribute a meaningful fraction to channel discharge lower in the network. Given the growing footprint of stream temperature observations, diel air–water temperature signals can provide distributed metrics sensitive to upstream groundwater discharge. Consequently, these metrics can support ongoing efforts by resource managers and researchers seeking to forecast the resilience of cold-water habitat to climate warming and changing precipitation regimes in mountain headwater streams.

New York, Virginia

Slope Unit Maker (SUMak): An efficient and parameter-free algorithm for delineating slope units to improve landslide modeling

Slope units are terrain partitions bounded by drainage and divide lines. In landslide modeling, including susceptibility modeling and event-specific modeling of landslide occurrence, slope units provide several advantages over gridded units, such as better capturing terrain geometry, improved incorporation of geospatial landslide-occurrence data in different formats (e.g., point and polygon), and better accommodating the varying data accuracy and precision in landslide inventories. However, the use of slope units in regional ( > 100 km 2 ) landslide studies remains limited due, in part, to the large computational costs and/or poor reproducibility with current delineation methods. We introduce a computationally efficient algorithm for the parameter-free delineation of slope units that leverages tools from within TauDEM and GRASS, using an R interface. The algorithm uses geomorphic laws to define the appropriate scaling of the slope units representative of hillslope processes, avoiding the often ambiguous determination of slope unit size. We then demonstrate how slope units enable more robust regional-scale landslide susceptibility and event-specific landslide occurrence maps.

Natural Hazards and Earth Systems Sciences (NHESS)

A glimpse into the future of tectonic tremor monitoring

Tectonic tremor is a weak, long-duration seismic signal often observed in subduction zones and on some other plate-bounding faults. Because of tremor's characteristically low amplitude (and low signal-to-noise) and lack of clear phase arrivals, detecting and locating tremor usually requires techniques distinct from those applied to typical earthquakes. Major advances in detection and understanding of tremor have derived in the past from a powerful combination of new data and new analysis techniques. In a recent study, Sagae et al. (2025, https://doi.org/10.1029/2025jb031348 ) exploit that combination again, developing a new machine-learning based workflow and applying it to the S-net cabled seismic network in the Japan trench offshore northern Honshu. Their approach, although complex, succeeds in detecting several times more tremor activity than earlier studies, resulting in new insights and providing a blueprint for similar approaches that could be applied elsewhere. As real-time earthquake monitoring adopts similar tools, it may present an opportunity to bring tremor monitoring into operational workflows. In turn, this could solidify tremor monitoring as a component of future operational earthquake forecasting.

JGR Solid Earth

Groundwater-storage change in the north Phoenix aquifer, Arizona, 2020–23

The city of Phoenix, Arizona, relies primarily on surface water for municipal water supply. The city also maintains wells to withdraw groundwater, particularly in times of drought and reduced surface-water supply, and to recharge groundwater when excess surface water is available. As of 2023, withdrawals from the aquifer in the northeastern part of the city are a small volume of water, less than 3,000 acre-feet in most years. Each year a similar volume of water is recharged through injection wells. The withdrawal wells are permitted to produce in total more than 20,000 acre-feet per year; increased pumping could lead to future groundwater-storage declines. To better understand groundwater-storage change in the north Phoenix aquifer, a repeat microgravity monitoring network was established in 2020. Measurements of changes in Earth’s gravity provide a direct, non-invasive measurement of subsurface mass change. Groundwater-storage changes were small during the 2020–23 study period, consistent with the relatively small volumes of pumping and recharge. Groundwater levels measured in monitoring wells were stable or increased slightly during this period, although the number of monitoring wells within the area of the gravity monitoring network is sparse. In total, about 15,000 acre-feet of water were pumped and 31,000 acre-feet recharged through injection in the north Phoenix aquifer during the 2020–22 period. Within the monitored area, groundwater storage increased by a small amount, about 1,000 acre-feet, in 2020, and decreased by a larger amount, about 6,000 acre-feet, each year in 2021 and 2022. Storage decreased at 89 of 102 stations from 2020 to 2023. Groundwater-storage decreases greater than the volume of net pumping indicate down-gradient subflow from the aquifer is greater than recharge plus incoming subflow, drying of the unsaturated zone resulting from decreased land-surface recharge, or both. At present (2023), the aquifer appears able to store and supply the relatively small amounts of water needed without excessive drawdown or harmful effects, such as land subsidence and surface inundation. If pumping or recharge increases in the future, the established repeat microgravity network is well suited to capture the spatial extent and magnitude of groundwater-storage changes in the aquifer.

Arizona

An early Holocene wet period in the southwestern United States

Multiple generations of spring-fed streams traversed ∼800 km 2 of the Las Vegas Valley in southern Nevada between ca. 10.9 ka and 8.5 ka, depositing an extensive tufa network. The scale of this network and diversity of tufa morphologies is novel in North America and offers an opportunity to obtain quantitative paleoclimate data for the region during the early Holocene. We determined isotopic compositions and estimated past temperatures using clumped isotope data from early Holocene tufa on the valley floor (698 m) as well as tufa forming today at higher elevation in the nearby Spring Mountains at Cold Creek Spring (1856 m). Modern and fossil tufa yielded comparably low δ 18 O values, implying that source waters for both were derived from high-elevation winter precipitation. Clumped isotope temperatures of modern tufa average 15.8 ± 2.5 °C, aligning with mean summer temperatures of the emergent spring water, and indicate equilibrium conditions of tufa formation. The early Holocene tufa yielded similar clumped isotope temperatures, averaging 15.2 ± 3.9 °C, meaning it precipitated at temperatures that occur at much higher elevations today. The Las Vegas tufa record, combined with nearby and temporally correlative paleospring and lacustrine records, suggest that cool/wet conditions prevailed throughout the Mojave Desert during the early Holocene. These records also demonstrate that spring ecosystems responded to millennial-scale hydroclimate variations that supersede climate change driven solely by insolation. The previously unrecognized pattern of ecosystem response to hydroclimate documented here may assist in understanding climate drivers for the early Holocene and provide critical information for the fate of groundwater-dependent ecosystems in the southwestern United States.

Nevada

Groundwater flow model for the Des Moines River alluvial aquifer near Des Moines, Iowa

Des Moines Water Works (DMWW) is a regional municipal water utility that provides residential and commercial water resources to about 600,000 customers in Des Moines, Iowa, and surrounding municipalities in central Iowa. DMWW has identified a need for increased water supply and is exploring the potential for expanding groundwater production capabilities in the Des Moines River alluvial aquifer, where it operates two radial collector wells (RCWs). The U.S. Geological Survey, in cooperation with DMWW, completed a study of the Des Moines River alluvial aquifer and interactions of the RCWs with the aquifer; no previously published model has included the existing well locations, which is the focus of this model. A conceptual and numerical groundwater flow model have been developed to characterize the Des Moines River alluvial aquifer under existing conditions, to simulate water levels observed in the RCWs, and to provide publicly accessible hydrologic data and research that advance understanding of the regional hydrologic system and can potentially be used in the future to evaluate groundwater production scenarios. Model performance was assessed by comparing observed and simulated groundwater levels that included water level elevations, water level changes, water level inequality observations, surface water streamflow, and change in surface water volume from upstream to downstream. Water table elevation in the aquifer layers is on average slightly overestimated with average absolute value error less than 1.5 meters at both RCWs and less than 2.5 meters for all observation wells in the alluvial aquifer layers. The model also accurately simulated water tables greater than the RCW design minimum (a water level threshold at which RCW pumping is reduced) in all timesteps for which water level observation data existed. Water table elevation error was higher in other model layers that were not the focus of the study, and the model did not accurately match streamflow targets.

Iowa

High frequency and region-scale simulations of large (Mw7+) earthquakes on the southern Whidbey Island fault, Washington, USA

We simulate ground shaking in western Washington State from hypothetical M w 7.0–7.5 earthquakes on the southern Whidbey Island fault (SWIF). Ground motions are modeled considering kinematic source distributions on a complex fault plane, a 3D seismic velocity model, and region‐specific soil velocity models. We run simulations with varying model resolutions, including regional‐scale simulations with a maximum‐modeled frequency of ∼1 Hz and local‐scale simulations with a maximum‐modeled frequency of ∼2.5 Hz. Additional local‐scale simulations are run considering high‐resolution surface topography. We explore how source parameters (i.e., magnitude, hypocenter location, and dip direction) and 3D velocity structure impact peak shaking intensity and its variability. In particular, we find that earthquakes on the SWIF would likely produce strong shaking throughout the populated Puget Lowland, including in the cities of Everett, Seattle, Bellevue, and Tacoma, Washington. Simulated short‐period ( T ≤ 2 s) spectral accelerations are strong throughout the Puget Lowland, and long‐period shaking ( T ≥ 5 s) is strong in the deep regional sedimentary basins, especially the Everett and Seattle basins. Source parameters strongly influence intra‐ and interevent variability in response, primarily through changes in source and site geometry, as well as rupture directivity. We also note a potential coupling between rupture directivity and basin effects, wherein directivity pulses are seemingly guided through the region’s deep, interconnected sedimentary basins. Overall, this work highlights the impacts of 3D source, path, and site effects on seismic hazard in the U.S. Pacific Northwest and substantially expands the catalog of simulated ground motions for Puget Sound area crustal faults.

Washington

Approximate inland extent of saltwater intrusion at the base of the Biscayne aquifer, Miami-Dade County, Florida, 2022

Miami-Dade County is part of a densely populated urban corridor in southeastern Florida. The Biscayne aquifer serves as Miami-Dade County’s primary drinking water source and is characterized by highly permeable karstic limestone and carbonate sand. The aquifer’s coastal location and permeable nature make it susceptible to saltwater intrusion. Monitoring the current inland extent and the rate of movement of the saltwater front in the aquifer can inform management strategies for conserving the long-term sustainability of the county’s water supply. In the 1950s, the U.S. Geological Survey published a map of the inland extent of saltwater intrusion in the Biscayne aquifer and has continued to update this map to monitor changes over time, with the most recent update published in 2018. An updated map has been created showing the approximate inland extent of saltwater intrusion in the Biscayne aquifer in eastern Miami-Dade County in 2022, with the 2018 extent shown for comparison. The inland extent of saltwater intrusion was mapped through the interpretation of borehole electromagnetic induction logs and measurements of chloride and specific conductance in groundwater samples. The location of the saltwater interface at the base of the Biscayne aquifer was represented by the 1,000-milligram-per-liter isochlor. This report describes changes in the location of the saltwater interface from 2018 to 2022. By 2022, the saltwater interface had moved farther inland in both the northern and southern parts of the county, advancing by as much as 0.3 kilometer in the north and up to 0.8 kilometer in the Model Land Area to the south. However, it remained relatively unchanged from its 2018 position in the east-central part of the county.

Florida

Colored shaded-relief bathymetric and acoustic-backscatter maps of Jenkinson Lake with orthomosaic of the Sly Park Creek and Hazel Creek area, California

The Caldor Fire was ignited on August 14, 2021, and burned almost 222,000 acres (898 square kilometers) in forested terrain of the central and western Sierra Nevada, California. During the subsequent two months, the fire burned nearly all of Sly Park Creek watershed in El Dorado County. The El Dorado Irrigation District manages the water supply for the area using storage in Jenkinson Lake, a 1.6-kilometer- (1.0-mile-) wide and 3.6-kilometer- (2.2-mile-) long reservoir, located south of the town of Pollock Pines. Several weeks after the fire, the U.S. Geological Survey began investigations into post-fire landscape responses, including sediment yield, by measuring new sediment deposition in Jenkinson Lake. This study focused on the collection and processing of bathymetric and acoustic-backscatter data, as well as onshore aerial imagery in and around Jenkinson Lake, to support wildfire science after the Caldor Fire. A colored shaded-relief bathymetric map (sheet 1) and an acoustic backscatter map (sheet 2) show the lake floor morphology and backscatter intensities.

California

Monitoring and assessment of urban stormwater best management practices at selected Chicago public schools in Chicago, Illinois, from September 1, 2016, to July 1, 2017

The Space to Grow program helps transform aging and neglected schoolyards of Chicago Public Schools into outdoor community spaces with the goal of promoting health and learning while addressing neighborhood flooding issues. Virgil I. Grissom Elementary School and Donald L. Morrill Math and Science School were selected in 2014 for schoolyard upgrades and the installation of various green infrastructure (GI) improvements. The U.S. Geological Survey installed sensors to measure precipitation, groundwater levels, and stormwater runoff volumes from September 1, 2016, to July 1, 2017. At Virgil I. Grissom Elementary School, about 933,000 gallons of water fell on the schoolyard during the monitoring period. No discharge was recorded coming from the GI sewer lines, but backflow indicated water was flowing from the sewer line draining the impervious running track into the combined manhole structure and backwards into the GI retention basins (as designed). This design allowed for a 100-percent capture rate. Native soil at Virgil I. Grissom Elementary School also was conducive to rapid infiltration. Soil borings at Virgil I. Grissom Elementary School indicated about 10.5 feet (ft) of fine sand overlying silty clay to a depth of at least 16 ft. At Donald L. Morrill Math and Science School, about 1,120,000 gallons of water fell on the schoolyard during the monitoring period. About 72.5 precent of this water was discharged into the sewer system, and the other 27.5 percent was captured by the GI. Unlike Virgil I. Grissom Elementary School, the soil profile at Donald L. Morrill Math and Science School consisted of about 5 ft of clay loam overlying stiff blue clay to a depth of at least 12 ft. The sewer line coming from the GI under the football field was at the bottom of the reservoir. This design seemed to allow water to flow out of the line before being absorbed by the retention basin.

Illinois

Water-resources inventory and assessment at Katahdin Woods and Waters National Monument

The U.S. Geological Survey, in cooperation with the National Park Service, prepared a water-resources inventory and assessment for Katahdin Woods and Waters National Monument (KAWW). This compilation includes published and publicly accessible hydrologic data and resource assessments of streams, rivers, ponds, lakes, wetlands, vernal pools, and groundwater in and near KAWW. It also includes reports and datasets summarizing attributes of KAWW’s hydrologic infrastructure, such as stream crossings, dams, wastewater discharge plants, groundwater monitoring wells, and U.S. Geological Survey streamflow-gaging stations. Descriptions of data and details of current limitations in available datasets are included. Wetland, groundwater, streamflow, and water-quality information are all limited. Hydrography data are available; however, there are limited ground-truth data. Accurate streamlines within KAWW were developed from light detection and ranging (lidar) as a part of this work. Hydrologic infrastructure information is available from multiple sources; however, differences exist among the datasets. Datasets are summarized in appendix 1.

Maine

Mechanics and statistics of postseismic shaking

Analysis of two weeks of continuous post-seismic shaking after the 2019 M7.1 Ridgecrest, CA earthquake sequence using 4 nearby borehole seismometers reveals that continuous ground motions decay as Omori’s law in time and follow the Gutenberg-Richter distribution in logarithmic amplitude. The measured temporal decay in amplitudes agrees with predictions of the rate-and-state framework and indicates shaking amplitudes are proportional to the velocity of afterslip. Our ground motion-based statistical framework provides a basis to forecast shaking intensity in the minutes to hours after a large earthquake.

California

Transport and dispersion of fluorescent tracer particles for the dune-bed condition, Atrisco Feeder Canal near Bernalillo, New Mexico

A tracer technique in which mineral particles were coated with fluorescent dyes was used to study the rates of transport and dispersion of sediment particles`of various diameters and specific gravities for a dune-bed condition in an alluvial channel. The experiment was conducted in the Atrisco Feeder Canal near Bernalillo, N. Mex., between May 1 and July 14, 1967. A continuous point source of tracers, approximated by injections at 5- or 10-minute intervals, was maintained for 7 days so that the steady-dilution procedure could be used to calculate the transport rate of bed material. After termination of the injection process, the spatial-integration procedure was used to follow the movement of the tracers downstream and to calculate the transport rate. Samples of the bed material in transport and the accompanying tracers moving along the surface of the dune bed were obtained periodically throughout the study with the "dustpan" sampler especially designed for fluorescent tracer studies. In addition, the spatial distributions of the tracers in the dune bed were determined three times during the study by core sampling.

New Mexico

Slow slip detectability in seafloor pressure records offshore Alaska

In subduction zones worldwide, seafloor pressure data are used to observe tectonic deformation, particularly from megathrust earthquakes and slow slip events (SSEs). However, such measurements are also sensitive to oceanographic circulation-generated pressures over a range of frequencies that conflate with tectonic signals of interest. Using seafloor pressure and temperature data from the Alaska Amphibious Community Seismic Experiment, and sea surface height data from satellite altimetry, we evaluate the efficacy of various seasonal and oceanographic pressure signal proxy corrections and conduct synthetic tests to determine their impact on the timing and amplitude prediction of ramp-like signals typical of SSEs. We find that subtracting out the first mode of the complex empirical orthogonal functions of the pressure records on either the shelf or slope yields signal root-mean-square error (RMS) reductions up to 73% or 80%, respectively. Additional correction with proxies that exploit the depth-dependent spatial coherence of pressure records provides cumulative variance reductions up to 83% and 93%, respectively. Our detectability tests show that the timing and amplitude of synthetic SSE-like ramps can be well constrained for ramp amplitudes ≥4 cm on the shelf and ≥2 cm on the slope, using a fully automated detector. The principal limits on detectability are residual abrupt changes in pressure that occur as part of the transition to and from summer to winter conditions but are not adequately characterized by our seasonal corrections, as well as the inability to properly account for instrumental drift, which is not readily separated from the seasonal signal.

Alaska

Greater sage-grouse habitat of Nevada and northeastern California—Integrating space use, habitat selection, and survival indices to guide areas for habitat management

Executive Summary Greater sage-grouse populations ( Centrocercus urophasianus ; hereafter sage-grouse) are threatened by a suite of disturbances and anthropogenic factors that have contributed to a net loss of sagebrush-dominant shrub cover in recent decades. Declines in sage-grouse populations are largely linked to habitat loss across their range. A key component of conservation and land use planning efforts for sage-grouse involves the continued monitoring and modeling of habitat requirements and suitability across its range. The Bureau of Land Management (BLM) is addressing the management of sage-grouse habitats on BLM-authorized public lands throughout the western United States through a land use planning amendment and associated environmental impact statement (86 FR 66331). More than 25 percent of the range-wide distribution of sage-grouse is within Nevada and northeastern California, and information on sage-grouse distribution and habitat requirements is important to guide appropriate management decisions. Therefore, the BLM has identified the need for updated spatially explicit information on sage-grouse habitat in Nevada and northeastern California to guide the land use planning amendment and associated management decisions. To address this need, researchers with the U.S. Geological Survey, in close cooperation with multiple State and Federal resource agency partners, including BLM, Nevada Department of Wildlife (NDOW) and California Department of Fish and Wildlife (CDFW), sought to map sage-grouse distribution and produce example habitat designations in these states. Herein, we report results of our primary study objective, which was to map sage-grouse habitat and create example habitat management areas, based on more than a decade of location and survival data collected from marked sage-grouse across the study region coupled with lek count survey data managed by the NDOW and the CDFW. We expanded on previously developed methodology to incorporate information on habitat selection and survival during reproductive life stages and specific seasons with updated sage-grouse location and known fate datasets, while also including brood-rearing areas that are understood to be threatened and important for population persistence. We combined predictive habitat map surfaces for each life stage and season with updated information on current occupancy patterns to classify habitat based on its suitability and probability of occupancy. We carried out additional steps to delineate specific example habitat management areas, specifically (1) incorporated corridors connecting key nesting and brood-rearing habitat, (2) corrected outputs for pre-wildfire habitat conditions within areas burned in the last 16 years, and (3) masked out areas of anthropogenic development. Our methodological example of deriving habitat management areas was intended to help inform decisions by BLM and other land managers regarding conservation and management of sage-grouse. Associated data products in the form of habitat maps provide updated, detailed, and comprehensive information about the status of habitats and can be useful to partner agencies in their efforts to designate and rank habitats for this species of high conservation concern in Nevada and California, with full recognition that on-the-ground field data and local sources of information and expertise should be used in conjunction with inferences from these models.

California, Idaho, Nevada

Geohydrology of the Cross-Florida Barge Canal area, with special reference to the Ocala vicinity

The Cross-Florida Barge Canal route commences at Palatka on the St. Johns River, about 75 miles upstream from the Atlantic Ocean, and extends 110 miles southwestward across Peninsular Florida into deep water in the Gulf of Mexico near Yankeetown. The canal will be equipped with five locks, each 600 feet long and 84 feet wide, and the channel will be a minimum of 12 feet deep and 150 feet wide. From near Ocala northeastward, the canal channel will replace much of the natural channel of the Oklawaha River, and will be excavated into beds of the so-called shallow sand aquifer of Miocene age and younger, which overlies limestone of the Floridan aquifer. Westward from Ocala, most of the canal will be excavated below the potentiometric surface into limestone and dolomite of the Floridan aquifer. Water levels of Rodman, Eureka, and Inglis Pools will be controlled by dams and spillways with the limited exchange of water between the pools and the aquifers. The water levels in the Summit Pools will fluctuate with the natural changes in the ground-water level of the Floridan aquifer, although the stage of the pool will be controlled partly by the stage held in the Eureka Pool. A dynamic inflow-outflow relationship will exist between the Summit Pool and the Floridan aquifer. The Floridan aquifer in the canal area is 1,000 to 1,200 feet thick and consists of limestone and dolomite of middle Eocene Miocene age, including from older to younger, the Lake City, Avon Park, and Ocala limestones plus permeable sandy, dolomitic limestone in the lower part of the Hawthorn Formation. It is possible that most of the flow to the two major springs in the area occurs in the upper 100 feet or so of the aquifer in the Ocala Limestone. The aquifer is underlain by the Oldsmar limestone of early Eocene age and is overlain by sand, clayey sand, clay and shell beds of Miocene through Holocene age, in thickness from a few feet to 300 feet. The permeable beds overlying the Floridan aquifer constitute the shallow aquifer, while the poorly permeable ones act as confining beds where the Floridan aquifer is under artesian conditions. A north-south line drawn separating the head of Silver Springs on the west from the Oklawaha River on the east marks the approximate western limit of a continuous blanket of materials of Miocene-Pliocene(?) age covering the rocks of the Floridan aquifer. East of the line, much of the aquifer is under artesian conditions, particularly in the Oklawaha River valley, although in some areas east of the valley, direct recharge through thick permeable Miocene-Pliocene(?) sands occurs. West of the line, only scattered remnants of a once continuous Miocene-Pliocene(?) cover remain. Lack of the cover is a result of erosion on the crest and flank of the Ocala Uplift, a broad northwest-southeast trending anticlinal upwarp, the axis of which is crossed by the canal route in the Dunnellon area. Over most of this area the Floridan aquifer is unconfined and receives direct recharge through a cover of a few tens of feet of sand and clayey sand of Quaternary age. Tensional stresses during the structural evolution of the Ocala Uplift produced an intersecting system of fractures and normal faults in rocks of the Floridan aquifer. The fractures and faults are important controls for orientation of solution channels and, therefore, for development of ground-water circulation patterns. When the system surface streams, which once drained the Barge Canal area, eroded the poorly permeable Miocene-Pliocene(?) cover from the flanks of the Ocala Uplift, surface runoff was reduced and precipitation began to directly infiltrate the underlying limestones. Now only principal streams remain, such as the Oklawaha and Withlacoochee Rivers and a few short tributaries, while one of the most highly developed subsurface drainage systems in the world has evolved in cavernous limestone of the Floridan aquifer. Two of the larger freshwater springs in the world now discharge from the Floridan aquifer in the canal area. Silver Springs near Ocala discharges an average 531 mgd (million gallons per day) down the 4-mile long Silver River, which flows on poorly permeable beds to the Oklawaha River. Rainbow Springs near Dunnellon discharges on average 468 mgd from numerous orifices in the bed of the 5-mile-long Rainbow River, which flows into the Withlacoochee River. The heads of the springs have migrated to their present positions partly because of a tendency of ground-water levels to decline as permeability in the aquifer is increased due to removal of limestone by solution, and because of mechanical erosion of the limestone in the vicinity of the spring heads. Also, points of principal spring discharges have shifted in the past due to changes in ground-water levels in response to changes in sea level. The subsurface drainage system is continuing to evolve today, as evidenced in part by frequent occurrence of new sinkholes and by the presence of significant amounts of calcium bicarbonate in the spring waters. Rodman Pool, at the east end of the canal, is separated from the Floridan aquifer by poorly permeable materials. The pool's operating water level will be only a few feet above the potentiometric surface at the downstream end, and at or slightly below the potentiometric surface at the upstream end. Little exchange of water between the Rodman Pool and the Floridan aquifer is expected. Eureka Pool, just upstream from Rodman Pool, will also be separated from the Floridan aquifer by poorly permeable beds. However, the stage of the pool will be about 15 feet higher than the natural potentiometric surface at the pool's downstream end, and some seepage into the Floridan aquifer is anticipated through faults and leaky parts of the poorly permeable beds, with a consequent rise in ground-water levels in areas adjacent to the lower end of the pool. Possibilities for particulate contamination of the aquifer will tend to be minimized because of the filtering capacity of the materials through which water must pass to reach the aquifer, although the natural filter will not preclude movement into the aquifer of contaminants which might become dissolved in the pool waters. No significant interchange of water between the pool and the aquifer is expected at the upstream end of the Eureka Pool. Present construction plans indicate an operating stage for Eureka Pool which will range between 38 and 40 feet above mean sea level, although it is possible to dredge the pool deep enough to permit a range in stage of 36 to 40 feet. A backwater effect extending up Silver River from the Eureka Pool is expected to regulate the stage at the head of Silver Springs between 39 and 44 feet above mean sea level if the pool ranges between 36 and 40 feet. If Eureka Pool ranges only between 38 and 40 feet, the range of stage at the head of the springs should be about 41 to 44 feet above mean sea level. From Inglis Lock west, the canal will have direct connection with the Gulf, and canal stage will fluctuate with the Gulf tide. Since the canal stage will be slightly lower than the adjacent ground-water levels along much of the reach, there will be some ground-water inflow to the canal. No significant changes in the existing ground-water regime are expected in the vicinity of Inglis Pool, the first step up in the canal east of the Gulf. Existing ground-water and surface-water levels in the area will not change appreciable, and the natural stage and flow of Rainbow Springs, which will flow by way of Rainbow River into Inglis Pool, should not be affected by canal operations. A possible adverse effect of the canal on the Inglis Pool area could result if sea water is locked up from the Gulf through Inglis Lock. However, the high step of 25 feet at the lock, flushing action of continuous flow from Inglis Pool to the lower reaches of the Withlacoochee River, and use of possible preventive locking procedures should minimize the problem. The potential for adverse effects on the ground-water regime is greatest in the area of the Summit Pool. Through most of the length of the pool, the canal channel will be excavated into limestone of the Floridan aquifer to depths of 12 to 27 feet below the potentiometric surface. Changes that will take place in the ground-water flow system in the Silver Springs drainage area, once the canal is completed, were estimated by flow-net analysis. Variation in aquifer transmissivity was determined by calculating transmissivity in 25 different flow cells surrounding Silver Springs. Transmissivity in the 25 cells averages about 15,600,000 gpd/ft [2,090,000 ft2/day (feet squared)], but transmissivity in the six cells through which the Summit Pool passes ranges from 9,000,000 to 44,000,000 gpd/ft (1,210,000 to 5,900,000 ft2/day). Transmissivity was used to compute static stage of the Summit Pool under given ground-water level conditions. Had the canal existed in May 1968 and had the stage of Eureka Pool been held at 36 feet at the time, the static stage in Summit Pool would have been about 42.1 feet above mean sea level. Thus, a conceptual model of the changes in the potentiometric surface wrought by the finished canal was drawn, and zones of ground-water inflow and outflow were delineated. Most outflow from the Summit Pool to the aquifer should be limited to one 4-mile-long zone along the north side of the pool, about 5 miles south of Silver Springs. It is estimated that a water volume equivalent to about 8 percent of the daily flow of Silver Springs will enter the Summit Pool each day from the southern one-third of the Silver Springs drainage area. A like amount will reenter the aquifer at the main zone of outflow and move toward Silver Springs at an estimated average velocity of about 200 feet per day, if something close to the natural static stage of the pool is maintained by return pumpage of the lockage losses. At a velocity of 200 feet per day, water from the Summit Pool would discharge at Silver Springs about 140 days later. However, any estimate of velocity in the highly cavernous limestone aquifer in the area should be used with caution, because difficult to measure changes in porosity and thickness of major zones of flow may cause large variations in velocity. If all lockage losses are returned to the Summit Pool by pumping from Eureka Pool, no net loss from the Silver Springs drainage area, except for some evaporation from the water surface in the canal and possible leakage around locks, will result from canal operations. The zone of outflow from the Summit Pool to the aquifer will be in a natural potentiometric trough, and the zone of inflow will be in a potentiometric ridge area. The equilibrium water level in the pool will tend to be about 1 foot higher than the altitude of the lowest level in the pre-canal potentiometric trough, and about 2 feet lower than the highest level on the pre-canal potentiometric ridge. West of the Silver Springs drainage area just east of Dunnellon Lock, in the area of a local potentiometric high, the water level in the Summit Pool is expected to be about 15 feet below the natural potentiometric surface. In most areas 2 to 3 miles away from the Summit Pool, effects of the canal on the natural potentiometric surface should be slight. The stage of the Summit Pool, judging from the 36-year record for ground-water level changes and the anticipated indirect effect of the controlled stage in Eureka Pool, should have a maximum of about 10.5 feet with a maximum water level of about 51.5 feet above mean sea level and a minimum of about 41.0 feet above mean sea level. Of particular importance in the Summit Pool is an implementation of well planned construction and operational procedures designed to minimize risks of ground-water contamination.

Florida

Declining groundwater storage expected to amplify mountain streamflow reductions in a warmer world

Groundwater interactions with mountain streams are often simplified in model projections, potentially leading to inaccurate estimates of streamflow response to climate change. Here, using a high-resolution, integrated hydrological model extending 400 m into the subsurface, we find groundwater an important and stable source of historical streamflow in a mountainous watershed of the Colorado River. In a warmer climate, increased forest water use is predicted to reduce groundwater recharge resulting in groundwater storage loss. Losses are expected to be most severe during dry years and cannot recover to historical levels even during simulated wet periods. Groundwater depletion substantially reduces annual streamflow with intermittent conditions predicted when precipitation is low. Expanding results across the region suggests groundwater declines will be highest in the Colorado Headwater and Gunnison basins. Our research highlights the tight coupling of vegetation and groundwater dynamics and that excluding explicit groundwater response to warming may underestimate future reductions in mountain streamflow.

Nature Water

Comparison of hydrologic data and water budgets between 2003–08 and 2018–23 for the eastern part of the Arbuckle-Simpson aquifer, south-central Oklahoma

The Arbuckle-Simpson aquifer is divided spatially into three parts (eastern, central, and western). The largest groundwater withdrawals are from the eastern part of the Arbuckle-Simpson aquifer, which provides water to approximately 39,000 people in Ada and Sulphur, Oklahoma, and surrounding areas. The Arbuckle-Simpson aquifer, including the eastern part, is designated a sole source aquifer for its service area. Based primarily on data collected between 2003 and 2008, a series of comprehensive hydrologic studies of the Arbuckle-Simpson aquifer was published to provide the information necessary to perform groundwater-flow model simulations so that the Oklahoma Water Resources Board could determine how much water could be withdrawn from the aquifer while maintaining flow to springs and streams. As part of the Phase 1 studies, an aquifer water budget was developed from a numerical model for the period 2003–08. For this report, Phase 1 refers to the 2003–08 data collection period, although for some of the analyses, data collected prior to 2003 were used to inform model development work. Allocation of water from this aquifer was then established by the Oklahoma Water Resources Board in 2013. Additional well-spacing rules were also established by the Oklahoma Water Resources Board for sensitive sole source groundwater basins. To determine how the water budget for the eastern part of the Arbuckle-Simpson aquifer has changed over time, recently collected hydrologic data (2018–23) were compared to data collected during 2003–08. The analysis of changes in the aquifer water budget from 2003–08 to 2018–23 could help resource managers better understand changes in the overall balance of water in storage and the potential effects on streamflow, changes in groundwater levels, and the effects of different water uses in the aquifer area on available water in the eastern part of the Arbuckle-Simpson aquifer and streams overlying the eastern part of the Arbuckle-Simpson aquifer.

Oklahoma