Search USGS⌕ Search

SEARCH · Search USGS

Results for “Natural Resource Modeling”

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

Snowpack signals in North American tree rings

Climate change has contributed to recent declines in mountain snowpack and earlier runoff, which in turn has intensified hydrological droughts in western North America. Climate model projections suggest that continued and severe snowpack reductions are expected over the 21st century, with profound consequences for ecosystems and human welfare. Yet the current understanding of trends and variability in mountain snowpack is limited by the relatively short and strongly temperature forced observational record. Motivated by the urgent need to better understand snowpack dynamics in a long-term, spatially coherent framework, here we examine snow-growth relationships in western North American tree-ring chronologies. We present an extensive network of snow-sensitive proxy data to support high space/time resolution paleosnow reconstruction, quantify and interpret the type and spatial density of snow related signals in tree-ring records, and examine the potential for regional bias in the tree-ring based reconstruction of different snow drought types (dry versus warm). Our results indicate three distinct snow-growth relationships in tree-ring chronologies: moisture-limited snow proxies that include a spring temperature signal, moisture-limited snow proxies lacking a spring temperature signal, and energy-limited snow proxies. Each proxy type is based on distinct physiological tree-growth mechanisms related to topographic and climatic site conditions, and provides unique information on mountain snowpack dynamics that can be capitalized upon within a statistical reconstruction framework. This work provides a platform and foundational background required for the accelerated production of high-quality annually-resolved snowpack reconstructions from regional to high (<12 km) spatial scales in western North America, and by extension, will support an improved understanding of the vulnerability of snowmelt-derived water resources to natural variability and future climate warming.

Environmental Research Letters↗

Social Values for Ecosystem Services (SolVES): using GIS to include social values information in ecosystem services assessments

Ecosystem services can be defined in various ways; simply put, they are the benefits provided by nature, which contribute to human well-being. These benefits can range from tangible products such as food and fresh water to cultural services such as recreation and esthetics. As the use of these benefits continues to increase, additional pressures are placed on the natural ecosystems providing them. This makes it all the more important when assessing possible tradeoffs among ecosystem services to consider the human attitudes and preferences that express underlying social values associated with their benefits. While some of these values can be accounted for through economic markets, other values can be more difficult to quantify, and attaching dollar amounts to them may not be very useful in all cases. Regardless of the processes or units used for quantifying such values, the ability to map them across the landscape and relate them to the ecosystem services to which they are attributed is necessary for effective assessments. To address some of the needs associated with quantifying and mapping social values for inclusion in ecosystem services assessments, scientists at the Rocky Mountain Geographic Science Center (RMGSC), in collaboration with Colorado State University, have developed a public domain tool, Social Values for Ecosystem Services (SolVES). SolVES is a geographic information system (GIS) application designed to use data from public attitude and preference surveys to assess, map, and quantify social values for ecosystem services. SolVES calculates and maps a 10-point Value Index representing the relative perceived social values of ecosystem services such as recreation and biodiversity for various groups of ecosystem stakeholders. SolVES output can also be used to identify and model relationships between social values and physical characteristics of the underlying landscape. These relationships can then be used to generate predicted Value Index maps for areas where survey data are not available. RMGSC will continue to develop more robust versions of SolVES by pursuing opportunities to work with land and resource managers as well as other researchers to apply SolVES to specific ecosystem management problems.

Fact Sheet↗

Comparison of two regression-based approaches for determining nutrient and sediment fluxes and trends in the Chesapeake Bay watershed

