USGS Science⌕ Search

SEARCH · USGS Science

Results for “Cooperator Report”

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 181 records · Page 10Linked to original sources

Groundwater quality for 75 domestic wells in Lycoming County, Pennsylvania, 2014

Groundwater is a major source of drinking water in Lycoming County and adjacent counties in north-central and northeastern Pennsylvania, which are largely forested and rural and are currently undergoing development for hydrocarbon gases. Water-quality data are needed for assessing the natural characteristics of the groundwater resource and the potential effects from energy and mineral extraction, timber harvesting, agriculture, sewage and septic systems, and other human influences. This report, prepared in cooperation with Lycoming County, presents analytical data for groundwater samples from 75 domestic wells sampled throughout Lycoming County in June, July, and August 2014. The samples were collected using existing pumps and plumbing prior to any treatment and analyzed for physical and chemical characteristics, including nutrients, major ions, metals and trace elements, volatile organic compounds, gross-alpha particle and gross beta-particle activity, uranium, and dissolved gases, including methane and radon-222. Results indicate groundwater quality generally met most drinking-water standards, but that some samples exceeded primary or secondary maximum contaminant levels (MCLs) for arsenic, iron, manganese, total dissolved solids (TDS), chloride, pH, bacteria, or radon-222. Arsenic concentrations were higher than the MCL of 10 micrograms per liter (µg/L) in 9 of the 75 (12 percent) well-water samples, with concentrations as high as 23.6 μg/L; arsenic concentrations were higher than the health advisory level (HAL) of 2 μg/L in 23 samples (31 percent). Total iron concentrations exceeded the secondary maximum contaminant level (SMCL) of 300 μg/L in 20 of the 75 samples. Total manganese concentrations exceeded the SMCL of 50 μg/L in 20 samples and the HAL of 300 μg/L in 2 of those samples. Three samples had chloride concentrations that exceeded the SMCL of 250 milligrams per liter (mg/L); two of those samples exceeded the SMCL of 500 mg/L for TDS. The pH ranged from 5.3 to 9.15 and did not meet the SMCL range of 6.5 to 8.5 in 22 samples, with 17 samples having a pH less than 6.5 and 8 samples having pH greater than 8.5. Generally, the samples that had elevated TDS, chloride, or arsenic concentrations had high pH. Total coliform bacteria were detected in 39 of 75 samples (52 percent), with Escherichia coli detected in 10 of those 39 samples. Radon-222 activities ranged from non-detect to 7,420 picocuries per liter (pCi/L), with a median of 863 pCi/L, and exceeded the proposed drinking-water standard of 300 pCi/L in 50 (67 percent) of the 75 samples; radon-222 activities were higher than the alternative proposed standard of 4,000 pCi/L in 3 samples. Water from 15 of 75 (20 percent) wells had concentrations of methane greater than the reporting level of 0.01 mg/L; detectable methane concentrations ranged from 0.04 to 16.8 mg/L. Two samples had methane concentrations (13.1 and 16.8 mg/L) exceeding the action level of 7 mg/L. Low levels of ethane (up to 0.12 mg/L) were present in the five samples with the highest methane concentrations (near or above 1 mg/L) that were analyzed for hydrocarbon compounds and isotopic composition. The isotopic composition of methane in four of these groundwater samples, from the Catskill and Lock Haven Formations and the Hamilton Group, have sample carbon isotopic ratio delta values (carbon-13/carbon-12) ranging from –42.36 to –36.08 parts per thousand (‰) and hydrogen isotopic ratio delta values (deuterium/protium) ranging from –212.0 to –188.4 ‰, which are consistent with the isotopic compositions reported for mud-gas logging samples from these geologic units and a thermogenic source of the methane. However, the isotopic composition and ratios of methane to ethane in a fifth sample indicate the methane in that sample may be of microbial origin that subsequently underwent oxidation. The fifth sample had the highest concentration of methane, 16.8 mg/L, with an carbon isotopic ratio delta values of -50.59 ‰ and a hydrogen isotopic ratio delta values of -209.7 ‰. The six well-water samples with the highest methane concentrations also had among the highest pH values (8.25 to 9.15) and elevated concentrations of sodium, lithium, boron, fluoride, arsenic, and bromide. Relatively elevated concentrations of some other constituents, such as barium, strontium, and chloride, commonly were present in, but not limited to, those well-water samples with elevated methane. Three of the six groundwater samples with the highest methane concentrations had chloride/bromide ratios that indicate mixing with a small amount of brine (0.02 percent or less) similar in composition to those reported at undetermined depth below the freshwater aquifer and for gas and oil well brines in Pennsylvania. The sample with the highest methane concentration and most other samples with low methane concentrations (less than about 1 mg/L) have chloride/bromide ratios that indicate predominantly anthropogenic sources of chloride, such as road-deicing salt, septic systems, and (or) animal waste. Brines that are naturally present may originate from deeper parts of the aquifer system, while anthropogenic sources are more likely to affect shallow groundwater because they occur on or near the land-surface. The spatial distribution of groundwater compositions generally indicate that (1) uplands along the western border of Lycoming County usually have dilute, slightly acidic oxygenated, calcium-bicarbonate type waters; (2) intermediate altitudes or areas of carbonate bedrock usually have water of near neutral pH, with highest amounts of hardness (calcium and magnesium); (3) stream valleys, low elevations where groundwater may be discharging, and deep wells in uplands usually have water with pH values greater than 8 and highest arsenic, sodium, lithium, bromide concentrations. Geochemical modeling indicated that for samples with elevated pH, sodium, lithium, bromide, and alkalinity, the water chemistry could have resulted by dissolution of calcite (calcium carbonate) combined with cation-exchange and mixing with a small amount of brine. Through cation-exchange reactions between water and bedrock, which are equivalent to processes in a water softener, calcium ions released by calcite dissolution are exchanged for sodium ions on clay minerals. Thus, the assessment of groundwater quality in Lycoming County indicates groundwater is generally of good quality, but various parts of Lycoming County can have groundwater with low to moderate concentrations of methane and other constituents that appear in naturally present brine and produced waters from gas and oil wells at high concentrations."

