SEARCH · USGS Science

Results for “Agricultural Water Management”

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

SourceResearch connectionFrom yearThrough year

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

1,439 records · Page 80Linked to original sources

USGS2024

Simulation of groundwater flow and brine discharge to the Dolores River in the Paradox Valley, Montrose County, Colorado

Salinity, or total dissolved solids (TDS), of the Colorado River affects agricultural, municipal, and industrial water users and is an important concern in the Western United States. In the Paradox Valley of southwestern Colorado, natural discharge of sodium-chloride brine to the Dolores River from the underlying core of a salt-valley anticline accounts for about 6 percent of the salinity load to the Colorado River. Formation of the Paradox Valley began during the Miocene, and subsequent erosion exposed the Pennsylvania Paradox Formation in the core of the anticline where a cap rock, collapse features, breccia, and sodium-chloride saturated brine developed at the top of the exposed salt diapir. The discharge of brine to the Dolores River is affected by these dissolution features, along with seasonal hydrologic conditions and density-dependent flow between older dense brine and the younger fresh groundwater in the overlying alluvial aquifer. To reduce TDS concentrations in the Dolores River through the Paradox Valley, the Bureau of Reclamation has pumped brine from a series of shallow wells adjacent to the river since July 1996. The pumped brine is collected and piped to a deep disposal well where it is injected into the Mississippian Leadville Limestone at a depth of about 4,570-meters below land surface. The pumping and injection operation is collectively known as the Paradox Valley Unit (PVU), and by 2015, the PVU had substantially reduced TDS concentrations in the Dolores River by about 70 percent. Since 2019, injection-pressure limits and related seismic activity have constrained deep-well injection and thus brine pumping at the PVU. In cooperation with the Bureau of Reclamation, the U.S. Geological Survey developed a MODFLOW-6 three-dimensional, variable-density groundwater flow and TDS transport model of the Paradox Valley to evaluate the effects of PVU pumping operations on brine discharge to the Dolores River and to guide additional research. The finite-difference model grid consists of 76 rows and 48 columns oriented from northwest to southeast in alignment with valley topography and groundwater-flow directions in the near-surface freshwater alluvial aquifer. A 7-layer hydrogeologic framework was developed from existing datasets to represent the alluvial aquifer, cap rock, collapse breccia, and groundwater flow and TDS transport from the underlying Paradox Formation salt to the Dolores River. The model represents a 33-year transient calibration period from 1987 through 2020 that includes pre-PVU conditions from 1987 through June 1996 and post-PVU conditions from July 1996 through 2020. A 1,000-year simulation of groundwater flow and coupled TDS transport computed the initial conditions for the subsequent 33-year transient simulation. Observations of precipitation, streamflow, evaporation, agricultural land use, and PVU brine pumping rates were used to specify appropriate boundary conditions to the model representing time-varying recharge, tributary streamflow, groundwater underflow, evapotranspiration (ET), and PVU pumping. Values for average monthly streamflow and TDS concentration at the upstream streamgage, the Dolores River at Bedrock (USGS streamgage 09169500), were specified as model input where the Dolores River enters Paradox Valley. Observed pumping from the PVU, water levels and TDS concentrations in groundwater, and streamflow and estimated TDS concentrations at the downstream streamgage, the Dolores River near Bedrock (USGS streamgage 09171100), were calibration targets that constrained the manual calibration of model parameters representing aquifer hydraulic conductivity, storage, streambed conductance, recharge, and (ET). Two primary model-calibration targets were the match between observed and simulated TDS mass flux from PVU pumping wells and the match between estimated and simulated TDS mass flux to the Dolores River. The simulated TDS mass withdrawn by pumping wells is calculated by the model as the product of the assigned pumping rate and simulated groundwater TDS concentrations. Because actual pumping rates were assigned as simulated values, the total simulated PVU pumping for the 33-year calibration is within 0.5 percent of the observed values. However, simulated concentrations and thus mass flux of TDS withdrawn by the PVU pumping wells were consistently about 26 percent less than observed values for all the simulated time periods (33-year simulation, pre-PVU, and post-PVU). The representation of brine inflow was explored through additional modeling to evaluate the effect of the simulated brine source on groundwater TDS concentrations. Results indicated that a saturated-salt constant-flux brine source best replicated the magnitude and transient pattern observed for TDS mass flux from PVU pumping wells. The simulated TDS mass flux to the Dolores River is compared to estimates based on observed streamflow and specific conductance (SC) data for the downstream streamgage. The calibrated model provided a close fit of simulated to measured streamflow at the downstream streamgage, and the calibrated model fit to estimated TDS concentrations at the downstream streamgage was reasonable. The greatest differences between simulated and estimated values occurred during drought periods from June 2000 to March 2003, May 2012 to June 2013, and October 2013 to October 2014, when simulated TDS concentrations in the river were greater than estimated concentrations. In general, simulated TDS mass flux to the river for the pre-PVU period is in good agreement with estimated values (2-percent difference), but the model overestimated TDS mass flux to the river by about 41 percent during the post-PVU period. The model uncertainty with respect to TDS mass flux to the river indicates other processes or model parameters not well represented by the model are affecting the system, especially during drought. During model calibration, the most sensitive parameters were identified as vertical hydraulic conductivity of the alluvial aquifer, conductance of the Dolores River streambed, ET extinction depth and rate, and recharge rate. Five 5-year scenarios of conditions for 2021–25 were simulated to assist evaluation of alternative strategies to manage the discharge of brine into the Dolores River. The first scenario simulates no PVU pumping and serves as a base case for comparison to the other scenarios. Two scenarios simulate the effects of varying withdrawal timing at an annual rate about one-third less than during 2010 through 2018. During high-flow spring snowmelt runoff periods when brine discharge is naturally minimized, PVU pumping does not substantially affect salinity in the Dolores River, and comparison of these two scenarios indicates that scheduling brine withdrawals during times of low river stage is nearly as effective at reducing TDS mass flux to the river as pumping brine year-round. Cessation of pumping during periods of high river stage may be advantageous for system maintenance, brine injection, and seismic-risk reduction. The fourth scenario tested the effect of reducing irrigation-return flow on brine discharge and predicted a slight reduction of TDS mass flux to the Dolores River, but not as great a reduction as that of using the PVU to remove brine. The fifth scenario simulated 5 years of drought conditions without PVU pumping and indicates brine discharge during drought about 15 percent greater than during average hydrologic conditions. Results from scenario 5 are consistent with the calibrated model results and indicate that aquifer properties and ET processes and parameters may be affecting simulation results during drought. The Paradox Valley groundwater model provides a reasonable overall match to observed conditions in the Dolores River. The model is useful for evaluating relative differences between brine management scenarios to inform PVU operational decisions and to identify gaps in data and process understanding. Representation of the brine source, hydraulic-conductivity parameters, and recharge and ET processes were identified as potential areas for additional field and modeling research. Additional research in the Paradox Valley might include field-data collection that provides additional information on the hydrogeologic framework, groundwater levels, groundwater TDS concentrations, stream characteristics, and aquifer properties. Additional modeling efforts could benefit from applying advanced tools for model development, calibration, and visualization including parameter-estimation and sensitivity analysis. Statistical evaluation of known model uncertainties such as hydraulic conductivity, streambed conductance, representations of the brine source, recharge, and ET could improve the match between simulated and estimated TDS mass flux from PVU pumping wells and to the Dolores River further informing model predictions and system understanding for the Paradox Valley.

Colorado↗

USGS2013

Estimating nitrate concentrations in groundwater at selected wells and springs in the surficial aquifer system and Upper Floridan aquifer, Dougherty Plain and Marianna Lowlands, Georgia, Florida, and Alabama, 2002-50

Groundwater from the surficial aquifer system and Upper Floridan aquifer in the Dougherty Plain and Marianna Lowlands in southwestern Georgia, northwestern Florida, and southeastern Alabama is affected by elevated nitrate concentrations as a result of the vulnerability of the aquifer, irrigation water-supply development, and intensive agricultural land use. The region relies primarily on groundwater from the Upper Floridan aquifer for drinking-water and irrigation supply. Elevated nitrate concentrations in drinking water are a concern because infants under 6 months of age who drink water containing nitrate concentrations above the U.S. Environmental Protection Agency maximum contaminant level of 10 milligrams per liter as nitrogen can become seriously ill with blue baby syndrome. In response to concerns about water quality in domestic wells and in springs in the lower Apalachicola–Chattahoochee–Flint River Basin, the Florida Department of Environmental Protection funded a study in cooperation with the U.S. Geological Survey to examine water quality in groundwater and springs that provide base flow to the Chipola River. A three-dimensional, steady-state, regional-scale groundwater-flow model and two local-scale models were used in conjunction with particle tracking to identify travel times and areas contributing recharge to six groundwater sites—three long-term monitor wells (CP-18A, CP-21A, and RF-41) and three springs (Jackson Blue Spring, Baltzell Springs Group, and Sandbag Spring) in the lower Apalachicola–Chattahoochee–Flint River Basin. Estimated nitrate input to groundwater at land surface, based on previous studies of nitrogen fertilizer sales and atmospheric nitrate deposition data, were used in the advective transport models for the period 2002 to 2050. Nitrate concentrations in groundwater samples collected from the six sites during 1993 to 2007 and groundwater age tracer data were used to calibrate the transport aspect of the simulations. Measured nitrate concentrations (as nitrogen) in wells and springs sampled during the study ranged from 0.37 to 12.73 milligrams per liter. Average apparent ages of groundwater calculated from measurements of chlorofluorocarbon, sulfur hexafluoride, and tritium from wells CP-18A, CP-21A,and RF-41 were about 23, 29, and 32 years, respectively. Average apparent ages of groundwater from Baltzell Springs Group, Sandbag Spring, and Jackson Blue Spring were about 16, 18, and 19 years, respectively. Simulated travel times of particles from the six selected sites ranged from less than 1 day to 511 years; both the minimum and maximum particle travel times were estimated for water from Jackson Blue Spring. Median simulated travel times of particles were about 30, 38, and 62 years for Jackson Blue Spring, Sandbag Spring, and Baltzell Springs Group, respectively. Study results indicated that travel times for approximately 50 percent of the particles from all spring sites were less than 50 years. The median simulated travel times of particles arriving at receptor wells CP-18A, CP-21A, and RF-41 were about 50, 35, and 36 years, respectively. All particle travel times were within the same order of magnitude as the tracer-derived average apparent ages for water, although slightly older than the measured ages. Travel time estimates were substantially greater than the measured age for groundwater reaching well CP-18A, as confirmed by the average apparent age of water determined from tracers. Local-scale particle-tracking models were used to predict nitrate concentrations in the three monitor wells and three springs from 2002 to 2050 for three nitrogen management scenarios: (1) fixed input of nitrate at the 2001 level, (2) reduction of nitrate inputs of 4 percent per year (from the previous year) from 2002 to 2050, and (3) elimination of nitrate input after 2001. Simulated nitrate concentrations in well CP-21A peaked at 7.82 milligrams per liter in 2030, and concentrations in background well RF-41 peaked at 1.10 milligrams per liter in 2020. The simulated particle travel times were longer than indicated by age dating analysis for groundwater in well CP-18A; to account for the poor calibration fit at this well, nitrate concentrations were shifted 21 years. With the shift, simulated nitrate concentrations in groundwater at CP-18A peaked at 13.76 milligrams per liter in 2026. For groundwater in Baltzell Springs Group, Jackson Blue Spring, and Sandbag Spring, simulated nitrate concentrations peaked at 3.77 milligrams per liter in 2006, 3.51 milligrams per liter in 2011, and 0.81 milligram per liter in 2018, respectively, under the three management scenarios. In management scenario 3 (elimination of nitrate input after 2001), simulated nitrate concentrations in Baltzell Springs Group decreased to less than background concentrations (0.10 milligram per liter) by 2033, and in Sandbag Spring concentrations decreased to less than background by 2041. Simulations using nitrate management scenarios 1 (fixed input of nitrate at 2001 levels) and 2 (reduction of 4.0 percent per year) indicate that nitrate concentrations in groundwater may remain above background concentrations through 2050 at all sites.

