Search USGS⌕ Search

SEARCH · Search USGS

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 721 records · Page 40Linked to original sources

Trends in streamflow and concentrations and flux of nutrients and total suspended solids in the Upper White River at Muncie, near Nora, and near Centerton, Indiana

The U.S. Geological Survey (USGS), in cooperation with The Nature Conservancy, completed a study to estimate and assess trends in streamflow and annual mean concentrations and flux of nutrients (nitrate plus nitrite, total Kjeldahl nitrogen, and total phosphorus) and total suspended solids at three USGS streamgages (hereafter referred to as “study gages”) on the Upper White River at Muncie (USGS station 03347000), near Nora (USGS station 03351000), and near Centerton (USGS station 03354000), Indiana. Water-quality data used in the analyses were collected by several agencies between calendar years 1991 and 2017, and streamflow (discharge) data were collected by the USGS. For most of the water-quality constituents, there were suitable data to facilitate an analysis of the 26-year period extending from calendar years 1991 to 2017 (water years 1992 to 2017); however, shorter analytical periods were necessary for total Kjeldahl nitrogen for the study gages at Muncie and near Centerton and for total suspended solids for the study gage near Centerton. Temporal trends in streamflows at the study gages for the period extending from water years 1978 to 2017 were assessed using Exploration and Graphics for RivEr Trends (EGRET) and Mann-Kendall and Pettitt tests. With just one exception, the annual maximum and mean daily streamflows and the annual minimum 7-day mean streamflows at the study gages demonstrated upward trends (increasing streamflows) in the EGRET analyses. The exception was the annual 7-day minimum streamflow at the study gage near Nora, which indicated no trend. Mann-Kendall tests also indicated that the average trend for the annual maximum daily, annual mean daily, and annual 7-day minimum streamflow statistics between water years 1978 and 2017 was upward at each of the study gages; however, only the trends in the annual mean daily streamflows at the study gage at Muncie and the annual maximum daily streamflows at the study gages near Nora and near Centerton were statistically significant at a 0.05 probability level. The Pettitt tests indicated that a statistically significant step trend (abrupt change) in annual mean daily streamflows occurred at each of the study gages around water year 2001. The seasonal distributions of total suspended solids, total phosphorus, nitrate plus nitrite, and total Kjeldahl nitrogen concentrations at the study gages were evaluated to identify patterns and other distinguishing characteristics by examining boxplots of concentrations as a function of month of the year. Seasonal distributions of nitrate plus nitrite concentrations and total suspended solids concentrations differed from each other but were generally similar among the three study gages for a given constituent. Median concentrations of nitrate plus nitrite were highest during the January–June months, whereas median concentrations of total suspended solids were highest during June and July. Seasonal distributions of total phosphorus concentrations were similar at the study gages near Nora and near Centerton, but the seasonal distribution was noticeably different at the study gage at Muncie, which had monthly median concentrations that were substantially lower than at the two downstream study gages (near Nora and near Centerton). The seasonal distribution of total Kjeldahl nitrogen concentrations differed in pattern among the three study gages; however, in general, some of the higher monthly median total Kjeldahl nitrogen concentrations at each study gage were associated with the late spring and summer periods. The Weighted Regressions on Time, Discharge, and Season (WRTDS) method implemented in EGRET was used to estimate water-year annual mean daily concentrations and flux of nutrients and total suspended solids, as well as estimates of concentrations and flux that were “normalized” to remove the effect of year-to-year variation in streamflow. The approximate coefficients of determination for the WRTDS regression models ranged from a high of 0.82 for total phosphorus for the study gage near Centerton to a low of 0.19 for nitrate plus nitrite for the study gage near Nora. Loads and yields of total suspended solids, total phosphorus, nitrate plus nitrite, and total Kjeldahl nitrogen were estimated for analytical periods consisting of the longest periods of concurrent record at the three study gages. Loads of each of the constituents increased sequentially from the most upstream study gage to the most downstream study gage; however, the same was not true for yields. The highest yields of total suspended solids, total phosphorus, and total Kjeldahl nitrogen occurred at the most upstream study gage (at Muncie); however, the highest yield of nitrate plus nitrite occurred at the most downstream study gage (near Centerton). WRTDS bootstrap tests were used to assess the magnitude, direction, and likelihood of changes in annual flow-normalized mean daily concentrations and flux of total suspended solids, total phosphorus, nitrate plus nitrite, and total Kjeldahl nitrogen at the study gages between water years 1997 and 2017. Changes in flow-normalized concentrations and flux of the constituents between water years 1997 and 2017 were mostly downward (decreasing). The exceptions were likely to highly likely upward (increasing) changes in (1) flow-normalized annual mean daily concentration and annual flux for total suspended solids and total phosphorus at the study gage at Muncie, (2) flow-normalized annual mean daily total phosphorus concentration at the study gage near Centerton, (3) flow-normalized annual flux of total phosphorus at the study gage near Centerton, and (4) flow-normalized annual mean daily nitrate plus nitrite concentration at the study gage near Centerton. Although an upward change in flow-normalized nitrate plus nitrite concentrations was likely at the study gage near Centerton, flow-normalized annual flux of nitrate plus nitrite at that study gage was determined to have a highly likely downward change. EGRET and Exploration and Graphics for RivEr Trends Confidence Intervals (EGRETci) analyses can be used to improve our understanding of how concentrations and flux change as functions of time and streamflow, as well as provide information on how the relations between streamflow and constituent concentrations have changed within the calendar year between any 2 years included in the analyses. Examples of those uses, illustrating changes between calendar years 1992 and 2017, were given for total suspended solids concentrations at the study gage near Nora and for nitrate plus nitrite concentrations at the study gage near Centerton.

Indiana↗

Hydrogeology and Aquifer Storage and Recovery Performance in the Upper Floridan Aquifer, Southern Florida

Well construction, hydraulic well test, ambient water-quality, and cycle test data were inventoried and compiled for 30 aquifer storage and recovery facilities constructed in the Floridan aquifer system in southern Florida. Most of the facilities are operated by local municipalities or counties in coastal areas, but five sites are currently being evaluated as part of the Comprehensive Everglades Restoration Plan. The relative performance of all sites with adequate cycle test data was determined, and compared with four hydrogeologic and design factors that may affect recovery efficiency. Testing or operational cycles include recharge, storage, and recovery periods that each last days or months. Cycle test data calculations were made including the potable water (chloride concentration of less than 250 milligrams per liter) recovery efficiency per cycle, total recovery efficiency per cycle, and cumulative potable water recovery efficiencies for all of the cycles at each site. The potable water recovery efficiency is the percentage of the total amount of potable water recharged for each cycle that is recovered; potable water recovery efficiency calculations (per cycle and cumulative) were the primary measures used to evaluate site performance in this study. Total recovery efficiency, which is the percent recovery at the end of each cycle, however, can be substantially higher and is the performance measure normally used in the operation of water-treatment plants. The Upper Floridan aquifer of the Floridan aquifer system currently is being used, or planned for use, at 29 of the aquifer storage and recovery sites. The Upper Floridan aquifer is continuous throughout southern Florida, and its overlying confinement is generally good; however, the aquifer contains brackish to saline ground water that can greatly affect freshwater storage and recovery due to dispersive mixing within the aquifer. The hydrogeology of the Upper Floridan varies in southern Florida; confinement between flow zones is better in southwestern Florida than in southeastern Florida. Vertical hydraulic conductivity in the upper part of the aquifer also may be higher in southeastern Florida because of unconformities present at formation contacts within the aquifer that may be better developed in this area. Recovery efficiencies per cycle varied widely. Eight sites had recovery efficiencies of less than about 10 percent for the first cycle, and three of these sites had not yet achieved recoveries exceeding 10 percent, even after three to five cycles. The highest recovery efficiency achieved per cycle was 94 percent. Three southeastern coastal sites and two southwestern coastal sites have achieved potable water recoveries per cycle exceeding 60 percent. One of the southeastern coastal sites and both of the southwestern coastal sites achieved good recoveries, even with long storage periods (from 174 to 191 days). The high recovery efficiencies for some cycles apparently resulted from water banking?an operational approach whereby an initial cycle with a large recharge volume of water is followed by cycles with much smaller recharge volume. This practice flushes out the aquifer around the well and builds up a buffer zone that can maintain high recovery efficiency in the subsequent cycles. The relative performance of all sites with adequate cycle test data was determined. Performance was arbitrarily grouped into ?high? (greater than 40 percent), ?medium? (between 20 and 40 percent), and ?low? (less than 20 percent) categories based primarily on their cumulative recovery efficiency for the first seven cycles, or projected to seven cycles if fewer cycles were conducted. The ratings of three sites, considered to be borderline, were modified using the overall recharge rate derived from the cumulative recharge volumes. A higher overall recharge rate (greater than 300 million gallons per year) can improve recovery efficiency because of the water-bankin

