Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Scientific Monthly”

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

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

At least 775 records · Page 43Linked to original sources
USGS2012

Computing daily mean streamflow at ungaged locations in Iowa by using the Flow Anywhere and Flow Duration Curve Transfer statistical methods

The U.S. Geological Survey (USGS) maintains approximately 148 real-time streamgages in Iowa for which daily mean streamflow information is available, but daily mean streamflow data commonly are needed at locations where no streamgages are present. Therefore, the USGS conducted a study as part of a larger project in cooperation with the Iowa Department of Natural Resources to develop methods to estimate daily mean streamflow at locations in ungaged watersheds in Iowa by using two regression-based statistical methods. The regression equations for the statistical methods were developed from historical daily mean streamflow and basin characteristics from streamgages within the study area, which includes the entire State of Iowa and adjacent areas within a 50-mile buffer of Iowa in neighboring states. Results of this study can be used with other techniques to determine the best method for application in Iowa and can be used to produce a Web-based geographic information system tool to compute streamflow estimates automatically. The Flow Anywhere statistical method is a variation of the drainage-area-ratio method, which transfers same-day streamflow information from a reference streamgage to another location by using the daily mean streamflow at the reference streamgage and the drainage-area ratio of the two locations. The Flow Anywhere method modifies the drainage-area-ratio method in order to regionalize the equations for Iowa and determine the best reference streamgage from which to transfer same-day streamflow information to an ungaged location. Data used for the Flow Anywhere method were retrieved for 123 continuous-record streamgages located in Iowa and within a 50-mile buffer of Iowa. The final regression equations were computed by using either left-censored regression techniques with a low limit threshold set at 0.1 cubic feet per second (ft3/s) and the daily mean streamflow for the 15th day of every other month, or by using an ordinary-least-squares multiple linear regression method and the daily mean streamflow for the 15th day of every other month. The Flow Duration Curve Transfer method was used to estimate unregulated daily mean streamflow from the physical and climatic characteristics of gaged basins. For the Flow Duration Curve Transfer method, daily mean streamflow quantiles at the ungaged site were estimated with the parameter-based regression model, which results in a continuous daily flow-duration curve (the relation between exceedance probability and streamflow for each day of observed streamflow) at the ungaged site. By the use of a reference streamgage, the Flow Duration Curve Transfer is converted to a time series. Data used in the Flow Duration Curve Transfer method were retrieved for 113 continuous-record streamgages in Iowa and within a 50-mile buffer of Iowa. The final statewide regression equations for Iowa were computed by using a weighted-least-squares multiple linear regression method and were computed for the 0.01-, 0.05-, 0.10-, 0.15-, 0.20-, 0.30-, 0.40-, 0.50-, 0.60-, 0.70-, 0.80-, 0.85-, 0.90-, and 0.95-exceedance probability statistics determined from the daily mean streamflow with a reporting limit set at 0.1 ft 3 /s. The final statewide regression equation for Iowa computed by using left-censored regression techniques was computed for the 0.99-exceedance probability statistic determined from the daily mean streamflow with a low limit threshold and a reporting limit set at 0.1 ft 3 /s. For the Flow Anywhere method, results of the validation study conducted by using six streamgages show that differences between the root-mean-square error and the mean absolute error ranged from 1,016 to 138 ft 3 /s, with the larger value signifying a greater occurrence of outliers between observed and estimated streamflows. Root-mean-square-error values ranged from 1,690 to 237 ft 3 /s. Values of the percent root-mean-square error ranged from 115 percent to 26.2 percent. The logarithm (base 10) streamflow percent root-mean-square error ranged from 13.0 to 5.3 percent. Root-mean-square-error observations standard-deviation-ratio values ranged from 0.80 to 0.40. Percent-bias values ranged from 25.4 to 4.0 percent. Untransformed streamflow Nash-Sutcliffe efficiency values ranged from 0.84 to 0.35. The logarithm (base 10) streamflow Nash-Sutcliffe efficiency values ranged from 0.86 to 0.56. For the streamgage with the best agreement between observed and estimated streamflow, higher streamflows appear to be underestimated. For the streamgage with the worst agreement between observed and estimated streamflow, low flows appear to be overestimated whereas higher flows seem to be underestimated. Estimated cumulative streamflows for the period October 1, 2004, to September 30, 2009, are underestimated by -25.8 and -7.4 percent for the closest and poorest comparisons, respectively. For the Flow Duration Curve Transfer method, results of the validation study conducted by using the same six streamgages show that differences between the root-mean-square error and the mean absolute error ranged from 437 to 93.9 ft 3 /s, with the larger value signifying a greater occurrence of outliers between observed and estimated streamflows. Root-mean-square-error values ranged from 906 to 169 ft 3 /s. Values of the percent root-mean-square-error ranged from 67.0 to 25.6 percent. The logarithm (base 10) streamflow percent root-mean-square error ranged from 12.5 to 4.4 percent. Root-mean-square-error observations standard-deviation-ratio values ranged from 0.79 to 0.40. Percent-bias values ranged from 22.7 to 0.94 percent. Untransformed streamflow Nash-Sutcliffe efficiency values ranged from 0.84 to 0.38. The logarithm (base 10) streamflow Nash-Sutcliffe efficiency values ranged from 0.89 to 0.48. For the streamgage with the closest agreement between observed and estimated streamflow, there is relatively good agreement between observed and estimated streamflows. For the streamgage with the poorest agreement between observed and estimated streamflow, streamflows appear to be substantially underestimated for much of the time period. Estimated cumulative streamflow for the period October 1, 2004, to September 30, 2009, are underestimated by -9.3 and -22.7 percent for the closest and poorest comparisons, respectively.

