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

Evaluating the spatial and temporal variability of groundwater uptake by riparian vegetation in a humid southeastern US catchment

In environments with shallow water tables, vegetation may use groundwater to support transpiration (T G ). This process has been carefully studied in some arid climates but rarely in humid climates—even those with severe droughts and seasonal water deficits. As such, the role of T G in humid-catchment hydrology is poorly constrained. We analysed water table fluctuations from nine monitoring wells along three transects in a second-order forested catchment to estimate T G at plot and whole-riparian zone scales. Average T G estimated around all well locations ranged from 1.06 to 4.95 mm d −1 and did not change systematically as a function of distance from stream channel or with plot-scale tree basal area. Counter to some previous studies, we found that T G was greater when the water table depth was deeper. Furthermore, the pattern of T G with water table depth was not monotonic at all locations. The ratio of T G to potential evapotranspiration tended to increase over the growing season, reflecting the progressive decrease in soil moisture storage and a greater reliance by vegetation on groundwater. Due to the lack of consistent spatial patterns in T G , we explored the number of monitoring wells needed to consistently estimate average T G within the 95% confidence bounds of the true mean. Based on this analysis, six or more wells were needed to consistently fall within the 95% confidence interval of the true mean. While this is based on the observed variability at a single site, it provides information for others considering this approach in similar upland forested catchments in humid regions.

Ecohydrology↗

Geochemical and geophysical indicators of oil and gas wastewater can trace potential exposure pathways following releases to surface waters

Releases of oil and gas (OG) wastewaters can have complex effects on stream-water quality and downstream organisms, due to sediment-water interactions and groundwater/surface water exchange. Previously, elevated concentrations of sodium (Na), chloride (Cl), barium (Ba), strontium (Sr), and lithium (Li), and trace hydrocarbons were determined to be key markers of OG wastewater releases when combined with Sr and radium (Ra) isotopic compositions. Here, we assessed the persistence of an OG wastewater spill in a creek in North Dakota using a combination of geochemical measurements and modeling, hydrologic analysis, and geophysical investigations. OG wastewater comprised 0.1 to 0.3% of the stream-water compositions at downstream sites in February and June 2015 but could not be quantified in 2016 and 2017. However, OG-wastewater markers persisted in sediments and pore water for 2.5 years after the spill and up to 7.2-km downstream from the spill site. Concentrations of OG wastewater constituents were highly variable depending on the hydrologic conditions. Electromagnetic measurements indicated substantially higher electrical conductivity under the bank adjacent to a seep 7.2 km downstream from the spill site. Geomorphic investigations revealed mobilization of sediment is an important contaminant transport process. Labile Ba, Ra, Sr, and ammonium (NH 4 ) concentrations extracted from sediments indicated sediments are a long-term reservoir of these constituents, both in the creek and on the floodplain. Using the drivers of ecological effects identified at this intensively studied site we identified 41 watersheds across the North Dakota landscape that may be subject to similar episodic inputs from OG wastewater spills. Effects of contaminants released to the environment during OG waste management activities remain poorly understood; however, analyses of Ra and Sr isotopic compositions, as well as trace inorganic and organic compound concentrations at these sites in pore-water provide insights into potentials for animal and human exposures well outside source-remediation zones.

North Dakota↗

Guidelines for determining flood flow frequency: Bulletin #17B of the Hydrology Subcommittee