Scientific Investigations Report↗

Estimating spatially and temporally varying recharge and runoff from precipitation and urban irrigation in the Los Angeles Basin, California

A daily precipitation-runoff model, referred to as the Los Angeles Basin watershed model (LABWM), was used to estimate recharge and runoff for a 5,047 square kilometer study area that included the greater Los Angeles area and all surface-water drainages potentially contributing recharge to a 1,450 square kilometer groundwater-study area underlying the greater Los Angeles area, referred to as the Los Angeles groundwater-study area. The recharge estimates for the Los Angeles groundwater-study area included spatially distributed recharge in response to the infiltration of precipitation, runoff, and urban irrigation, as well as mountain-front recharge from surface-water drainages bordering the groundwater-study area. The recharge and runoff estimates incorporated a new method for estimating urban irrigation, consisting of residential and commercial landscape watering, based on land use and the percentage of pervious land area. The LABWM used a 201.17-meter gridded discretization of the study area to represent spatially distributed climate and watershed characteristics affecting the surface and shallow sub-surface hydrology for the Los Angeles groundwater study area. Climate data from a local network of 201 monitoring sites and published maps of 30-year-average monthly precipitation and maximum and minimum air temperature were used to develop the climate inputs for the LABWM. Published maps of land use, land cover, soils, vegetation, and surficial geology were used to represent the physical characteristics of the LABWM area. The LABWM was calibrated to available streamflow records at six streamflow-gaging stations. Model results for a 100-year target-simulation period, from water years 1915 through 2014, were used to quantify and evaluate the spatial and temporal variability of water-budget components, including evapotranspiration (ET), recharge, and runoff. The largest outflow of water from the LABWM was ET; the 100-year average ET rate of 362 millimeters per year (mm/yr) accounted for 66 percent of the combined water inflow of 551 mm/yr, including 488 mm/yr from precipitation and 63 mm/yr from urban irrigation. The simulated ET rate varied from a minimum of 0 mm/yr for impervious areas to high values of more than 1,000 mm/yr for many areas, including the south-facing slopes of the San Gabriel Mountains, stream channels underlain by permeable soils and thick root zones, and pervious locations receiving inflows both from urban irrigation and surface water. Runoff was the next largest outflow, averaging 145 mm/yr for the 100-year period, or 26 percent of the combined precipitation and urban-irrigation inflow. Recharge averaged 45 mm/yr, or about 8 percent of the combined inflow from precipitation and urban irrigation. Simulation results indicated that recharge in response to urban irrigation was an important component of spatially distributed recharge, contributing an average of 56 percent of the total recharge to the eight LABWM subdomains containing the Los Angeles groundwater study area. The 100‑year average recharge rate for the eight subdomains was 41 mm/yr, or 8,473 hectare-meters per year (ha-m/yr), with urban irrigation included in the simulation compared to a recharge rate of 18 mm/yr, or 3,741 ha-m/yr, with urban irrigation excluded. In contrast to recharge, the effect of urban irrigation on runoff was slight; runoff was 72,667 ha-m/yr with urban irrigation included compared to 72,618 ha-m/yr with urban irrigation excluded, an increase of only 48 ha-m/yr (about 0.1 percent). Simulation results also indicated that potential recharge from hilly drainages outside of, but bordering and tributary to, the lower-lying area of the Los Angeles groundwater study area, in this study referred to as mountain-front recharge, could provide an important contribution to the total recharge for the groundwater basins. The time-averaged recharge rate was similar to the combined direct and mountain-front recharge components estimated in a previous study and used as input for a calibrated groundwater model. The annual (water year) recharge estimates simulated in this study, however, indicated much greater year-to-year variability, which was dependent on year-to-year variability in the magnitude and distribution of daily precipitation, compared to the previous estimates.

California↗

Modeling Escherichia coli in the Missouri River near Omaha, Nebraska, 2012–16

The city of Omaha, Nebraska, has a combined sewer system in some areas of the city. In Omaha, Nebr., a moderate amount of rainfall will lead to the combination of stormwater and untreated sewage or wastewater being discharged directly into the Missouri River and Papillion Creek and is called a combined sewer overflow (CSO) event. In 2009, the city of Omaha began the implementation of their Long Term Control Plan (LTCP) to mitigate the effects of CSOs on the Missouri River and Papillion Creek. As part of the LTCP, the city partnered with the U.S. Geological Survey (USGS) in 2012 to begin monitoring in the Missouri River. Since 2012, monthly discrete water-quality samples for many constituents have been collected from the Missouri River at four sites. At 3 of the 4 sites, water quality has been monitored continuously for selected constituents and physical properties. These discrete water-quality samples and continuous water-quality monitoring data (from July 2012 to 2020) have been collected to better understand the water quality of the Missouri River, how it is changing with time, how it changes upstream from the city of Omaha to downstream, and how it varies during base-flow conditions and during periods of runoff. The purpose of this report is to document the development of Escherichia coli ( E. coli ) concentration models for these four Missouri River sites. Analysis was completed using the first 5 years of data (through 2016) to determine if the current approach is sufficient to meet future analysis goals and to understand if proposed models such as Load Estimator (LOADEST) models will be able to represent water-quality changes in the Missouri River. Multiple linear regression models were developed to estimate E. coli concentration using LOADEST as implemented in the rloadest package in the R statistical software program. A set of explanatory variables, including streamflow and streamflow anomalies, precipitation, information about CSOs, and continuous water quality, were evaluated for potential inclusion in regression models. The best model at Missouri River at NP Dodge Park at Omaha, Nebr. (USGS station 412126095565201; hereafter “NP Dodge”) included basin explanatory variables of upstream antecedent precipitation index measured at Tekamah, Nebr.; decimal time; season; and turbidity. The best model at Missouri River at Freedom Park Omaha, Nebr. (USGS station 411636095535401; hereafter “Freedom Park”) included the same explanatory variables as the NP Dodge model with the addition of turbidity anomalies and flow anomalies. The best models at the two downstream sites (Missouri River near Council Bluffs, Iowa, USGS station 06610505 and Missouri River near La Platte, Nebr., USGS station 410333095530101) included the same explanatory variables as the Freedom Park model with the addition of local antecedent precipitation index as measured at Eppley Airport in Omaha, Nebr., and additional turbidity and flow anomalies. The final selected models were the best models given our modeling design constraint in which explanatory variables included in the model for the upstream site were included in the downstream models. Explanatory variables currently (2020) being collected and included in the selected models through 2016 explained 64–75 percent of the variability of E. coli concentration in the Missouri River. Explaining 64–75 percent of the variability might be considered low when working with physical constituents (total nitrogen or sediment), but with the natural variability of biological constituents such as E. coli , the uncertainty of E. coli laboratory measurements, and the added complexity of modeling in a large drainage basin with multiple sources, these results are adequate and indicate that the explanatory variables being collected and models such as LOADEST can represent water-quality changes in the Missouri River for E. coli concentration from 2012 to 2016.

Nebraska↗

Variability in stream chemistry in relation to urban development and biological condition in seven metropolitan areas of the United States, 1999-2004

Beginning in 1999, the U.S. Geological Survey National Water Quality Assessment Program investigated the effects of urban development on stream ecosystems in nine metropolitan study areas across the United States. In seven of these study areas, stream-chemistry samples were collected every other month for 1 year at 6 to 10 sites. Within a study area, the sites collectively represented a gradient of urban development from minimally to highly developed watersheds, based on the percentage of urban land cover; depending on study area, the land cover before urban development was either forested or agricultural. The stream-chemistry factors measured in the samples were total nitrogen, total phosphorus, chloride, and pesticide toxicity. These data were used to characterize the stream-chemistry factors in four ways (hereafter referred to as characterizations)—seasonal high-flow value, seasonal low-flow value, the median value (representing a single integrated value of the factor over the year), and the standard deviation of values (representing the variation of the factor over the year). Aquatic macroinvertebrate communities were sampled at each site to infer the biological condition of the stream based on the relative sensitivity of the community to environmental stressors. A Spearman correlation analysis was used to evaluate relations between (1) urban development and each characterization of the stream-chemistry factors and (2) the biological condition of a stream and the different characterizations of chloride and pesticide toxicity. Overall, the study areas where the land cover before urban development was primarily forested had a greater number of moderate and strong relations compared with the study areas where the land cover before urban development was primarily agriculture; this was true when urban development was correlated with the stream-chemistry factors (except chloride) and when chloride and pesticide toxicity was correlated with the biological condition. Except for primarily phosphorus in two study areas, stream-chemistry factors generally increased with urban development, and among the different characterizations, the median value typically indicated the strongest relations. The variation in stream-chemistry factors throughout the year generally increased with urban development, indicating that water quality became less consistent as watersheds were developed. In study areas with high annual snow fall, the variation in chloride concentrations throughout the year was particularly strongly related to urban development, likely a result of road salt applications during the winter. The relations of the biological condition to chloride and pesticide toxicity were calculated irrespective of urban development, but the overall results indicated that the relations were still stronger in the study areas that had been forested before urban development. The weaker relations in the study areas that had been agricultural before urban development were likely the results of biological communities having been degraded from agricultural practices in the watersheds. Collectively, these results indicated that, compared with sampling a stream at a single point in time, sampling at regular intervals during a year may provide a more representative measure of water quality, especially in the areas of high urban development where water quality fluctuated more widely between samples. Furthermore, the use of "integrated" values of stream chemistry factors may be more appropriate when assessing relations to the biological condition of a stream because the taxa composition of a biological community typically reflects the water-quality conditions over time.