Pennsylvania↗

Temporal changes in nutrient concentrations in the Lower Grand River and selected drainage basins, Missouri and Iowa, during the Mississippi River Basin Healthy Watersheds Initiative (2010–23)

This report describes a cooperative study by the U.S. Geological Survey and Missouri Department of Natural Resources that evaluated temporal changes in total nitrogen (TN) and total phosphorus (TP) concentrations in the Lower Grand River hydrologic unit. The study focused on trends since 2010, when the basin was designated as a priority drainage basin of the Mississippi River Basin Healthy Watersheds Initiative (MRBI). At three local drainage basins within the Lower Grand hydrological unit (MRBI sites), stream nutrient trends were evaluated using flow-adjusted (FA) TN and TP concentrations for water years 2011 through 2023. FATN concentration trends were not statistically significant for any MRBI site. One site (site 2) showed a statistically significant increasing trend in FATP concentration, indicating a possible increase in phosphorus sources in parts of the basin. Overall, streamflow variability appeared to be the dominant factor affecting nutrient concentrations at MRBI sites. At five regional drainage basins, including the Grand River and nearby rivers with data from 1994 through 2023 (long-term sites), annual flow-normalized (FN) TN and TP concentrations were evaluated for trends before (water years 2000–10) and during (water years 2010–23) the MRBI. For water years 2010 through 2023, annual FNTN and FNTP concentrations decreased in the Grand River, as well as in the Nodaway and Chariton Rivers, which were not targeted by the MRBI. The Grand River (site 9) reversed from increasing to decreasing FNTP concentrations after 2010, with a 26-percent reduction. Annual FNTN and FNTP concentrations also decreased at the Missouri River sites. While nutrient reductions in the Grand River may reflect the effects of implemented conservation practices, similar trends in nearby, nontargeted rivers and the absence of strong decreasing trends at MRBI sites suggest that broader regional factors, instead of or in addition to MRBI efforts, may have contributed to nutrient reductions in the Grand River.

Iowa, Missouri↗

A review of surface energy balance models for estimating actual evapotranspiration with remote sensing at high spatiotemporal resolution over large extents

Many approaches have been developed for measuring or estimating actual evapotranspiration ( ET a ), and research over many years has led to the development of remote sensing methods that are reliably reproducible and effective in estimating ET a . Several remote sensing methods can be used to estimate ET a at the high spatial resolution of agricultural fields and the large extent of river basins. More complex remote sensing methods apply an analytical approach to ET a estimation using physically based models of varied complexity that require a combination of ground-based and remote sensing data, and are grounded in the theory behind the surface energy balance model. This report, funded through cooperation with the International Joint Commission, provides an overview of selected remote sensing methods used for estimating water consumed through ET a and focuses on Mapping Evapotranspiration at High Resolution with Internalized Calibration (METRIC) and Operational Simplified Surface Energy Balance (SSEBop), two energy balance models for estimating ET a that are currently applied successfully in the United States. The METRIC model can produce maps of ET a at high spatial resolution (30 meters using Landsat data) for specific areas smaller than several hundred square kilometers in extent, an improvement in practice over methods used more generally at larger scales. Many studies validating METRIC estimates of ET a against measurements from lysimeters have shown model accuracies on daily to seasonal time scales ranging from 85 to 95 percent. The METRIC model is accurate, but the greater complexity of METRIC results in greater data requirements, and the internalized calibration of METRIC leads to greater skill required for implementation. In contrast, SSEBop is a simpler model, having reduced data requirements and greater ease of implementation without a substantial loss of accuracy in estimating ET a . The SSEBop model has been used to produce maps of ET a over very large extents (the conterminous United States) using lower spatial resolution (1 kilometer) Moderate Resolution Imaging Spectroradiometer (MODIS) data. Model accuracies ranging from 80 to 95 percent on daily to annual time scales have been shown in numerous studies that validated ET a estimates from SSEBop against eddy covariance measurements. The METRIC and SSEBop models can incorporate low and high spatial resolution data from MODIS and Landsat, but the high spatiotemporal resolution of ET a estimates using Landsat data over large extents takes immense computing power. Cloud computing is providing an opportunity for processing an increasing amount of geospatial “big data” in a decreasing period of time. For example, Google Earth Engine TM has been used to implement METRIC with automated calibration for regional-scale estimates of ET a using Landsat data. The U.S. Geological Survey also is using Google Earth Engine TM to implement SSEBop for estimating ET a in the United States at a continental scale using Landsat data.

Scientific Investigations Report↗

Water quality of the Little Arkansas River and Equus Beds Aquifer before and concurrent with large-scale artificial recharge, south-central Kansas, 1995-2012