Illinois;Iowa;Minnesota;Missouri;Nebraska;Wisconsi↗
USGS2021

Assessment of runoff volume reduction associated with soil amendments added to portions of highway median-strip catchments in Ohio, 2018–20

The U.S. Geological Survey installed 10 rain gages and 12 calibrated H-flumes to measure rainfall and runoff volumes at 10 locations in Ohio Department of Transportation highway median-strip catchments. Data were collected to facilitate comparisons of rainfall and runoff volumes at study sites before and after stormwater best management practices (BMPs) were installed and between sites with different BMPs. The BMP treatments comprised removing the top layer of the existing soil, rototilling the remaining soil to a 6-inch depth, mixing the soils with one of two soil amendments (compost with sand or shale) at one of two thicknesses (4 inches or 6 inches), topping with a compost blanket, seeding, and installing erosion control matting. The overall treatment used at a given study site is referred to as “BMP.” At two locations where soil amendments were installed, a second “control” site was installed to measure runoff from an adjacent catchment in the same median strip where no soil amendment was installed. This no-treatment option (no soil amendment) was considered its own class of BMP. Rainfall and runoff data were collected during periods when air temperatures were above freezing (including all months except January, February, and parts of December and March) from 2018 to 2020. The data collection period for each study site was divided into “pre-BMP” and “post-BMP” periods. Equipment to measure rainfall and runoff was installed and data were collected from April to December 2018 before installation of soil amendments (the pre-BMP period). The post-BMP period started between April and May of 2019 at the first measured rainfall after soil amendments were installed. Rainfall and runoff monitoring continued through September 2020. For control sites, the post-BMP periods were assigned to start with the first measured rainfall in the 2019 data collection season. A rainfall-runoff “event” was defined as beginning at the time of the first measured rainfall and ending when rainfall and runoff (if any) ceased and remained ceased for at least 3 hours. A value referred to as “event runoff percentage,” defined as the total volume of runoff during an event expressed as a percentage of the total volume of rainfall falling over the catchment, was computed for each event. The distribution of rainfall totals associated with events was similar between the pre-BMP and post-BMP periods; however, there were appreciable between-site differences in the distribution of event runoff percentages during the pre-BMP and post-BMP periods. Empirical distribution function (EDF) tests were performed with and without data from events that resulted in no runoff to determine whether the distribution of event runoff percentages changed from the pre-BMP period to the post-BMP period. The null hypothesis that the EDFs of event runoff percentages were equal in the pre-BMP and post-BMP periods was rejected ( α =0.05) in at least one of the two tests for four sites (one site with a shale amendment and three sites with sand amendments). Mean event runoff percentages at each of those four sites decreased from the pre-BMP period to the post-BMP period. The null hypothesis that the EDFs of event runoff percentages were equal was not rejected for the other six sites’ draining catchments with soil amendments or the two control sites. EDF tests performed on event rainfall totals indicated no statistically significant changes between the pre-BMP and post-BMP period distributions for any of the sites. Double-mass analyses of cumulative runoff were performed for two pairs of closely spaced sites (each pair located in a common median strip): one site in each pair drained a catchment where soil amendments were installed, and the other (a control) drained a catchment without soil amendments. Those double-mass analyses indicated a small reduction in runoff from the pre-BMP to post-BMP period at the site whose catchment received the sand and compost amendment, but no perceptible reduction in runoff at the site whose catchment received the shale and compost amendment. Regression analyses indicated that (a) three rainfall factors (event rainfall totals, total rainfall for the previous 7 days, and a cross product of the factors) and the intercept term were the four most important factors explaining event runoff percentages, (b) the effect of amendment type on event runoff percentage was small in comparison to the rainfall and intercept terms, (c) event runoff percentages tended to be lower for sites with shale amendments than sites with sand amendments; however, event runoff percentages tended to be lower for control sites than for sites with shale or sand amendments, and (d) event runoff percentages increased with increasing amendment thickness. The counterintuitive results that event runoff percentages increased with increasing amendment thickness and that control sites tended to have lower event runoff percentages than sites draining soil-amended catchments likely reflects unmeasured factors that existed at the sites before BMPs were installed rather than the effect of the BMP treatments. Although not definitive, some support for the conclusion that the sand amendment was generally more effective at reducing runoff than the shale amendment was provided by results from the EDF tests, double-mass analyses, and runoff statistics.

