USGS Science⌕ Search

SEARCH · USGS Science

Results for “Hydrological Processes”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,099 records · Page 61Linked to original sources
USGS1989

Hydrologic and geochemical monitoring in Long Valley caldera, Mono County, California, 1986

The U.S. Geological Survey continued to monitor hydrologic and geochemical conditions in the Long Valley caldera during 1986. The monitoring is directed toward detecting changes in the hydrologic system caused by tectonic or magmatic processes. Data collected during 1986 include chemical and isotopic composition of water from selected streams sites, springs, and wells; pumpage from four geothermal wells; flow rates of selected springs and stream sites; mean daily water or gas temperatures at selected sites; mean daily atmospheric pressures and water level at selected wells, and precipitation records for two sites. Seismicity within the caldera persisted at a relatively low level compared with the more active periods of 1978-84. The most significant events of seismicity that affected hydrologic monitoring sites in Long Valley during 1986 occurred during July , in response to the Chalfant Valley earthquakes, centered about 20 miles southeast of the caldera. Water level records for three wells show distinct responses to the Chalfant Valley earthquakes. (USGS)

Water-Resources Investigations Report↗
USGS2012

Modeling thermal dynamics of active layer soils and near-surface permafrost using a fully coupled water and heat transport model

Thawing and freezing processes are key components in permafrost dynamics, and these processes play an important role in regulating the hydrological and carbon cycles in the northern high latitudes. In the present study, we apply a well-developed soil thermal model that fully couples heat and water transport, to simulate the thawing and freezing processes at daily time steps across multiple sites that vary with vegetation cover, disturbance history, and climate. The model performance was evaluated by comparing modeled and measured soil temperatures at different depths. We use the model to explore the influence of climate, fire disturbance, and topography (north- and south-facing slopes) on soil thermal dynamics. Modeled soil temperatures agree well with measured values for both boreal forest and tundra ecosystems at the site level. Combustion of organic-soil horizons during wildfire alters the surface energy balance and increases the downward heat flux through the soil profile, resulting in the warming and thawing of near-surface permafrost. A projection of 21st century permafrost dynamics indicates that as the climate warms, active layer thickness will likely increase to more than 3 meters in the boreal forest site and deeper than one meter in the tundra site. Results from this coupled heat-water modeling approach represent faster thaw rates than previously simulated in other studies. We conclude that the discussed soil thermal model is able to well simulate the permafrost dynamics and could be used as a tool to analyze the influence of climate change and wildfire disturbance on permafrost thawing.

Journal of Geophysical Research D: Atmospheres↗
USGS2013

Hydrologic connectivity of floodplains, northern Missouri: implications for management and restoration of floodplain forest communities in disturbed landscapes

Hydrologic connectivity between the channel and floodplain is thought to be a dominant factor determining floodplain processes and characteristics of floodplain forests. We explored the role of hydrologic connectivity in explaining floodplain forest community composition along streams in northern Missouri, USA. Hydrologic analyses at 20 streamgages (207–5827 km 2 area) document that magnitudes of 2-year return floods increase systematically with increasing drainage area whereas the average annual number and durations of floodplain-connecting events decrease. Flow durations above the active-channel shelf vary little with increasing drainage area, indicating that the active-channel shelf is in quasi-equilibrium with prevailing conditions. The downstream decrease in connectivity is associated with downstream increase in channel incision. These relations at streamflow gaging stations are consistent with regional channel disturbance patterns: channel incision increases downstream, whereas upstream reaches have either not incised or adjusted to incision by forming new equilibrium floodplains. These results provide a framework to explain landscape-scale variations in composition of floodplain forest communities in northern Missouri. Faust ( 2006 ) had tentatively explained increases of flood-dependent tree species, and decreases of species diversity, with a downstream increase in flood magnitude and duration. Because frequency and duration of floodplain-connecting events do not increase downstream, we hypothesize instead that increases in relative abundance of flood-dependent trees at larger drainage area result from increasing size of disturbance patches. Bank-overtopping floods at larger drainage area create large, open, depositional landforms that promoted the regeneration of shade-intolerant species. Higher tree species diversity in floodplains with small drainage areas is associated with non-incised floodplains that are frequently connected to their channels and therefore subject to greater effective hydrologic variability compared with downstream floodplains. Understanding the landscape-scale geomorphic and hydrologic controls on floodplain connectivity provides a basis for more effective management and restoration of floodplain forest communities.

