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.

SourceResearch connectionFrom yearThrough year

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

At least 1,513 records · Page 84Linked to original sources

USGS1998

Application of the top specified boundary layer (TSBL) approximation to initial characterization of an inland aquifer mineralization: 2. Seepage of saltwater through semi-confining layers

This paper presents a generalized basic study that addresses practical needs for an understanding of the major mechanisms involved in the mineralization of groundwater in the Great Bend Prairie aquifer in south- central Kansas. This Quaternary alluvial aquifer and associated surface waters are subject to contamination by saltwater, which in some areas seeps from the deeper Permian bedrock formation into the overlying freshwater aquifer through semiconfining layers. A simplified conceptual model is adopted. It incorporates the freshwater aquifer whose bottom is comprised of a semiconfining layer through which a hydrologically minor but geochemically important saline water discharge seeps into the aquifer. A hierarchy of approximate approaches is considered to analyze the mineralization processes taking place in the aquifer. The recently developed top specified boundary layer (TSBL) approach is very convenient to use for the initial characterization of these processes, and is further adapted to characterization of head-driven seepage through semi-confining layers. TSBL calculations indicate that the seeping saline water may create two distinct new zones in the aquifer: (1) a completely saline zone (CSZ) adjacent to the semiconfining bottom of the aquifer, and (2) a transition zone (TZ) which develops between the CSZ and the freshwater zone. Some possible scenarios associated with the various mineralization patterns are analyzed and discussed.

Journal of Contaminant Hydrology↗

USGS2018

Mercury on a landscape scale—Balancing regional export with wildlife health

The Cosumnes River watershed requires a 57–64 percent reduction in loads to meet the new Delta methylmercury (MeHg) total maximum daily load allocation, established by the Central Valley Regional Water Quality Control Board. Because there are no large point sources of MeHg in the watershed, the focus of MeHg load reductions will fall upon non-point sources, particularly the expansive wetlands considered to be a primary source of MeHg in the region. Few management practices have been implemented and tested in order to meet load reductions in managed wetlands, but recent efforts have shown promise. This project examines a treatment approach to reduce MeHg loads to the Sacramento-San Joaquin River Delta by creating open-water deep cells with a small footprint at the downstream end of wetlands to promote net demethylation of MeHg and to minimize MeHg and Hg loads exiting wetlands at the Cosumnes River Preserve. Specifically, the deep cells were were located immediately up gradient of the wetland’s outflow weir and were deep enough (75–91 centimeter depth) to be vegetation-free. The topographic and hydrologic structure of each treatment wetland was modified to include open-water deep cells so that the removal of aqueous MeHg might be enhanced through (1) particle settling, (2) photo-degradation, and (3) benthic microbial demethylation. These deep cells were, therefore, expected to clean MeHg from surface water prior to its discharge to the Cosumnes River and the downstream Delta. Our goal was to test whether the implementation of the deep cells within wetlands would minimize MeHg and total Hg export. Further, we sought to test whether continuous flow-through hydrology, would lower MeHg concentrations in resident biota, compared to traditional wetland management operations. The dominant practice in seasonal wetlands management is the “fill-and-maintain” approach, in which wetlands are filled with water and the water levels maintained without substantial draining until drawdown. Our approach was to create and characterize replicate treatment wetland complexes, in conjunction with monitoring of hydrologic, biologic, and chemical indicators of MeHg exposure for two full annual cycles within winter-spring flooded seasonal wetlands. In addition to the creation of deep cells within treatment wetlands, hydrology was manipulated so that there was a constant flow-through of water, while the control wetlands utilized the fill-and-maintain approach. Specifically, the treatment wetlands were maintained in a flow-through manner, while the control wetlands were maintained in a fill-and-maintain manner from September through May, to test the hypothesis that the flow of water through the seasonal wetland can lower fish bioaccumulation through dilution of MeHg-concentrated water within the wetland by constant inflows of water into the wetland. The major tasks of this study included: (1) field design and implementation, (2) water and wetland management, (3) hydrologic monitoring and water quality sampling, (4) MeHg export and load estimates, (5) caged fish experiments for examining MeHg bioaccumulation, (6) site and process characterization to improve understanding and transferability of results, (7) adaptive management, transferability, and outreach, and (8) reporting of results and conclusions. This report summarizes the key findings of this study, which focuses on MeHg load estimates from control and treatment wetlands, quantification of three MeHg removal mechanisms (particulate settling, benthic demethylation, and photo-demethylation) in the deep cells within the treatment wetlands, and MeHg bioaccumulation in wetland fishes. Key findings include: Over two years of study, mean whole-water MeHg load decreased 37 percent in deep cells, when comparing inlet of check weir flows to outlet. Of the 37 percent MeHg load removed within the deep cell, photodegradation accounted for 7 percent and particle flux to the benthos accounted for 24 percent of the mass removed, with the remaining 6 percent apparent MeHg loss unexplained. Benthic MeHg degradation did not appear to be a major MeHg removal process in the deep cells, as changes in the ambient MeHg pool over 7-day bottle incubations showed that the surface sediment exhibited net MeHg production in the majority (87 percent) of incubation experiments. In only 13 percent of the incubations (3 out of 24) was net MeHg degradation observed. Estimates of benthic diffusive flux of MeHg across the sediment/water interface were small relative to particulate flux and variable (positive or negative), suggesting this is likely a minor term in the overall MeHg budget within the deep cells. Although the deep cells served as net MeHg sink overall, MeHg export from the flow-through treatment wetlands (shallow and deep combined) exceeded export from the fill-and-maintain managed control wetlands, because of the differences in hydrologic management between the two wetland types. Shallow wetlands under flow-through conditions generated a net export of MeHg. Most of the annual MeHg export from the treatment wetlands occurred within the first 3 months of flood up (September to November), shortly after hydrologic management began. Despite the effectiveness of the deep cell in lowering MeHg export concentrations, total mercury (THg) concentration did not decrease in biosentinel fish ( Gambusia affinis , Mosquitofish) between the deep cell inlet and outlet. Mosquitofish THg concentrations were higher in treatment wetlands than in control wetlands during the first year of study, likely because of an associated increase in MeHg availability immediately following wetland construction activities. Mosquitofish THg concentrations declined in the treatment wetlands during the second year of study, and fish THg concentrations in treatment wetlands were no different from those in the control. Similarly, the increased hydrologic flow rates in the treatment wetlands did not lower fish THg concentrations nor aqueous MeHg concentrations in the shallow cells, suggesting that MeHg flux from the sediment to water column exceeded the flow-through flushing rate in the shallow portion of the treatment wetlands. Reductions in MeHg concentrations of surface water and fish may require higher flow rates than used in the study to achieve the region’s regulatory goals. However, the flow rates necessary may not be feasible for these managed wetlands because of limited water supply and the associated costs for water and pumping. The use of deep cells in seasonal wetlands were effective in lowering MeHg exports under continuous water flow-through hydrology. However, fill-and-maintain hydrology had lower exports overall, because of a single major drainage event at the end of the flood season. Future studies focused on limiting MeHg export should consider combining deep cells with the fill-and-maintain or fill-and-trickle hydrologic management approach.