Scientific Investigations Report↗

Hydrology and simulation of ground-water flow in Cedar Valley, Iron County, Utah

Cedar Valley, located in the eastern part of Iron County in southwestern Utah, is experiencing rapid population growth. Cedar Valley traditionally has supported agriculture, but the growing population needs a larger share of the available water resources. Water withdrawn from the unconsolidated basin fill is the primary source for public supply and is a major source of water for irrigation. Water managers are concerned about increasing demands on the water supply and need hydrologic information to manage this limited water resource and minimize flow of water unsuitable for domestic use toward present and future public-supply sources. Surface water in the study area is derived primarily from snowmelt at higher altitudes east of the study area or from occasional large thunderstorms during the summer. Coal Creek, a perennial stream with an average annual discharge of 24,200 acre-feet per year, is the largest stream in Cedar Valley. Typically, all of the water in Coal Creek is diverted for irrigation during the summer months. All surface water is consumed within the basin by irrigated crops, evapotranspiration, or recharge to the ground-water system. Ground water in Cedar Valley generally moves from primary recharge areas along the eastern margin of the basin where Coal Creek enters, to areas of discharge or subsurface outflow. Recharge to the unconsolidated basin-fill aquifer is by seepage of unconsumed irrigation water, streams, direct precipitation on the unconsolidated basin fill, and subsurface inflow from consolidated rock and Parowan Valley and is estimated to be about 42,000 acre-feet per year. Stable-isotope data indicate that recharge is primarily from winter precipitation. The chloride mass-balance method indicates that recharge may be less than 42,000 acre-feet per year, but is considered a rough approximation because of limited chloride concentration data for precipitation and Coal Creek. Continued declining water levels indicate that recharge is not sufficient to meet demand. Water levels in many areas are at or close to historic lows. In 2000, withdrawal from wells was estimated to be 36,000 acre-feet per year. About 4,000 acre-feet per year are estimated to discharge to evapotranspiration or as subsurface outflow. Prior to large-scale ground-water development, ground-water discharge by evapotranspiration and discharge to springs was much larger. Ground water along the eastern margin of the valley between Cedar City and Enoch is unsuitable for domestic use because of high dissolved-solids and nitrate concentrations. The predominant ions of Ca and SO4 in this area indicate dissolution of gypsum in the Markagunt Plateau to the east. Data collected during this study were compared to historic data; there is no evidence to indicate deterioration in ground-water quality. The spatial distribution of ground water with high nitrate concentration does not appear to be migrating beyond its previously known extent. No single source can be identified as the cause for elevated nitrate concentrations in ground water. Low nitrogen-15 values north of Cedar City indicate a natural geologic source. Higher nitrogen-15 values toward the center of the basin and associated hydrologic data indicate probable recharge from waste-water effluent. Excess dissolved nitrogen gas and low nitrate concentrations in shallow ground water indicate that denitrification is occurring in some areas. A computer ground-water flow model was developed to simulate flow in the unconsolidated basin fill. The method of determining recharge from irrigation was changed during the calibration process to incorporate more areal and temporal variability. In general, the model accurately simulates water levels and water-level fluctuations and can be considered an adequate tool to help determine the valley-wide effects on water levels of additional ground-water withdrawals and changes in water use. The model was used to simulated water-level changes caused by projecting current withdrawal rates, increased withdrawal rates, and a 10-year drought. Water levels declined 20 to 275 feet in the southern and central parts of the valley and less than 20 feet north of Enoch

Utah↗

Organic wastewater compounds in water and sediment in and near restored wetlands, Great Marsh, Indiana Dunes National Lakeshore, 2009–11

A cooperative investigation between the U.S. Geological Survey and the National Park Service was completed from 2009 through 2011 to understand the occurrence, distribution, and environmental processes affecting concentrations of organic wastewater compounds in water and sediment in and near Great Marsh at the Indiana Dunes National Lakeshore in Beverly Shores, Indiana. Sampling sites were selected to represent hydrologic inputs to the restored wetlands from adjacent upstream residential and less developed areas and to represent discharge points of cascading cells within the restored wetland. A multiphase approach was used for the investigation. Discrete water samples and time-integrated passive samples were analyzed for 69 organic wastewater compounds. Continuous water-level information and periodic streamflow measurements characterized flow conditions at discharge points from restored wetland cells. Wetland sediments were collected and analyzed for sorptive losses of organic wastewater compounds and to evaluate of the potential for wetland sediments to biotransform organic wastewater compounds. A total of 52 organic wastewater compounds were detected in discrete water samples at 1 or more sites. Detections of organic wastewater compounds were widespread, but concentrations were generally low and 95 percent were less than 2.1 micrograms per liter. Six compounds were detected at concentrations greater than 2.1 micrograms per liter—four fecal sterols (beta-sitosterol, cholesterol, beta-stigmastanol, and 2-beta coprostanol), one plasticizer (bis-2-ethylhex ylphthalate), and a non-ionic detergent (4-nonylphenol diethoxylate). Two 1-month deployments of time-integrative passive samplers, called polar organic chemical integrative samplers, detected organic wastewater compounds at lower concentrations than were possible with discrete water samples. Isopropyl benzene (solvent), caffeine (plant alkaloid, stimulant), and hexahydrohexamethyl cyclopentabenzopyran (fragrance) were detected in more than half of the extracts from passive samplers, but they were not detected in any discrete water sample. The Yeast Estrogen Screen assay identified measurable estrogenicity in one passive sampler extract from the most downstream wetland site in both the April and November–December 2011 deployments and in passive sampler extracts from one residential and one upstream site in the November–December 2011 deployment only. Surface-water levels in the restored wetland cells were monitored continuously using submersible pressure transducers in hand-driven well points screened in the surface water. Surface-water levels in the wetland cells responded quickly to precipitation and substantially receded within 2 days following the largest rainfall events. Seasonal patterns in water levels generally showed higher and more variable surface-water levels in the wetland cells during spring and early summer. Water levels in the wetland cells fell below the elevation of the control structures and ceased to flow over the spillways during extended dry periods (primarily late summer and early fall). Daily loads of seven organic wastewater compounds, as indicators of septic system effluent, were estimated for samples collected at wetland outlet spillways when flow measurements could be made. Median daily loads of the indicator organic wastewater compounds increased in downstream order, and the largest median loads were measured at the most downstream site. Median daily loads were higher for samples collected in spring and summer than those collected in fall, as the higher seasonal water levels increased streamflow at the wetland outlet spillways. Wetland sediment samples were analyzed for 84 organic wastewater compounds, polycyclic aromatic hydrocarbons, and semivolatile organic compounds to investigate the fate of contaminants in Great Marsh. The top five detected compounds by total mass in wetland sediment samples were beta-sitosterol, beta-stigmastanol, cholesterol, bis(2-ethylhexyl) phthalate, and phenol. Polycyclic aromatic hydrocarbons also were frequently detected in wetland sediment samples. Source apportionment of polycyclic aromatic hydrocarbon detections indicated atmospheric sources of pyrogenic compounds, rather than residential sources. Comparisons of polycyclic aromatic hydrocarbon concentrations in wetland sediment samples to sediment quality target guidelines indicated the potential for harmful effects on sediment-dwelling organisms at several sites. Biodegradation of select endocrine-disrupting compounds (17α-ethinylestradiol, 4-nonylphenol, triclocarban, and bisphenol A) in shallow wetland sediments was evaluated in laboratory experiments by using carbon-14 radiolabeled model contaminants. Substantial biodegradation of certain organic wastewater compounds were demonstrated, primarily in oxic (oxygen containing) environments. One of four modeled compounds, bisphenol A, was biodegraded in anoxic (oxygen free) environments. Only sediments collected nearest residential areas exhibited degradation of the synthetic birth control pharmaceutical, 17α-ethinylestradiol, possibly owing to adaptation and acclimation of the indigenous microbial community to septic discharge and the resultant selection of a microbial capability for biodegradation of 17α-ethinylestradiol.