Missouri↗
USGS2016

Lateral and subsurface flows impact arctic coastal plain lake water budgets

Arctic thaw lakes are an important source of water for aquatic ecosystems, wildlife, and humans. Many recent studies have observed changes in Arctic surface waters related to climate warming and permafrost thaw; however, explaining the trends and predicting future responses to warming is difficult without a stronger fundamental understanding of Arctic lake water budgets. By measuring and simulating surface and subsurface hydrologic fluxes, this work quantified the water budgets of three lakes with varying levels of seasonal drainage, and tested the hypothesis that lateral and subsurface flows are a major component of the post-snowmelt water budgets. A water budget focused only on post-snowmelt surface water fluxes (stream discharge, precipitation, and evaporation) could not close the budget for two of three lakes, even when uncertainty in input parameters was rigorously considered using a Monte Carlo approach. The water budgets indicated large, positive residuals, consistent with up to 70% of mid-summer inflows entering lakes from lateral fluxes. Lateral inflows and outflows were simulated based on three processes; supra-permafrost subsurface inflows from basin-edge polygonal ground, and exchange between seasonally drained lakes and their drained margins through runoff and evapotranspiration. Measurements and simulations indicate that rapid subsurface flow through highly conductive flowpaths in the polygonal ground can explain the majority of the inflow. Drained lakes were hydrologically connected to marshy areas on the lake margins, receiving water from runoff following precipitation and losing up to 38% of lake efflux to drained margin evapotranspiration. Lateral fluxes can be a major part of Arctic thaw lake water budgets and a major control on summertime lake water levels. Incorporating these dynamics into models will improve our ability to predict lake volume changes, solute fluxes, and habitat availability in the changing Arctic.

Alaska↗
USGS2026

Groundwater–surface water exchange at a New England salt marsh: The interplay between tidal pumping at multiple frequencies and terrestrial groundwater discharge

Groundwater–surface water exchange plays a fundamental role in the transport of energy and solutes in coastal salt marshes, which are positioned at an intersection between terrestrial and marine hydrologic environments. Hydraulic gradients in shallow marsh sediments continuously change due to tidal water level variations that alternately flood and drain the marsh platform, often superimposed on terrestrial groundwater drainage. We use timeseries data from Vertical Temperature Profiler and Conductivity-Temperature-Depth instruments to investigate the space–time patterns of groundwater–surface water exchange at a hillslope-adjacent fringing salt marsh on the Herring River in Wellfleet, MA. Analyses of data acquired from September to December 2023 demonstrate that time-averaged exchange rates systematically vary from recharge near the tidal channel to saline groundwater discharge towards the upland fringe and systematically decrease during the transition from summer to winter conditions. Exchange rates were continuously modulated by tides, including high-frequency compound and harmonic components generated by non-linear interactions as the astronomical tide propagated into the estuary. The exchange rate sensitivity to tidal forcing decreased with both frequency and distance from the tidal channel, consistent with the lateral diffusion of pore pressure perturbations through the high permeability layer underlying the surficial, low permeability peat. We find that tidal loading of the marsh sediments is a two-stage process consisting of direct, instantaneous loading during inundation and indirect, time delayed loading via lateral diffusion of pore pressure variations in-between inundation intervals. The combination of terrestrial groundwater entering the marsh through the high permeability layer beneath the peat, seawater flooding the marsh platform during inundation, and tidally induced pore pressure variations diffusing into the marsh generates a complex space–time exchange pattern, including bi-directional exchange at tidal periods, with poorly understood but important implications for biogeochemical processes in the marsh sediments and constituent fluxes to the coastal ocean.

Massachusetts↗
USGS2009

