SEARCH · USGS Science

Results for “Scientific Monthly”

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

SourceResearch connectionFrom yearThrough year

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

At least 775 records · Page 43Linked to original sources

USGS2024

Groundwater hydrology, groundwater and surface-water interactions, aquifer testing, and groundwater-flow simulations for the Fountain Creek alluvial aquifer, near Colorado Springs, Colorado, 2018–20

From 2018 through 2020, the U.S. Geological Survey, in cooperation with the Air Force Civil Engineering Center, conducted an integrated study of the Fountain Creek alluvial aquifer located near Colorado Springs, Colorado. The objective of the study was to characterize hydrologic conditions for the alluvial aquifer pertinent to the potential for transport of solutes. Specific goals of this report were to characterize the groundwater hydrology of the area, to quantify groundwater and surface-water interactions, to estimate hydraulic properties of the aquifer using aquifer testing, and to complete numerical simulations of groundwater flow. Synoptic groundwater-level elevation measurements completed throughout this study, and as part of other U.S. Geological Survey programs between 1994 and 2020, indicate groundwater-level elevations fluctuate on annual and interannual timeframes. Groundwater-level fluctuations likely were caused by temporally variable groundwater recharge and discharge components in the area, with many wells showing maximum groundwater-level elevations during the winter months (November through March). From an interannual perspective, groundwater-level fluctuations appear to have reached maximum values during 2000 to 2003, decreased during 2003 to 2006, and remained relatively constant since that time, with the exception of several wells which have displayed rising groundwater-level elevations since 2018. Spatial evaluation of groundwater-level elevations indicates groundwater flow is generally from northeast to southwest within the vicinity of several alluvial paleochannels occurring along the northeastern margin of the aquifer. Within the center of the aquifer along Fountain Creek, groundwater flow is generally from north to south, approximately paralleling surface-water flow. To quantitatively understand the potential effect of groundwater recharge and groundwater pumping on fluctuations in groundwater-level elevation, a statistical transfer-function-noise model was applied. Results of the statistical model indicate throughout most of the aquifer, fluctuations were primarily the result of recharge seasonality. In the main stem of the aquifer where groundwater pumping wells were more concentrated, however, groundwater-level elevation fluctuations were also attributable to groundwater pumping through time. Three-dimensional evaluation of the aquifer geometry near Fountain Creek was combined with synoptic streamflow measurement and accounting of stream gains and losses to evaluate groundwater and surface-water interactions in the study area. Streamflow gain or loss calculations indicate Fountain Creek both gains from and loses flow to the alluvial aquifer, and gaining or losing reaches of the stream may be partially controlled by the depth to bedrock near the stream. Reaches with streamflow gains tend to coincide with areas where the estimated depth to bedrock is decreasing, meaning the alluvial aquifer is likely thinning in these areas and groundwater-flow paths may be converging and discharging groundwater to the stream. Losing reaches tended to coincide with locally greater depth to bedrock where the alluvial aquifer is likely thicker and has greater storage potential for surface water lost from Fountain Creek. Results of aquifer testing indicate hydraulic conductivity, estimated from slug tests and single-well pumping tests, ranged from 0.32 to 1,410 feet per day (ft/d) and 4.13 to 664 ft/d, respectively. These results are similar to the range of values from previous aquifer tests in the study area. Hydraulic conductivities from aquifer testing for this study were generally greater than the estimates of previous slug tests and had a mean value less than the estimates from previous pumping tests. Spatial evaluation of aquifer testing results indicates hydraulic conductivity tends to be greater in the main stem of the alluvial aquifer and lower in paleochannels upgradient from the main stem of the aquifer. The spatial variation in hydraulic conductivity may be attributed to the geomorphologic processes that formed the alluvial aquifer. Compacted sediment in the paleochannels has not been potentially transported sufficient distance to cause grain-size sorting, resulting in a poorly sorted deposit and lower hydraulic conductivities. In the central portion of the alluvial aquifer, near Fountain Creek, the sediments have been transported farther from their source areas and are likely better sorted, removing finer grained sediments that would cause lower hydraulic conductivity. A numerical groundwater-flow model was calibrated for the Fountain Creek alluvial aquifer for 2000–19 to simulate water-budget components, groundwater-flow directions, and groundwater-flow paths. The model simulated precipitation recharge, groundwater and surface-water interactions, evapotranspiration, high-volume groundwater pumping by pumping wells, and external inflows and outflows occurring along the boundaries of the alluvial aquifer. Model calibration was completed using manual and automated approaches, the latter of which assisted in quantifying model results sensitivity to input parameters. The calibrated model corresponds well with groundwater-level elevation observations, with a mean residual (observed minus simulated groundwater-level elevation) equal to −0.60 feet. Simulated groundwater base flow to streams was typically within 10 percent of base flow estimated by independent methods. Groundwater and surface-water interactions represented the largest water-budget components of the aquifer, with the second largest groundwater discharge component coming from pumping wells. Groundwater and surface-water interactions represent both the largest gain and loss terms in the water budget, because these interactions differ spatially, meaning in some areas of the model domain groundwater is being recharged by streams, whereas in other areas, groundwater is discharged to streams. Estimates of advective groundwater-flow paths indicate pumping wells may capture groundwater recharged from losing streams and groundwater that flows into the main stem of the alluvial aquifer from paleochannels.

Colorado↗

USGS2014

Hydrosalinity studies of the Virgin River, Dixie Hot Springs, and Littlefield Springs, Utah, Arizona, and Nevada

The Virgin River contributes a substantial amount of dissolved solids (salt) to the Colorado River at Lake Mead in the lower Colorado River Basin. Degradation of Colorado River water by the addition of dissolved solids from the Virgin River affects the suitability of the water for municipal, industrial, and agricultural use within the basin. Dixie Hot Springs in Utah are a major localized source of dissolved solids discharging to the Virgin River. The average measured discharge from Dixie Hot Springs during 2009–10 was 11.0 cubic feet per second (ft 3 /s), and the average dissolved-solids concentration was 9,220 milligrams per liter (mg/L). The average dissolved-solids load—a measurement that describes the mass of salt that is transported per unit of time—from Dixie Hot Springs during this period was 96,200 tons per year (ton/yr). Annual dissolved-solids loads were estimated at 13 monitoring sites in the Virgin River Basin from streamflow data and discrete measurements of dissolved-solids concentrations and (or) specific conductance. Eight of the sites had the data needed to estimate annual dissolved-solids loads for water years (WYs) 1999 through 2010. During 1999–2010, the smallest dissolved-solids loads in the Virgin River were upstream of Dixie Hot Springs (59,900 ton/yr, on average) and the largest loads were downstream of Littlefield Springs (298,200 ton/yr, on average). Annual dissolved-solids loads were smallest during 2002–03, which was a period of below normal precipitation. Annual dissolved-solids loads were largest during 2005—a year that included a winter rain storm that resulted in flooding throughout much of the Virgin River Basin. An average seepage loss of 26.7 ft 3 /s was calculated from analysis of monthly average streamflow from July 1998 to September 2010 in the Virgin River for the reach that extends from just upstream of the Utah/Arizona State line to just above the Virgin River Gorge Narrows. Seepage losses from three river reaches in the Virgin River Gorge containing known fault zones accounted for about 48 percent of this total seepage loss. An additional seepage loss of 6.7 ft 3 /s was calculated for the reach of the Virgin River between Bloomington, Utah, and the Utah/Arizona State line. This loss in flow is small compared to total flow in the river and is comparable to the rated error in streamflow measurements in this reach; consequently, it should be used with caution. Littlefield Springs were studied to determine the fraction of its discharge that originates as upstream seepage from the Virgin River and residence time of this water in the subsurface. Geochemical and environmental tracer data from groundwater and surface-water sites in the Virgin River Gorge area suggest that discharge from Littlefield Springs is a mixture of modern (post-1950s) seepage from the Virgin River upstream of the springs and older groundwater from a regional carbonate aquifer. Concentrations of the chlorofluorocarbons (CFCs) CFC-12 and CFC-113, chloride/fluoride and chloride/bromide ratios, and the stable isotope deuterium indicate that water discharging from Littlefield Springs is about 60 percent seepage from the Virgin River and about 40 percent discharge from the regional carbonate aquifer. The river seepage component was determined to have an average subsurface traveltime of about 26 ±1.6 years before discharging at Littlefield Springs. Radiocarbon data for Littlefield Springs suggest groundwater ages from 1,000 to 9,000 years. Because these are mixed waters, the component of discharge from the carbonate aquifer is likely much older than the groundwater ages suggested by the Littlefield Springs samples. If the dissolved-solids load from Dixie Hot Springs to the Virgin River were reduced, the irrigation water subsequently applied to agricultural fields in the St. George and Washington areas, which originates as water from the Virgin River downstream of Dixie Hot Springs, would have a lower dissolved-solids concentration. Dissolved-solids concentrations in excess irrigation water draining from the agricultural fields are about 1,700 mg/L higher than the concentrations in the Virgin River water that is currently (2014) used for irrigation that contains inflow from Dixie Hot Springs; this increase results from evaporative concentration and dissolution of mineral salts in the irrigated agricultural fields. The water samples collected from drains downgradient from the irrigated areas are assumed to include the dissolution of all available minerals precipitated in the soil during the previous irrigation season. Based on this assumption, a change to more dilute irrigation water will not dissolve additional minerals and increase the dissolved-solids load in the drain discharge. Following the hypothetical reduction of salts from Dixie Hot Springs, which would result in more dilute Virgin River irrigation water than is currently used, the dissolution of minerals left in the soil from the previous irrigation season would result in a net increase in dissolved-solids concentrations in the drain discharge, but this increase should only last one irrigation season. After one (or several) seasons of irrigating with more dilute irrigation water, mineral precipitation and subsequent re-dissolution beneath the agricultural fields should be greatly reduced, leading to a reduction in dissolved-solids load to the Virgin River below the agricultural drains. A mass-balance model was used to predict changes in the dissolved-solids load in the Virgin River if the salt discharging from Dixie Hot Springs were reduced or removed. Assuming that 33.4 or 26.7 ft 3 /s of water seeps from the Virgin River to the groundwater system upstream of the Virgin River Gorge Narrows, the immediate hypothetical reduction in dissolved-solids load in the Virgin River at Littlefield, Arizona is estimated to be 67,700 or 71,500 ton/yr, respectively. The decrease in dissolved-solids load in seepage from the Virgin River to the groundwater system is expected to reduce the load discharging from Littlefield Springs in approximately 26 years, the estimated time lag between seepage from the river and discharge of the seepage water, after subsurface transport, from Littlefield Springs. At that time, the entire reduction in dissolved solids seeping from the Virgin River is expected to be realized as a reduction in dissolved solids discharging from Littlefield Springs, resulting in an additional reduction of 24,700 ton/yr (based on 33.4 ft 3 /s of seepage loss) or 21,000 ton/yr (based on 26.7 ft 3 /s of seepage loss) in the river’s dissolved-solids load at Littlefield.

Arizona, Nevada, Utah↗

USGS2014

Concentrations and transport of suspended sediment, nutrients, and pesticides in the lower Mississippi-Atchafalaya River subbasin during the 2011 Mississippi River flood, April through July