The city of Wichita artificially recharged about 1 billion gallons of water into the Equus Beds aquifer during 2007–2012 as part of Phase I recharge of the Artificial Storage and Recovery project. This report, prepared in cooperation by the U.S. Geological Survey and the city of Wichita, Kansas, summarizes Little Arkansas River (source-water for artificial recharge) and Equus Beds aquifer water quality before (1995–2006) and during (2007–2012) Artificial Storage and Recovery Phase I recharge. Additionally, aquifer water-quality distribution maps are presented and water-quality changes associated with Phase I recharge timing are described. Computed chloride concentrations in the Little Arkansas River exceeded the Federal secondary maximum contaminant level (SMCL) about 20 percent of the time during 1999 through 2012, primarily during low-flow conditions. Groundwater chloride concentrations during 2001 through 2012 exceeded the SMCL in about 6 percent of shallow wells and 7 percent of deep wells, primarily near Burrton, Kansas and along the Arkansas River. Nearly all surface water nitrate plus nitrite concentrations during 1995 through 2012 were less than the Federal maximum contaminant level (MCL); groundwater nitrate plus nitrite concentrations exceeded the MCL in about 16 percent of shallow groundwater samples and were minimal in the deeper parts of the aquifer. Several trace elements frequently exceeded drinking water criteria, including arsenic, iron, and manganese. Recharge activities at Phase I recharge wells have not resulted in substantial effects on groundwater quality in the area, likely because the total amount of water recharged is relatively small (1 billion gallons) compared to aquifer storage volume (greater than 990 billion gallons in winter 2012). The eastward movement of the Burrton chloride plume is likely being slowed by a line of recharge locations associated with Phase I; however, chloride concentrations in deep groundwater still advanced to less than one half mile from the central part of the study area. Water-quality constituents of concern (major ions, nutrients, trace elements, triazine herbicides, and fecal indicator bacteria) have not increased substantially and are likely more affected by climatological (natural recharge by precipitation) and natural (geochemical oxidation/reduction, metabolic and decay rates) processes than artificial recharge. Arsenic remains a water-quality constituent of concern because of natural and continued persistence of concentrations exceeding the Federal maximum contaminant level of 10 micrograms per liter, especially in the deeper parts of the Equus Beds aquifer.

Kansas↗

Simulating the effects of climate-related changes to air temperature and precipitation on streamflow and water temperature in the Meduxnekeag River watershed, Maine

Responsible stewardship of native fish populations and riparian plants in the Meduxnekeag River watershed in northeastern Maine is a high priority for the Houlton Band of Maliseet Indians. Understanding the potential changes in hydrology and water temperature as a result of climate change is important to this priority for evaluating future habitat conditions in the watershed. This report, prepared in cooperation with the Houlton Band of Maliseet Indians, documents and presents the results of a model using the Precipitation-Runoff Modeling System (PRMS), a hydrologic model designed to provide streamflow and water temperature simulations under predicted changes in precipitation and air temperature during the next century. To estimate streamflows and water temperature in the Meduxnekeag River watershed, a PRMS model was developed and calibrated. By using the calibrated PRMS model, simulations were made for projected scenarios of 0, 5, 10, and 15 percent increases in precipitation and for increases in air temperature of 0.0, 3.6, 7.0, and 10.4 degrees Fahrenheit (°F). The increases in precipitation and temperature were applied to all the daily input values uniformly. These scenarios were based upon the results from 30 climate change models summarized in the National Climate Change Viewer. Streamflows and water temperatures modeled for different climate scenarios were compared with streamflows and water temperatures modeled with unadjusted climate inputs. Overall, streamflow increased with increasing precipitation and decreased with increasing air temperature. Water temperature increased with increasing air temperature. At the outlet of the studied Meduxnekeag River watershed, with both a 15 percent increase in precipitation and a 10.4 °F increase in air temperature, the mean annual streamflow increased by 17 percent from 489 cubic feet per second (ft 3 /s) to 572 ft 3 /s, and the mean annual maximum streamflow decreased by 8.3 percent from 3,870 ft 3 /s to 3,550 ft 3 /s. At the same location and under the same scenario, the mean annual water temperature increased by 17.5 percent from 47.4 °F to 55.7 °F. Significant changes in mean monthly streamflows were found with increasing air temperature. The PRMS model results showed that when air temperature was increased, there was an increase in mean monthly streamflow during the winter months and a decrease in mean monthly streamflow during the spring months. In addition, with a 10.4 °F increase in the air temperature, the month with the greatest monthly streamflow changed from April to December. In addition, the PRMS model estimated that the mean annual maximum snowpack in snow water equivalent for the watershed would decrease from 7.67 inches to 1.26 inches, and the mean annual date of the maximum snowpack would change from March 21 to January 28 with a 15 percent increase in precipitation and a 10.4 °F increase in air temperature.

Maine↗

American and Sacramento Rivers, California, erodibility measurements and model