Estimating 3D variation in active-layer thickness beneath arctic streams using ground-penetrating radar

We acquired three-dimensional (3D) ground-penetrating radar (GPR) data across three stream sites on the North Slope, AK, in August 2005, to investigate the dependence of thaw depth on channel morphology. Data were migrated with mean velocities derived from multi-offset GPR profiles collected across a stream section within each of the 3D survey areas. GPR data interpretations from the alluvial-lined stream site illustrate greater thaw depths beneath riffle and gravel bar features relative to neighboring pool features. The peat-lined stream sites indicate the opposite; greater thaw depths beneath pools and shallower thaw beneath the connecting runs. Results provide detailed 3D geometry of active-layer thaw depths that can support hydrological studies seeking to quantify transport and biogeochemical processes that occur within the hyporheic zone.

Journal of Hydrology↗
USGS2006

Watershed Data Management (WDM) Database for Salt Creek Streamflow Simulation, DuPage County, Illinois

The U.S. Geological Survey (USGS), in cooperation with DuPage County Department of Engineering, Stormwater Management Division, maintains a database of hourly meteorologic and hydrologic data for use in a near real-time streamflow simulation system, which assists in the management and operation of reservoirs and other flood-control structures in the Salt Creek watershed in DuPage County, Illinois. The majority of the precipitation data are collected from a tipping-bucket rain-gage network located in and near DuPage County. The other meteorologic data (wind speed, solar radiation, air temperature, and dewpoint temperature) are collected at Argonne National Laboratory in Argonne, Illinois. Potential evapotranspiration is computed from the meteorologic data. The hydrologic data (discharge and stage) are collected at USGS streamflow-gaging stations in DuPage County. These data are stored in a Watershed Data Management (WDM) database. This report describes a version of the WDM database that was quality-assured and quality-controlled annually to ensure the datasets were complete and accurate. This version of the WDM database contains data from January 1, 1997, through September 30, 2004, and is named SEP04.WDM. This report provides a record of time periods of poor data for each precipitation dataset and describes methods used to estimate the data for the periods when data were missing, flawed, or snowfall-affected. The precipitation dataset data-filling process was changed in 2001, and both processes are described. The other meteorologic and hydrologic datasets in the database are fully described in the annual U.S. Geological Survey Water Data Report for Illinois and, therefore, are described in less detail than the precipitation datasets in this report.

Open-File Report↗
USGS1983

Hydrologic and geomorphic studies of the Platte River basin

The channels of the Platte River and its major tributaries, the South Platte and North Platte Rivers in Colorado, Wyoming, and Nebraska, have undergone major changes in hydrologic regime and morphology since about 1860, when the water resources of the basin began to be developed for agricultural, municipal, and industrial uses. These water uses have continued to increase with growth in population and land development. Diversion of flow from channels, storage of water in reservoirs, and increased use of ground water have affected the distribution and timing of streamflows and the transport of fluvial sediments. All these factors have contributed to changes in channel geometry and the riverine environment. In 1979, the U.S. Geological Survey began investigations in the Platte River basin to determine the effects of water use on the hydrology and morphology of the Platte River and its major tributaries. These investigations also considered the relationship of hydrologic regime to factors that control or affect the habitat of migratory waterfowl in the Platte River valley. This volume brings together the results of several research studies on historical changes in channel morphology, surface-water hydrology, hydraulic geometry, sediment-transport and bedform processes, ground-water and surface-water relations, stochastic models of streamflow and precipitation, and methods for estimating discharge required to maintain channel width. In each of the studies, data on some segment of the Platte River hydrologic system were collected and interpreted. All the studies are interrelated; together they provide some degree of understanding of regime changes that are occurring. The hydrologic research described in the following chapters will be useful in decision-making pertaining to the management of water resources and migratory waterfowl habitats.

Colorado, Nebraska, Wyoming↗
USGS2017

Vertical datum conversion process for the inland and coastal gage network located in the New England, Mid-Atlantic, and South Atlantic-Gulf hydrologic regions