Alabama, Florida, Georgia↗

USGS2013

Groundwater and surface-water interaction within the upper Smith River Watershed, Montana 2006-2010

The 125-mile long Smith River, a tributary of the Missouri River, is highly valued as an agricultural resource and for its many recreational uses. During a drought starting in about 1999, streamflow was insufficient to meet all of the irrigation demands, much less maintain streamflow needed for boating and viable fish habitat. In 2006, the U.S. Geological Survey, in cooperation with the Meagher County Conservation District, initiated a multi-year hydrologic investigation of the Smith River watershed. This investigation was designed to increase understanding of the water resources of the upper Smith River watershed and develop a detailed description of groundwater and surface-water interactions. A combination of methods, including miscellaneous and continuous groundwater-level, stream-stage, water-temperature, and streamflow monitoring was used to assess the hydrologic system and the spatial and temporal variability of groundwater and surface-water interactions. Collectively, data are in agreement and show: (1) the hydraulic connectedness of groundwater and surface water, (2) the presence of both losing and gaining stream reaches, (3) dynamic changes in direction and magnitude of water flow between the stream and groundwater with time, (4) the effects of local flood irrigation on groundwater levels and gradients in the watershed, and (5) evidence and timing of irrigation return flows to area streams. Groundwater flow within the alluvium and older (Tertiary) basin-fill sediments generally followed land-surface topography from the uplands to the axis of alluvial valleys of the Smith River and its tributaries. Groundwater levels were typically highest in the monitoring wells located within and adjacent to streams in late spring or early summer, likely affected by recharge from snowmelt and local precipitation, leakage from losing streams and canals, and recharge from local flood irrigation. The effects of flood irrigation resulted in increased hydraulic gradients (increased groundwater levels relative to stream stage) or even reversed gradient direction at several monitoring sites coincident with the onset of nearby flood irrigation. Groundwater-level declines in mid-summer were due to groundwater withdrawals and reduced recharge from decreased precipitation, increased evapotranspiration, and reduced leakage in some area streams during periods of low flow. Groundwater levels typically rebounded in late summer, a result of decreased evapotranspiration, decreased groundwater use for irrigation, increased flow in losing streams, and the onset of late-season flood irrigation at some sites. The effect of groundwater and surface-water interactions is most apparent along the North and South Forks of the Smith River where the magnitude of streamflow losses and gains can be greater than the magnitude of flow within the stream. Net gains consistently occurred over the lower 15 miles of the South Fork Smith River. A monitoring site near the mouth of the South Fork Smith River gained (flow from the groundwater to the stream) during all seasons, with head gradients towards the stream. Two upstream sites on the South Fork Smith River exhibited variable conditions that ranged from gaining during the spring, losing (flowing from the stream to the groundwater) during most of the summer as groundwater levels declined, and then approached or returned to gaining conditions in late summer. Parts of the South Fork Smith River became dry during periods of losing conditions, thus classifying this tributary as intermittent. The North Fork Smith River is highly managed at times through reservoir releases. The North Fork Smith River was perennial throughout the study period although irrigation diversions removed a large percentage of streamflow at times and losing conditions persisted along a lower reach. The lowermost reach of the North Fork Smith River near its mouth transitioned from a losing reach to a gaining reach throughout the study period. Groundwater and surface-water interactions occur downstream from the confluence of the North and South Fork Smith Rivers, but are less discernible compared to the overall magnitude of the main-stem streamflow. The Smith River was perennial throughout the study. Monitoring sites along the Smith River generally displayed small head gradients between the stream and the groundwater, while one site consistently showed strongly gaining conditions. Synoptic streamflow measurements during periods of limited irrigation diversion in 2007 and 2008 consistently showed gains over the upper 41.4 river miles of the main stem Smith River where net gains ranged from 13.0 to 28.9 cubic feet per second. Continuous streamflow data indicated net groundwater discharge and small-scale tributary inflow contributions of around 25 cubic feet per second along the upper 10-mile reach of the Smith River for most of the 2010 record. A period of intense irrigation withdrawal during the last two weeks in May was followed by a period (early June 2010 to mid-July 2010) with the largest net increase (an average of 71.1 cubic feet per second) in streamflow along this reach of the Smith River. This observation is likely due to increased groundwater discharge to the Smith River resulting from irrigation return flow. By late July, the apparent effects of return flows receded, and the net increase in streamflow returned to about 25 cubic feet per second. Two-dimensional heat and solute transport VS2DH models representing selected stream cross sections were used to constrain the hydraulic properties of the Quaternary alluvium and estimate temporal water-flux values through model boundaries. Hydraulic conductivity of the Quaternary alluvium of the modeled sections ranged from 3x10-6 to 4x10-5 feet per second. The models showed reasonable approximations of the streambed and shallow aquifer environment, and the dynamic changes in water flux between the stream and the groundwater through different model boundaries.

Montana↗

USGS2018

Monitoring the water-quality response of agricultural conservation practices in the Bucks Branch watershed, Sussex County, Delaware, 2014–16

The purpose of this study was to evaluate the effects of irrigation and cover crops as conservation practices on water quality in groundwater and streams. Bucks Branch, a stream in the Nanticoke River watershed in southwestern Delaware, was identified as having one of the highest concentrations of nitrate in all surface-water sites sampled by the Delaware Department of Natural Resources and Environmental Control (DNREC). The study site is on two adjacent fields bordering Bucks Branch, one that has used irrigation since 2000 and one with dryland farming; both under conservation tillage and long-term rotation of corn, soybean, and small grain crops. A streamgage was installed near the study site fields to measure streamflow and water quality. The study area is typical of farming practices and environmental conditions throughout much of the intensively farmed agricultural land of the Coastal Plain of Delaware and surrounding parts of Maryland. Monitoring was conducted from January 2014 through June 2016. Corn was grown on both fields during the two growing seasons of the study period, and cover crops were planted before or shortly after harvest on both fields. During the second year of data collection, the effects of radish and rye grass cover crops on nutrient transport were studied. The combined results from data collected for this study show that water and nitrate moved below the root zone year round when soil moisture was high, especially after significant rainfall and frequently after irrigation. Soil water sampled 2 to 3 weeks after nutrients were applied had nitrate concentrations greater than 50 milligrams per liter as nitrogen (mg/L as N) and may be a significant source of nitrate to groundwater. Whereas recharge containing elevated nitrate concentrations also occurred under the dryland field, it was less frequent and of lower concentration than recharge under the irrigated field. Nitrate was present in all groundwater samples from these sites. Groundwater estimated to have recharged within 10 years or less had higher median concentrations of nitrate than in older water samples. The oldest groundwater encountered was over 30 years old, and had traveled along the longest, deep flowpaths from upland fields to the stream. The median nitrate concentration was 18 mg/L as N in younger water (less than 10 years old) beneath the irrigated field, compared to about 10 mg/L as N in younger water beneath the dryland field. Samples from the shallow upland wells in both study fields showed little, if any, evidence of denitrification. Several samples from deeper wells and from wells near forested riparian zone wetlands that border both fields did show partial denitrification. A mixing model estimated that between 12 and 22 percent of the nitrate discharging to the stream was lost through uptake and denitrification upstream of the streamgage on Bucks Branch. Continuous data collected at this site and evidence of denitrification in the surface-water samples showed a greater potential for loss of nitrate during the warmer months than the colder months. This pattern was similar to that seen below the streamgage at the most downstream site in the watershed. A mixed cover crop of radishes and rye was planted prior to removal (radishes) and just after harvest (rye) of the corn crop on the irrigated field. Rye grass was planted shortly after crop harvest on the dryland field. Cover crop biomass samples collected while radishes were growing and after they were killed by freezing temperatures indicates that the early planted radish crop effectively scavenged available nitrogen from the soil. Whereas radish biomass initially held more nitrogen than rye, at 55 to 8 pounds per acre, respectively, leaching of inorganic nitrate following radish die-off was minimal. Soil-water nitrate concentrations during the cover-crop growing period were lower than during the growing season prior to planting of the cover crop. There also was an increase in soil fertility and dissolved organic nitrate in samples of soil water that was likely related to increased soil microbial metabolism. Results indicate that cover crops stored plant nutrients over the winter and did not increase shallow groundwater concentrations of nitrate. Although conservation practices such as cover crops and nutrient management have been applied to these fields, there was still significant leaching of nitrate to groundwater, especially under the irrigated field. This will likely continue to be a challenge in this area and other parts of the Coastal Plain where soil moisture capacity is relatively low and managing irrigation around rainfall is difficult. Cover crops, when planted in standing corn, are one practice that can effectively pull nitrate from below the root zone to the top layer of soil, thus limiting the amount of potential nitrate leaching to groundwater. Irrigation management that would lower average soil moisture conditions during the growing season also could potentially limit nitrogen transport.

Delaware↗

USGS2013

Phosphorus losses from agricultural watersheds in the Mississippi Delta

Phosphorus (P) loss from agricultural fields is of environmental concern because of its potential impact on water quality in streams and lakes. The Mississippi Delta has long been known for its fish productivity and recreational value, but high levels of P in fresh water can lead to algal blooms that have many detrimental effects on natural ecosystems. Algal blooms interfere with recreational and aesthetic water use. However, few studies have evaluated P losses from agricultural watersheds in the Mississippi Delta. To better understand the processes influencing P loss, rainfall, surface runoff, sediment, ortho-P (orthophosphate, PO 4 –P), and total P (TP) were measured (water years 1996–2000) for two subwatersheds (UL1 and UL2) of the Deep Hollow Lake Watershed and one subwatershed of the Beasley Lake Watershed (BL3) primarily in cotton production in the Mississippi Delta. Ortho-P concentrations ranged from 0.01 to 1.0 mg/L with a mean of 0.17 mg/L at UL1 (17.0 ha), 0.36 mg/L at UL2 (11.2 ha) and 0.12 mg/L at BL3 (7.2 ha). The TP concentrations ranged from 0.14 to 7.9 mg/L with a mean of 0.96 mg/L at UL1, 1.1 mg/L at UL2 and 1.29 mg/L at BL3. Among the three sites, UL1 and UL2 received P application in October 1998, and BL3 received P applications in the spring of 1998 and 1999. At UL1, ortho-P concentrations were 0.36, 0.25 and 0.16 for the first, second and third rainfall events after P application, respectively; At UL2, ortho-P concentrations were 1.0, 0.66 and 0.65 for the first, second and third rainfall events after P application, respectively; and at BL3, ortho-P concentrations were 0.11, 0.22 and 0.09 for the first, second and third rainfall events after P application, respectively. P fertilizer application did influence P losses, but high P concentrations observed in surface runoff were not always a direct result of P fertilizer application or high rainfall. Application of P in the fall (UL1 and UL2) resulted in more ortho-P losses, likely because high rainfall often occurred in the winter months soon after application. The mean ortho-P concentrations were higher at UL1 and UL2 than those at BL3, although BL3 received more P application during the monitoring period, because P was applied in spring at BL3. However, tillage associated with planting and incorporating applied P in the spring (BL3) may have resulted in more TP loss in sediment, thus the mean TP concentration was the highest at BL3. Ortho-P loss was correlated with surface runoff; and TP loss was correlated with sediment loss. These results indicate that applying P fertilizer in the spring may be recommended to reduce potential ortho-P loss during the fallow winter season; in addition, conservation practices may reduce potential TP loss associated with soil loss.

Mississippi↗

USGS2014

Simulation of the effects of rainfall and groundwater use on historical lake water levels, groundwater levels, and spring flows in central Florida