In December 1967, Bulletin No. 15, "A Uniform Technique for Determining Flood Flow Frequencies," was issued by the Hydrology Committee of the Water Resources Council. The report recommended use of the Pearson Type III distribution with log transformation of the data (log-Pearson Type III distribution) as a base method for flood flow frequency studies. As pointed out in that report, further studies were needed covering various aspects of flow frequency determinations. In March 1976, Bulletin 17, "Guidelines for Determining Flood Flow Frequency" was issued by the Water Resources Council. The guide was an extension and update of Bulletin No. 15. It provided a more complete guide for flood flow frequency analysis incorporating currently accepted technical methods with sufficient detail to promote uniform application. It was limited to defining flood potentials in terms of peak discharge and exceedance probability at locations where a systematic record of peak flood flows is available. The recommended set of procedures was selected from those used or described in the literature prior to 1976, based on studies conducted for this purpose at the Center for Research in Water Resources of the University of Texas at Austin (summarized in Appendix 14) and on studies by the Work Group on Flood Flow Frequency. The "Guidelines" were revised and reissued in June 1977 as Bulletin 17A. Bulletin 17B is the latest effort to improve and expand upon the earlier publications. Bulletin 17B provides revised procedures for weighting a station skew value with the results from a generalized skew study, detecting and treating outliers, making two station comparisons, and computing confidence limits about a frequency curve. The Work Group that prepared this revision did not address the suitability of the original distribution or the generalized skew map. Major problems are encountered when developing guides for flood flow frequency determinations. There is no procedure or set of procedures that can be adopted which, when rigidly applied to the available data, will accurately define the flood potential of any given watershed. Statistical analysis alone will not resolve all flood frequency problems. As discussed in subsequent sections of this guide, elements of risk and uncertainty are inherent in any flood frequency analysis. User decisions must be based on properly applied procedures and proper interpretation of results considering risk and uncertainty. Therefore, the judgment of a professional experienced in hydrologic analysis will enhance the usefulness of a flood frequency analysis and promote appropriate application. It is possible to standarize many elements of flood frequency analysis. This guide describes each major element of the process of defining the flood potential at a specific location in terms of peak discharge and exceedance probability. Use is confined to stations where available records are adequate to warrant statistical analysis of the data. Special situations may require other approaches. In those cases where the procedures of this guide are not followed, deviations must be supported by appropriate study and accompanied by a comparison of results using the recommended procedures. As a further means of achieving consistency and improving results, the Work Group recommends that studies be coordinated when more than one analyst is working currently on data for the same location. This recommendation holds particularly when defining exceedance probabilities for rare events, where this guide allows more latitude. Flood records are limited. As more years of record become available at each location, the determination of flood potential may change. Thus, an estimate may be outdated a few years after it is made. Additional flood data alone may be sufficient reason for a fresh assessment of the flood potential. When making a new assessment, the analyst should incorporate in his study a review of earlier estimates. Where differences appear, they should be acknowledged and explained.

Bulletin↗

A one-dimensional, steady-state, dissolved-oxygen model and waste-load assimilation study for Little Laughery Creek, Ripley and Franklin counties, Indiana

The Indiana State Board of Health is developing a State water-quality management plan that includes establishing Limits for wastewater effluents discharged into Indiana streams. A digital model calibrated to conditions in Little Laughery Creek tributary and Little Laughery Creek was used to predict alternatives for future waste loadings that would be compatible with Indiana stream water-quality standards defined for two critical hydrologic conditions, summer and winter low flows. Natural streamflow during the summer and annual 7-day, 10-year low flow is zero. Headwater flow upstream from the wastewater-treatment facilities consists solely of process cooling water from an industrial discharger. This flow is usually Less than 0.5 cubic foot per second. Consequently, benefits from dilution are minimal. As a result, current and projected ammonia-nitrogen concentrations from the municipal discharges will result in in-stream ammonia-nitrogen concentrations that exceed the Indiana ammonia-nitrogen toxicity standards (maximum stream ammonia-nitrogen concentrations of 2.5 and 4.0 milligrams per liter during summer and winter low flows, respectively). This waste-Load assimilation study is not based on a verified model. The changes in stream water quality predicted by the model represent only possible stream responses to differing effluent conditions. Benthic-oxygen demand is probably the most significant factor affecting Little Laughery Creek and is probably responsible for the in-stream dis-solved-oxygen concentration being less than the Indiana stream dissolved-oxygen standard (5.0 milligrams per liter) during two water-quality surveys. After municipal dischargers complete advanced waste-treatment facilities, benthic-oxygen demand should be less significant in the stream dissolved-oxygen dynamics.

Indiana↗

Precipitation and runoff simulations of select perennial and ephemeral watersheds in the middle Carson River basin, Eagle, Dayton, and Churchill Valleys, west-central Nevada

