SEARCH · Search USGS

Results for “Environmental Monitoring and Assessment”

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

SourceResearch connectionFrom yearThrough year

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

1,343 records · Page 75Linked to original sources

USGS2010

Effects of selected low-impact-development (LID) techniques on water quality and quantity in the Ipswich River Basin, Massachusetts: Field and modeling studies

During the months of August and September, flows in the Ipswich River, Massachusetts, dramatically decrease largely due to groundwater withdrawals needed to meet increased residential and commercial water demands. In the summer, rates of groundwater recharge are lower than during the rest of the year, and water demands are higher. From 2005 to 2008, the U.S. Geological Survey, in a cooperative funding agreement with the Massachusetts Department of Conservation and Recreation, monitored small-scale installations of low-impact-development (LID) enhancements designed to diminish the effects of storm runoff on the quantity and quality of surface water and groundwater. Funding for the studies also was contributed by the U.S. Environmental Protection Agency’s Targeted Watersheds Grant Program through a financial assistance agreement with Massachusetts Department of Conservation and Recreation. The monitoring studies examined the effects of (1) replacing an impervious parking lot surface with a porous surface on groundwater quality, (2) installing rain gardens and porous pavement in a neighborhood of 3 acres on the quantity and quality of stormwater runoff, and (3) installing a 3,000-square foot (ft2) green roof on the quantity and quality of stormwater runoff. In addition, the effects of broad-scale implementation of LID techniques, reduced water withdrawals, and water-conservation measures on streamflow in large areas of the basin were simulated using the U.S. Geological Survey’s Ipswich River Basin model. From June 2005 to 2007, groundwater quality was monitored at the Silver Lake town beach parking lot in Wilmington, MA, prior to and following the replacement of the conventional, impervious-asphalt surface with a porous surface consisting primarily of porous asphalt and porous pavers. Changes in the concentrations of the water-quality constituents, phosphorus, nitrogen, cadmium, chromium, copper, lead, nickel, zinc, and total petroleum hydrocarbons, were monitored. Increased infiltration of precipitation did not result in discernible increases in concentrations of these potential groundwater contaminants. Concentrations of dissolved oxygen increased slightly in groundwater profiles following the removal of the impervious asphalt parking lot surface. In Wilmington, MA, in a 3-acre neighborhood, stormwater runoff volume and quality were monitored to determine the ability of selected LID enhancements (rain gardens and porous paving stones) to reduce flows and loads of the above constituents to Silver Lake. Flow-proportional water-quality samples were analyzed for nutrients, metals, total petroleum hydrocarbons, and total-coliform and Escherichia coli bacteria. In general, when all storms were considered, no substantial decreases were observed in runoff volume as a result of installing LID enhancements. However, the relation between rainfall and runoff did provide some insight into how the LID enhancements affected the effective impervious area for the neighborhood. A decrease in runoff was observed for storms of 0.2 inches (in.) or less of precipitation, which indicated a reduction in effective impervious area from approximately 10 percent to about 4.5 percent for the 3-acre area. Water-quality-monitoring results were inconclusive; there were no statistically significant differences in concentrations or loads when the pre- and post-installation-period samples were compared. Three factors were probably most important in minimizing differences: (1) the small decrease in effective impervious area, (2) the differences in the size of storms sampled for water-quality constituents before and after installation of the infiltration enhancing measures, and (3) small sample sizes. In a third field study, the characteristics of runoff from a vegetated “green” roof and a conventional, rubber-membrane roof were compared. The amount of precipitation and the length of the antecedent dry period were the two primary factors affecting the green roof’s water-storage capacity. The green roof retained more than 50 percent of the precipitation from storms with 0.04 to 1.0 in. of rain. Approximately 95 percent of the precipitation from one storm of nearly 2 in. was retained by the green roof. On the rubber-membrane roof, only a small, shallow puddle of insubstantial volume ever remained after a storm. Bulk precipitation from 10 storms was monitored for the same constituents (nutrients, metals, and total petroleum hydrocarbons) as the roof runoff, and the results were compared with those for roof-runoff samples. The use of fertilizers to help establish the vegetation during the study probably distorted any effect the plants and growing medium may have had on the retention of target analytes. As a result of the fertilizer and growing medium chemistry, median concentrations of total nitrogen, total phosphorus, cadmium, copper, and nickel in runoff from the green roof were greater than in the runoff from the conventional roof or in bulk precipitation. Concentrations of lead and zinc were greater in runoff from the conventional roof, probably a result of passage through the old, metal drainpipes. Simulations of the effects of LID on streamflow in the Ipswich River Basin were conducted with a previously calibrated Hydrological Simulation Program-FORTRAN (HSPF) precipitation-runoff model. Simulations were conducted at multiple spatial scales to evaluate the effects of (1) updated water withdrawals for the towns of Reading and Wilmington; (2) potential land-use changes at buildout (potential future development); (3) effective impervious area reductions upstream from the South Middleton streamgage to represent the effects of widespread implementation of LID retrofit techniques; (4) basin-scale water withdrawal reductions scaled up (expanded to the town level) from water-conservation pilot programs conducted by the Massachusetts Department of Conservation and Recreation; and (5) land-use change and LID techniques at a local scale, which is smaller than the HSPF subbasin. Effects on streamflow generally were evaluated by comparing results of two or more related simulations for selected reaches in the basin; thus, relative rather than absolute changes in simulated flow were the focus of the assessment. Simulations indicated that reduced withdrawals for the towns of Reading and Wilmington led to substantially higher medium and low flows in most of the reaches upstream from the South Middleton streamgage. Simulations of water-conservation measures resulted in negligible effects on streamflow. Overall, simulations indicated that spatial scale is an important factor in determining the effects of land-use change and LID practices on streamflow. Potential land-use changes at buildout had modest (percent differences of less than 20 percent) effects on streamflow in most subbasins because relatively little land in the basin was available for development (about 17 percent); moreover, most of the available open land is zoned for low-density residential development, and this land-use category was simulated to contain relatively little effective impervious area and to be similar hydrologically to the forested land in place prior to development. Results of the simulations conducted to evaluate widespread effective impervious area reductions upstream from the South Middleton streamgage indicated that the percentage of urban land use and associated effective impervious area was too small for a 50-percent reduction of effective impervious area to appreciably affect streamflow (percent differences of less than 20 percent) in most subbasins. In contrast, the results of the hypothetical local-scale simulations indicated that for smaller streams, where the percentage of urban land use and associated effective impervious area in the drainage area may be substantially higher, land-use change, development patterns, and LID practices potentially have much greater effects on streamflow. Modeling results also indicated that LID was potentially most beneficial for minimizing streamflow alteration when applied to dense urban development, largely because larger tracts of effective impervious area were available for reduction than were available for other land-use categories. For example, commercial-industrial-transportation land use is composed of 37 percent pervious area and 63 percent effective impervious area in the HSPF model, whereas low-density residential area is composed of 97.5 percent pervious area and only 2.5 percent effective impervious area. Field and modeling studies concurred in the assessment that LID enhancements would likely have the greatest effect on decreasing stormwater runoff when broadly applied to highly impervious urban areas. A measurable effect for small rainfall events (less than 0.25 inch) was determined in the small, highly pervious area that was monitored in this study, but the volume difference was not great.

Massachusetts↗

USGS2020

Research, monitoring, and evaluation of emerging issues and measures to recover the Snake River Fall Chinook Salmon ESU: January 2019 - December 2019