High streamflow associated with the April–July 2011 Mississippi River flood forced the simultaneous opening of the three major flood-control structures in the lower Mississippi-Atchafalaya River subbasin for the first time in history in order to manage the amount of water moving through the system. The U.S. Geological Survey (USGS) collected samples for analysis of field properties, suspended-sediment concentration, particle-size, total nitrogen, nitrate plus nitrite, total phosphorus, orthophosphate, and up to 136 pesticides at 11 water-quality stations and 2 flood-control structures in the lower Mississippi-Atchafalaya River subbasin from just above the confluence of the upper Mississippi and Ohio Rivers downstream from April through July 2011. Monthly fluxes of suspended sediment, suspended sand, total nitrogen, nitrate plus nitrite, total phosphorus, orthophosphate, atrazine, simazine, metolachlor, and acetochlor were estimated at 9 stations and 2 flood-control structures during the flood period. Although concentrations during the 2011 flood were within the range of what has been observed historically, concentrations decreased during peak streamflow on the lower Mississippi River. Prior to the 2011 flood, high concentrations of suspended sediment and nitrate were observed in March 2011 at stations downstream of the confluence of the upper Mississippi and Ohio Rivers, which probably resulted in a loss of available material for movement during the flood. In addition, the major contributor of streamflow to the lower Mississippi-Atchafalaya River subbasin during April and May was the Ohio River, whose water contained lower concentrations of suspended sediment, pesticides, and nutrients than water from the upper Mississippi River. Estimated fluxes for the 4-month flood period were still quite high and contributed approximately 50 percent of the estimated annual suspended sediment, nitrate, and total phosphorus fluxes in 2011; the largest fluxes were estimated at the water-quality station located at Vicksburg, Mississippi. The majority of the suspended-sediment flux introduce into the lower Mississippi-Atchafalaya River subbasin during the 2011 flood was in the form of fine-grained particles from the upper Mississippi River—77 percent of the suspended-sediment flux compared to 23 percent from the Ohio River. As water moved downstream along the lower Mississippi River, there were losses in suspended-sediment flux because of deposition and backwater areas. Fluxes showed a greater response to increased streamflow in the Atchafalaya River than in the lower Mississippi River. The result was a gain in suspended-sediment flux with distance downstream in the Atchafalaya River because of resuspension of previously deposited materials—particularly sand particles. Overall, 13 percent less suspended sediment left the lower Mississippi-Atchafalaya River subbasin than entered it from the confluence of the upper Mississippi and Ohio Rivers during the flood. The loss in suspended-sediment flux during the flood accounted for 14 percent of the 2011 annual suspended-sediment flux loss within the lower Mississippi-Atchafalaya River subbasin. Nitrate composed approximately 70 percent of the total nitrogen flux at all of the sampled water-quality stations, excluding the Arkansas River. Almost 2.4 times more nitrate flux entered the lower Mississippi-Atchafalaya River subbasin from the upper Mississippi River than from the Ohio River. As nitrate moved down the lower Mississippi River and the Atchafalaya River, there were no substantial losses or gains in flux, indicating that nitrate moved conservatively within the subbasin during the 2011 flood. Although streamflow was the largest on record, nitrate flux during the flood period resulted in a zone of hypoxia in the Gulf of Mexico that was only the tenth largest on record. The flux of total phosphorus in the lower Mississippi-Atchafalaya River subbasin during the 2011 flood was strongly related to suspended-sediment flux at most of the stations. There were significant gains in total phosphorus flux in the Atchafalaya River during the flood period and losses between the stations along the lower Mississippi River. Overall, however, the amount of total phosphorus flux that left the lower Mississippi-Atchafalaya River subbasin was only 1.7 percent less than the flux that entered it from the upper Mississippi River and the Ohio River, indicating that total phosphorus flux within the subbasin during the flood was conservative. As streamflow was decreasing within the lower Mississippi-Atchafalaya River subbasin, orthophosphate composed an increasing percentage of the total phosphorus concentration, probably because of the return of waters low in oxygen concentration from areas such as inundated lands, backwater streams, and floodways. Poorly oxygenated waters promote the release of sediment-bound phosphorus into the more-readily available dissolved form (measured as orthophosphate in this study). Because of processing within the subbasin during the flood period, there was a 25-percent gain in orthophosphate flux between the confluence of the upper Mississippi and Ohio Rivers and the outlet of the subbasin. Of the 136 pesticide compounds and degradates that were analyzed, only 18 were detected above the method reporting level. The 18 compounds that were detected fell into three categories: (1) compounds that were frequently detected and showed a response in concentration to the flood; (2) compounds that were detected in almost every sample at every station but at low concentrations; and (3) compounds that were infrequently detected. Fluxes for the most frequently detected pesticides having the highest concentrations (atrazine, metolachlor, acetochlor, and simazine) were within the low-to-middle range of historic fluxes. An average of 66,450 cubic feet per second of streamflow was diverted from the lower Mississippi River through the Morganza Floodway into the Atchafalaya River from May 14 through July 7, 2011. Dissolved oxygen concentrations in the floodway decreased with the amount of time that the flood control structure was open, which affected nitrate and orthophosphate concentrations. As dissolved oxygen concentrations decreased in the floodway, nitrate concentrations decreased and orthophosphate concentrations increased. Oil and gas samples were also collected at 1 station upstream and 1 station downstream from the outlet of the Morganza Floodway into the Atchafalaya River. There were no detections of petroleum hydrocarbons in the upstream or downstream samples. All concentrations of oil and grease were relatively low, and the effect of water from the floodway on water quality in the Atchafalaya River could not be determined because oil and grease samples were not collected from the floodway.

Lower Mississippi-atchafalaya River Subbasin;Missi↗

USGS2012

Simulation of groundwater conditions and streamflow depletion to evaluate water availability in a Freeport, Maine, watershed

In order to evaluate water availability in the State of Maine, the U.S. Geological Survey (USGS) and the Maine Geological Survey began a cooperative investigation to provide the first rigorous evaluation of watersheds deemed "at risk" because of the combination of instream flow requirements and proportionally large water withdrawals. The study area for this investigation includes the Harvey and Merrill Brook watersheds and the Freeport aquifer in the towns of Freeport, Pownal, and Yarmouth, Maine. A numerical groundwater- flow model was used to evaluate groundwater withdrawals, groundwater-surface-water interactions, and the effect of water-management practices on streamflow. The water budget illustrates the effect that groundwater withdrawals have on streamflow and the movement of water within the system. Streamflow measurements were made following standard USGS techniques, from May through September 2009 at one site in the Merrill Brook watershed and four sites in the Harvey Brook watershed. A record-extension technique was applied to estimate long-term monthly streamflows at each of the five sites. The conceptual model of the groundwater system consists of a deep, confined aquifer (the Freeport aquifer) in a buried valley that trends through the middle of the study area, covered by a discontinuous confining unit, and topped by a thin upper saturated zone that is a mixture of sandy units, till, and weathered clay. Harvey and Merrill Brooks flow southward through the study area, and receive groundwater discharge from the upper saturated zone and from the deep aquifer through previously unknown discontinuities in the confining unit. The Freeport aquifer gets most of its recharge from local seepage around the edges of the confining unit, the remainder is received as inflow from the north within the buried valley. Groundwater withdrawals from the Freeport aquifer in the study area were obtained from the local water utility and estimated for other categories. Overall, the public-supply withdrawals (105.5 million gallons per year (Mgal/yr)) were much greater than those for any other category, being almost 7 times greater than all domestic well withdrawals (15.3 Mgal/yr). Industrial withdrawals in the study area (2.0 Mgal/yr) are mostly by a company that withdraws from an aquifer at the edge of the Merrill Brook watershed. Commercial withdrawals are very small (1.0 Mgal/yr), and no irrigation or other agricultural withdrawals were identified in this study area. A three-dimensional, steady-state groundwater-flow model was developed to evaluate stream-aquifer interactions and streamflow depletion from pumping, to help refine the conceptual model, and to predict changes in streamflow resulting from changes in pumping and recharge. Groundwater levels and flow in the Freeport aquifer study area were simulated with the three-dimensional, finite-difference groundwater-flow modeling code, MODFLOW-2005. Study area hydrology was simulated with a 3-layer model, under steady-state conditions. The groundwater model was used to evaluate changes that could occur in the water budgets of three parts of the local hydrologic system (the Harvey Brook watershed, the Merrill Brook watershed, and the buried aquifer from which pumping occurs) under several different climatic and pumping scenarios. The scenarios were (1) no pumping well withdrawals; (2) current (2009) pumping, but simulated drought conditions (20-percent reduction in recharge); (3) current (2009) recharge, but a 50-percent increase in pumping well withdrawals for public supply; and (4) drought conditions and increased pumping combined. In simulated drought situations, the overall recharge to the buried valley is about 15 percent less and the total amount of streamflow in the model area is reduced by about 19 percent. Without pumping, infiltration to the buried valley aquifer around the confining unit decreased by a small amount (0.05 million gallons per day (Mgal/d)), and discharge to the streams increased by about 8 percent (0.3 Mgal/d). A 50-percent increase in pumping resulted in a simulated decrease in streamflow discharge of about 4 percent (0.14 Mgal/d). Streamflow depletion in Harvey Brook was evaluated by use of the numerical groundwater-flow model and an analytical model. The analytical model estimated negligible depletion from Harvey Brook under current (2009) pumping conditions, whereas the numerical model estimated that flow to Harvey Brook decreased 0.38 cubic feet per second (ft 3 /s) because of the pumping well withdrawals. A sensitivity analysis of the analytical model method showed that conducting a cursory evaluation using an analytical model of streamflow depletion using available information may result in a very wide range in results, depending on how well the hydraulic conductivity variables and aquifer geometry of the system are known, and how well the aquifer fits the assumptions of the model. Using the analytical model to evaluate the streamflow depletion with an incomplete understanding of the hydrologic system gave results that seem unlikely to reflect actual streamflow depletion in the Freeport aquifer study area. In contrast, the groundwater-flow model was a more robust method of evaluating the amount of streamflow depletion that results from withdrawals in the Freeport aquifer, and could be used to evaluate streamflow depletion in both streams. Simulations of streamflow without pumping for each measurement site were compared to the calibratedmodel streamflow (with pumping), the difference in the total being streamflow depletion. Simulations without pumping resulted in a simulated increase in the steady-state flow rate of 0.38 ft 3 /s in Harvey Brook and 0.01 ft 3 /s in Merrill Brook. This translates into a streamflow-depletion amount equal to about 8.5 percent of the steady-state base flow in Harvey Brook, and an unmeasurable amount of depletion in Merrill Brook. If pumping was increased by 50 percent and recharge reduced by 20 percent, the amount of streamflow depletion in Harvey Brook could reach 1.41 ft 3 /s.

Maine↗

USGS2016

Temporal trends and stationarity in annual peak flow and peak-flow timing for selected long-term streamflow-gaging stations in or near Montana through water year 2011: Chapter B in Montana StreamStats

A large-scale study by the U.S. Geological Survey, in cooperation with the Montana Department of Transportation and the Montana Department of Natural Resources and Conservation, was done to investigate general patterns in peak-flow temporal trends and stationarity through water year 2011 for 24 long-term streamflow-gaging stations (hereinafter referred to as gaging stations) in Montana. Hereinafter, all years refer to water years; a water year is the 12-month period from October 1 through September 30 and is designated by the year in which it ends. The primary focus of the study was to identify general patterns in peak-flow temporal trends and stationarity that are relevant to application of peak-flow frequency analyses within a statewide gaging-station network. Temporal trends were analyzed for two hydrologic variables: annual peak flow and peak-flow timing. Annual peak flow is the maximum instantaneous discharge, in cubic feet per second, recorded each year a gaging station was operated. Peak-flow timing is the day of the annual peak flow (hereinafter referred to as day of peak), recorded each year a gaging station was operated. Study results provide evidence that annual peak flow for most of the long-term gaging stations can be reasonably considered as stationary for application of peak-flow frequency analyses within a statewide gaging station network. Upward trends in annual peak flow during 1930–76 generally were stronger than downward trends during 1967–2011 for most long-term gaging stations. Statistical distributions of annual peak flow generally were similar among three summary time periods (1930–78, 1979–2011, and the entire period of record). However, for two low-elevation gaging stations in eastern Montana (Poplar River at international boundary [gaging station 06178000] and Powder River at Moorhead, Montana [gaging station 06324500]), substantial downward trends in annual peak flow during 1967–2011 were of similar or stronger magnitude than the upward trends during 1930–76, and the annual-peak-flow medians for 1979–2011 were substantially lower than the medians for the entire period of record. For peak-flow timing for most long-term gaging stations, differences in trends between 1930–76 and 1967–2011 are variable and not particularly strong. Statistical distributions generally are similar among the summary time periods. However, for two high-elevation gaging stations on a Missouri River headwater tributary (Gallatin River near Gallatin Gateway, Montana [gaging station 06043500] and Gallatin River at Logan, Montana [gaging station 06052500]) and for five high-elevation gaging stations in the Yellowstone River Basin (Yellowstone River at Corwin Springs, Montana [gaging station 06191500], Yellowstone River near Livingston, Montana [gaging station 06192500], Clarks Fork Yellowstone River near Belfry, Montana [gaging station 06207500], Clarks Fork Yellowstone River at Edgar, Montana [gaging station 06208500], and Yellowstone River at Billings, Montana [gaging station 06214500]) downward trends in peak-flow timing during 1967–2011 generally were stronger than upward trends during 1930–1976, and day-of-peak medians for 1979–2011 were considerably less than medians for 1930–78. The downward trends in peak-flow timing for 1967–2011 indicate that the timing of annual peak flows changed from later in the year to earlier in the year. For the seven high-elevation gaging stations, the mean change during 1967–2011 was about 13 days (range of 8 to 22 days). For most of the high-elevation gaging stations in the Missouri River headwaters, Yellowstone River Basin, and Columbia River Basin, there was general correspondence between trend patterns for annual peak flow and trend patterns for peak-flow timing; that is, during periods when there were upward trends in annual peak flow, there generally also were upward trends in peak-flow timing. Conversely, during periods when there were downward trends in annual peak flow, there generally also were downward trends in peak-flow timing. The two low-elevation gaging stations in eastern Montana (Poplar River at international boundary [gaging station 06178000] and Powder River at Moorhead, Montana [gaging station 06324500]) had considerable changes in annual-peakflow characteristics after the mid-1970s, which might provide evidence of potential nonstationarity in the peak-flow records. The two low-elevation gaging stations that have potential nonstationarity are located in drainage basins that are strongly affected by agricultural activities that potentially affect the hydrologic regimes. Primary agricultural activities that might alter natural hydrologic conditions include construction of small impoundments (primarily for stock-watering purposes) and irrigation diversions. Temporal variability in these activities might contribute to the potential nonstationarity issues. Changes in climatic characteristics after the mid-1970s also possibly contribute to the potential nonstationarity issues. Lack of considerable indication of potential nonstationarity in annual peak flow for the other long-term gaging stations in this study might indicate that climatic changes have been more pronounced with respect to effects on peak flows in low elevation areas in eastern Montana than in areas represented by the other long-term gaging stations. Another possibility is that climatic changes after the mid-1970s are exacerbated in low-elevation areas where small-impoundment development and potential effects of irrigation diversions might be more extensive.

Montana↗

USGS2015

Climate and streamflow characteristics for selected streamgages in eastern South Dakota, water years 1945–2013