Executive Summary A previous report by the authors described sediment sampling and drilling by the U.S. Geological Survey (USGS) beside the American and Sacramento Rivers near Sacramento, California, in support of a U.S. Army Corps of Engineers project focused on regional flood control. The drilling was performed to define lithology, extract samples for laboratory testing, and perform borehole erosion tests (BETs). The U.S. Department of Agriculture (USDA) performed jet erodibility tests (JETs) near each drilling site, and a team from Texas A&M University performed laboratory tests with an erosion function apparatus (EFA). Collectively, the effort was intended to reveal spatial variations in sediment erodibility and provide data for use in a model to simulate morphological response to a major flood. The data collected by the USGS are available in a public data release. This report, developed in cooperation with the U.S. Army Corps of Engineers, provides comparisons of the three types of measurements of the erodibility of riverbed sediments. The BET is performed in the field and reveals erodibility of sediments below the bed surface. The JET is likewise performed in the field but reveals only erodibility of exposed sediments. The EFA test is done in the laboratory and was performed on soils extracted from different depths beneath the bed surface, in many cases reconstituted for laboratory testing. Tests were performed at nominally similar locations but differed by meters to tens of meters in horizontal locations. The comparison was undertaken to investigate differences among results obtained by the individual measurement approaches and to elucidate pros and cons of each method. The critical shear stress to initiate erosion and the rate of change of erosion rate per unit increase of excess shear stress, sometimes referred to as the erosion coefficient, served as the primary basis for comparison. The three test methods in some cases resulted in order of magnitude differences in estimates of these parameters. Some differences could be attributed to variances in site location or result from testing surface sediment versus a deeper layer, but systematic differences are also evident in the results. The tests performed in the laboratory using the EFA resulted in much lower values of critical shear stress and much higher values of the erosion coefficient compared to the JET tests performed by the USDA team on surface sediments. Critical shear stress was poorly resolved in the BET results because of the limited number of results per site, but the erosion coefficients derived from BET results were systematically lower than those obtained using the EFA. A new, simplified approach is also proposed to estimate the increase in channel cross-sectional area during a large flood, given data describing the initial river cross section, riverbed erodibility parameters, and peak flood discharge and duration. The model runs until the cross section erodes to an equilibrium condition or the flood ends. Output describes the area of the cross section at the end of the simulation and the time required to reach equilibrium if it was reached within the simulated period. The model assumes unique, constant values for both the critical shear stress and the erosion coefficient and represents the fluid mechanics in a simplified way, making it of limited value for quantitative predictions. It does, however, provide an indication of which cross sections are most likely to undergo the greatest change in the design event and can be used to investigate sensitivity of erosion predictions to variability in sediment erodibility measurements.

California↗

Results of a preimpoundment water-quality study of Swatara Creek, Pennsylvania

The water quality of Swatara Creek prior to impoundment by the proposed Swatara Creek Reservoir in south-central Pennsylvania was studied from July 1981 through October 1982. The report, done in cooperation with the Pennsylvania Department of Environmental Resources (PaDER), Bureau of State Parks, presents information on existing water-quality conditions. A discussion of possible water-quality conditions in and downstream from the planned impoundment is also included. Precipitation measured near the study area at Lebanon, Pennsylvania from October 1981 through September 1982 was 8 percent below normal. Streamflow for the same period at Swatara Creek at Harper Tavern just downstream from the study area was 15 percent below the average annual flow. Swatara Creek above Highway 895 has been degraded by acid mine drainage. The main inflow to the planned impoundment has 2.1 times the discharge of Lower Little Swatara Creek--a forested and agricultural basin that is also tributary to the proposed impoundment. During the 1982 water year, 17,400 tons of suspended sediment were transported from the study area. About 46 percent of the annual load was transported during 3 days of high flow. Inflows to the planned impoundment from both Lower Little Swatara Creek and Swatara Creek above Highway 895 were poorly buffered. Measured concentrations of alkalinity and acidity were usually less than 10 mg/L (milligrams per liter) and 5 mg/L as CaCO3, respectively. The inflows contain high concentrations of nutrients and metals that would probably stratify in a reservoir. Maximum concentrations of dissolved nitrate and total phosphorus were 2.6 mg/L and 0.31 mg/L, respectively. At Lower Little Swatara Creek; these concentrations are well above those needed for growth of algae. Maximum observed concentrations for total recoverable iron, aluminum, and manganese at Swatara Creek above Highway 895 at Pine Grove were 100,000 ug/L (micrograms per liter), 66,000 ug/L, and 2,300 ug/L, respectively. Large increases in metal concentrations along with simultaneous decreases in pH and increases in acidity confirm that mine drainage continues to degrade the water quality of Swatara Creek and may have a large impact on water quality of the planned impoundment. Iron, lead, copper, and zinc concentrations periodically exceeded the U.S. Environmental Protection Agency (U.S. EPA) criteria for freshwater aquatic life. Concentrations of manganese and lead also exceeded the U.S. EPA criteria for domestic water supplies and human health, respectively. The water quality of the Swatara Creek Reservoir will depend on characteristics such as (1) the detention time of water in the lake, (2) the timing and extent of thermal and chemical stratification, (3) sedimentation, and (4) the chemical loading and concentrations in the lake. Each of these characteristics may depend in part, on streamflow. The impoundment will act as a sediment trap and thus reduce the concentrations of total phosphorus, iron, aluminum, lead, copper, and zinc immediately downstream from the impoundment. Large storm discharges and releases from the hypolimnion of the reservoir to attain the winter-pool level may contain low oxygen concentrations and elevated concentrations of iron, aluminum, lead, copper, and zinc. Unless conservation releases from the multi-level release gates are carefully controlled, low dissolved-oxygen levels and high metal concentrations may degrade the downstream water quality and be detrimental to the aquatic community.

Pennsylvania↗

Effects of historical coal mining and drainage from abandoned mines on streamflow and water quality in Bear Creek, Dauphin County, Pennsylvania — March 1999–December 2002