The urbanization of central Florida has progressed substantially in recent decades, and the total population in Lake, Orange, Osceola, Polk, and Seminole Counties more than quadrupled from 1960 to 2010. The Floridan aquifer system is the primary source of water for potable, industrial, and agricultural purposes in central Florida. Despite increases in groundwater withdrawals to meet the demand of population growth, recharge derived by infiltration of rainfall in the well-drained karst terrain of central Florida is the largest component of the long-term water balance of the Floridan aquifer system. To complement existing physics-based groundwater flow models, artificial neural networks and other data-mining techniques were used to simulate historical lake water level, groundwater level, and spring flow at sites throughout the area. Historical data were examined using descriptive statistics, cluster analysis, and other exploratory analysis techniques to assess their suitability for more intensive data-mining analysis. Linear trend analyses of meteorological data collected by the National Oceanic and Atmospheric Administration at 21 sites indicate 67 percent of sites exhibited upward trends in air temperature over at least a 45-year period of record, whereas 76 percent exhibited downward trends in rainfall over at least a 95-year period of record. Likewise, linear trend analyses of hydrologic response data, which have varied periods of record ranging in length from 10 to 79 years, indicate that water levels in lakes (307 sites) were about evenly split between upward and downward trends, whereas water levels in 69 percent of wells (out of 455 sites) and flows in 68 percent of springs (out of 19 sites) exhibited downward trends. Total groundwater use in the study area increased from about 250 million gallons per day (Mgal/d) in 1958 to about 590 Mgal/d in 1980 and remained relatively stable from 1981 to 2008, with a minimum of 559 Mgal/d in 1994 and a maximum of 773 Mgal/d in 2000. The change in groundwater-use trend in the early 1980s and the following period of relatively slight trend is attributable to the concomitant effects of increasing public-supply withdrawals and decreasing use of water by the phosphate industry and agriculture. On the basis of available historical data and exploratory analyses, empirical lake water-level, groundwater-level, and spring-flow models were developed for 22 lakes, 23 wells, and 6 springs. Input time series consisting of various frequencies and frequency-band components of daily rainfall (1942 to 2008) and monthly total groundwater use (1957 to 2008) resulted in hybrid signal-decomposition artificial neural network models. The final models explained much of the variability in observed hydrologic data, with 43 of the 51 sites having coefficients of determination exceeding 0.6, and the models matched the magnitude of the observed data reasonably well, such that models for 32 of the 51 sites had root-mean-square errors less than 10 percent of the measured range of the data. The Central Florida Artificial Neural Network Decision Support System was developed to integrate historical databases and the 102 site-specific artificial neural network models, model controls, and model output into a spreadsheet application with a graphical user interface that allows the user to simulate scenarios of interest. Overall, the data-mining analyses indicate that the Floridan aquifer system in central Florida is a highly conductive, dynamic, open system that is strongly influenced by external forcing. The most important external forcing appears to be rainfall, which explains much of the multiyear cyclic variability and long-term downward trends observed in lake water levels, groundwater levels, and spring flows. For most sites, groundwater use explains less of the observed variability in water levels and flows than rainfall. Relative groundwater-use impacts are greater during droughts, however, and long-term trends in water levels and flows were identified that are consistent with historical groundwater-use patterns. The sensitivity of the hydrologic system to rainfall is expected, owing to the well-drained karst terrain and relatively thin confinement of the Floridan aquifer system in much of central Florida. These characteristics facilitate the relatively rapid transmission of infiltrating water from rainfall to the water table and contribute to downward leakage of water to the Floridan aquifer system. The areally distributed nature of rainfall, as opposed to the site-specific nature of groundwater use, and the generally high transmissivity and low storativity properties of the semiconfined Floridan aquifer system contribute to the prevalence of water-level and flow patterns that mimic rainfall patterns. In general, the data-mining analyses demonstrate that the hydrologic system in central Florida is affected by groundwater use differently during wet periods, when little or no system storage is available (high water levels), compared to dry periods, when there is excess system storage (low water levels). Thus, by driving the overall behavior of the system, rainfall indirectly influences the degree to which groundwater use will effect persistent trends in water levels and flows, with groundwater-use impacts more prevalent during periods of low water levels and spring flows caused by low rainfall and less prevalent during periods of high water levels and spring flows caused by high rainfall. Differences in the magnitudes of rainfall and groundwater use during wet and dry periods also are important determinants of hydrologic response. An important implication of the data-mining analyses is that rainfall variability at subannual to multidecadal timescales must be considered in combination with groundwater use to provide robust system-response predictions that enhance sustainable resource management in an open karst aquifer system. The data-driven approach was limited, however, by the confounding effects of correlation between rainfall and groundwater use, the quality and completeness of the historical databases, and the spatial variations in groundwater use. The data-mining analyses indicate that available historical data when used alone do not contain sufficient information to definitively quantify the related individual effects of rainfall and groundwater use on hydrologic response. The knowledge gained from data-driven modeling and the results from physics-based modeling, when compared and used in combination, can yield a more comprehensive assessment and a more robust understanding of the hydrologic system than either of the approaches used separately.

Florida↗

USGS2012

Sediment transport to and from small impoundments in northeast Kansas, March 2009 through September 2011

The U.S. Geological Survey, in cooperation with the Kansas Water Office, investigated sediment transport to and from three small impoundments (average surface area of 0.1 to 0.8 square miles) in northeast Kansas during March 2009 through September 2011. Streamgages and continuous turbidity sensors were operated upstream and downstream from Atchison County, Banner Creek, and Centralia Lakes to study the effect of varied watershed characteristics and agricultural practices on sediment transport in small watersheds in northeast Kansas. Atchison County Lake is located in a predominantly agricultural basin of row crops, with wide riparian buffers along streams, a substantial amount of tile drainage, and numerous small impoundments (less than 0.05 square miles; hereafter referred to as “ponds”). Banner Creek Lake is a predominantly grassland basin with numerous small ponds located in the watershed, and wide riparian buffers along streams. Centralia Lake is a predominantly agricultural basin of row crops with few ponds, few riparian buffers along streams, and minimal tile drainage. Upstream from Atchison County, Banner Creek, and Centralia Lakes 24, 38, and 32 percent, respectively, of the total load was transported during less than 0.1 percent (approximately 0.9 days) of the time. Despite less streamflow in 2011, larger sediment loads during that year indicate that not all storm events transport the same amount of sediment; larger, extreme storms during the spring may transport much larger sediment loads in small Kansas watersheds. Annual sediment yields were 360, 400, and 970 tons per square mile per year at Atchison County, Banner, and Centralia Lake watersheds, respectively, which were less than estimated yields for this area of Kansas (between 2,000 and 5,000 tons per square mile per year). Although Centralia and Atchison County Lakes had similar percentages of agricultural land use, mean annual sediment yields upstream from Centralia Lake were about 2.7 times those at Atchison County or Banner Creek Lakes. These data indicate larger yields of sediment from watersheds with row crops and those with fewer small ponds, and smaller yields in watersheds which are primarily grassland, or agricultural with substantial tile drainage and riparian buffers along streams. These results also indicated that a cultivated watershed can produce yields similar to those observed under the assumed reference (or natural) condition. Selected small ponds were studied in the Atchison County Lake watershed to characterize the role of small ponds in sediment trapping. Studied ponds trapped about 8 percent of the sediment upstream from the sediment-sampling site. When these results were extrapolated to the other ponds in the watershed, differences in the extent of these ponds was not the primary factor affecting differences in yields among the three watersheds. However, the selected small ponds were both 45 years old at the time of this study, and have reduced capacity because of being filled in with sediments. Additionally, trapping efficiency of these small ponds decreased over five observed storms, indicating that processes that suspended or resuspended sediments in these shallow ponds, such as wind and waves, affected their trapping efficiencies. While small ponds trapped sediments in small storms, they could be a source of sediment in larger or more closely spaced storm events. Channel slope was similar at all three watersheds, 0.40, 0.46, and 0.31 percent at Atchison County, Banner Creek, and Centralia Lake watersheds, respectively. Other factors, such as increased bank and stream erosion, differences in tile drainage, extent of grassland, or riparian buffers, could be the predominant factors affecting sediment yields from these basins. These results show that reference-like sediment yields may be observed in heavily agricultural watersheds through a combination of field-scale management activities and stream channel protection. When computing loads using published erosion rates obtained by single-point survey methodology, streambank contributions from the main stem of Banner Creek are three times more than the sediment load observed by this study at the sediment sampling site at Banner Creek, 2.6 times more than the sediment load observed by this study at the sediment sampling site at Clear Creek (upstream from Atchison County Lake), and are 22 percent of the load observed by this study at the sediment sampling site at Black Vermillion River above Centralia Lake. Comparisons of study sites to similarly sized urban and urbanizing watersheds in Johnson County, Kansas indicated that sediment yields from the Centralia Lake watershed were similar to those in construction-affected watersheds, while much smaller sediment yields in the Atchison County and Banner Creek watersheds were comparable to stable, heavily urbanized watersheds. Comparisons of study sites to larger watersheds upstream from Tuttle Creek Lake indicate the Black Vermillion River watershed continues to have high sediment yields despite 98 percent of sediment from the Centralia watershed (a headwater of the Black Vermillion River) being trapped in Centralia Lake. Estimated trapping efficiencies for the larger watershed lakes indicated that Banner Creek and Centralia Lakes trapped 98 percent of incoming sediment, whereas Atchison County Lake trapped 72 percent of incoming sediment during the 3-year study period.

Kansas↗

USGS2021

Global cropland-extent product at 30-m resolution (GCEP30) derived from Landsat satellite time-series data for the year 2015 using multiple machine-learning algorithms on Google Earth Engine cloud