Nutrient and sediment fluxes and changes in fluxes over time are key indicators that water resource managers can use to assess the progress being made in improving the structure and function of the Chesapeake Bay ecosystem. The U.S. Geological Survey collects annual nutrient (nitrogen and phosphorus) and sediment flux data and computes trends that describe the extent to which water-quality conditions are changing within the major Chesapeake Bay tributaries. Two regression-based approaches were compared for estimating annual nutrient and sediment fluxes and for characterizing how these annual fluxes are changing over time. The two regression models compared are the traditionally used ESTIMATOR and the newly developed Weighted Regression on Time, Discharge, and Season (WRTDS). The model comparison focused on answering three questions: (1) What are the differences between the functional form and construction of each model? (2) Which model produces estimates of flux with the greatest accuracy and least amount of bias? (3) How different would the historical estimates of annual flux be if WRTDS had been used instead of ESTIMATOR? One additional point of comparison between the two models is how each model determines trends in annual flux once the year-to-year variations in discharge have been determined. All comparisons were made using total nitrogen, nitrate, total phosphorus, orthophosphorus, and suspended-sediment concentration data collected at the nine U.S. Geological Survey River Input Monitoring stations located on the Susquehanna, Potomac, James, Rappahannock, Appomattox, Pamunkey, Mattaponi, Patuxent, and Choptank Rivers in the Chesapeake Bay watershed. Two model characteristics that uniquely distinguish ESTIMATOR and WRTDS are the fundamental model form and the determination of model coefficients. ESTIMATOR and WRTDS both predict water-quality constituent concentration by developing a linear relation between the natural logarithm of observed constituent concentration and three explanatory variables—the natural log of discharge, time, and season. ESTIMATOR uses two additional explanatory variables—the square of the log of discharge and time-squared. Both models determine coefficients for variables for a series of estimation windows. ESTIMATOR establishes variable coefficients for a series of 9-year moving windows; all observed constituent concentration data within the 9-year window are used to establish each coefficient. Conversely, WRTDS establishes variable coefficients for each combination of discharge and time using only observed concentration data that are similar in time, season, and discharge to the day being estimated. As a result of these distinguishing characteristics, ESTIMATOR reproduces concentration-discharge relations that are closely approximated by a quadratic or linear function with respect to both the log of discharge and time. Conversely, the linear model form of WRTDS coupled with extensive model windowing for each combination of discharge and time allows WRTDS to reproduce observed concentration-discharge relations that are more sinuous in form. Another distinction between ESTIMATOR and WRTDS is the reporting of uncertainty associated with the model estimates of flux and trend. ESTIMATOR quantifies the standard error of prediction associated with the determination of flux and trends. The standard error of prediction enables the determination of the 95-percent confidence intervals for flux and trend as well as the ability to test whether the reported trend is significantly different from zero (where zero equals no trend). Conversely, WRTDS is unable to propagate error through the many (over 5,000) models for unique combinations of flow and time to determine a total standard error. As a result, WRTDS flux estimates are not reported with confidence intervals and a level of significance is not determined for flow-normalized fluxes. The differences between ESTIMATOR and WRTDS, with regard to model form and determination of model coefficients, have an influence on the determination of nutrient and sediment fluxes and associated changes in flux over time as a result of management activities. The comparison between the model estimates of flux and trend was made for combinations of five water-quality constituents at nine River Input Monitoring stations. The major findings with regard to nutrient and sediment fluxes are as follows: (1)WRTDS produced estimates of flux for all combinations that were more accurate, based on reduction in root mean squared error, than flux estimates from ESTIMATOR; (2) for 67 percent of the combinations, WRTDS and ESTIMATOR both produced estimates of flux that were minimally biased compared to observed fluxes(flux bias = tendency to over or underpredict flux observations); however, for 33 percent of the combinations, WRTDS produced estimates of flux that were considerably less biased (by at least 10 percent) than flux estimates from ESTIMATOR; (3) the average percent difference in annual fluxes generated by ESTIMATOR and WRTDS was less than 10 percent at 80 percent of the combinations; and (4) the greatest differences related to flux bias and annual fluxes all occurred for combinations where the pattern in observed concentration-discharge relation was sinuous (two points of inflection) rather than linear or quadratic (zero or one point of inflection). The major findings with regard to trends are as follows: (1) both models produce water-quality trends that have factored in the year-to-year variations in flow; (2) trends in water-quality condition are represented by ESTIMATOR as a trend in flow-adjusted concentration and by WRTDS as a flow normalized flux; (3) for 67 percent of the combinations with trend estimates, the WRTDS trends in flow-normalized flux are in the same direction and magnitude to the ESTIMATOR trends in flow-adjusted concentration, and at the remaining 33 percent the differences in trend magnitude and direction are related to fundamental differences between concentration and flux; and (4) the majority (85 percent) of the total nitrogen, nitrate, and orthophosphorus combinations exhibited long-term (1985 to 2010) trends in WRTDS flow-normalized flux that indicate improvement or reduction in associated flux and the majority (83 percent) of the total phosphorus (from 1985 to 2010) and suspended sediment (from 2001 to 2010) combinations exhibited trends in WRTDS flow-normalized flux that indicate degradation or increases in the flux delivered.

Delaware, Maryland, New York, Pennsylvania, Virgin↗

Analysis of ambient conditions and simulation of hydrodynamics, constituent transport, and water-quality characteristics in Lake Maumelle, Arkansas, 1991-92

Lake Maumelle is the major drinking-water source for the Little Rock metropolitan area in central Arkansas. Urban and agricultural development has increased in the Lake Maumelle Basin and information is needed related to constituent transport and water quality response to changes in constituent loading or hydrologic regime. This report characterizes ambient conditions in Lake Maumelle and its major tributary, Maumelle River; describes the calibration and verification of a numerical model of hydrodynamics and water quality; and provides several simulations that describe constituent transport and water quality response to changes in constituent loading and hydrologic regime. Ambient hydrologic and water-quality conditions demonstrate the relatively undisturbed nature of Lake Maumelle and the Maumelle River. Nitrogen and phosphorus concentrations were low, one to two orders of magnitude lower than estimates of national background nutrient concentrations. Phosphorus and chlorophyll a concentrations in Lake Maumelle demonstrate its oligotrophic/mesotrophic condition. However, concentrations of chlorophyll a appeared to increase since 1990 within the upper and middle reaches of the reservoir. A two-dimensional, laterally averaged hydrodynamic and water-quality model developed and calibrated for Lake Maumelle simulates water level, currents, heat transport and temperature distribution, conservative material transport, and the transport and transformation of 11 chemical constituents. Simulations included the movement and dispersion of spills or releases in the reservoir during stratified and unstratified conditions, release of the fish nursery pond off the southern shore of Lake Maumelle, and algal responses to changes in external loading. The model was calibrated using 1991 data and verified using 1992 data. Simulated temperature and dissolved oxygen concentrations related well when compared to measured values. Simulated nutrient and algal biomass also related reasonably well when compared to measured values. A simulated spill of conservative material at the upper end of Lake Maumelle during a major storm event took less than 102 hours to disperse the entire length of the reservoir. Simulation of a nursery pond release into a tributary to Lake Maumelle demonstrated how the released water plunges within the receiving embayment and enters the main stem of the reservoir at mid depths. Simulations of algal response to increases of nitrogen and phosphorus loads demonstrate the phosphorus limiting condition in Lake Maumelle. Results from this study will provide water-resource management with information to better understand how changes in hydrology and water quality in the basin affects water quality in the reservoir. With this information, managers will be able to more effectively manage their drinking-water source supply.

Arkansas↗

Inventories and mobilization of unsaturated zone sulfate, fluoride, and chloride related to land use change in semiarid regions, southwestern United States and Australia