Ohio↗
USGS2014

California State Waters Map Series: Offshore of Coal Oil Point, California

In 2007, the California Ocean Protection Council initiated the California Seafloor Mapping Program (CSMP), designed to create a comprehensive seafloor map of high-resolution bathymetry, marine benthic habitats, and geology within the 3-nautical-mile limit of California’s State Waters. The CSMP approach is to create highly detailed seafloor maps through collection, integration, interpretation, and visualization of swath sonar data, acoustic backscatter, seafloor video, seafloor photography, high-resolution seismic-reflection profiles, and bottom-sediment sampling data. The map products display seafloor morphology and character, identify potential marine benthic habitats, and illustrate both the surficial seafloor geology and shallow (to about 100 m) subsurface geology. The Offshore of Coal Oil Point map area lies within the central Santa Barbara Channel region of the Southern California Bight. This geologically complex region forms a major biogeographic transition zone, separating the cold-temperate Oregonian province north of Point Conception from the warm-temperate California province to the south. The map area is in the southern part of the Western Transverse Ranges geologic province, which is north of the California Continental Borderland. Significant clockwise rotation—at least 90°—since the early Miocene has been proposed for the Western Transverse Ranges province, and geodetic studies indicate that the region is presently undergoing north-south shortening. Uplift rates (as much as 2.0 mm/yr) that are based on studies of onland marine terraces provide further evidence of significant shortening. The cities of Goleta and Isla Vista, the main population centers in the map area, are in the western part of a contiguous urban area that extends eastward through Santa Barbara to Carpinteria. This urban area is on the south flank of the east-west-trending Santa Ynez Mountains, on coalescing alluvial fans and uplifted marine terraces underlain by folded and faulted Miocene bedrock. In the map area, the relatively low-relief, elevated coastal bajada narrows from about 2.5 km wide in the east to less than 500 m wide in the west. Several beaches line the actively utilized coastal zone, including Isla Vista County Park beach, Coal Oil Point Reserve, and Goleta Beach County Park. The beaches are subject to erosion each winter during storm-wave attack, and then they undergo gradual recovery or accretion during the more gentle wave climate of the late spring, summer, and fall months. The Offshore of Coal Oil Point map area lies in the central part of the Santa Barbara littoral cell, which is characterized by littoral drift to the east-southeast. Longshore drift rates have been reported to range from about 160,000 to 800,000 tons/yr, averaging 400,000 tons/yr. Sediment supply to the western and central parts of the littoral cell, including the map area, is largely from relatively small transverse coastal watersheds. Within the map area, these coastal watersheds include (from east to west) Las Llagas Canyon, Gato Canyon, Las Varas Canyon, Dos Pueblos Canyon, Eagle Canyon, Tecolote Canyon, Winchester Canyon, Ellwood Canyon, Glen Annie Canyon, and San Jose Creek. The Santa Ynez and Santa Maria Rivers, the mouths of which are about 100 to 140 km northwest of the map area, are not significant sediment sources because Point Conception and Point Arguello provide obstacles to downcoast sediment transport and also because much of their sediment load is trapped in dams. The Ventura and Santa Clara Rivers, the mouths of which are about 45 to 55 km southeast of the map area, are much larger sediment sources. Still farther east, eastward-moving sediment in the littoral cell is trapped by Hueneme and Mugu Canyons and then transported to the deep-water Santa Monica Basin. The offshore part of the map area consists of a relatively flat and shallow continental shelf, which dips gently seaward (about 0.8° to 1.0°) so that water depths at the shelf break, roughly coincident with the California’s State Waters limit, are about 90 m. This part of the Santa Barbara Channel is relatively well protected from large Pacific swells from the north and northwest by Point Conception and from the south and southwest by offshore islands and banks. The shelf is underlain by variable amounts of upper Quaternary marine and fluvial sediments deposited as sea level fluctuated in the late Pleistocene. The large (130 km2) Goleta landslide complex lies along the shelf break in the southern part of the map area. This compound slump complex may have been initiated more than 200,000 years ago, but it also includes three recent failures that may have been generated between 8,000 to 10,000 years ago. A local, 5- to 10-m-high tsunami may have been generated from these failure events. The map area has had a long history of hydrocarbon development, which began in 1928 with discovery of the Ellwood oil field. Subsequent discoveries in the offshore include South Ellwood offshore oil field, Coal Oil Point oil field, and Naples oil and gas field. Development of South Ellwood offshore field began in 1966 from platform “Holly,” the last platform to be installed in California’s State Waters. The area also is known for “the world’s most spectacular marine hydrocarbon seeps,” and large tar seeps are exposed on beaches east of the mouth of Goleta Slough. Offshore seeps adjacent to South Ellwood oil field release about 40 tons per day of methane and about 19 tons per day of ethane, propane, butane, and higher hydrocarbons. Seafloor habitats in the broad Santa Barbara Channel region consist of significant amounts of soft sediment and isolated areas of rocky habitat that support kelp-forest communities nearshore and rocky-reef communities in deep water. The potential marine benthic habitat types mapped in the Offshore of Coal Oil Point map area are directly related to its Quaternary geologic history, geomorphology, and active sedimentary processes. These potential habitats, which lie primarily within the Shelf (continental shelf) but also partly within the Flank (basin flank or continental slope) megahabitats, range from soft, unconsolidated sediment to hard sedimentary bedrock. This heterogeneous seafloor provides promising habitat for rockfish, groundfish, crabs, shrimp, and other marine benthic organisms.