California↗

USGS2015

Groundwater recharge assessment at local and episodic scale in a soil mantled perched karst aquifer in southern Italy

Groundwater recharge assessment of karst aquifers, at various spatial and temporal scales, is a major scientific topic of current importance, since these aquifers play an essential role for both socio-economic development and fluvial ecosystems. In this study, groundwater recharge was estimated at local and episodic scales in a representative perched karst aquifer in a region of southern Italy with a Mediterranean climate. The research utilized measurements of precipitation, air temperature, soil water content, and water-table depth, obtained in 2008 at the Acqua della Madonna test area (Terminio Mount karst aquifer, Campania region). At this location the aquifer is overlain by ash-fall pyroclastic soils. The Episodic Master Recession (EMR) method, an improved version of the Water Table Fluctuation (WTF) method, was applied to estimate the amount of recharge generated episodically by individual rainfall events. The method also quantifies the amount of precipitation generating each recharge episode, thus permitting calculation of the Recharge to the Precipitation Ratio (RPR) on a storm-by-storm basis. Depending on the seasonally varying air temperature, evapotranspiration, and precipitation patterns, calculated values of RPR varied between 35% and 97% among the individual episodes. A multiple linear correlation of the RPR with both the average intensity of recharging rainfall events and the antecedent soil water content was calculated. Given the relatively easy measurability of precipitation and soil water content, such an empirical model would have great hydrogeological and practical utility. It would facilitate short-term forecasting of recharge in karst aquifers of the Mediterranean region and other aquifers with similar hydrogeological characteristics. By establishing relationships between the RPR and climate-dependent variables such as average storm intensity, it would facilitate prediction of climate-change effects on groundwater recharge. The EMR methodology could further be applied to other aquifers for evaluating the relationship of recharge to various hydrometeorological and hydrogeological processes.

Journal of Hydrology↗

USGS2010

An optimum reduction of gauges to meet data program constraints

Budget or manpower constraints may force a reduction in data collection activities. However, information may be transferred from continued gauge sites to discontinued gauge sites provided there is ‘sufficient’ correlation between flow sequences. Information defined as the reciprocal of variance (of the parameter estimator for which the gauge has been established) is provided at ungauged sites by regression techniques. A set of gauge sites is chosen to satisfy the budget constraint and to maximize the information content. The information content based on regression is a function of the correlations between flow sequences and is subject to the errors in their estimators. Post optimal analysis gives estimates of correlation error effects on the decision process as to which gauges to continue and which to discontinue.

Hydrological Sciences Bulletin↗

USGS2015

A 2-D process-based model for suspended sediment dynamics: A first step towards ecological modeling

In estuaries suspended sediment concentration (SSC) is one of the most important contributors to turbidity, which influences habitat conditions and ecological functions of the system. Sediment dynamics differs depending on sediment supply and hydrodynamic forcing conditions that vary over space and over time. A robust sediment transport model is a first step in developing a chain of models enabling simulations of contaminants, phytoplankton and habitat conditions. This works aims to determine turbidity levels in the complex-geometry delta of the San Francisco estuary using a process-based approach (Delft3D Flexible Mesh software). Our approach includes a detailed calibration against measured SSC levels, a sensitivity analysis on model parameters and the determination of a yearly sediment budget as well as an assessment of model results in terms of turbidity levels for a single year, water year (WY) 2011. Model results show that our process-based approach is a valuable tool in assessing sediment dynamics and their related ecological parameters over a range of spatial and temporal scales. The model may act as the base model for a chain of ecological models assessing the impact of climate change and management scenarios. Here we present a modeling approach that, with limited data, produces reliable predictions and can be useful for estuaries without a large amount of processes data.

Hydrology and Earth System Sciences↗

USGS2016