The effect that land use may have on streamflow in the Carson River, and ultimately its impact on downstream users can be evaluated by simulating precipitation-runoff processes and estimating groundwater inflow in the middle Carson River in west-central Nevada. To address these concerns, the U.S. Geological Survey, in cooperation with the Bureau of Reclamation, began a study in 2008 to evaluate groundwater flow in the Carson River basin extending from Eagle Valley to Churchill Valley, called the middle Carson River basin in this report. This report documents the development and calibration of 12 watershed models and presents model results and the estimated mean annual water budgets for the modeled watersheds. This part of the larger middle Carson River study will provide estimates of runoff tributary to the Carson River and the potential for groundwater inflow (defined here as that component of recharge derived from percolation of excess water from the soil zone to the groundwater reservoir). The model used for the study was the U.S. Geological Survey's Precipitation-Runoff Modeling System, a physically based, distributed-parameter model designed to simulate precipitation and snowmelt runoff as well as snowpack accumulation and snowmelt processes. Models were developed for 2 perennial watersheds in Eagle Valley having gaged daily mean runoff, Ash Canyon Creek and Clear Creek, and for 10 ephemeral watersheds in the Dayton Valley and Churchill Valley hydrologic areas. Model calibration was constrained by daily mean runoff for the 2 perennial watersheds and for the 10 ephemeral watersheds by limited indirect runoff estimates and by mean annual runoff estimates derived from empirical methods. The models were further constrained by limited climate data adjusted for altitude differences using annual precipitation volumes estimated in a previous study. The calibration periods were water years 1980-2007 for Ash Canyon Creek, and water years 1991-2007 for Clear Creek. To allow for water budget comparisons to the ephemeral models, the two perennial models were then run from 1980 to 2007, the time period constrained somewhat by the later record for the high-altitude climate station used in the simulation. The daily mean values of precipitation, runoff, evapotranspiration, and groundwater inflow simulated from the watershed models were summed to provide mean annual rates and volumes derived from each year of the simulation. Mean annual bias for the calibration period for Ash Canyon Creek and Clear Creek watersheds was within 6 and 3 percent, and relative errors were about 18 and -2 percent, respectively. For the 1980-2007 period of record, mean recharge efficiency and runoff efficiency (percentage of precipitation as groundwater inflow and runoff) averaged 7 and 39 percent, respectively, for Ash Canyon Creek, and 8 and 31 percent, respectively, for Clear Creek. For this same period, groundwater inflow volumes averaged about 500 acre-feet for Ash Canyon and 1,200 acre-feet for Clear Creek. The simulation period for the ephemeral watersheds ranged from water years 1978 to 2007. Mean annual simulated precipitation ranged from 6 to 11 inches. Estimates of recharge efficiency for the ephemeral watersheds ranged from 3 percent for Eureka Canyon to 7 percent for Eldorado Canyon. Runoff efficiency ranged from 7 percent for Eureka Canyon and 15 percent at Brunswick Canyon. For the 1978-2007 period, mean annual groundwater inflow volumes ranged from about 40 acre-feet for Eureka Canyon to just under 5,000 acre-feet for Churchill Canyon watershed. Watershed model results indicate significant interannual variability in the volumes of groundwater inflow caused by climate variations. For most of the modeled watersheds, little to no groundwater inflow was simulated for years with less than 8 inches of precipitation, unless those years were preceded by abnormally high precipitation years with significant subsurface storage carryover.

Nevada↗

Assessment of long-term changes in surface-water extent within Klamath Marsh, south-central Oregon, 1985–2021