California↗
USGS2008

Occurrence and distribution of algal biomass and Its relation to nutrients and selected basin characteristics in Indiana streams, 2001-2005

Algal biomass and nutrient data were gathered at 322 randomly selected sites on 261 streams in the West Fork White River, Whitewater River, East Fork White River, Upper Wabash River, Kankakee River, Lower Wabash River, Tributaries to the Great Lakes, and Tributaries to the Ohio River Basins in Indiana from May through October for years 2001 through 2005. Basin characteristics (land use and drainage area), substrate, turbidity, and nutrient concentrations were determined for the basin and sampling sites. The relations of the seasonal algal biomass parameters periphyton chlorophyll a (CHLa), ash-free dry mass (AFDM), seston CHLa, and particulate organic carbon (POC) to concentrations of the seasonal nutrients nitrate, total Kjeldahl nitrogen (TKN), total nitrogen (TN), and total phosphorus (TP) were determined using Spearman’s rho. The effects of streamflow were determined using data collected at U.S. Geological Survey (USGS) streamflow-gaging stations spatially located throughout the study basins. Throughout the 5-year study, the magnitude and frequency of stream discharge varied monthly and annually and greatly influenced algal biomass concentrations through scour and algal drift. Algal biomass median concentrations in Indiana streams consisted of periphyton CHLa, 41.2 milligrams per square meter (mg/m 2 ); AFDM, 52.1 grams per square meter (g/m 2 ); seston CHLa, 2.44 micrograms per liter (µg/L); and POC, 0.75 milligrams per liter (mg/L). Approximately 32 percent of the periphyton CHLa and 6 percent of the seston CHLa samples would be considered eutrophic (nutrient enriched). To ascertain seasonal variability, samples were collected in the spring (May), summer (June through August), and fall (September through October). The highest median concentration of periphyton CHLa was in the spring, 63.2 mg/m 2 , while the highest median concentrations of AFDM, seston CHLa, and POC were in the summer 55.4 g/m 2 , 2.96 µg/L, and 0.81 mg/L respectively. There were no significant differences among seasons for periphyton CHLa and AFDM; there were significant differences among seasons for seston CHLa and POC. There were no significant relations with nutrients and periphyton or seston CHLa parameters. The only significant positive relations were observed between summer POC and summer TP as well as summer POC and summer TKN. Positive relations also related spring POC and spring TP. These significant relations with TP are most likely related to phosphorus associated within seston algal cells and attached to sediment. Drainage area and land use were analyzed to understand the effect of site location on algal growth. Study basins varied in size (headwater streams, 0–51 km 2 ; wadable streams, 52–2,590 km 2 ; and boatable streams, 2,591–38,900 km 2 ) and were dominated by agricultural land use. Basin characteristics (land use, drainage area) as well as substrate type, turbidity, and nutrients, affected the concentration of algal biomass parameters. Of the eight basins in which samples were collected during the 5-year study, the Whitewater River Basin (2002) had the highest median concentration of periphyton CHLa (63.1 mg/m 2 ), the Tributaries to the Great Lakes (2005) exhibited the highest median concentration for AFDM (160 g/m 2 ), the East Fork White River Basin (2002) had nearly twice the median concentration (4.01 µg/L) of seston CHLa as the other basins, and the West Fork White River Basin (2001) exhibited the highest median concentration of POC (1.10 mg/L). Of the eight major basins sampled, 15–45 percent of the periphyton CHLa and up to 20 percent of the seston CHLa samples were eutrophic. Samples collected at headwater and wadable streams were the most eutrophic for periphyton CHLa (31–36 percent) and 28 percent of samples collected at boatable streams were eutrophic for seston CHLa. As basin size increased, seston CHLa and POC concentrations increased while periphyton CHLa and AFDM concentrations decreased. The median turbidity values ranged from 6.95 NTU for headwater streams to 8.27 NTU for wadable streams, and 17.0 NTU for boatable streams. In addition, the types and availability of periphytic substrates (epilithic, epipsammic, or epidendric) were an important factor when comparing periphyton CHLa and AFDM concentrations in the study due to the periphytic substrates individual susceptibility to bed movement and scouring. Periphyton CHLa median concentrations ranged from 53.8 mg/m 2 for epilithic substrates, to 41.8 mg/m 2 for epipsammic substrates, and 17.2 mg/m 2 for epidendric substrates. Higher AFDM concentrations were collected from epipsammic substrates during years of low stream discharge velocity, which enhanced the settling of organic matter on epipsammic substrates. AFDM median concentrations ranged from 141 g/m 2 for epipsammic to 28.8 g/m 2 for epilithic and 22.9 g/m 2 for epidendric substrates. The seasonal values for nutrients (nitrate, TKN, TN, and TP) and algal biomass (periphyton CHLa, AFDM, seston CHLa, and POC) were compared to published U. S. Environmental Protection Agency (USEPA) values for their respective ecoregions. Algal biomass values either were greater than the 25th percentile published USEPA values or extended the range of data in Aggregate Nutrient Ecoregions VI, VII, IX and USEPA Level III Ecoregions 54, 55, 56, 71, and 72. If the values for the 25th percentile proposed by the USEPA were adopted as nutrient water-quality criteria, then about 71 percent of the nutrient samples and 57 percent of the CHLa samples within the eight study basins would be considered nutrient enriched.