Groundwater and surface-water interaction, water quality, and processes affecting loads of dissolved solids, selenium, and uranium in Fountain Creek, near Pueblo, Colorado, 2012–2014

In 2012, the U.S. Geological Survey, in cooperation with the Arkansas River Basin Regional Resource Planning Group, initiated a study of groundwater and surface-water interaction, water quality, and loading of dissolved solids, selenium, and uranium to Fountain Creek near Pueblo, Colorado, to improve understanding of sources and processes affecting loading of these constituents to streams in the Arkansas River Basin. Fourteen monitoring wells were installed in a series of three transects across Fountain Creek near Pueblo, and temporary streamgages were established at each transect to facilitate data collection for the study. Groundwater and surface-water interaction was characterized by using hydrogeologic mapping, groundwater and stream-surface levels, groundwater and stream temperatures, vertical hydraulic-head gradients and ratios of oxygen and hydrogen isotopes in the hyporheic zone, and streamflow mass-balance measurements. Water quality was characterized by collecting periodic samples from groundwater, surface water, and the hyporheic zone for analysis of dissolved solids, selenium, uranium, and other selected constituents and by evaluating the oxidation-reduction condition for each groundwater sample under different hydrologic conditions throughout the study period. Groundwater loads to Fountain Creek and in-stream loads were computed for the study area, and processes affecting loads of dissolved solids, selenium, and uranium were evaluated on the basis of geology, geochemical conditions, land and water use, and evapoconcentration. During the study period, the groundwater-flow system generally contributed flow to Fountain Creek and its hyporheic zone (as a single system) except for the reach between the north and middle transects. However, the direction of flow between the stream, the hyporheic zone, and the near-stream aquifer was variable in response to streamflow and stage. During periods of low streamflow, Fountain Creek generally gained flow from groundwater. However, during periods of high streamflow, the hydraulic gradient between groundwater and the stream temporarily reversed, causing the stream to lose flow to groundwater. Concentrations of dissolved solids, selenium, and uranium in groundwater generally had greater spatial variability than surface water or hyporheic-zone samples, and constituent concentrations in groundwater generally were greater than in surface water. Constituent concentrations in the hyporheic zone typically were similar to or intermediate between concentrations in groundwater and surface water. Concentrations of dissolved solids, selenium, uranium, and other constituents in groundwater samples collected from wells located on the east side of the north monitoring well transect were substantially greater than for other groundwater, surface-water, and hyporheic-zone samples. With one exception, groundwater samples collected from wells on the east side of the north transect exhibited oxic to mixed (oxic-anoxic) conditions, whereas most other groundwater samples exhibited anoxic to suboxic conditions. Concentrations of dissolved solids, selenium, and uranium in surface water generally increased in a downstream direction along Fountain Creek from the north transect to the south transect and exhibited an inverse relation to streamflow with highest concentration occurring during periods of low streamflow and lowest concentrations occurring during periods of high streamflow. Groundwater loads of dissolved solids, selenium, and uranium to Fountain Creek were small because of the small amount of groundwater flowing to the stream under typical low-streamflow conditions. In-stream loads of dissolved solids, selenium, and uranium in Fountain Creek varied by date, primarily in relation to streamflow at each transect and were much larger than computed constituent loads from groundwater. In-stream loads generally decreased with decreases in streamflow and increased as streamflow increased. In-stream loads of dissolved solids and selenium increased between the north and middle transects but generally decreased between the middle and south transects. By contrast, uranium loads generally decreased between the north and middle transects but increased between the middle and south transects. In-stream load differences between transects appear primarily to be related to differences in streamflow. However, because groundwater typically flows to Fountain Creek under low-flow conditions, and groundwater has greater concentrations of dissolved solids, selenium, and uranium than surface water in Fountain Creek, increases in loads between transects likely are affected by inflow of groundwater to the stream, which can account for a substantial proportion of the in-stream load difference between transects. When loads decreased between transects, the primary cause likely was decreased streamflow as a result of losses to groundwater and flow through the hyporheic zone. However, localized groundwater inflow likely attenuated the magnitude by which the in-stream loads decreased. The combination of localized soluble geologic sources and oxic conditions likely is the primary reason for the occurrence of high concentrations of dissolved solids, selenium, and uranium in groundwater on the east side of the north monitoring well transect. To evaluate conditions potentially responsible for differences in water quality and redox conditions, physical characteristics such as depth to water, saturated thickness, screen depth below the water table, screen height above bedrock, and aquifer hydraulic conductivity were compared by using Wilcoxon rank-sum tests. Results indicated no significant difference between depth to water, screen height above bedrock, and hydraulic conductivity for groundwater samples collected from wells on the east side of the north transect and groundwater samples from all other wells. However, saturated thickness and screen depth below the water table both were significantly smaller for groundwater samples collected from wells on the east side of the north transect than for groundwater samples from other wells, indicating that these characteristics might be related to the elevated constituent concentrations found at that location. Similarly, saturated thickness and screen depth below the water table were significantly smaller for groundwater samples under oxic or mixed (oxic-anoxic) conditions than for those under anoxic to suboxic conditions. The greater constituent concentrations at wells on the east side of the north transect also could, in part, be related to groundwater discharge from an unnamed alluvial drainage located directly upgradient from that location. Although the quantity and quality of water discharging from the drainage is not known, the drainage appears to collect water from a residential area located upgradient to the east of the wells, and groundwater could become concentrated in nitrate and other dissolved constituents before flowing through the drainage. High levels of nitrate, whether from anthropogenic or natural geologic sources, could promote more soluble forms of selenium and other constituents by affecting the redox condition of groundwater. Whether oxic conditions at wells on the east side of the north transect are the result of physical characteristics or of groundwater inflow from the alluvial drainage, the oxic conditions appear to cause increased dissolution of minerals from the shallow shale bedrock at that location. Because ratios of hydrogen and oxygen isotopes indicate evaporation likely has not had a substantial effect on groundwater, constituent concentrations at that location likely are not the result of evapoconcentration.