More than 100 years of anthracite coal mining has changed surface- and ground-water hydrology and contaminated streams draining the Southern Anthracite Coal Field in east-central Pennsylvania. Bear Creek drains the western prong of the Southern Anthracite Coal Field and is affected by metals in drainage from abandoned mines and streamwater losses. Total Maximum Daily Loads (TMDL) developed for dissolved iron of about 5 lb/d (pounds per day) commonly are exceeded in the reach downstream of mine discharges. Restoration of Bear Creek using aerobic ponds to passively remove iron in abandoned mine drainage is under consideration (2004) by the Dauphin County Conservation District. This report, prepared in cooperation with the Dauphin County Conservation District, evaluates chemical and hydrologic data collected in Bear Creek and its receiving waters prior to implementation of mine-drainage treatment. The data collected represent the type of baseline information needed for documentation of water-quality changes following passive treatment of mine drainage in Pennsylvania and in other similar hydrogeologic settings. Seven surface-water sites on Bear Creek and two mine discharges were monitored for nearly three years to characterize the chemistry and hydrology of the following: (1) Bear Creek upstream of the mine discharges (BC-UMD), (2) water draining from the Lykens-Williamstown Mine Pool at the Lykens Water-Level Tunnel (LWLT) and Lykens Drift (LD) discharges, (3) Bear Creek after mixing with the mine discharges (BC-DMD), and (4) Bear Creek prior to mixing with Wiconisco Creek (BCM). Two sites on Wiconisco Creek, upstream and downstream of Bear Creek (WC-UBC and WC-DBC, respectively), were selected to evaluate changes in streamflow and water quality upon mixing with Bear Creek. During periods of below-normal precipitation, streamwater loss was commonly 100 percent upstream of site BC-UMD (streamflow range = 0 to 9.7 ft3/s (cubic feet per second)) but no loss was detected downstream owing to sustained mine water drainage from the Lykens Water-Level Tunnel (range = 0.41 to 3.7 ft3/s), Lykens Drift (range = 0.40 to 6.1 ft3/s), and diffuse zones of seepage. Collectively, mine water inputs contributed about 84 percent of base flow and 53 percent of stormflow measured in the downstream reach. An option under consideration by the Dauphin County Conservation District for treatment of the discharge from the LWLT requires the source of the discharge to be captured and rerouted downstream, bypassing approximately 1,000 feet of stream channel. Because streamwater loss upstream of the tunnel was commonly 100 percent, rerouting the discharge from the LWLT may extend the reach of Bear Creek that is subject to dryness. Differences in the chemistry of water discharging from the LWLT compared to the LD suggest that the flow path through the Lykens-Williamstown Mine Pool to each mine discharge is unique. The LWLT is marginally alkaline (median net acid neutralizing capacity (ANC) = 9 mg/L (milligrams per liter) as CaCO3; median pH = 5.9), commonly becomes acidic (minimum net ANC = -74 mg/L as CaCO3) at low flow, and may benefit from alkaline amendments prior to passive treatment. Water discharging from the LD provides excess ANC (median net ANC = 123 mg/L as CaCO3; median pH = 6.5) to the downstream reach and is nearly anoxic at its source (median dissolved oxygen = 0.5 mg/L). Low dissolved oxygen water with relatively high ANC and metals concentrations discharging from the LD is characteristic of a deeper flow path and longer residence time within the mine pool than the more acidic, oxygenated water discharging from the LWLT. TMDLs for iron have been developed for dissolved species only. Consequently, distinguishing between dissolved and suspended iron in Bear Creek is important for evaluating water-quality improvement through TMDL attainment. Median total iron concentration increased from 550 mg/L (micrograms per liter) at site BC-UM

Pennsylvania↗

Estimating daily public supply water use by drinking water service area in New Jersey

This report, prepared in cooperation with the New Jersey Department of Environmental Protection, presents a method for estimating daily public supply water use by drinking water service area systems for New Jersey. The ability to accurately estimate daily public supply water use could help water supply planners in New Jersey better understand and manage the state’s limited water resources and balance the competing needs for freshwater resources. Data sources for this work include daily public supply water-use data from 2016 through 2020 acquired from New Jersey American Water for 15 drinking water service areas and monthly data exported from the New Jersey Department of Environmental Protection’s online water transfer data model database (known as NJWaTr). The two datasets were compared by aggregating the daily data to a monthly timescale. Statistical regression analysis was applied to the daily data, along with climate data, to evaluate what factors are influential in estimating daily fluctuations and trends in daily public supply water use. Fifteen regression equations were developed, one for each of the 15 drinking water service area systems for which daily data were acquired. Regression equations for systems that had seasonal patterns performed better than equations for non-seasonal systems. For the test year (2020), the average adjusted coefficient of determination for the linear regression with autoregressive errors model among systems with seasonality was 0.78; the average adjusted coefficient of determination for the linear regression with autoregressive errors model among systems with little or no seasonality was 0.25. The effects of anomalous data in the regression analysis were examined by comparing adjusted coefficient of determination values when the atypical data points were removed versus when they were retained in the analysis. Overall, including the anomalous data did not have a large effect on the results, and thus the data were retained for this study. In addition to developing regression equations, all 589 unique drinking water service area systems in New Jersey were characterized based on socio-economic data and monthly water-use data from NJWaTr. Systems that are located near the New Jersey coast, serve populations larger than 1,970 people, or serve areas that have median property values over $256,250 tended to demonstrate seasonal water-use behaviors. Systems that have mostly urban residential land use tended to show little to no seasonal water-use behaviors. Finally, a method was developed to disaggregate monthly data to a daily timescale and was tested against systems for which daily data were not available. Two regression equation forms were developed to be applied to systems beyond the 15 systems from which the original equations were developed; one equation was developed for use when all drinking water service area systems showed little to no seasonality, and the other equation was developed for use when systems displayed seasonal behavior. To the extent possible, uncertainty and possible sources of error were identified and examined in relation to the regression model equations developed. Additional daily data from these 15 systems (over different years) and daily data from different systems could be used to further evaluate the results of the disaggregation through a comprehensive assessment of error. Further adjustments to the regression equations could be made, ultimately enhancing their accuracy.

New Jersey↗

Numerical groundwater-flow model of the Minnelusa and Madison hydrogeologic units in the Rapid City area, South Dakota