Scientific Investigations Report↗
USGS2018

Groundwater contributions to excessive algal growth in the East Fork Carson River, Carson Valley, west-central Nevada, 2010 and 2012

Excessive algal growth and low dissolved oxygen concentrations were observed during low streamflow conditions during summer months along a 5,800-foot reach of the East Fork Carson River in Carson Valley, west-central Nevada. Algal growth from nutrient enrichment of a stream reduces aquatic diversity, threatens fish ecology and stream health, and can be a recreational nuisance. In response to concerns that groundwater discharging to the 5,800-foot reach of the East Fork Carson River may be a source of nutrients to the stream, the U.S. Geological Survey, in cooperation with the Carson Water Subconservancy District and the Nevada Division of Environmental Protection, conducted studies during the summers of 2010 and 2012 to gain an improved understanding of the contributions of nutrients to the stream from groundwater, characterize algal conditions and algal effects on water quality, assess potential sources of nitrate in groundwater discharging to the stream, and evaluate nitrate reduction in groundwater from denitrification. A reconnaissance study in the summer of 2010 along the 5,800-foot study reach located a subreach with clear evidence of nutrient-rich groundwater discharging to the stream. At the subreach, nitrate plus nitrite (referred to hereafter as nitrate) concentrations in groundwater discharging to the stream were high (average 2.75 milligrams per liter as nitrogen) along the right bank. The stream at this location had the highest stream nitrate concentrations (average 0.056 milligrams per liter as nitrogen) compared to other locations upstream and downstream of the subreach. As a result, the 2012 study focused on a 405-foot subreach of the East Fork Carson River centered where results from the 2010 study found the highest stream and groundwater concentrations of nitrate, as well as the greatest observed contributions of groundwater discharge to the stream. Groundwater nutrient concentrations were much higher than stream nutrient concentrations during the summer of 2012 during low streamflow conditions at the 405-foot subreach of the East Fork Carson River. Average groundwater nitrate and orthophosphate concentrations along the right bank of the 405‑foot subreach were 9 and 12 times higher, respectively, than in the stream at this subreach. Groundwater discharge rates to the study reach based on different methods varied from 0.09 to 1.2 cubic feet per second per mile. Estimated groundwater discharge rates to the right bank of the study subreach were used to calculate groundwater nutrient load estimates to the subreach right bank, which were found to be low when compared to stream nutrient loads. Elevated algal biomass levels above nuisance thresholds were observed during the summers of 2010 and 2012. The study reach was characterized as mesotrophic-eutrophic during the 2010 study and eutrophic during the 2012 study. The presence of algae caused daily dissolved oxygen and pH fluctuations in the stream, resulting in exceedances of the State of Nevada water-quality standards owing to low dissolved oxygen concentrations and high pH concentrations, although the standards might not have been applicable during 2012 because of extremely low streamflow. The addition of nutrients to the stream from the constant supply in groundwater discharge sustains the growth of algae during low streamflow conditions. In the summer when streamflow is low or very low, nutrient-rich groundwater discharge enters the stream through the sediment-water interface at the streambed. Because the attached algae is thick and stream velocity is low, the nutrient-rich water pools at the sediment-water interface. Higher nutrient concentrations at the streambed create a favorable microenvironment for algae attached to the substrate to consume available nutrients from the groundwater before the groundwater mixes with overlying stream water. The source of nitrate in groundwater in this subreach is anthropogenic because nitrate concentrations are greater than background groundwater nitrate concentrations in Douglas County, high groundwater nitrate concentrations are only found at the right bank of the stream near a housing development, and organic wastewater compounds indicative of human-derived sources were also detected in groundwater wells on the right bank of the stream. Nitrogen and oxygen isotope concentrations of nitrate in shallow groundwater were used to determine the specific source of the nitrate, but the isotopic values indicated denitrification was occurring. Further investigation is needed to determine the specific anthropogenic source of the nitrate in the groundwater because the denitrification present in all samples obscures the original source of nitrogen.