Colorado↗

USGS2008

Use of an integrated flow model to estimate ecologically relevant hydrologic characteristics at stream biomonitoring sites

We developed an integrated hydroecological model to provide a comprehensive set of hydrologic variables representing five major components of the flow regime at 856 aquatic-invertebrate monitoring sites in New Jersey. The hydroecological model simulates streamflow by routing water that moves overland and through the subsurface from atmospheric delivery to the watershed outlet. Snow accumulation and melt, evapotranspiration, precipitation, withdrawals, discharges, pervious- and impervious-area runoff, and lake storage were accounted for in the water balance. We generated more than 78 flow variables, which describe the frequency, magnitude, duration, rate of change, and timing of flow events. Highly correlated variables were filtered by principal component analysis to obtain a non-redundant subset of variables that explain the majority of the variation in the complete set. This subset of variables was used to evaluate the effect of changes in the flow regime on aquatic-invertebrate assemblage structure at 856 biomonitoring sites. We used non-metric multidimensional scaling (NMS) to evaluate variation in aquatic-invertebrate assemblage structure across a disturbance gradient. We employed multiple linear regression (MLR) analysis to build a series of MLR models that identify the most important environmental and hydrologic variables driving the differences in the aquatic-invertebrate assemblages across the disturbance gradient. The first axis of NMS ordination was significantly related to many hydrologic, habitat, and land-use/land-cover variables, including the average number of annual storms producing runoff, ratio of 25-75% exceedance flow (flashiness), diversity of natural stream substrate, and the percentage of forested land near the stream channel (forest buffer). Modifications in the hydrologic regime as the result of changes in watershed land use appear to promote the retention of highly tolerant aquatic species; in contrast, species that are sensitive to hydrologic instability and other anthropogenic disturbance become much less prevalent. We also found strong relations between an index of invertebrate-assemblage impairment, its component metrics, and the primary disturbance gradient. The process-oriented watershed modeling approach used in this study provides a means to evaluate how natural landscape features interact with anthropogenic factors and assess their effects on flow characteristics and stream ecology. By combining watershed modeling and indirect ordination techniques, we were able to identify components of the hydrologic regime that have a considerable effect on aquatic-assemblage structure and help in developing short- and long-term management measures that mitigate the effects of anthropogenic disturbance in stream systems.

New Jersey↗

USGS2014

Control on groundwater flow in a semiarid folded and faulted intermountain basin

The major processes controlling groundwater flow in intermountain basins are poorly understood, particularly in basins underlain by folded and faulted bedrock and under regionally realistic hydrogeologic heterogeneity. To explore the role of hydrogeologic heterogeneity and poorly constrained mountain hydrologic conditions on regional groundwater flow in contracted intermountain basins, a series of 3-D numerical groundwater flow models were developed using the South Park basin, Colorado, USA as a proxy. The models were used to identify the relative importance of different recharge processes to major aquifers, to estimate typical groundwater circulation depths, and to explore hydrogeologic communication between mountain and valley hydrogeologic landscapes. Modeling results show that mountain landscapes develop topographically controlled and predominantly local-scale to intermediate-scale flow systems. Permeability heterogeneity of the fold and fault belt and decreased topographic roughness led to permeability controlled flow systems in the valley. The structural position of major aquifers in the valley fold and fault belt was found to control the relative importance of different recharge mechanisms. Alternative mountain recharge model scenarios showed that higher mountain recharge rates led to higher mountain water table elevations and increasingly prominent local flow systems, primarily resulting in increased seepage within the mountain landscape and nonlinear increases in mountain block recharge to the valley. Valley aquifers were found to be relatively insensitive to changing mountain water tables, particularly in structurally isolated aquifers inside the fold and fault belt.

Colorado↗

USGS2016

Completion summary for boreholes TAN-2271 and TAN‑2272 at Test Area North, Idaho National Laboratory, Idaho