Executive Summary Global food and water security analysis and management require precise and accurate global cropland-extent maps. Existing maps have limitations, in that they are (1) mapped using coarse-resolution remote-sensing data, resulting in the lack of precise mapping location of croplands and their accuracies; (2) derived by collecting and collating national statistical data that are often subjective, leading to substantial uncertainties in cropland-area estimates, as well as their locations; and (3) extracted from one or more classes of a land use–land cover product in which cropland classes are not the focus of mapping, leading to their mixing with other classes and creating significant errors of omission and commission. These limitations can be overcome by producing high-resolution cropland-extent maps using satellite-sensor data, such as Landsat 30-m resolution or higher. The most fundamental cropland product is the high-resolution cropland-extent map because all higher level cropland products, such as crop-watering method (that is, whether crops are irrigated or rainfed), crop types, cropping intensities, cropland fallows, crop productivity, and crop-water productivity, are dependent on a precise and accurate cropland-extent product. Given these realities, the overarching goal of this study was to produce a Landsat satellite-derived global cropland-extent product at 30-m resolution. The work, which involved a paradigm shift in how global cropland-extent maps are produced, involved the following five key steps: (1) petabyte-scale computing that involved multiyear, 8- to 16-day, time-series Landsat 30-m resolution data for the global land surface; (2) composition of analysis-ready data (ARD) cubes; (3) creation of a large global-reference data hub for machine learning; (4) use of multiple machine-learning algorithms (MLAs) by writing software and computing in the cloud; and (5) Google Earth Engine (GEE) cloud computing. The five key steps involved nine distinct phases. First, the world was segmented into 74 agroecological zones (AEZs). Second, Landsat 8- to 16-day data were used to time-composite 10-band (blue, green, red, near-infrared, short-wave infrared band 1, short-wave infrared band 2, thermal infrared, enhanced vegetation index, normalized difference water index, and normalized difference vegetation index) Landsat 30-m resolution data cubes for every 2- to 4-month time period during 3- to 4-year periods (stated as nominal-year 2015 or, simply, 2015), along with two additional 30-m resolution bands (Shuttle Radar Topography Mission elevation, and slope) in each of the 74 AEZs. Third, more than 100,000 reference-training data samples were collected using ground data (some of which were collected using a mobile application), as well as submeter- to 5-m-resolution, very high-resolution imagery sourced from other reliable sources. Fourth, reference-training data were used to create a knowledge base for separating cropland from noncropland. Fifth, MLAs such as the pixel-based supervised random forest and support-vector machines were written on the GEE using Python and JavaScript. Sixth, object-based recursive hierarchical segmentation algorithm was used, in addition to MLAs, to overcome uncertainties. Seventh, MLAs used the knowledge base to classify and separate cropland from noncropland. Eighth, accuracy assessment was conducted by generating error matrices for each of the 74 AEZs using 19,171 independent validation-data samples. Ninth, cropland areas were computed for all countries of the world and compared with United Nation’s (UN’s) Food and Agricultural Organization (FAO) and other national statistics. The outcome was a Landsat-derived global cropland-extent product at 30-m resolution (GCEP30), which has an overall accuracy of 91.7 percent. For the cropland class, producer’s accuracy was 83.4 percent, and user’s accuracy was 78.3 percent. GCEP30 calculated (using direct pixel count) the global net-cropland area (GNCA) for the year 2015 as 1.873 billion hectares (~12.6 percent of the Earth’s terrestrial area). The continental cropland distribution as a percentage of GNCA was Asia, 33 percent; Europe, 25.5 percent; Africa, 16.7 percent; North America, 14.4 percent; South America, 8.1 percent; and Australia and Oceania, 2.4 percent. The worldwide cropland areas in GCEP30 for 2015 were higher by 236 to 299 million hectares (Mha) compared to national statistics reported elsewhere for the same year (for example, in Food and Agriculture Organization’s corporate statistical database [FAOSTAT] and in the monthly irrigated and rainfed crop areas [MIRCA] database). The global cropland area reported for 2015 increased by 344 Mha (22.5 percent), compared to the year 2000. During the same period (2000–2015), the world’s population increased by 20 percent. Whereas some of these areal increases are real increases in cropland areas, others are due to the types of data, methods, and approaches used. Using the highest known resolution (compared to previous coarse-resolution global products) enabled this study to capture fragmented croplands. Coarse-resolution data compute areas on the basis of subpixels, which, for a large proportion of certain land use–land cover classes, will show only a certain percentage of the total pixel area as actual area. Subpixel areas can lead to substantial uncertainties in area computation, as determining the exact fraction of cropland areas within a coarse-resolution pixel is resource intensive and subject to errors. Other innovations in GCEP30 include reference-data hubs, machine learning, and cloud computing. Cropland areas in 214 countries, territories, departments, and regions were calculated for the year 2015 using GCEP30, on the basis of UN’s global administrative unit layers (GAUL) boundaries. The 10 leading countries in terms of cropland area (as a percentage of the GNCA) were India (9.6 percent), United States (8.95 percent), China (8.82 percent), Russia (8.32 percent), Brazil (3.42 percent), Ukraine (2.32 percent), Canada (2.29 percent), Argentina (2.05 percent), Indonesia (2 percent), and Nigeria (1.91 percent). Together, these 10 countries occupy 50 percent of the global cropland, and they have 52 percent of the global population. Their combined cropland area increased by 2 percent between 2000 and 2015, compared to the substantial increase in population of 517 million (15.5 percent). Together, India, United States, China, and Russia encompass 36 percent of the total area. In the United States and Canada, from 2000 to 2015, cropland decreased by about 2 percent, whereas their populations increased by 14 and 13 percent, respectively. The additional food requirements in these 10 countries, which are caused by increased populations, as well as increasing nutritional demands, are met by production increases in existing cropland or through virtual food trade, or both. More than 18 countries, territories, departments, or regions had 60 percent or more of their geographic area as cropland: Republic of Moldova, San Marino, and Hungary had more than 80 percent of the country’s area as cropland; Denmark, Ukraine, Ireland, and Bangladesh, 70 to 80 percent; and Uruguay, Netherlands, United Kingdom, Spain, Lithuania, Poland, Gaza Strip, Czechia, Italy, India, and Azerbaijan, 60 to 70 percent. Europe and South Asia can be considered agricultural capitals of the world, on the basis of their percentages of geographic area as cropland. United States, China, and Russia, which all have high cropland areas, are ranked second, third, and fourth in the world; India is ranked first. However, the amount of cropland as a percentage of the country’s geographic area is relatively very low for United States (18.3 percent), China (17.7 percent), and Russia (9.5 percent), whereas it is 60.5 percent for India. Most African and South American countries, territories, departments, or regions have less than 15 percent of their geographic area as cropland. China and India together house 36 percent of the world’s population; however, between 2000 and 2015, the amount of China’s cropland area fell by 18.9 percent, owing to urban expansion and the abandonment of farmlands caused by demographic changes (that is, the movement of population from villages to cities). In contrast, China’s population grew by 10 percent. The amount of India’s cropland increased by 8.5 percent, whereas its population grew by 20 percent. This study showed that, out of the 10 leading cropland countries, Ukraine, Nigeria, Russia, and Indonesia showed an 18 to 31 percent increase in cropland areas, on the basis of GCEP30 by the year 2015, compared to 2000. Nigeria’s cropland area increased by 25 percent, and its population increased by 31 percent in the same period. In these countries, food security is maintained by cropland expansion, productivity increases, and virtual food trade. Nevertheless, this trend of increasing net-cropland area and productivity will likely become difficult to maintain, owing to diminishing arable lands and plateauing of 50 years of continual yield increases, requiring policymakers to explore novel and data-supported approaches to solving future food security issues. The GCEP30 product, which can be browsed at full resolution at www.croplands.org , has been released for public download and use through U.S. Geological Survey (USGS)–National Aeronautics and Space Administration (NASA) Land Processes Distributed Active Archive Center (see https://lpdaac.usgs.gov/news/release-of-gfsad-30-meter-cropland-extent-products/ ).

Professional Paper↗

USGS1970

Hydrology of the Upper Malad River basin, southeastern Idaho

The report area comprises 485 square miles in the Basin and Range physiographic province. It includes most of eastern' Oneida County and parts of Franklin, Bannock, and Power Counties of southeastern Idaho. Relief is about 5,000 feet; the floor of the Malad Valley is at an average altitude of about 4,400 feet. Agriculture is, by far, ,the principal economic .activity. In 1960 the population of the upper Malad River basin was about 3,600, of which about 60 percent resided in Malad City, the county seat of Oneida County. The climate is semiarid throughout the Malad Valley and its principal tributary valleys; ,above 6,500 feet the climate is subhumid. Annual precipitation ranges from about 13 inches in the lower Malad Valley to more than 30 inches on the highest peaks of the Bannock and Malad ranges. Owing to ,the normally clear atmospheric conditions, large daily and seasonal temperature fluctuations are common. Topography, distance from the Pacific Ocean, .and the general atmospheric circulation are the principal factors governing the climate of the Malad River basin. The westerlies transport moisture from the P.acific Ocean toward southeastern Idaho. The north-south tren4ing mountains flanking the basin are oriented orthogonally to the moisture flux so that they are very effective in removing precipitable water from the air. A minimum uplift of 6,000 feet is required to transport moisture from the Pacific source region; accordingly, most air masses are desiccated long before they reach the Malad basin. Heaviest precipitation is generally associated with steep pressure gradients in the midtroposphere that are so oriented as to cause a deep landward penetration of moisture from the Pacific Ocean. Annual water yields in the project area range from about 0.8 inch in the, lower Malad Valley to more than 19 inches on the high peaks north and east of Malad City. The mean annual water yield for the entire basin is 4 inches, or about 115,000 acre-feet. Evaporation is greatest in July when about 7 inches is lost from lakes, reservoirs, and waterlogged areas; losses from free-water surfaces may be as much .as 38 inches annually. An extensive ground-water reservoir consisting of sand and gravel interbedded with relatively impermeable beds of silt .and clay underlies much of the Malad Valley. Wells near the center of the valley exceeding 700 feet in depth do not reach bedrock. The Woodruff fault, which transects the constricted lower Malad Valley, is one of the main factors creating artesian conditions south of the latitude of Malad City. Recharge is obtained principally from mountain runoff which flows onto highly permeable alluvial fans surrounding the valley and from streams that flow across the valley floor. On the basis of a water balance analysis, under flow from the project area was estimated to be 28,000 acre-feet annually, surface-water outflow was 51,000 acre-feet, and transbasin imports were about 4,000 acre-feet. The principal tributaries of the Malad River are perennial along their upper and middle reaches and have well-sustained low flows. During the growing season, all surface water entering the Malad Valley is used for irrigation. Spine irrigation is practiced in the principal tributary valleys; however, a shortage of suitable reservoir sites has hampered surface-water development in these areas. The highly porous deposits underlying the Malad Valley tend to attenuate flood peaks. An unusual combination of meteorologic events early in 1962 effectively counteracted the high absorptive capacity of the valley and predisposed the basin to high flood risk. Subsequent rapid snowmelt combined with frozen ground produced the extraordinary flood of February 12, 1962. Calcium and bicarbonate commonly are the most abundant ions in the surface waters of the upper Malad River basin. In August 1967, the dissolved-solids content of streamflow ranged from 200 to 350 milligrams per liter in the middle and upper parts of the basin; however, much greater values were measured in the Malad River between Woddruff and Cherry Creek Lane. With the exception of that reach, the surface water of the project area is suitable for irrigating all but the most sensitive crops. The total water yield is not sufficient to meet all the water needs of the basin. A comprehensive water-management plan is required to ensure optimal use of the water resource.

Idaho↗

USGS2021

Groundwater quality and age of secondary bedrock aquifers in the glaciated portion of eastern Nebraska, 2016–18

The Eastern Nebraska Water Resources Assessment (ENWRA) project was initiated in 2006 to assist water managers by developing a hydrogeologic framework and water budget for the glaciated portion of eastern Nebraska. Within the ENWRA area, the primary groundwater sources for municipal, domestic, and irrigation water needs are provided by withdrawals from alluvial, buried paleovalley, and the High Plains aquifer (where present). Generally, other bedrock aquifers are considered a secondary water source. However, in some areas, such as parts of Sarpy and Nemaha Counties, these secondary bedrock aquifers are the only source of water within glaciated upland areas. To improve the understanding of the quality, geochemistry, and age of groundwater from bedrock aquifers, the U.S. Geological Survey (USGS), in cooperation with the ENWRA group, which includes the Lewis and Clark, Lower Elkhorn, Lower Platte North, Lower Platte South, Nemaha, and Papio-Missouri River Natural Resources Districts, designed a study to sample 31 wells completed in the secondary bedrock aquifers and analyze samples for major ions, physical properties, nutrients, stable isotopes, and selected age tracers. Of the 31 samples collected for this report, 22 samples were collected from the Dakota aquifer contained in the Dakota Sandstone, 3 from the Niobrara aquifer contained in the Niobrara Formation of Colorado Group, and 6 from Paleozoic aquifers contained in undifferentiated Paleozoic-age units. The results of this study indicate that major ion data collected from the Dakota aquifer can be used for assessing the quality, recharge source, and age of groundwater. Calcium bicarbonate dominant samples were characterized as modern or mixed, indicating that, in these areas, groundwater is unconfined and is recharged by precipitation and (or) surface water. If groundwater extraction rates exceed recharge rates, total dissolved solid concentrations may increase as a result of upwelling of groundwater from deeper units or formations, which can adversely affect groundwater quality. Sampling results presented in this report indicate water quality is good, but that groundwater in the Dakota aquifer with calcium bicarbonate water type may be vulnerable to surface contamination. In contrast, groundwater sampled from the Dakota aquifer, having a dominant water type other than calcium bicarbonate, generally has low dissolved oxygen and nitrate concentrations, and higher concentrations of total dissolved solids and trace elements, including iron and strontium. The geochemical characteristics of noncalcium bicarbonate samples from the Dakota aquifer indicated confining conditions and limited groundwater recharge from local precipitation. Apparent groundwater ages estimated from radiocarbon (carbon-14) sampling of noncalcium bicarbonate samples from the Dakota aquifer indicated that the time of groundwater recharge to the Dakota aquifer occurred during Pleistocene time. Depleted stable isotopes results indicate recharge during a colder climate. Groundwater under confined conditions is not easily recharged from precipitation or surface water. Future groundwater-level monitoring in locations where the Dakota aquifer appears to be confined could provide information to evaluate whether groundwater supplies remain sufficient to meet future municipal, domestic, and irrigation needs. For the Niobrara aquifer and Paleozoic aquifers, the dominant water type was not a diagnostic indicator of recharge source, age, and groundwater quality as with the Dakota aquifer. Most likely this is because the host formation was dominated by calcium-carbonate-rich rocks; however, few samples were collected from these aquifers to be able to confirm this interpretation. Samples collected from wells completed in the Niobrara aquifer and Paleozoic aquifers and characterized as calcium sulfate water type have statistically significantly higher concentrations of total dissolved solids compared to other samples from the Niobrara aquifer and Paleozoic aquifers characterized as calcium bicarbonate. Given that six of the nine of samples collected from the Niobrara and Paleozoic aquifers indicated modern recharge, these secondary bedrock aquifers are reliant on precipitation to sustain groundwater levels and may be vulnerable to a multiyear drought. Well yields of the Niobrara and Paleozoic aquifers are dependent on the presence of secondary porosity and these units offer little storage. Samples collected from wells completed in Paleozoic aquifers were the most isotopically enriched and similar to modern precipitation and had the highest concentrations of nitrate, indicating that groundwater is affected by agricultural activities. Future groundwater sampling would be beneficial to characterize groundwater-quality changes within the Niobrara and Paleozoic aquifers over time.

Nebraska↗

USGS2016

Sediment load and distribution in the lower Skagit River, Skagit County, Washington

The Skagit River delivers about 40 percent of all fluvial sediment that enters Puget Sound, influencing flood hazards in the Skagit lowlands, critically important estuarine habitat in the delta, and some of the most diverse and productive agriculture in western Washington. A total of 175 measurements of suspended-sediment load, made routinely from 1974 to 1993, and sporadically from 2006 to 2009, were used to develop and evaluate regression models of sediment transport (also known as “sediment-rating curves”) for estimating suspended-sediment load as a function of river discharge. Using a flow-range model and 75 years of daily discharge record (acquired from 1941 to 2015), the mean annual suspended-sediment load for the Skagit River near Mount Vernon, Washington, was estimated to be 2.5 teragrams (Tg, where 1 Tg = 1 million metric tons). The seasonal model indicates that 74 percent of the total annual suspended‑sediment load is delivered to Puget Sound during the winter storm season (from October through March), but also indicates that discharge is a poor surrogate for suspended‑sediment concentration (SSC) during the summer low-flow season. Sediment-rating curves developed for different time periods revealed that the regression model slope of the SSC-discharge relation increased 66 percent between the periods of 1974–76 and 2006–09 when suspended-sediment samples were collected, implying that changes in sediment supply, channel hydraulics, and (or) basin hydrology occurred between the two time intervals. In the relatively wet water year 2007 (October 1, 2006, through September 30, 2007), an automated sampler was used to collect daily samples of suspended sediment from which an annual load of 4.5 Tg was calculated, dominated by a single large flood event that contributed 1.8 Tg, or 40 percent of the total. In comparison, the annual load calculated for water year 2007 using the preferred flow-range model was 4.8 Tg (+6.7 percent), in close agreement with the measured value. Particle size affects sediment transport, fate and distribution across watersheds, and therefore is important for predicting how coastal environments, particularly deltas and beaches, will respond to changes in climate and sea-level. Particle-size analysis of winter storm samples indicated that about one-half of the suspended-sediment load consisted of fines (that is, silt- and clay-sized particles smaller than 0.0625 mm in diameter), and the remainder consisted of mostly fine- to medium-sized sand (0.0625–0.5 mm), whereas bedload during winter storm flows (about 1–3 percent of total sediment load) was predominantly composed of medium to coarse sand (0.25–1 mm). A continuous turbidity record from the Anacortes Water Treatment Plant (water years 1999–2013), used as a surrogate for the concentration of fines ( R 2 = 0.93, p = 4.2E-10, n = 17), confirms that about one-half of the mean annual suspended-sediment load is composed of fines. The distribution of flow through the delta distributaries (that is, the channels into which the main stem splits as it approaches the delta) is dynamic, with twice as much flow through the North Fork of the Skagit River relative to the South Fork during low-flow conditions, and close to equal flows in the two channels during high-flow conditions. Turbidity, monitored at several locations in the lower river in spring 2009, was essentially uniform among sites, indicating that fines are well mixed in the lower Skagit River system (defined as the Skagit River and all its distributaries downstream of the Mount Vernon streamgage). A strong relation ( R 2 = 0.95, p = 3.2E-14, n = 21; linear regression) between the concentration of fines and turbidity measured at various locations in summer 2009 indicates that turbidity is an effective surrogate for the concentration of fines, independent of location in the river, under naturally well-mixed fluvial conditions. This relation is especially useful for monitoring suspended sediment in western Washington rivers that are seasonally dominated by glacier meltwater because glacial melting typically produces suspended-sediment concentrations that are not well correlated with discharge. These results provide a comprehensive set of tools to estimate sediment delivery and delta responses of interest to scientists and resource managers including decision-makers examining options for flood hazard mitigation, estuary restoration, and climate change adaptation.

Washington↗

USGS2015

Water-quality trends in the Scituate reservoir drainage area, Rhode Island, 1983-2012

The Scituate Reservoir is the primary source of drinking water for more than 60 percent of the population of Rhode Island. Water-quality and streamflow data collected at 37 surface-water monitoring stations in the Scituate Reservoir drainage area, Rhode Island, from October 2001 through September 2012, water years (WYs) 2002-12, were analyzed to determine water-quality conditions and constituent loads in the drainage area. Trends in water quality, including physical properties and concentrations of constituents, were investigated for the same period and for a longer period from October 1982 through September 2012 (WYs 1983-2012). Water samples were collected and analyzed by the Providence Water Supply Board, the agency that manages the Scituate Reservoir. Streamflow data were collected by the U.S. Geological Survey. Median values and other summary statistics for pH, color, turbidity, alkalinity, chloride, nitrite, nitrate, total coliform bacteria, Escherichia coli ( E. coli ), and orthophosphate were calculated for WYs 2003-12 for all 37 monitoring stations. Instantaneous loads and yields (loads per unit area) of total coliform bacteria and E. coli , chloride, nitrite, nitrate, and orthophosphate were calculated for all sampling dates during WYs 2003-12 for 23 monitoring stations with streamflow data. Values of physical properties and concentrations of constituents were compared with State and Federal water-quality standards and guidelines and were related to streamflow, land-use characteristics, varying classes of timber operations, and impervious surface areas. Tributaries in the Scituate Reservoir drainage area for WYs 2003-12 were slightly acidic (median pH of all stations equal to 6.1) and contained low median concentrations of chloride (22 milligrams per liter [mg/L]), nitrate (0.01 mg/L as nitrogen), nitrite (0.001 mg/L as nitrogen), and orthophosphate (0.02 milligrams per liter as phosphorus [mg/L as P]). Turbidity and alkalinity values also were low with medians of 0.57 nephelometric turbidity units and 5.1 mg/L as calcium carbonate, respectively. Total coliform bacteria and E. coli were detected in most samples from all stations, but median concentrations were generally low-43 colony-forming units per 100 milliliters (mL) and 15 colony-forming units per 100 milliliters, respectively. Median values of several physical properties and median concentrations of several constituents correlated positively with the percentages of developed land and negatively with the percentages of forest cover in the drainage areas above the monitoring stations. Median concentrations of chloride correlated positively with the percentages of impervious land use in the subbasins of monitoring stations, likely reflecting the effects of deicing compounds applied to roadways during winter maintenance. Median concentrations of alkalinity also correlated positively with the percentage of impervious land use, which may be related to the deterioration of fabricated structures containing calcium carbonate. Median values of color correlated positively with the percentage of wetland area in the subbasins of monitoring stations, reflecting the natural sources of color in tributaries. Streamflows were negatively correlated with turbidity and concentrations of total coliform bacteria and E. coli, possibly reflecting seasonal patterns in which relatively high values of these properties and constituents occur during warmer low-flow conditions late in the water year. Similar seasonal patterns were observed for pH, alkalinity, and color. Negative correlations between concentrations of chloride and streamflow also were significant, indicating that deicing salts from roadways and other impervious surfaces that lack direct connection to the tributaries are likely infiltrating to the groundwater and discharging to some of the tributaries late in the water year. While salt-laden runoff directly enters some of the tributaries at roadway crossings, most of the roadway runoff infiltrates into the adjacent berms throughout the drainage area. Statistically significant correlations were not identified between various degrees of tree-canopy reduction caused by timber operations in the subbasins and median values or concentrations of water-quality properties. Loads and yields of chloride, nitrate, nitrite, orthophosphate, and bacteria varied at monitoring stations in the Scituate Reservoir drainage area in WYs 2003-12. Loads generally were greater at stations in the Barden Reservoir and the Regulating Reservoir Subbasins that have larger drainage areas than in subbasins with smaller drainage areas. Subbasin yields of fecal-indicator bacteria and orthophosphate generally were largest in the Westconnaug Reservoir Subbasin, and subbasin yields for chloride, nitrate, and nitrite were largest in the Moswansicut Reservoir Subbasin in the northeastern part of the drainage area. Upward trends in pH were identified for nearly half of the monitoring stations for WYs 1983-2012 and may reflect regional reductions in acid precipitation. Many upward trends in alkalinity also were identified for both the WYs 1983-2012 and for WYs 2003-12 periods and are likely related to the natural weathering of structures containing concrete or, in some cases, the application of lime or fertilizers on agriculture lands. Significant trends in chloride concentrations at most stations during WYs 1983-2012 were upward; however, results for WYs 2003-12 substantiate few significant upward trends and, in a few cases, downward trends were identified in several tributary drainage areas.

Rhode Island↗

USGS1967

Swatara Creek basin of southeastern Pennsylvania: An evaluation of its hydrologic system

Local concentrations of population in the Swatara Creek basin of Pennsylvania find it necessary to store, transport, and treat water because local supplies are either deficient or have been contaminated by disposal of wastes in upstream areas. Water in the basin is available for the deficient areas and for dilution of the coal-mine drainage in the northern parts and the sewage wastes in the southern parts. Swatara Creek drains 576 square miles just east of Harrisburg, Pa., and is the largest tributary to the Susquehanna River from the north side below Harrisburg. It rises in the southern Pocono Mountains and flows southwestward across the Lebanon Plateau. On an average day Swatara Creek discharges more than 630 million gallons into the Susquehanna River at Middletown, Pa. In a year this amounts to about 23 inches of water over the entire basin and is the residual from an average annual precipitation of 45.5 inches. During an average year the flow in Swatara Creek from the upper third of the basin above Harper Tavern is always greater than 1,300 mgd (million gallons per day) for at least 15 days and is always greater than 25 mgd for at least 350 days. The daily streamflow from the basin averages 1.1 mgd per sq mi, but yields from different areas range from 0.97 to 1.22 mgd per sq mi. These variations are caused chiefly by differences in precipitation and land cover. The area of lowest yield is in the valleys west of Tremont, and the highest yields are in the Upper and Lower Little Swatara Creek subbasins. At high and medium stages the chemical character of the water in the streams is suitable for public and private supplies. At lower stages, defending on the areas and the amounts of contamination by coal-mine drainage and sewage pollution, the natural flow may require some treatment. At low stages the chemical characteristics of the natural flow not affected by man is almost identical with that of the ground water in the area drained by the stream. In general, the total dissolved solids range from about 25 to 400 parts per million and the hardness is as much as about 300 parts per million. The ground-water increment to the base flow of Swatara Creek averages about 240 mgd, or about 8.8 inches annually, for the basin. Generally, ground-water supplies in amounts of less than 0.5 mgd can be developed south of Blue Mountain. Supplies of several million gallons per day have been developed for industrial use from the permeable limestones in the south-central part of the basin. More intensive investigation in other parts of the basin would indicate areas where supplies of more than 0.5 mgd could be developed from properly spaced wells. The chemical character of water from wells depends largely on the host rock. In highly soluble rocks water contains large amount of dissolved solids; in more resistant rocks concentrations are lower. The chemical character of unpolluted ground water generally reflects the composition of the more readily soluble minerals in the local geologic environment. Areas contaminated by septic- tank effluent may have above normal amounts of nitrate and detergent products. Except where polluted, most ground water is suitable for public and industrial uses without extensive treatment. Sites for storage of surface water exist in the part of the basin lying in the valley and ridge area. As much as 30 to 40 percent of the annual flow could be impounded for release as low-flow augmentation for dilution of mine drainage and other wastes in the basin. Low sediment yields of supplying drainage areas would ensure a long life expectancy of reservoirs at these sites. Overbank flooding of the main stem of the Swatara Creek and its tributaries has occurred many times in the past. However, it has not been a hazard because urban development has not encroached on the flood plain. An inundation map of the August 1933 flood provides a basis that urban planners may use to avoid future damage. As water in the Swatara Creek moves downstream to the Susquehanna River, the flow is influenced consecutively by a large annual rainfall on the northern valley and ridge area, the wastes of surface and subsurface coal-mining activities, and less annual rainfall on the part of the basin lying in the Lebanon Plateau area; the flow is supplemented and further influenced by many tributaries and by the industrial and domestic wastes that are carried by these secondary streams. The annual precipitation ranges from 52 inches at the east edge and 49 inches at the west edge of the mountainous part of the basin to about 41 inches at the southwestern part at Middletown. The rainfall generally is adequate during the growing season to mature the crops. The mean annual temperature at Lebanon is about 52°F, and the growing season is about 180 days. In this report the basin has been divided into eight hydrologic zones, leased on runoff, natural use of water, and chemical character of water. Four zones lie in the valley and ridge area, three lie in the Lebanon Plateau area, and one lies in the highland along the southeastern basin boundary. In each of the zones the hydrologic characteristics are virtually the same, but they may be completely different from those in adjacent zones. The boundaries of the zones generally coincide with boundaries between geologic formations, and the areas in each zone include rocks of similar influence on water. Streams in zone 4 at the northeast edge of the plateau have the highest average surface runoff from 1.2 to 1.1 mgd per sq mi whereas those in zone 2 at the northwest edge of the valley and ridge area have the lowest, about 1.0 mgd. Streams in zone 8, along the southeast edge of the basin, have the largest sustained low-flow yield, about 0.26 to 0.19 mgd per sq mi; those in zone 5 overlying the Martinsburg Shale east of Harrisburg have the smallest sustained low-flow yields, 0.03 to 0.01 mgd. Streams in the limestone area of zone 7 have the greatest range in low-flow yields in any one zone from 0.60 to 0 mgd per sq mi. Low-flow yields in zones 1 through 4 range from 0.13 to 0.03 mgd per sq mi. Surface flows from zones 1 and 2 are generally acidic and contain high concentrations of sulfate, iron, and total dissolved solids especially where contaminated with mine wastes. Surface flows from zones 3 and 4 are dilute, slightly alkaline, and suitable for public water supplies. Surface flows from zones 5, 6, and 7 are alkaline and contain moderate concentrations of dissolved solids with waters of highest hardness occurring in zone 7. Surface flows from zone 8 are dilute to moderately mineralized and are relatively high in silica concentration. Nitrate concentrations are high in surf Fee flows below sewage outfalls and in ground water contaminated by septic tank effluent and industrial wastes. Average annual sediment yields of 550 to 650 tons per square mile are characteristic of zones 1 and 2 where strip mining has destroyed the forest cover and coal culm is carried into the streams. From agricultural lands on the Martinsburg Shale in zones 5 and 6, annual sediment yields range from 300 to 350 tons per square mile; but from agricultural lands on the siliceous rocks in zone 8 and zones 3 and 4 in the valley and ridge area, the sediment yield ranges from 200 to 250 tons annually per square mile. Lowest annual sediment yields in the basin are in the forested areas of siliceous rocks in zones 2, 3, 4, and 5, and in the sinkhole topography of the limestones in zone 7 where the yield ranges from 30 to 35 tons and 50 to 60 tons per square mile, respectively. The amount of ground water that can be developed in the basin is dependent on the ability of the underlying rocks to yield water to wells. More than 300 gpm (gallons per minute) can be obtained from wells in alluvial materials in the valley bottoms and in some of the limestones where large solution channels and fractures are penetrated by the wells. From 50 to 300 gpm can be obtained from wells in loosely cemented sandstones and in fractured limestones. From 10 to 50 gpm can be developed from wells in the shales and harder sandstones. The most dense rocks will yield from 1 to 10 gpm from fractures and crevices. Most wells yield water from the upper 350 feet of the formation, for this part contains the most fractures or solution channels. Studies show that the velocity at which a contaminant will move downstream in the basin is related to the discharge of the stream at the time. At a stream discharge of about 400 mgd at Pine Grove, a contaminant in Swatara Creek would require about 40 hours to move from Pine Grove to Middletown. As a result of dispersion and dilution, the maximum concentration of the contaminant at Middletown would be less than 20 percent the concentration at Pine Grove under these conditions. An evaluation of the availability of water in the basin indicates that about I,239 mgd enters as precipitation, 630 mgd leaves as streamflow, 580 mgd is evaporated and transpired, and 56 mgd is diverted for use by man. Not all the diversions for man's use are lost to the basin, as about 27 mgd is returned as sewage for reuse. About one-fourth of the waste water is returned to the ground and the remainder to stream drainageways. Of that diverted by man, 11.6 mgd is used for public supply and 44.4 mgd for industrial and private supplies. Diversions of streamflow furnish 86 percent of the public supply and 27 percent of the industrial supply, and ground-water sources yield the remainder. Municipal and private sewage treatment plants are upgrading the waste water in many places, but no provisions are being made for treatment other than natural dilution and assimilation for the 15 mgd of coal-nine drainage in the northern part of the basin. Technology for economic treatment of mine water is not available at this time, although research in this field is being done. Urbanization eastward from Harrisburg and around Lebanon has increased the population density of the basin. Densities of 500 people per square mile and water use exceeding 2.0 mgd per sq mi can be expected in the future. By the year 2000 the population of the basin may increase 60 percent; and if the per capita rate of use increases 0.5 percent per year the domestic requirements for water will be about two times the present use, or 23 mgd. Similarly, if the present 1:4 ratio of domestic use to industrial use of water continues, at least 89 mgd will be needed for industry in the future. Although an increase to twice the present use of water can be foreseen, or 112 mgd, water for the dilution and assimilation of wastes from treatment systems are not included. Providing water for dilution of wastes from treatment plants has not been a problem, but in the future the amounts needed for this purpose will be greater as the population increases. As water becomes more valuable, treatment of sewage wastes to reduce the biochemical-oxygen-demand load by at least 80 to 90 percent will be necessary to conserve water for more productive uses. As much as 100 mgd may be needed for waste dilution in the basin by year 2000. The present trends in suburban and light industrial development will probably persist in the basin. Problems arising through changes in economic value of water, conflicts in use, and alternatives in development are typical of those confronting the manager of a water-resource system.

Pennsylvania↗

USGS2004

Hydrology and cycling of nitrogen and phosphorus in Little Bean Marsh: A remnant riparian wetland along the Missouri River in Platte County, Missouri, 1996–97

The lack of concurrent water-quality and hydrologic data on riparian wetlands in the Midwestern United States has resulted in a lack of knowledge about the water-quality functions that these wetlands provide. Therefore, Little Bean Marsh, a remnant riparian wetland along the Missouri River, was investigated in 1996 and 1997 primarily to determine the magnitude and character of selected water-quality benefits that can be produced in such a wetland and to identify critical processes that can be managed in remnant or restored riparian wetlands for amelioration of water quality. Little Bean Marsh averages 69 hectares in size, has a maximum depth of about 1 meter, and the majority of the marsh is covered by macrophytes. In 1997, 41 percent of the water received by Little Bean Marsh was from direct precipitation, 14 percent was from ground-water seepage, 30 percent from watershed runoff, and 15 percent was backflow from Bean Lake. Although, Little Bean Marsh was both a ground-water recharge and discharge area, discharge to the marsh was three times the recharge to ground water. Ground-water levels closely tracked marsh water levels indicating a strong hydraulic connection between ground water and the marsh. Reduced surface runoff and ground-water availability are stabilizing influences on marsh hydrology and probably contribute to the persistence of emergent vegetation. The rapid hydraulic connection between Little Bean Marsh and ground water indicates that the hydrologic regime of most wetlands along the lower Missouri River is largely a function of the altitude of the marsh bottom relative to the altitude of the water table. More water was lost from the marsh through evapotranspiration (59 percent) than all other pathways combined. This is partially because the transpiration process of abundant macrophytes can greatly contribute to the evapotranspiration above that lost from open water surfaces. Surface outflow accounted for 36 percent and ground-water seepage accounted for only 5 percent of the losses. Large residence times allows the marsh to greatly affect water quality before water escapes as ground-water recharge or surface outflow. The shallowness of Little Bean Marsh and ion exclusion during ice formation caused the highest specific conductances of 1,100 to 1,300 microsiemens per centimeter at 25 degrees Celsius to occur during the winter. This concentration of dissolved solutes under ice can make wetlands more vulnerable to toxic contaminants than deeper surface-water bodies. Dissolved oxygen was less than 5 mg/L (milligrams per liter) for 3 to 4 months and near 0 mg/L for about 1 month in summer. Despite depths of less than 1 meter, temperature stratification persisted more than 3 months during the summers of 1996 and 1997, preventing mixing and contributing to periods of anoxia. Shallow depths and extended periods of anoxia in the marsh limit the ability of some organisms to escape high-temperature stress. Turbidity in Little Bean Marsh usually was low for several reasons: sediment loadings from the largely flood-plain drainage were low, emergent vegetation shade out algae and shield the water from wind, and high concentrations of bivalent cations increase flocculation rates of inorganic suspended material. The high concentrations of bivalent cations was largely because of a substantial amount of ground-water seepage into the marsh. Dissolved organic nitrogen was the dominant nitrogen species in Little Bean Marsh. Denitrification and biotic uptake kept more than 62 percent of nitrate (NO3) and 43 percent of ammonium (NH4) concentrations in marsh samples less than a detection limit of 0.005 mg/L. This contrasts with the Missouri River where inorganic NO3 dominates. Consequently, artificial flood-plain drainage that bypasses riparian wetlands likely deliver substantially more biotically available inorganic nitrogen to receiving waters than surface water that has been routed through a remnant wetland. Average total nitrogen concentrations in Little Bean Marsh were substantially less than those at other Missouri River wetlands, roughly one-half the mean concentrations in the Missouri River, but roughly twice the average nitrogen values in reservoirs of the glaciated plains of Missouri. The largest concentrations of nearly all species of nitrogen and phosphorus and the most intense period of hypereutrophy coincided with a phytoplankton bloom and senescence of River Bulrush (Scirpus fluviatilis) and common cattail (Typha latifolia) in September 1997. The rapid leaching of nitrogen that occurs soon after macrophyte senescence combined with a recent destratification of the marsh probably provided nitrogen to the nitrogen-limited open-water areas and triggered a phytoplankton bloom. Despite the rarity of runoff events, surface runoff from the watershed, combined with atmospheric deposition, contributed more than seven times the 530 kg (kilograms) of nitrogen that escaped Little Bean Marsh in surface outflow during 1997. Atmospheric deposition alone was more than 530 kg. Seepage to ground water contained less than 1.5 percent of the nitrogen leaving the marsh in surface outflow. The slow decay rate of Scirpus fluviatilis and reducing conditions in bottom sediments make burial of organic nitrogen a substantial sink of nitrogen. Denitrification experiments indicate that denitrification rates were limited by NO3 in the water column. Consequently, decomposition and nitrification of NH4 and organic nitrogen are the rate limiting steps of nitrogen removal in Little Bean Marsh. The NO3-limited rates of denitrification also indicate that Little Bean Marsh has a large unused capacity for nitrogen removal. These data indicate that the vast extent of riparian marshes along the Missouri and Mississippi Rivers may have had a substantial role in limiting NO3 loads to the Gulf of Mexico before agricultural development of flood plains. Drainage and removal of riparian marshes may be a major cause of the increased NO3 loads to the Gulf of Mexico. Periods of anoxia had much larger effects on phosphorus release than the other variables. The largest concentrations of phosphorus occurred in late summer and corresponded with senescing macrophytes, periods of anoxia, and a large algal bloom in Little Bean Marsh. Low water levels prevented the escape of phosphorus in surface outflow during these periods of highest phosphorus concentrations. Dry weather in late summer is typical and probably makes the correspondence of low water levels, anoxia, and consequent low phosphorus release a common occurrence in marshes along the Missouri River. Little Bean Marsh retained more than 95 percent of the phosphorus it received. The amount of phosphorus in surface inflows to the marsh were more than one order of magnitude greater than that escaping in surface outflows. The long hydraulic residence time of the marsh and large contributions of iron from ground water (that provide many sorption sites for phosphorus) make the marsh an effective sediment and phosphorus trap.

Missouri↗

USGS2002

Water quality and environmental isotopic analyses of ground-water samples collected from the Wasatch and Fort Union Formations in areas of coalbed methane development — Implications to recharge and ground-water flow, eastern Powder River Basin, Wyoming

Chemical analyses of ground-water samples were evaluated as part of an investigation of lower Tertiary aquifers in the eastern Powder River Basin where coalbed methane is being developed. Ground-water samples were collected from two springs discharging from clinker, eight monitoring wells completed in the Wasatch aquifer, and 13 monitoring or coalbed methane production wells completed in coalbed aquifers. The ground-water samples were analyzed for major ions and environmental isotopes (tritium and stable isotopes of hydrogen and oxygen) to characterize the composition of waters in these aquifers, to relate these characteristics to geochemical processes, and to evaluate recharge and ground-water flow within and between these aquifers. This investigation was conducted in cooperation with the Wyoming State Engineer's Office and the Bureau of Land Management. Water quality in the different aquifers was characterized by major-ion composition. Samples collected from the two springs were classified as calcium-sulfate-type and calcium-bicarbonate-type waters. All ground-water samples from the coalbed aquifers were sodium-bicarbonate-type waters as were five of eight samples collected from the overlying Wasatch aquifer. Potential areal patterns in ionic composition were examined. Ground-water samples collected during this and another investigation suggest that dissolved-solids concentrations in the coalbed aquifers may be lower south of the Belle Fourche River (generally less than 600 milligrams per liter). As ground water in coalbed aquifers flows to the north and northwest away from an inferred source of recharge (clinker in the study area), dissolved-solids concentrations appear to increase. Variation in ionic composition in the vertical dimension was examined qualitatively and statistically within and between aquifers. A relationship between ionic composition and well depth was noted and corroborates similar observations by earlier investigators in the Powder River Basin in both Wyoming and Montana. This relationship results in two different water-quality zones with different characteristics - a shallow zone, comprising the upper part of the Wasatch aquifer, characterized by mixed cation composition and either sulfate or bicarbonate as the dominant anion; and a deeper zone, comprising the lower (deeper) part of the Wasatch aquifer and the underlying coalbed aquifers, characterized by sodium-bicarbonate-type waters. The zonation appears to be related to geochemical processes described by earlier investigators such as dissolution and precipitation of minerals, ion exchange, sulfate reduction, and mixing of waters. Qualitative and statistically significant differences were observed in sulfate concentrations between the coalbed aquifers and the overlying Wasatch aquifer. Ionic composition suggests that bacterially mediated redox processes such as sulfate reduction were probably the dominant geochemical processes in the anaerobic coalbed aquifers. Tritium was used to qualitatively estimate the time of ground-water recharge. Tritium concentrations in both springs suggests that both were recharged after 1952 and contain modern water. Tritium was not detected at concentrations suggestive of modern water in any ground-water samples collected from the coalbed aquifers or in six of eight ground-water samples collected from the overlying Wasatch aquifer. Tritium concentrations in the remaining two wells from the Wasatch aquifer suggest a mixture between submodern (recharged before 1952) and modern water, although the low concentrations suggest that ground water in these two wells have very little modern water. The relative absence of modern water in all aquifers in the study area suggests that recharge processes to these aquifers are probably very slow. Paired d 2 H (deuterium/hydrogen isotopic ratio) and d 18 O (oxygen-18/oxygen-16 isotopic ratio) values for samples collected from the springs and all aquifers are close to the Global Meteoric Water Line, a meteoric water line for North American continental precipitation, and an estimated local meteoric water line, suggesting the water in the aquifers is of meteoric origin. The d 2 H and d 18 O values suggest that the waters were recharged in a colder climate or temperature, mid-latitudes, and mid-continent. In general, the samples do not form separate groups based on aquifer origin; this suggests either intermixing of the waters in the aquifers or that the different aquifers are subject to similar recharge and/or evolutional paths for the water. However, examination of the differences in the values of d 2 H and d 18 O, in combination with major-ion chemistry at three monitoring-well clusters, suggest that changes in the values with depth may represent different timing or sources of recharge to the different aquifers. The areal distribution of d 2 H was examined and an apparent break in the d 2 H along a northwest to southeast trend was observed. In the coalbed aquifers, all but one ground-water sample (collected from the Big George coal bed), show a pattern where the d 2 H values become more negative towards the center of the Powder River Basin and values greater (less negative) than an arbitrary reference value of -140 ‰ (per mil or parts per thousand) were observed near the outcrop area of the Wyodak-Anderson coal zone. In the overlying Wasatch aquifer, the d 2 H values became less negative towards the center of the basin. The values more negative than -140 ‰ are near the outcrop area and the values that are less negative than -140 ‰ are closer to the basin center. It is unclear if this pattern is a result of sample size, different recharge mechanisms, geochemical processes, or if the processes producing these differences are independent or unrecognized. Results of water-quality sampling were compared with selected regulatory and non-regulatory standards as well as commonly-used guidelines for proposed water uses. Dissolved solids was the measure that most frequently exceeded U.S. Environmental Protection Agency public water-supply standards and State of Wyoming domestic-use standards in ground-water samples collected from the Wasatch aquifer and coalbed aquifers. The State of Wyoming agricultural standards (irrigation) for sulfate and dissolved solids were exceeded in some samples collected from the Wasatch aquifer and coalbed aquifers. The State of Wyoming livestock standard for pH was exceeded in some samples collected from the Wasatch aquifer. Water from the Wasatch aquifer ranged from soft to very hard, and water from the coalbed aquifers ranged from moderately hard to very hard. Samples collected from wells completed in both the Wasatch aquifer and coalbed aquifers plotted in a wide range of both sodium- and salinity-hazard classes, but most samples clustered in or near the combined medium-sodium-hazard—high-salinity-hazard classes.

Wyoming↗

USGS2013

Simulation of groundwater flow, effects of artificial recharge, and storage volume changes in the Equus Beds aquifer near the city of Wichita, Kansas well field, 1935–2008

The Equus Beds aquifer is a primary water-supply source for Wichita, Kansas and the surrounding area because of shallow depth to water, large saturated thickness, and generally good water quality. Substantial water-level declines in the Equus Beds aquifer have resulted from pumping groundwater for agricultural and municipal needs, as well as periodic drought conditions. In March 2006, the city of Wichita began construction of the Equus Beds Aquifer Storage and Recovery project to store and later recover groundwater, and to form a hydraulic barrier to the known chloride-brine plume near Burrton, Kansas. In October 2009, the U.S. Geological Survey, in cooperation with the city of Wichita, began a study to determine groundwater flow in the area of the Wichita well field, and chloride transport from the Arkansas River and Burrton oilfield to the Wichita well field. Groundwater flow was simulated for the Equus Beds aquifer using the three-dimensional finite-difference groundwater-flow model MODFLOW-2000. The model simulates steady-state and transient conditions. The groundwater-flow model was calibrated by adjusting model input data and model geometry until model results matched field observations within an acceptable level of accuracy. The root mean square (RMS) error for water-level observations for the steady-state calibration simulation is 9.82 feet. The ratio of the RMS error to the total head loss in the model area is 0.049 and the mean error for water-level observations is 3.86 feet. The difference between flow into the model and flow out of the model across all model boundaries is -0.08 percent of total flow for the steady-state calibration. The RMS error for water-level observations for the transient calibration simulation is 2.48 feet, the ratio of the RMS error to the total head loss in the model area is 0.0124, and the mean error for water-level observations is 0.03 feet. The RMS error calculated for observed and simulated base flow gains or losses for the Arkansas River for the transient simulation is 7,916,564 cubic feet per day (91.6 cubic feet per second) and the RMS error divided by (/) the total range in streamflow (7,916,564/37,461,669 cubic feet per day) is 22 percent. The RMS error calculated for observed and simulated streamflow gains or losses for the Little Arkansas River for the transient simulation is 5,610,089 cubic feet per day(64.9 cubic feet per second) and the RMS error divided by the total range in streamflow (5,612,918/41,791,091 cubic feet per day) is 13 percent. The mean error between observed and simulated base flow gains or losses was 29,999 cubic feet per day (0.34 cubic feet per second) for the Arkansas River and -1,369,250 cubic feet per day (-15.8 cubic feet per second) for the Little Arkansas River. Cumulative streamflow gain and loss observations are similar to the cumulative simulated equivalents. Average percent mass balance difference for individual stress periods ranged from -0.46 to 0.51 percent. The cumulative mass balance for the transient calibration was 0.01 percent. Composite scaled sensitivities indicate the simulations are most sensitive to parameters with a large areal distribution. For the steady-state calibration, these parameters include recharge, hydraulic conductivity, and vertical conductance. For the transient simulation, these parameters include evapotranspiration, recharge, and hydraulic conductivity. The ability of the calibrated model to account for the additional groundwater recharged to the Equus Beds aquifer as part of the Aquifer Storage and Recovery project was assessed by using the U.S. Geological Survey subregional water budget program ZONEBUDGET and comparing those results to metered recharge for 2007 and 2008 and previous estimates of artificial recharge. The change in storage between simulations is the volume of water that estimates the recharge credit for the aquifer storage and recovery system. The estimated increase in storage of 1,607 acre-ft in the basin storage area compared to metered recharge of 1,796 acre-ft indicates some loss of metered recharge. Increased storage outside of the basin storage area of 183 acre-ft accounts for all but 6 acre-ft or 0.33 percent of the total. Previously estimated recharge credits for 2007 and 2008 are 1,018 and 600 acre-ft, respectively, and a total estimated recharge credit of 1,618 acre-ft. Storage changes calculated for this study are 4.42 percent less for 2007 and 5.67 percent more for 2008 than previous estimates. Total storage change for 2007 and 2008 is 0.68 percent less than previous estimates. The small difference between the increase in storage from artificial recharge estimated with the groundwater-flow model and metered recharge indicates the groundwater model correctly accounts for the additional water recharged to the Equus Beds aquifer as part of the Aquifer Storage and Recovery project. Small percent differences between inflows and outflows for all stress periods and all index cells in the basin storage area, improved calibration compared to the previous model, and a reasonable match between simulated and measured long-term base flow indicates the groundwater model accurately simulates groundwater flow in the study area. The change in groundwater level through recent years compared to the August 1940 groundwater level map has been documented and used to assess the change of storage volume of the Equus Beds aquifer in and near the Wichita well field for three different areas. Two methods were used to estimate changes in storage from simulation results using simulated change in groundwater levels in layer 1 between stress periods, and using ZONEBUDGET to calculate the change in storage in the same way the effects of artificial recharge were estimated within the basin storage area. The three methods indicate similar trends although the magnitude of storage changes differ. Information about the change in storage in response to hydrologic stresses is important for managing groundwater resources in the study area. The comparison between the three methods indicates similar storage change trends are estimated and each could be used to determine relative increases or decreases in storage. Use of groundwater level changes that do not include storage changes that occur in confined or semi-confined parts of the aquifer will slightly underestimate storage changes; however, use of specific yield and groundwater level changes to estimate storage change in confined or semi-confined parts of the aquifer will overestimate storage changes. Using only changes in shallow groundwater levels would provide more accurate storage change estimates for the measured groundwater levels method. The value used for specific yield is also an important consideration when estimating storage. For the Equus Beds aquifer the reported specific yield ranges between 0.08 and 0.35 and the storage coefficient (for confined conditions) ranges between 0.0004 and 0.16. Considering the importance of the value of specific yield and storage coefficient to estimates of storage change over time, and the wide range and substantial overlap for the reported values for specific yield and storage coefficient in the study area, further information on the distribution of specific yield and storage coefficient within the Equus Beds aquifer in the study area would greatly enhance the accuracy of estimated storage changes using both simulated groundwater level, simulated groundwater budget, or measured groundwater level methods.

Kansas↗

USGS2013

Plankton communities and summertime declines in algal abundance associated with low dissolved oxygen in the Tualatin River, Oregon

Phytoplankton populations in the Tualatin River in northwestern Oregon are an important component of the dissolved oxygen (DO) budget of the river and are critical for maintaining DO levels in summer. During the low-flow summer period, sufficient nutrients and a long residence time typically combine with ample sunshine and warm water to fuel blooms of cryptophyte algae, diatoms, green and blue-green algae in the low-gradient, slow-moving reservoir reach of the lower river. Algae in the Tualatin River generally drift with the water rather than attach to the river bottom as a result of moderate water depths, slightly elevated turbidity caused by suspended colloidal material, and dominance of silty substrates. Growth of algae occurs as if on a “conveyor belt” of streamflow, a dynamic system that is continually refreshed with inflowing water. Transit through the system can take as long as 2 weeks during the summer low-flow period. Photosynthetic production of DO during algal blooms is important in offsetting oxygen consumption at the sediment-water interface caused by the decomposition of organic matter from primarily terrestrial sources, and the absence of photosynthesis can lead to low DO concentrations that can harm aquatic life. The periods with the lowest DO concentrations in recent years (since 2003) typically occur in August following a decline in algal abundance and activity, when DO concentrations often decrease to less than State standards for extended periods (nearly 80 days). Since 2003, algal populations have tended to be smaller and algal blooms have terminated earlier compared to conditions in the 1990s, leading to more frequent declines in DO to levels that do not meet State standards. This study was developed to document the current abundance and species composition of phytoplankton in the Tualatin River, identify the possible causes of the general decline in algae, and evaluate hypotheses to explain why algal blooms diminish in midsummer. Plankton and water-quality sample data from 2006 to 2008 were combined with parts of a larger discrete-sample and continuous water-quality monitoring dataset and examined to identify patterns in water-quality and algal conditions since 1991, with a particular emphasis on 2003–08. Longitudinal plankton surveys were conducted in 2006–08 at six sites between river miles (RM) 24.5 and 3.4 at 2- to 3-week intervals, or 5–6 per season, and in-situ bioassay experiments were conducted in 2008 to examine the potential effects of wastewater treatment facility (WWTF) effluent and phosphorus additions on phytoplankton biomass and algal photosynthesis. Phytoplankton and zooplankton community composition, streamflow, and water-quality data were analyzed using multivariate statistical techniques to gain insights into plankton dynamics to determine what factors might be most tied to the abundance and characteristics of the phytoplankton assemblages, and identify possible causes of their declines. The connection between low-DO events and algal declines was clearly evident, as bloom crashes were nearly always followed by periods of low DO. Algal blooms occurred each year during 2006–08, producing maximum chlorophyll-a (Chl-a) values in June or July generally in the range of 50–80 micrograms per liter (µg/L). Bloom crashes and absence of sufficient algal photosynthesis in mid- to late-summer contributed to minimum DO concentrations that were less than the State standard of 6.5 milligrams per liter (mg/L) based on the 30-day mean daily concentration, for 62–74 days each year. At times, the absolute minimum State standard (4 mg/L DO) also was not met. To learn more about why low-DO events occurred, specific algal declines during 2003–08 were scrutinized to determine their likely causal factors. From this information, a series of hypotheses were formulated and evaluated in terms of their ability to explain recent declines in algal populations in the river in late summer. Meteorological, streamflow, turbidity, water temperature, and conductance conditions in the Tualatin River during the 2006–08 summer seasons were not atypical. Natural flow comprised the majority (70–80 percent) of flow each year during spring, but then reduced to 38–40 percent during midsummer when WWTF effluent—which contributed as much as 36 percent—and flow augmentation releases comprised a greater fraction of the flow. Summer 2008 was unusual, however, in the prolonged influence from the Wapato Lake agricultural area near Gaston in the upper part of the basin. The previous winter flooding and levee breach at Wapato Lake caused a much greater area of inundation. As a result, drainage from this area continued into July, much later than normal. A subsequent algal bloom in Wapato Lake then seeded the upper Tualatin River, and this drainage had a profound effect on the downstream plankton community. A large blue-green algae bloom developed—the largest in recent memory—prompting a public health advisory for recreational contact for about two weeks. Algal growths and surface blooms are a common feature of the Tualatin River. Most of the dominant algae have growth forms and morphologies that are well suited for planktonic life, employing spines and gas vacuoles to resist settling, forming colonies, and producing mucilage (or toxins) to resist zooplankton grazing. In 2006–08, 143 algal taxa were identified in 117 main-stem samples; diatoms and green algae were more diverse than blue-green, golden, and cryptophyte algae, although these later groups sometimes dominated the overall volumetric abundance (biovolume). The most frequently occurring taxa, occurring in 97–99 percent of samples, were flagellated cryptophytes Cryptomonas erosa and Rhodomonas minuta. Other important algal taxa included centric diatoms Stephanodiscus, Cyclotella, and Melosira species and colonial green algae Scenedesmus and Actinastrum. These taxa comprised the majority of the algal biovolume during much of the growing season. A general seasonal trend in the phytoplankton assemblages was observed, with dominance by filamentous centric diatoms Stephanodiscus and Melosira in spring and early summer, and flagellated cryptophytes and green algae, particularly Chlamydomonas sp., in late-summer; or, in 2008, dominance by blue-green algae Anabaena flos-aquae and Aphanizomenon flos-aquae during the Wapato Lake bloom event. There were 99 zooplankton taxa identified from the Tualatin River in 2006–08, composed primarily of cladocerans, copepods, and rotifers. A seasonal increase in zooplankton abundance was observed in early summer just as or shortly after the phytoplankton population began to increase, with populations growing to 15,000−120,000 organisms per cubic meter in the lower river. Zooplankton abundance showed a predictable and distinct longitudinal downstream increase, particularly downstream of Highway 99W (RM 11.6). Although grazing rates were not measured, the data suggest that, at times, zooplankton grazing may affect algal abundance and species composition in the Tualatin River, with diatoms becoming relatively less abundant and flagellated cryptophytes and green algae relatively more abundant during periods when zooplankton densities were highest. Multivariate statistical analyses identified soluble reactive phosphorus (SRP), natural flow, flow augmentation, and WWTF effluent as important factors influencing Tualatin River phytoplankton populations, with zooplankton density (particularly rotifers and copepods), specific conductance, chloride, and water temperature also having an important influence. Although SRP was highly correlated with the plankton communities, that correlation was likely the result of high or low algal activity (uptake) as SRP concentrations were often reduced to low levels during blooms. While previous studies have already established that phosphorus, among other factors such as flow, places a theoretical cap on the size of the phytoplankton population in the river, sometimes algal declines occur when SRP concentrations are apparently sufficient. To identify alternative causal factors, additional analyses were performed without SRP to focus on other water-quality parameters, zooplankton density, and flow factors. Considering data for all 3 years and including just those samples from the lower Tualatin River not affected by the 2008 Wapato Lake drainage event, three factors (percentage of reservoir flow augmentation, total natural flow, and rotifer density) best explained variations in the phytoplankton assemblages. Analyses focusing on the possible causes of algal declines included the above multivariate analyses, scrutiny of 10 specific instances of declines in algal populations during 2003–08 including several bloom–crash sequences, and analyses of historic routine watershed monitoring data from Clean Water Services. Six factors were hypothesized to be important in causing bloom crashes or impeding blooms from rebounding in August: (1) light limitation from cloudy weather, (2) a reduction in the plankton inocula or “seed” entering the lower river from upstream sources, (3) increased summer streamflows, (4) changes in the dominant sources of flow as the percentage of flow augmentation and WWTF discharges have increased, (5) zooplankton grazing, and (6) low concentrations of bioavailable phosphorus (<0.015 milligram per liter). All of these hypotheses are supported in some fashion by the available data and statistical analyses. Zooplankton grazing, short-term declines in photosynthesis from cloudy weather, total flow as it affects residence time, and the dominant source of flow are primary factors responsible for the low-DO events caused by declines in algae in the lower Tualatin River during late summer. Cloudy weather and increased turbidity are known to inhibit algal growth in the Tualatin River, and slight increases in turbidity in recent years may be a problem. Upstream sources of algae are critical in determining the characteristics and size of downstream populations, as illustrated by the Wapato Lake bloom in 2008, but more data are needed from upstream to fully define the importance of this connection. The sources of flow, through their differential contribution of plankton inocula (quality and amount), were, at times, important factors affecting phytoplankton populations. While SRP concentrations were often most highly correlated with phytoplankton species community, the bioavailability of phosphorus is still somewhat unknown and there are several sources to consider. Preliminary bioassay tests suggested that while treated wastewater effluent may stimulate algae at 30 percent concentrations, negative effects (or decreased stimulation) on Chl-a and DO production may occur at concentrations of 50 percent. Targeted data collection and future research will be needed to further understand the importance of these factors on Tualatin River phytoplankton. While the data and analysis completed for this report provide insights into future research and monitoring that would be useful to continue, additional monitoring of turbidity, Chl-a, and plankton abundance and species composition in the upper part of the basin would enhance our understanding of plankton dynamics and factors affecting phytoplankton abundance in the lower river. Assessment of the key upstream sources of algal inocula via surveys of the major flow sources as well as tributaries and wetlands would provide useful information for the management of river water quality. Other studies that could prove useful for developing management strategies include targeted experiments to evaluate the bioavailability of phosphorus from a variety of sources. New research on phytoplankton–zooplankton interactions, and studies of planktivorous fish, might also provide insight about food web dynamics and potential “top-down” effects of fish predation on the plankton communities. In addition, further development of neural-network or other water-quality models would help to evaluate management strategies and provide forecasts of water-quality conditions. Finally, periodic future reassessments of the available data with the multivariate statistical tools used in this study would be helpful to assess whether and how plankton communities are changing, and to continue to shed light on the importance of factors shaping the plankton. Although certain types and sizes of algal blooms are undesirable, minimum phytoplankton populations are an important part of aquatic food webs and are needed to maintain healthy levels of DO in the river. By understanding the sources, characteristics, causal factors, and responses of the plankton communities, management strategies can be developed to improve DO conditions in the lower Tualatin River during the important summer low-flow period.

Oregon↗