The annual maximum extent of surface water in Klamath Marsh has naturally fluctuated in response to periods of wet and dry conditions in the surrounding basin. Field observations during the 2010s indicate that the annual maximum extent of surface water has been declining and the marsh is not responding to hydrologic inputs as it had historically. This report describes the results of a hydrologic evaluation of Klamath Marsh to characterize and understand multi-year declines in the surface-water extent and increased intermittency of streamflow exiting the marsh. Landsat imagery collected during 1985–2021 was processed to create a time series of annual maximum surface-water extent to assess changes in surface-water inundation within the marsh. A 50-percent decrease in the mean surface area of annual total open-water extent (TOWE) during the latter half of the study period (2003–21) compared to the first half (1985–2003) was observed in this 37-year time-series dataset. The change in open-water extent was offset by a corresponding increase in dry land in the marsh. Time series of streamflow, groundwater level, total annual precipitation, annual mean temperature, and anthropogenic water use and water management were compiled and evaluated to improve understanding of the factors affecting TOWE. Statistically significant downward trends in the regional groundwater table and streamflow into and out of the marsh were identified as well as statistically significant upward trends in annual mean temperature. Statistically significant correlations among TOWE, streamflow, and groundwater level also were identified. The decreasing trends could not be attributed to changes in total annual precipitation or changing anthropogenic groundwater use within the study area. Declines in the open-water extent of Klamath Marsh since 2000 principally are due to a decoupling of the groundwater and surface-water system beneath the marsh because of regional declines in groundwater level. Regional increases in air temperature and the reestablishment of more than 55,000 acres of forested land within the study area have likely contributed to increasing evapotranspiration, leaving less water available for groundwater recharge and stream base flow and resulting in basin-wide declines in streamflow and groundwater levels.

Oregon↗

A new groundwater energy transport model for the MODFLOW hydrologic simulator

Heat transport in the subsurface is an important aspect of research related to the effects of a warming climate on ecological services (i.e., cold-water refugia); the development of geothermal resources for energy banking schemes (i.e., aquifer thermal energy storage [ATES]); and the effects of temperature on other aspects of groundwater quality, such as nutrient cycling. Historically, simulation of heat transport using the MODFLOW groundwater simulator and related codes was performed by scaling the input parameters of a solute-transport model to emulate heat transport. However, that approach required additional pre- and post-processing of input and output and could not account for the variation in effective thermal storage and transport properties during transient, unsaturated flow, for example. True heat-transport capabilities in the context of MODFLOW were first introduced in a variant called USG-Transport. More recently, a new groundwater energy-transport (GWE) model type has been added to MODFLOW 6, the core version of the MODFLOW hydrologic simulator. GWE supports the simulation of heat transport on structured or unstructured grids as well as within and between features of advanced packages that represent streams, lakes, multi-aquifer wells, and the unsaturated zone. GWE is integrated within MODFLOW 6 and is accessible through the FloPy Python package and the MODFLOW 6 application programming interface (API). An example simulation demonstrates conduction between grid cells through both the water and the solid aquifer material, including thermal bleeding from saturated overburden cells into a groundwater flow field.

Groundwater↗

Processing of next generation weather radar-multisensor precipitation estimates and quantitative precipitation forecast data for the DuPage County streamflow simulation system

The U.S. Geological Survey, in cooperation with DuPage County Stormwater Management Department, is testing a near real-time streamflow simulation system that assists in the management and operation of reservoirs and other flood-control structures in the Salt Creek and West Branch DuPage River drainage basins in DuPage County, Illinois. As part of this effort, the U.S. Geological Survey maintains a database of hourly meteorological and hydrologic data for use in this near real-time streamflow simulation system. Among these data are next generation weather radar-multisensor precipitation estimates and quantitative precipitation forecast data, which are retrieved from the North Central River Forecasting Center of the National Weather Service. The DuPage County streamflow simulation system uses these quantitative precipitation forecast data to create streamflow predictions for the two simulated drainage basins. This report discusses in detail how these data are processed for inclusion in the Watershed Data Management files used in the streamflow simulation system for the Salt Creek and West Branch DuPage River drainage basins.

Illinois↗

A method for quantifying cloud immersion in a tropical mountain forest using time-lapse photography