Nevada↗
USGS2008

Evaluation of four structural best management practices for highway runoff in Beaufort and Colleton Counties, South Carolina, 2005–2006

From 2005 to 2006, the U.S. Geological Survey worked cooperatively with the South Carolina Department of Transportation in Beaufort and Colleton Counties, South Carolina, to assess the performance of four different structural devices that served as best management practices (BMPs). These structural devices were installed to mitigate the effects of stormwater runoff on waterways near State roads. The South Carolina Department of Transportation is required to address the quality of stormwater runoff from State-maintained roadways as part of the National Pollutant Discharge Elimination System stormwater program mandated in the Clean Water Act. The performance assessment of the four structural best management practices was based on stormflow measurements and chemical analyses of stormwater-quality samples collected during a 20-month period from March 2005 through October 2006, which represented a range of seasons and rainfall intensities. A total of 49 sample sets that included stormwater from the inlet and outlet of each of the four structural devices were collected as flow-weighted composites to provide event-mean concentrations of suspended sediment, nutrients, and trace metals. In addition, each set included grab samples that were collected to provide the first flush concentrations of oil and grease and fecal-indicator bacteria. A tiered statistical approach was used in the data analysis. Performances of the four structural BMPs were assessed individually based on how well the BMPs were able to reduce the selected constituents. Descriptive statistics and nonparametric Wilcoxon signed rank tests were applied to event-mean concentrations and loads in the water entering the inlet and in the water leaving the outlet of each BMP for each constituent to identify if significant reductions occurred. If significant reductions occurred, the BMP was considered efficient at reducing the constituent. To quantify efficiency, a simplistic approach was applied to compute mean and geometric mean efficiency ratios for the significantly reduced constituents in each BMP. Each BMP performance was ranked based on its computed efficiency ratios, however, the computed efficiency ratios were not sufficient to determine if statistical differences occurred among the performances of the four BMPs. Consequently, a more complex approach was used to apply statistical comparison tests to reduction percentages computed for individual storms (a modified removal efficiency of individual storm-load approach) to determine if differences in event-mean concentrations, loads, and reduction percentages for significantly reduced constituents occurred among the four structural BMPs. Overall, the four BMPs were efficient in reducing suspended-sediment event-mean concentrations and loads in the stormwater entering the inlets of the BMPs to significantly lower event-mean concentrations before discharging the stormwater from the outlets. The cumulative suspended-sediment event-mean load in stormwater entering the BMPs from the storms sampled during the data-collection period was 1,026 kilograms (1.13 tons). The BMPs removed a cumulative suspended-sediment load of 558 kilograms (0.62 ton). The BMPs tended to preferentially trap the sand-size fraction of the sediment, thereby releasing a greater percentage of fine-grained (silt and clay) sediment in the water discharging from the outlet. The preferential trapping of fine-grained sediment by the BMPs could explain, at least in part, why the BMPs were not successful at significantly reducing these constituents. In general, the four BMPs were not successful at significantly reducing fecal bacteria, nutrients, and total organic carbon (including associated properties of biochemical oxygen demand and chemical oxygen demand). Three of the four BMPs significantly lowered oil and grease concentrations before the stormwater discharged from the outlet. Additionally, only one BMP was effective at reducing all total and particulate trace-metal event-mean concentrations and particulate trace-metal event-mean loads in stormwater entering the inlet. With respect to trace-metal event-mean concentrations, however, minimal or no improvement in outlet water quality was observed for the four BMPs, and the majority of the outlet concentrations were above the established acute and chronic aquatic-life criteria by the South Carolina Department of Health and Environmental Control. No statistical differences among the removal-efficiency of the four BMPs were determined for suspended-sediment event-mean concentrations, total suspended solids event-mean concentrations, or oil and grease concentrations. These statistical findings indicated that differences among the mean efficiency ratios were not significant among the BMPs for these properties. Additionally, one BMP generally had statistically greater removal efficiency for total and particulate cadmium, copper, lead, and zinc than one or more of the other three BMPs. Statistical correlation tests were unable to identify a single major factor that would explain the high variability in inlet and outlet water concentrations and in removal efficiencies estimated by reduction percentage. Highly variable inlet and outlet concentrations for each BMP that produced highly variable reduction percentages were probably the result of multiple interacting factors, particularly rainfall intensity, the amount of rainfall between sampling events, traffic density, and the period of time since the last maintenance (clean out) of the BMP.

South Carolina↗
USGS2018

Simulation of groundwater flow and analysis of projected water use for the Rush Springs aquifer, western Oklahoma