In 2015, the U.S. Geological Survey, in cooperation with the U.S. Department of Energy, drilled and constructed boreholes TAN-2271 and TAN-2272 for stratigraphic framework analyses and long-term groundwater monitoring of the eastern Snake River Plain aquifer at the Idaho National Laboratory in southeast Idaho. Borehole TAN-2271 initially was cored to collect continuous geologic data, and then re-drilled to complete construction as a monitor well. Borehole TAN-2272 was partially cored between 210 and 282 feet (ft) below land surface (BLS) then drilled and constructed as a monitor well. Boreholes TAN-2271 and TAN-2272 are separated by about 63 ft and have similar geologic layers and hydrologic characteristics based on geologic, geophysical, and aquifer test data collected. The final construction for boreholes TAN-2271 and TAN-2272 required 10-inch (in.) diameter carbon-steel well casing and 9.9-in. diameter open-hole completion below the casing to total depths of 282 and 287 ft BLS, respectively. Depth to water is measured near 228 ft BLS in both boreholes. Following construction and data collection, temporary submersible pumps and water-level access lines were placed to allow for aquifer testing, for collecting periodic water samples, and for measuring water levels. Borehole TAN-2271 was cored continuously, starting at the first basalt contact (about 33 ft BLS) to a depth of 284 ft BLS. Excluding surface sediment, recovery of basalt and sediment core at borehole TAN-2271 was better than 98 percent. Based on visual inspection of core and geophysical data, material examined from 33 to 211ft BLS primarily consists of two massive basalt flows that are about 78 and 50 ft in thickness and three sediment layers near 122, 197, and 201 ft BLS. Between 211 and 284 ft BLS, geophysical data and core material suggest a high occurrence of fractured and vesicular basalt. For the section of aquifer tested, there are two primary fractured aquifer intervals: the first between 235 and 255 ft BLS and the second between 272 and 282 ft BLS. Basalt texture for borehole TAN-2271 generally was described as aphanitic, phaneritic, and porphyritic. Sediment layers, starting near 122 ft BLS, generally were composed of fine-grained sand and silt with a lesser amount of clay. Basalt flows generally ranged in thickness from 2 to 78 ft and varied from highly fractured to dense with high to low vesiculation. Geophysical data and limited core material collected from TAN-2272 show similar lithologic sequences to those reported for TAN-2271. Geophysical and borehole video logs were collected during certain stages of the drilling and construction process at boreholes TAN-2271 and TAN-2272. Geophysical logs were examined synergistically with available core material to confirm geologic and hydrologic similarities and suggest possible fractured network interconnection between boreholes TAN-2271 and TAN-2272. Natural gamma log measurements were used to assess the completeness of the vapor port lines behind 10-in. diameter well casing. Electromagnetic flow meter results were used to identify downward flow conditions that exist for boreholes TAN-2271 and TAN-2272. Furthermore, gyroscopic deviation measurements were used to measure horizontal and vertical displacement at all depths in boreholes TAN-2271 and TAN-2272. After borehole construction was completed, single‑well aquifer tests were done within wells TAN-2271 and TAN ‑ 2272 to provide estimates of transmissivity and hydraulic conductivity. The transmissivity and hydraulic conductivity were estimated for the pumping well and observation well during the aquifer tests conducted on August 25 and August 27, 2015. Estimates for transmissivity range from 4.1 . 10 3 feet squared per day (ft 2 /d) to 8.1 . 10 3 ft 2 /d; estimates for hydraulic conductivity range from 5.8 to 11.5 feet per day (ft/d). Both TAN-2271 and TAN ‑ 2272 show sustained pumping rates of about 30 gallons per minute (gal/min) with measured drawdown in the pumping well of 1.96 ft and 1.14 ft, respectively. The transmissivity estimates for wells tested were within the range of values determined from previous aquifer tests in other wells near Test Area North. Groundwater samples were collected from both wells and were analyzed for cations, anions, metals, nutrients, volatile organic compounds, stable isotopes, and radionuclides. Groundwater samples for most of the inorganic constituents showed similar water chemistry in both wells. Groundwater samples for strontium-90, trichloroethene, and vinyl chloride exceeded maximum contaminant levels for public drinking water supplies in one or both wells.

Idaho↗

USGS1977

Hydrologic characteristics of the Madison Limestone, the Minnelusa Formation, and equivalent rocks as determined by well-logging formation evaluation, Wyoming, Montana, South Dakota, and North Dakota

Geophysical logs from wells distributed throughout the Powder River Basin were digitized, processed, and interpreted to get a regional understanding of the lithologic and ground-water characteristics of aquifers in the Madison Limestone and Minnelusa Formation. The percentage of sand, porosity, and apparent ground-water resistivity of the Minnelusa closely follow structural trends in the basin. Values increase from the structural axis of the basin toward the flanks. Close correlations between the water resistivities in the Madison and the Minnelusa indicate possible hydrologic connection between the two units. The low primary porosity from geophysical logs in the Madison indicates that secondary porosity is necessary to develop water supplies of economic importance from the Madison. Because of increased porosity and better water quality, the flanks of the basin appear to be the most favorable places for water production from either the Madison or the Minnelusa.

Wyoming, Montana, South Dakota, North Dakota↗

USGS2013

Extending airborne electromagnetic surveys for regional active layer and permafrost mapping with remote sensing and ancillary data, Yukon Flats ecoregion, central Alaska

Machine-learning regression tree models were used to extrapolate airborne electromagnetic resistivity data collected along flight lines in the Yukon Flats Ecoregion, central Alaska, for regional mapping of permafrost. This method of extrapolation (r = 0.86) used subsurface resistivity, Landsat Thematic Mapper (TM) at-sensor reflectance, thermal, TM-derived spectral indices, digital elevation models and other relevant spatial data to estimate near-surface (0–2.6-m depth) resistivity at 30-m resolution. A piecewise regression model (r = 0.82) and a presence/absence decision tree classification (accuracy of 87%) were used to estimate active-layer thickness (ALT) (< 101 cm) and the probability of near-surface (up to 123-cm depth) permafrost occurrence from field data, modelled near-surface (0–2.6 m) resistivity, and other relevant remote sensing and map data. At site scale, the predicted ALTs were similar to those previously observed for different vegetation types. At the landscape scale, the predicted ALTs tended to be thinner on higher-elevation loess deposits than on low-lying alluvial and sand sheet deposits of the Yukon Flats. The ALT and permafrost maps provide a baseline for future permafrost monitoring, serve as inputs for modelling hydrological and carbon cycles at local to regional scales, and offer insight into the ALT response to fire and thaw processes.

Alaska↗

USGS2021