Indiana↗

Movement of the saltwater interface in the surficial aquifer system in response to hydrologic stresses and water-management practices, Broward County, Florida

A study was conducted to evaluate the relation between water-level fluctuations and saltwater intrusion in Broward County, Florida. The objective was achieved through data collection at selected wells in Broward County and through the development of a variable-density ground-water flow model. The numerical model is representative of many locations in Broward County that contain a well field, control structure, canal, the Intracoastal Waterway, and the Atlantic Ocean. The model was used to simulate short-term movement (from tidal fluctuations to monthly changes) and long-term movement (greater than 10 years) of the saltwater interface resulting from changes in rainfall, well-field withdrawals, sea-level rise, and upstream canal stage. The SEAWAT code, which is a combined version of the computer codes, MODFLOW and MT3D, was used to simulate the complex variable-density flow patterns. Model results indicated that the canal, control structure, and sea level have major effects on ground-water flow. For periods greater than 10 years, the upstream canal stage controls the movement and location of the saltwater interface. If upstream canal stage is decreased by 1 foot (0.3048 meter), the saltwater interface takes 50 years to move inland and stabilize. If the upstream canal stage is then increased by 1 foot (0.3048 meter), the saltwater interface takes 90 years to move seaward and stabilize. If sea level rises about 48 centimeters over the next 100 year as predicted, then inland movement of the saltwater interface may cause well-field contamination. For periods less than 10 years, simulation results indicated that a 3-year drought with increased well-field withdrawals probably will not have long-term effects on the position of the saltwater interface in the Biscayne aquifer. The saltwater interface returns to its original position in less than 10 years. Model results, however, indicated that the interface location in the lower part of the surficial aquifer system takes longer than 10 years to recover from a drought. Additionally, rainfall seems to have the greatest effect on saltwater interface movement in areas some distance from canals, but the upstream canal stage has the greatest effect on the movement of the saltwater interface near canals. Field data indicated that saltwater interface movement includes short-term fluctuations caused by tidal fluctuations and long-term seasonal fluctuations. Statistical analyses of daily-averaged data indicated that the saltwater interface moves in response to pumpage, rainfall, and upstream canal stage. In areas near the canal, the saltwater interface is most affected by canal stage because water-management structures control the stage in the upstream part of the canal and allow movement of the saltwater interface. In areas away from the canal, the saltwater interface is most affected by pumpage and rainfall, depending on the location of well fields. Data analyses also revealed that rainfall changes the vertical flow direction in the Biscayne aquifer. Results from the study indicated that upstream canal stage substantially affects the long-term position of the saltwater interface in the surficial aquifer system. The saltwater interface moves faster inland than seaward because of changes in upstream canal stage. For short-term problems, such as drought, the threat of saltwater intrusion in the Biscayne aquifer does not appear to be severe if the well-field withdrawal is increased; however, this conclusion is based on the assumption that well-field withdrawals will decrease once the drought is over. Sea-level rise may be a potential threat to the water supply in Broward County as the saltwater interface moves inland toward well fields.

Florida↗

Occurrence of volatile organic compounds in selected urban streams in the United States, 1995-2003

As part of the U.S. Geological Survey's (USGS) National Water-Quality Assessment (NAWQA) Program, urban indicator sites were monitored to (1) characterize the stream quality from drainage basins with predominantly residential and commercial land use, and (2) determine which selected natural and anthropogenic factors affect stream quality. A total of 869 water samples were collected from 37 urban streams during 1995-2003 and were analyzed for 87 volatile organic compounds (VOCs). The occurrence of VOCs in urban streams is described in this report for (1) all samples as a single dataset, (2) all samples grouped by streamflow pentiles, and (3) all samples grouped by warmer (April through September) and cooler (October through March) months by the detection frequency and (or) concentration of (a) any VOC, (b) VOC groups, and (c) individual compounds. An assessment level of 0.02 microgram per liter (ug/L) was used to compute the detection frequencies and concentrations of VOCs. Concentrations of VOCs were compared to (1) U.S. Environmental Protection Agency's (USEPA) drinking-water Maximum Contaminant Levels (MCLs) or Drinking Water Advisories, (2) Health-Based Screening Levels (HBSLs) developed by the USGS in collaboration with the USEPA and other agencies, and (3) USEPA and Canadian aquatic-life criteria. One or more VOCs were detected in 97.1 percent of 869 samples, and one or more VOCs were detected frequently (greater than 80 percent) at all sites. The median total VOC concentration for all samples was 0.57 ug/L, and total VOC concentrations in a single sample ranged from not detected to 698 ug/L. About 85 percent of the samples contained two or more VOCs, and about one-half contained five or more VOCs. The gasoline hydrocarbons were the most frequently occurring VOC group followed by solvents, trihalomethanes (THMs), gasoline oxygenates, organic synthesis compounds, fumigants, and refrigerants. Concentration ranges for most VOC groups were distributed over at least two orders of magnitude. Fifty-seven of the 87 VOCs analyzed were detected in at least one sample at an assessment level of 0.02 ug/L. More than one-half of the 30 VOCs not detected in samples were organic synthesis compounds. Fifteen compounds had detection frequencies greater than or equal to 10 percent. With the exception of toluene and chloroform, the median concentration of each VOC for all samples was less than the assessment level. Furthermore, the median concentrations of detections for the 15 most frequently occurring VOCs ranged from 0.03 to 3.9 ug/L, and typically were less than or equal to 0.10 ug/L. The 869 samples from the 37 sites were stratified into five streamflow pentiles (less than 20, 20-less than 40, 40-less than 60, 60-less than 80, and greater than or equal to 80 percent of estimated long-term streamflow statistics) for comparison of the occurrence of VOCs. The detection frequency of one or more VOCs by streamflow pentile varied only slightly from 96.7 to 97.7 percent. The median total VOC concentrations in samples for the five streamflow pentiles ranged from 0.39 to 1.0 ug/L. Two or more VOCs were present in more than 80 percent of samples in each of the five pentiles. The gasoline hydrocarbons, solvents, THMs, and gasoline oxygenates occurred frequently (greater than 30 percent) in all streamflow pentiles, in contrast to the organic synthesis compounds, fumigants, and refrigerants that occurred less frequently in urban streams under all streamflow conditions. The median total VOC concentrations for gasoline hydrocarbons, solvents, gasoline oxygenates, and organic synthesis compounds generally increased as streamflow increased. In contrast, the median total VOC concentrations for THMs and fumigants generally decreased as streamflow increased. The median total VOC concentrations for refrigerants showed no pattern as streamflow increased. Because differences between VOC occurrence and streamflow pentiles were small for most compariso

Scientific Investigations Report↗

Trends in groundwater levels in and near the Rosebud Indian Reservation, South Dakota, water years 1956–2017