Unsaturated zone salt reservoirs are potentially mobilized by increased groundwater recharge as semiarid lands are cultivated. This study explores the amounts of pore water sulfate and fluoride relative to chloride in unsaturated zone profiles, evaluates their sources, estimates mobilization due to past land use change, and assesses the impacts on groundwater quality. Inventories of water‐extractable chloride, sulfate, and fluoride were determined from borehole samples of soils and sediments collected beneath natural ecosystems ( N = 4), nonirrigated (“rain‐fed”) croplands ( N = 18), and irrigated croplands ( N = 6) in the southwestern United States and in the Murray Basin, Australia. Natural ecosystems contain generally large sulfate inventories (7800–120,000 kg/ha) and lower fluoride inventories (630–3900 kg/ha) relative to chloride inventories (6600–41,000 kg/ha). Order‐of‐magnitude higher chloride concentrations in precipitation and generally longer accumulation times result in much larger chloride inventories in the Murray Basin than in the southwestern United States. Atmospheric deposition during the current dry interglacial climatic regime accounts for most of the measured sulfate in both U.S. and Australian regions. Fluoride inventories are greater than can be accounted for by atmospheric deposition in most cases, suggesting that fluoride may accumulate across glacial/interglacial climatic cycles. Chemical modeling indicates that fluorite controls fluoride mobility and suggests that water‐extractable fluoride may include some fluoride from mineral dissolution. Increased groundwater drainage/recharge following land use change readily mobilized chloride. Sulfate displacement fronts matched or lagged chloride fronts by up to 4 m. In contrast, fluoride mobilization was minimal in all regions. Understanding linkages between salt inventories, increased recharge, and groundwater quality is important for quantifying impacts of anthropogenic activities on groundwater quality and is required for remediating salinity problems.

Water Resources Research↗

Interlaboratory comparability, bias, and precision for four laboratories measuring constituents in precipitation, November 1982-August 1983

Four laboratories were evaluated in their analysis of identical natural and simulated precipitation water samples. Interlaboratory comparability was evaluated using analysis of variance coupled with Duncan 's multiple range test, and linear-regression models describing the relations between individual laboratory analytical results for natural precipitation samples. Results of the statistical analyses indicate that certain pairs of laboratories produce different results when analyzing identical samples. Analyte bias for each laboratory was examined using analysis of variance coupled with Duncan 's multiple range test on data produced by the laboratories from the analysis of identical simulated precipitation samples. Bias for a given analyte produced by a single laboratory has been indicated when the laboratory mean for that analyte is shown to be significantly different from the mean for the most-probable analyte concentrations in the simulated precipitation samples. Ion-chromatographic methods for the determination of chloride, nitrate, and sulfate have been compared with the colorimetric methods that were also in use during the study period. Comparisons were made using analysis of variance coupled with Duncan 's multiple range test for means produced by the two methods. Analyte precision for each laboratory has been estimated by calculating a pooled variance for each analyte. Analyte estimated precisions have been compared using F-tests and differences in analyte precisions for laboratory pairs have been reported. (USGS)

Water-Resources Investigations Report↗

Geohydrology and simulation of ground-water flow in the carbonate rocks of the Valley Creek basin, eastern Chester County, Pennsylvania

Sixty-eight percent of the 22.6-square-mile Valley Creek basin is underlain by Cambrian and Ordovician limestone and dolomite. Ground water flows through a network of interconnected secondary openings; primary porosity is virtually nonexistent. Some of these openings have been enlarged by solution. Secondary porosity and permeability exhibit great spatial variability, and the yield and specific capacity of wells are highly variable. The number of water-bearing zones decreases with depth. Fifty percent of water-bearing zones are encountered within 100 feet of the land surface, and 81 percent are within 200 feet. Most ground-water flow in the Valley Creek basin is local and discharges to nearby streams. Ground-water discharge comprised an average of 76 percent of the flow of Valley Creek during 1983--87, including both natural ground- water discharge and quarry pumpage discharged to Valley Creek. Discharge from the Cedar Hollow quarry comprised 21 to 26 percent of the base flow of Valley Creek; the average was 23 percent. The average natural base flow of Valley Creek would be 8 percent lower if the quarry were not operating. Regional ground-water flow is to the northeast to the Schuylkill River. On the western side of the Valley Creek basin, the ground-water divide is 1/2 mile west of the surface-water divide. An estimated 0.75 million gallons per day of ground water flows from the adjacent West Valley Creek basn eastward into the Valley Creek basin. A ground-water divide is not present on the eastern side of the basin; the water table slopes gently eastward toward the Schuylkill River. On the northeaster side, an estimated 1.76 million gallons per day of ground water flows northeastward out of the basin to the Schuylkill River beneath the surface-water divide. On the southeaster side, an estimated 0.85 million gallons per day of ground water flows beneath the surface-water divide into the basin. Annual water budgets and an average water budget were calculated for 1983-87 for the 20.8-square-mile area bove the streamflow-gaging station. Annual precipitation for 1983-87 ranged from 40.61 to 56.55 inches and averaged 47.25 inches; annual streamflow ranged from 15.55 to 28.57 inches and averaged 22.31 inches; annual evapotranspiration ranged from 18.21 to 24.83 inches and averaged 22.90 inches; and annual recharge ranged from 15.89 to 26.84 inches and averaged 21.04 inches. The Valley Creek basin was modeled as a two-dimensional water-table aquifer. Recharge to, ground-water flow through, and discharge from the rocks of Chester valley were simulated. In order to include the natural hydrologic boundaries of the ground-water-flwo system, the 66.4-square-mile area between the Brandywine Creek and the Schuylkill River was modeled. The model was calibrated under stead-state conditions using avareage recharge and evapotranspiration rates. Aquifer hydraulic conductivity was estimated from specific-capacity and quifer-test data. The average (1983-87) annual water budget for hte Valley Creek basin was simualted. The effect of increased ground-water development on base flow and underflow was simulated by locating a hypothetical well field produceing 4 million gallons per day in different parts of the basin. Pumpage from a well field near surface-water divides would induce as much as an additional 1.41 inches per year of underflow from an adjacent surface-water basin. Pumpage from a well field near the center of the basin would affect base flow more than underflow. Increased seepage of ground water into quarries as a result of their expansion was simulated as increased withdrawal by pumping. A 100-percent increase in the pumping rate of the Cedar Hollow quarry, from 3.93 to 7.86 million gallons per day, owuld reduce the natural base flow of Valley Creek by 18 percent. However, the quarry pumpage would be discharged to Valley Creek, thereby increasing the base flow at the gaging station by