Upland burning and grazing as strategies to offset climate-change effects on wetlands

Wetland ecosystems perform a multitude of services valued by society and provide critical habitat for migratory birds and other wildlife. Despite their importance, wetlands have been lost to different local, regional, and global drivers. Remaining wetlands are extremely sensitive to changing temperature and precipitation regimes. Management of grassland areas in wetland catchments may be an effective strategy for counteracting potentially negative impacts of climate change on wetlands. Our objective was to estimate the effects of climate changes on wetland hydrology, and to explore strategies for increasing surface-water inputs to wetlands. We coupled a field study with process-based simulation modeling of wetland-water levels. We found that climate change could decrease the number of wetlands that hold ponded water during the waterfowl breeding season by 14% under a hot wet scenario or 29% under a hot dry scenario if no upland-management actions were taken. Upland burning reduced pond losses to 9% (hot wet) and 26% (hot dry). Upland grazing resulted in the smallest loss of ponded wetlands, 6% loss under the hot-and-wet scenario and 22% loss under the hot-and-dry scenario. Overall, water inputs could be increased by either burning or grazing of upland vegetation thereby reducing pond losses during the waterfowl breeding season. While field results suggest that both grazing and burning can reduce the vegetative structure that could lead to increases in runoff in grassland catchments, our model simulations indicated that additional actions may be needed for managers to minimize future meteorologically driven water losses.

Wetlands Ecology and Management↗

USGS2020

Linking mesoscale meteorology with extreme landscape response: Effects of narrow cold frontal rainbands (NCFR)

Landscapes evolve in response to prolonged and/or intense precipitation resulting from atmospheric processes at various spatial and temporal scales. Whereas synoptic (large‐scale) features (e.g., atmospheric rivers and hurricanes) govern regional‐scale hydrologic hazards such as widespread flooding, mesoscale features such as thunderstorms or squall lines are more likely to trigger localized geomorphic hazards such as landslides. Thus, to better understand relations between hydrometeorological drivers and landscape response, a knowledge of mesoscale meteorology and its impacts is needed. Here we investigate the extreme geomorphic response associated with one type of mesoscale meteorological feature, the narrow cold frontal rainband (NCFR). Resulting from low‐level convergence and shallow convection along a cold front, NCFRs are narrow bands of high‐intensity rainfall that occur in midlatitude areas of the world. Our study examines an NCFR impacting the Sierra Nevada foothills (California, USA) that initiated over 500 landslides, mobilized ~360,000 metric tons of sediment to the fluvial system (as much as 16 times the local annual sediment yield), and severely damaged local infrastructure and regional water transport facilities. Coupling geomorphological field investigations with meteorological analyses, we demonstrate that precipitation associated with the NCFR was both intense (maximum 15 min intensity of 70 mm/hr) and localized, resulting in a highly concentrated band of shallow landsliding. This meteorological phenomenon likely plays an important role in landscape evolution and hazard initiation. Other types of mesoscale meteorological features also occur globally and offer new avenues for understanding the effects of storms on landscapes.

California↗

USGS2024

The effect of drying boreal lakes on plants, soils, and microbial communities in lake margin habitats

Decadal scale lake drying in interior Alaska results in lake margin colonization by willow shrub and graminoid vegetation, but the effects of these changes on plant production, biodiversity, soil properties, and soil microbial communities are not well known. We studied changes in soil organic carbon (SOC) and nitrogen (N) storage, plant and microbial community composition, and soil microbial activities in drying and non-drying lakes in the Yukon Flats National Wildlife Refuge. Historic changes in lake area were determined using Landsat imagery. Results showed that SOC storage in drying lake margins declined by 0.13 kg C m −2 yr −1 over 30 years of exposure of lake sediments, with no significant change in soil N. Lake drying resulted in an increase in graminoid and shrub aboveground net primary production (ANPP, +3% yr −1 ) with little change in plant functional composition. Increases in ANPP were similar in magnitude (but opposite in sign) to losses in SOC over a 30-year drying trend. Potential decomposition rates and soil enzyme activities were lower in drying lake margins compared to stable lake margins, possibly due to high salinities in drying lake margin soils. Microbial communities shifted in response to changing plant communities, although they still retained a legacy of the previous plant community. Understanding how changing lake hydrology impacts the ecology and biogeochemistry of lake margin terrestrial ecosystems is an underexamined phenomenon with large impacts to landscape processes.

Alaska↗

USGS2011

Simulating daily water temperatures of the Klamath River under dam removal and climate change scenarios