The portion of the Snake River fall Chinook salmon Oncorhynchus tshawytscha evolutionary significant unit (ESU) that spawns upstream of Lower Granite Dam transitioned from low to high abundance during 19922019 in association with U.S. Endangered Species Act recovery efforts and other federally mandated actions. This annual report focuses on (1) numeric and habitat use responses by natural- and hatchery-origin spawners, (2) phenotypic and numeric responses by natural-origin juveniles, and (3) use of a small unmanned aerial system (sUAS) to search for fall Chinook salmon redds and carcasses. Spawners have located and used most of the available spawning habitat and that habitat is gradually approaching redd capacity. Timing of spawning and fry emergence have been relatively stable, but effects of density dependence are evident in juvenile life stages. Apparent abundance of juvenile fall Chinook salmon has increased and we noted the following responses: parr dispersal from riverine rearing habitat into Lower Granite Reservoir has become earlier; growth rate (g/d) and dispersal size of parr declined; and passage timing of smolts from the two Snake River reaches has become earlier and downstream movement rate faster. These findings coupled with stock-recruitment analyses presented in this report provide evidence for density-dependence in the Snake River reaches and in Lower Granite Reservoir that was influenced by the expansion of the recovery program. The long-term goal is to use this information in a comprehensive modeling effort to conduct action-effectiveness and uncertainty research and to inform Fish Population, Hydrosystem, Harvest, Hatchery, and Predation and Invasive Species Management Research, Monitoring, and Evaluation (RM&E) progams. In 2019, the U.S. Geological Survey (USGS) shifted survey efforts in the Snake River toward deepwater redd searches and fish collection for parentage-based tagging (PBT) analyses because all unmanned aerial system (UAS) activities were suspended by the Department of the Interior two weeks into the spawning season. We counted 81 deepwater redds at 17 of the 29 sites surveyed. Redd depths averaged 3.6 m. We collected a total of 123 live fall Chinook salmon from 16 unique geographic locations that spanned 55 river kilometers. Forty-six fish were recovered at Eureka Bar (rkm 307.1) and Kirby Creek (rkm 352.0), which accounted for 37% of all collected fish in 2019. Most (73 fish) post-spawned salmon were collected from early to mid-November just after peak spawning. A summary of 2019 PBT results can be found in Appendix A.1. In 2019, we PIT tagged subyearling fall Chinook salmon in both the Snake and Clearwater rivers. In the Snake River, we tagged 410 fish with 8-mm tags, 666 fish with 9-mm tags, and 1,082 fish 12-mm tags. During seining, our recapture rate of previously tagged fish was slightly higher in the lower reach at 11.4% than in the upper reach at 10.6%. In an effort to represent more of the population through tagging, we tagged fish as small as 45 mm with 8-mm tags at one site in the upper reach. This allowed us to increase the number of fish tagged in that reach by 17.7%. An additional 10.6% of collected fish could have been tagged in the lower reach had we used 8-mm tags in that reach. In the Clearwater River, we tagged 2,451 subyearlings and recaptured 260 (10.6%) fish in the river and 66 fish (48 tagged by USGS, 18 tagged by the Nez Perce Tribe) at Lower Granite Dam during October to provide information for growth estimation. Within riverine habitats, growth in both length and mass were higher for fish tagged with 8-mm tags than with 9- and 12-mm tags. Estimated growth in length and mass of subyearlings was higher in Lower Granite Reservoir than in riverine habitats. We adapted existing statistical models used to estimate abundance of steelhead and spring/summer Chinook salmon for fall Chinook salmon passing Lower Granite Dam. Run reconstruction efforts to date at Lower Granite Dam for Snake River fall Chinook salmon have provided estimates of the number of returning adults but with no measure of uncertainty about the estimates. The objective of this study was to estimate the abundance, with uncertainty, of marked (coded-wire tagged CWT or adipose clipped) and unmarked fall Chinook salmon past Lower Granite Dam for return years 20032018. Estimating uncertainty is important for informing the state-space life cycle model (Chapter 5), which incorporates both observation and process uncertainty into parameter estimates. The coefficient of variation (CV) for log- abundance was 1.0% or less in all years, whereas the CV for abundance averaged 4.2% and ranged from 1.4% to 10.4% among years. Over the past five years, we have been developing a two-stage state-space life-cycle model for naturally produced fall Chinook salmon in the Snake River basin. Initial efforts focused on generating juvenile and adult abundance estimates, with estimates of uncertainty, for informing the life-cycle model. In this report we 1) describe the statistical life-cycle model and improvements made since our last report to the Independent Scientific Advisory Board (ISAB), 2) estimate the effects of covariates on key demographic parameters, and 3) use the fitted life- cycle model to simulate population trajectories under hydrosystem actions proposed for the NOAA 2020 Biological Opinion (hereafter, the Proposed Action). Major recent advancements to the model include revised juvenile abundance estimates, the ability to estimate smolt-to-adult return rates (SAR) separately for subyearling and yearling juvenile fall Chinook salmon, and improvements in the observation model for estimating age, sex, and outmigration structure in adult returns. We examined the effect of numerous environmental, hydrosystem, and ocean covariates on key demographic parameters but only a few covariates were significant. For the adult-to-juvenile transition, we found the maximum weekly river flows during the winter egg- incubation period had a significant negative effect on the juvenile outmigrant abundance from that brood year. For subyearling outmigrants, percent spill during the summer had a significant positive effect on SAR and the mean winter PDO (Pacific Decadal Oscillation) had a significant negative effect on SAR. For yearling outmigrants, NPGO (North Pacific Gyre Oscillation) had a significant positive effect on SAR. We used the fitted model to simulate population trajectories under the Proposed Action, and the median 10-year geometric mean abundance was 8,222 female spawners (interquartile range: 2,592 26,714). Overall, the probability of quasi- extinction (probability of falling below 50 female spawners for 4 consecutive years) was low, with only 1.6% of all simulations having a quasi-extinction probability >0.95. Although quasi-extinction probability was low, we did not assess the additional effect of climate change, which would be expected to increase quasi-extinction probability.

Idaho, Oregon, Washington↗

USGS2017

Water quality and quantity and simulated surface-water and groundwater flow in the Laurel Hill Creek Basin, southwestern Pennsylvania, 1991–2007