The city of Rapid City and other water users in the Rapid City area obtain water supplies from the Minnelusa and Madison aquifers, which are contained in the Minnelusa and Madison hydrogeologic units. A numerical groundwater-flow model of the Minnelusa and Madison hydrogeologic units in the Rapid City area was developed to synthesize estimates of water-budget components and hydraulic properties, and to provide a tool to analyze the effect of additional stress on water-level altitudes within the aquifers and on discharge to springs. This report, prepared in cooperation with the city of Rapid City, documents a numerical groundwater-flow model of the Minnelusa and Madison hydrogeologic units for the 1,000-square-mile study area that includes Rapid City and the surrounding area. Water-table conditions generally exist in outcrop areas of the Minnelusa and Madison hydrogeologic units, which form generally concentric rings that surround the Precambrian core of the uplifted Black Hills. Confined conditions exist east of the water-table areas in the study area. The Minnelusa hydrogeologic unit is 375 to 800 feet (ft) thick in the study area with the more permeable upper part containing predominantly sandstone and the less permeable lower part containing more shale and limestone than the upper part. Shale units in the lower part generally impede flow between the Minnelusa hydrogeologic unit and the underlying Madison hydrogeologic unit; however, fracturing and weathering may result in hydraulic connections in some areas. The Madison hydrogeologic unit is composed of limestone and dolomite that is about 250 to 610 ft thick in the study area, and the upper part contains substantial secondary permeability from solution openings and fractures. Recharge to the Minnelusa and Madison hydrogeologic units is from streamflow loss where streams cross the outcrop and from infiltration of precipitation on the outcrops (areal recharge). MODFLOW-2000, a finite-difference groundwater-flow model, was used to simulate flow in the Minnelusa and Madison hydrogeologic units with five layers. Layer 1 represented the fractured sandstone layers in the upper 250 ft of the Minnelusa hydrogeologic unit, and layer 2 represented the lower part of the Minnelusa hydrogeologic unit. Layer 3 represented the upper 150 ft of the Madison hydrogeologic unit, and layer 4 represented the less permeable lower part. Layer 5 represented an approximation of the underlying Deadwood aquifer to simulate upward flow to the Madison hydrogeologic unit. The finite-difference grid, oriented 23 degrees counterclockwise, included 221 rows and 169 columns with a square cell size of 492.1 ft in the detailed study area that surrounded Rapid City. The northern and southern boundaries for layers 1-4 were represented as no-flow boundaries, and the boundary on the east was represented with head-dependent flow cells. Streamflow recharge was represented with specified-flow cells, and areal recharge to layers 1-4 was represented with a specified-flux boundary. Calibration of the model was accomplished by two simulations: (1) steady-state simulation of average conditions for water years 1988-97 and (2) transient simulations of water years 1988-97 divided into twenty 6-month stress periods. Flow-system components represented in the model include recharge, discharge, and hydraulic properties. The steady-state streamflow recharge rate was 42.2 cubic feet per second (ft3/s), and transient streamflow recharge rates ranged from 14.1 to 102.2 ft3/s. The steady-state areal recharge rate was 20.9 ft3/s, and transient areal recharge rates ranged from 1.1 to 98.4 ft3/s. The upward flow rate from the Deadwood aquifer to the Madison hydrogeologic unit was 6.3 ft3/s. Discharge included springflow, water use, flow to overlying units, and regional outflow. The estimated steady-state springflow of 32.8 ft3/s from seven springs was similar to the simulated springflow of 31.6 ft3/s, which included 20.5 ft3

South Dakota↗

Simulated water budgets and ground-water/surface-water interactions in Bushkill and parts of Monocacy Creek watersheds, Northampton County, Pennsylvania: A preliminary study with identification of data needs

This report, prepared in cooperation with the Department of Environmental Protection, Office of Mineral Resources Management, provides a preliminary analysis of water budgets and generalized ground-water/surface-water interactions for Bushkill and parts of Monocacy Creek watersheds in Northampton County, Pa., by use of a ground-water flow model. Bushkill Creek watershed was selected for study because it has areas of rapid growth, ground-water withdrawals from a quarry, and proposed stream-channel modifications, all of which have the potential for altering ground-water budgets and the interaction between ground water and streams. Preliminary 2-dimensional, steady-state simulations of ground-water flow by the use of MODFLOW are presented to show the status of work through September 2005 and help guide ongoing data collection in Bushkill Creek watershed. Simulations were conducted for (1) predevelopment conditions, (2) a water table lowered for quarry operations, and (3) anthropogenic changes in hydraulic conductivity of the streambed and aquifer. Preliminary results indicated under predevelopment conditions, the divide between the Bushkill and Monocacy Creek ground-water basins may not have been coincident with the topographic divide and as much as 14 percent of the ground-water discharge to Bushkill Creek may have originated from recharge in the Monocacy Creek watershed. For simulated predevelopment conditions, Schoeneck Creek and parts of Monocacy Creek were dry, but Bushkill Creek was gaining throughout all reaches. Simulated lowering of the deepest quarry sump to an altitude of 147 feet for quarry operations caused ground-water recharge and streamflow leakage to be diverted to the quarry throughout about 14 square miles and caused reaches of Bushkill and Little Bushkill Creeks to change from gaining to losing streams. Lowering the deepest quarry sump to an altitude of 100 feet caused simulated ground-water discharge to the quarry to increase about 4 cubic feet per second. Raising the deepest sump to an altitude of 200 feet caused the simulated discharge to the quarry to decrease about 14 cubic feet per second.Decreasing the hydraulic conductivity of the streambed of Bushkill Creek in the reach of large losses of flow caused simulated ground-water levels to decline and ground-water discharge to a quarry to decrease from 74 to 45 cubic feet per second. Decreasing the hydraulic conductivity of a hypothesized highly transmissive zone with a plug of relatively impermeable material caused ground-water levels to increase east of the plug and decline west of the plug, and decreased the discharge to a quarry from 74 to 53 cubic feet per second. Preliminary results of the study have significant limitations, which need to be recognized by the user. The results demonstrated the usefulness of ground-water modeling with available data sets, but as more data become available through field studies, a more complete evaluation could be conducted of the preliminary assumptions in the conceptual model, model sensitivity, and effects of boundary conditions. Additional streamflow and ground-water-level measurements would be needed to better quantify recharge and aquifer properties, particularly the anisotropy of carbonate rocks. Measurements of streamflow losses at average, steady-state hydrologic conditions could provide a more accurate estimate of ground-water recharge from this source, which directly affects water budgets and contributing areas simulated by the model.