Datum conversions from the National Geodetic Vertical Datum of 1929 to the North American Vertical Datum of 1988 among inland and coastal gages throughout the hydrologic regions of New England, the Mid-Atlantic, and the South Atlantic-Gulf have implications among river and storm surge forecasting, general commerce, and water-control operations. The process of data conversions may involve the application of a recovered National Geodetic Vertical Datum of 1929–North American Vertical Datum of 1988 offset, a simplistic datum transformation using VDatum or VERTCON software, or a survey, depending on a gaging network datum evaluation, anticipated uncertainties for data use among the cooperative water community, and methods used to derive the conversion. Datum transformations from National Geodetic Vertical Datum of 1929 to North American Vertical Datum of 1988 using VERTCON purport errors of ± 0.13 foot at the 95 percent confidence level among modeled points, claiming more consistency along the east coast. Survey methods involving differential and trigonometric leveling, along with observations using Global Navigation Satellite System technology, afford a variety of approaches to establish or perpetuate a datum during a survey. Uncertainties among leveling approaches are generally < 0.1 foot, and and Global Navigation Satellite System approaches may be categorized with uncertainties of ≤0.1 foot for a Level I quality category and ≥0.1 foot for Level II or III quality categories (defined by the U.S. Geological Survey) by observation and review of experienced practice. The conversion process is initiated with an evaluation of the inland and coastal gage network datum, beginning with altitude datum components and the history of those components queried through the U.S. Geological Survey Groundwater Site Inventory database. Subsequent edits to the Groundwater Site Inventory database may be required and a consensus reached among the U.S. Geological Survey Water Science Centers to identify the outstanding workload categorized as in-office datum transformations or offset applications versus out-of-office survey efforts. Datum conversions or datum establishment for the inland or coastal gaging network should meet datum uncertainty requirements among other Federal agencies. Datum uncertainty requirements are ±0.25 foot for U.S. Army Corps of Engineers water-control or construction projects and ±0.16 foot for Federal Emergency Management Agency field surveys and checkpoint surveys used for mapping. River level forecasts generally are defined as ± 0.10 foot among the National Oceanic and Atmospheric Administration–National Weather Service. Collaboration and communication among the cooperative water community is necessary during a datum conversion or datum change. Datum notification time-change requirements set by the National Oceanic and Atmospheric Administration–National Weather Service vary from 30 to 120 days, depending on datum conversion or datum-change case scenarios. Notification times associated with these case scenarios may be useful to the National Oceanic and Atmospheric Administration–National Weather Service and U.S. Army Corps of Engineers, because their daily operations are time sensitive, unlike the notification time change requirements of other entities that make up the cooperative water community. At the time of this writing, a future geopotential datum resulting from Gravity for the Redefinition of the American Vertical Datum is anticipated in 2022. A future version of VDatum and VERTCON is anticipated to provide a transformation among North American Vertical Datum of 1988 elevations to the new geopotential datum.

Techniques and Methods↗
USGS2011

Digital hydrologic networks supporting applications related to spatially referenced regression modeling

Digital hydrologic networks depicting surface-water pathways and their associated drainage catchments provide a key component to hydrologic analysis and modeling. Collectively, they form common spatial units that can be used to frame the descriptions of aquatic and watershed processes. In addition, they provide the ability to simulate and route the movement of water and associated constituents throughout the landscape. Digital hydrologic networks have evolved from derivatives of mapping products to detailed, interconnected, spatially referenced networks of water pathways, drainage areas, and stream and watershed characteristics. These properties are important because they enhance the ability to spatially evaluate factors that affect the sources and transport of water-quality constituents at various scales. SPAtially Referenced Regressions On Watershed attributes (SPARROW), a process-based &frasl; statistical model, relies on a digital hydrologic network in order to establish relations between quantities of monitored contaminant flux, contaminant sources, and the associated physical characteristics affecting contaminant transport. Digital hydrologic networks modified from the River Reach File (RF1) and National Hydrography Dataset (NHD) geospatial datasets provided frameworks for SPARROW in six regions of the conterminous United States. In addition, characteristics of the modified RF1 were used to update estimates of mean-annual streamflow. This produced more current flow estimates for use in SPARROW modeling.