Upward trends in precipitation and streamflow have been observed in the northeastern Missouri River Basin during the past century, including the area of eastern South Dakota. Some of the identified upward trends were anomalously large relative to surrounding parts of the northern Great Plains. Forcing factors for streamflow trends in eastern South Dakota are not well understood, and it is not known whether streamflow trends are driven primarily by climatic changes or various land-use changes. Understanding the effects that climate (specifically precipitation and temperature) has on streamflow characteristics within a region will help to better understand additional factors such as land-use alterations that may affect the hydrology of the region. To aid in this understanding, a study was completed by the U.S. Geological Survey, in cooperation with the East Dakota Water Development District and James River Water Development District, to assess trends in climate and streamflow characteristics at 10 selected streamgages in eastern South Dakota for water years (WYs) 1945–2013 (69 years) and WYs 1980–2013 (34 years). A WY is the 12-month period, October 1 through September 30, and is designated by the calendar year in which it ends. One streamgage is on the Whetstone River, a tributary to the Minnesota River, and the other streamgages are in the James, Big Sioux, and Vermillion River Basins. The watersheds for two of the James River streamgages extend into North Dakota, and parts of the watersheds for two of the Big Sioux River streamgages extend into Minnesota and Iowa. The objectives of this study were to document trends in streamflow and precipitation in these watersheds, and characterize the residual streamflow variability that might be attributed to factors other than precipitation. Residuals were computed as the departure from a locally-weighted scatterplot smoothing (LOWESS) model. Significance of trends was based on the Mann-Kendall nonparametric test at a 0.10 significance level. Of the 10 streamgages selected, only the Elm River at Westport (in the upper part of James River Basin) did not have a significant upward trend in annual mean streamflow for WYs 1945–2013, whereas only one-half of the streamgages had significant upward trends in annual mean streamflow for WYs 1980–2013. Mean and 7-day minimum streamflows also had upward trends for the spring runoff period (March–May) for most of the streamgages during WYs 1945–2013 and for one streamgage during WYs 1980–2013. Magnitudes of increases in streamflow were as great as 30 cubic feet per second per year for the streamgage on the James River near Scotland during WYs 1980–2013. Precipitation trends for WYs 1945–2013 were not necessarily significant for the watersheds of streamgages with a significant streamflow trend. Annual total precipitation had a significant upward trend for the watersheds of 4 of the 10 streamgages during WYs 1945–2013 and no significant trends for WYs 1980–2013. The most widespread precipitation increase was for September–November, with significant upward trends for the watersheds of 8 of the 10 streamgages during WYs 1945–2013; however, no trends in September– November precipitation were significant for WYs 1980–2013. The greatest magnitude of increase in precipitation was for the December–May season during WYs 1980–2013, which had a mean increase of 0.106 inch per year in the watersheds of streamgages with significant trends. The correlation between streamflow and precipitation metrics was low as indicated by the mean coefficient of determination (R 2 ) of 0.18 for all pairs considered. The highest locally-weighed scatterplot smoothing (LOWESS) correlation was between annual precipitation (by water year) and annual mean streamflow (by water year), which had a mean R 2 of 0.47 for all streamgages and was as high as 0.72 for one streamgage. The correlation between annual precipitation and March–May mean streamflow had a mean R 2 of 0.33 for all streamgages and was as high as 0.52 for one streamgage. Other metrics had R 2 values for LOWESS correlations that were less than 0.3 and were not further considered for analyses of residuals. For annual precipitation as a predictor of annual mean flow, precipitation-removed streamflow had significant upward trends during WYs 1945–2013 for one-half of the streamgages. Upward trends in residual annual mean streamflow were indicated for the Whetstone River and lower part of the Big Sioux River Basin, indicating that other factors are contributors to streamflow variability during WYs 1945–2013. In contrast, most of the streamgages in the James and Vermillion River Basins had no trends in residual annual mean streamflow, indicating that streamflow trends can be explained primarily by precipitation. Precipitation-removed streamflow had fewer trends during the more recent analysis period of WYs 1980–2013 than WYs 1945–2013 for all streamgages in eastern South Dakota. Upward trends in residuals for March– May mean streamflow were indicated for Skunk Creek at Sioux Falls and the Big Sioux River at Akron, but trends in residuals were not significant at the remaining streamgages. For the streamgages with significant trends in residual streamflow (such as the streamgage on the Whetstone River and streamgages in the Big Sioux River Basin), land-use changes likely are minor factors, with the main factors probably being changes in the timing and frequency of large precipitation events and persistently wetter antecedent conditions. Changes in the relation between precipitation and streamflow since 1945 were evident when considering the runoff efficiency of the watershed. For example, the streamflow response to annual precipitation of 25 inches for the James River near Scotland increased from approximately 1,000 cubic feet per second for WYs 1945–1990 to about 2,500 cubic feet per second for WYs 1991–2013. The importance of antecedent conditions on annual mean streamflow also was indicated by the significance of the multiple linear regression coefficients of annual mean streamflow and precipitation from preceding water years for all but one streamgage. In addition, rising groundwater levels are present in wells in eastern South Dakota, particularly since the 1980s.

South Dakota↗

USGS2012

Simulation of groundwater flow and hydrologic effects of groundwater withdrawals from the Kirkwood-Cohansey aquifer system in the Pinelands of southern New Jersey

The Kirkwood-Cohansey aquifer system is an important source of present and future water supply in southern New Jersey. Because this unconfined aquifer system also supports sensitive wetland and aquatic habitats within the New Jersey Pinelands (Pinelands), water managers and policy makers need up-to-date information, data, and projections that show the effects of potential increases in groundwater withdrawals on these habitats. Finite-difference groundwater flow models (MODFLOW) were constructed for three drainage basins (McDonalds Branch Basin, 14.3 square kilometers (km 2 ); Morses Mill Stream Basin, 21.63 km 2 ; and Albertson Brook Basin, 52.27 km 2 ) to estimate the effects of potential increases in groundwater withdrawals on water levels and the base-flow portion of streamflow, in wetland and aquatic habitats. Three models were constructed for each drainage basin: a transient model consisting of twenty-four 1-month stress periods (October 2004 through September 2006); a transient model to simulate the 5- to 10-day aquifer tests that were performed as part of the study; and a high-resolution, steady-state model used to assess long-term effects of increased groundwater withdrawals on water levels in wetlands and on base flow. All models were constructed with the same eight-layer structure. The smallest horizontal cell dimensions among the three model areas were 150 meters (m) for the 24-month transient models, 10 m for the steady-state models, and 3 m for the transient aquifer-test models. Boundary flows of particular interest to this study and represented separately are those for wetlands, streams, and evapotranspiration. The final variables calibrated from both transient models were then used in steady-state models to assess the long-term effects of increased groundwater withdrawals on water levels in wetlands and on base flow. Results of aquifer tests conducted in the three study areas illustrate the effects of withdrawals on water levels in wetlands and on base flow. Pumping stresses at aquifer-test sites resulted in measurable drawdown in each observation well installed for the tests. The magnitude of drawdown in shallow wetland observation wells at the end of pumping ranged from 5.5 to 16.7 centimeters (cm). The stresses induced by the respective tests reduced the flow of the smallest stream (McDonalds Branch) by 75 percent and slightly reduced flow in a side channel of Morses Mill Stream, but did not measurably affect the flow of Morses Mill Stream or Albertson Brook. Results of aquifer-test simulations were used to refine the estimates of hydraulic properties used in the models and to confirm the ability of the model to replicate observed hydrologic responses to pumping. Steady-state sensitivity simulation results for a variety of single well locations and depths were used to define overall “best-case” (smallest effect on wetland water levels and base flow) and “worst-case” (greatest effect on wetland water levels and base flow) groundwater withdrawal configurations. “Best-case” configurations are those for which the extent of the wetland areas within a 1-kilometer (km) radius of the withdrawal well is minimized, the well is located at least 100 m and as far from wetland boundaries as possible, and the withdrawal is from a deep well (50–90 m deep). “Worst-case” configurations are those for which the extent of wetlands within a 1-km radius of the withdrawal well is maximized, the well is located 100 m or less from a wetland boundary, and the withdrawal is from a relatively shallow well (30–67 m deep). “Best-” and “worst-case” simulations were applied by locating hypothetical wells across the study areas and assigning groundwater withdrawals so that the sum of the withdrawals for the basin is equal to 5, 10, 15, and 30 percent of overall recharge. The results were compared to the results of simulations of no groundwater withdrawals. Results for withdrawals of 5 percent of recharge show that the area of wetland water-level decline that exceeded 15 cm was as much as 1.5 percent of the total wetland area for the “best-case” simulations and as much as 9.7 percent of the total wetland area for the “worst-case” simulations. For the same withdrawals, base-flow reduction was as much as 5.1 percent for the “best-case” simulations and as much as 8.6 percent for the “worst-case” simulations. Results for withdrawals of 30 percent of recharge show that the area of wetland water-level decline that exceeded 15 cm was as much as 70 percent of the total wetland area for the “best-case” simulations and as much as 84 percent of the total wetland area for the “worst-case” simulations. For the same withdrawals, base-flow reduction was as much as 30 percent for the “best-case” simulations and as much as 51 percent for the “worst-case” simulations. Results for withdrawals of 10 and 15 percent of recharge show decreased water levels and base flow that are intermediate between those simulated for 5 and 30 percent of recharge. Several approaches for applying the results of this study to other parts of the Pinelands were explored. An analytical-modeling technique based on the Thiem equation and image-well theory was developed to estimate local drawdown distributions resulting from withdrawals in other areas within the Pinelands. Results of example applications of this technique were compared with those of the numerical simulations used in this study and were shown to be useful. Differences among the three basins in the simulated percentage of basin wetlands affected by drawdown were found to be related to the proximity of wetlands to streams, the proximity of wetlands to pumped wells, and the vertical conductance of the aquifer system. These factors formed the basis for an index of wetland vulnerability to drawdown. An empirically-derived model based on the Gompertz function and the wetland vulnerability index was developed, tested, and shown to be an effective means to evaluate potential drawdown in wetlands at a basin scale throughout the Pinelands. Base-flow reduction can be estimated from generalized results of the numerical models, estimates of evapotranspiration reduction, or available regional groundwater flow models. These approaches could be used to evaluate alternative water-supply strategies and, in conjunction with ecological-modeling results, to determine maximum basin withdrawal rates within the limits of acceptable ecological change.

New Jersey↗

USGS1976

Summary appraisals of the nation's ground-water resources – California region

Most people in the California Region live in a semiarid or arid climate, with precipitation less than the potential evapotranspiration- environments of perennial water deficiency. The deficiency becomes most onerous during the characteristically rainless summers and during recurrent droughts that may continue for 10--20 years. However, water from winter rain and snow can be stored for use during the dry summer months, and water stored during a wet climatic period can be used in a succeeding dry period; moreover, perennial deficiency can be overcome by bringing water from areas of perennial surplus. Ground-water reservoirs have especial significance in arid and semiarid regions as repositories where water is stored or can be stored with minimum loss by evaporation. Nearly all the ground-water reservoirs of the California Region are in alluvial sediments of valleys and plains that flank the mountain ranges. The largest, underlying the vast Central Valley, occupies 10 percent of the area of the region, has an estimated usable capacity exceeding 100 million acre-feet (125 cubic kilometres), and has an annual pumpage from wells of about 13 million acre-feet (16 cubic kilometres). Another 10 percent of the region is occupied by 55 developed ground-water reservoirs that are widely distributed; aggregate annual pumpage from them is about 3 1/2 million acre-feet (4 cubic kilometres). In the southeastern desert about 60 ground-water reservoirs occupy still another 10 percent of the region; these have been explored only enough to show that most have some usable water, but current use is negligible. In northeastern California and adjacent Oregon and Nevada, ground-water reservoirs are identified only in valleys and lowlands where wells are feasible, but basaltic rocks of the Cascade Range and Modoc Plateau are excellent aquifers distributed over an area constituting about 15 percent of the region. In sum, slightly less than half the California Region is underlain by ground-water reservoirs, either in valley fill or in volcanic rocks, which can yield significant quantities of water to wells. The rest of the California Region includes the mountains, canyons, slopes, and foothills of the Sierra Nevada, Coast Ranges, and Basin Ranges, whose consolidated rocks and products of their weathering may be permeable locally but are not generally so. Here, the prevailing method of ground-water development is still mostly trial and error, and while in many places a well can yield enough water for a family, some families might have to do without amenities such as flush toilets and automatic washers. For more than half a century the California Region has led all others in North America in pumping of ground water as well as in the area, variety, yield, and export of crops irrigated by water from wells. It has led in the development and use of deep-well turbine pumps for large yield and in the drilling of water wells to great depths. At the same time, such developments have resulted in the elimination of artesian pressures that produced thousands of flowing wells in the 19th century and led to the wide distribution of "falling water tables." Also, California was first to induce encroachment of seawater into wells (in 1906); first to recognize subsidence of land caused by pumping from wells (in 1933), generating news about land sinking in San Jose, Long Beach, and along the Delta-Mendota and Friant-Kern Canals; and first to experience pollution of ground-water reservoirs by brines, chemicals, industrial wastes, and petroleum byproducts including gasoline. The region has led in research in several fields leading to solution of many of these problems. Ground-water problems developed rapidly after World War II with booming population, agriculture, industry, and water demand during several years of regionwide drought. Water levels in wells trended downward almost everywhere as a natural effect of the drought and at accelerated rates in areas of pumping for new enterprises or to supplement subnormal surface-water supplies. The declines in many pumping areas exceeded 100 feet (30 metres), and in some confined aquifers the potentiometric surface was drawn down more than 330 feet (100 metres). The depletion of ground-water storage has had "permanent" side effects, including subsidence of the land exceeding 10 feet (3 metres) in extensive areas, and seawater intrusion that ended the useful lives of many wells along the coast and as much as 6 miles (10 kilometres) inland. Some problems have been solved, but these solutions have at times created other problems. Many ground-water reservoirs have gone through one or more stages - exploration for productive aquifers, exploitation and development for use of the water, restriction to the perennial supply or "safe" yield, importation of surface water, artificial recharge of ground water, conjunctive use of surface and ground water, protection of water quality, and integrated management of use and disposal of water. This evolutionary sequence is unique for each reservoir, and so generalizations become difficult in a regional appraisal; also, the changes with time are significant and varied, and knowledge of prior events is a prerequisite in an appraisal of the resource in a specific year. As of 1970, water levels in many wells had risen significantly from the minimum levels of record reached during the 1960's or earlier; only in areas of new development and in desert areas of long-continued "mining" of nonreplenished water was ground-water storage still being depleted. Land subsidence has continued at diminishing rates and practically has come to a halt in some areas; invading seawater has been stopped or nudged back in most places where problems were significant. The current, favorable situation has been helped by climatic variations, from drought in 1945-52 and exceedingly dry years in 1959 and 1961 to above-normal precipitation in 1969 and 1970; but most of the serious problems have been solved by human efforts, including especially the implementation of the California Water Plan, transporting water from areas of perennial surplus to areas where it is used in lieu of ground water or for ground-water replenishment. All major urban areas now import water to supplement or replace the water pumped from wells. Extensive agricultural areas that formerly were irrigated solely by ground water now obtain some of their water from surface reservoirs and canals, especially in the Central Valley. With surface water available as an alternative supply, well owners can view their ground water with complacency. But complacency can lead to neglect and carelessness and consequent deterioration of the ground-water resource by pollution. Claiming heritage from the English Common Law, the existing California law grants to the landowner (riparian) and private enterprise (appropriator) rights to the water stored in ground-water reservoirs or discharged from them, including the base flow of streams. Ground-water development has been by private and local enterprise, and the California legislature has protected and encouraged local responsibility, control, and management of ground water. As to surface water, a constitutional amendment in 1928 limited riparian rights to the quantities of water that were "reasonably required for the beneficial use to be served." The surpluses have become public waters which are collected, stored, transported, and delivered under various contracts by Federal, State, and other agencies. The agencies have not stored water underground because of uncertainty as to their rights, but some local agencies have been encouraged with favorable pricing schedules to undertake the artificial recharge and management of ground-water reservoirs. Thus, conjunctive use of surface and ground water has become a matter of interagency negotiation. Of all the constraints on effective use of ground-water reservoirs, the most formidable may be the attitudes of people. Assurance of water supply is vital in areas of water deficiency, and Government has assumed increasing responsibility for the welfare of people in these areas. Unfortunately, when Government provides this assurance, most beneficiaries demand continued subsidy to the exclusion of perhaps cheaper private development. Indeed, as the water resources are presently segregated-with private rights predominant in ground water and public interest dominant in surface water-ground-water development has suffered for lack of public concern. The region has the scientific and technologic capability for effective use of groundwater reservoirs, as shown by the achievements and programs of several districts, but many districts are not organized or staffed for such comprehensive management and will need assistance and scientific expertise available from State and Federal agencies. Those agencies, in turn, may not have the scientific data that are essential to prevent haphazard activities and to enable programs to be organized for the most effective and attractive utilization of the water resources. In these days of increasing concern over pollution, existing data are generally inadequate to assess the natural deterioration of ground waters as a basis for defining pollution.