The U.S. Geological Survey (USGS), in cooperation with the Rosebud Sioux Tribe, completed a study to characterize water-level fluctuations in observation wells to examine driving factors that affect water levels in and near the Rosebud Indian Reservation, which comprises all of Todd County. The study investigates concerns regarding potential effects of groundwater withdrawals and climate conditions on groundwater levels within an area that includes Todd County and a surrounding area that extends 10 miles north, east, and west of the county border. Characterization of water-level fluctuations in observation wells and relative driving factors was accomplished by statistical trend analysis. Two statistical methods were used for analysis of temporal trends for climatic and hydrologic data. To determine which trend analysis to use, applicable datasets were tested for statistically significant short-term persistence (STP). In the absence of significant STP, existence of statistical trends was determined using the standard Mann-Kendall test for probability values less than or equal to 0.10 (90-percent confidence level); however, a modified Mann-Kendall test was used for datasets where statistically significant STP was detected. Trend magnitudes were computed using the Sen’s slope estimator. Monthly data from the Parameter-elevation Regressions on Independent Slopes Model (PRISM) were aggregated to obtain annual and seasonal datasets for total precipitation, minimum air temperature ( T min ), and maximum air temperature ( T max ) for the study area and a surrounding buffer area. Trend tests for total precipitation, T min , and T max were completed for annual and seasonal time series for water years 1956–2017, which is about 2 years before the earliest available water-level measurements. A 2-year offset was arbitrarily selected because scrutiny of water-level and precipitation data indicated that responses of groundwater levels for many of the observation wells lagged major changes in precipitation patterns by about 2 years. Statistically significant upward trends were detected for annual precipitation and annual T min for almost all of the study area and the surrounding buffer area. Statistically significant downward trends in T max were detected for a very small part of the study area; however, the sparse spatial coverage reduces confidence that these are true trends. Spatial distributions of statistically significant trends in seasonal climate data were generally similar to the annual trends, but with substantial differences in the spatial density of the trends. Groundwater trends for 58 observation wells were analyzed for three separate water-level parameters (minimum, median, and maximum) because wells are measured sporadically and data are biased towards more frequent measurements during periods of heaviest irrigation demand. Trends in the time series of annual precipitation (from PRISM) starting 2 years earlier than for the associated water-level trend also were analyzed for the location of each individual observation well. Sen’s slope and Mann-Kendall probability values (p-values) were computed for the three water-level parameters and for the annual precipitation time series. Graphs showing results of trend analyses for each observation well also showed changes over time in the sum of licensed groundwater withdrawals within six specified radii (0.5, 1, 2, 3, 4, and 5 miles) of each well as a qualitative indicator of proximal groundwater demand. Of all 58 observation wells considered, 28 wells had significant upward trends for at least one of the three water-level parameters, 11 wells had significant downward trends for at least one water-level parameter, and 19 wells did not have any significant trends. Significant upward trends in annual precipitation were detected for 48 of the 58 wells. Results of trend analyses likely show the effects of groundwater withdrawals on water levels in the Ogallala aquifer in areas of substantial demand. Precipitation trends are significantly upward for 43 of the 48 wells completed in the Ogallala aquifer that were analyzed. Of the 48 Ogallala aquifer wells, 24 had significant upward trends for at least one water-level parameter (17 with all 3); however, 10 wells had statistically significant downward trends for at least one water-level parameter (8 with all 3 parameters). All but one of the wells with significant downward trends are located in the south-central part of the study area where licensed irrigation withdrawals are concentrated.

South Dakota↗

Diffuse-flow conceptualization and simulation of the Edwards aquifer, San Antonio region, Texas

A numerical ground-water-flow model (hereinafter, the conduit-flow Edwards aquifer model) of the karstic Edwards aquifer in south-central Texas was developed for a previous study on the basis of a conceptualization emphasizing conduit development and conduit flow, and included simulating conduits as one-cell-wide, continuously connected features. Uncertainties regarding the degree to which conduits pervade the Edwards aquifer and influence ground-water flow, as well as other uncertainties inherent in simulating conduits, raised the question of whether a model based on the conduit-flow conceptualization was the optimum model for the Edwards aquifer. Accordingly, a model with an alternative hydraulic conductivity distribution without conduits was developed in a study conducted during 2004-05 by the U.S. Geological Survey, in cooperation with the San Antonio Water System. The hydraulic conductivity distribution for the modified Edwards aquifer model (hereinafter, the diffuse-flow Edwards aquifer model), based primarily on a conceptualization in which flow in the aquifer predominantly is through a network of numerous small fractures and openings, includes 38 zones, with hydraulic conductivities ranging from 3 to 50,000 feet per day. Revision of model input data for the diffuse-flow Edwards aquifer model was limited to changes in the simulated hydraulic conductivity distribution. The root-mean-square error for 144 target wells for the calibrated steady-state simulation for the diffuse-flow Edwards aquifer model is 20.9 feet. This error represents about 3 percent of the total head difference across the model area. The simulated springflows for Comal and San Marcos Springs for the calibrated steady-state simulation were within 2.4 and 15 percent of the median springflows for the two springs, respectively. The transient calibration period for the diffuse-flow Edwards aquifer model was 1947-2000, with 648 monthly stress periods, the same as for the conduit-flow Edwards aquifer model. The root-mean-square error for a period of drought (May-November 1956) for the calibrated transient simulation for 171 target wells is 33.4 feet, which represents about 5 percent of the total head difference across the model area. The root-mean-square error for a period of above-normal rainfall (November 1974-July 1975) for the calibrated transient simulation for 169 target wells is 25.8 feet, which represents about 4 percent of the total head difference across the model area. The root-mean-square error ranged from 6.3 to 30.4 feet in 12 target wells with long-term water-level measurements for varying periods during 1947-2000 for the calibrated transient simulation for the diffuse-flow Edwards aquifer model, and these errors represent 5.0 to 31.3 percent of the range in water-level fluctuations of each of those wells. The root-mean-square errors for the five major springs in the San Antonio segment of the aquifer for the calibrated transient simulation, as a percentage of the range of discharge fluctuations measured at the springs, varied from 7.2 percent for San Marcos Springs and 8.1 percent for Comal Springs to 28.8 percent for Leona Springs. The root-mean-square errors for hydraulic heads for the conduit-flow Edwards aquifer model are 27, 76, and 30 percent greater than those for the diffuse-flow Edwards aquifer model for the steady-state, drought, and above-normal rainfall synoptic time periods, respectively. The goodness-of-fit between measured and simulated springflows is similar for Comal, San Marcos, and Leona Springs for the diffuse-flow Edwards aquifer model and the conduit-flow Edwards aquifer model. The root-mean-square errors for Comal and Leona Springs were 15.6 and 21.3 percent less, respectively, whereas the root-mean-square error for San Marcos Springs was 3.3 percent greater for the diffuse-flow Edwards aquifer model compared to the conduit-flow Edwards aquifer model. The root-mean-square errors for San Antonio and San Pedro Springs were appreciably greater, 80.2 and 51.0 percent, respectively, for the diffuse-flow Edwards aquifer model. The simulated water budgets for the diffuse-flow Edwards aquifer model are similar to those for the conduit-flow Edwards aquifer model. Differences in percentage of total sources or discharges for a budget component are 2.0 percent or less for all budget components for the steady-state and transient simulations. The largest difference in terms of the magnitude of water budget components for the transient simulation for 1956 was a decrease of about 10,730 acre-feet per year (about 2 per-cent) in springflow for the diffuse-flow Edwards aquifer model compared to the conduit-flow Edwards aquifer model. This decrease in springflow (a water budget discharge) was largely offset by the decreased net loss of water from storage (a water budget source) of about 10,500 acre-feet per year.

Texas↗

Groundwater levels, geochemistry, and water budget of the Tsala Apopka Lake system, west-central Florida, 2004–12