Pennsylvania↗

Water resources inventory of Connecticut Part 10: Lower Connecticut River basin

The lower Connecticut River basin study area in south-central Connecticut includes 639 square miles and is drained principally by the Connecticut River and by seven smaller streams that flow directly to Long Island Sound between the West River on the west and the Connecticut River on the east. The population in 1979 was estimated to be 210,380. Much of the industrial development and population centers are in the Mattabesset River basin in the northwestern part, and the largest water use is also in the Mattabesset River basin. Precipitation averages 47 inches per year and provides an abundant supply of water. About 20 inches returns to the atmosphere as evapotranspiration, and the remainder either flows directly to streams or percolates to the water table, eventually discharging to Long Island Sound. Small quantities of water are exported from the basin by the New Haven and Meridan Water Departments, and small quantities are imported by the New Britain Water Department and Metropolitan Direct Commission. Precipitation during 1931-60 resulted in an average annual runoff of 302 billion gallons. In inflow from the Connecticut River is added to the average annual runoff, the 4,370 billion gallon s per year is potentially available for water ue. The domestic, institutional, commercial, and industrial (other than cooling water) water use for 1970 was 7 billion gallons, which is only 3 percent of the total water used, whereas 97 percent of the total is cooling water for power plants. Approximately 60 percent of the 7 billion gallons is treated before being discharged back to the streams. The total amount of fresh water used during 1970 was estimated to be 256,000 million gallons (Mgal), of which 247,000 Mgal was used for cooling water at stream electric-generating plants. The quantity for domestic, commercial, industrial, and agricultural used was 9,000 Mgal, which was approximately 120 gallons a day per person. Public water systems providing 70 percent of these requirement and all the systems supplying water met the drinking water standards of the Connecticut General Assembly (1975). Till is widespread and generally provides only small amounts of water. Wells in till normally yield only a few hundred gallons of water daily and may be inadequate during dry periods. The thickness of of till ranges from 0 to 15 feet; a median thickness of 26 feet is estimated from information provided in drillers' logs of 467 wells penetrating underlying bedrock. The till is generally used only as an emergency or secondary source of water. Bedrock aquifers underlie the entire area and include sedimentary and crystalline (igneous and metamorphic) rock types. These aquifers supply small and usually reliable quantities of water to wells and are the chief source of water for many rural homes and farms., About 90 percent of the wells tapping bedrock yield at least 2 gal/min. The median yields from wells tapping aquifers in sedimentary, igneous, and metamorphic rocks are 11, 8, and 6.5 gal/min, respectively. The quantity of water potentially available from stratified drift was estimated on the basis of hydraulic characteristics of the aquifers, mathematical modeling of the aquifer system, and evaluation of natural and induced recharge. Long-term yields estimated or ten areas underlain by significant thickness of stratified drift range from 0.4 to 4.4 million gallons per day (Mgal/d). A change in well spacing or numbering could increase the long-term yields, but detailed modeling verification studies are needed to confirm optimal well locations. The chemical and physical (turbidity, color, taste, and sediment load) quality of water is good. The water if generally low in dissolved solids and is classified as soft to hard. Surface water is less mineralized than ground water, especially during high flow, when it is primarily derived from surface runoff rather than groundwater runoff. A median dissolved-solids concentration of 42 milligrams per liter (mg/L) and median hardness of 18 mg/L were determined from water samples collected from 26 streams during the high-flow period. During the low-flow period, median dissolved-solids concentration of 61 mg/L and median hardness of 27 mg/L were determined from sample from the same streams. The quality of water in stratified-drift and crystalline-rock aquifers is generally better than that in the sedimentary-rock aquifers. Water from 32 wells tapping stratified drift had median dissolved-solids concentrations of 116 mg/L; and 33 wells tapping stratified drift and 42 tapping crystalline rock had median hardness of 73 mg/L and 68 mg/L, respectively. Water from 32 wells tapping sedimentary rock had median dissolved concentrations of 231 and 156 mg/L, respectively. Sedimentary rock generally yields the hardest water. Iron and manganese occur objectionable concentrations in places, particularly in water from streams draining swamps and in water from aquifers either rich in iron and manganese-nearing minerals or where the reducing environment for solution of these minerals is favorable. Concentrations of iron in excess of 0.3 mg/L were found in 35 percent of the high streamflow samples, and in 45 percent of the ground-water samples. Most of the high iron and manganese concentration in streams and aquifers are found east of the Connecticut River. Human activities and tidal influence along the coast have modified the quality of water in much f the study area. The greatest influence from human activities has been in the northwestern part, in the Mattabesset River basin. There, the quality of water has been affected by domestic and animal wastes, which cause high dissolved-solids concentrations, high nitrate and phosphate loads and high bacterial counts. In the entire area, high nitrate in groundwater occur only locally, and its presence in an individual water supply is chiefly a function of its proximity to sources of contamination, of well construction, and of thickness of overburden. Thirty public-supply wells did have water that had high sodium concentrations or objectionable iron and manganese concentrations, but these are not considered health hazards in the concentrations found in the water samples. Streams, wetlands, and some aquifers along the sough boundary of the basin contain salty water because of tidal movement or extensive ground-water withdrawals. High sediment concentrations also occur as a result of tidal influence in this area.