California↗

USGS2020

Hydrogeology, numerical simulation of groundwater flow, and effects of future water use and drought for reach 1 of the Washita River alluvial aquifer, Roger Mills and Custer Counties, western Oklahoma, 1980–2015

The Washita River alluvial aquifer is a valley-fill and terrace alluvial aquifer along the valley of the Washita River in western Oklahoma that provides a productive source of groundwater for agricultural irrigation and water supply. The Oklahoma Water Resources Board (OWRB) has designated the westernmost section of the aquifer in Roger Mills and Custer Counties, Okla., as reach 1 of the Washita River alluvial aquifer; reach 1 is the focus of this report. The OWRB issued an order on November 13, 1990, that established the maximum annual yield (MAY; 120,320 acre-feet per year [acre-ft/yr]) and equal-proportionate-share (EPS) pumping rate (2.0 acre-feet per acre per year [(acre-ft/acre)/yr]) for reach 1 of the Washita River alluvial aquifer. The MAY and EPS were based on hydrologic investigations that evaluated the effects of potential groundwater withdrawals on groundwater availability in the Washita River alluvial aquifer. Every 20 years, the OWRB is statutorily required to update the hydrologic investigation on which the MAY and EPS were based. Because 30 years have elapsed since the last order was issued, the U.S. Geological Survey, in cooperation with the OWRB, conducted a new hydrologic investigation and evaluated the effects of potential groundwater withdrawals on groundwater flow and availability in the Washita River alluvial aquifer. The Washita River is the primary source of inflow to Foss Reservoir, a Bureau of Reclamation reservoir constructed in 1961 for flood control, water supply, and recreation. Foss Reservoir provides water for Bessie, Clinton, New Cordell, and Hobart, Okla. Nearly 98 percent of the total groundwater use from the Washita River alluvial aquifer during 1967 to 2015 was for irrigation; other uses of groundwater in the study area include public supply, mining, and agriculture. A hydrogeologic framework was developed for the Washita River alluvial aquifer and included the physical characteristics of the aquifer, the geologic setting, the hydraulic properties of hydrogeologic units, the potentiometric surface (water table), and groundwater-flow directions at a scale that captures the regional controls on groundwater flow. The Washita River alluvial aquifer consists of alluvium and terrace deposits that were transported primarily by water and range from clay to gravel in size. The terrace includes windblown deposits of silt size and, in some cases, contains gravel laid down at several levels along former courses of present-day rivers. A conceptual flow model is a simplified description of the aquifer system that includes hydrologic boundaries, major inflow and outflow sources of the groundwater-flow system, and a conceptual water budget with the estimated mean flows between those hydrologic boundaries. During the study period 1980–2015, mean annual groundwater withdrawals, predominantly used for agricultural irrigation, totaled 5,502 acre-ft/yr, or 14 percent of aquifer outflows. When applied across the 132-square-mile aquifer area used for modeling purposes (84,366 acres), mean annual recharge of 3.15 inches per year corresponds to a mean annual recharge volume of 22,169 acre-ft/yr, or 56 percent of aquifer inflows. The annual saturated-zone evapotranspiration outflow was 11,828 acre-ft/yr for the Washita River alluvial aquifer, or about 30 percent of aquifer outflows. For the Washita River alluvial aquifer, lateral flow was 17,157 acre-ft/yr, or 44 percent of the aquifer inflows. The conceptual flow model and hydrogeologic framework were used to conceptualize, design, and build the numerical groundwater-flow model. A numerical groundwater-flow model of the Washita River alluvial aquifer was constructed by using MODFLOW-2005. The Washita River alluvial aquifer groundwater-model grid was spatially discretized into 350-foot (ft) cells and two layers. Layer 1 represented the undifferentiated alluvium and terrace deposits of Quaternary age, and layer 2 represented the bedrock of Permian age, which was given a uniform nominal thickness of 100 ft. The groundwater-simulation period was temporally discretized into 433 monthly transient stress periods, representing January 1980 to December 2015. An initial 365-day steady-state stress period was configured to represent mean annual inflows and outflows from the Washita River alluvial aquifer for the study period. The groundwater-flow model was calibrated manually and by automated adjustment of model inputs by using PEST++. Calibration targets for the Washita River alluvial aquifer model included groundwater-level observations and reservoir-stage observations, as well as base-flow and stream-seepage estimates. Three groundwater-availability scenarios were used in the calibrated groundwater model to (1) estimate the EPS pumping rate that retains the saturated thickness that meets the minimum 20-year life of the aquifer, (2) quantify the effects of projected pumping rates on groundwater storage over a 50-year period, and (3) evaluate how projected pumping rates extended 50 years into the future and sustained hypothetical drought conditions over a 10-year period affect base flow and groundwater in storage. The results of the groundwater-availability scenarios could be used by the OWRB to reevaluate the established MAY of groundwater from the Washita River alluvial aquifer. EPS scenarios for the Washita River alluvial aquifer were run for periods of 20, 40, and 50 years. The 20-, 40-, and 50-year EPS pumping rates under normal recharge conditions were 1.7, 1.6, and 1.6 (acre-ft/acre)/yr, respectively. Given the aquifer area used for modeling purposes (84,366 acres), these rates correspond to annual yields of 142,579, 134,986, and 134,986 acre-ft/yr, respectively. Groundwater storage at the end of the 20-year EPS scenario was about 281,000 acre-feet (acre-ft), or about 306,000 acre-ft (52 percent) less than the starting storage. Considering the land-surface area of the Washita River alluvial aquifer and using a specific yield of 0.12, this decrease in storage was equivalent to a mean groundwater-level decline of about 30 ft. The Washita River downstream from Foss Reservoir and most of the streams in the study area were dry at the end of the 20-year EPS scenario. Foss Reservoir stage was below the dead-pool stage of 1,597 ft after about 7 years of pumping in the 20-year EPS scenario. Four projected 50-year groundwater-use scenarios were used to simulate the effects of selected well withdrawal rates on groundwater storage in the Washita River alluvial aquifer. These four scenarios used (1) no groundwater use, (2) groundwater use at the 2015 pumping rate, (3) mean groundwater use for the simulation period, and (4) increasing groundwater use. Groundwater storage after 50 years with no groundwater use was 545,249 acre-ft, or 693 acre-ft (0.1 percent) greater than the initial groundwater storage; this groundwater storage increase is equivalent to a mean groundwater-level increase of 0.1 ft. Groundwater storage at the end of the 50-year period with 2015 pumping rates was 543,831 acre-ft, or 723 acre-ft (0.1 percent) less than the initial storage; this groundwater storage decrease is equivalent to a mean groundwater-level decrease of 0.1 ft. Groundwater storage after 50 years with the mean pumping rate for the study period was 543,202 acre-ft, or 1,349 acre-ft (0.2 percent) less than the initial groundwater storage; this groundwater storage decrease is equivalent to a mean groundwater-level decrease of 0.1 ft. Groundwater storage at the end of the 50-year period with an increasing demand groundwater-pumping rate, which was 38 percent greater than the 2015 groundwater-pumping rate, was 542,584 acre-ft, or 1,967 acre-ft (0.4 percent) less than the initial storage; this groundwater storage decrease is equivalent to a mean groundwater-level decrease of 0.2 ft. A hypothetical 10-year-drought scenario was used to simulate the effects of a prolonged period of reduced recharge on groundwater storage in the Washita River alluvial aquifer and Foss Reservoir stage and storage. To simulate the hypothetical drought, recharge in the calibrated model was reduced by 50 percent during the simulated drought period (1983–1992). Groundwater storage at the end of the drought period in December 1992 was 562,000 acre-ft, or 36,000 acre-ft (6 percent) less than the groundwater storage of the calibrated groundwater model (598,000 acre-ft). At the end of the hypothetical drought, the largest changes in saturated thickness (as great as 43.5 ft) were in the area upgradient from Foss Reservoir, particularly in the terrace at the model boundary. Substantial decreases in the Foss Reservoir stage began during the fall of 1985 in conjunction with base-flow decreases of up to 100 percent at U.S. Geological Survey streamgage 07324200 Washita River near Hammon, Okla. These lake-stage declines outpaced groundwater-level declines in the surrounding aquifer. The minimum Foss Reservoir storage simulated during the drought period was 77,954 acre-ft, which was a decrease of 46 percent from the nondrought storage.

Oklahoma↗

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↗

USGS2016

Effects of aquifer storage and recovery activities on water quality in the Little Arkansas River and Equus Beds Aquifer, south-central Kansas, 2011–14