Tsala Apopka Lake is a complex system of lakes and wetlands, with intervening uplands, located in Citrus County in west-central Florida. It is located within the 2,100 square mile watershed of the Withlacoochee River, which drains north and northwest towards the Gulf of Mexico. The lake system is managed by the Southwest Florida Water Management District as three distinct “pools,” which from upstream to downstream are referred to as the Floral City Pool, Inverness Pool, and Hernando Pool. Each pool contains a mixture of deep-water lakes that remain wet year round, ephemeral (seasonal) ponds and wetlands, and dry uplands. Many of the major deep-water lakes are interconnected by canals. Flow from the Withlacoochee River, when conditions allow, can be diverted into the lake system. Flow thorough the canals can be used to control the distribution of water between the three pools. Flow in the canals is controlled using structures, such as gates and weirs. Hydrogeologic units in the study area include a surficial aquifer consisting of Quaternary-age sediments, a discontinuous intermediate confining unit consisting of Miocene- and Pliocene-age sediments, and the underlying Upper Floridan aquifer, which consists of Eocene- and Oligocene-age carbonates. The fine-grained quartz sands that constitute the surficial aquifer are generally thin, typically less than 25 feet thick, within the vicinity of Tsala Apopka Lake. A thin, discontinuous, sandy clay layer forms the intermediate confining unit. The Upper Floridan aquifer is generally unconfined in the vicinity of Tsala Apopka Lake because the intermediate confining unit is discontinuous and breached by numerous karst features. In the study area, the Upper Floridan aquifer includes the upper Avon Park Formation and Ocala Limestone. The Ocala Limestone is the primary source of drinking water and spring flow in the area. The objectives of this study are to document the interaction of Tsala Apopka Lake, the surficial aquifer, and the Upper Floridan aquifer; and to estimate an annual water budget for each pool and for the entire lake system for 2004–12. The hydrologic interactions were evaluated using hydraulic head and geochemical data. Geochemical data, including major ion, isotope, and age-tracer data, were used to evaluate sources of water and to distinguish flow paths. Hydrologic connection of the surficial environment (lakes, ponds, wetlands, and the surficial aquifer) was quantified on the basis of a conceptualized annual water-budget model. The model included the change in surface water and groundwater storage, precipitation, evapotranspiration, surface-water inflow and outflow, and net groundwater exchange with the underlying Upper Floridan aquifer. The control volume for each pool extended to the base of the surficial aquifer and covered an area defined to exceed the maximum inundated area for each pool during 2004–12 by 0.5 foot. Net groundwater flow was computed as a lumped value and was either positive or negative, with a negative value indicating downward or lateral leakage from the control volume and a positive value indicating upward leakage to the control volume. The annual water budget for Tsala Apopka Lake was calculated using a combination of field observations and remotely sensed data for each of three pools and for the composite three pool area. A digital elevation model at a 5-foot grid spacing and bathymetric survey data were used to define the land-surface elevation and volume of each pool and to calculate the changes in inundated area with change in lake stage. Continuous lake-stage and groundwater-level data were used to define the change in storage for each pool. The rainfall data used in the water-budget calculations were based on daily radar reflectance data and measured rainfall from weather stations. Evapotranspiration was computed as a function of reference evapotranspiration, adjusted to actual evapotranspiration using a monthly land-cover coefficient (based on evapotranspiration measurements at stations located in representative landscapes). Surface-water inflows and outflows were determined using stage data collected at a series of streamgages installed primarily at the water-control structures. Discharge was measured under varying flow regimes and ratings were developed for the water-control structures. The discharge data collected during the study period were used to calibrate a surface-water flow model for 2004–12. Flows predicted by the model were used in the water-budget analysis. Net groundwater flow was determined as the residual term in the water-budget equation. The results of the water-budget analysis indicate that rainfall was the largest input of water to Tsala Apopka Lake, whereas evapotranspiration was the largest output. For the 2004–12 analysis period, surface-water inflow accounted for 11 percent of the inputs, net groundwater inflow accounted for 1 percent of inputs (annual periods with positive net groundwater flow were included as inputs, while annual periods with negative net groundwater flow were counted as outputs), and rainfall accounted for the remaining 88 percent. For the same period, the outputs consisted of 2 percent surface-water outflow, 12 percent net groundwater outflow, and 86 percent evapotranspiration. Net groundwater inflows and surface-water/groundwater storage were negligible during the water-budget period but could be important components of the budget in individual years. The net groundwater flow was negative (downward) for 8 out of the 9 years modeled (2004–12), indicating that the Tsala Apopka Lake study area was primarily a recharge area for the underlying Upper Floridan aquifer during this time period. Groundwater-level elevation in paired wells (adjacent wells completed in the surficial aquifer and Upper Floridan aquifer) typically was higher in the surficial aquifer than the Upper Floridan aquifer. However, hydraulic head data indicate that the surficial aquifer often has discharge potential to the surface-water system, especially in the low lying areas near the major lakes. Surficial-aquifer water levels were often higher than lake stages, especially during wet periods, which is likely an indication of aquifer-to-lake seepage in these areas. East of the major lakes, hydraulic head data were nearly equal in the surficial aquifer and Upper Floridan aquifer, which is an indication that the Upper Floridan aquifer is unconfined. Based on deuterium and oxygen stable isotope data collected in December 2011 and December 2012, there was no evidence of recharge to the Upper Floridan aquifer from the wetlands east of the major lakes; aquifer isotopic ratios did not indicate an enriched source, which is typical of lake and wetland sources. West of the major lakes, there was evidence of enriched isotopic ratios in water samples from the Upper Floridan aquifer. Differences in hydraulic head at paired wells in the surficial aquifer and Upper Floridan aquifer indicated that the surficial aquifer has the potential to recharge the Upper Floridan aquifer in the western part of the pools and west of the major lakes.

Florida↗

Conceptualization and simulation of the Edwards aquifer, San Antonio region, Texas

A new numerical ground-water-flow model (Edwards aquifer model) that incorporates important components of the latest information and plausible conceptualization of the Edwards aquifer was developed. The model includes both the San Antonio and Barton Springs segments of the Edwards aquifer in the San Antonio region, Texas, and was calibrated for steady-state (1939–46) and transient (1947–2000) conditions, excluding Travis County. Transient simulations were conducted using monthly recharge and pumpage (withdrawal) data. The model incorporates conduits simulated as continuously connected (other than being separated in eastern Uvalde and southwestern Medina Counties), one-cell-wide (1,320 feet) zones with very large hydraulic-conductivity values (as much as 300,000 feet per day). The locations of the conduits were based on a number of factors, including major potentiometric-surface troughs in the aquifer, the presence of sinking streams, geochemical information, and geologic structures (for example, faults and grabens). The simulated directions of flow in the Edwards aquifer model are most strongly influenced by the presence of simulated conduits and barrier faults. The simulated flow in the Edwards aquifer is influenced by the locations of the simulated conduits, which tend to facilitate flow. The simulated subregional flow directions generally are toward the nearest conduit and subsequently along the conduits from the recharge zone into the confined zone and toward the major springs. Structures simulated in the Edwards aquifer model influencing groundwater flow that tend to restrict flow are barrier faults. The influence of simulated barrier faults on flow directions is most evident in northern Medina County. A water budget is an accounting of inflow to, outflow from, and storage change in the aquifer. For the Edwards aquifer model steady-state simulation, recharge (from seepage losses from streams and infiltration of rainfall) accounts for 93.5 percent of the sources of water to the Edwards aquifer, and inflow through the northern and northwestern model boundaries contributes 6.5 percent. The largest discharges are springflow (73.7 percent) and ground-water withdrawals by wells (25.7 percent). The principal source of water to the Edwards aquifer for the Edwards aquifer model transient simulation was recharge, constituting about 60 percent of the sources of water (excluding change in storage) to the Edwards aquifer during 1956, a drought period, and about 97 percent of the sources (excluding change in storage) during 1975, a period of above-normal rainfall and recharge. The principal discharges from the Edwards aquifer for the transient simulation were springflow and withdrawals by wells. During 1956, representing drought conditions, the change in storage (net water released from storage) was much greater than recharge, accounting for 75.9 percent of the total flow compared to 14.5 percent for recharge. Conversely, during 1975, representing above-normal rainfall and recharge conditions, recharge constituted 79.9 percent of the total flow, compared to 7.1 percent for the change in storage (net water added to storage). A series of sensitivity tests was made to ascertain how the model results were affected by variations greater than and less than the calibrated values of input data. Simulated hydraulic heads in the Edwards aquifer model were most sensitive to recharge, withdrawals, hydraulic conductivity of the conduit segments, and specific yield and were comparatively insensitive to spring-orifice conductance, northern boundary inflow, and specific storage. Simulated springflow in the Edwards aquifer model was most sensitive to recharge, withdrawals, hydraulic conductivity of the conduit segments, specific yield, and increases in northern boundary inflow and was comparatively insensitive to spring-orifice conductance and specific storage.

Texas↗

Hydrologic data and groundwater-flow simulations in the Brown Ditch Watershed, Indiana Dunes National Lakeshore, near Beverly Shores and Town of Pines, Indiana