Connecticut↗

The impacts of human recreation on brown bears (Ursus arctos): A review and new management tool

Increased popularity of recreational activities in natural areas has led to the need to better understand their impacts on wildlife. The majority of research conducted to date has focused on behavioral effects from individual recreations, thus there is a limited understanding of the potential for population-level or cumulative effects. Brown bears (Ursus arctos) are the focus of a growing wildlife viewing industry and are found in habitats frequented by recreationists. Managers face difficult decisions in balancing recreational opportunities with habitat protection for wildlife. Here, we integrate results from empirical studies with expert knowledge to better understand the potential population-level effects of recreational activities on brown bears. We conducted a literature review and Delphi survey of brown bear experts to better understand the frequencies and types of recreations occurring in bear habitats and their potential effects, and to identify management solutions and research needs. We then developed a Bayesian network model that allows managers to estimate the potential effects of recreational management decisions in bear habitats. A higher proportion of individual brown bears in coastal habitats were exposed to recreation, including photography and bear-viewing than bears in interior habitats where camping and hiking were more common. Our results suggest that the primary mechanism by which recreation may impact brown bears is through temporal and spatial displacement with associated increases in energetic costs and declines in nutritional intake. Killings in defense of life and property were found to be minimally associated with recreation in Alaska, but are important considerations in population management. Regulating recreation to occur predictably in space and time and limiting recreation in habitats with concentrated food resources reduces impacts on food intake and may thereby, reduce impacts on reproduction and survival. Our results suggest that decisions managers make about regulating recreational activities in time and space have important consequences for bear populations. The Bayesian network model developed here provides a new tool for managers to balance demands of multiple recreational activities while supporting healthy bear populations.

PLoS ONE↗

Drain tiles and groundwater resources: Understanding the relations

Executive Summary Drainage for agricultural production over the past 150 years has been an integral component of human-driven change to Minnesota’s rural landscapes. Benefits of drainage Historically, poorly drained soils across much of the State would often remain saturated or flooded after spring snowmelt, preventing timely farm operations such as tilling and planting crops (Arneman, 1963). Installation of agricultural drainage, both surface ditches and subsurface drainage, accelerated transport of water off farm fields and imparted producers higher crop yields (Beauchamp, 1987; Stoner and others, 1993). Agricultural drainage offered many other benefits such as preventing crop drown out, aerating the soil profile for improved plant growth, limiting surface runoff and soil erosion, and allowing farmers better access to croplands (Fausey and others, 1987). Without agricultural drainage on much of Minnesota’s croplands, it would have been difficult to realize high enough crop yields to remain economically viable. Environmental concerns While drainage of Minnesota croplands provided the benefits mentioned above, several environmental concerns result. These include wetland loss, degradation of downstream water quality, and reduced [potential for] recharge. Early agricultural drainage efforts (pre-20th century) led to the disappearance of much of Minnesota’s natural wetlands. Increased focus on preventing or mitigating wetland loss over the last 50 years has helped curtail further losses, even as agricultural drainage proceeds. Prior to establishment of Minnesota statehood, wetlands accounted for more than 10 million acres in Minnesota, including prairie wetlands, peatlands, and forest wetlands that comprised approximately 19 percent of the total land area (Palmer, 1915; King, 1980). In 2018, only half of Minnesota’s pre-settlement wetlands remain, mostly in parts of the State that have not experienced widespread drainage, such as northern Minnesota. Water-quality monitoring has shown that agricultural drainage, in particular the practice of subsurface drainage, provides a direct flow path for nutrient (nitrogen and soluble phosphorus) losses to surface water resources. The negative consequences of agricultural drainage on surface water quality are well documented (for example, Dinnes and others, 2002; Kladivko and others, 2004; Richards and others, 2008; Rozemeijer and others, 2010; Schottler and others, 2013). Agricultural basins with a high percentage of agricultural drainage have been implicated as part of the cause of the Gulf of Mexico hypoxia zone due to excessive nitrogen export (Goolsby and Battaglin, 2001; Randall and Mulla, 2001). The connection of hydrological effects of agricultural subsurface drainage on groundwater recharge and aquifers, on the other hand, has not been well-established. Agricultural subsurface drainage intercepts infiltrating water below croplands and directly discharges the water to nearby surface waters. However, the size of the water balance shift from drained water that would have evapotranspired or run off the land to drained water that would recharge underlying aquifers has been poorly characterized (Schuh, 2008). Drain Tiles and Groundwater Given the poor accounting of subsurface drainage effects on groundwater resources, the Minnesota Ground Water Association (MGWA) deemed it imperative that we document these effects so that groundwater resources in agricultural regions with substantial drainage can be effectively managed. This white paper documents the relations of drain tiles and groundwater resources and discusses the historical significance of agricultural drainage practices, the recognized positive benefits and potential negative consequences of agricultural drainage practices, and the gaps in understanding of the connections between agricultural drainage and groundwater resources. The major messages emerged from the findings of this white paper are: Complex history. Minnesota has a long history of agricultural drainage, spanning over 150 years. Agricultural drainage, and the eventual widespread usage of subsurface drainage, can be separated into at least four distinct periods of time: (a) early drainage to get water off the land, pre-20th century; (b) the boom and bust era (1900-1945); (c) postwar resurgence of subsurface drainage and early conservation efforts (1945-1960); and, (d) emergence of the environmental movement (1960 to present). The State’s regulatory framework that both allowed for drainage and controlled its usage during these periods is complicated and has been governed by a patchwork of both State and Federal statutes. Drainage Provinces. This white paper advances the concept of tile drainage provinces to aid in the discussion of regional differences in subsurface drainage and its overall effect on groundwater resources. Built upon the concept of groundwater provinces (Minnesota Department of Natural Resources, 2001), three distinct tile drainage provinces were conceptualized: (1) the Southeastern Province; (2) the SouthCentral Province; and, (3) the Western Province. The distinct geology and the soils that developed in these regions have implications for each region’s subsurface drainage density and the potential implications for groundwater. Knowledge gaps. Several critical knowledge gaps are identified in this paper, creating opportunities for further research to improve our understanding for better managed water resources: Extent of drainage is unknown. Direct estimates of the extent of subsurface drainage do not exist in Minnesota. However, several indirect methods have been utilized to estimate subsurface drainage, from the field-scale to countylevel through the use of geographic information system (GIS) analysis and aerial photography. Based on a 2012 U.S. Geological Survey estimate of subsurface drainage extent (Nakagaki and Wieczorek, 2016), about 21% of the land area in Minnesota has some density of subsurface drainage. Effect of drainage on underlying aquifers is unknown. A basic understanding of unconfined and confined aquifers and their recharge is necessary to connect any hydrological effects from agricultural drainage to groundwater. The basic goal of subsurface drainage to efficiently drain saturated soils clearly alters the water balance in croplands. However, its overall effect on groundwater resources has been poorly characterized, and is in large part determined by the geology below drained areas and the arrangement of underlying aquifers. Water balance shifts. An improved understanding of historical water balance shifts from pre- to post-drainage periods is necessary to understand long-term implications on net groundwater recharge. Also, more direct field-scale studies and indirect modeling studies are needed to characterize water budgets for fields with subsurface drainage.