Quantifying the frequency, duration, and elevation range of fog or cloud immersion is essential to estimate cloud water deposition in water budgets and to understand the ecohydrology of cloud forests. The goal of this study was to develop a low-cost and high spatial-coverage method to detect occurrence of cloud immersion within a mountain cloud forest by using time-lapse photography. Trail cameras and temperature/relative humidity sensors were deployed at five sites covering the elevation range from the assumed lifting condensation level to the mountain peaks in the Luquillo Mountains of Puerto Rico. Cloud-sensitive image characteristics (contrast, the coefficient of variation and the entropy of pixel luminance, and image colorfulness) were used with a k -means clustering approach to accurately detect cloud-immersed conditions in a time series of images from March 2014 to May 2016. Images provided hydrologically meaningful cloud-immersion information while temperature-relative humidity data were used to refine the image analysis using dew point information and provided temperature gradients along the elevation transect. Validation of the image processing method with human-judgment based classification generally indicated greater than 90% accuracy. Cloud-immersion frequency averaged 80% at sites above 900 m during nighttime hours and 49% during daytime hours, and was consistent with diurnal patterns of cloud immersion measured in a previous study. Results for the 617 m site demonstrated that cloud immersion in the Luquillo Mountains rarely occurs at the previously-reported cloud base elevation of about 600 m (11% during nighttime hours and 5% during daytime hours). The framework presented in this paper will be used to monitor at a low cost and high spatial resolution the long-term variability of cloud-immersion patterns in the Luquillo Mountains, and can be applied to ecohydrology research at other cloud-forest sites or in coastal ecosystems with advective sea fog.

El Yunque National Forest, Puerto Rico↗

Can the planetary health concept save freshwater biodiversity and ecosystems?

People clearly need and benefit from healthy freshwater ecosystems; Given the precarious state of these important systems and services, current efforts to address the freshwater biodiversity crisis remain insufficient. Planetary health is an emerging framework that aims to secure the state of natural systems within environmental limits that ensure humanity can flourish. The planetary health concept is tied to the planetary boundaries framework in which various ecological thresholds are identified with the goal of constraining human activity to within those boundaries (so-called safe operating spaces). Freshwater systems are influenced by some planetary-scale processes like the climate systems and phosphorus and nitrogen cycles. Nonetheless, safe boundaries to guide the conservation and management of freshwater ecosystems need to consider their uneven distribution around the globe, and their ecological and hydrologic limits, which are often site and context dependent. Efforts to down-scale planetary boundaries concepts to the management of freshwater recreational fisheries at the lake scale, suggest that there are opportunities for rethinking planetary health as a nested cross-scale approach from the planet to the watershed.

Lancet Planetary Health↗

Reservoir stratification modulates the influence of impoundments on fish mercury concentrations along an arid land river system

Impoundment is among the most common hydrologic alterations with impacts on aquatic ecosystems that can include effects on mercury (Hg) cycling. However, landscape-scale differences in Hg bioaccumulation between reservoirs and other habitats are not well characterized nor are the processes driving these differences. We examined total Hg (THg) concentrations of Smallmouth Bass ( Micropterus dolomieu ) collected from reservoir, tailrace, and free-flowing reaches along an 863 km segment of the Snake River, USA, a semiarid river with 22 impoundments along its course. Across three size-classes (putative 1-year-old, first reproductive, and harvestable sized fish), THg concentrations in reservoirs and tailraces averaged 76% higher than those in free-flowing segments. Among reservoirs, THg concentrations were highest in reservoirs with inconsistent stratification patterns, 47% higher than annually stratified, and 144% higher than unstratified reservoirs. Fish THg concentrations in tailraces immediately downstream of stratified reservoirs were higher than those below unstratified (38–130%) or inconsistently stratified (32–79%) reservoirs. Stratification regimes influenced the exceedance of fish and human health benchmarks, with 52–80% of fish from stratifying reservoirs and downstream tailraces exceeding a human consumption benchmark, compared to 6–17% where stratification did not occur. These findings suggest that impoundment and stratification play important roles in determining the patterns of Hg exposure risk across the landscape.

Environmental Science & Technology↗

Snowmelt-triggered earthquake swarms at the margin of Long Valley Caldera, California

Fluids are well known to influence earthquakes, yet rarely are earthquakes convincingly linked to precipitation. Weak modulation or limited data often leads to ambiguous interpretations. In contrast, here we find that shallow seismicity in the Sierra Nevada range near Long Valley Caldera is strongly modulated by snowmelt. Over 33 years, shallow seismicity rates were ~37 times higher during very wet periods versus very dry periods. Relative earthquake relocations from a swarm in 2017 reveal downward migration from ~1- to 3-km depth along a steeply inclined plane. Steeply dipping strata may provide high-permeability pathways and faulting plane. Here we combine the correlated seismicity and hydrologic time series with the propagation observed in the relatively relocated earthquakes. From this combined evidence, we infer that pressure diffusion from groundwater recharge dramatically accelerated shallow seismicity rates, causing seismic swarms unrelated to volcanic processes.