A one-dimensional daily averaged water temperature model was used to simulate Klamath River temperatures for two management alternatives under historical climate conditions and six future climate scenarios. The analysis was conducted for the Secretarial Determination on removal of four hydroelectric dams on the Klamath River. In 2012, the Secretary of the Interior will determine if dam removal and implementation of the Klamath Basin Restoration Agreement (KBRA) (Klamath Basin Restoration Agreement, 2010) will advance restoration of salmonid fisheries and is in the public interest. If the Secretary decides dam removal is appropriate, then the four dams are scheduled for removal in 2020. Water temperature simulations were conducted to compare the effect of two management alternatives: the no-action alternative where dams remain in place, and the action alternative where dam removal occurs in 2020 along with habitat restoration. Each management alternative was simulated under historical climate conditions (1961-2010) and six 50-year (2012-2061) climate scenarios. The model selected for the study, River Basin Model-10 (RBM10), was used to simulate water temperatures over a 253-mile reach of the Klamath River located in south-central Oregon and northern California. RBM10 uses a simple equilibrium flow model, assuming discharge in each river segment on each day is transmitted downstream instantaneously. The model uses a heat budget formulation to quantify heat flux at the air-water interface. Inputs for the heat budget were calculated from daily-mean meteorological data, including net shortwave solar radiation, net longwave atmospheric radiation, air temperature, wind speed, vapor pressure, and a psychrometric constant needed to calculate the Bowen ratio. The modeling domain was divided into nine reaches ranging in length from 10.8 to 42.4 miles, which were calibrated and validated separately with measured water temperature data collected irregularly from 1961 to 2010. Calibration root mean square errors of observed versus simulated water temperatures for the nine reaches ranged from 0.8 to 1.5 degrees C. Mean absolute errors ranged from 0.6 to 1.2 degrees C. For model validation, a k-fold cross-validation technique was used. Validation root mean square error and mean absolute error for the nine reaches ranged from 0.8 to 1.4 degrees C and 0.8 to 1.2 degrees C, respectively. Input data for the six future climate scenarios (2012-2061) were derived from historical hydrological and meteorological data and simulated meteorological output from five Global Circulation Models. Total Maximum Daily Loads or other regulatory processes that might reduce future water temperatures were not included in the simulations. Under the current climate conditions scenario, impacts of dam removal on water temperatures were greatest near Iron Gate Dam (near Yreka, California) and were attenuated in the lower reaches of the Klamath River. May and October simulated mean water temperatures increased and decreased by approximately 1-2 degrees C and 2-4 degrees C, respectively, downstream of Iron Gate Dam after dam removal. Dam removal also resulted in an earlier annual temperature cycle shift of 18 days, 5 days, and 2 days, near Iron Gate Dam, Scott River, and Trinity River, respectively. Although the magnitude of precipitation and air temperature change predicted by the five Global Circulation Models varied, all five models resulted in progressive incremental increases in water temperatures with each decade from 2012 to 2061. However, dam removal under KBRA appeared to delay the effects of climate change to some extent near Iron Gate Dam. With dam removal under KBRA, annual-mean water temperatures exceeded the 49-year historical mean temperature beginning in 2045; whereas with dams, annual-mean temperatures exceeded the historical mean beginning in 2025. Potential changes in seasonal water temperatures resulting from dam removal, with or without future climat

Open-File Report↗

USGS2006

Methods Used to Assess the Susceptibility to Contamination of Transient, Non-Community Public Ground-Water Supplies in Indiana