Laurel Hill Creek is considered one of the most pristine waterways in southwestern Pennsylvania and has high recreational value as a high-quality cold-water fishery; however, the upper parts of the basin have documented water-quality impairments. Groundwater and surface water are withdrawn for public water supply and the basin has been identified as a Critical Water Planning Area (CWPA) under the State Water Plan. The U.S. Geological Survey, in cooperation with the Somerset County Conservation District, collected data and developed modeling tools to support the assessment of water-quality and water-quantity issues for a basin designated as a CWPA. Streams, springs, and groundwater wells were sampled for water quality in 2007. Streamflows were measured concurrent with water-quality sampling at main-stem sites on Laurel Hill Creek and tributaries in 2007. Stream temperatures were monitored continuously at five main-stem sites from 2007 to 2010. Water usage in the basin was summarized for 2003 and 2009 and a Water-Analysis Screening Tool (WAST) developed for the Pennsylvania State Water Plan was implemented to determine whether the water use in the basin exceeded the “safe yield” or “ the amount of water that can be withdrawn from a water resource over a period of time without impairing the long-term utility of a water resource .” A groundwater and surface-water flow (GSFLOW) model was developed for Laurel Hill Creek and calibrated to the measured daily streamflow from 1991 to 2007 for the streamflow-gaging station near the outlet of the basin at Ursina, Pa. The CWPA designation requires an assessment of current and future water use. The calibrated GSFLOW model can be used to assess the hydrologic effects of future changes in water use and land use in the basin. Analyses of samples collected for surface-water quality during base-flow conditions indicate that the highest nutrient concentrations in the main stem of Laurel Hill Creek were at sites in the northeastern part of the basin where agricultural activity is prominent. All of the total nitrogen (N) and a majority of the total phosphorus (P) concentrations in the main stem exceeded regional nutrient criteria levels of 0.31 and 0.01 milligrams per liter (mg/L), respectively. The highest total N and total P concentrations in the main stem were 1.42 and 0.06 mg/L, respectively. Tributary sites with the highest nutrient concentrations are in subbasins where treated wastewater is discharged, such as Kooser Run and Lost Creek. The highest total N and total P concentrations in subbasins were 3.45 and 0.11 mg/L, respectively. Dissolved chloride and sodium concentrations were highest in the upper part of the basin downstream from Interstate 76 because of road deicing salts. The mean base-flow concentrations of dissolved chloride and sodium were 117 and 77 mg/L, respectively, in samples from the main stem just below Interstate 76, and the mean concentrations in Clear Run were 210 and 118 mg/L, compared to concentrations less than 15 mg/L in tributaries that were not affected by highway runoff. Water quality in forested tributary subbasins underlain by the Allegheny and Pottsville Formations was influenced by acidic precipitation and, to a lesser extent, the underlying geology as indicated by pH values less than 5.0 and corresponding specific conductance ranging from 26 to 288 microsiemens per centimeter at 25 degrees Celsius for some samples; in contrast, pH values for main stem sites ranged from 6.6 to 8.5. Manganese (Mn) was the only dissolved constituent in the surface-water samples that exceeded the secondary maximum contaminant level (SMCL). More than one-half the samples from the main stem had Mn concentrations exceeding the SMCL level of 50 micrograms per liter (μg/L), whereas only 19 percent of samples from tributaries exceeded the SMCL for Mn. Stream temperatures along the main stem of Laurel Hill Creek became higher moving downstream. During the summer months of June through August, the daily mean temperatures at the five sites exceeded the limit of 18.9 degrees Celsius (°C) for a cold-water fishery. The maximum instantaneous values for each site ranged from 27.2 to 32.8 °C. Water-quality samples collected at groundwater sites (wells and springs) indicate that wells developed within the Mauch Chunk Formation had the best water quality, whereas wells developed within the Allegheny and Pottsville Formations yielded the poorest water quality. Waters from the Mauch Chunk Formation had the highest median pH (7.6) and alkalinity (80 mg/L calcium carbonate) values. The lowest pH and alkalinity median values were in waters from the Allegheny and Pottsville Formations. Groundwater samples collected from wells in the Allegheny and Pottsville Formations also had the highest concentrations of dissolved iron (Fe) and dissolved Mn. Seventy-eight percent of the groundwater samples collected from the Allegheny Formation exceeded the SMCL of 300 μg/L for Fe and 50 μg/L for Mn. Forty-three and 62 percent of the groundwater samples collected from the Pottsville Formation exceeded the SMCL for iron and Mn, respectively. The highest Fe and Mn concentrations for surface waters were measured for tributaries draining the Pottsville Formation. The highest median Fe concentration for tributaries was in samples from streams draining the Allegheny Formation. During base-flow conditions, the streamflow per unit area along the main stem of Laurel Hill Creek was lowest in the upper parts of the basin [farthest upstream site 0.07 cubic foot per second per square mile (ft 3 /s/mi 2 )] and highest (two sites averaging about 0.20 (ft 3 /s/mi 2 ) immediately downstream from Laurel Hill Lake in the center of the basin. Tributaries with the highest streamflow per unit area were those subbasins that drain the western ridge of the Laurel Hill Creek Basin. The mean streamflow per unit area for tributaries draining areas that extend into the western ridge and draining eastern or central sections was 0.24 and 0.05 ft 3 /s/mi 2 , respectively. In general, as the drainage area increased for tributary basins, the streamflow per unit area increased. Criteria established by the Pennsylvania Department of Environmental Protection indicate that the safe yield of water withdrawals from the Laurel Hill Creek Basin is 1.43 million gallons per day (Mgal/d). Water-use data for 2009 indicate that net (water withdrawals subtracted by water discharges) water withdrawals from groundwater and surface-water sources in the basin were approximately 1.93 Mgal/d. Water withdrawals were concentrated in the upper part of the basin with approximately 80 percent of the withdrawals occurring in the upper 36 mi 2 of the basin. Three subbasins—Allen Creek, Kooser Run, and Shafer Run— in the upper part were affected the most by water withdrawals such that safe yields were exceeded by more than 1,000 percent in the first two and more than 500 percent in the other. In the subbasin of Shafer Run, intermittent streamflow characterizes sections that historically have been perennial. The GSFLOW model of the Laurel Hill Creek Basin is a simple one-layer representation of the groundwater flow system. The GSFLOW model was primarily calibrated to reduce the error term associated with base-flow periods. The total amount of observed streamflow at the Laurel Hill Creek at Ursina, Pa. streamflow-gaging station and the simulated streamflow were within 0.1 percent over the entire modeled period; however, annual differences between simulated and observed streamflow showed a range of -27 to 24 percent from 1992 to 2007 with nine of the years having less than a 10-percent difference. The primary source of simulated streamflow in the GSFLOW model was the subsurface (interflow; 62 percent), followed by groundwater (25 percent) and surface runoff (13 percent). Most of the simulated subsurface flow that reached the stream was in the form of slow flow as opposed to preferential (fast) interflow.

Pennsylvania↗

USGS2014

Dissolved-solids sources, loads, yields, and concentrations in streams of the conterminous United States

Recent studies have shown that excessive dissolved-solids concentrations in water can have adverse effects on the environment and on agricultural, domestic, municipal, and industrial water users. Such effects motivated the U.S. Geological Survey’s National Water Quality Assessment Program to develop a SPAtially-Referenced Regression on Watershed Attributes (SPARROW) model that has improved the understanding of sources, loads, yields, and concentrations of dissolved solids in streams of the conterminous United States. Using the SPARROW model, long-term mean annual dissolved-solids loads from 2,560 water-quality monitoring stations were statistically related to several spatial datasets that are surrogates for dissolved-solids sources and land-to-water delivery processes. Specifically, sources in the model included variables representing geologic materials, road deicers, urban lands, cultivated lands, and pasture lands. Transport of dissolved solids from these sources was modulated by land-to-water delivery variables that represent precipitation, streamflow, soil, vegetation, terrain, population, irrigation, and artificial drainage characteristics. Where appropriate, the load estimates, source variables, and transport variables were statistically adjusted to represent conditions for the base year 2000. The nonlinear least-squares estimated SPARROW model was used to predict long-term mean annual conditions for dissolved-solids sources, loads, yields, and concentrations in a digital hydrologic network representing nearly 66,000 stream reaches and their corresponding incremental catchments that drain the Nation. Nationwide, the predominant source of dissolved solids yielded from incremental catchments and delivered to local streams is geologic materials in 89 percent of the catchments, road deicers in 5 percent of the catchments, pasture lands in 3 percent of the catchments, urban lands in 2 percent of the catchments, and cultivated lands in 1 percent of the catchments. Whereas incremental catchments with dissolved solids that originated predominantly from geologic sources or from urban lands are found across much of the Nation, incremental catchments with dissolved solids yields that originated predominantly from road deicers are largely found in the Northeast, and incremental catchments with dissolved solids that originated predominantly from cultivated or pasture lands are largely found in the West. The total amount of dissolved solids delivered to the Nation’s streams is 271.9 million metric tons (Mt) annually, of which 194.2 million Mt (71.4%) come from geologic sources, 37.7 million Mt (13.9%) come from road deicers, 18.2 million Mt (6.7%) come from pasture lands, 13.9 million Mt (5.1%) come from urban lands, and 7.9 million Mt (2.9%) come from cultivated lands. Nationwide, the median incremental-catchment yield delivered to local streams is 26 metric tons per year per square kilometer [(Mt/yr)/km 2 ]. Ten percent of the incremental catchments yield less than 4 (Mt/yr)/km 2 , and 10 percent yield more than 90 (Mt/yr)/km 2 . Incremental-catchment yields greater than 50 (Mt/yr)/km 2 mostly occur along the northern part of the West Coast and in a crescent shaped band south of the Great Lakes. For example, the median incremental-catchment yield is 81 (Mt/yr)/km 2 for the Great Lakes, 78 (Mt/yr)/km 2 for the Ohio, and 74 (Mt/yr)/km 2 for the Upper Mississippi water-resources regions. Incremental-catchment yields less than 10 (Mt/yr)/km 2 mostly occur in a wide band across the arid lowland of the interior West that excludes areas along the coast and the extensive, higher mountain ranges. For example, the median incremental-catchment yield is 3 (Mt/yr)/km 2 for the Lower Colorado, 5 (Mt/yr)/km 2 for the Rio Grande, and 8 (Mt/yr)/km 2 for the Great Basin water-resources regions. Predicted incremental loads were cascaded down through the reach network, with loads accumulating from reach to reach. For most stream reaches, the entire incremental load of dissolved solids delivered to the reach was transported to either the ocean or to one of the large streams flowing along the U.S. international boundary without losses occurring along the way. The exceptions to this include streams in the southwestern part of the country, such as the Colorado River, Rio Grande, and streams of internally drained drainages in the Great Basin, where dissolved-solids loads decreased through streamflow diversion for off-stream use, or by infiltration through the streambed. Long-term mean annual flow-weighted concentrations were derived from the predicted accumulated-load and stream-discharge data. Widespread low concentrations, generally less than 100 milligrams per liter (mg/L), occur in many reaches of the New England, South Atlantic-Gulf, and Pacific Northwest water-resources regions as a result of moderate dissolved-solids yields and high runoff rates. Widespread moderate concentrations, generally between 100 and 500 mg/L, occur in many reaches of the Great Lakes, Ohio, and Upper Mississippi River water-resources regions. Whereas dissolved-solids yields are generally high in these regions, runoff rates are also high, which helps moderate concentrations in these regions. Widespread higher concentrations, generally greater than 500 mg/L, occur across a belt of reaches that extends almost continuously from Canada to Mexico in the Midwest, cutting through the Souris-Red-Rainy, Missouri, Arkansas-White-Red, Texas-Gulf, and Rio Grande water-resources regions. Although dissolved-solids yields are moderate to low in these areas, low runoff rates result in the high concentrations for these areas. In 12.6 percent of the Nation’s stream reaches, predicted concentrations of dissolved solids exceed 500 mg/L, the U.S. Environmental Protection Agency’s secondary, nonenforceable drinking water standard. While this standard provides a metric for evaluating predicted concentrations in the context of drinking-water supplies, it should be noted that it only applies to drinking water actually served to customers by water utilities, and it does not apply to all stream reaches in the Nation nor does it apply during times when water is not being withdrawn for use. Exceedance of 500 mg/L is more pronounced in certain water-resources regions than others. For example, about half of the reaches in the Souris-Red-Rainy region have concentrations predicted to exceed 500 mg/L, and between 25 and 37 percent of the reaches in the Missouri, Arkansas-White-Red, Texas-Gulf, Rio Grande, and Lower Colorado regions are predicted to exceed 500 mg/L. Development of stream-load data for use in the SPARROW model also provided long-term temporal trend information in dissolved-solids concentrations at the monitoring stations for their period of record, which was constrained between 1980 and 2009. For the 2,560 monitoring stations used in this study, long-term trends in flow-adjusted dissolved-solids concentrations increased over time at 23 percent of the stations, decreased at 18 percent of the stations, and did not change over time at 59 percent of the stations. Long-term trends show a strong regional spatial pattern where from the western parts of the Great Plains to the West Coast, concentrations mostly either did not change or decreased over time, and from the eastern parts of the Great Plains to the East Coast, concentrations mostly either did not change or increased over time. Results from the trend analysis and from the SPARROW model indicate that, compared to monitoring stations with no trends or decreasing trends, stations with increasing trends are associated with a smaller percentage of the predicted dissolved-solids load originating from geologic sources, and a larger percentage originating from urban lands and road deicers. Conversely, compared to stations with increasing trends or no trends, stations with decreasing trends have a larger percentage of the predicted dissolved-solids load originating from geologic sources and a smaller percentage originating from urban lands and road deicers. Stations with decreasing trends also have larger percentages of predicted dissolved-solids load originating from cultivated lands and pasture lands, compared to stations with increasing trends or no trends.

Scientific Investigations Report↗

USGS2016

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

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

Kansas↗

USGS2021

Capacity assessment for Earth Monitoring, Analysis, and Prediction (EarthMAP) and future integrated monitoring and predictive science at the U.S. Geological Survey

Executive Summary Managers of our Nation’s resources face unprecedented challenges driven by the convergence of increasing, competing societal demands and a changing climate that affects the stability, vulnerability, and predictability of those resources. To help meet these challenges, the scientific community must take advantage of all available technologies, data, and integrative Earth systems modeling capacity to better inform resource and risk management decisions. This is the overarching goal of the U.S. Geological Survey (USGS) Earth Monitoring, Analysis, and Prediction (EarthMAP) vision: “By 2030, the USGS will deliver well integrated observations and predictions of the future state of natural systems—water, ecosystems, energy, minerals, hazards—at regional and national scales, working primarily with federal, state, and academic partners to develop and operate the capability” (U.S. Geological Survey, 2021). Providing more integrated Earth systems science and actionable information to decision makers, stakeholders, and the public requires a better understanding of the depth and distribution of existing capacity (capabilities, tools, and techniques) across the Bureau. Identifying existing capacity is also a critical first step toward gap analysis and targeted investments to increase capacity over time. The USGS formed a Capacity Assessment Team (CAT) and charged it with (1) conducting a Request for Information (RFI) to identify existing USGS expertise and activities supportive of integrated and predictive science to inform decision making, (2) developing a strategy and proof-of-concept for a continuously updated capacity assessment capability, and (3) identifying lessons learned to inform development of best practices for future capacity assessment efforts. The RFI took the form of a survey, with content guided by the science and technology needs identified in a USGS report titled “Grand Challenges for Integrated U.S. Geological Survey Science—A Workshop Report” (Jenni and others, 2017). The 44-question survey provided respondents the ability to rate their level of experience with a suite of priority disciplines, analysis and modeling approaches, technologies, and stakeholder engagement strategies and to enter optional narrative text for supporting context. An introductory portion focused on general science capacity assessment, followed by three sections targeting capabilities related to the foundational components of EarthMAP: (1) data and information integration, (2) integrated predictive science, and (3) actionable information. The survey results provided a high-level snapshot of USGS capacity in the targeted areas. Respondents (1,035 individuals) represented approximately 13 percent of the USGS across all mission areas and regions. Seventy-four percent of the respondents held a science-focused position title and the remainder had position titles in information technology, computer science, management, administrative, or other (contractors, volunteers, emeritus, and unknown). To provide greater insight into respondent capabilities and activities, information from the U.S. Department of the Interior and USGS enterprise information systems were used to further characterize topical expertise and organizational associations of survey respondents. To address the ongoing need to assess the Bureau’s capacity to address integrated predictive science priorities, the CAT developed a software-based proof-of-concept called the Integrated Science Assessment Information Database (iSAID) for assembling various information sources together toward making the full extent of USGS capabilities and scientific assets available for routine capacity assessment. This proof-of-concept is intended to serve as a catalyst for further development. The process of implementing the EarthMAP capacity assessment survey, analyzing survey responses, and developing the proof-of-concept resulted in lessons learned, findings, and recommendations. Example scenarios throughout the report demonstrate how capacity assessment data can inform science planning. Three overarching findings and recommendations are: (1) Finding: Capacity is limited in some critical disciplines, skills, and technology applications, but “sufficient” depends on the question and the need relative to availability at a given point in time. Recommendation: Develop an on-demand capacity assessment framework that enables rapid identification and evaluation of existing and available expertise to support decision needs as they arise. (2) Finding: Institutional barriers and lack of awareness constrain the ability of USGS staff to adopt new technologies, collaborate across administrative boundaries, and deliver actionable information to stakeholders in a timely manner. However, these barriers are not universally experienced. Recommendation: Pursue more targeted inquiries to clarify which institutional barriers are obstructing the adoption of new technologies and approaches or the sharing of expertise and equipment across organizational and regional boundaries. These inquiries should inform USGS leadership, mission areas, and regions whether policies can be revised or whether a lack of understanding is creating perceived obstacles. Highlight cases when staff have successfully adopted new technologies and approaches to advance EarthMAP priorities and provide actionable information in a timely manner to spread awareness of how perceived obstacles can be navigated and overcome when appropriate. (3) Finding: Examples of people and projects integrating across disciplines and scales and applying advanced approaches to meet complex stakeholder needs exist. Such examples provide transfer value across the spectrum from approach to decision making. Many projects, already underway, appear to meet elements of the EarthMAP vision, and the USGS has people who can provide leadership in multiple types of specific integrated science efforts. Recommendation: Use these findings as a starting point for near-term strategic planning for integrated science. Highlight, incentivize, and build on existing interdisciplinary predictive science and information delivery activities across the USGS to advance toward further realization of an EarthMAP capacity. The CAT efforts to develop and assess existing USGS capacity to advance the EarthMAP vision revealed a fundamental challenge for not only this effort but any effort to assess existing capacity: A considerable amount of thought, time, and effort is required to survey and assess capabilities and tools available to support a given need, yet best results are still likely to provide an incomplete assessment. To better meet the frequent need to assess capabilities, tools, products, and projects that address an expressed strategic priority, the CAT proposes the concept of an on-demand capacity assessment framework supported by a software package that dynamically pulls and integrates information from existing USGS information systems and public domain registries. Although existing USGS enterprise information systems currently lack the structure, cross-system consistency, interoperability, and stability to support a continuously updated capacity assessment capability, we identify reasonable near-term steps to improve the utility of information gathered on expertise and project capacity and to improve the consistency and completeness of information and the ability of USGS systems to share that information. The ability to search and characterize this information will make future assessments of capacity faster, more complete, more efficient, and more targeted. This approach would grow the Bureau’s capacity knowledge over time, iteratively improving the ability to access, leverage, and synthesize existing capabilities and assets as well as identify and fill critical gaps. The greatest promise for developing integrated science could lie in linking across existing projects and expertise to create a multi-project capacity for addressing large, complex environmental issues.

Open-File Report↗

USGS2014

Elk monitoring in Lewis and Clark National Historical Park: 2008-2012 synthesis report

Maintaining elk (Cervus elaphus roosevelti) herds that frequent Lewis and Clark National Historical Park (NHP) is central to the park’s purpose of preserving the historic, cultural, scenic, and natural resources associated with the winter encampment of the Lewis and Clark expedition. Elk were critically important to the Lewis and Clark expedition in providing food and hides that sustained the expedition during the winter of 1805-06 and supplied them for their return east during 1806. Today, elk remain a key component of interpreting the Lewis and Clark story to over 200,000 park visitors each year at the Fort Clatsop visitor center. In 2008, the US Geological Survey (USGS) began collaborating with Lewis and Clark NHP and the NPS North Coast and Cascades Network to develop a protocol for monitoring long-term changes in the magnitude and spatial patterns of elk use within and adjacent to Lewis and Clark NHP (Griffin et al. 2011). Specific objectives of the monitoring program were to measure trends in (1) relative use of the Fort Clatsop unit by elk during winter; (2) the proportion of areas where elk sign is present in the Fort Clatsop unit in winter; and (3) the frequency of elk sightings from roads in and around the Fort Clatsop unit. This report synthesizes the results of the first four years of monitoring elk distribution and use in Lewis and Clark NHP from 2008-2012. We also present data from FY2012 (Appendix 1), in lieu of an annual report for that year. We used fecal pellet group surveys as the cornerstone for monitoring trends in both relative use of the Fort Clatsop Unit by elk and the proportion of areas where elk sign was present at the end of winter. We estimated pellet group density based on data collected from a network of fecal pellet plots distributed systematically throughout the unit. We developed a double observer sampling scheme that enabled us to estimate detection biases and improve the accuracy of pellet group density estimates. We computed the estimated detection probability for any pellet group observed; this probability was a function of the pellet group size and stage of decay, as well as lighting and vegetation conditions, and the number of observers (one or two) searching for pellets in that subplot. We then used these estimated detection probabilities to adjust the raw counts of the detected pellet groups to account for groups that likely went undetected under similar pellet and environmental conditions (each observed pellet group was weighted by the inverse of its estimated detection probability). We also used results from the late winter fecal pellet surveys to quantify the proportion of areas where elk pellets occurred (PAO), which was based on the presence of fecal pellet groups and estimation of detection biases (i.e., accounting for pellet groups that likely went undetected by both observers). In this synthesis, we report temporal trends in both pellet group density and PAO from 2008-2012, based on weighted linear regression analyses as well as spatial variation of pellet group densities over time. We completed late winter fecal pellet surveys at 61-66 plots annually, depending on yearly variation in access. We cleared fecal pellets at survey points in late October / early November each year and returned in late February / early March to count pellet groups left by elk over the winter. The estimated probability that a pellet group was detected by any one observer during late winter was affected most by the pellet group size and was less affected by decay class and lighting conditions. Per-observer detection probabilities ranged from as low as ~10-15% for single pellets to ~85-90% for pellet groups with 50 pellets. Average pellet group density in the Fort Clatsop unit ranged annually from 0.58 (+/- 1.43 standard error [SE]) to 0.93 (+/- 2.25 SE) pellet groups per 3-m radius subplot. Pellet group density declined over time, at approximately 8.8% per year (+/- 2.5% SE), but that slope was not statistically distinguishable from zero (2- tailed P=0.16). Following correction for detection biases, the proportion of surveyed points used by elk (i.e., PAO) ranged from 0.44 (+/- 0.07 SE) to 0.53 (+/- 0.07 SE) during the 4 winters. The estimated proportion of areas where elk pellets occurred (PAO) declined at a rate of 2.6% per year (+/- 1.2% per year SE), but that trend also was not statistically distinguishable from zero (2- tailed P=0.17). Statistical significance of a measure’s trend depends on both the magnitude (i.e., slope) of the observed trend and the number of years the trend continues in the same increasing or decreasing direction. Through simulation modeling we determined how many additional years of surveys would be required to reveal a statistically significant trend, based on the same trends in pellet group density and PAO, and associated variation, observed from 2009-2012. Assuming the same trends persist in the future, simulations indicated that there is a 70% probability that a statistically significant trend would be detected after two more years of conducting pellet group surveys. Relative use by elk during winter, as indexed by elk pellet group density, was generally greatest in the southeast region of the Fort Clatsop unit in or near the large freshwater marsh at the mouth of Colewort Creek and adjacent upland areas. Pellet group density was also higher than average in the north-central forested area, not far from a privately-owned pasture north of the park boundary. This spatial pattern in pellet group densities across the Fort Clatsop unit was consistent across all four years, although specific pellet group densities varied from year to year. Pellet group density declined significantly over time at two points in the southeast of the Fort Clatsop unit, even though pellet group density at those points remained higher than the unit average. Pellet group density increased significantly over time at one point in the north-central region, and at one point in the south-central region of the unit, indicating a slight shift in the distribution of elk use within the Fort Clatsop Unit over the four years. As an index of visitors’ opportunities to see elk in and around the Fort Clatsop Unit, we conducted replicated roadside elk surveys 3-5 times monthly during February, April, June, August, October and December 2008-2012. During each morning of survey, we searched for elk along four routes that totaled 32 km. We examined bimonthly trends in the numbers of elk groups seen, the total number of elk seen, and the observed composition ratios for those six months of the year. The average number of elk groups seen per survey ranged from 0.75 (+/- 0.32 SE) during February to a peak of 1.95 (+/- 0.36 SE) during June. Despite this seasonal variation in numbers of elk groups seen, the average total number of elk seen per morning was less variable. The average ratios of antlered elk to antlerless adult elk (i.e., bulls:cows) and calves to antlerless adult elk (i.e. calves:cows) varied seasonally, with the highest of both average ratios observed in August. We detected no significant trends in the average number of elk groups and total numbers of elk seen per survey from 2008-2012. Similarly, ratios of calves and antlered elk per antlerless elk did not differ over time. Elk groups were frequently seen from January to August in the southeast region of the Fort Clatsop unit, in the vicinity of Colewort Creek. Outside of NPS lands, we observed elk most frequently in open areas near the Astoria regional airport, in the pastures and forests immediately north of the Fort Clatsop unit and, prior to the construction of a residential development, in a pasture northwest of the Fort Clatsop unit. Elk monitoring at Lewis and Clark NHP is still in its initial years and additional monitoring will be required to verify trends that appear to be emerging. For example, the initial monitoring suggested incipient declining trends in both pellet group density and proportion of plots with pellets present, as well as, potentially, a small shift in elk distribution away from a new trail that was recently constructed in the southeast portion of the Fort Clatsop unit. Continued monitoring will aid in determining whether this local change in distribution persists (or, alternatively, resulted from short-term random variation), and whether there will be any positive or negative effect in the northern portion of the unit where a new trail has been constructed. High variability in road counts prevented our ability to find any clear trend in numbers or composition of elk observed in and near Fort Clatsop, but changes in the patterns of observations of elk from roadways suggest that residential development outside the park has reduced the available habitat for elk in some of the areas surveyed, and may have affected spatial use patterns of elk adjacent to some areas of the park. In addition to monitoring future effects of land use changes outside the park, continued monitoring may also prove useful for assessing elk responses to natural succession in forests disturbed by windthrow in December 2007 and to NPS vegetation management activities such as variable density thinning in the forest, trail development, and restoration at Otter Point tidal area and Colewort Creek Slough.

Oregon;Washington↗

USGS2023

Evaluating management alternatives for Wyoming elk feedgrounds in consideration of chronic wasting disease

Executive Summary The authors used decision and modeling analyses to evaluate management alternatives for a decision on whether to permit Cervus canadensis (elk) feeding on two sites on Bridger-Teton National Forest, Dell Creek and Forest Park. Supplemental feeding of elk could increase the transmission of chronic wasting disease (CWD) locally and disease spread regionally, potentially impacting elk populations over time with wider implications for Odocoileus hemionus (mule deer) and Odocoileus virginianus (white-tailed deer) populations and hunting, tourism, and regional revenue. Supplemental feeding is thought to improve overwinter elk survival and reduce the commingling of elk with cattle during months when brucellosis transmission risk is highest. We worked with the U.S. Department of Agriculture Forest Service to identify their fundamental objectives and associated performance metrics related to this feedground decision. We then developed disease and habitat selection models to quantify the effect of four management alternatives on select performance metrics. The four alternatives were to continue to permit feeding, phaseout permits to feed in three years, permit feeding on an emergency basis, or stop permitting feeding. In this report, we present methods and summarized results on disease and habitat selection models and summaries of other performance metrics analyzed by BIO-WEST, Inc. and Cirrus Ecological Solutions as part of an Environmental Impact Statement. Data from Wyoming Game and Fish Department (WGFD) supported the assumption that supplemental elk feeding allows for larger elk populations in a region. We documented that herd units (HU) without feedgrounds had 23 percent lower densities of elk per area of winter range when compared against HUs with feedgrounds, after accounting for differences in sightability of elk during counts on and off feedgrounds. Thus, throughout our analyses, we assumed feedground closures would reduce elk carrying capacity resulting in an average decline of previously fed elk population segments by 23 percent (5th and 95th percentiles = [11 percent, 35 percent]) by year 20. Most of that decline occurred within the first few years after a feedground ceases to operate. We used a panel of CWD experts to help estimate CWD trans-mission in fed and unfed elk population segments. In aggregate, the expert panel estimated that median values of direct and indirect transmission of CWD are expected to be 1.9 and 4 times higher, respectively, in fed elk populations compared to unfed elk. We used these disease transmission estimates in combination with local elk demographic rates and carrying capacity estimates to project disease and population dynamics. In year 20, we predicted CWD prevalence would increase to 42 percent (5th and 95th percentiles = [29 percent, 55 percent]), and 13 percent (5th and 95th percentiles = [4 percent, 26 percent]) on average for fed and unfed elk population segments, respectively, given a starting prevalence of 1.6 percent. The prevalence estimates for the unfed elk population segments are in the range of previous observations of CWD in elk in the western United States. The average CWD prevalence from 2016 to 2018 in the unfed elk population of Wind Cave National Park in South Dakota was 18 percent overall but up to 30 percent in some regions (Sargeant and others, 2021). Meanwhile, CWD prevalence in the Iron Mountain and Laramie Peak elk herds in Wyoming from 2016 to 2018 was 14 percent and 7 percent, respectively, despite being present since at least 2002 (Wyoming Game and Fish Department, 2020b). From 2016 to 2020, elk that were fed at Dell Creek and Forest Park constituted on average 12–20 percent of the total elk on their respective HUs. As a result, the differences between management alternatives are modest when considering the closure of only one feedground on a HU. The no feeding alternative for Forest Park resulted in a CWD prevalence of 17 percent (SD = 7 percent) in the Afton HU compared to 20 percent (SD = 7 percent) with continued feeding by year 20. In the Upper Green River HU, no feeding on Dell Creek resulted in a CWD prevalence of 27 percent (SD = 6 percent) compared to 30 percent (SD = 5 percent) with continued feeding. In terms of disease-associated mortality, we predicted the closure of Forest Park and Dell Creek feedgrounds would reduce the total number of CWD mortalities by 9 percent in the Upper Green River HU and 26 percent in the Afton HU during the 20-year timespan. Our spatial analyses predicted that management alternative effects vary by HU as a function of private property and other wildlife winter ranges proximity relative to feedground location. The predicted number of elk abortions on private land, as a proxy for brucellosis risk to cattle, may increase by 8–21 percent in the absence of feeding at Dell Creek and Forest Park. Eight feedgrounds are located on Bridger-Teton National Forest, all of which have permits that have expired or will expire prior to 2028. In addition, WGFD could change their management of feedgrounds given new information; therefore, we also assessed the cumulative effects of continued feeding, phaseout, and no feeding management alternatives across five HUs south of Jackson, Wyoming (Afton HU, Fall Creek HU, Piney HU, Pinedale HU, and Upper Green River HU). These five HUs ranged from about 41 to 85 percent of the elk herd using feedgrounds, which corresponded to a CWD prevalence at year 20 of 23–34 percent if all feedgrounds in those five HUs remained open relative to 12 to 14 percent if all feedgrounds were closed. We predicted feedground closures may result in immediate reductions in population size relative to alternatives that continue feeding (for example, continued feeding and emergency feeding alternatives); however, over longer periods of time, CWD-associated mortality leads to larger population reductions. The no feeding alternative resulted in higher elk population sizes compared to the continued feeding alternative after about 10 years of implementation. Delayed action under a phaseout alternative resulted in increasing the CWD prevalence to 20 percent relative to 12 to 14 percent, on average, without feeding on HUs with a large population of fed elk such as the Upper Green River HU. Summarizing our cumulative results across all five of the analyzed HUs, we predicted continued feeding will lead to fewer elk by year 20 (mean = 8,300, standard deviation [SD] = 740) compared to no feeding at U.S. Department of Agri-culture Forest Service sites (10,700, SD = 890). The closure of all feedgrounds was projected to result in the largest elk populations at year 20 (12,500, SD = 980). No feeding at all sites also resulted in the largest cumulative harvest of 57,700 (SD = 2,600) compared to 51,100 (SD = 3,800) for continued feeding at all current feedground sites on the five HUs. Continued feeding also resulted in the lowest brucellosis costs to producers ($194,600, SD = $11,500) compared to no feeding on all feedgrounds ($243,000, SD = $13,700). Assuming moderate reductions in hunter interest because of increasing CWD prevalence in elk, we predicted that no feeding resulted in regional revenues generated by hunting activities of $190 million (SD = $10 million) compared to $173 million (SD = $10 million) for continued feeding over the 20-year timeframe. Recent CWD detections in mule deer and elk in Grand Teton National Park has elevated the importance of the cur-rent decision on whether, and how, to permit elk feeding on Dell Creek and Forest Park and the management of the other feedgrounds. Aggressive male harvest has slowed, but not stopped, the increasing prevalence of CWD in mule deer (Conner and others, 2021). It is unclear whether harvest management can be an effective tool to slow the spread of CWD in elk. There are also no effective treatments or vaccines for CWD, and it is unlikely that any will be developed that can be easily deployed in the near future. Thus, reducing artificial aggregations is one of the few management approaches suggested by the Western Association of Fish and Wildlife Agencies (Almberg and others, 2017). Future surveillance and monitoring can be designed to resolve uncertainties that can improve future decision-making. If feedgrounds close, research could quantify elk population reductions in the absence of feeding, the redistribution of fed elk to other places, or the consequences of elk movement on private property. If feedgrounds remain open, research could assess how rapidly CWD spreads in artificial aggregations of elk; however, surveillance programs would need to be designed with sufficient power to detect initial changes of CWD prevalence. Delaying action on feedground management was projected to be costly. Results of the phaseout alternative relative to the no feeding alternative suggested a 3-year delay was enough for substantial long-term changes in CWD prevalence. The long-term persistence of infectious CWD prions in the environment suggests that feedground management decisions may have long-lasting consequences. Our results indicated tradeoffs in the ability of a management agency to achieve all their objectives, and all management alternatives resulted in significant reductions in elk population size. This report contains the foundational elements for formal decision analysis methods, which can be implemented to help decision makers transparently evaluate the consequences of decision alternatives and identify the set of actions that best achieve agency and stakeholder priorities.

Wyoming↗

USGS2011

Borehole geophysical investigation of a formerly used defense site, Machiasport, Maine, 2003-2006

The U.S. Geological Survey, in cooperation with the U.S. Army Corps of Engineers, collected borehole geophysical logs in 18 boreholes and interpreted the data along with logs from 19 additional boreholes as part of an ongoing, collaborative investigation at three environmental restoration sites in Machiasport, Maine. These sites, located on hilltops overlooking the seacoast, formerly were used for military defense. At each of the sites, chlorinated solvents, used as part of defense-site operations, have contaminated the fractured-rock aquifer. Borehole geophysical techniques and hydraulic methods were used to characterize bedrock lithology, fractures, and hydraulic properties. In addition, each geophysical method was evaluated for effectiveness for site characterization and for potential application for further aquifer characterization and (or) evaluation of remediation efforts. Results of borehole geophysical logging indicate the subsurface is highly fractured, metavolcanic, intrusive, metasedimentary bedrock. Selected geophysical logs were cross-plotted to assess correlations between rock properties. These plots included combinations of gamma, acoustic reflectivity, electromagnetic induction conductivity, normal resistivity, and single-point resistance. The combined use of acoustic televiewer (ATV) imaging and natural gamma logs proved to be effective for delineating rock types. Each of the rock units in the study area could be mapped in the boreholes, on the basis of the gamma and ATV reflectivity signatures. The gamma and mean ATV reflectivity data were used along with the other geophysical logs for an integrated interpretation, yielding a determination of quartz monzonite, rhyolite, metasedimentary units, or diabase/gabbro rock types. The interpretation of rock types on the basis of the geophysical logs compared well to drilling logs and geologic mapping. These results may be helpful for refining the geologic framework at depth. A stereoplot of all fractures intersecting the boreholes indicates numerous fractures, a high proportion of steeply dipping fractures, and considerable variation in fracture orientation. Low-dip-angle fractures associated with unloading and exfoliation are also present, especially at a depth of less than 100 feet below the top of casing. These sub-horizontal fractures help to connect the steeply dipping fractures, making this a highly connected fracture network. The high variability in the fracture orientations also increases the connectivity of the fracture network. A preliminary comparison of all fracture data from all the boreholes suggests fracturing decreases with depth. Because all the boreholes were not drilled to the same depth, however, there is a clear sampling bias. Hence, the deepest boreholes are analyzed separately for fracture density. For the deepest boreholes in the study, the intensity of fracturing does not decline significantly with depth. It is possible the fractures observed in these boreholes become progressively tighter or closed with depth, but this is difficult to verify with the borehole methods used in this investigation. The fact that there are more sealed fractures at depth (observed in optical televiewer logs in some of the boreholes) may indicate less opening of the sealed fractures, less water moving through the rock, and less weathering of the fracture infilling minerals. Although the fracture orientation remained fairly constant with depth, differences in the fracture patterns for the three restoration sites indicate the orientation of fractures varies across the study area. The fractures in boreholes on Miller Mountain predominantly strike northwest-southeast, and to a lesser degree they strike northeast. The fractures on or near the summit of Howard Mountain strike predominantly east-west and dip north and south, and the fractures near the Transmitter Site strike northeast-southwest and dip northwest and southeast. The fracture populations for the boreholes on or near the summit of Howard Mountain show more variation than at the other two sites. This variation may be related to the proximity of the fault, which is northeast of the summit of Howard Mountain. In a side-by-side comparison of stereoplots from selected boreholes, there was no clear correspondence between fracture orientation and proximity to the fault. There is, however, a difference in the total populations of fractures for the boreholes on or near the summit of Howard Mountain and the boreholes near the Transmitter Site. Further to the southwest and further away from the fault, the fractures at the Transmitter Site predominantly strike northeast-southwest and northwest-southeast.Heat-pulse flowmeter (HPFM) logging was used to identify transmissive fractures and to estimate the hydraulic properties along the boreholes. Ambient downflow was measured in 13 boreholes and ambient upflow was measured in 9 boreholes. In nine other bedrock boreholes, the HPFM did not detect measurable vertical flow. The observed direction of vertical flow in the boreholes generally was consistent with the conceptual flow model of downward movement in recharge locations and upward flow in discharge locations or at breaks in the slope of land surface. Under low-rate pumping or injection rates [0.25 to 1 gallon per minute (gal/min)], one to three inflow zones were identified in each borehole. Two limitations of HPFM methods are (1) the HPFM can only identify zones within 1.5 to 2 orders of magnitude of the most transmissive zone in each borehole, and (2) the HPFM cannot detect flow rates less than 0.010 + or - 0.005 gal/min, which corresponds to a transmissivity of about 1 foot squared per day (ft2/d). Consequently, the HPFM is considered an effective tool for identifying the most transmissive fractures in a borehole, down to its detection level. Transmissivities below that cut-off must be measured with another method, such as packer testing or fluid-replacement logging. Where sufficient water-level and flowmeter data were available, HPFM results were numerically modeled. For each borehole model, the fracture location and measured flow rates were specified, and the head and transmissivity of each fracture zone were adjusted until a model fit was achieved with the interpreted ambient and stressed flow profiles. The transmissivities calculated by this method are similar to the results of an open-hole slug test; with the added information from the flowmeter, however, the head and transmissivity of discrete zones also can be determined. The discrete-interval transmissivities ranged from 0.16 to 330 ft2/d. The flowmeter-derived open-hole transmissivity, which is the combined total of each of the transmissive zones, ranged from 1 to 511 ft2/d. The whole-well open-hole transmissivity values determined with HPFM methods were compared to the results of open-hole hydraulic tests. Despite the fact that the flowmeter-derived transmissivities consistently were lower than the estimates derived from open-hole hydraulic tests alone, the correlation was very strong (with a coefficient of determination, R2, of 0.9866), indicating the HPFM method provides a reasonable estimate of transmissivities for the most transmissive fractures in the borehole. Geologic framework, fracture characterization, and estimates of hydraulic properties were interpreted together to characterize the fracture network. The data and interpretation presented in this report should provide information useful for site investigators as the conceptual site groundwater flow model is refined. Collectively, the results and the conceptual site model are important for evaluating remediation options and planning or implementing the design of a well field and borehole completions that will be adequate for monitoring flow, remediation efforts, groundwater levels, and (or) water quality. Similar kinds of borehole geophysical logging (specifically the borehole imaging, gamma, fluid logs, and HPFM) should be conducted in any newly installed boreholes and integrated with interpretations of any nearby boreholes. If boreholes are installed close to existing or other new boreholes, cross-hole flowmeter surveys may be appropriate and may help characterize the aquifer properties and connections between the boreholes.

Maine↗

USGS2016

Groundwater quality, age, and susceptibility and vulnerability to nitrate contamination with linkages to land use and groundwater flow, Upper Black Squirrel Creek Basin, Colorado, 2013

The Upper Black Squirrel Creek Basin is located about 25 kilometers east of Colorado Springs, Colorado. The primary aquifer is a productive section of unconsolidated deposits that overlies bedrock units of the Denver Basin and is a critical resource for local water needs, including irrigation, domestic, and commercial use. The primary aquifer also serves an important regional role by the export of water to nearby communities in the Colorado Springs area. Changes in land use and development over the last decade, which includes substantial growth of subdivisions in the Upper Black Squirrel Creek Basin, have led to uncertainty regarding the potential effects to water quality throughout the basin. In response, the U.S. Geological Survey, in cooperation with Cherokee Metropolitan District, El Paso County, Meridian Service Metropolitan District, Mountain View Electric Association, Upper Black Squirrel Creek Groundwater Management District, Woodmen Hills Metropolitan District, Colorado State Land Board, and Colorado Water Conservation Board, and the stakeholders represented in the Groundwater Quality Study Committee of El Paso County conducted an assessment of groundwater quality and groundwater age with an emphasis on characterizing nitrate in the groundwater. Groundwater-quality samples were collected from 50 randomly selected wells between May and June 2013. The samples were analyzed for major ions, nutrients, dissolved gases, tritium ( 3 H), chlorofluorocarbons (CFC-11, CFC-12, and CFC-113), and fuel products (such as benzene, toluene, ethylbenzene, and xylenes). None of the groundwater samples exceeded the U.S. Environmental Protection Agency (EPA) National Primary Drinking Water Regulations for primary maximum contaminant levels (MCL) for major ions. Secondary maximum contaminant levels, which are not health concerns and affect mainly taste, color, or odor of the water, were observed in rare instances for pH (2 samples), chloride (1 sample), iron (3 samples), and manganese (8 samples). The secondary maximum contaminant level for total dissolved solids was also exceeded for two samples. Nitrate (nitrite plus nitrate as nitrogen in groundwater) was elevated above the estimated background concentration of natural recharge waters of 1 milligram per liter (mg/L) in 44 of the 50 wells sampled and showed a median concentration of 5.4 mg/L. Nitrate concentrations were above the MCL of 10 mg/L in 5 of the 50 wells sampled and above half of the EPA MCL (5 mg/L) in 27 of the 50 wells sampled, which included samples above the MCL. Dissolved-oxygen concentrations exceeded 0.5 mg/L in 95 percent of reported values (40 of 42 samples) and exceeded 2.0 mg/L in 90 percent of reported values (38 of 42 samples). The oxidized conditions observed in most areas indicate that nitrate from fertilizers and animal or human waste was geochemically stable and could persist in the groundwater for decades or perhaps longer. A historical analysis of median nitrate concentrations over nearly three decades showed an increase in nitrate of approximately 1 mg/L from 4.3 to 5.4 mg/L, although the increase was not determined to be significantly different using nonparametric statistical methods. Major-ion data indicate that groundwater representative of the primary aquifer was classified as calcium-sodium bicarbonate type water. Other water samples from wells located mainly along the periphery of the primary aquifer had cation-anion compositions consistent with distinct water sources, including groundwater contributions from the underlying bedrock aquifers. The areas with differentiable water sources were located mainly where alluvial deposits were thin and geologic contacts to the underlying bedrock aquifers were relatively shallow. Nitrate concentrations in the groundwater were evaluated for relations to land use. An agricultural region was defined using a sequence of land satellite imagery. Groundwater flow directions interpreted from median water-table elevations measured from 2000 to 2013 were used in conjunction with cropland locations to define the agricultural region boundaries by encompassing potential pathways of nitrate transport in the groundwater from nitrogen-based fertilizers. A statistically significant higher median nitrate concentration was observed for areas inside the agricultural region (6.7 mg/L) compared to areas outside the agricultural region (2.3 mg/L), although median concentrations in both areas were below the MCL (10 mg/L). Median nitrate concentration was also significantly greater in land parcels with septic use (4.9 mg/L) compared to nonseptic parcels (1.7 mg/L). In general, agriculture or septic use was identified as the primary source of nitrate, depending on location, while commercial, county, grazing, and residential land uses were generally secondary sources of nitrate. Apparent groundwater ages were estimated from chlorofluorocarbons (CFC-11, CFC-12, and CFC-113) and tritium ( 3 H) data using models that assumed piston flow and binary mixing (dilution of a young component with old, tracer-free water). The mean and median groundwater ages were about 30 years and the standard deviation was 6 years, indicating that most groundwater in the primary aquifer was “young” water that had recharged to the aquifer over the last few decades (post-1950s). The median fraction of young water was about 71 percent, and the standard deviation was 29 percent. The remaining water predated the 1950s, which may have originated from deeper geologic formations or may represent slow moving groundwater within the primary aquifer. Some of the oldest groundwater ages (older than 30 years) were observed in the upper reaches of the aquifer to the northwest where the primary aquifer is thin and intersects bedrock, supporting the hypothesis of geochemically distinct groundwater entering the primary aquifer from below. Groundwater that had reached the central part of the aquifer from upgradient areas of the basin was variable in age because of differences in flow paths and travel velocities. The groundwater age analysis showed that current (2013) land-use practices could affect water quality over decades to come, and that responses to remedial actions could be slow, especially for constituents, such as nitrate, that are stable under oxidized conditions. Fuel products (including acetone, benzene, diisopropyl ether, ethylbenzene, methyl acetate, methyl tertiary butyl ether (MTBE), methyl tert-pentyl ether, m- + p-xylene, o-xylene, tert-amyl alcohol, tert-butyl alcohol, tert-butyl ethyl ether, and toluene) were analyzed in groundwater from 49 of the 50 wells. Water from seven sites had detections for fuel compounds; all concentrations were below MCL. The results provided assurance of water quality and a valuable baseline to evaluate future trends of fuel constituents as the region is further developed. Probability maps were developed from logistic regression models to examine the likelihood that nitrate concentrations in groundwater exceeded specified levels. Susceptibility analysis examined relations between mid-level (5.0 mg/L) nitrate concentrations and climatic, hydrologic, and geologic variables; the significant variables were identified as depth to groundwater, soil organic matter, and soil water storage to 25-centimeter (cm) depth. The vulnerability assessments included natural factors driving susceptibility but also human factors related to land use and septic use. Vulnerability to low-level (2.5 mg/L) nitrate was related to depth to groundwater, septic zoning, and soil organic matter. The results highlighted that septic zoning affected low-level nitrate concentrations. Vulnerability to mid-level (5.0 mg/L) nitrate was examined using all 50 samples and also with two data outliers removed, which showed relatively high nitrate concentrations but also anomalous water chemistry or were located beyond the primary study area. Vulnerability to mid-level (5.0 mg/L) nitrate using all 50 samples was related to depth to groundwater, land use, septic use within a 500-meter (m) radius, soil water storage to a 25-cm depth, soil organic matter, and whether a location was within the agricultural region. The mid-level (5.0 mg/L) vulnerability model using 48 samples (two outliers removed) produced the best overall fit and was related to the same variables as when using all samples except septic use. The results for mid-level vulnerability provided additional support that septic use was associated with low levels of nitrate in the groundwater. Soil properties and land use were identified as the main drivers of moderate nitrate concentrations. Probabilities of exceeding low-level nitrate concentrations were high in most areas with the lowest probabilities usually to the northwest along thin geologic deposits in the upper part of the basin. The results of this investigation offer the foundational information needed for developing best management practices to mitigate nitrate contamination, basic concepts on water quality to aid public education, and information to guide regulatory measures if policy makers determine this is warranted. Science-based decision making will require continued monitoring and analysis of water quality in the future.

Colorado↗

USGS2013

Investigations of groundwater system and simulation of regional groundwater flow for North Penn Area 7 Superfund site, Montgomery County, Pennsylvania

Groundwater in the vicinity of several industrial facilities in Upper Gwynedd Township and vicinity, Montgomery County, in southeast Pennsylvania has been shown to be contaminated with volatile organic compounds (VOCs), the most common of which is the solvent trichloroethylene (TCE). The 2-square-mile area was placed on the National Priorities List as the North Penn Area 7 Superfund site by the U.S. Environmental Protection Agency (USEPA) in 1989. The U.S. Geological Survey (USGS) conducted geophysical logging, aquifer testing, and water-level monitoring, and measured streamflows in and near North Penn Area 7 from fall 2000 through fall 2006 in a technical assistance study for the USEPA to develop an understanding of the hydrogeologic framework in the area as part of the USEPA Remedial Investigation. In addition, the USGS developed a groundwater-flow computer model based on the hydrogeologic framework to simulate regional groundwater flow and to estimate directions of groundwater flow and pathways of groundwater contaminants. The study area is underlain by Triassic- and Jurassic-age sandstones and shales of the Lockatong Formation and Brunswick Group in the Mesozoic Newark Basin. Regionally, these rocks strike northeast and dip to the northwest. The sequence of rocks form a fractured-sedimentary-rock aquifer that acts as a set of confined to partially confined layers of differing permeabilities. Depth to competent bedrock typically is less than 20 ft below land surface. The aquifer layers are recharged locally by precipitation and discharge locally to streams. The general configuration of the potentiometric surface in the aquifer is similar to topography, except in areas affected by pumping. The headwaters of Wissahickon Creek are nearby, and the stream flows southwest, parallel to strike, to bisect North Penn Area 7. Groundwater is pumped in the vicinity of North Penn Area 7 for industrial use, public supply, and residential supply. Results of field investigations by USGS at the site and results from other studies support, and are consistent with, a conceptual model of a layered leaky aquifer where the dip of the beds has a strong control on hydraulic connections in the groundwater system. Connections within and (or) parallel to bedding tend to be greater than across bedding. Transmissivities of aquifer intervals isolated by packers ranged over three orders of magnitude [from about 2.8 to 2,290 square feet per day (ft 2 /d) or 0.26 to 213 square meters per day (m 2 /d)], did not appear to differ much by mapped geologic unit, but showed some relation to depth being relatively smaller in the shallowest and deepest intervals (0 to 50 ft and more than 250 ft below land surface, respectively) compared to the intermediate depth intervals (50 to 250 ft below land surface) tested. Transmissivities estimated from multiple-observation well aquifer tests ranged from about 700 to 2,300 ft 2 /d (65 to 214 m 2 /d). Results of chemical analyses of water from isolated intervals or monitoring wells open to short sections of the aquifer show vertical differences in concentrations; chloride and silica concentrations generally were greater in shallow intervals than in deeper intervals. Chloride concentrations greater than 100 milligrams per liter (mg/L), combined with distinctive chloride/bromide ratios, indicate a different source of chloride in the western part of North Penn Area 7 than elsewhere in the site. Groundwater flow at a regional scale under steady-state conditions was simulated by use of a numerical model (MODFLOW-2000) for North Penn Area 7 with different layers representing saprolite/highly weathered rock near the surface and unweathered competent bedrock. The sedimentary formations that underlie the study area were modeled using dipping model layers for intermediate and deep zones of unweathered, fractured rock. Horizontal cell model size was 100 meters (m) by 100 meters (328 ft by 328 ft), and model layer thickness ranged from 6 m (19.7 ft) representing shallow weathered rock and saprolite up to 200 m (656 ft) representing deeper dipping bedrock. The model did not include detailed structure to account for local-scale differences in hydraulic properties, with the result that local-scale groundwater flow may not be well simulated. Additional detailed multi-well aquifer tests would be needed to establish the extent of interconnection between intervals at the local scale to address remediation of contamination at each source area. This regional groundwater-flow model was calibrated against measured groundwater levels (1996, 2000, and 2005) and base flow estimated from selected streamflow measurements by use of nonlinear-regression parameter-estimation algorithms to determine hydraulic conductivity and anisotropy of hydraulic conductivity, streambed hydraulic conductivity, and recharge during calibration periods. Results of the simulation using the calibrated regional model indicate that the aquifer appears to be anisotropic where hydraulic conductivity is greatest parallel to the orientation of bedding of the formations underlying the area and least in the cross-bed direction. The maximum hydraulic conductivity is aligned with the average regional strike of the formations, which is “subhorizontal” in the model because the altitudes of the beds and model cells vary in the strike, as well as dip, direction. Estimated subhorizontal hydraulic conductivities (in strike direction parallel to dipping beds) range from 0.001 to 1.67 meters per day (0.0032 to 5.5 feet per day). The ratio of minimum (dip direction) to maximum (strike direction) subhorizontal hydraulic conductivity ranges from 1/3.1 to 1/8.6, and the ratio of vertical to horizontal hydraulic conductivity ranges from 1/1 to 1/478. However, limited available field data precluded rigorous calibration of vertical anisotropy in the model. Estimated recharge rates corresponding to calibration periods in 1996, 2000, and 2005 are 150, 109, and 124 millimeters per year (5.9, 4.3, and 4.9 inches per year), respectively. The calibrated groundwater-flow model was used to simulate groundwater flow under steady-state conditions during periods of relatively high withdrawals (pumpage) (1990) and relatively low withdrawals (2000 and 2005). Groundwater-flow paths originating from recharge areas near known areas of soil contamination (sources) were simulated. Pumped industrial and production wells captured more groundwater from several of these sources during 1990 than after 1990 when pumping declined or ceased and greater amounts of contaminated groundwater moved away from North Penn Area 7 Superfund site to surrounding areas. Uncertainty in simulated groundwater-flow paths from contaminant sources and contributing areas, resulting from uncertainty in estimated hydraulic properties of the model, was illustrated through Monte Carlo simulations. The effect of uncertainty in the vertical anisotropy was not included in the Monte Carlo simulations. Contributing areas indicating the general configuration of groundwater flow towards production well MG-202 (L-22) in the study area also were simulated for the different time periods; as simulated, the flow paths do not pass through any identified contaminant source in North Penn Area 7. However, contributing areas to wells, such as MG-202, located near many pumped wells are particularly complex and, in some cases, include areas that contribute flow to streams that subsequently recharge the aquifer through stream loss. In these cases, water-quality constituents, including contaminants that are present in surface water may be drawn into the aquifer to nearby pumped wells. Results of a simulated shutdown of well MG-202 under steady-state 2005 conditions showed that the area contributing recharge for nearby production well MG-76 (L-17), when MG-202 is not pumping, shifts downstream and is similar to the area contributing recharge for MG-202 when both wells are pumping. Concentrations of constituents in groundwater samples collected in fall 2005 or spring 2006 were compared to simulated groundwater-flow paths for the year 2005 to provide a qualitative assessment of model results. The observed spatial distribution of selected constituents, including TCE, CFC-11, and CFC-113 in groundwater in 2005 and the chloride/bromide mass ratios in 2006, generally were consistent with the model results of the simulated 2005 groundwater-flow paths at North Penn Area 7, indicating the presence of several separate sources of contaminants within North Penn Area 7.

Pennsylvania↗