The Equus Beds aquifer in south-central Kansas is aprimary water source for the city of Wichita. The Equus Beds aquifer storage and recovery (ASR) project was developed to help the city of Wichita meet increasing current (2016) and future water demands. The Equus Beds ASR project pumps water out of the Little Arkansas River during above-base flow conditions, treats it using drinking-water quality standards as a guideline, and recharges it into the Equus Beds aquifer for later use. Phase II of the Equus Beds ASR project currently (2016) includes a river intake facility and a surface-water treatment facility with a 30 million gallon per day capacity. Water diverted from the Little Arkansas River is delivered to an adjacent presedimentation basin for solids removal. Subsequently, waste from the surface-water treatment facility and the presedimentation basin is returned to the Little Arkansas River through a residuals return line. The U.S. Geological Survey, in cooperation with the city of Wichita, developed and implemented a hydrobiological monitoring program as part of the ASR project to characterize and quantify the effects of aquifer storage and recovery activities on the Little Arkansas River and Equus Beds aquifer water quality. Data were collected from 2 surface-water sites (one upstream and one downstream from the residuals return line), 1 residuals return line site, and 2 groundwater well sites (each having a shallow and deep part): the Little Arkansas River upstream from the ASR facility near Sedgwick, Kansas (upstream surface-water site 375350097262800), about 0.03 mile (mi) upstream from the residuals return line site; the Little Arkansas River near Sedgwick, Kans. (downstream surface-water site 07144100), about 1.68 mi downstream from the residuals return line site; discharge from the Little Arkansas River ASR facility near Sedgwick, Kansas (residuals return line site 375348097262800); 25S 01 W 07BCCC01 SMW–S11 near CW36 (MW–7 shallow groundwater well site 375327097285401); 25S01 W 07BCCC02 DMW–S10 near CW36 (MW–7 deep groundwater well site 375327097285402); 25S 01W 07BCCA01 SMW–S13 near CW36 (MW–8 shallow groundwater well site 375332097284801); and 25S 01W 07BCCA02 DMW–S14 near CW36 (MW–8 deep groundwater well site 375332097284802). The U.S. Geological Survey, in cooperation with the city of Wichita, assessed the effects of the ASR Phase II facility residuals return line discharges on stream quality of the Little Arkansas River by measuring continuous physicochemical properties and collecting discrete water-quality and sediment samples for about 2 years pre- (January 2011 through April 2013) and post-ASR (May 2013 through December 2014) Phase II facility operation upstream and downstream from the ASR Phase II facility. Additionally, habitat variables were quantified and macroinvertebrate and fish communities were sampled upstream and downstream from the ASR Phase II facility during the study period. To assess the effects of aquifer recharge on Equus Beds groundwater quality, continuous physicochemical properties were measured and discrete water-quality samples were collected before and during the onset of Phase II aquifer recharge in two (shallow and deep) groundwater wells. Little Arkansas River streamflow was about 10 times larger after the facility began operating because of greater rainfall. Residuals return line release volumes were a very minimal proportion (0.06 percent) of downstream streamflow volume during the months the ASR facility was operating. Upstream and downstream continuously measured water temperature and dissolved oxygen median differences were smaller post-ASR than pre-ASR. Turbidity generally was smaller at the downstream site throughout the study period and decreased at both sites after the ASR Phase II facility began discharging despite a median residuals return line turbidity that was about an order of magnitude larger than the median turbidity at the downstream site. Upstream and downstream continuously measured turbidity median differences were larger post-ASR than pre-ASR. Median post-ASR continuously measured nitrite plus nitrate and continuously computed total suspended solids and suspended-sediment concentrations were smaller than pre-ASR likely because of higher streamflows and dilution; whereas, median continuously computed dissolved and total organic carbon concentrations were larger likely because of higher streamflows and runoff conditions. None of the discretely measured water-quality constituents (dissolved and suspended solids, primary ions, suspended sediment, nutrients, carbon, trace elements, viral and bacterial indicators, and pesticides) in surface water were significantly different between the upstream and downstream sites after the ASR Phase II facility began discharging; however, pre-ASR calcium, sodium, hardness, manganese, and arsenate concentrations were significantly larger at the upstream site, which indicates that some water-quality conditions at the upstream and downstream sites were more similar post-ASR. Most of the primary constituents that make up dissolved solids decreased at both sites after the ASR Phase II facility began operation. Discretely collected total suspended solids concentrations were similar between the upstream and downstream sites before the facility began operating but were about 27 percent smaller at the downstream site after the facility began operating, despite the total suspended solids concentrations in the residuals return line being 15 times larger than the downstream site. Overall habitat scores were indicative of suboptimal conditions upstream and downstream from the ASR Phase II facility throughout the study period. Substrate fouling and sediment deposition mean scores indicated marginal conditions at the upstream and downstream sites during the study period, demonstrating that sediment deposition was evident pre- and post-ASR and no substantial changes in these habitat characteristics were noted after the ASR Phase II facility began discharging. Macroinvertebrate community composition (evaluated using functional feeding, behavioral, and tolerance metrics) generally was similar between sites during the study period. Fewer macroinvertebrate metrics were significant between the upstream and downstream sites post-ASR (6) than pre-ASR (14), which suggests that macroinvertebate communities were more similar after the ASR facility began discharging. Upstream-downstream comparisons in macroinvertebrate aquatic-life-support metrics had no significant differences for the post-ASR time period and neither site was fully supporting for any of the Kansas Department of Health and Environment aquatic-life-support metrics (Macroinvertebrate Biotic Index; Kansas Biotic Index with tolerances for nutrients and oxygen-demanding substances; Ephemeroptera, Plecoptera, and Trichoptera [EPT] richness; and percentage of EPT species). Overall, using macroinvertebrate aquatic life-support criteria from the Kansas Department of Health and Environment, upstream and downstream sites were classified as partially supporting before and after the onset of ASR facility operations. Fish community trophic status and tolerance groups generally were similar among sites during the study period. Fish community Little Arkansas River Basin Index of Biotic Integrity scores at the upstream and downstream sites were indicative of fair-to-good conditions before the facility began operating and decreased to fair conditions after the facility began operating. Groundwater physicochemical changes concurrent with the beginning of recharge operations at the Sedgwick basin were more pronounced in shallow groundwater. No constituent concentrations in the pre-recharge period in comparison to the post-recharge period increased to concentrations exceeding drinking water regulations; however, nitrate decreased significantly from a pre-recharge exceedance of the U.S. Environmental Protection Agency maximum contaminant level to a post recharge nonexceedance. Shallow groundwater chemical concentrations or rates of detection increased after artificial recharge began for the ions potassium, chloride, and fluoride; phosphorus and organic carbon species; trace elements barium, manganese, nickel, arsenate, arsenic, and boron; agricultural pesticides atrazine, metolachlor, metribuzin, and simazine; organic disinfection byproducts bromodichloromethane and trichloromethane; and gross beta levels. Additionally, water temperature, and pH were larger after recharge began; and total solids and slime-forming bacteria concentrations and densities were smaller. Total solids, nitrate, and selenium significantly decreased; and potassium, chloride, nickel, arsenic, fluoride, phosphorus and carbon species, and gross beta levels significantly increased in shallow groundwater after artificial recharge. Results of biological activity reaction tests indicated that water quality microbiology was different before and after artificial recharge began; at times, these differences may lead to changes in dominant bacterial populations that, in turn, may lead to formation and expansion in populations that may cause bioplugging and other unwanted effects. Calcite, iron (II) hydroxide, hydroxyapatite, and similar minerals, had shifts in saturation indices that generally were from undersaturation toward equilibrium and, in some cases, toward oversaturation. These shifts toward neutral saturation indices might suggest reduced weathering of the minerals present in the Equus Beds aquifer. Chemical weathering in the shallow parts of the aquifer may be accelerated because of the increased water temperatures and the system is more vulnerable to clogged pores and mineral dissolution as the equilibrium state is affected by recharge and withdrawal. When oversaturation is indicated for iron minerals, plugging of aquifer materials may happen.

Kansas↗

USGS2017

Hydrogeology and simulated groundwater flow and availability in the North Fork Red River aquifer, southwest Oklahoma, 1980–2013

On September 8, 1981, the Oklahoma Water Resources Board established regulatory limits on the maximum annual yield of groundwater (343,042 acre-feet per year) and equal-proportionate-share (EPS) pumping rate (1.0 acre-foot per acre per year) for the North Fork Red River aquifer. The maximum annual yield and EPS were based on a hydrologic investigation that used a numerical groundwater-flow model to evaluate the effects of potential groundwater withdrawals on groundwater availability in the North Fork Red River aquifer. The Oklahoma Water Resources Board is statutorily required (every 20 years) to update the hydrologic investigation on which the maximum annual yield and EPS were based. Because 20 years have elapsed since the final order was issued, the U.S. Geological Survey, in cooperation with the Oklahoma Water Resources Board, conducted an updated hydrologic investigation and evaluated the effects of potential groundwater withdrawals on groundwater flow and availability in the North Fork Red River aquifer in Oklahoma. This report describes a hydrologic investigation of the North Fork Red River aquifer that includes an updated summary of the aquifer hydrogeology. As part of this investigation, groundwater flow and availability were simulated by using a numerical groundwater-flow model. The North Fork Red River aquifer in Beckham, Greer, Jackson, Kiowa, and Roger Mills Counties in Oklahoma is composed of about 777 square miles (497,582 acres) of alluvium and terrace deposits along the North Fork Red River and tributaries, including Sweetwater Creek, Elk Creek, Otter Creek, and Elm Fork Red River. The North Fork Red River is the primary source of surface-water inflow to Lake Altus, which overlies the North Fork Red River aquifer. Lake Altus is a U.S. Bureau of Reclamation reservoir with the primary purpose of supplying irrigation water to the Lugert-Altus Irrigation District. A hydrogeologic framework was developed for the North Fork Red River aquifer and included a definition of the aquifer extent and potentiometric surface, as well as a description of the textural and hydraulic properties of aquifer materials. The hydrogeologic framework was used in the construction of a numerical groundwater-flow model of the North Fork Red River aquifer described in this report. A conceptual model of aquifer inflows and outflows was developed for the North Fork Red River aquifer to constrain the construction and calibration of a numerical groundwater-flow model that reasonably represented the groundwater-flow system. The conceptual-model water budget estimated mean annual inflows to and outflows from the North Fork Red River aquifer for the period 1980–2013 and included a sub-accounting of mean annual inflows and outflows for the portions of the aquifer that were upgradient and downgradient from Lake Altus. The numerical groundwater-flow model simulated the period 1980–2013 and was calibrated to water-table-altitude observations at selected wells, monthly base flow at selected streamgages, net streambed seepage as estimated for the conceptual model, and Lake Altus stage. Groundwater-availability scenarios were performed by using the calibrated numerical groundwater-flow model to (1) estimate the EPS pumping rate that guarantees a minimum 20-, 40-, and 50-year life of the aquifer, (2) quantify the potential effects of projected well withdrawals on groundwater storage over a 50-year period, and (3) simulate the potential effects of a hypothetical (10-year) drought on base flow and groundwater storage. The results of the groundwater-availability scenarios could be used by the Oklahoma Water Resources Board to reevaluate the maximum annual yield of groundwater from the North Fork Red River aquifer. EPS scenarios for the North Fork Red River aquifer were run for periods of 20, 40, and 50 years. The 20-, 40-, and 50-year EPS pumping rates under normal recharge conditions were 0.59, 0.52, and 0.52 acre-foot per acre per year, respectively. Given the 497,582-acre aquifer area, these rates correspond to annual yields of about 294,000, 259,000, and 259,000 acre-feet per year, respectively. Groundwater storage at the end of the 20-year EPS scenario was about 951,000 acre-feet, or about 1,317,000 acre-feet (58 percent) less than the starting EPS scenario storage. This decrease in storage was equivalent to a mean water-level decline of about 22 feet. Most areas of the active alluvium near the North Fork Red River, Elk Creek, and Elm Fork Red River remained partially saturated through the end of the EPS scenario because of streambed seepage. Lake Altus storage was reduced to zero after 6–7 years of EPS pumping in each scenario. Projected 50-year pumping scenarios were used to simulate the effects of selected well withdrawal rates on groundwater storage of the North Fork Red River aquifer and base flows in the North Fork Red River upstream from Lake Altus. The effects of well withdrawals were evaluated by comparing changes in groundwater storage and base flow between four 50-year scenarios using (1) no groundwater pumping, (2) mean pumping rates for the study period (1980–2013), (3) 2013 pumping rates, and (4) increasing demand pumping rates. The increasing demand pumping rates assumed a 20.4-percent increase in pumping over 50 years based on 2010–60 demand projections for southwest Oklahoma. Groundwater storage after 50 years with no pumping was about 2,606,000 acre-feet, or 137,000 acre-feet (5.5 percent) greater than the initial groundwater storage; this groundwater storage increase is equivalent to a mean water-level increase of 2.3 feet. Groundwater storage after 50 years with the mean pumping rate for the study period (1980–2013) was about 2,476,000 acre-feet, or about 7,000 acre-feet (0.3 percent) greater than the initial groundwater storage; this groundwater storage increase is equivalent to a mean water-level increase of 0.1 foot. Groundwater storage at the end of the 50-year period with 2013 pumping rates was about 2,398,000 acre-feet, or about 70,000 acre-feet (2.8 percent) less than the initial storage; this groundwater storage decrease is equivalent to a mean water-level decline of 1.2 feet. Groundwater storage at the end of the 50-year period with increasing demand pumping rates was about 2,361,000 acre-feet, or about 107,000 acre-feet (4.3 percent) less than the initial storage; this groundwater storage decrease is equivalent to a mean water-level decline of 1.8 feet. Mean annual base flow simulated at the Carter streamgage (07301500) on North Fork Red River increased by about 4,000 acre-feet (10 percent) after 50 years with no pumping and decreased by about 5,400 acre-feet (13 percent) after 50 years with increasing demand pumping rates. Mean annual base flow simulated at the North Fork Red River inflow to Lake Altus increased by about 7,400 acre-feet (15 percent) after 50 years with no pumping and decreased by about 5,800 acre-feet (12 percent) after 50 years with increasing demand pumping rates. A hypothetical 10-year drought scenario was used to simulate the effects of a prolonged period of reduced recharge on groundwater storage and Lake Altus stage and storage. Drought effects were quantified by comparing the results of the drought scenario to those of the calibrated numerical model (no drought). To simulate the hypothetical drought, recharge in the calibrated numerical model was reduced by 50 percent during the simulated drought period (1984–1993). Groundwater storage at the end of the drought period was about 2,271,000 acre-feet, or about 426,000 acre-feet (15.8 percent) less than the groundwater storage of the calibrated numerical model. This decrease in groundwater storage is equivalent to a mean water-table-altitude decline of 7.1 feet. At the end of the 10-year hypothetical drought period, base flows at the Sweetwater (07301420), Carter (07301500), Headrick (07305000), and Snyder (07307010) streamgages had decreased by about 37, 61, 44, and 45 percent, respectively. The minimum Lake Altus storage simulated during the drought period was 403 acre-feet, which was a decline of 92 percent from the nondrought storage. Reduced base flows in the North Fork Red River were the primary cause of Lake Altus storage declines.

Oklahoma↗

USGS2020

Precipitation, temperature, groundwater-level elevation, streamflow, and potential flood storage trends within the Brazos, Colorado, Big Cypress, Guadalupe, Neches, Sulphur, and Trinity River basins in Texas through 2017