Pennsylvania↗

Water budgets and groundwater volumes for abandoned underground mines in the Western Middle Anthracite Coalfield, Schuylkill, Columbia, and Northumberland Counties, Pennsylvania: Preliminary estimates with identification of data needs

This report, prepared in cooperation with the Pennsylvania Department of Environmental Protection (PaDEP), the Eastern Pennsylvania Coalition for Abandoned Mine Reclamation, and the Dauphin County Conservation District, provides estimates of water budgets and groundwater volumes stored in abandoned underground mines in the Western Middle Anthracite Coalfield, which encompasses an area of 120 square miles in eastern Pennsylvania. The estimates are based on preliminary simulations using a groundwater-flow model and an associated geographic information system that integrates data on the mining features, hydrogeology, and streamflow in the study area. The Mahanoy and Shamokin Creek Basins were the focus of the study because these basins exhibit extensive hydrologic effects and water-quality degradation from the abandoned mines in their headwaters in the Western Middle Anthracite Coalfield. Proposed groundwater withdrawals from the flooded parts of the mines and stream-channel modifications in selected areas have the potential for altering the distribution of groundwater and the interaction between the groundwater and streams in the area. Preliminary three-dimensional, steady-state simulations of groundwater flow by the use of MODFLOW are presented to summarize information on the exchange of groundwater among adjacent mines and to help guide the management of ongoing data collection, reclamation activities, and water-use planning. The conceptual model includes high-permeability mine voids that are connected vertically and horizontally within multicolliery units (MCUs). MCUs were identified on the basis of mine maps, locations of mine discharges, and groundwater levels in the mines measured by PaDEP. The locations and integrity of mine barriers were determined from mine maps and groundwater levels. The permeability of intact barriers is low, reflecting the hydraulic characteristics of unmined host rock and coal. A steady-state model was calibrated to measured groundwater levels and stream base flow, the latter at many locations composed primarily of discharge from mines. Automatic parameter estimation used MODFLOW-2000 with manual adjustments to constrain parameter values to realistic ranges. The calibrated model supports the conceptual model of high-permeability MCUs separated by low-permeability barriers and streamflow losses and gains associated with mine infiltration and discharge. The simulated groundwater levels illustrate low groundwater gradients within an MCU and abrupt changes in water levels between MCUs. The preliminary model results indicate that the primary result of increased pumping from the mine would be reduced discharge from the mine to streams near the pumping wells. The intact barriers limit the spatial extent of mine dewatering. Considering the simulated groundwater levels, depth of mining, and assumed bulk porosity of 11 or 40 percent for the mined seams, the water volume in storage in the mines of the Western Middle Anthracite Coalfield was estimated to range from 60 to 220 billion gallons, respectively. Details of the groundwater-level distribution and the rates of some mine discharges are not simulated well using the preliminary model. Use of the model results should be limited to evaluation of the conceptual model and its simulation using porous-media flow methods, overall water budgets for the Western Middle Anthracite Coalfield, and approximate storage volumes. Model results should not be considered accurate for detailed simulation of flow within a single MCU or individual flooded mine. Although improvements in the model calibration were possible by introducing spatial variability in permeability parameters and adjusting barrier properties, more detailed parameterizations have increased uncertainty because of the limited data set. The preliminary identification of data needs includes continuous streamflow, mine discharge rate, and groundwater levels in the mines and adjacent areas. Data collected whe

Pennsylvania↗

A workshop model simulating fate and effect of drilling muds and cuttings on benthic communities