The U.S. Geological Survey (USGS) collected data and simulated groundwater flow to increase understanding of the hydrology and the effects of drainage alterations on the water table in the vicinity of Great Marsh, near Beverly Shores and Town of Pines, Indiana. Prior land-management practices have modified drainage and caused changes in the distribution of open water, streams and ditches, and groundwater abundance and flow paths. Collected hydrologic data indicate that the majority of water entering Great Marsh flows from the southern dune ridge beneath Town of Pines, Indiana. Groundwater flow is intercepted by Brown Ditch in the eastern portion of the study area and Derby Ditch in the western portion of the study area. A smaller amount of groundwater from the northern dune ridge beneath Beverly Shores also contributed water to Great Marsh. Continuous groundwater-level data collected indicate that the predominant north-south groundwater-flow gradients vary during the course of the year due to increased levels of precipitation or during periods of drainage obstructions. Continuous surface-water discharge and surface-water elevation were measured at three USGS streamgages, one each on Brown, Kintzele and Derby Ditches. The monthly mean discharge statistics indicate that during the period of record— June 2012 to September 2013—streamflow in Kintzele Ditch was lowest during July 2012 and highest during April 2013. In Derby Ditch, streamflow also was lowest during July 2012 and highest during April 2013. Periods of relatively high and low groundwater levels during August 1982, March 2013, and April 2014 were examined and simulated by using MODFLOW and companion software. Results from the simulation of conditions during March 2013 include that nearly 100 percent of all water entering the area simulating Town of Pines is from recharge. Of all the water simulated to enter the eastern and western portions of Great Marsh, nearly 20 and 18 percent, respectively, flows from Town of Pines to the western and eastern portions of Great Marsh. The dune ridges beneath Town of Pines and to a lesser extent beneath Beverly Shores are a major source of recharge to the surficial aquifer and Great Marsh. Results from the simulation of the conditions of April 2014 include that, despite increases in the amount of water entering Great Marsh due to a beaver-dam-modified hydrologic condition, there is still virtually zero simulated groundwater flow from Great Marsh to Town of Pines. The volume of water simulated to be entering the zone representing Beverly Shores decreased by 0.43 cubic foot per second from the results of the March 2013 simulation. This simulated difference in water budgets can be attributed to increased simulated recharge in Great Marsh and Town of Pines. Effects of the inclusion of the beaver dam included the increase of the simulated water table and simulated inundated area upstream of the beaver dam due to the effects of ponding surface water. Results from the simulation scenario that includes six proposed pool-riffle control structures in Brown Ditch under the hydrologic conditions of March 2013 indicate areas inundated by water are larger, including areas just to the north of the entrance of Brown Ditch into Great Marsh, and areas north of the confluence of Brown and Kintzele Ditches. Results from the scenario simulating the increase of the Lake Michigan water level to the historical high of May 31, 1998, showed inundated areas of Great Marsh south of Beverly Shores enlarged on both sides of Lakeshore County Road with the greatest enlargement simulated to be southeast of the intersection of Lakeshore County Road and Beverly Drive. For the scenario simulating the decrease of the Lake Michigan water level to the historical low of December 23, 2007, results show little change from the original March 2013 inundated area. The results of this study can be used by water-resource managers to understand how surrounding ditches affect water levels in Great Marsh and other inland wetlands and residential areas. The groundwater model developed can be applied to answer questions about how alterations to the drainage system in the area affects water levels in the public and residential areas surrounding Great Marsh. The modeling methods developed in this study provide a template for other studies of groundwater flow and groundwater/surface-water interactions within the shallow surficial aquifer in northern Indiana, and in similar hydrologic settings that include surficial sand aquifers in coastal areas.

Indiana↗

Use of set blanks in reporting pesticide results at the U.S. Geological Survey National Water Quality Laboratory, 2001-15

Executive Summary Background .—Pesticide results from the U.S. Geological Survey (USGS) National Water Quality Laboratory (NWQL) are used for water-quality assessments by many agencies and organizations. The USGS is committed to providing data of the highest possible quality to the consumers of its data. A cooperator’s inquiries about specific pesticide detections in water revealed potential laboratory contamination issues for some results. Consequently, the USGS conducted an extensive evaluation of potential low-level contamination related to processing or analysis of water-quality samples at NWQL for 21 pesticide compounds of interest to the cooperator. This is the most comprehensive study of NWQL pesticide quality-control (QC) results to date. Purpose and scope .—The purpose of this study was to document protocols used by the NWQL to censor pesticide results and to determine the effects of laboratory contamination—as determined from detections in laboratory set blanks—on pesticide detections in groundwater and surface-water samples. More than 30,000 pesticide results from 113 selected batches of samples (2 percent or less of total batches) analyzed by the NWQL during the 15 years from 2001 to 2015 were reviewed. All laboratory results from the selected batches, including results from environmental (surface water and groundwater) and QC (set-blank, blind-blank, and blind-spike) samples, were evaluated. The study includes results for 21 pesticide compounds analyzed in groundwater and surface-water samples collected across the United States. Eleven pesticide compounds were analyzed by a gas chromatography/mass spectrometry method and 10 compounds by a liquid chromatography/mass spectrometry method. Objectives and methods .—The objectives of this study were to (1) determine the characteristics of laboratory contamination over time, (2) compare distributions of pesticide results in set blanks with distributions in environmental samples, (3) evaluate the potential for false-positive and false-negative reporting of results, and (4) evaluate the effects of reevaluating historical pesticide results using 2017 compound identification protocols on detections of pesticides in groundwater and surface-water samples. The 113 instrument batches selected for this study contained detections of one or more of the 21 pesticide compounds in set blanks or were among those batches with the highest pesticide detection frequencies in set blanks. As a result, the dataset for this study was targeted toward pesticides and batches with laboratory contamination. The objectives were addressed by statistically comparing environmental and set-blank results; computing moving averages of set-blank detection frequencies to identify periods of episodic contamination; and using summary statistics, tabular summaries, and graphical approaches, such as time-series plots and cumulative distribution functions. Results .— Objective 1 : Laboratory contamination, as determined by pesticide detections in set blanks, was found in 13 percent of set-blank results from the 113 targeted batches included in this study (as compared to 6 percent of set-blank results from all 7,620 batches analyzed during the study period). It is estimated that 92 percent of the laboratory contamination during the study period was episodic, meaning that it occurred during discrete periods of time. All 21 of the targeted pesticide compounds had periods of episodic contamination, with most episodes ranging in duration from about 1 to 8 months. The remaining 8 percent of laboratory contamination was random or from a known source (deterministic). Objective 2 : For some compounds, graphs of cumulative distribution functions of the entire distributions of set-blank and environmental samples overlap, suggesting that there is no difference in the distributions of the two types of samples. However, time-series graphs show that detections in set blanks often occur at different times (sometimes separated by years) than detections in environmental samples, indicating clear differences in those distributions, and indicating the importance of evaluating the timing of detections in all sample types. For most compounds detected in set-blank and environmental samples, detection frequencies were significantly greater in set blanks than in groundwater or surface-water samples ( p <0.05). There are several explanations for this finding, including that the 113 batches of samples chosen for this study targeted batches with detections in set blanks or that detections in set-blank samples were historically determined with less stringent identification criteria than for environmental samples (groundwater and surface-water samples). Objective 3 : The false-positive and false-negative rates from blind samples submitted during the study period by the USGS Quality Systems Branch generally were less than 1 and 5 percent, respectively, for the 21 pesticides. The only compound with a false-positive rate greater than 1 percent was flumetsulam (2.6 percent), indicating that there is a higher likelihood of flumetsulam being reported as a detection when it is not present in an environmental sample compared with the reporting of other compounds. Objective 4 : Altogether, for data in targeted batches, NWQL would have reported 0.1 percent of results from groundwater samples and 1.4 percent of results from surface-water samples differently if 2017 identification protocols were applied to historical pesticide results. In most of these cases, detections observed in historical results would change to nondetections. The small percentages of changes that would occur if historical data were reevaluated indicate that historical protocols used by the NWQL to identify detections in environmental samples were robust and produced results that are predominantly consistent with current [2017] practices. Conclusions .—The NWQL produces high-quality pesticide results at environmentally relevant concentrations. NWQL identification protocols and censoring practices are largely effective at minimizing the reporting of false-positive and false-negative results. Laboratory contamination, when it occurred, tended to occur in episodes; thus, evaluating the timing and magnitude of detections in set blanks relative to detections in environmental samples was determined to be an important consideration for analysis of environmental results. Because NWQL censoring practices do not address all types and occurrences of laboratory contamination, options for additional censoring practices are provided for data users with more specific or stringent data-quality objectives. The methods used to analyze the 21 compounds for this report can similarly be applied to all 173 pesticide compounds that were analyzed by the NWQL during the same time period. This study also has helped to identify potential improvements in reporting USGS data, such as conducting more frequent review of set-blank datasets.

Scientific Investigations Report↗

Evaluation of Street Sweeping as a Stormwater-Quality-Management Tool in Three Residential Basins in Madison, Wisconsin