The U.S. Geological Survey (USGS), in cooperation with the U.S. Army Corps of Engineers (USACE), analyzed streamflow trends and streamflow-related variables through 2017 in seven important water-supply basins to provide information that can help water managers with the USACE and river authorities make future water management decisions. The primary purpose of this report is to document trends in long-term streamflow data at 114 selected USGS streamflow-gaging stations and 36 simulated reservoir-inflow stations in 7 river basins primarily in Texas: Brazos, Colorado, Big Cypress, Guadalupe, Neches, Sulphur, and Trinity. In this report, trends were considered statistically significant if their p -values were less than or equal to 0.05 ( p -value ≤0.05). Streamflow data selected for temporal trend analyses included annual minimum streamflow, annual peak streamflow, and streamflow volume. Precipitation, air temperature, and groundwater-level-elevation data were analyzed for trends that may help to explain changes observed in the streamflow statistics. Basins were divided into sections along county lines for precipitation analyses. Streamflow volumes were analyzed for associations with potential flood storage. The potential flood storage, defined as the difference between maximum storage and normal storage, was computed for each dam from the National Inventory of Dams database and accumulated over time based on the completion date of the dam. Precipitation and air temperature trends were analyzed for each of the eight climate divisions (High Plains, Trans-Pecos, Low Rolling Hills, Edwards Plateau, North Central Texas, South Central Texas, East Texas, and Upper Coast). Results of precipitation trend analyses indicated moderate upward trends in the Upper Coast and East Texas Climate Divisions analyzed on an annual time step from 1900 through 2017. These two climate divisions are in the eastern and southeastern parts of the State, and they receive more mean annual precipitation (45.88 and 46.09 inches, respectively) than the other climate divisions. The results of air temperature analyses indicated upward trends in annual mean air temperature within all climate divisions, with a mean slope of 0.02 degree Fahrenheit per year, or 1 degree every 50 years. Within the Brazos River Basin, results of precipitation trend analyses on an annual time step indicated that precipitation amounts are most likely increasing in the lower and middle sections of the basin. Downward trends in annual streamflow and in the ratio of streamflow volume to precipitation volume were indicated at 7 of the 15 stations in the upper sections of the basin. The lower sections of the basin had mostly downward trends in annual minimum streamflow, whereas upward trends in annual minimum streamflow were indicated in the upper sections of the basin. Downward trends in annual peak streamflow were indicated at many of the stations in the upper sections of the basin. At the same seven stations in the upper sections of the basin where there were downward trends in annual streamflow, there were also downward trends in the ratio of streamflow volume to precipitation volume. The data from the same seven stations indicated negative associations between potential flood storage volume and annual streamflow volume and downward trends in the ratio of annual streamflow volume to potential flood storage volume. With the known addition of 13,006,394 acre-feet of potential flood storage between 1900 and 2010 in the subbasins analyzed, streamflow volumes have decreased in the upper sections of the Brazos River Basin. Within the Colorado River Basin, results of precipitation trend analyses on an annual time step indicated no trends in the basin. Downward trends in annual streamflow were indicated at 16 stations in the upper sections of the basin, whereas no trends in annual streamflow were indicated in the lower section of the basin. In the lower section of the basin, one station that was operated as a continuous streamflow-gaging station through 2017 had a downward trend in annual minimum streamflow, and another station (operated through 2007) had an upward trend in annual minimum streamflow. In the upper sections of the basin, data from seven stations indicated upward trends in annual minimum streamflow, and data from six stations indicated downward trends. Data from 18 stations in the upper sections of the basin indicated downward trends in annual peak streamflow. Thirteen of the 16 stations in the upper sections of the basin with data that indicated downward trends in annual streamflow also have data that indicated downward trends in the ratio of streamflow volume to precipitation volume. Data from the same 13 stations indicated negative associations between potential flood storage volume and annual streamflow volume and downward trends in the ratio of annual streamflow volume to potential flood storage volume. With the known addition of 7,193,147 acre-feet of potential flood storage between 1891 and 2014 in the subbasins analyzed, streamflow volumes have decreased in the upper sections of the Colorado River Basin. Within the Big Cypress Basin, results of precipitation trend analyses on annual, seasonal, and monthly time steps indicated almost no trends in the basin as defined in this report. However, the annual precipitation p -value only slightly exceeded the p -value threshold for a statistically significant trend. Given the upward trend in precipitation in the East Texas Climate Division, which includes the Big Cypress Basin, and the low p -value for annual precipitation within the basin, precipitation in the basin may be increasing over time. Two annual streamflow trends, one upward and one downward, were in the upper parts of the basin. Data from USGS streamflow-gaging station 07346000 Big Cypress Bayou near Jefferson, Texas, indicated an upward trend in annual minimum streamflow and a downward trend in annual peak streamflow. The station is immediately downstream from Lake O’ the Pines; presumably, minimums have increased because of regulated releases, and annual peaks have decreased because of storage from the lake for flood control. Despite the known addition of 2,737,154 acre-feet of potential flood storage between 1898 and 2011 in the subbasins analyzed, there have not been widespread reductions in streamflow volumes in the Big Cypress Basin, except for within the drainage area for the farthest upstream station on the main stem downstream from Mount Pleasant, Texas. Within the Guadalupe River Basin, results of precipitation trend analyses on an annual time step indicated an upward trend in the lower section of the basin, but no trends in annual streamflow were indicated in the lower section of the basin. In the upper section of the basin, data from 1 of the 13 stations indicated an upward trend in annual streamflow. Data from 6 of the 13 stations in the upper section of the basin indicated a trend in annual minimum streamflow with 4 upward and 2 downward trends. Data from 2 of the 13 stations in the upper section of the basin indicated downward trends in annual peak streamflow. Despite the known addition of 2,016,534 acre-feet of potential flood storage between 1849 and 2013 in the subbasins analyzed, streamflow volumes have not decreased in the Guadalupe River Basin. Within the Neches River Basin, results of precipitation trend analyses on an annual time step indicated upward trends in the basin. None of the data from stations analyzed in the Neches River Basin indicated annual trends in streamflow despite upward trends in annual precipitation within the basin. Data from 9 of the 19 stations analyzed in the basin indicated upward trends in annual minimum streamflow. Data from one of the simulated-inflow stations indicated a downward trend in annual minimum streamflow into Sam Rayburn Reservoir. Data from two stations indicated downward trends in annual peak streamflow, and data from one small subbasin indicated an upward trend in annual peak streamflow. Despite the known addition of 4,839,609 acre-feet of potential flood storage between 1888 and 2008 in the subbasins analyzed, there have not been widespread reductions in streamflow volumes in the Neches River Basin. Within the Sulphur River Basin, results of precipitation trend analyses on an annual time step indicated a moderate upward trend within the basin. Data from only one of the stations, the simulated inflow to Jim Chapman Lake, indicated an annual upward trend in streamflow despite an upward trend in annual precipitation throughout the basin. Data from three of the six stations in the Sulphur River Basin indicated upward trends in annual minimum streamflow, and data from one of the six stations indicated a downward trend in annual peak streamflow. Despite the known addition of 6,933,361 acre-feet of potential flood storage between 1904 and 2006 in the subbasins analyzed, streamflow volumes have not decreased in the Sulphur River Basin. Within the Trinity River Basin, results of precipitation trend analyses on an annual time step indicated upward trends in most sections of the basin. Data from 8 of the 36 stations analyzed for trends in annual streamflow indicated upward trends, and all 8 stations are in the upper sections of the basin. None of the data from stations in the lower sections of the basin indicated trends in annual streamflow. Data from 16 of the 36 stations indicated upward trends in annual minimum streamflow. Upward trends in annual minimum streamflow could be the result of managed reservoir releases in combination with wastewater treatment plant releases in the large Dallas-Fort Worth metroplex in the upper sections of the basin. All the trends in annual peak streamflow were in the sections of the basin that include the Dallas-Fort Worth metroplex. Data from two stations, one USGS streamflow-gaging station and one simulated-inflow station, indicated upward trends in annual peak streamflow, and data from one streamflow-gaging station indicated a downward trend in annual peak streamflow. Of the basins included in this study, the Trinity River Basin has the second largest amount of potential flood storage of 8,947,349 acre-feet from dams added between 1890 and 2013. Eleven stations in the Trinity River Basin had positive associations between potential flood storage volume and annual streamflow volume, indicating that annual streamflow increases as potential flood storage increases. Data from 7 of the 11 stations also indicated upward trends in annual streamflow. The positive associations may be the result of increases in minimum streamflow, which could be the result of any combination of managed reservoir releases, wastewater treatment plant releases, or increased runoff from urbanized areas, particularly in the urbanized area of the Dallas-Fort Worth metroplex.

Texas↗

USGS2021

Hydrogeology and model-simulated groundwater availability in the Salt Fork Red River aquifer, southwestern Oklahoma, 1980–2015

The 1973 Oklahoma Water Law (82 OK Stat § 82-1020.5) requires that the Oklahoma Water Resources Board (OWRB) conduct hydrologic investigations of the State’s groundwater basins to support a determination of the maximum annual yield for each groundwater basin (hereinafter referred to as an “aquifer”). The maximum annual yield allocated per acre of land is known as the equal-proportionate-share (EPS) pumping rate. At present (2021), the OWRB has not yet established a maximum annual yield and EPS pumping rate for the Salt Fork Red River aquifer. To provide updated information to the OWRB that could support evaluation and determination of an appropriate maximum annual yield, the U.S. Geological Survey (USGS), in cooperation with the OWRB, conducted a hydrologic investigation and evaluated the effects of potential groundwater withdrawals on groundwater availability in the Salt Fork Red River aquifer. The Salt Fork Red River aquifer in Greer, Harmon, and Jackson Counties of southwestern Oklahoma is composed of about 274.5 square miles of alluvium and terrace deposits associated with the Salt Fork Red River. The mean annual recharge rate to the Salt Fork Red River aquifer for the period 1980–2015 was estimated to be about 2.94 inches per year, or 10.0 percent of the mean annual precipitation for the same period (29.4 inches per year). This 1980–2015 mean annual recharge rate is equivalent to a mean annual recharge rate of about 38,000 acre-feet per year (acre-ft/yr) for the Salt Fork Red River aquifer excluding about 19,764 acres comprising the Mulberry Creek and Horse Creek terraces. The mean annual recharge rates upgradient and downgradient from USGS streamgage 07300500 Salt Fork Red River at Mangum, Okla. (hereinafter referred to as the “Mangum gage”), apportioned by aquifer area (41.5 and 58.5 percent, respectively), were about 16,000 and 22,000 acre-ft/yr, respectively. Mean annual groundwater use for the study period (1980–2015) was 3,532.7 acre-ft/yr; about 77 percent of that groundwater use was for irrigation, and about 23 percent was for public supply. Most groundwater use for irrigation was associated with wells in the Martha terrace. A hydrogeologic framework was developed for the Salt Fork Red River aquifer and included a definition of the aquifer extent and potentiometric surface, as well as a description of the textural and hydraulic properties of aquifer materials. The hydrogeologic framework was used in the construction of the numerical groundwater-flow model of the Salt Fork Red River aquifer described in this report. A conceptual model for the Salt Fork Red River aquifer that reasonably represents the groundwater-flow system was developed to constrain the construction and calibration of the numerical model. The conceptual-model water budget estimated mean annual inflows to, and outflows from, the Salt Fork Red River aquifer for the period 1980–2015 and included a subaccounting of mean annual inflows and outflows for the portions of the aquifer that were upgradient and downgradient from the Mangum gage. The numerical groundwater-flow model of the Salt Fork Red River aquifer was constructed by using MODFLOW-2005 with the Newton formulation solver. The model of the Salt Fork Red River aquifer was spatially discretized into 1,050 rows, 1,125 columns, about 170,000 active cells measuring 200 by 200 feet (ft), and a single convertible layer. The model was temporally discretized into 432 monthly transient stress periods (each with two time steps to improve model stability). An initial steady-state stress period represented mean annual inflows to, and outflows from, the aquifer and produced a solution that was used as the initial condition for subsequent transient stress periods as well as some groundwater-availability scenarios. The model was calibrated to water-table-altitude observations at selected wells and base-flow observations at selected streamgages. The simulated saturated thickness of the Salt Fork Red River aquifer was determined by subtracting the altitude of the aquifer base from the simulated water-table altitude at the end of the numerical-model period (2015). The simulated saturated thickness was more than 75 ft in a paleochannel in the Dodson terrace near the Texas border. The mean aquifer thickness (sum of saturated and unsaturated) was 49.62 ft, and the mean saturated thickness was 28.55 ft. A simulated mean transmissivity of 1,024 feet squared per day was computed from the calibrated hydraulic conductivity and saturated thickness of each cell. The simulated available water in storage at the end of the numerical-model period (2015) was 526,117 acre-feet (acre-ft); about 42 percent of that total was available upgradient from the Mangum gage, and about 58 percent of that total was available downgradient from the Mangum gage (including the Mangum terrace). Three types of groundwater-availability scenarios were run using the calibrated numerical model. These scenarios were used to (1) estimate the EPS pumping rate that ensures a minimum 20-, 40-, and 50-year life of the aquifer, (2) quantify the potential effects of projected well withdrawals on groundwater storage over a 50-year period, and (3) simulate the potential effects of a hypothetical 10-year drought on base flow and groundwater storage. The 20-, 40-, and 50-year EPS pumping rates under normal recharge conditions were about 0.51, 0.48, and 0.48 acre-foot per acre per year, respectively. Given the 155,929-acre modeled aquifer area, these rates correspond to annual yields of about 78,800, 74,900, and 74,700 acre-ft/yr, respectively. For the 20-year EPS scenario, decreasing and increasing recharge by 10 percent resulted in a 6-percent change in the EPS pumping rate in both cases; for the 40- and 50-year EPS scenarios, decreasing and increasing recharge by 10 percent resulted in a 7-percent change in the EPS pumping rate in both cases. Projected 50-year pumping scenarios were used to simulate the effects of selected well withdrawal rates on groundwater storage of the Salt Fork Red River aquifer and base flows in the Salt Fork Red River. The effects of well withdrawals were evaluated by quantifying differences in groundwater storage and base flow in four 50-year scenarios, which applied (1) no groundwater pumping, (2) mean pumping rates for the study period (1980–2015), (3) 2015 pumping rates, and (4) increasing demand pumping rates at simulated wells. The increasing demand pumping rates assumed a cumulative 20.4-percent increase in pumping over 50 years based on 2010–60 demand projections for southwestern Oklahoma. Groundwater storage after 50 years with no pumping was 535,000 acre-ft, or 8,900 acre-ft (1.7 percent) greater than the initial groundwater storage; this groundwater storage increase is equivalent to a mean water-table-altitude increase of 0.48 ft. Groundwater storage after 50 years of pumping at the mean rate for the study period (1980–2015) was 519,900 acre-ft, or 6,200 acre-ft (1.2 percent) less than the initial groundwater storage; this groundwater storage decrease is equivalent to a mean water-table-altitude decline of 0.34 ft. Groundwater storage at the end of the 50-year period with 2015 pumping rates was 513,100 acre-ft, or 13,000 acre-ft (2.5 percent) less than the initial storage; this groundwater storage decrease is equivalent to a mean water-table-altitude decline of 0.71 ft. Groundwater storage at the end of the 50-year period with increasing demand pumping rates was 509,700 acre-ft, or 16,500 acre-ft (3.1 percent) less than the initial storage; this groundwater storage decrease is equivalent to a mean water-table-altitude decline of 0.89 ft. A hypothetical 10-year drought scenario was used to simulate the effects of a prolonged period of reduced recharge on groundwater storage. The period January 1983–December 1992 was chosen as the simulated drought period. Drought effects were quantified by comparing the results of the drought scenario to those of the calibrated numerical model (no drought) at the end of the simulated drought period (1992). To simulate the hypothetical drought, recharge in the calibrated numerical model was reduced by 50 percent during the simulated drought period (1983–92). Upstream inflows from the Salt Fork Red River, Turkey Creek, and Bitter Creek were reduced by 75 percent. Groundwater storage at the end of the drought period (1992) was 479,200 acre-ft, or 53,200 acre-ft (10.0 percent) less than the groundwater storage of the calibrated numerical model at the end of the drought period. This decrease in groundwater storage is equivalent to a mean water-table-altitude decline of 2.9 ft. At the end of the 10-year hypothetical drought period, simulated base flows at the Mangum gage and USGS streamgage 07301110 Salt Fork Red River near Elmer, Okla., had decreased by about 80 and 70 percent, respectively.