California↗

Capturing the transient hydrological response in sandy soils during a rare cloudburst associated with shallow slope failures; A case study in the Atlantic Highlands, New Jersey, USA

A cloudburst on 7 August 2018 in the coastal bluffs of the Atlantic Highlands, New Jersey, induced flooding, erosion and multiple shallow slope failures that adversely affected the surrounding hillside residential area. Historically, short-duration deluges are rare in the New York Bay region, with only eight cloudbursts of greater magnitude documented since 1948. The coastal bluffs consist of a variably thick, sandy surficial material overlying flat-lying, mostly non-indurated Cretaceous and Tertiary sediments, including some low-permeability glauconitic units. The bluffs have been affected by both historical deep-seated and shallow landslide movement, the latter typically related to heavy, relatively long-duration rainfall associated with tropical cyclones and nor'easters. The shallow hydrological response during the rare cloudburst was captured at two hydrological monitoring sites and yielded insights into rapidly changing moisture conditions resulting in slope failure. Additional information is provided on historical cloudbursts that have affected the region, antecedent moisture conditions, and documented landslide types and processes.

New Jersey↗

Hierarchical stochastic modeling of large river ecosystems and fish growth across spatio-temporal scales and climate models: the Missouri River endangered pallid sturgeon example

We present a hierarchical series of spatially decreasing and temporally increasing models to evaluate the uncertainty in the atmosphere – ocean global climate model (AOGCM) and the regional climate model (RCM) relative to the uncertainty in the somatic growth of the endangered pallid sturgeon (Scaphirhynchus albus). For effects on fish populations of riverine ecosystems, cli- mate output simulated by coarse-resolution AOGCMs and RCMs must be downscaled to basins to river hydrology to population response. One needs to transfer the information from these climate simulations down to the individual scale in a way that minimizes extrapolation and can account for spatio-temporal variability in the intervening stages. The goal is a framework to determine whether, given uncertainties in the climate models and the biological response, meaningful inference can still be made. The non-linear downscaling of climate information to the river scale requires that one realistically account for spatial and temporal variability across scale. Our down- scaling procedure includes the use of fixed/calibrated hydrological flow and temperature models coupled with a stochastically parameterized sturgeon bioenergetics model. We show that, although there is a large amount of uncertainty associated with both the climate model output and the fish growth process, one can establish significant differences in fish growth distributions between models, and between future and current climates for a given model.

Missouri River basin↗

Assessing carbon stocks, carbon sequestration, and greenhouse-gas fluxes in ecosystems of the United States under present conditions and future scenarios

The Energy Independence and Security Act of 2007 (EISA) requires the U.S. Department of the Interior (DOI) to develop a methodology and conduct an assessment of carbon storage, carbon sequestration, and greenhouse-gas (GHG) fluxes in the Nation's ecosystems. The U.S. Geological Survey (USGS) has developed and published the methodology (U.S. Geological Survey Scientific Investigations Report 2010-5233) and has assembled an interdisciplinary team of scientists to conduct the assessment over the next three to four years, commencing in October 2010. The assessment will fulfill specific requirements of the EISA by (1) quantifying, measuring, and monitoring carbon sequestration and GHG fluxes using national datasets and science tools such as remote sensing, and biogeochemical and hydrological models, (2) evaluating a range of management and restoration activities for their effects on carbon-sequestration capacity and the reduction of GHG fluxes, and (3) assessing effects of climate change and other controlling processes (including wildland fires) on carbon uptake and GHG emissions from ecosystems.

Fact Sheet↗

Modflow-setup: Robust automation of groundwater model construction