Journal of the American Water Resources Associatio↗
USGS2018

Implementation of MOVE.1, censored MOVE.1, and piecewise MOVE.1 low-flow regressions with applications at partial-record streamgaging stations in New Jersey

The U.S. Geological Survey (USGS) uses Maintenance of Variance Extension Type 1 (MOVE.1) regression to transfer streamflows measured at long-term continuous-record streamgaging stations to partial-record (PR) streamgaging stations where intermittent base-flow measurements are available. MOVE.1 regression is used widely throughout the hydrologic community to extend historic low flows and low-flow statistics at continuous-record streamgaging stations to streamgaging stations that have access to only a partial record of low flows. The method correlates base-flow measurements at PR streamgaging stations with daily mean streamflows measured at index stations that exhibit similar streamflow characteristics. Following changes in the computing platform for storing, processing, retrieving, and publishing National Water Information System (NWIS) hydrologic data, legacy Statistical Analysis System (SAS) code developed by the USGS to implement the MOVE.1 regression was no longer suitable for reading and processing NWIS streamflow data. To migrate the MOVE.1 program so that it could continue to read streamflow data using the new hydrologic data platform, the SAS code was re-written in R, an open source programming language and software environment for statistical computing and graphics supported by the R Foundation for Statistical Computing. The work described in this report was performed in a study conducted by USGS in cooperation with the New Jersey Department of Environmental Protection. During migration from SAS to R, graphical and tabular output generated by the R script was compared to output produced by the legacy SAS code to ensure that equations used to perform the MOVE.1 regression remained the same. An option to perform censored MOVE.1 regression was added to extend the MOVE.1 methodology to cases where one or more measured continuous-record or PR streamgaging station flows are zero valued. In addition to permitting censored regression, the new R script includes an option to perform piecewise MOVE.1 regression when the relation between PR station and index station low flows varies significantly across the range of index station streamflows. Together with traditional MOVE.1 regression, censored, and piecewise MOVE.1 regression methods implemented by the R script offer less biased estimates than ordinary least squares regression for the annual 7-day 10-year and other low-flow statistics at PR stations for a range of base-flow conditions. The R script is used to implement the MOVE.1 regression methods across a variety of computing platforms.

New Jersey↗
USGS2020

Time scales of arsenic variability and the role of high-frequency monitoring at three water-supply wells in New Hampshire, USA

Groundwater geochemistry, redox process classification, high-frequency physicochemical and hydrologic measurements, and climate data were analyzed to identify controls on arsenic (As) concentration changes. Groundwater was monitored in two public-supply wells (one glacial aquifer and one bedrock aquifer), and one bedrock-aquifer domestic well in New Hampshire, USA, from 2014 to 2018 to identify time scales of and controls on As concentration changes. Concentrations of As and other geochemical constituents were measured bimonthly. Specific conductance (SC), pH, dissolved oxygen, and pumping rate/water level were measured at high frequency (every 5 to 15 min). Median (and 95% confidence interval) As concentrations at the three wells were 4.1 (3.7–4.6), 18.9 (17.2–23.6), and 37.5 (30.4–42.9) μg/L. Arsenic variability in each of the three wells, in relative standard deviation, ranged from 9 to 12%. Median quarterly As concentrations were highest in all wells in the spring. The bedrock-aquifer public-supply well As concentration increased over the period of study while pumping rate decreased. In the public-supply wells, As variability was correlated with SC and pH, and As species were related to SC, pH, pumping, precipitation, and changes in redox process. Specific conductance also had a seasonal pattern in the two public-supply wells and was correlated with Na and Cl. Excess Na in water samples suggests possible ion exchange with dissolved Ca, creating more capacity to dissolve CaCO 3 from calcareous rocks, which can increase pH and in turn, As concentrations in wells. High-frequency monitoring data are cost effective to collect, which could be advantageous in other parts of the United States and in the many parts of the world where glacial aquifers are in direct contact with other water supply aquifers or where water from different aquifers have potential to mix.