Minnesota↗

Variability of bed mobility in natural, gravel‐bed channels and adjustments to sediment load at local and reach scales

Local variations in boundary shear stress acting on bed‐surface particles control patterns of bed load transport and channel evolution during varying stream discharges. At the reach scale a channel adjusts to imposed water and sediment supply through mutual interactions among channel form, local grain size, and local flow dynamics that govern bed mobility. In order to explore these adjustments, we used a numerical flow model to examine relations between model‐predicted local boundary shear stress (т j ( and measured surface particle size ( D 50 ) at bank‐full discharge in six gravel‐bed, alternate‐bar channels with widely differing annual sediment yields. Values of т j and D 50 were poorly correlated such that small areas conveyed large proportions of the total bed load, especially in sediment‐poor channels with low mobility. Sediment‐rich channels had greater areas of full mobility; sediment‐poor channels had greater areas of partial mobility; and both types had significant areas that were essentially immobile. Two reach‐mean mobility parameters (Shields stress and Q *) correlated reasonably well with sediment supply. Values which can be practicably obtained from carefully measured mean hydraulic variables and particle size would provide first‐order assessments of bed mobility that would broadly distinguish the channels in this study according to their sediment yield and bed mobility.

Water Resources Research↗

Effects of low-impact-development (LID) practices on streamflow, runoff quantity, and runoff quality in the Ipswich River Basin, Massachusetts: A summary of field and modeling studies

Low-impact-development (LID) approaches are intended to create, retain, or restore natural hydrologic and water-quality conditions that may be affected by human alterations. Wide-scale implementation of LID techniques may offer the possibility of improving conditions in river basins, such as the Ipswich River Basin in Massachusetts, that have run dry during the summer because of groundwater withdrawals and drought. From 2005 to 2008, the U.S. Geological Survey, in a cooperative funding agreement with the Massachusetts Department of Conservation and Recreation, monitored small-scale installations of LID enhancements designed to diminish the effects of storm runoff on the quantity and quality of surface water and groundwater. Funding for the studies also was contributed by the U.S. Environmental Protection Agency’s Targeted Watersheds Grant Program through a financial assistance agreement with Massachusetts Department of Conservation and Recreation. The monitoring studies examined the effects of replacing an impervious parking-lot surface with a porous surface on groundwater quality, installing rain gardens and porous pavement in a neighborhood of 3 acres on the quantity and quality of stormwater runoff, and installing a 3,000-ft2 (square-foot) green roof on the quantity and quality of rainfall-generated roof runoff. In addition to these small-scale installations, the U.S. Geological Survey’s Ipswich River Basin model was used to simulate the basin-wide effects on streamflow of several changes: broad-scale implementation of LID techniques, reduced water-supply withdrawals, and water-conservation measures. Water-supply and conservation scenarios for application in model simulations were developed with the assistance of two technical advisory committees that included representatives of State agencies responsible for water resources, the U.S. Environmental Protection Agency, the U.S. Geological Survey, water suppliers, and non-governmental organizations. From June 2005 to June 2007, groundwater quality was monitored at the Silver Lake town beach parking lot in Wilmington, Massachusetts, prior to and following the replacement of the conventional, impervious-asphalt surface with a porous surface consisting primarily of porous asphalt and porous pavers designed to enhance rainfall infiltration into the groundwater and to minimize runoff to Silver Lake. Concentrations of phosphorus, nitrogen, cadmium, chromium, copper, lead, nickel, zinc, and total petroleum hydrocarbons in groundwater were monitored. Enhancing infiltration of precipitation did not result in discernible increases in concentrations of these potential groundwater contaminants. Concentrations of dissolved oxygen increased slightly in groundwater profiles following the removal of the impervious asphalt parking-lot surface. In Wilmington, Massachusetts, in a 3-acre neighborhood, stormwater runoff volume and quality were monitored to determine the ability of selected LID enhancements (rain gardens and porous paving stones) to reduce flows and loads of the selected constituents to Silver Lake. Water-quality samples were analyzed for nutrients, metals, total petroleum hydrocarbons, and total-coliform and E. coli bacteria. A decrease in runoff quantity was observed for storms of 0.25 inch or less of precipitation. Water-quality-monitoring results were inconclusive; there were no statistically significant differences in concentrations or loads when the pre- and post-installation-period samples were compared. In a third field study, the characteristics of runoff from a vegetated "green" roof and a conventional, rubber-membrane roof were compared. The two primary factors affecting the green roof’s water-storage capacity were the amount of precipitation and antecedent dry period. Although concentrations of many of the chemicals in roof runoff were higher from the green roof than from the conventional roof, the ability of the green roof to retain water generally resulted in decreased differences between the total amounts (loads) of the chemicals that ran off the roofs. Land-use and water-management changes associated with LID implementation were investigated at multiple spatial scales, using the U.S. Geological Survey’s Ipswich River Basin model, to evaluate the effects of updated water-supply withdrawals for the towns of Reading and Wilmington (representing new baseline conditions for all simulations), potential land-use changes at buildout (potential future development), widespread implementation of retrofitting LID techniques, basin-scale water withdrawal reductions based on water-conservation pilot programs conducted by the Massachusetts Department of Conservation and Recreation, and land-use change and LID applications at a local scale. The new baseline simulation indicated that reduced water-supply withdrawals for the towns of Reading and Wilmington led to substantially higher medium and low flows in most of the reaches upstream from the South Middleton streamgage in the upper Ipswich River basin. Overall, simulations pointed to the importance of spatial scale in determining the effects of land-use change and LID practices on streamflow. Potential land-use changes at buildout had modest effects on streamflow in most subbasins (percent differences of less than 20 percent) because relatively little land in the basin was available for development. Results of the simulations conducted to evaluate widespread effective-impervious-area reductions upstream from the South Middleton streamgage indicated that the percentages of urban land use and associated effective impervious area were too small for even a 50-percent reduction of effective impervious area to appreciably affect streamflow in most subbasins. In contrast, the results of the hypothetical local-scale simulations indicated that for smaller streams, with high percentages of urban land use and associated effective impervious area, land-use change, development patterns, and LID practices may have substantial effects on streamflow. Modeling studies concurred with the results of fieldwork in the assessment that LID enhancements would likely have the greatest effect on decreasing stormwater runoff when broadly applied to highly impervious urban areas.

Massachusetts↗

Over, under, and through: Hydrologic connectivity and the future of coastal landscape salinization

Seawater intrusion (SWI) affects coastal landscapes worldwide. Here we describe the hydrologic pathways through which SWI occurs - over land via storm surge or tidal flooding, under land via groundwater transport, and through watersheds via natural and artificial surface water channels—and how human modifications to those pathways alter patterns of SWI. We present an approach to advance understanding of spatiotemporal patterns of salinization that integrates these hydrologic pathways, their interactions, and how humans modify them. We use examples across the East Coast of the United States that exemplify mechanisms of salinization that have been reported around the planet to illustrate how hydrologic connectivity and human modifications alter patterns of SWI. Finally, we suggest a path for advancing SWI science that includes (a) deploying standardized and well-distributed sensor networks at local to global scales that intentionally track SWI fronts, (b) employing remote sensing and geospatial imaging techniques targeted at integrating above and belowground patterns of SWI, and (c) continuing to develop data analysis and model-data fusion techniques to measure the extent, understand the effects, and predict the future of coastal salinization.

Water Resources Research↗

Methods for predicting peak discharge of floods caused by failure of natural and constructed earthen dams

Floods from failures of natural and constructed dams constitute a widespread hazard to people and property. Expeditious means of assessing flood hazards are necessary, particularly in the case of natural dams, which may form suddenly and unexpectedly. We revise statistical relations (derived from data for past constructed and natural dam failures) between peak discharge ( Q p ) and water volume released ( V 0 ) or drop in lake level ( d ) but assert that such relations, even when cast into a dimensionless form, are of limited utility because they fail to portray the effect of breach-formation rate. We then analyze a simple, physically based model of dam-breach formation to show that the hydrograph at the breach depends primarily on a dimensionless parameter η= kV 0 / g l/2 d 7/2 , where k is the mean erosion rate of the breach and g is acceleration due to gravity. The functional relationship between Q p and η takes asymptotically distinct forms depending on whether η ≪ 1 (relatively slow breach formation or small lake volume) or η ≫ 1 (relatively fast breach formation or large lake volume). Theoretical predictions agree well with data from dam failures for which k , and thus η, can be estimated. The theory thus provides a rapid means of predicting the plausible range of values of peak discharge at the breach in an earthen dam as long as the impounded water volume and the water depth at the dam face can be estimated.

Water Resources Research↗

The international river interface cooperative: Public domain flow and morphodynamics software for education and applications

This paper describes a new, public-domain interface for modeling flow, sediment transport and morphodynamics in rivers and other geophysical flows. The interface is named after the International River Interface Cooperative (iRIC), the group that constructed the interface and many of the current solvers included in iRIC. The interface is entirely free to any user and currently houses thirteen models ranging from simple one-dimensional models through three-dimensional large-eddy simulation models. Solvers are only loosely coupled to the interface so it is straightforward to modify existing solvers or to introduce other solvers into the system. Six of the most widely-used solvers are described in detail including example calculations to serve as an aid for users choosing what approach might be most appropriate for their own applications. The example calculations range from practical computations of bed evolution in natural rivers to highly detailed predictions of the development of small-scale bedforms on an initially flat bed. The remaining solvers are also briefly described. Although the focus of most solvers is coupled flow and morphodynamics, several of the solvers are also specifically aimed at providing flood inundation predictions over large spatial domains. Potential users can download the application, solvers, manuals, and educational materials including detailed tutorials at www.-i-ric.org . The iRIC development group encourages scientists and engineers to use the tool and to consider adding their own methods to the iRIC suite of tools.

Advances in Water Resources↗

The severity of the 2014–2015 snow drought in the Oregon Cascades in a multicentury context

The western United States (US) is a hotspot for snow drought. The Oregon Cascade Range is highly sensitive to warming and as a result has experienced the largest mountain snowpack losses in the western US since the mid-20th century, including a record-breaking snow drought in 2014–2015 that culminated in a state of emergency. While Oregon Cascade snowpacks serve as the state's primary water supply, short instrumental records limit water managers' ability to fully constrain long-term natural snowpack variability prior to the influence of ongoing and projected anthropogenic climate change. Here, we use annually-resolved tree-ring records to develop the first multi-century reconstruction of Oregon Cascade April 1st Snow Water Equivalent (SWE). The model explains 58% of observed snowpack variability and extends back to 1688 AD, nearly quintupling the length of the existing snowpack record. Our reconstruction suggests that only one other multiyear event in the last three centuries was as severe as the 2014–2015 snow drought. The 2015 event alone was more severe than nearly any other year in over three centuries. Extreme low-to-high snowpack “whiplash” transitions are a consistent feature throughout the reconstructed record. Multi-decadal intervals of persistent below-the-mean peak SWE are prominent features of pre-instrumental snowpack variability, but are generally absent from the instrumental period and likely not fully accounted for in modern water management. In the face of projected snow drought intensification and warming, our findings motivate adaptive management strategies that address declining snowpack and increasingly variable precipitation regimes.

Oregon↗

Porphyry copper deposit tract definition - A global analysis comparing geologic map scales

Geologic maps are a fundamental data source used to define mineral-resource potential tracts for the first step of a mineral resource assessment. Further, it is generally believed that the scale of the geologic map is a critical consideration. Previously published research has demonstrated that the U.S. Geological Survey porphyry tracts identified for the United States, which are based on 1:500,000-scale geology and larger scale data and published at 1:1,000,000 scale, can be approximated using a more generalized 1:2,500,000-scale geologic map. Comparison of the USGS porphyry tracts for the United States with weights-of-evidence models made using a 1:10,000,000-scale geologic map, which was made for petroleum applications, and a 1:35,000,000-scale geologic map, which was created as context for the distribution of porphyry deposits, demonstrates that, again, the USGS US porphyry tracts identified are similar to tracts defined on features from these small scale maps. In fact, the results using the 1:35,000,000-scale map show a slightly higher correlation with the USGS US tract definition, probably because the conceptual context for this small-scale map is more appropriate for porphyry tract definition than either of the other maps. This finding demonstrates that geologic maps are conceptual maps. The map information shown in each map is selected and generalized for the map to display the concepts deemed important for the map maker's purpose. Some geologic maps of small scale prove to be useful for regional mineral-resource tract definition, despite the decrease in spatial accuracy with decreasing scale. The utility of a particular geologic map for a particular application is critically dependent on the alignment of the intention of the map maker with the application. ?? International Association for Mathematical Geology 2007.

Natural Resources Research↗

Model analysis of effects on water levels at Indiana Dunes National Lakeshore caused by construction dewatering

Two computer models were developed to investigate possible hydrologic effects within the Indiana Dunes National Lakeshore caused by planned dewatering at the adjacent Bailly Nuclear Generator construction site. The dewatering, which is scheduled to last for about 18 months, would cause ground-water levels to be drawn down 20 to 30 feet (6 to 9 metres) to an elevation of 4 ft (1.2 m) below Lake Michigan datum under the approximately 7-acre (2.8-square hectometre) construction site, which is about 800 ft (245 m) west of the Lakeshore property. The results of this study have been used by the National Park Service to help evaluate the environmental impact of the dewatering, particularly the effects on the ecosystem of the lakeshore. The aquifer underlying the area is about 50 ft (15 m) thick and is composed of fine sand with layers of medium sand and gravel. All natural ponds in the area are separated from the aquifer by a low-permeability layer of silty organic muck and clay that ranges from about 1 to 4 ft (0.3 to 1.2 m) in thickness. This pond-bottom material acts only to retard the flow of water either to or from the aquifer beneath the ponds, depending upon the relative head in the pond and the aquifer. Elsewhere, the aquifer is under water-table conditions. The model analysis indicates that the planned dewatering would cause a drawdown of about 4 ft (1.2 m) under the westernmost pond of the Lakeshore and that this drawdown would cause the pond to go almost dry--less than 0.5 ft (0.15 m) of water remaining in about 1 percent of the pond--under average conditions during the 18-month dewatering period. When water levels are below average, as during late July and early August 1974, the pond would go dry in about 5 1/2 months. However, the pond may not have to go completely dry to damage the ecosystem. If the National Park Service's independent study determines the minimum pond level at which ecosystem damage would be minimized, the models developed in this study could be used to predict the hydrologic conditions necessary to maintain that level.

Indiana↗