Oklahoma↗

USGS2015

Water Quality, Cyanobacteria, and Environmental Factors and Their Relations to Microcystin Concentrations for Use in Predictive Models at Ohio Lake Erie and Inland Lake Recreational Sites, 2013-14

Harmful cyanobacterial “algal” blooms (cyanoHABs) and associated toxins, such as microcystin, are a major water-quality issue for Lake Erie and inland lakes in Ohio. Predicting when and where a bloom may occur is important to protect the public that uses and consumes a water resource; however, predictions are complicated and likely site specific because of the many factors affecting toxin production. Monitoring for a variety of environmental and water-quality factors, for concentrations of cyanobacteria by molecular methods, and for algal pigments such as chlorophyll and phycocyanin by using optical sensors may provide data that can be used to predict the occurrence of cyanoHABs. To test these monitoring approaches, water-quality samples were collected at Ohio recreational sites during May–November in 2013 and 2014. In 2013, samples were collected monthly at eight sites at eight lakes to facilitate an initial assessment and select sites for more intensive sampling during 2014. In 2014, samples were collected approximately weekly at five sites at three lakes. Physical water-quality parameters were measured at the time of sampling. Composite samples were preserved and analyzed for dissolved and total nutrients, toxins, phytoplankton abundance and biovolume, and cyanobacterial genes by molecular methods. Molecular assays were done to enumerate (1) general cyanobacteria, (2) general Microcystis and Dolichospermum ( Anabaena ), (3) mcyE genes for Microcystis , Dolichospermum ( Anabaena ), and Planktothrix targeting deoxyribonucleic acid (DNA), and (4) mcyE transcripts for Microcystis , Dolichospermum ( Anabaena ), and Planktothrix targeting ribonucleic acid (RNA).The DNA assays for the mcyE gene provide data on cyanobacteria that have the potential to produce microcystin, whereas the RNA assays provide data on cyanobacteria that are actively transcribing the toxin gene. Environmental data were obtained from available online sources. Quality-control (QC) samples were collected and analyzed for all constituents to characterize bias and variability; however, QC data for molecular assays were examined in more detail than for the other constituents. The QC data for molecular assays suggested that sampling variability and qPCR variability were small in comparison with the combined variability associated with sample filtering, extraction and purification, and the matrix itself. A total of 46 water-quality samples were collected during 2013 at 8 beach sites—Buck Creek, Buckeye Crystal, Deer Creek, Harsha Main, Maumee Bay State Park (MBSP) Inland (negative control site), MBSP Lake Erie, Port Clinton, and Sandusky Bay. Microcystin was detected in 67–100 percent of samples at all sites except for MBSP Inland, where microcystin was detected in only 20 percent of samples. Microcystin concentrations ranged from <0.10 to 48 micrograms per liter (µ/L), with the widest range found at MBSP Lake Erie and the highest concentrations found at Buckeye Crystal. Saxitoxin was detected in five samples, and cylindrospermopsin was not detected in any samples. A total of 65 water-quality samples were collected during 2014 at 5 sites on 3 lakes—Buckeye Fairfield and Onion Island, Harsha Main and Campers, and MBSP Lake Erie beach. Four of the sites were bathing beaches and one site, Onion Island, was an offshore boater swim area. Concentrations of microcystin ranged from <0.10 to 240 µ/L and, as in 2013, the widest range was found at MBSP Lake Erie. At Buckeye Lake, microcystin concentrations were consistently high (greater than 20 µ/L), ranging from 23 to 81 µ/L. At Harsha Main and Campers, microcystin concentrations ranged from <0.10 to 15 µ/L. Saxitoxin was detected in four samples collected at MBSP Lake Erie. Throughout the 2014 season, the cyanobacterial community, as determined by molecular and microscopy methods, and the dominance associated with the highest microcystin concentrations were unique to individual lakes. At Buckeye Lake, Planktothrix dominated the cyanobacterial community throughout the season and Planktothrix DNA and RNA were found in 100 percent of samples; Microcystis mcyE DNA was found in low concentrations. At Harsha Lake, Dolichospermum and Microcystis were a substantial percentage of the community from late May through August, and the highest microcystin concentrations occurred in June and July. At MBSP Lake Erie, Microcystis generally dominated from mid-July through early November, and the highest microcystin concentrations occurred in August. Spearman’s correlation coefficient (rho) was computed to determine the relations between environmental and water-quality factors and microcystin concentrations at four sites—Buckeye Fairfield, Buckeye Onion Island, Harsha Main, and MBSP Lake Erie. Factors were evaluated for use as potential independent variables in two types of predictive models—daily and long-term models. Easily or continuously measured water-quality factors and available environmental data are used for daily predictions that do not require a site visit. Data from factors used in daily predictions and results from samples collected and analyzed in a laboratory are used for long-term predictions (a few days to several weeks). A few statistically significant correlations (p ≤ 0.05) between microcystin concentrations and factors for both daily and long-term predictions were found at Buckeye Onion Island, and many were found at Harsha Main and MBSP Lake Erie. There were only a few statistically significant factors for daily predictions at Buckeye Fairfield, likely because of the lack of variability in microcystin concentrations. Among factors for daily predictions, phycocyanin had the highest Spearman’s correlation to microcystin concentrations (rho = 0.79 to 0.93) at all sites except for Buckeye Fairfield. Turbidity, pH, algae category, and Secchi depth were significantly correlated to microcystin concentrations at Harsha Main and MBSP Lake Erie. Algae categories were observational categories from 0 (none) to 4 (extreme). Several discharge variables (Maumee River at Waterville, river mouth is approximately 3.5 miles from the beach) at MBSP Lake Erie were promising environmental factors for daily predictions. In addition to discrete water-quality measurements recorded at Harsha Main at the time of sampling, many manipulated measurements (factors derived from mathematical manipulation of time-series data) available from a nearby continuous monitor were strongly correlated to microcystin concentrations; the highest correlation was found for the relation between microcystin concentrations and the antecedent 7-day average phycocyanin (rho = 0.98). For long-term predictions, the most highly correlated molecular assays were Planktothrix mcyE DNA at Buckeye Onion Island and Microcystis mcyE DNA at Harsha Main and MBSP Lake Erie. Concentrations of several nutrient constituents were significantly correlated to microcystin concentrations including total nitrogen at Buckeye Onion Island, ammonia and nitrate plus nitrite (both negatively correlated) at Harsha Main and MBSP Lake Erie, and total phosphorus at MBSP Lake Erie. The results of this study showed that water-quality and environmental variables are promising for use in site-specific daily or long-term predictive models. In order to develop more accurate models to predict toxin concentrations at freshwater lake sites, data need to be collected more frequently and for consecutive days in future studies.

Ohio↗

USGS2019

Spatial and temporal variability of harmful algal blooms in Milford Lake, Kansas, May through November 2016

The U.S. Geological Survey, in cooperation with the Kansas Department of Health and Environment (KDHE), completed a study to quantify the spatial and temporal variability of cyanobacterial blooms in Milford Lake, Kansas, over a range of environmental conditions at various time scales (hours to months). A better understanding of the spatial and temporal variability of cyanobacteria and microcystin will inform sampling and management strategies for Milford Lake and for other lakes with cyanobacterial harmful algal bloom (CyanoHAB) issues throughout the Nation. Spatial and temporal variability were assessed in the upstream one-third of Milford Lake (designated as “Zone C” by KDHE) during May through November 2016 using a combination of time-lapse photography, continuous water-quality monitors, discrete phytoplankton, chlorophyll, and microcystin samples, and spatially dense near-surface data. Combined, these data were used to characterize variability of cyanobacterial abundance, algal biomass, and microcystin concentrations in Zone C of Milford Lake before, during, and after cyanobacterial blooms in 2016. Temporal patterns were evaluated during May through November 2016 using time-lapse photography at six locations in Zone C and at a single point location (the Wakefield site) using a combination of discrete and continuously measured water-quality data (including the cyanobacterial pigment phycocyanin). Based on time-lapse photography, CyanoHABs developed in Zone C of Milford Lake in early July and persisted through the end of November. Bloom accumulations at individual sites were dependent on wind direction. After a change in wind direction, it would take about 1 day for accumulations to become visible at different locations. During periods with low wind, accumulations were widespread and visible at all sites. Cyanobacteria were absent from the algal community at the Wakefield site in late May and were a minor component of the community in June; however, by mid-July the cyanobacteria were dominant and remained dominant until early November. Chlorophyll and microcystin concentrations at the Wakefield site were estimated using sensor-measured phycocyanin based on regression models developed for Zone C. Regression-estimated concentrations likely are more indicative of seasonal patterns in algal biomass (as indicated by chlorophyll concentrations) and microcystin than discretely collected samples because regression-estimated data have a much higher temporal resolution. Based on regression estimates, algal biomass and microcystin concentrations at the Wakefield site steadily increased from May through August. After August, concentrations decreased but remained relatively high compared to May and June. Daily chlorophyll maxima were as much as 400 times higher than daily minima, and daily microcystin maxima were as many as several orders of magnitude higher than daily minima. The extreme variability in algal biomass and microcystin concentrations at the Wakefield site reflects the development and dissipation of blooms, as indicated by the time-lapse cameras. Based on regression-estimated microcystin concentrations, the KDHE watch and warning thresholds for microcystin were exceeded during mid-June through late November. Exceedance of KDHE advisory thresholds often changed from no advisory to watch or warning over the course of the day because of the variability in algal biomass and microcystin concentrations caused by bloom development and dissipation. Continuous water-quality monitors may be useful in informing public-health decisions in lakes with variable CyanoHAB conditions; however, site-specific models need to be developed, and best practices for using continuous water-quality monitors to inform CyanoHAB management strategies need to be established. Spatial data were collected on May 26, July 21, and September 15, 2016, using a combination of a boat-mounted array and discrete water-quality samples analyzed for phytoplankton community composition and chlorophyll and microcystin concentrations. Spatial patterns were described using regression-estimated chlorophyll and microcystin concentrations. During the May 26, 2016, spatial surveys, cyanobacterial abundances were relatively low throughout Zone C and did not exceed KDHE guidance values compared to spatial surveys on July 21 and September 15. Regression-estimated chlorophyll concentrations were indicative of higher algal biomass uplake in Zone C, and decreases in the downlake direction towards Zone B. Regression-estimated chlorophyll concentrations also were more variable uplake than downlake. Based on regression estimates, microcystin concentrations did not exceed KDHE guidance values anywhere in Zone C on May 26. Spatial patterns in microcystin throughout Zone C did not match patterns in regression-estimated chlorophyll concentrations, likely because the algal community was not dominated by cyanobacteria at most locations in May. During the July 21, 2016, spatial surveys, cyanobacterial abundances in Zone C exceeded KDHE guidance values in 50 percent of samples. The algal community in Zone C was dominated by cyanobacteria at all locations except two, where cyanobacteria codominated with diatoms. Both locations where cyanobacteria and diatoms codominated were north of the causeway. Regression-estimated chlorophyll concentrations were indicative of higher algal biomass north of the causeway and on the eastern shore of Zone C. On July 21, algal biomass did not always decrease in the downlake direction. There was a decrease just south of the causeway but an increase shortly after with higher concentrations into Zone B. Spatial maps indicated changes in algal distribution at a 0.5-meter depth, with algae moving to the central part of the lake north of the causeway and along the eastern shore south of the causeway. Most regression-estimated microcystin concentrations on July 21 exceeded KDHE guidance values, reflecting the pervasive bloom conditions in Zone C during this period. Spatial patterns in regression-estimated microcystin concentrations throughout Zone C were similar to patterns seen in discrete samples and regression-estimated chlorophyll concentrations, with higher concentrations north of the causeway and on the east shore of Zone C. During the September 15, 2016, spatial surveys, cyanobacterial abundances did not exceed KDHE guidance values. The algal community north of the causeway was dominated by diatoms. The algal community throughout the rest of Zone C was dominated by cyanobacteria. Of regression-estimated microcystin concentrations on September 15, 80 percent did not exceed KDHE guidance values. Spatial patterns indicated northward movement of the cyanobacterial bloom consistent with a wind shift noted the previous day. On September 14, winds were generally from the north to northwest, shifting to the south by September 15. There was a northward progression of chlorophyll and microcystin during the spatial surveys. These data, along with the camera data and spatial and wind data from May and July, indicate that wind can be a major driver of the spatial and temporal variability of cyanobacterial blooms in Milford Lake and likely plays a role in the extent and duration of near-shore accumulations.