New Hampshire↗
USGS2021

Rapid-response unsaturated zone hydrology: Small-scale data, small-scale theory, big problems

The unsaturated zone (UZ) extends across the Earth’s terrestrial surface and is central to many problems related to land and water resource management. Flow of water through the UZ is typically thought to be slow and diffusive, such that it could attenuate fluxes and dampen variability between atmospheric inputs and underlying aquifer systems. This would reduce water resource vulnerability to contaminants and water-related hazards. Reducing or negating that effect, however, spatially concentrated and rapid flow and transport through the unsaturated zone is surprisingly common and becoming more so with the increasing frequency and magnitude of extreme hydroclimatic events. Arising from the wide range in the rates and complex modes of nonlinear flow processes, these effects are among the most poorly characterized hydrologic phenomena. Issues of scale present additional difficulties. Equations representing unsaturated processes have been developed and tested on the basis of field and laboratory measurements typically made at scales from pore size to plot size. In contrast, related problems of significant interest to society, including floods, aquifer recharge, landslides, and groundwater contamination, range from watershed to regional scales. The disparity between the scale of our understanding and the scale of interest for societal problems has spurred application of these model equations at increasingly coarse resolutions over larger areas than can be justified by existing measurements or theory. This mismatch in scales requires an assumption that spatially averaging slow diffusive flow and rapid preferential flow can effectively represent the influence of both processes across vast areas. Given the currently inadequate recognition and quantitative characterization of focused and rapid processes in unsaturated flow, these phenomena are critically in need of expanded attention and effort.

Frontiers in Earth Science↗
USGS2007

Nutrient vectors and riparian processing: A review with special reference to African semiarid savanna ecosystems

This review article describes vectors for nitrogen and phosphorus delivery to riparian zones in semiarid African savannas, the processing of nutrients in the riparian zone and the effect of disturbance on these processes. Semiarid savannas exhibit sharp seasonality, complex hillslope hydrology and high spatial heterogeneity, all of which ultimately impact nutrient fluxes between riparian, upland and aquatic environments. Our review shows that strong environmental drivers such as fire and herbivory enhance nitrogen, phosphorus and sediment transport to lower slope positions by shaping vegetative patterns. These vectors differ significantly from other arid and semiarid ecosystems, and from mesic ecosystems where the impact of fire and herbivory are less pronounced and less predictable. Also unique is the presence of sodic soils in certain hillslopes, which substantially alters hydrological flowpaths and may act as a trap where nitrogen is immobilized while sediment and phosphorus transport is enhanced. Nutrients and sediments are also deposited in the riparian zone during seasonal, intermittent floods while, during the dry season, subsurface movement of water from the stream into riparian soils and vegetation further enrich riparian zones with nutrients. As is found in mesic ecosystems, nutrients are immobilized in semiarid riparian corridors through microbial and plant uptake, whereas dissimilatory processes such as denitrification may be important where labile nitrogen and carbon are in adequate supply and physical conditions are suitable-such as in seeps, wallows created by animals, ephemeral wetlands and stream edges. Interaction between temporal hydrologic connectivity and spatial heterogeneity are disrupted by disturbances such as large floods and extended droughts, which may convert certain riparian patches from sinks to sources for nitrogen and phosphorus. In the face of increasing anthropogenic pressure, the scientific challenges are to provide a basic understanding of riparian biogeochemistry in semiarid African savannas to adequately address the temporal and spatial impact of disturbances, and to apply this knowledge to better regional land and water management. An integrated, multidisciplinary approach applied in protected as well as human-disturbed ecosystems in southern Africa is essential for underpinning a strong environmental basis for sustainable human-related expansion.

Ecosystems↗
USGS2008

Technical Manual for the Geospatial Stream Flow Model (GeoSFM)