The Safe Water Drinking Act of 1974 as amended in 1996 gave each State the responsibility of developing a Source-Water Assessment Plan (SWAP) that is designed to protect public-water supplies from contamination. Each SWAP must include three elements: (1) a delineation of the source-water protection area, (2) an inventory of potential sources of contaminants within the area, and (3) a determination of the susceptibility of the public-water supply to contamination from the inventoried sources. The Indiana Department of Environmental Management (IDEM) was responsible for preparing a SWAP for all public-water supplies in Indiana, including about 2,400 small public ground-water supplies that are designated transient, non-community (TNC) supplies. In cooperation with IDEM, the U.S. Geological Survey compiled information on conditions near the TNC supplies and helped IDEM complete source-water assessments for each TNC supply. The delineation of a source-water protection area (called the assessment area) for each TNC ground-water supply was defined by IDEM as a circular area enclosed by a 300-foot radius centered at the TNC supply well. Contaminants of concern (COCs) were defined by IDEM as any of the 90 contaminants for which the U.S. Environmental Protection Agency has established primary drinking-water standards. Two of these, nitrate as nitrogen and total coliform bacteria, are Indiana State-regulated contaminants for TNC water supplies. IDEM representatives identified potential point and nonpoint sources of COCs within the assessment area, and computer database retrievals were used to identify potential point sources of COCs in the area outside the assessment area. Two types of methods-subjective and subjective hybrid-were used in the SWAP to determine susceptibility to contamination. Subjective methods involve decisions based upon professional judgment, prior experience, and (or) the application of a fundamental understanding of processes without the collection and analysis of data for a specific condition. Subjective hybrid methods combine subjective methods with quantitative hydrologic analyses. The subjective methods included an inventory of potential sources and associated contaminants, and a qualitative description of the inherent susceptibility of the area around the TNC supply. The description relies on a classification of the hydrogeologic and geomorphic characteristics of the general area around the TNC supply in terms of its surficial geology, regional aquifer system, the occurrence of fine- and coarse-grained geologic materials above the screen of the TNC well, and the potential for infiltration of contaminants. The subjective hybrid method combined the results of a logistic regression analysis with a subjective analysis of susceptibility and a subjective set of definitions that classify the thickness of fine-grained geologic materials above the screen of a TNC well in terms of impedance to vertical flow. The logistic regression determined the probability of elevated concentrations of nitrate as nitrogen (greater than or equal to 3 milligrams per liter) in ground water associated with specific thicknesses of fine-grained geologic materials above the screen of a TNC well. In this report, fine-grained geologic materials are referred to as a geologic barrier that generally impedes vertical flow through an aquifer. A geologic barrier was defined to be thin for fine-grained materials between 0 and 45 feet thick, moderate for materials between 45 and 75 feet thick, and thick if the fine-grained materials were greater than 75 feet thick. A flow chart was used to determine the susceptibility rating for each TNC supply. The flow chart indicated a susceptibility rating using (1) concentrations of nitrate as nitrogen and total coliform bacteria reported from routine compliance monitoring of the TNC supply, (2) the presence or absence of potential sources of regulated contaminants (nitrate as nitrogen and coliform bac

Indiana↗

USGS2016

Simulating groundwater flow in karst aquifers with distributed parameter models—Comparison of porous-equivalent media and hybrid flow approaches

Understanding karst aquifers, for purposes of their management and protection, poses unique challenges. Karst aquifers are characterized by groundwater flow through conduits (tertiary porosity), and (or) layers with interconnected pores (secondary porosity) and through intergranular porosity (primary or matrix porosity). Since the late 1960s, advances have been made in the development of numerical computer codes and the use of mathematical model applications towards the understanding of dual (primary [matrix] and secondary [fractures and conduits]) porosity groundwater flow processes, as well as characterization and management of karst aquifers. The Floridan aquifer system (FAS) in Florida and parts of Alabama, Georgia, and South Carolina is composed of a thick sequence of predominantly carbonate rocks. Karst features are present over much of its area, especially in Florida where more than 30 first-magnitude springs occur, numerous sinkholes and submerged conduits have been mapped, and numerous circular lakes within sinkhole depressions are present. Different types of mathematical models have been applied for simulation of the FAS. Most of these models are distributed parameter models based on the assumption that, like a sponge, water flows through connected pores within the aquifer system and can be simulated with the same mathematical methods applied to flow through sand and gravel aquifers; these models are usually referred to as porous-equivalent media models. The partial differential equation solved for groundwater flow is the potential flow equation of fluid mechanics, which is used when flow is dominated by potential energy and has been applied for many fluid problems in which kinetic energy terms are dropped from the differential equation solved. In many groundwater model codes (basic MODFLOW), it is assumed that the water has a constant temperature and density and that flow is laminar, such that kinetic energy has minimal impact on flow. Some models have been developed that incorporate the submerged conduits as a one-dimensional pipe network within the aquifer rather than as discrete, extremely transmissive features in a porous-equivalent medium; these submerged conduit models are usually referred to as hybrid models and may include the capability to simulate both laminar and turbulent flow in the one-dimensional pipe network. Comparisons of the application of a porous-equivalent media model with and without turbulence (MODFLOW-Conduit Flow Process mode 2 and basic MODFLOW, respectively) and a hybrid (MODFLOW-Conduit Flow Process mode 1) model to the Woodville Karst Plain near Tallahassee, Florida, indicated that for annual, monthly, or seasonal average hydrologic conditions, all methods met calibration criteria (matched observed groundwater levels and average flows). Thus, the increased effort required, such as the collection of data on conduit location, to develop a hybrid model and its increased computational burden, is not necessary for simulation of average hydrologic conditions (non-laminar flow effects on simulated head and spring discharge were minimal). However, simulation of a large storm event in the Woodville Karst Plain with daily stress periods indicated that turbulence is important for matching daily springflow hydrographs. Thus, if matching streamflow hydrographs over a storm event is required, the simulation of non-laminar flow and the location of conduits are required. The main challenge in application of the methods and approaches for developing hybrid models relates to the difficulty of mapping conduit networks or having high-quality datasets to calibrate these models. Additionally, hybrid models have long simulation times, which can preclude the use of parameter estimation for calibration. Simulation of contaminant transport that does not account for preferential flow through conduits or extremely permeable zones in any approach is ill-advised. Simulation results in other karst aquifers or other parts of the FAS may differ from the comparison demonstrated herein.

Florida↗

USGS2001

Updating flood maps efficiently using existing hydraulic models, very-high-accuracy elevation data, and a geographic information system: A pilot study on the Nisqually River, Washington

A method of updating flood inundation maps at a fraction of the expense of using traditional methods was piloted in Washington State as part of the U.S. Geological Survey Urban Geologic and Hydrologic Hazards Initiative. Large savings in expense may be achieved by building upon previous Flood Insurance Studies and automating the process of flood delineation with a Geographic Information System (GIS); increases in accuracy and detail result from the use of very-high-accuracy elevation data and automated delineation; and the resulting digital data sets contain valuable ancillary information such as flood depth, as well as greatly facilitating map storage and utility. The method consists of creating stage-discharge relations from the archived output of the existing hydraulic model, using these relations to create updated flood stages for recalculated flood discharges, and using a GIS to automate the map generation process. Many of the effective flood maps were created in the late 1970?s and early 1980?s, and suffer from a number of well recognized deficiencies such as out-of-date or inaccurate estimates of discharges for selected recurrence intervals, changes in basin characteristics, and relatively low quality elevation data used for flood delineation. FEMA estimates that 45 percent of effective maps are over 10 years old (FEMA, 1997). Consequently, Congress has mandated the updating and periodic review of existing maps, which have cost the Nation almost 3 billion (1997) dollars. The need to update maps and the cost of doing so were the primary motivations for piloting a more cost-effective and efficient updating method. New technologies such as Geographic Information Systems and LIDAR (Light Detection and Ranging) elevation mapping are key to improving the efficiency of flood map updating, but they also improve the accuracy, detail, and usefulness of the resulting digital flood maps. GISs produce digital maps without manual estimation of inundated areas between cross sections, and can generate working maps across a broad range of scales, for any selected area, and overlayed with easily updated cultural features. Local governments are aggressively collecting very-high-accuracy elevation data for numerous reasons; this not only lowers the cost and increases accuracy of flood maps, but also inherently boosts the level of community involvement in the mapping process. These elevation data are also ideal for hydraulic modeling, should an existing model be judged inadequate.

Washington↗