Kansas↗

USGS2016

Water-quality trends and constituent-transport analysis for selected sampling sites in the Milltown Reservoir/Clark Fork River Superfund Site in the upper Clark Fork Basin, Montana, water years 1996–2015

During the extended history of mining in the upper Clark Fork Basin in Montana, large amounts of waste materials enriched with metallic contaminants (cadmium, copper, lead, and zinc) and the metalloid trace element arsenic were generated from mining operations near Butte and milling and smelting operations near Anaconda. Extensive deposition of mining wastes in the Silver Bow Creek and Clark Fork channels and flood plains had substantial effects on water quality. Federal Superfund remediation activities in the upper Clark Fork Basin began in 1983 and have included substantial remediation near Butte and removal of the former Milltown Dam near Missoula. To aid in evaluating the effects of remediation activities on water quality, the U.S. Geological Survey began collecting streamflow and water-quality data in the upper Clark Fork Basin in the 1980s. Trend analysis was done on specific conductance, selected trace elements (arsenic, copper, and zinc), and suspended sediment for seven sampling sites in the Milltown Reservoir/Clark Fork River Superfund Site for water years 1996–2015. The most upstream site included in trend analysis is Silver Bow Creek at Warm Springs, Montana (sampling site 8), and the most downstream site is Clark Fork above Missoula, Montana (sampling site 22), which is just downstream from the former Milltown Dam. Water year is the 12-month period from October 1 through September 30 and is designated by the year in which it ends. Trend analysis was done by using a joint time-series model for concentration and streamflow. To provide temporal resolution of changes in water quality, trend analysis was conducted for four sequential 5-year periods: period 1 (water years 1996–2000), period 2 (water years 2001–5), period 3 (water years 2006–10), and period 4 (water years 2011–15). Because of the substantial effect of the intentional breach of Milltown Dam on March 28, 2008, period 3 was subdivided into period 3A (October 1, 2005–March 27, 2008) and period 3B (March 28, 2008–September 30, 2010) for the Clark Fork above Missoula (sampling site 22). Trend results were considered statistically significant when the statistical probability level was less than 0.01. In conjunction with the trend analysis, estimated normalized constituent loads (hereinafter referred to as “loads”) were calculated and presented within the framework of a constituent-transport analysis to assess the temporal trends in flow-adjusted concentrations (FACs) in the context of sources and transport. The transport analysis allows assessment of temporal changes in relative contributions from upstream source areas to loads transported past each reach outflow. Trend results indicate that FACs of unfiltered-recoverable copper decreased at the sampling sites from the start of period 1 through the end of period 4; the decreases ranged from large for one sampling site (Silver Bow Creek at Warm Springs [sampling site 8]) to moderate for two sampling sites (Clark Fork near Galen, Montana [sampling site 11] and Clark Fork above Missoula [sampling site 22]) to small for four sampling sites (Clark Fork at Deer Lodge, Montana [sampling site 14], Clark Fork at Goldcreek, Montana [sampling site 16], Clark Fork near Drummond, Montana [sampling site 18], and Clark Fork at Turah Bridge near Bonner, Montana [sampling site 20]). For period 4 (water years 2011–15), the most notable changes indicated for the Milltown Reservoir/Clark Fork River Superfund Site were statistically significant decreases in FACs and loads of unfiltered-recoverable copper for sampling sites 8 and 22. The period 4 changes in FACs of unfiltered-recoverable copper for all other sampling sites were not statistically significant. Trend results indicate that FACs of unfiltered-recoverable arsenic decreased at the sampling sites from period 1 through period 4 (water years 1996–2015); the decreases ranged from minor (sampling sites 8–20) to small (sampling site 22). For period 4 (water years 2011–15), the most notable changes indicated for the Milltown Reservoir/Clark Fork River Superfund Site were statistically significant decreases in FACs and loads of unfiltered-recoverable arsenic for sampling site 8 and near statistically significant decreases for sampling site 22. The period 4 changes in FACs of unfiltered-recoverable arsenic for all other sampling sites were not statistically significant. Trend results indicate that FACs of suspended sediment decreased at the sampling sites from period 1 through period 4 (water years 1996–2015); the decreases ranged from moderate (sampling site 8) to small (sampling sites 11–22). For period 4 (water years 2011–15), the changes in FACs of suspended sediment were not statistically significant for any sampling sites. The reach of the Clark Fork from Galen to Deer Lodge is a large source of metallic contaminants and suspended sediment, which strongly affects downstream transport of those constituents. Mobilization of copper and suspended sediment from flood-plain tailings and the streambed of the Clark Fork and its tributaries within the reach results in a contribution of those constituents that is proportionally much larger than the contribution of streamflow from within the reach. Within the reach from Galen to Deer Lodge, unfiltered-recoverable copper loads increased by a factor of about 4 and suspended-sediment loads increased by a factor of about 5, whereas streamflow increased by a factor of slightly less than 2. For period 4 (water years 2011–15), unfiltered-recoverable copper and suspended-sediment loads sourced from within the reach accounted for about 41 and 14 percent, respectively, of the loads at Clark Fork above Missoula (sampling site 22), whereas streamflow sourced from within the reach accounted for about 4 percent of the streamflow at sampling site 22. During water years 1996–2015, decreases in FACs and loads of unfiltered-recoverable copper and suspended sediment for the reach generally were proportionally smaller than for most other reaches. Unfiltered-recoverable copper loads sourced within the reaches of the Clark Fork between Deer Lodge and Turah Bridge near Bonner (just upstream from the former Milltown Dam) were proportionally smaller than contributions of streamflow sourced from within the reaches; these reaches contributed proportionally much less to copper loading in the Clark Fork than the reach between Galen and Deer Lodge. Although substantial decreases in FACs and loads of unfiltered-recoverable copper and suspended sediment were indicated for Silver Bow Creek at Warm Springs (sampling site 8), those substantial decreases were not translated to downstream reaches between Deer Lodge and Turah Bridge near Bonner. The effect of the reach of the Clark Fork from Galen to Deer Lodge as a large source of copper and suspended sediment, in combination with little temporal change in those constituents for the reach, contributes to this pattern. With the removal of the former Milltown Dam in 2008, substantial amounts of contaminated sediments that remained in the Clark Fork channel and flood plain in reach 9 (downstream from Turah Bridge near Bonner) became more available for mobilization and transport than before the dam removal. After the removal of the former Milltown Dam, the Clark Fork above Missoula (sampling site 22) had statistically significant decreases in FACs of unfiltered-recoverable copper in period 3B (March 28, 2008, through water year 2010) that continued in period 4 (water years 2011–15). Also, decreases in FACs of unfiltered-recoverable arsenic and suspended sediment were indicated for period 4 at this site. The decrease in FACs of unfiltered-recoverable copper for sampling site 22 during period 4 was proportionally much larger than the decrease for the Clark Fork at Turah Bridge near Bonner (sampling site 20). Net mobilization of unfiltered-recoverable copper and arsenic from sources within reach 9 are smaller for period 4 than for period 1 when the former Milltown Dam was in place, providing evidence that contaminant source materials have been substantially reduced in reach 9.

Montana↗

USGS2012

Water quality and the composition of fish and macroinvertebrate communities in the Devils and Pecos Rivers within and upstream from the Amistad National Recreation Area, Texas, 2005-7

To gain a better understanding of the water quality and status of fish and macroinvertebrate communities, the U.S. Geological Survey, in cooperation with the National Park Service and Amistad National Recreation Area, completed a reconnaissance-level survey of the water quality and fish and macroinvertebrate communities in the Devils and Pecos Rivers in and upstream from the Amistad National Recreation Area in southwest Texas during 2005–7. Water-quality conditions during the spring and summer months of 2005 in the Devils and Pecos Rivers were assessed at locations just upstream from the Amistad National Recreation Area, and the composition of fish and macroinvertebrate communities were assessed during 2006 and 2007 in and upstream from the Amistad National Recreation Area and Amistad Reservoir. Water-quality samples were collected at one site on both the Devils and Pecos Rivers. Fish and macroinvertebrates were collected at the water-quality sampling site on each river and at three additional sites on each river. The water-quality constituents of primary concern were total dissolved solids, chloride, sulfate, ammonia plus organic nitrogen, nitrate plus nitrite, orthophosphate, phosphorus, selenium, and selected pesticides. During the spring and summer of 2005, the concentrations of total dissolved solids ranged from 208 to 232 milligrams per liter (mg/L) in samples from the Devils River compared to 1,460 to 2,390 mg/L in samples from the Pecos River. Total dissolved solid concentrations measured in samples collected from the Devils River and Pecos River did not exceed the proposed State of Texas water-quality standard applicable for the segments of each river where samples were collected. During the spring and summer of 2005, chloride concentrations measured in samples collected in 2005 from the Devils River ranged from 11.6 to 12.9 mg/L, compared to chloride concentrations measured in samples collected from the Pecos River, which ranged from 519 to 879 mg/L. Chloride concentrations in samples collected from the Devils River in 2005 were less than the lower quartile (25th percentile) value of 14.0 mg/L reported for chloride concentrations in water-quality samples collected at the same sampling location during 1978–95 by the U.S. Geological Survey as part of the Hydrologic Benchmark Network program. The chloride concentrations measured in samples collected from the Pecos River during the spring and summer of 2005 represented a range of values similar to the interquartile range of 548 to 942 mg/L reported for samples collected during 1974–2007 at the same sampling location by the U.S. Geological Survey as part of the National Stream Quality Accounting Network program. None of the chloride concentrations measured in samples collected from the Devils or Pecos Rivers in 2005 exceeded applicable proposed State of Texas water-quality standards for chloride. Sulfate concentrations ranged from 7.55 to 8.20 mg/L in samples from the Devils River compared to 298 to 503 mg/L in samples from the Pecos River. Concentrations of sulfate did not exceed applicable proposed State of Texas water-quality standards. Ammonia plus organic nitrogen concentrations were reported as nitrogen ranged from 0.12 to 0.14 mg/L of nitrogen in samples collected from the Devils River compared to 0.15 to 0.32 mg/L of nitrogen in samples collected from the Pecos River. Ammonia plus organic nitrogen concentrations measured in samples collected from the Devils River in 2005 were less than the lower quartile (25th percentile) value of 0.23 mg/L of nitrogen for concentrations in water-quality samples collected at the same sampling location during 1978–95 by the U.S. Geological Survey as part of the Hydrologic Benchmark Network program. Ammonia plus nitrogen concentrations measured in samples collected from the Pecos River were similar to the range of historical concentrations measured in samples collected from the same Pecos River sampling location by the U.S. Geological Survey National Stream Quality Accounting Network program. Nitrate plus nitrite concentrations in samples from the Devils River and Pecos Rivers were within the historical range of concentrations for samples collected at the same locations on each river. Total phosphorous and orthophosphate concentrations were less than the laboratory reporting levels in the water samples from the Devils and Pecos Rivers. None of the selenium concentrations measured in samples collected during the spring and summer of 2005 from the Devils or Pecos Rivers exceeded the Texas Surface Water Quality Standards (chronic criterion of 5 μg/L or the acute criterion of 20 μg/L) established by the State for the protection of aquatic life. Concentrations of pesticides in the samples collected from the Devils and Pecos Rivers during March–August 2005 were very low and not present in detectable amounts (all reported concentrations were below laboratory reporting levels). The total number of fish species collected was the same in the Devils River and Pecos River, but the species found in the two rivers varied slightly. The number of fish species generally increased from the site farthest upstream to the site farthest downstream in the Devils River, and decreased between the site farthest upstream and site farthest downstream in the Pecos River. The redbreast sunfish was the most abundant species collected in the Devils River, and the blacktail shiner was the most abundant species collected in the Pecos River. Comparing the species from each river, the percentage of omnivorous fish species was larger at the more downstream sites closer to Amistad Reservoir, and the percentage of species tolerant of environmental stressors was larger in the Pecos River. The fish community, assessed on the basis of the number of shared species among the sites sampled, was more similar to the fish community at the other sites on the same river than it was to the fish community from any other site in the other river. More macroinvertebrate taxa were collected in the Devils River than in the Pecos River. The largest number of macroinvertebrate taxa were from the site second farthest downstream on the Devils River, and the smallest numbers of macroinvertebrate taxa were from the farthest downstream site on the Pecos River. Mayflies were more common in the Devils River, and caddisflies were less common than mayflies at most sites. Net-spinning caddisflies were more common at the Devils River sites. The combined percent of mayfly, caddisfly, and stonefly taxa was generally larger at the Pecos River sites. Riffle beetles were the most commonly collected beetle taxon among all sites, and water-penny beetles were only collected at the Pecos River sites. A greater number of true midge taxa were collected more than any other taxa at the genus and species taxonomic level. Non-insect macroinvertebrate taxa were more common at the Devils River sites. Corbicula sp. (presumably the introduced Asian clam) was found at sites in both rivers, and amphipods were more abundant in the Devils River. The Margalef species richness index, based on aquatic insect taxa only, was larger at the Devils River sites than at the Pecos River sites. The Hilsenhoff's biotic index was largest at the site farthest downstream in the Devils River and smallest at the site second farthest downstream in the Pecos River. Overall similarity among sites based on the number of shared macroinvertebrate taxa indicated that each site is more similar to other sites on the same river than to sites on the other river.

Texas↗