The monitoring of wide-area hydrologic events requires the use of geospatial and time series data available in near-real time. These data sets must be manipulated into information products that speak to the location and magnitude of the event. Scientists at the U.S. Geological Survey Earth Resources Observation and Science (USGS EROS) Center have implemented a hydrologic modeling system which consists of an operational data processing system and the Geospatial Stream Flow Model (GeoSFM). The data processing system generates daily forcing evapotranspiration and precipitation data from various remotely sensed and ground-based data sources. To allow for rapid implementation in data scarce environments, widely available terrain, soil, and land cover data sets are used for model setup and initial parameter estimation. GeoSFM performs geospatial preprocessing and postprocessing tasks as well as hydrologic modeling tasks within an ArcView GIS environment. The integration of GIS routines and time series processing routines is achieved seamlessly through the use of dynamically linked libraries (DLLs) embedded within Avenue scripts. GeoSFM is run operationally to identify and map wide-area streamflow anomalies. Daily model results including daily streamflow and soil water maps are disseminated through Internet map servers, flood hazard bulletins and other media.

Open-File Report↗
USGS2024

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↗
USGS2004

Sensitivity to acidification of subalpine ponds and lakes in north-western Colorado

Although acidifying deposition in western North America is lower than in many parts of the world, many high-elevation ecosystems there are extremely sensitive to acidification. Previous studies determined that the Mount Zirkel Wilderness Area (MZWA) has the most acidic snowpack and aquatic ecosystems that are among the most sensitive in the region. In this study, spatial and temporal variability of ponds and lakes in and near the MZWA were examined to determine their sensitivity to acidification and the effects of acidic deposition during and after snowmelt. Within the areas identified as sensitive to acidification based on bedrock types, there was substantial variability in acid-neutralizing capacity (ANC), which was related to differences in hydrological flowpaths that control delivery of weathering products to surface waters. Geological and topographic maps were of limited use in predicting acid sensitivity because their spatial resolution was not fine enough to capture the variability of these attributes for lakes and ponds with small catchment areas. Many of the lakes are sensitive to acidification (summer and autumn ANC < 100 µeq L −1 ), but none of them appeared to be threatened immediately by episodic or chronic acidification. In contrast, 22 ponds had minimum ANC < 30 µeq L −1 , indicating that they are extremely sensitive to acidic deposition and could be damaged by episodic acidification, although net acidity (ANC < 0) was not measured in any of the ponds during the study. The lowest measured pH value was 5·4, and pH generally remained less than 6·0 throughout early summer in the most sensitive ponds, indicating that biological effects of acidification are possible at levels of atmospheric deposition that occurred during the study. The aquatic chemistry of lakes was dominated by atmospheric deposition and biogeochemical processes in soils and shallow ground water, whereas the aquatic chemistry of ponds was also affected by organic acids and biogeochemical processes in the water column and at the sediment–water interface. These results indicate that conceptual and mechanistic acidification models that have been developed for lakes and streams may be inadequate for predicting acidification in less-understood systems such as ponds.

Hydrological Processes↗
USGS2011

Direction of unsaturated flow in a homogeneous and isotropic hillslope

The distribution of soil moisture in a homogeneous and isotropic hillslope is a transient, variably saturated physical process controlled by rainfall characteristics, hillslope geometry, and the hydrological properties of the hillslope materials. The major driving mechanisms for moisture movement are gravity and gradients in matric potential. The latter is solely controlled by gradients of moisture content. In a homogeneous and isotropic saturated hillslope, absent a gradient in moisture content and under the driving force of gravity with a constant pressure boundary at the slope surface, flow is always in the lateral downslope direction, under either transient or steady state conditions. However, under variably saturated conditions, both gravity and moisture content gradients drive fluid motion, leading to complex flow patterns. In general, the flow field near the ground surface is variably saturated and transient, and the direction of flow could be laterally downslope, laterally upslope, or vertically downward. Previous work has suggested that prevailing rainfall conditions are sufficient to completely control these flow regimes. This work, however, shows that under time-varying rainfall conditions, vertical, downslope, and upslope lateral flow can concurrently occur at different depths and locations within the hillslope. More importantly, we show that the state of wetting or drying in a hillslope defines the temporal and spatial regimes of flow and when and where laterally downslope and/or laterally upslope flow occurs.

Water Resources Research↗