The U.S. Geological Survey, in cooperation with the Bureau of Reclamation and the Oklahoma Water Resources Board, (1) quantified the groundwater resources of the Rush Springs aquifer in western Oklahoma by developing a numerical groundwater-flow model, (2) evaluated the effects of estimated equal-proportionate-share (EPS) pumping rates on aquifer storage and streamflow for time periods of 20, 40, and 50 years into the future, (3) assessed the uncertainty in the EPS scenario results, and (4) evaluated the effects of (a) projected groundwater-use rates extended 50 years into the future and (b) sustained hypothetical drought conditions over a 10-year period on stream base flow and groundwater in storage. The Rush Springs aquifer is an important source of water for municipal and irrigation use by many communities and agricultural users in the study area. The study area is composed of about 4,970 square miles (3,181,003 acres) of Rush Springs aquifer bedrock deposits located in 14 counties. The study area also includes the alluvium and terrace deposits of the Canadian and Washita Rivers, as well as alluvium along the Little Washita River, Deer Creek, and a number of smaller tributaries of the Washita River that overlie the bedrock. A numerical groundwater-flow model of the Rush Springs aquifer was constructed by using MODFLOW with the Newton solver. Groundwater flow was simulated for January 1979–December 2015 by using monthly stress periods, and an initial steady-state stress period was configured to represent mean annual inflows and outflows. The model was calibrated to groundwater-level observations at selected wells, monthly base flow at nine streamgages, stream seepage as estimated for the conceptual water budget, and Fort Cobb Reservoir stage. The EPS scenarios for the Rush Springs 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.82, 0.49, and 0.43 acre-foot per acre per year, respectively. Given the 2,954,545-acre aquifer area used for the EPS scenarios, the 20-year rate corresponds to an annual yield of about 2,422,727 acre-feet per year. Groundwater storage at the end of the 20-year EPS scenario was about 13,321,000 acre-feet, or about 31,516,437 acre-feet (70 percent) less than the starting EPS scenario storage. This decrease in storage was equivalent to a mean groundwater-level decline of about 152 feet. Water availability under the EPS pumping rate was primarily from the western area of the model. Saturation was sustained though the entire EPS scenario where the aquifer was sufficiently thick or a shallow hydraulic gradient was present. Fort Cobb Reservoir stage was below the dead-pool stage after about 5 years of 20-year EPS pumping. An uncertainty analysis was conducted to assess the uncertainty in the EPS scenario results. An ensemble of 400 random sets of possible parameter values was performed for the uncertainty analysis by using a multivariate normal distribution centered on the calibrated parameter values. The parameter bounds for the uncertainty analysis were determined by using the posterior covariance matrix, which allows for the incorporation of knowledge gained during the calibration process as well as observation uncertainty and the correlation between estimated parameters. The uncertainty results indicate a 95-percent confidence interval for the 20-year EPS pumping rate between 0.73 and 0.95 acre-foot per acre per year. Projected 50-year pumping scenarios were used to simulate the effects of selected well withdrawal rates on groundwater storage of the Rush Springs aquifer. The effects of well withdrawals were evaluated by comparing changes in groundwater storage between four 50-year scenarios using (1) no groundwater use, (2) mean groundwater use for the study period (1979–2015), (3) increasing groundwater use, and (4) groundwater use at the 2015 rate. The increasing-use scenario assumed a 38-percent increase in pumping over 50 years on the basis of 2010–60 demand projections for western Oklahoma. Simulated groundwater storage changes ranged between an increase of 6.3 percent for the scenario with no groundwater use, and 0.9 percent for the scenario with 2015 groundwater-use rates. For the Fort Cobb Reservoir surface watershed, simulated groundwater storage changes ranged between an increase of 23.6 percent for the scenario with no groundwater use and a decrease of 4.0 percent for the increasing groundwater-use scenario. Groundwater-level changes were generally greater in areas with a large concentration of groundwater wells and groundwater use such as the Fort Cobb Reservoir surface watershed. A hypothetical 10-year drought scenario was used to simulate the effects of a prolonged period of reduced recharge on the Rush Springs aquifer groundwater storage and Fort Cobb Reservoir stage and storage. Drought effects were quantified by comparing the results of the drought scenario to those of the calibrated numerical model. To simulate the hypothetical drought, recharge in the calibrated numerical model was reduced by 50 percent during the simulated drought period (1983–1992), and upstream inflows to the Canadian and Washita Rivers and associated tributaries were reduced by 37 percent. Groundwater storage at the end of the hypothetical drought period in December 1992 was about 42,983,000 acre-feet, or about 3,525,000 acre-feet (7.6 percent) less than the groundwater storage of the calibrated numerical model. This change in groundwater storage is equivalent to a mean groundwater-level decline of 15.8 feet. Simulated mean base-flow declines at the Canadian and Washita River streamgages were between 39 and 59 percent during the drought period. The minimum stage in Fort Cobb Reservoir at the end of the hypothetical drought period was 1,311 feet, indicating a storage capacity of only 10 percent of active conservation pool storage. The Fort Cobb Reservoir storage declines mostly resulted from reduced base flows in Cobb, Lake, and Willow Creeks upstream from the reservoir.