Recent technological improvements have increased the ability of street sweepers to remove sediment and other debris from street surfaces; the effect of these technological advancements on stormwater quality is largely unknown. The U.S. Geological Survey, in cooperation with the City of Madison and the Wisconsin Department of Natural Resources, evaluated three street-sweeper technologies from 2002 through 2006. Regenerative-air, vacuum-assist, and mechanical-broom street sweepers were operated on a frequency of once per week (high frequency) in separate residential basins in Madison, Wis., to measure each sweeper's ability to not only reduce street-dirt yield but also improve the quality of stormwater runoff. A second mechanical-broom sweeper operating on a frequency of once per month (low frequency) was also evaluated to measure reductions in street-dirt yield only. A paired-basin study design was used to compare street-dirt and stormwater-quality samples during a calibration (no sweeping) and a treatment period (weekly sweeping). The basis of this paired-basin approach is that the relation between paired street-dirt and stormwater-quality loads for the control and tests basins is constant until a major change is made at one of the basins. At that time, a new relation will develop. Changes in either street-dirt and/or stormwater quality as a result of street sweeping could then be quantified by use of statistical tests. Street-dirt samples collected weekly during the calibration period and twice per week during the treatment period, once before and once after sweeping, were dried and separated into seven particle-size fractions ranging from less than 63 micrometers to greater than 2 millimeters. Street-dirt yield evaluation was based on a computed mass per unit length of pounds per curb-mile. An analysis of covariance was used to measure the significance of the effect of street sweeping at the end of the treatment period and to quantify any reduction in street-dirt yield. Both the regenerative-air and vacuum-assist sweepers produced reductions in street-dirt yield at the 5-percent significance level. Street-dirt yield was reduced by an average of 76, 63, and 20 percent in the regenerative-air, vacuum-assist, and high-frequency broom basins, respectively. The low-frequency broom basin showed no significant reductions in street-dirt yield. Sand-size particles (greater than 63 micrometers) recorded the greatest overall reduction. Street-sweeper pickup efficiency was determined by computing the difference between weekly street-dirt yields before and after sweeping cleaning. The regenerative-air and vacuum-assist sweepers had similar pickup efficiencies of 25 and 30 percent, respectively. The mechanical broom sweeper operating at high frequency was considerably less efficient, removing an average of 5 percent of street-dirt yield. The effects of street sweeping on stormwater quality were evaluated by use of statistical tests to compare event mean concentrations and loads computed for individual storms at the control and test basins. Loads were computed by multiplying the event mean concentrations by storm-runoff volumes. Only ammonia-nitrogen for the test basin with the vacuum-assist sweeper showed significant load increases over the control basin, at the 10-percent significance level, of 63 percent. Difficulty in detecting significant changes in constituent stormwater-quality loads could be due, in part, to the large amount of variability in the data. Coefficients of variation for the majority of constituent loads were greater than 1, indicating substantial variability. The ability to detect changes in constituent stormwater-quality loads was likely hampered by an inadequate number of samples in the data set. However, sediment transport in the storm-sewer pipe, sediment washing onto the street from other source areas, winter sand application, and sampling challenges were additional sources of variability within each study ba

Scientific Investigations Report↗

Hydrologic drought of water year 2011 compared to four major drought periods of the 20th century in Oklahoma

Water year 2011 (October 1, 2010, through September 30, 2011) was a year of hydrologic drought (based on streamflow) in Oklahoma and the second-driest year to date (based on precipitation) since 1925. Drought conditions worsened substantially in the summer, with the highest monthly average temperature record for all States being broken by Oklahoma in July (89.1 degrees Fahrenheit), June being the second hottest and August being the hottest on record for those months for the State since 1895. Drought conditions continued into the fall, with all of the State continuing to be in severe to exceptional drought through the end of September. In addition to effects on streamflow and reservoirs, the 2011 drought increased damage from wildfires, led to declarations of states of emergency, water-use restrictions, and outdoor burning bans; caused at least $2 billion of losses in the agricultural sector and higher prices for food and other agricultural products; caused losses of tourism and wildlife; reduced hydropower generation; and lowered groundwater levels in State aquifers. The U.S. Geological Survey, in cooperation with the Oklahoma Water Resources Board, conducted an investigation to compare the severity of the 2011 drought with four previous major hydrologic drought periods during the 20th century – water years 1929–41, 1952–56, 1961–72, and 1976–81. The period of water years 1925–2011 was selected as the period of record because few continuous record streamflow-gaging stations existed before 1925, and gaps in time existed where no streamflow-gaging stations were operated before 1925. In water year 2011, statewide annual precipitation was the 2d lowest, statewide annual streamflow was 16th lowest, and statewide annual runoff was 42d lowest of those 87 years of record. Annual area-averaged precipitation totals by the nine National Weather Service climate divisions from water year 2011 were compared to those during four previous major hydrologic drought periods to show how precipitation deficits in Oklahoma varied by region. The nine climate divisions in Oklahoma had precipitation in water year 2011 ranging from 43 to 76 percent of normal annual precipitation, with the Northeast Climate Division having the closest to normal precipitation and the Southwest Climate Division having the greatest percentage of annual deficit. Based on precipitation amounts, water year 2011 ranked as the second driest of the 1925–2011 period, being exceeded only in one year of the 1952 to 1956 drought period. Regional streamflow patterns for water year 2011 indicate that streamflow in the Arkansas-White-Red water resources region, which includes all of Oklahoma, was relatively large, being only the 26th lowest since 1930, primarily because of normal or above-normal streamflow in the northern part of the region. Twelve long-term streamflow-gaging stations with periods of record ranging from 67 to 83 years were selected to show how streamflow deficits varied by region in Oklahoma. Statewide, streamflow in water year 2011 was greater than streamflows measured in years during the drought periods of 1929–41, 1952–56, 1961–72, and 1976–81. The hydrologic drought worsened going from the northeast toward the southwest in Oklahoma, ranging from 140 percent (above normal streamflow) in the northeast, to 13 percent of normal streamflow in southwestern Oklahoma. The relatively low streamflow in 2011 resulted in 83.3 percent of the statewide conservation storage being available at the end of the water year in major reservoirs, similar to conservation storage in the preceding severe drought year of 2006. The ranking of streamflow as the 16th smallest for the 1925–2011 period, despite precipitation being ranked the 2d smallest, may have been caused, in part, by the relatively large streamflow in northeastern Oklahoma during water year 2011.

Oklahoma↗

Changes in phosphorus concentrations and loads in the Assabet River, Massachusetts, October 2008 through April 2014

Treated effluent discharged from municipal wastewater-treatment plants to the Assabet River in central Massachusetts includes phosphorus, which leads to increased growth of nuisance aquatic plants that decrease the river’s water quality and aesthetics in impounded reaches during the growing season. To improve the river’s water quality and aesthetics, the U.S. Environmental Protection Agency approved a total maximum daily load for phosphorus in 2004 that directed the wastewater-treatment plants to reduce the amount of total phosphorus discharged to the river by 2012. The permitted total phosphorus monthly average of 0.75 milligrams per liter during the aquatic plant growing season (April 1 through October 31) was reduced by the total maximum daily load to a target of 0.1 milligrams per liter by 2012, and the nongrowing-season limit was unchanged at 1.0 milligrams per liter. From October 2008 through April 2014, the U.S. Geological Survey, in cooperation with the Massachusetts Department of Environmental Protection, measured streamflow and collected weekly flow-proportional, composite samples of water from the Assabet River for analysis of concentrations of total phosphorus and orthophosphate. Streamflow and concentration data were used to estimate total phosphorus and orthophosphate loads in the river. The purpose of this monitoring effort was to evaluate phosphorus concentrations and loads in the river before, during, and after the wastewater-treatment-plant upgrades and to assess the effects of seasonal differences in permitted discharges. The locations of water-quality-monitoring stations, with respect to the Hudson and Ben Smith impoundments, enabled examination of effects of phosphorus entering and leaving the impoundments. Annual median concentrations of total phosphorus in wastewater-treatment plants were reduced by more than 80 percent with the plant upgrades. Measured instream annual median concentrations of total phosphorus in the Assabet River decreased by about 38 to 50 percent at three of the four monitoring stations following the wastewater-treatment-plant upgrades. At the station farthest upstream, the median total phosphorus concentration remained unchanged throughout the study; this may be attributed to the site location and potential resuspension of particulate organic matter during periods of increased streamflow. Annual median loads from the wastewater-treatment plants were reduced by up to 91 percent following the upgrades, instream annual median total phosphorus loads at the three downstream stations decreased by 71 to 76 percent, and instream orthophosphate loads at the three downstream stations decreased by 79 to 87 percent. Seasonal fluctuations (growing versus nongrowing) of total phosphorus and orthophosphate were observed instream before the upgrades. However, after the upgrades, fluctuations in phosphorus released from the treatment plants were slight and seasonal changes were typically not observed instream. Annual loads entering and leaving the two impoundments were inconclusive in determining whether the impoundments were sources or sinks of total phosphorus during the study. Total phosphorus loads entering the Hudson impoundment were consistently greater than those leaving; however, there was uncertainty about the loads at the monitoring station upstream from this impoundment. At the Ben Smith impoundment, total phosphorus and orthophosphate loads downstream were slightly greater than those upstream from the impoundment, but the differences may reflect additions from tributaries and overland runoff. Estimated instream total phosphorus concentrations and loads indicated that the decreases in total phosphorus in wastewater-treatment-plant discharges were accompanied by reductions measured in the Assabet River. A statistical analysis which incorporates the effect of varying flow conditions demonstrated significant reductions in total phosphorus concentrations after the wastewater-treatment-plant upgrades at three of the four instream monitoring stations. No significant change was observed at the most upstream location, the Assabet River at Port Street at Hudson, Massachusetts (station number 01096835), which may have been affected by flow-related resuspension of particulate phosphorus.

Massachusetts↗