In an age of both big data and increasing strain on water resources, sound management decisions often rely on numerical models. Numerical models provide a physics-based framework for assimilating and making sense of information that by itself only provides a limited description of the hydrologic system. Often, numerical models are the best option for quantifying even intuitively obvious connections between human activities and water resource impacts. However, despite many recent advances in model data assimilation and uncertainty quantification, the process of constructing numerical models remains laborious, expensive, and opaque, often precluding their use in decision making. Modflow-setup aims to provide rapid and consistent construction of MODFLOW groundwater models through robust and repeatable automation. Common model construction tasks are distilled in an open-source, online code base that is tested and extensible through collaborative version control. Input to Modflow-setup consists of a single configuration file that summarizes the workflow for building a model, including source data, construction options, and output packages. Source data providing model structure and parameter information including shapefiles, rasters, NetCDF files, tables, and other (geolocated) sources to MODFLOW models are read in and mapped to the model discretization, using Flopy and other general open-source scientific Python libraries. In a few minutes, an external array-based MODFLOW model amenable to parameter estimation and uncertainty quantification is produced. This paper describes the core functionality of Modflow-setup, including a worked example of a MODFLOW 6 model for evaluating pumping impacts to a lake in central Wisconsin, United States.

Wisconsin↗

Tundra landform and vegetation productivity trend maps for the Arctic Coastal Plain of northern Alaska

Arctic tundra landscapes are composed of a complex mosaic of patterned ground features, varying in soil moisture, vegetation composition, and surface hydrology over small spatial scales (10–100 m). The importance of microtopography and associated geomorphic landforms in influencing ecosystem structure and function is well founded, however, spatial data products describing local to regional scale distribution of patterned ground or polygonal tundra geomorphology are largely unavailable. Thus, our understanding of local impacts on regional scale processes (e.g., carbon dynamics) may be limited. We produced two key spatiotemporal datasets spanning the Arctic Coastal Plain of northern Alaska (~60,000 km 2 ) to evaluate climate-geomorphological controls on arctic tundra productivity change, using (1) a novel 30 m classification of polygonal tundra geomorphology and (2) decadal-trends in surface greenness using the Landsat archive (1999–2014). These datasets can be easily integrated and adapted in an array of local to regional applications such as (1) upscaling plot-level measurements (e.g., carbon/energy fluxes), (2) mapping of soils, vegetation, or permafrost, and/or (3) initializing ecosystem biogeochemistry, hydrology, and/or habitat modeling.

Alaska↗

Simulation of ground-water flow and solute transport in the Glen Canyon aquifer, East-Central Utah

The extraction of methane from coal beds in the Ferron coal trend in central Utah started in the mid-1980s. Beginning in 1994, water from the extraction process was pressure injected into the Glen Canyon aquifer. The lateral extent of the aquifer that could be affected by injection is about 7,600 square miles. To address regional-scale effects of injection over a decadal time frame, a conceptual model of ground-water movement and transport of dissolved solids was formulated. A numerical model that incorporates aquifer concepts was then constructed and used to simulate injection. The Glen Canyon aquifer within the study area is conceptualized in two parts—an active area of ground-water flow and solute transport that exists between recharge areas in the San Rafael Swell and Desert, Waterpocket Fold, and Henry Mountains and discharge locations along the Muddy, Dirty Devil, San Rafael, and Green Rivers. An area of little or negligible ground-water flow exists north of Price, Utah, and beneath the Wasatch Plateau. Pressurized injection of coal-bed methane production water occurs in this area where dissolved-solids concentrations can be more than 100,000 milligrams per liter. Injection has the potential to increase hydrologic interaction with the active flow area, where dissolved-solids concentrations are generally less than 3,000 milligrams per liter. Pressurized injection of coal-bed methane production water in 1994 initiated a net addition of flow and mass of solutes into the Glen Canyon aquifer. To better understand the regional scale hydrologic interaction between the two areas of the Glen Canyon aquifer, pressurized injection was numerically simulated. Data constraints precluded development of a fully calibrated simulation; instead, an uncalibrated model was constructed that is a plausible representation of the conceptual flow and solute-transport processes. The amount of injected water over the 36-year simulation period is about 25,000 acre-feet. As a result, simulated water levels in the injection areas increased by 50 feet and dissolved-solids concentrations increased by 100 milligrams per liter or more. These increases are accrued into aquifer storage and do not extend to the rivers during the 36-year simulation period. The amount of change in simulated discharge and solute load to the rivers is less than the resolution accuracy of the numerical simulation and is interpreted as no significant change over the considered time period.

Utah↗