Oil and gas exploration and production at marine sites has generated concern over potential environmental impacts resulting from the discharge of spent drilling muds and cuttings. This concern has led to a broad array of publicly and privately sponsored research. This report described a cooperative modeling effort designed to focus information resulting from this research through construction of explicit equations that simulate the potential impacts of discharge drilling fluids (muds) and cuttings on marine communities. The model is the result of collaboration among more than 30 scientists. The principal cooperating organizations were the E.S. Environmental Protection Agency, the U.S. Minerals Management Service, the Offshore Operators Committee, and the Alaska Oil and Gas Association. The overall simulation model can be conceptualized as three connected submodels: Discharge and Plume Fate, Sediment Redistribution, and Benthic Community Effects. On each day of simulation, these submodels are executed in sequence, with flows of information between submodels. The Benthic Community Effects submodel can be further divided into sections that calculate mortality due to burial, mortality due to toxicity, mortality due to resuspension disturbance, and growth of the community. The model represents a series of seven discrete 1-m 2 plots at specified distances along a transect in one direction away from a discharge point. It consists of coupled difference equations for which parameter values can easily be set to evaluate different conditions or to examine the sensitivity of output to various assumptions. Sets of parameter values were developed to represent four general cases or scenarios: (1) a shallow (5 m), cold environment with ice cover during a substantial fraction of the year, such as might be encountered in the Beaufort Sea, Alaska; (2) a shallow (20 m), temperate environment, such as might be encountered in the Gulf of Mexico; (3) a deeper (80 m), temperate environment, such as might be encountered in the Gulf of Mexico; and (4) a very deep (1,000 m) environment, such as might be encountered on the Atlantic slope. The focus of the modeling effort was on the connection of a reasonable representation of physical fate to the biological responses of populations, rather than on highly detailed representations of individual processes. For example, the calculations of physical fate are not as detailed as those in the recently published model of Brandsma et al. (1983). The value of the model described herein is in the broad scope of processes that are explicitly represented and linked together. The model cannot be considered to produce reliable predictions of the quantitative impacts of discharged drilling fluids and cuttings on biological populations at a particular site. Limitations of the model in predicting integrated fate and effects can be traced to three general areas: level of refinement of the algorithms used in the model; lack of understanding of the processes determining fate and effects; and parameter and data values. Despite the limitations, several qualitative conclusions concerning both potential impacts and the importance of various remaining data gaps can be drawn from the modeling effort. These include: (1) Simple, unequivocal conclusions about fate and effects across geographical regions and drilling operations are difficult, if not misleading, due to the large amount of variability in characteristics of discharged materials (e.g., oil content and toxicity), discharge conditions (e.g., duration of drilling operations), physical environments (e.g., water depth, current direction, and sediment disturbance regimes), and biological communities (e.g., intrinsic growth rates). Different combinations of these characteristics can result in substantial differences in simulated environmental fate and biological effects. For examples, simulated recovery in some high-energy environments occurs within months after the cessation of discharge operations, even at heavily impacted sites, whereas simulated recover in some low-energy environments takes years at heavily impacted sites. <2) Considerable difficulties remain in the reliable extrapolation of results from laboratory toxicity experiments to predictions of population effects in the field. (3) The volume of material discharged and duration of operations in the production drilling operations simulated by the model are sufficient to produce substantial simulated biological impacts at some plots, both in terms of differences from a control plot during the period of discharge operations, and in terms of the recovery period following the perturbations. Evaluation of the significance of potential effects involves the following factors: • Definition of a specific spatial and temporal reference frame (e.g., What is the natural variation? Is 1 year to be considered a "long" or "short" time? Is 50 m to be considered a "large" or "trivial" distance? • Consideration of rare or unique resources and particularly sensitive biotic assemblages. • Consideration of the potential for long term, cumulative effects. Some of these aspects are clearly beyond the scope of this modeling efforts (e.g., the model does not simulate the long term fate of resuspended material). The model does, however, contain an internal "reference frame" by comparison to simulated behavior at a control plot. The model, in general, simulates substantial "natural" variation at the reference or control plots, both over time, due to sediment disturbance events in medium to high energy environments, and over space, due to geographically varying conditions, such as water depth and current regime.

Report↗

Summary of water-resources activities of the U.S. Geological Survey in Colorado: fiscal year 1984

Water-resources investigations of the U.S. Geological Survey in Colorado consist of collecting water-resources data and conducting interpretive hydrologic investigations. The water-resources data and the results of the investigations are published or released by either the U.S. Geological Survey or by cooperating agencies. This report describes the water- resources investigations in Colorado for the 1984 fiscal year (October 1, 1983, to September 30, 1984).

Colorado↗

Water-resources investigations of the U.S. Geological Survey in Colorado, fiscal year 1977

Water-resources investigations of the U.S. Geological Survey in Colorado consist of collecting water-resources data and conducting interpretive hydro-logic investigations. The water-resources data and the results of the investigations are published or released by either the U.S. Geological Survey or by cooperating agencies. This report describes the water-resources investigations in Colorado for the 1977 fiscal year (October 1, 1976, to September 30, 1977). The U.S. Geological Survey's investigations of the water-resources of Colorado are under the direction of James E. Biesecker, District Chief. The Colorado District office is located in Building 53, Denver Federal Center, Lakewood, Colo. (fig. 1). The Colorado District has four subdistrict offices located in Grand Junction, Lakewood, Meeker, and Pueblo (fig. 1).

Colorado↗

Programs and activities of the Missouri District, Water Resources Division, U.S. Geological Survey, fiscal year 1979

Water-resources investigations of the U.S. Geological Survey in Missouri consist of collecting hydrologic data and conducting interpretive investigations. The data and the results of the investigations are published or released by either the U.S. Geological Survey or by cooperating agencies. This report describes the data-collection activities and investigations in Missouri for the 1979 fiscal year and provides an extensive list of water-resources references for the State.

Missouri↗

Water-resources investigations of the U.S. Geological Survey in Missouri : fiscal year 1980

Water-resources investigations of the U.S. Geological Survey in Missouri consist of collecting hydrologic data and conducting interpretative investigations. The data and the results of the investigations are published or released by either the U.S. Geoloogical Survey or by cooperating agencies. The report describes the data-collection activities and investigations in Missouri for the 1980 fiscal year and provides an extensive list of water-resources references for the State. (USGS)

Open-File Report↗

Water-resources investigations of the USGS in Texas, fiscal year 1981

Water-resources investigations in Texas consist of the collection of basic records through the hydrologic-data network, interpretive studies, and research projects. The hydrologic records and the results of investigations are published by the Geological Survey or by cooperating agencies. This report describes the water-resources projects and activities of the Geological Survey in Texas for the 1981 fiscal year (October 1, 1980 to September 30, 1981). The Geological Survey's investigations of the water resources of Texas are under the general direction of Jack Rawson, Acting District Chief. The Texas District office is in the Federal Building, 300 East Eighth Street, Austin, Texas 78701.

Texas↗