Oklahoma↗
USGS2016

Sediment loads and transport at constructed chutes along the Missouri River - Upper Hamburg Chute near Nebraska City, Nebraska, and Kansas Chute near Peru, Nebraska

The U.S. Geological Survey, in cooperation with the U.S. Army Corps of Engineers, monitored suspended sediment within constructed Missouri River chutes during March through October 2012. Chutes were constructed at selected river bends by the U.S. Army Corps of Engineers to help mitigate aquatic habitat lost through the creation and maintenance of the navigation channel on the Missouri River. The restoration and development of chutes is one method for creating shallow-water habitat within the Missouri River to meet requirements established by the amended 2000 Biological Opinion. Understanding geomorphic channel-evolution processes and sediment transport is important for the design of chutes, monitoring and maintenance of existing chutes, and characterizing the habitat that the chutes provide. This report describes the methods used to monitor suspended sediment at two Missouri River chutes and presents the results of the data analysis to help understand the suspended-sediment characteristics of each chute and the effect the chutes have on the Missouri River. Upper Hamburg chute, near Nebraska City, Nebraska, and Kansas chute, near Peru, Nebraska, were selected for monitoring. At each study site, monthly discrete samples were collected from April through October in the Missouri River main-channel transects upstream from the chute inlet, downstream from the chute outlet, at the outlet (downstream transect) of both chutes, and at the inlet (upstream transect) of Kansas chute. In addition, grab samples from all chute sampling locations were collected using autosamplers. Suspended-sediment concentration (SSC) and grain-size metrics were determined for all samples (discrete and grab). Continuous water-quality monitors recorded turbidity and water temperature at 15-minute intervals at the three chute sampling locations. Two acoustic Doppler velocimeters, one within each chute, measured water depth and current velocities continuously. The depth and velocity data were used to estimate streamflow within each chute. The sampling design was developed to understand the suspended-sediment differences within each chute and between the chute and the Missouri River main channel during discrete sampling. The sampling design also allowed for site-specific surrogate relations between SSC and turbidity to be developed, which could be used to compute real-time estimates of SSC and sediment loads within the chutes. Real-time estimates of SSC and sediment loads enable a better understanding of sediment transport within the chutes during times when physical samples are not collected, including periods of high flow. High flows during the summer of 2011 resulted in substantial alterations to both studied chutes; therefore, the U.S. Army Corps of Engineers repaired and modified both chutes during 2012. These unforeseen repairs and modifications within the chutes added uncertainty to the analysis because concentrations were altered by construction equipment and flow alteration. Daily suspended-sediment and suspended-silt loads were estimated based on surrogate relations with turbidity. A linear regression was used to estimate equal-width increment (EWI)-equivalent SSC from autosampler SSC before using the model-calibration dataset to determine the best-fit model for prediction of SSC from the turbidity and, in some cases, discharge. Correlation between suspended-sand concentration ( SSandC ) in EWI samples and concurrent samples collected by an autosampler was low; therefore, SSandC was excluded from development of surrogate relations because a large part of the calibration dataset was from autosamples. Instead, SSandC was estimated as SSC minus suspended-silt-clay concentration ( SSiltC ). At all sites, the best-fit models included the base-10 logarithm of concentration and turbidity, and at Kansas chute upstream, the base-10 logarithm of streamflow was also included in the best-fit models. These surrogate models were used to estimate continuous time series of SSC and SSiltC . Estimated concentrations of suspended sediment were used to estimate instantaneous and daily loads for total suspended sediment, suspended silt-clay, and suspended sand. Estimated daily suspended-sediment loads were not significantly different between upstream and downstream transects within the Kansas chute, and most individual daily loads within the chute were not significantly different between upstream and downstream transects when evaluated using overlap in daily 95-percent confidence intervals. The comparison of daily load values for upstream and downstream chute transects, as estimated from turbidity-based surrogate models for Kansas chute, documents the daily dynamic nature of sediment transport within the chute with a temporal resolution that is not practical with discrete suspended-sediment sampling alone. Comparisons of concentrations and loads from EWI samples collected from different transects within a study site resulted in few significant differences, but comparisons are limited by small sample sizes and large within-transect variability. When comparing the Missouri River upstream transect to the chute inlet transect, similar results were determined in 2012 as were determined in 2008—the chute inlet affected the amount of sediment entering the chute from the main channel. In addition, the Kansas chute is potentially affecting the sediment concentration within the Missouri River main channel, but small sample size and construction activities within the chute limit the ability to fully understand either the effect of the chute in 2012 or the effect of the chute on the main channel during a year without construction. Finally, some differences in SSC were detected between the Missouri River upstream transects and the chute downstream transects; however, the effect of the chutes on the Missouri River main-channel sediment transport was difficult to isolate because of construction activities and sampling variability.

Nebraska↗
USGS2013

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

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

Alabama, Florida, Georgia↗
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↗