Search USGS⌕ Search

SEARCH · Search USGS

Results for “Regional Environmental Change”

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

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

1,404 records · Page 78Linked to original sources

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↗

Geochemical conditions and nitrogen transport in nearshore groundwater and the subterranean estuary at a Cape Cod embayment, East Falmouth, Massachusetts, 2013–14

Nitrogen transport and transformation were studied during 2013 to 2014 by the U.S. Geological Survey, in cooperation with the U.S. Environmental Protection Agency, in a subterranean estuary beneath onshore locations on the Seacoast Shores peninsula, a residential area in Falmouth, Massachusetts, served by septic systems and cesspools, and adjacent offshore locations in the Eel River, a saltwater embayment connected to the ocean. The field investigation included installation and sampling of clusters of wells and temporary sampling points near a transect extending from about 35 meters (m) onshore to 18 m offshore. The fresh groundwater at the study site formed a lens about 11 m thick at the shoreline that was underlain by saline groundwater. Groundwater flow in the water-table aquifer was oriented northwestward toward the embayment. Nitrate concentrations in the fresh groundwater at a site about 35 m onshore increased in the downward direction from less than 500 micromoles per liter near the water table to about 1,700 micromoles per liter just above the freshwater/saltwater transition zone. Dissolved oxygen was largely absent in the onshore fresh groundwater. Distributions of salinity, dissolved oxygen, and nitrate at the shoreline and offshore generally were similar to those onshore; at some locations, however, shallow saline water was present above the freshwater, and there were scattered occurrences of elevated dissolved oxygen concentrations. Geochemical indicators of nitrate reduction, including concentrations of the reaction product nitrogen gas, stable isotope ratios of nitrate and nitrogen gas, and changes in alkalinity, provided evidence for nitrate reduction in two zones separated vertically by a zone 7–8 m thick with no evidence of nitrate reduction. The shallow nitrate-reduction zone was near the water table in fresh groundwater onshore, where nitrate reduction may be related to particular recharge conditions at nearby sources. The shallow nitrate-reduction zone also may be related to an interval of fine-grained sediments at about the same altitude (−1 to −6 m relative to the National Geodetic Vertical Datum of 1929), where flow is slower and reactive electron donors such as solid organic carbon, iron, or sulfide phases may be present to drive the reduction. The deep nitrate-reduction zone was near the freshwater/saltwater transition zone, where nitrate reduction may be related to mixing of freshwater containing nitrate and saltwater containing dissolved organic carbon and ammonium, or to fine-grained sediments near the transition zone. The maximum amount of nitrate converted to nitrogen gas was estimated to be less than or equal to 300 micromoles per liter in both nitrate-reduction zones. The presence of nitrate and low dissolved oxygen concentrations in the 7–8-meter-thick zone between the shallow and deep nitrate-reduction zones are conditions that could permit nitrate reduction. The absence of evidence of nitrate reduction in the high-nitrate zone may have resulted from the lack of reactive electron donors in that depth interval. The high-nitrate zone dissipated somewhat in the offshore direction, but the current study did not extend far enough to encompass the fresh groundwater discharge area or determine how much of the nitrate was removed prior to discharge. A shallow intertidal saltwater cell was formed during a spring tide by saltwater infiltration during tidal run-up on the beach. Nitrate reduction might have occurred if nitrate-containing fresh groundwater discharging to the estuary mixed with the saltwater containing dissolved organic carbon in this zone, but samples collected from the intertidal saltwater cell during this study were not analyzed for indicators of nitrate reduction. Elevated dissolved oxygen concentrations in fresh groundwater 9 m offshore may indicate that groundwater flow was partly oblique to the sampling transect or that groundwater from a regional flow system was converging under the river near the study area. Flow directions also may have been affected by aquifer heterogeneity such as the shallow fine-grained sediments onshore and at the bottom of the Eel River. Improved understanding of the fate of nitrate in this type of complex setting might be gained by including additional characterization of aquifer heterogeneity and groundwater flow and extending investigations of nitrate reduction to the shallow sediments in the intertidal saltwater cell and adjacent subtidal zone and to locations farther offshore beneath the estuary.

Massachusetts↗

Recycling ground water in Waushara County, Wisconsin : resource management for cold-water fish hatcheries

Recycling water within the local ground-water system can increase the quantity of water available for use, control or avoid environmental pollution, and control temperature of the water supply. Pumped ground water supplied a fish-rearing facility for 15 months, and the waste water recharged the local ground-water system through an infiltration pond. Eighty-three percent of the recharged water returned to the well (recycled). Make-up water from the ground-water system provided the remaining 17 percent. Pumping 300 gallons per minute (20 litres per second), combined with recycling, resulted in water-level declines equivalent to a pumping rate of approximately 50 gallons per minute (3 litres per second). Using this effective pumping rate in the Theis nonequilibrium equation resulted in predicted drawdowns within 0.5 foot (0.2 metres) of those observed throughout the 15-month period. The concentration of nitrate in the water supply increased only slightly during the 15 months of recycling. Nitrate levels in a closed recycling system (100 percent recycling efficiency) would have reached 9 milligrams per litre, but observed levels did not exceed 4 milligrams per litre, and at the end of the recycling period they were lower than the initial levels. Mass-balance equations relate observed nitrate levels to the loading imposed on the system, the pumping rate, the volume of ground water affected by recycling, and the recycling efficiency. The equations predict nitrate concentrations (or other ions not attenuated by movement through the unsaturated part of the aquifer) within 1 milligram per litre of observed concentrations except during periods when nutrients are used in plant growth. The method does not account for nutrients utilized by aquatic vegetation in the infiltration pond, so that observed levels would usually be even lower than predicted. The predicted response of the local groundwater system to a nutrient loading equivalent to that generated by a hatchery producing approximately 100,000 pounds (50,000 kilograms) of trout and salmon per year indicates that maximum nitrate levels would remain significantly below the limit established by the State of Wisconsin (and the U.S. Public Health Service, 1962) for drinking water. The water-supply temperature can be maintained within the optimum range for trout and salmon rearing (10.0° to 15.5°C or 50° to 60°F) during recycling. Continuous recycling during the study period resulted in watersupply temperatures ranging between 7.0° and 14.0°C (45° and 57°F). Longterm continuous recycling would result in water-supply temperatures ranging between 7.0° and 14.5°C (45° and 58°F). Selective recycling (recycling water for only 8 months of the year) would provide water-supply temperatures ranging from 9.5° to 15.5°C (49° to 60°F). A permanent recharge pond with supplementary ponds would fulfill the needs of a hatchery development at the Greenwood Wildlife Refuge site, Waushara County, Wisconsin, and insure protection of the ground-water system. Eighty percent or more of the water pumped could be recycled by recharging waste water near the supply well. Selective recycling (recharging water to the ground-water system outside the zone of recycling during approximately 4 months of each year) could maintain optimum water-supply temperatures, reduce water-level declines by 50 percent (compared to no recycling), maintain nitrate levels in the water supply below limits established for drinking water supplies, and minimize the effect of water-supply development on the regional ground-water system. Other recharge-recycling schemes can also be evaluated. Estimating the recycling efficiency (of recharge ponds, trenches, spreading areas, or irrigated fields) provides a basis for predicting water-level declines, the concentration of conservative ions (conservative in the sense that no reaction other than mixing occurs to change the character of the ion being considered) in the water supply and in the regional ground-water system, and the temperature of the water supply. Hatchery development and management schemes can be chosen to optimize hatchery productivity or minimize operation costs while protecting the ground-water system.

Wisconsin↗

Hydrogeology and hydrology of the Punta Cabullones wetland area, Ponce, southern Puerto Rico, 2007-08

The U.S. Geological Survey, in cooperation with the Municipio Autónomo de Ponce and the Puerto Rico Department of Natural and Environmental Resources, conducted a study of the hydrogeology and hydrology of the Punta Cabullones area in Ponce, southern Puerto Rico. (Punta Cabullones is also referred to as Punta Cabullón.) The Punta Cabullones area is about 9 square miles and is an ecological system made up of a wetland, tidal flats, saltflats, mangrove forests, and a small fringing reef located a short distance offshore. The swales or depressions between successive beach ridges became development avenues for saline to hypersaline wetlands. The Punta Cabullones area was designated by the U.S. Fish and Wildlife Service as a coastal barrier in the 1980s because of its capacity to act as a buffer zone to ameliorate the impacts of natural phenomenon such as storm surges. Since 2003, Punta Cabullones has been set aside for preservation as part of the mitigation effort mandated by Federal and State laws to compensate for the potential environmental effects that might be caused by the construction of the Las Américas Transshipment Port. Total rainfall measured during 2008 within the Punta Cabullones area was 36 inches, which is slightly greater than the long-term annual average of 32 inches for the coastal plain near Ponce. Two evapotranspiration estimates, 29 and 37 inches, were obtained for the subarea of the Punta Cabullones area that is underlain by fan-delta and alluvial deposits by using two variants of the Penman semi-empirical equation. The long-term water stage and chemical character of the wetland in Punta Cabullones are highly dependent on the seasonal and annual variations of both rainfall and sea-wave activity. Also, unseasonal short-term above-normal rainfall and sea-wave events resulting from passing storms may induce substantial changes in the water stage and the chemical character of the wetland. In general, tidal fluctuations exert a minor role in modifying the water quality and stage of the wetland in Punta Cabullones. The role of the tidal fluctuations becomes important during those times when the outlets/inlets to the sea are not blocked by a sand bar and is allowed to freely flow into the wetland interior. The salinity of the wetland varies from brackish to hypersaline. The hypersaline conditions, including the occurrence of saltflats, within the Punta Cabullones wetland area result from a high evapotranspiration rate. The hypersaline conditions are further enhanced by a sand bar that blocks the inlet/outlet of the wetland’s easternmost channel, particularly during the dry season. Groundwater in Punta Cabullones mostly is present within beds of silisiclastic sand and gravel. During the study period, the depth to groundwater did not exceed 4 feet below land surface. The movement and direction of the groundwater flow in Punta Cabullones are driven by density variations that in turn result from the wide range of salinities in the groundwater. The salinity of the groundwater decreases within the first 60 to 100 feet of depth and decreases outward from a mound of hypersaline groundwater centered on piezometer nest PN2. The main groundwater types within the Punta Cabullones area vary from calcium-bicarbonate type in the northernmost part of the study area to a predominantly sodium-potassium-chloride groundwater type southward. According to stable-isotope data, groundwater within the study area is both modern meteoric water and seawater highly affected by evaporation. The chemical and stable-isotopic character of local groundwater is highly influenced by evapotranspiration because of its shallow depth. Equivalent freshwater heads indicate groundwater moves away from a mound centered on piezometer nest PN2, in a pattern similar to the spatial distribution of groundwater salinity. Vertical groundwater flow occurs in Punta Cabullones due to local differences in density. In the wetland subarea of Punta Cabullones, groundwater and surface water are hydraulically coupled. Locally, surface-hypersaline water sinks into the aquifer, providing recharge and serving as a mechanism to redistribute salt throughout the study area. The evapotranspiration in the wetland subarea is estimated at about 11 million gallons per day (Mgal/d) that is equivalent to about 12,586 acre-feet per year. The balance of evapotranspiration, in excess of the about 0.5 Mgal/d of groundwater flow within the wetland, is supplied by saline to hypersaline surface water that may include seawater and meteoric water highly affected by evaporation with dissolved salts. In one of the extreme scenarios in which no groundwater is intercepted by pumpage at the Restaurada well field, the amount of saline to hypersaline water in the wetland consumed by evapotranspiration is about 10.5 Mgal/d. In the opposite extreme in which the entire regional groundwater flow is intercepted by pumpage in the Restaurada well field, the entire evapotranpiration requirement is met by saline to hypersaline water. Hydrologic, isotopic, and chemical data indicate that all of, or a large portion of, the historical groundwater flow to Punta Cabullones is being captured by the Puerto Rico Aqueducts and Sewer Authority pumpage at the Restaurada well field at a rate of about 2 Mgal/d. As a consequence, seawater intrusion into the aquifer at the Punta Cabullones area seems to be occurring, while the current pumpage at the Restaurada well field is sustained by storage depletion of the aquifer.

Puerto Rico↗

Hydrologic analyses in support of the Navajo Generating Station–Kayenta Mine Complex environmental impact statement

Introduction The U.S. Department of Interior’s Bureau of Reclamation, Lower Colorado Region (Reclamation) is preparing an environmental impact statement (EIS) for the Navajo Generating Station-Kayenta Mine Complex Project (NGS-KMC Project). The proposed project involves various Federal approvals that would facilitate continued operation of the Navajo Generating Station (NGS) from December 23, 2019 through 2044, and continued operation of the Kayenta Mine and support facilities (collectively called the Kayenta Mine Complex, or KMC) to supply coal to the NGS for this operational period. The EIS will consider several project alternatives that are likely to produce different effects on the Navajo (N) aquifer; the N aquifer is the principal water resource in the Black Mesa area used by the Navajo Nation, Hopi Tribe, and Peabody Western Coal Company (PWCC). The N aquifer is composed of three hydraulically connected formations—the Navajo Sandstone, the Kayenta Formation, and the Lukachukai Member of the Wingate Sandstone—that function as a single aquifer. The N aquifer is confined under most of Black Mesa, and the overlying stratigraphy limits recharge to this part of the aquifer. The N aquifer is unconfined in areas surrounding Black Mesa, and most recharge occurs where the Navajo Sandstone is exposed in the area near Shonto, Arizona. Overlying the N aquifer is the D aquifer, which includes the Dakota Sandstone, Morrison Formation, Entrada Sandstone, and Carmel Formation. The aquifer is named for the Dakota Sandstone, which is the primary water-bearing unit. The NGS is located near Page, Arizona on the Navajo Nation. The KMC, which delivers coal to NGS by way of a dedicated electric railroad, is located approximately 83 miles southeast of NGS (about 125 miles northeast of Flagstaff, Arizona). The Kayenta Mine permit area is located on about 44,073 acres of land leased within the boundaries of the Hopi and Navajo Indian Reservations. KMC has been conducting mining and reclamation operations within the Kayenta Mine permit boundary since 1973. The KMC part of the proposed project requires approval by the Office of Surface Mining (OSM) of a significant revision of the mine’s permit to operate in accordance with the Surface Mine Control and Reclamation Act (Public Law 95-87, 91 Stat. 445 [30 U.S.C. 1201 et seq. ]) . The revision will identify coal resource areas that may be used to continue extracting coal at the present rate of approximately 8.2 million tons per year. The Kayenta Mine Complex uses water pumped from the D and N aquifers beneath PWCC’s leasehold to support mining and reclamation activities. Prior to 2006, water from the PWCC well field also was used to transport coal by way of a coal-slurry pipeline to the now-closed Mohave Generating Station. Water usage at the leasehold was approximately 4,100 acre-feet per year (acre-ft/yr) during the period the pipeline was in use, and declined to an average 1,255 acre-ft/yr from 2006 to 2011. The Probable Hydrologic Consequences (PHC) section of the mining and reclamation permit must be modified to project the consequences of extended water use by the mine for the duration of the KMC part of the project, including a post-mining reclamation period. Since 1971, the U.S. Geological Survey (USGS) has conducted the Black Mesa Monitoring Program, which consists of monitoring water levels and water quality in the N aquifer, compiling information on water use by PWCC and tribal communities, maintaining several stream-gaging stations, measuring discharge at selected springs, conducting special studies, and reporting findings. These data are useful in evaluating the effects on the N aquifer from PWCC and community pumping, and the effects of variable precipitation. The EIS will assess the impacts of continued pumping on the N aquifer, including changes in storage, water quality, and effects on spring and baseflow discharge, by proposed mining through 2044, and during the reclamation process to 2057. Several groundwater models exist for the area and Reclamation concluded it would conduct a peer review of the groundwater flow model that will be used to assess the direct, reasonably foreseeable indirect, and cumulative effects of future groundwater withdrawals on the D and N aquifers in the Black Mesa area. Reclamation made this determination because of the level of controversy around the effects of continued water use and the comments received from the 2014 draft EIS scoping meetings. Reclamation requested assistance from the USGS in evaluating existing groundwater flow models of the Black Mesa Basin that can be used to predict the effects of different project alternatives on the D and N aquifers.

Arizona↗

Impacts of exploratory drilling for oil and gas on the benthic environment of Georges Bank

A 3-year monitoring program was performed to assess the impacts of exploratory drilling for oil and gas on the benthic environment of Georges Bank, an important commercial fishery region in the North Atlantic east of Massachusetts, USA. Surficial sediments were sampled for chemical and benthic infaunal analysis and bottom still photographs were taken to document bottom microtopography and epifauna at 46 stations during 12 field surveys. The surveys were performed quarterly from just before drilling began, during drilling, and for nearly 2 years after completion of drilling. Two of the eight drilling sites were selected for monitoring. Twenty-nine stations were positioned in a tight radial array around a drilling site in 80 m of water. A second group of three stations was positioned near another drilling site in 140 m of water. The remaining stations covered a broad expanse of the Bank and adjacent suspected sites of deposition of fine-grained sediments.\ Of the 12 elements analyzed in bulk sediments, only barium increased in concentration during the period when drilling was taking place (July 1981 to September 1982). The concentration of barium in surficial sediment increased 4·7-fold from 28 ppm before drilling to 131·7 ppm after drilling at the station closest to the drilling site in 80 m of water and 5·9-fold from 32 ppm before drilling to 172 ppm after drilling at the station closest to the drilling site in 140 m of water. The concentrations of both barium and chromium increased in the fine (silt/clay) fraction (usually less than 5% by weight of sediment from most stations) of sediments from several stations around one or both rig sites monitored during the period of drilling. Elevated concentrations of chromium (about two-fold) occurred only in sediments near the drilling site in 140 m of water. Statistically significant increases in the concentration of barium in the fine fraction to sediment were detected approximately 65 km west (downcurrent) and 35 km east of the drilling site in 80 m of water after drilling was completed. The benthic fauna were abundant and diverse throughout the study area. At most stations, the dominant species remained nearly the same from one season to another over the 3 years of sampling. Polychaetes were the most abundant, followed by crustaceans. The number of individuals of some species, particularly the amphipods Erichthonius fasciatus and Unciola inermis , showed large seasonal variations. Cluster analysis revealed a strong relationship between community structure and both sediment type and water depth. Little seasonal variation was detected, but some interannual differences were revealed by cluster analysis and correspondence analysis. The replicates from a station always resembled each other more than they resembled any replicates from other stations. In addition, the combined replicates from a station always clustered with samples from that station taken on other cruises. This excellent replication and uniformity of the benthic infaunal community at a station over time made it possible to detect very subtle changes in community parameters that might be related to discharges of drilling fluid and drill cuttings. Nevertheless, no changes were detected in benthic communities of Georges Bank that could be attributed to drilling activities.

Massachusetts↗

Water-quality assessment of the Connecticut, Housatonic, and Thames River Basins study unit: Analysis of available data on nutrients, suspended sediments, and pesticides, 1972-92

This retrospective report examines available nutrient, suspended sediment, and pesticide data in surface and ground water in the Connecticut, Housatonic and Thames Rivers Study Unit of the National Water-Quality Assessment Program. The purpose of this study is to improve the understanding of natural and anthropogenic factors affecting water quality in the study unit. Waterquality data were acquired from various sources, primarily, the U.S. Geological Survey and the U.S. Environmental Protection Agency. The report examines data for water years 1972-92, focusing on 1980-92, although it also includes additional data from as early as 1905. The study unit lies within the New England Physiographic Province and altitudes range from sea level in coastal Connecticut to 6,288 feet above sea level at Mount Washington, New Hampshire. Two major aquifer types underlie the study unit unconsolidated glacial deposits and fractured bedrock. The climate generally is temperate and humid, with four distinct seasons. Average annual precipitation ranges from 34 to 65 inches. The study unit has a population of about 4.5 million, which is most highly concentrated in southwestern Connecticut and along the south-central region of the Connecticut River Valley. Surface-water-quality data were screened to provide information about sites with adequate numbers of analyses (50) over sufficiently long periods (1980-90) to enable valid statistical analyses. In order to compare effects of different types of land use on surface-water quality, examination of data required application of several statistical and graphical techniques, including mapping, histograms, boxplots, concentration-discharge plots, trend analysis, and load estimation. Spatial and temporal analysis of surface-water-quality data indicated that, with a single exception, only stations in the Connecticut water-quality network had sufficient data collected over adequately long time periods to use in detailed analyses. Ground-water nutrient and pesticide data were compiled from several Federal and State agencies, primarily the U.S. Geological Survey, U.S. Environmental Protection Agency, and Connecticut Department of Health Services. Nutrient data were available for several thousand wells; nitrite plus nitrate as nitrogen was the most commonly reported constituent. Most wells with nutrient data are in Massachusetts and Connecticut. Relative to nutrient data in ground and surface water, pesticide data are less common. Pesticide data were available for slightly more than 200 surface-water sites and less than 500 wells; about 95 percent of the wells are completed in stratified-drift or till aquifers. Data for 81 pesticide compounds were available in various data bases. 2,4-D and silvex were the most commonly detected herbicides in surface water and dieldrin and diazinon were the most commonly detected insecticides. Most surface-water pesticide samples and detections are from bed sediment, but much of the data are not recent. Ethylene dibromide (EDB), a soil fumigant used in tobacco farming was detected in 268 wells in a 50 square-mile area of north-central Connecticut; EDB contamination also was detected in wells in Massachusetts. Atrazine, an herbicide commonly used in corn farming, commonly was detected in wells installed in tilled agricultural fields. Corn herbicides were commonly detected in the northern part of the study unit, although the sampling has been less frequent than in the southern part of the study unit. Pesticides were seldom detected in public-supply wells in Connecticut. Urban sites with relatively high population densities and high concentrations of dischargers were characterized by having the highest nutrient concentrations and loads when adjusted for differences in drainage area or population. Particularly, the Pequabuck, Naugatuck, and Quinnipiac River Basins were characterized by high nutrient concentrations median total nitrogen concentrations ranged from 3.3 to 4.2 mg/L (milligrams per liter) and median total phosphorus concentrations ranged from 0.42 to 0.8 mg/L. In contrast, the predominantly forested and low density residential land-use sites, such as Saugatuck and Salmon River Basins, were characterized by low nutrient concentrations median total nitrogen ranged from 0.50 to 0.60 mg/L and median total phosphorus concentrations ranged from 0.01 to 0.02 mg/L. Estimated total nitrogen loadings in median discharge years ranged from 940 kilograms per square mile at the Salmon River near East Hampton, Conn., to 5,800 kilograms per square mile at the Naugatuck River at Beacon Falls, Conn. Water quality, in terms of nutrient concentrations and areally adjusted loadings, for sites with large drainage basins integrating a wide variety of land-use categories fell between the extremes of the urban and forested sites total nitrogen was 1,400 kilograms per square mile per year at the Connecticut River at Thompsonville, Conn. Nitrate concentrations in ground water occasionally exceeded the safe drinking-water standard of 10 mg/L as nitrogen. The greatest number of detections exceeding the standard, however, were not in public-water supplies but in shallow observation wells in agricultural settings (the most frequently sampled type of well). None of the public-supply wells in Massachusetts exceeded the standard. Although nitrate concentrations for Vermont and New Hampshire generally were low, few data were available and those were seldom reported on the basis of drainage basin, making analysis difficult. Trend analysis indicated that flow-adjusted concentrations of total and dissolved phosphorus generally decreased during the period of analysis, however, total nitrogen did not change substantially. Decreases in ammonia concentrations with time were usually accompanied by increases in nitrate, suggesting improvements in sewage treatment. The lack of adequate data from more or less exclusively agricultural areas points to the need for further study of the effects of fanning on surface-water quality in the study unit. Furthermore, additional information is needed on the rates, transformations, and movements of nutrients and other materials through and between the aquatic and terrestrial components of the study unit.

Connecticut, Massachusetts, New Hampshire, Quebec,↗

National Acid Precipitation Assessment Program Report to Congress: An integrated assessment

Acid deposition, more commonly known as acid rain, occurs when emissions of sulfur dioxide (SO2) and nitrogen oxides (NOx) react in the atmosphere (with water, oxygen, and oxidants) to form various acidic compounds. These acidic compounds then fall to earth in either a wet form (rain, snow, and fog) or a dry form (gases, aerosols, and particles). Prevailing winds transport the acidic compounds hundreds of miles, often across state and national borders. At certain levels the acidic compounds, including small particles such as sulfates and nitrates can cause many negative human health and environmental effects. While ecosystems are subject to many stresses, including land-use changes, climate change, and variations in hydrologic and meteorologic cycles, the scientific literature has clearly demonstrated that these pollutants can: • Degrade air quality, • Impair visibility, • Damage public health, • Acidify lakes and streams, • Harm sensitive forests, • Harm sensitive coastal ecosystems, and • Accelerate the decay of building materials, paints, and cultural artifacts such as buildings, statues, and sculptures. Title IV was passed by Congress as part of the 1990 Clean Air Act Amendments to reduce emissions of SO 2 and NO x from fossil fuel-burning power plants in order to protect ecosystems suffering damage from acid deposition and to improve air quality. At the same time, the National Acid Precipitation Assessment Program (NAPAP) was asked to periodically assess and report to Congress on the implementation of the Acid Rain Program, recent scientific knowledge surrounding acid deposition and its effects, and the reduction in acid deposition necessary to prevent adverse ecological effects. This NAPAP Report focuses primarily on emission reductions from power plants, both in terms of assessing past reductions under the Acid Rain Program and in projecting the ecological effects of additional reductions of SO 2 and NO x . It should be noted that power generation currently contributes approximately 69% of the SO 2 emissions and 22% of the NO x emissions nationwide. This contribution is decreasing as emissions from power generation continue to decrease, making the other sources of these pollutants more prominent. Modeling suggests that even if SO 2 emissions from power generation were reduced to zero, some lakes and streams would remain acidic due to acid deposition. However, there are several other regulations that reduce emissions of SO 2 and NO x from these non-power generation sources, such as the Tier II mobile source standards, the Heavy Duty Diesel standards, and the Clean Air Non-Road Diesel Rule, that have also been promulgated since 1990. These regulations, primarily designed to bring counties into attainment with fine particle and ozone air quality standards, also incidentally reduce emissions that contribute to acid deposition. Implementation of Title IV has successfully and substantially reduced emissions of SO 2 and NO x from power generation at a significantly lower cost than expected: • In 2002, SO 2 emissions were 10.2 million tons, 35% lower than 1990 emissions and 40% lower than 1980 emissions.* • In 2002, NO x emissions were 4.5 million tons, 33% lower than 1990 emissions. In addition, SO 2 emissions from all sources have decreased by 32% since 1990 and emissions of NO x from all sources have decreased by 12% since 1990. Power generating sources continue to close in on the goal of reducing power plant SO 2 emissions from 1980 levels by 50% (to 8.95 million tons) as required by the 1990 Clean Air Act. Power generating sources have also exceed­ed the goal of a two million ton reduction in NO x emissions from projected 2000 levels as required by the 1990 Clean Air Act. These emission reductions have contributed to measurable improvements in air quality, reductions in acid deposition, and the beginnings of recovery of acid-sensitive waters in some areas: • SO2 concentrations in the atmosphere (a precursor to fine particles and acid deposition) have decreased since 1990. Average annual SO2 concentrations in the Northeast in 2000–2002 were 40% lower than they were in 1989–1991, concentrations in the mid-Atlantic were 30% lower, concentrations in the Southeast were 35% lower, and concentrations in the Midwest were 45% lower. • Sulfate concentrations in the atmosphere (a major component of fine particles, especially in the East) have decreased since 1990 as well. Average annual sulfate concentrations in the Northeast and Midwest in 2000–2002 were approximately 30% lower than they were in 1989–1991, and concentrations in the mid-Atlantic and Southeast were 25% lower. • Wet sulfate deposition, a major component of acid rain, has also decreased since 1990. Average annual sulfate deposition in the Northeast in 2000–2002 was 40% lower than it was in 1989–1991, deposition in the mid-Atlantic and Midwest was 35% lower, and deposition in the Southeast was 25% lower. • Wet nitrate deposition has not decreased regionally from historical levels because of the relatively moderate NO x reduction from power plants and the continuing large contribution (over 50% of total NOx emissions) from other sources of NO x such as vehicles and nonroad vehicles. • Although visibility has begun to improve in some parts of the U.S., there is still significant impairment of visibility in many national parks and other Class I areas throughout the U.S. • Acid neutralizing capacity is beginning to rise in some surface waters in the Northeast, including lakes in the Adirondack Mountains (see graphic below). This is an indication that recovery from acidification is occurring in those areas.

Report to Congress↗

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↗

U.S. Geological Survey Karst Interest Group Proceedings, Carlsbad, New Mexico, April 29-May 2, 2014

Karst aquifer systems are present throughout parts of the United States and some of its territories, and have developed in carbonate rocks (primarily limestone and dolomite) that span an interval of time encompassing more than 550 million years. The depositional environments, diagenetic processes, post-depositional tectonic events, and geochemical weathering processes that form karst aquifers are varied and complex, and involve biological, chemical, and physical changes. These factors, combined with the diverse climatic regimes under which karst development in these rocks has taken place, result in the unique dual- or triple-porosity nature of karst aquifers. These complex hydrogeologic systems typically represent challenging and unique conditions to scientists attempting to study groundwater flow and contaminant transport in these terrains. The dissolution of carbonate rocks and the subsequent development of distinct and beautiful landscapes, caverns, and springs has resulted in the most exceptional karst areas of the United States being designated as national or state parks; commercial caverns and known privately owned caves number in the tens of thousands. Both public and private properties provide access for scientists to study the flow of groundwater in situ. Likewise, the range and complexity of landforms and groundwater flow systems associated with karst terrains are enormous, perhaps more than for any other aquifer type. Karst aquifers and landscapes that form in tropical areas, such as the cockpit karst along the north coast of Puerto Rico, differ greatly from karst landforms in more arid climates, such as the Edwards Plateau in west-central Texas or the Guadalupe Mountains near Carlsbad, New Mexico, where hypogenic processes have played a major role in speleogenesis. Many of these public and private lands also contain unique flora and fauna associated with these karst hydrogeologic systems. As a result, numerous federal, state, and local agencies have a strong interest in the study of karst terrains. Many of the major springs and aquifers in the United States have developed in carbonate rocks, such as the Floridan aquifer system in Florida and parts of Alabama, Georgia, and South Carolina; the Ozark Plateaus aquifer system in parts of Arkansas, Kansas, Missouri, and Oklahoma; and the Edwards-Trinity aquifer system in west-central Texas. These aquifers, and the springs that discharge from them, serve as major water-supply sources and as unique ecological habitats. Competition for the water resources of karst aquifers is common, and urban development and the lack of attenuation of contaminants in karst areas can impact the ecosystem and water quality of these aquifers. The concept for developing a platform for interaction among scientists within the U.S. Geological Survey (USGS) working on karst-related studies evolved from the November 1999 National Ground-Water Meeting of the USGS. As a result, the Karst Interest Group (KIG) was formed in 2000. The KIG is a loose-knit, grass-roots organization of USGS and non-USGS scientists and researchers devoted to fostering better communication among scientists working on, or interested in, karst science. The primary mission of the KIG is to encourage and support interdisciplinary collaboration and technology transfer among scientists working in karst areas. Additionally, the KIG encourages collaborative studies between the different mission areas of the USGS as well as other federal and state agencies, and with researchers from academia and institutes. The KIG also encourages younger scientists by participation of students in the poster and oral sessions. To accomplish its mission, the KIG has organized a series of workshops that are held near nationally important karst areas. To date (2014) six KIG workshops, including the workshop documented in this report, have been held. The workshops typically include oral and poster sessions on selected karst-related topics and research, as well as field trips to local karst features. Proceedings of the workshops are published by the USGS and are available online at http://water.usgs.gov/ogw/karst/kig . The first KIG workshop was held in St. Petersburg, Florida, February 13–16, 2001, in the vicinity of the large springs and other karst features of the Floridan aquifer system. The second KIG workshop was held August 20–22, 2002, in Shepherdstown, West Virginia, in proximity to the carbonate aquifers of the northern Shenandoah Valley and highlighted an invited presentation on karst literature by the late Barry F. Beck of P.E. LaMoreaux and Associates. The third KIG workshop was held September 12–15, 2005, in Rapid City, South Dakota, nearby to karst features in evaporites and limestones of the Madison Group in the Black Hills of South Dakota, including Wind Cave National Park and Jewel Cave National Monument. The workshop also included a featured presentation by Thomas Casadevall, Central Region Director, USGS, on the status of earth science at the USGS and evening trips to Jewel Cave led by Mike Wiles, National Park Service (NPS) and Wind Cave led by Rod Horrocks, NPS. The fourth KIG workshop was held May 27–29, 2008, and hosted by the Hoffman Environmental Research Institute and Center for Cave and Karst Studies at Western Kentucky University in Bowling Green, Kentucky, near Mammoth Cave National Park and karst features of the Chester Upland and Pennyroyal Plateau. The workshop featured a late-night field trip into Mammoth Cave with Rickard Toomey and Rick Olsen, NPS. The fifth workshop was held April 26–29, 2011, and was a joint meeting of the USGS KIG and University of Arkansas HydroDays, hosted by the Department of Geosciences at the University of Arkansas in Fayetteville. The workshop featured an outstanding field trip to the unique karst terrain along the Buffalo National River of the southern Ozarks and a keynote presentation on paleokarst in the United States by Art and Peggy Palmer. This sixth and current 2014 KIG workshop is hosted by the National Cave and Karst Research Institute (NCKRI) in Carlsbad, New Mexico, with Director of NCKRI, George Veni, serving as co-chair of the workshop with Eve Kuniansky, USGS. The session planning committee for this sixth workshop includes Van Brahana, USGS retired and University of Arkansas Professor Emeritus; Tom Byl, USGS and Tennessee State University; Zelda Bailey, former Director of NCKRI and retired Director, National Institute of Standards and Technology, Boulder Laboratory, Colorado; Patrick Tucci, USGS retired; and Mike Bradley, Allan Clark, Geoff Delin, Daniel Doctor, James Kaufmann, Eve Kuniansky, Randy Orndorff, Larry Spangler, and Dave Weary of the USGS. The karst hydrology field trip on Thursday will be led by Lewis Land (NCKRI karst hydrologist) and the optional Friday field trip on the geology of Carlsbad Caverns National Park will be led by George Veni. The keynote speaker is Dr. Penelope Boston, Director of Cave and Karst Studies at New Mexico Tech, Socorro, and Academic Director at NCKRI, who will address the future of karst research. Additionally, there is a featured presentation “Irish karst and its management,” by Caoimhe Hickey, The Geological Survey of Ireland, preceding a panel discussion on “Collaboration During Times of Limited Resources.” The extended abstracts of USGS authors were peer reviewed and approved for publication by the U.S. Geological Survey. Articles submitted by university researchers and other federal and state agencies did not go through the formal USGS peer review and approval process, and therefore may not adhere to our editorial standards or stratigraphic nomenclature and is not research conducted or data collected by the USGS. However, all articles had at a minimum of two peer reviews, and all articles were edited for consistency of appearance in the published Proceedings. The use of trade, firm or product names in any article is for descriptive purposes only and does not imply endorsement by the U.S. Government. The USGS, Office of Groundwater, provides technical support for the Karst Interest Group website and public availability of the Proceedings from these workshops, and the USGS Groundwater Resources Program funds the publication costs. Finally, the cover illustration is the work of Ann Tihansky, USGS, used since the first KIG workshop in 2000.

Scientific Investigations Report↗

U.S. Geological Survey Karst Interest Group Proceedings, San Antonio, Texas, May 16–18, 2017

Introduction and Acknowledgments Karst aquifer systems are present throughout parts of the United States and some of its territories, and have developed in carbonate rocks (primarily limestone and dolomite) and evaporites (gypsum, anhydrite, and halite) that span an interval of time encompassing more than 550 million years. The depositional environments, diagenetic processes, post-depositional tectonic events, and geochemical weathering processes that form karst aquifers are varied and complex. These factors involve biological, chemical, and physical changes that when combined with the diverse climatic regimes in which karst development has taken place, result in the unique dual- or triple-porosity nature of karst aquifers. These complex hydrogeologic systems typically represent challenging and unique conditions to scientists attempting to study groundwater flow and contaminant transport in these terrains. The dissolution of carbonate rocks and the subsequent development of distinct and beautiful landscapes, caverns, and springs have resulted in the most exceptional karst areas being designated as national or state parks. Tens of thousands of similar areas in the United States have been developed into commercial caverns and known privately owned caves. Both public and private properties provide access for scientists to study the flow of groundwater in situ . Likewise, the range and complexity of landforms and groundwater flow systems associated with karst terrains are enormous, perhaps more than for any other aquifer type. Karst aquifers and landscapes that form in tropical areas, such as the cockpit karst along the north coast of Puerto Rico, differ greatly from karst landforms in more arid climates, such as the Edwards Plateau in west-central Texas or the Guadalupe Mountains near Carlsbad, New Mexico, where hypogenic processes have played a major role in speleogenesis. Many of these public and private lands also contain unique flora and fauna associated with these karst hydrogeologic systems. As a result, numerous federal, state, and local agencies have a strong interest in the study of karst terrains. Many of the major springs and aquifers in the United States have developed in carbonate rocks, such as the Floridan aquifer system in Florida and parts of Alabama, Georgia, and South Carolina; the Ozark Plateaus aquifer system in parts of Arkansas, Kansas, Missouri, and Oklahoma; and the Edwards-Trinity aquifer system in west-central Texas. These aquifers, and the springs that discharge from them, serve as major water-supply sources and form unique ecological habitats. Competition for the water resources of karst aquifers is common, and urban development and the lack of attenuation of contaminants in karst areas due to dissolution features that form direct pathways into karst aquifers can impact the ecosystem and water quality associated with these aquifers. The concept for developing a platform for interaction among scientists within the U.S. Geological Survey (USGS) working on karst-related studies evolved from the November 1999 National Groundwater Meeting of the USGS. As a result, the Karst Interest Group (KIG) was formed in 2000. The KIG is a loose-knit, grass-roots organization of USGS and non-USGS scientists and researchers devoted to fostering better communication among scientists working on, or interested in, karst science. The primary mission of the KIG is to encourage and support interdisciplinary collaboration and technology transfer among scientists working in karst areas. Additionally, the KIG encourages collaborative studies between the different mission areas of the USGS as well as with other federal and state agencies, and with researchers from academia and institutes. To accomplish its mission, the KIG has organized a series of workshops that have been held near nationally important karst areas. To date (2017) seven KIG workshops, including the workshop documented in this report, have been held. The workshops typically include oral and poster sessions on selected karst-related topics and research, as well as field trips to local karst areas. To increase non-USGS participation an effort was made for the workshops to be held at a university or institute beginning with the fourth workshop. Proceedings of the workshops are published by the USGS and are available online at the USGS publications warehouse https://pubs.er.usgs.gov/ by using the search term “karst interest group.” The first KIG workshop was held in St. Petersburg, Florida, in 2001, in the vicinity of the large springs and other karst features of the Floridan aquifer system. The second KIG workshop was held in 2002, in Shepherdstown, West Virginia, in proximity to the carbonate aquifers of the northern Shenandoah Valley, and highlighted an invited presentation on karst literature by the late Barry F. Beck of P.E. LaMoreaux and Associates. The third KIG workshop was held in 2005, in Rapid City, South Dakota, near evaporite karst features in limestones of the Madison Group in the Black Hills of South Dakota. The Rapid City KIG workshop included field trips to Wind Cave National Park and Jewel Cave National Monument, and featured a presentation by Thomas Casadevall, then USGS Central Region Director, on the status of Earth science at the USGS. The fourth KIG workshop in 2008 was hosted by the Hoffman Environmental Research Institute and Center for Cave and Karst Studies at Western Kentucky University in Bowling Green, Kentucky, near Mammoth Cave National Park and karst features of the Chester Upland and Pennyroyal Plateau. The workshop featured a late-night field trip into Mammoth Cave led by Rickard Toomey and Rick Olsen, National Park Service. The fifth KIG workshop in 2011 was a joint meeting of the USGS KIG and University of Arkansas HydroDays, hosted by the Department of Geosciences at the University of Arkansas in Fayetteville. The workshop featured an outstanding field trip to the unique karst terrain along the Buffalo National River in the southern Ozarks, and a keynote presentation on paleokarst in the United States was delivered by Art and Peggy Palmer. The sixth KIG workshop was hosted by the National Cave and Karst Research Institute (NCKRI) in 2014, in Carlsbad, New Mexico. George Veni, Director of the NCKRI, served as a co-chair of the workshop with Eve Kuniansky of the USGS. The workshop featured speaker Dr. Penelope Boston, Director of Cave and Karst Studies at New Mexico Tech, Socorro, and Academic Director at the NCKRI, who addressed the future of karst research. The field trip on evaporite karst of the lower Pecos Valley was led by Lewis Land (NCKRI karst hydrologist), and the field trip on the geology of Carlsbad Caverns National Park was led by George Veni. This current seventh KIG workshop is being held in San Antonio at the University of Texas at San Antonio (UTSA). This 2017 workshop is being hosted by the Department of Geological Sciences’ Student Geological Society (SGS), and student chapters of the American Association of Petroleum Geologists (AAPG) and Association of Engineering Geologists (AEG), with support by the UTSA Department of Geological Sciences and Center for Water Research. The UTSA student chapter presidents, Jose Silvestre (SGS), John Cooper (AAPG), and Tyler Mead (AEG) serve as co-chairs of the 2017 workshop with Eve Kuniansky of the USGS. The technical session committee is chaired by Eve Kuniansky, USGS, and includes Michael Bradley, Tom Byl, Rebecca Lambert, John Lane, and James Kaufmann, all USGS, and Patrick Tucci, retired USGS. The logistics committee includes Amy Clark, Yongli Gao, and Lance Lambert (Department Chair), UTSA Department of Geological Sciences; and Ryan Banta and Allan Clark, USGS, San Antonio, Texas. The field trip committee is chaired by Allan Clark and includes Amy Clark, Yongli Gao, and Keith Muehlestein, UTSA; Marcus Gary, Edwards Aquifer Authority and University of Texas at Austin; Ron Green, Southwest Research Institute; Geary Schindel, Edwards Aquifer Authority; and George Veni, NCKRI. Additionally, two organizations have assisted the UTSA student chapters in hosting the meeting by donating funds to the chapters: the Edwards Aquifer Authority, San Antonio, Texas, and the Barton Springs Edwards Aquifer Authority, Austin, Texas. Additionally, Yongli Gao, Center for Water Research and Department of Geological Sciences, UTSA, helped develop sessions on cave and karst research in China for this workshop. These proceedings could not have been accomplished without the assistance of Lawrence E. Spangler as co-editor who not only has subject matter expertise, but also serves as an editor with the USGS Science Publishing Network. We sincerely hope that this workshop continues to promote future collaboration among scientists of varied and diverse backgrounds, and improves our understanding of karst aquifer systems in the United States and its territories. The extended abstracts of USGS authors were peer reviewed and approved for publication by the USGS. Articles submitted by university researchers and other federal and state agencies did not go through the formal USGS peer review and approval process, and therefore may not adhere to USGS editorial standards or stratigraphic nomenclature. However, all articles had a minimum of two peer reviews and were edited for consistency of appearance in the proceedings. The use of trade, firm or product names is for descriptive purposes only and does not imply endorsement by the U.S. Government. The USGS Water Availability and Use Science Program funded the publication costs of the proceedings.

Scientific Investigations Report↗

Revised groundwater-flow model of the glacial aquifer system north of Aberdeen, South Dakota, through water year 2015

The city of Aberdeen, in northeastern South Dakota, requires an expanded and sustainable supply of water to meet current and future demands. Conceptual and numerical models of the glacial aquifer system in the area north of Aberdeen were developed by the U.S. Geological Survey in cooperation with the City of Aberdeen in 2012. The U.S. Geological Survey, in cooperation with the City of Aberdeen, completed a study to revise the original numerical groundwater-flow model using data through water year (WY) 2015 to aid the City of Aberdeen in their development of plans and strategies for a sustainable water supply and to increase understanding of the glacial aquifer system and groundwater-flow system near Aberdeen. The original model was revised to improve the fit between model-simulated values and observed (measured or estimated) data, provide greater insight into surface-water interactions, and improve the usefulness of the model for water-supply planning. The revised groundwater-flow model (hereafter referred to as the “revised model”) presented in this report supersedes the original model. The purpose of this report is to describe a revised groundwater-flow model including data collection, model calibration, and model results for the glacial aquifer system including the Elm, Middle James, and Deep James aquifers north of Aberdeen, South Dakota, using updated hydrologic data through WY 2015. The original numerical model was revised in several ways. The model was modified by adding four new layers, which included a surficial layer, two intervening confining layers, and a shale bedrock layer. The revised model provides an improved understanding of the groundwater-flow system in comparison to the original model. The principal aquifers of the model area include portions of the Elm, Middle James, and Deep James aquifers. The lithologic information used to define and describe the aquifers in the model area was unaltered; however, aquifer properties and boundary conditions were reviewed and updated using geological information reported by the South Dakota Department of Environmental and Natural Resources and information obtained from geophysical investigations for this study. The horizontal extent of the Elm, Middle James, and Deep James aquifers was unaltered from the original model. The thickness of the Deep James aquifer was modified based on interpretations from the geophysical investigations. In general, groundwater in the Elm aquifer flowed from northwest to southeast and locally towards rivers and streams. Similarly, in the Middle James and Deep James aquifers, groundwater also typically flowed southeast. The revisions made to the original model include use of the following MODFLOW stress packages: Recharge, Evapotranspiration, Time-Variant Specified Head, Wells, Drains, and Stream Flow Routing, all of which were updated from the original model except for the Stream Flow Routing Package, which replaced the River Package used in the original model. Model calibration is the process of estimating model parameters to minimize the differences, or residuals, between observed data and simulated values; therefore, Parameter ESTimation (PEST) software was used to optimize model input parameters by matching model-simulated values to observed data. Calibration parameters included horizontal hydraulic conductivity, vertical hydraulic conductivity, specific yield, specific storage, and vertical streambed conductance for stream and drain cells. Multipliers were used to calibrate the recharge and evapotranspiration stresses. Evapotranspiration extinction depth also was adjusted during model calibration. Comparisons to the original model are described to highlight the changes made in the revised model. In general, the revised model adequately simulates the natural system and compares favorably with observed hydrologic data. Simulated water levels were evaluated by comparing them to single water-level observations at selected well locations. The selected wells were the same wells used in the original model. The coefficient of determination value between simulated and observed water levels for the revised model was 0.89 and included simulated and observed values from October 1, 1974 (WY 1975), through September 30, 2015 (WY 2015). The coefficient of determination value for the original model was 0.94 and included simulated and observed values from October 1, 1974, through September 30, 2009. The difference may indicate that the original model could have been overfit to hydraulic head observations because base flow was not simulated. The additional data used in the revised model included some climatically wetter, more extreme periods, such as 2011, in which annual precipitation was 30.9 inches. Average annual precipitation for the original model timeframe, which included data from WYs 1975–2009, was 20.26 inches. Additional precipitation data for WYs 2010–15, included in the revised model timeframe, resulted in an average annual precipitation for WYs 1975–2015 in the model area of 20.6 inches. The larger variability in climate data coupled with the additional water-level data could explain the lower coefficient of determination for water levels in the revised model. The revised model was used to calculate various groundwater-budget components for steady-state and transient conditions for WYs 1975–2015. The time-variant specified-head cells in the revised model had the largest change when compared to the original steady-state model for inflows and outflows. Comparing the transient budget components between the original and the revised models indicated that inflow from recharge and time-variant specified-head cells had the greatest effect on groundwater inflows, and outflow from storage had the greatest effect on groundwater outflows. The simulated potentiometric contours from the revised model were compared with (1) the observed (interpreted) potentiometric surface (layer 2) and the hydraulic head values (layers 4 and 6) and (2) the simulated contours from the original model. The simulated hydraulic gradients and general direction of groundwater flow in the Elm aquifer in the revised model generally matched the observed potentiometric contours, the simulated potentiometric contours from the original model, and general flow directions interpreted to be perpendicular to the contours. Minor discrepancies between simulated potentiometric contours from the revised model and the observed potentiometric contours may be due to the lack of observed data in the model area. The revised model was designed to reduce the limitations of the original model. The revisions were validated by comparing the results of the original model with the revised model. A primary benefit of the revised model is the inclusion of the surficial deposits and the confining units as explicit layers in the model. The addition of the surficial layer was beneficial for three primary reasons: (1) more accurate representation of recharge from precipitation, (2) more accurate representation of groundwater evapotranspiration, and (3) more accurate representation of groundwater and surface-water interactions. The groundwater model is a numeric approximation of a complex physical hydrologic system, and the revised model data were interpolated in regions with sparse data. Additionally, model discretization included averaged and interpolated values for water use, withdrawal rates, and hydraulic conductivity. The revised model provides a useful estimate for hydraulic gradients, groundwater-flow directions, and aquifer response to groundwater withdrawals.

South Dakota↗

Monitoring Avian Productivity and Survivorship (MAPS) 6-year summary, Naval Outlying Landing Field, Imperial Beach, southwestern San Diego County, California, 2014–20

Executive Summary From 2014 to 2020, a Monitoring Avian Productivity and Survivorship (MAPS) banding station (station) was operated at the Naval Outlying Landing Field (NOLF), Imperial Beach, in southwestern San Diego County, California. The station was established as part of a long-term monitoring program of Neotropical migratory bird populations on NOLF and helps Naval Base Coronado (NOLF is a component) meet the goals and objectives of the Department of Defense Partners in Flight program and the Birds and Migratory Birds Management Strategies of the Naval Base Coronado Integrated Natural Resources Management Plan. The station was established in 2009 and has been in operation during the spring and summer since 2009 except for 2016 when it was not funded. The station was operated by AMEC Earth and Environmental, Inc., from 2009 to 2011, by the U.S. Geological Survey from 2012 to 2015, the San Diego Natural History Museum in 2017, and the U.S. Geological Survey again from 2018 to 2023. This report synthesizes results from 2014 to 2020. A prior report presents summaries and analyses from 2009 to 2013. The banding station at NOLF was operated according to the standard MAPS protocol with some exceptions. Ten mist nets used to capture birds were erected in fixed locations that remained consistent between and within years, with few minor relocations. Nets were open for 6 hours per day, once every 10 days (a netting period) for 13 netting periods starting April 1 each year. Occasionally, poor weather conditions (for example, rain, wind, or excessive heat) prevented net operation or forced nets to be closed early (or, rarely, late). Nets were checked periodically throughout the day and birds were removed, processed (leg bands affixed, measurements recorded), and released. From 2014 to 2020, we had 3,543 captures (including initial captures and recaptures) of a maximum of 3,264 year-unique captures (543±143 year-unique captures [the total number of individual birds captured for the first time each year]). The count of year-unique captures included 2,702 newly banded birds, 258 individuals that were recaptured from previous years, and 304 birds that were released unbanded (218 hummingbirds and 86 other birds that were intentionally released unbanded [game birds, and so forth] or escaped before banding). Individuals of 68 species were captured, 39 of which breed at or in the immediate vicinity of the MAPS banding station. Bird capture rate averaged 43±30 captures per day (corrected to account for variation in effort) for all years (range 7–163 effort-corrected captures per day) and species richness per year averaged 43±4. Bushtit ( Psaltriparus minimus ) was the most abundant species captured, followed by Orange-crowned Warbler ( Leiothlypis celata ), Wilson’s Warbler ( Cardellina pusilla ), House Finch ( Haemorhous mexicanus ), Song Sparrow ( Melospiza melodia ), and Common Yellowthroat ( Geothlypis trichas ). The mean adult sex ratio of all species combined across all years was 54:46 male:female. Adults averaged 73±12 percent of known age captures per year (range 59–94 percent), and juveniles averaged 27±12 percent (range 6–41 percent). Nineteen sensitive species were detected at NOLF (12 captured and 7 observed only). During 2014–20, we captured one State and federally endangered species, Least Bell’s Vireo ( Vireo bellii pusillus ); one federally threatened species, California Gnatcatcher ( Polioptila californica ); one State endangered species, Willow Flycatcher ( Empidonax traillii ); and two State species of concern, Yellow-breasted Chat ( Icteria virens ) and Yellow Warbler ( Setophaga petechia ). One additional State species of concern, Northern Harrier ( Circus hudsonius ), was observed at the MAPS banding station but not captured. Peregrine Falcon ( Falco peregrinus ) and White-tailed Kite ( Elanus leucurus ), California State fully protected species, also were observed at the MAPS banding station. Seven federal bird species of conservation concern—Calliope Hummingbird ( Selasphorus calliope ), Rufous Hummingbird ( Selasphorus rufus ), Allen’s Hummingbird ( Selasphorus sasin ), Nuttall’s Woodpecker ( Dryobates nuttallii ), Wrentit ( Chamaea fasciata ), California Thrasher ( Toxostoma redivivum ), and Lawrence’s Goldfinch ( Spinus lawrencei )—also were captured, and four additional federal bird species of conservation concern—Willet ( Tringa semipalmata ), Western Gull ( Larus occidentalis ), California Gull ( Larus californicus ), and Bullock’s Oriole ( Icterus bullockii )—were observed but not captured. Local population trends varied among species and years. From 2012 to 2019, year-round residents Bushtit, Song Sparrow, and Common Yellowthroat significantly decreased, whereas the migrant Least Bell’s Vireo increased. The total number of captures for all species except Least Bell’s Vireo was lowest in 2017, corresponding to the habitat damage caused by Kuroshio shot hole borer beetle ( Euwallacea kuroshio ) in the Tijuana River Valley. Annual productivity and annual adult survival were calculated for seven breeding species based on criteria used by the Institute for Bird Populations (Least Bell’s Vireo, Bushtit, Wrentit, House Wren [ Troglodytes aedon ], Song Sparrow, Orange-crowned Warbler, and Common Yellowthroat). Productivity was highest for most species in 2010 and 2019, years with high precipitation, and lowest in 2014 and 2018, years with low precipitation. Song Sparrow demonstrated the highest productivity among species and Least Bell’s Vireo had the lowest productivity. Annual adult survival was generally high from 2011 to 2012 and from 2018 to 2019. Bushtit had higher annual survival with lower late winter precipitation. Either temperature or precipitation was associated with productivity for all species except Wrentit, and with survival for all species except Least Bell’s Vireo and Common Yellowthroat. For most species, productivity was positively associated with precipitation, and both productivity and survival were negatively associated with temperature. Other studies have found that higher temperatures led to increased predation by snakes and birds and also increased vector-borne disease transmission, such as West Nile virus. Predicted regional increases in temperature over the next 30 years will likely affect the demographics of these species. The Song Sparrow population increased with higher breeding productivity during the previous year, and the Bushtit population increased with higher annual survival and higher productivity during the previous year. Aside from a possible positive association between survivorship and Common Yellowthroat population growth, productivity and survival rates did not appear to influence population change for other focal species.

California↗

The effects of artificial recharge on groundwater levels and water quality in the west hydrogeologic unit of the Warren subbasin, San Bernardino County, California

Between the late 1940s and 1994, groundwater levels in the Warren subbasin, California, declined by as much as 300 feet because pumping exceeded sparse natural recharge. In response, the local water district, Hi-Desert Water District, implemented an artificial-recharge program in early 1995 using imported water from the California State Water Project. Subsequently, the water table rose by as much as 250 feet; however, a study done by the U.S. Geological Survey found that the rising water table entrained high-nitrate septic effluent, which caused nitrate (as nitrogen) concentrations in some wells to increase to more than the U.S. Environmental Protection Agency maximum contaminant level of 10 milligrams per liter.. A new artificial-recharge site (site 3) was constructed in 2006 and this study, which started in 2004, was done to address concerns about the possible migration of nitrates in the unsaturated zone. The objectives of this study were to: (1) characterize the hydraulic, chemical, and microbiological properties of the unsaturated zone; (2) monitor changes in water levels and water quality in response to the artificial-recharge program at site 3; (3) determine if nitrates from septic effluent infiltrated through the unsaturated zone to the water table; (4) determine the potential for nitrates within the unsaturated zone to mobilize and contaminate the groundwater as the water table rises in response to artificial recharge; and (5) determine the presence and amount of dissolved organic carbon because of its potential to react with disinfection byproducts during the treatment of water for public use. Two monitoring sites were installed and instrumented with heat-dissipation probes, advanced tensiometers, suction-cup lysimeters, and wells so that the arrival and effects of recharging water from the State Water Project through the 250 to 425 foot-thick unsaturated zone and groundwater system could be closely observed. Monitoring site YVUZ-1 was located between two recharge ponds in the middle of site 3, and YVUZ-2 was located approximately 1,200 feet down-gradient and to the southeast in an area where septic systems have been in use since about 1960. Site YVUZ-3 only went to a depth of 42 feet and was used to sample the upper part of the unsaturated zone near a golf course. Prior to the start of artificial recharge at site 3, nitrate concentrations reported as nitrogen from the soil leachate below YVUZ-1 did not exceed 1.58 milligrams per kilogram. Nitrate-reducing bacteria concentrations of 4,300 most probable number were found at about 220 feet below land surface and at the top of the water table at YVUZ-1. Nitrate concentrations at YVUZ-2 reached a maximum concentration of about 25 milligrams per kilogram between about 100 and 121 feet below land surface; concentrations of nitrate-reducing or denitrifying bacteria were as high as 21,000 most probable number at 36 feet below land surface but did not exceed 40 most probable number below about 150 feet below land surface. Between June 2006 and September 2009, more than 9,800 acre feet of water from the State Water Project was released to site 3 ponds. The infiltration of the recharge water was predominantly vertical with limited lateral spreading to a depth of about 200 feet below land surface at YVUZ-1. Lateral spreading of the recharge water with depth was caused by geologic heterogeneities within the unsaturated zone, and resulted in varied arrival times of the recharge water to the instruments and slower rates of vertical movement with depth. No abrupt changes in soil moisture were observed at YVUZ-2, indicating that the recharge water had not reached that site by September 2009. Water levels from the monitoring wells at both sites and from five production wells nearby showed that the water table rose at a mean rate of about 0.08 feet per day between June 2006 and January 2009. The arrival of the water from the State Water Project caused relatively rapid changes in the stable-isotopic ratios from the lysimeters at YVUZ-1. The estimated average rate of infiltration of the recharge water through the unsaturated zone ranged from 3.7 to 25 feet per day. The recharge water arrived at the monitoring well below the recharge ponds between August 2007 and March 2008; the rate of vertical movement to the monitoring well was between 0.6 and 0.9 feet per day. By September 2008, a production well located 375 feet west of site 3 was producing almost 100 percent infiltrated recharge water. By contrast, the stable-isotope data from the lysimeters at YVUZ-2 showed that the recharge water had not reached this site by September 2009, but that septic effluent in the unsaturated zone likely had mixed with the native pore water to at least 154 feet below land surface. Assuming vertical infiltration, the minimum rate of infiltration of septic effluent since 1960 was about 3 feet per year. The isotopic data from the lysimeters at YVUZ-3 indicated two different sources of water to the upper 43 feet–irrigation-return flow and precipitation. Nitrate concentrations of the water from the State Water Project did not exceed 1 milligram per liter. Prior to artificial recharge, nitrate concentrations of the pore water at YVUZ-1 ranged between 6 to 18.2 milligrams per liter. After the arrival of the recharge water, the nitrate concentrations from the lysimeters and well at YVUZ-1 decreased to less than 1 milligram per liter, with the exception of samples collected at 205.5 feet, which did not exceed 4.12 milligrams per liter. The decrease in nitrate concentrations after artificial recharge indicated that the rising water table did not result in an increase of nitrates below YVUZ-1. At YVUZ-2, nitrate concentrations ranged between 12 to 479 milligrams per liter. The highest nitrate concentrations were at 92 feet below land surface and were almost seven times that of samples collected from a nearby septic tank. Nitrate concentrations from the lysimeter at 273 feet below land surface increased from 6 to almost 58 milligrams per liter after it was saturated by the rising water table in December 2007. These increases could be the result of the mobilization of high-nitrate water from regional sources of septic effluent after saturation, or the result of high-nitrate water present at the top of the water table that may be diluted deeper in the aquifer. Nitrate concentrations in groundwater from five nearby production wells and from both monitoring wells were less than 5 milligrams per liter before artificial recharge started. Nitrate concentrations decreased to less than 3 milligrams per liter in three of the production wells and the monitoring well below the recharge ponds after artificial recharge. Dissolved organic carbon concentrations were measured in the recharge water and groundwater because of the potential for dissolved organic carbon to react with chlorine to form trihalomethanes during the water-treatment process. The dissolved organic carbon concentrations of the recharge water were 3.1 milligrams per liter or less, and dissolved organic carbon concentrations of the groundwater were less than 1 milligram per liter. Even though recharge water was present in some of the wells by September 2008, the concentrations of both dissolved organic carbon and trihalomethane formation potential in the groundwater did not increase. Interpretation of these data suggests that the dissolved organic carbon from the recharge water is altered or metabolized in the unsaturated zone, either by absorption to the grain particles in the soil or by microbiological processes.

California↗

Fluvial sediment dynamics in the Shoshone River and tributaries around Willwood Dam, Park County, Wyoming

Sedimentation affects many of the aging reservoirs in the United States. Dams and water diversions from rivers have been central elements of infrastructure supporting agricultural irrigation in the arid and semiarid regions of the Western United States for more than a century. The Willwood Irrigation District diversion dam (hereafter referred to as “Willwood Dam”) in Park County, Wyoming, is approximately 12 miles northeast of Cody, Wyo.; has a structural height of 70 feet; and impounds the Shoshone River for diversion into the Willwood Canal. Willwood Dam is part of a larger irrigation scheme supported by water storage in the much larger Buffalo Bill Dam, which is approximately 20 miles upstream. In October 2016, renovation construction activities at Willwood Dam and the Willwood Canal caused an unplanned evacuation of nearly 96,000 cubic yards of fine sediment. The fine sediment release in 2016 raised concerns that ongoing sediment management at Willwood Dam could impose limits on the long-term health of the aquatic ecosystem and fish populations. The U.S. Geological Survey, in cooperation with Wyoming Department of Environmental Quality and Willwood Work Groups 2 and 3, initiated an investigation of the dynamics of sediment transport in the Shoshone River and selected tributaries between Buffalo Bill Dam and Willwood Dam. The goal of the study was to quantify sediment transport into and out of Willwood Dam on an annual, seasonal, and event basis to better understand the relative quantities of sediment coming from natural sources and human activities on the landscape. The study ran from March 2019 through October 2021 and used observations of streamflow, turbidity, and acoustic backscatter collected at streamgages upstream and downstream from Willwood Dam to quantify suspended-sediment loads into and out of the dam during irrigation and fallow seasons, precipitation-runoff events, and deliberate sediment releases. Each tributary’s relative contribution to the sediment load upstream from Willwood Dam was examined using discrete measurements of suspended-sediment concentration and bedload during irrigation and fallow seasons, precipitation events, and stable conditions. Analysis of daily precipitation and temperature data indicated that conditions in the study area during the 2019 agricultural year were wetter and colder than period of record normal, and drier and near normal temperatures for the 2020 and 2021 agricultural years. Not all sediment load records between 2019 and 2021 are complete because of rejected observations (outliers), instrument failures or fouling, and instrument removal for calibrations. Statistical modeling of suspended-sediment concentration using paired values of turbidity and acoustic backscatter produced four models that, after refinement, had coefficients of determination indicating that more than 84 percent of the variance was explained by either turbidity or acoustic backscatter. A system of rules was developed to select the model predictions based on the seasonal operations of Willwood Dam, assumptions about the grain sizes mobilized during these operations, and assumed accuracy of the models at the downstream streamgage (Shoshone River below Willwood Dam, near Ralston, Wyo. [streamgage 06284010]) under different operational conditions. The sediment budget between upstream and downstream estimates of loads was interpreted using the mean predicted values bound by their respective model prediction intervals. When mean predicted loads of one streamgage were contained in the prediction intervals of the other streamgage, and vice-versa, difference in the sediment budget were interpreted as “indeterminate.” Modeled sediment load balances demonstrated the depositional and erosional behaviors expected from the conceptual model of dam operations whereby sediment tends to accumulate during irrigation seasons when the dam is spilling over the top, and sediment tends to evacuate during the fallow seasons when it is flowing through the sluice gates at the base of the dam. The sediment load calculations using the rules-based model criteria indicated that between 14,200 and 380,000 tons of suspended sediment moved through the Shoshone River around Willwood Dam during the irrigation seasons of 2019, 2020, and 2021; 380,000 tons of suspended sediment were transported during the cool, wet year of 2019, and 14,200 tons of suspended sediment were transported in 2020, which was relatively dry. During fallow seasons 2019, 2020, and 2021, which had fewer complete records, between 1,140 and 106,000 tons of suspended sediment was estimated to have moved through the river. For all seasons except fallow season 2022, the models estimated that more sediment was released from the dam than entered the dam, but the modeled mean loads at each streamgage were nearly always within the prediction intervals of each other, making the sediment balance indeterminant. Examination of suspended-sediment loads during irrigation seasons indicated that between 65 and 85 percent of fine sediment was transported during annual high flows and storm events, with the remainder transported during steady, lower streamflows. Examination of suspended loads during fallow seasons indicated that deliberate sediment releases through Willwood Dam accounted for between 39 and 67 percent of the total sediment moved during the fallow seasons. Deliberate sediment releases from Willwood Dam had estimated net exports of between 1,360 and 22,400 tons. Between August 2017 and July 2023, suspended-sediment concentration and bedload sediment samples were collected from 9 tributaries to the Shoshone River during 137 sampling events, including stable and precipitation-runoff conditions. During irrigation season precipitation events, the mean total sediment yields ranged from 0.33 to 9.51 tons per day per square mile; during fallow season precipitation events, the mean total yields ranged from 0.04 to 0.95 ton per day per square mile. The mean total sediment yield per unit area across all samples at each tributary site ranged from 0.26 to 3.08 tons per day per square mile. Bedload was a minor fraction of the total load, constituting a mean of 4 percent across all samples; 3 and 6 percent for events and nonevents, respectively, during irrigation season; and 3 and 1 percent for events and nonevents, respectively, during the fallow season. With the exception of one tributary, Dry Creek, these mean yield values were within the range of watershed-scale background sediment yield values estimated from reservoir surveys and previous suspended-sediment studies. Imagery from irrigation seasons 2012, 2015, 2017, 2019, and 2022 was used to determine the planimetric backwater extent of the pool area in the Shoshone River behind Willwood Dam to identify any changes in sediment storage. Active river channel widths in the Shoshone River upstream from Willwood Dam were all similar between years except 2015, which was determined to be statistically different from all other years. Bathymetric data taken in the pool behind Willwood Dam during three different surveys between November 2017 and April 2022 indicated no statistically significant differences in bed elevations between the years. Results from the planimetric and bathymetric survey data provide multiple lines of evidence indicating that sediment did not accumulate behind the dam within the error of the methods used. Examination of how precipitation affects sediment transport in the Shoshone River upstream from Willwood Dam indicated that accumulated rainfall from the natural runoff events captured during the study period varied from a trace to as much as 4.26 inches, with associated predicted suspended-sediment loads varying from 112 to 232,000 tons of suspended sediment. The behavior of the sediment loads relative to accumulated precipitation did not appear to change depending on irrigation or fallow season. A model of suspended-sediment concentrations relative to the 2-day accumulated precipitation indicated that suspended-sediment concentrations in the Shoshone River upstream from Willwood Dam increased exponentially for accumulations of 0.3 inch or more; such storms accounted for 10 percent or less of precipitation events observed during the 1981 to 2018 period of record. The gaps in records, precision of the instrumentation, and large variation in grain sizes in suspended-sediment mixtures downstream from the dam made closing the sediment budgets for most seasons unattainable. The biggest recent change in sediment storage measured using the planimetric area of deposits behind Willwood Dam took place between 2015 and 2017. The main event between these two measurements was the installation of new Willwood Canal gates in October 2016, which resulted in the large unplanned sediment release. Because the sediment budgets were nearly always indeterminate and the planimetric and bathymetric data indicated little change in the bed and bank material, it is likely that the change in sediment storage behind the dam during the study period was small relative to the precision of the statistical models and other uncertainties. This body of evidence suggests that, averaged during the 3-year study period, no major changes in storage took place, and that the current operations may be keeping storage at near-equilibrium. This condition could have been initiated because the middle sluice gate has now been operational since 2014, and the sediment release in October 2016 evacuated a large amount of legacy sediment from storage. Although the uncertainties are large, sluicing events allow for controlled releases of sediment that contributed to the near equilibrium conditions observed over an annual basis during this study.

Wyoming↗

Selenium in ecosystems within the mountaintop coal mining and valley-fill region of southern West Virginia-assessment and ecosystem-scale modeling

Coal and associated waste rock are among environmental selenium (Se) sources that have the potential to affect reproduction in fish and aquatic birds. Ecosystems of southern West Virginia that are affected by drainage from mountaintop coal mines and valleys filled with waste rock in the Coal, Gauley, and Lower Guyandotte watersheds were assessed during 2010 and 2011. Sampling data from earlier studies in these watersheds (for example, Upper Mud River Reservoir) and other mining-affected watersheds also are included to assess additional hydrologic settings and food webs for comparison. Basin schematics give a comprehensive view of sampled species and Se concentration data specific to location and date. Food-web diagrams document the progression of Se trophic transfer across suspended particulate material, invertebrates, and fish for each site to serve as the basis for developing an ecosystem-scale model to predict Se exposure within the hydrologic conditions and food webs of southern West Virginia. This approach integrates a site-specific predator’s dietary exposure pathway into modeling to ensure an adequate link to Se toxicity and, thus, to species vulnerability. Site-specific fish abundance and richness data in streams documented various species of chub, shiner, dace, darters, bass, minnow, sunfish, sucker, catfish, and central stoneroller (Campostoma anomalum), mottled sculpin (Cottus bairdii), and least brook lamprey (Lampetra aepyptera). However, Se assessment species for streams, and hence, model species for streams, were limited to creek chub (Semotilus atromaculatus) and central stoneroller. Both of these species of fish are generally considered to have a high tolerance for environmental stress based on traditional comparative fish community assessment, with creek chub being present at all sites. Aquatic insects (mayfly, caddisfly, stonefly, dobsonfly, chironomid) were the main invertebrates sampled in streams. Collection of suspended particulate material acted as an integrator of organic-rich, fine-grained biomass present in streams. The base-case food web modeled for streams was suspended particulate material to aquatic insect to creek chub, with comparative modeling of a direct particulate-to-stoneroller food web. Model species for a reservoir setting were based on an earlier study of bluegill sunfish (Lepomis macrochirus), green sunfish (Lepomis cyanellus), and largemouth bass (Micropterus salmoides). Several reservoir food webs were considered based on a variety of invertebrates (insect, snail, clam). For stream and reservoir settings, predicted Se concentrations in exposure scenarios showed a high degree of correlation (r 2 = 0.91 for invertebrates and 0.75 for fish) with field observations of Se concentrations when modeling was initiated from suspended-particulate-material Se concentrations and model transfer parameters defined previously in the literature were used. These strong correlations validate the derived site-specific model and establish sufficient confidence that the predictions from the developed model can be quantitatively applied to the ecosystems in southern West Virginia. An application of modeling used a metric describing the partitioning of Se between particulate material and dissolved phases (K d ) to allow determination of a dissolved Se concentration that would be necessary to attain a site-specific Se fish body burden. The operationally defined K d quantifies the complex process of transformation at the base of a food web on a site-specific basis. The magnitude of this metric is known to vary with such factors as Se speciation, particulate-material type, and hydrology. This application (1) ties dissolved Se concentrations to fish tissue concentrations; (2) allows consideration of different choices for intervening site-specific exposure steps that set Se bioaccumulation, partitioning, and bioavailability; and (3) generates implications for management decisions that define protection through different regulatory pathways and guidelines. The range of model outcomes accounts for critical sources of variability and establishes whether site and food-web characterization were adequate to represent the dynamics of the system with certainty. This is especially true in terms of particulate-material phases at the base of the food web and utilization of K d in different hydrologic settings. For streams, a range of field-derived K d ds were applied to food-web exposure scenarios within a framework of locational and hydrologic variables (area of stream basin; stream gradient and discharge) that may affect the magnitude of K d . Overlaying even a coarse temporal scale that acknowledges variability in stream dissolved Se and Se speciation, such as through seasonal derivation of K d , can substantially narrow model uncertainty. Modeling that constrains the place and time of greatest ecosystem Se sensitivity within a specified food web gives insight into Se risk and identifies controlling management alternatives within a watershed or stream basin. If there is a range of hydrologic settings, specificity is needed to establish a hierarchy of in-stream and off-stream habitats for a watershed approach that takes into account Se-enriched water moving through different K d and food web environments. If there is a range of predator vulnerabilities (measured as a combination of food-web Se biodynamics and response in Se toxicity tests) within the site-specific community of fish species to be protected, then choice of fish species is critical to protection because it determines the food web and, hence, the magnitude of biotransfer through which Se is modeled. Whether creek chub is representative of the vulnerability to Se of all fish species encountered within the study-site ecosystems will require additional species-specific data and analysis. A range of site-specific scenarios illustrated here set model outcomes, but the final quantitative evaluation of alternatives and their implications will be those generated through choices and guidance formulated by state and other agencies in their decisionmaking processes. Proposed additions and refinements to the ecosystem-scale site-specific approach developed here include consideration of: measurement of temporally matched pairs of dissolved and suspended-particulate-material Se concentrations across a broader range of stream sites to expand the stream K d database and to test the representativeness of a suspended-particulate-material sample within a stream; characterization of different phases of particulate material across seasons to better define the base of the food web and connect to invertebrate feeding; refinement of model assumptions concerning dietary preferences and composition for fish to develop additional trophic transfer factors (TTFs) (for example, calculation of TTFinvertebrate composite for mixed diets); expansion of modeling of fish species and their food webs to include Se-vulnerable species; temporal characterization of a predator’s life cycle and habitat use as additional model layers to integrate with Se biodynamics in streams; investigation of the effect of stream gradient on K d based on a finer scale than presented here in terms of such variables as residence time, watershed dilution, and physical habitat attributes (for example, amount of ponding versus run or riffle within a stream); and linkage to discharge through use of stream gaging to record variability and enable model organization within water-year types and discharge seasons. Investigating the presence and variability of prey and predator species in demographically open systems such as streams also is key to model outcomes given the overall environmental stressors (for example, general landscape change, food-web disruption, recolonization potential) imposed on the composition of biological communities in coal mining and valley-fill affected watersheds

West Virginia↗

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↗

Davis Pond freshwater prediversion biomonitoring study: freshwater fisheries and eagles

In January 2001, the construction of the Davis Pond freshwater diversion structure was completed by the U.S. Army Corps of Engineers. The diversion of freshwater from the Mississippi River is intended to mitigate saltwater intrusion from the Gulf of Mexico and to lessen the concomitant loss of wetland areas. In addition to the freshwater inflow, Barataria Bay basin would receive nutrients, increased flows of sediments, and water-borne and sediment-bound compounds. The purpose of this biomonitoring study was, therefore, to serve as a baseline for prediversion concentrations of selected contaminants in bald eagle ( Haliaeetus leucocephalus ) nestlings (hereafter referred to as eaglets), representative freshwater fish, and bivalves. Samples were collected from January through June 2001. Two similarly designed postdiversion studies, as described in the biological monitoring program, are planned. Active bald eagle nests targeted for sampling eaglet blood (n = 6) were generally located southwest and south of the diversion structure. The designated sites for aquatic animal sampling were at Lake Salvador, at Lake Cataouatche, at Bayou Couba, and along the Mississippi River. Aquatic animals representative of eagle prey were collected. Fish were from three different trophic levels and have varying feeding strategies and life histories. These included herbivorous striped mullet ( Mugil cephalus ), omnivorous blue catfish ( Ictalurus furcatus ), and carnivorous largemouth bass ( Micropterus salmoides ). Three individuals per species were collected at each of the four sampling sites. Freshwater Atlantic rangia clams ( Rangia cuneata ) were collected at the downstream marsh sites, and zebra mussels ( Dreissena spp.) were collected on the Mississippi River. The U.S. Geological Survey (USGS) Biomonitoring of Environmental Status and Trends (BEST) protocols served as guides for fish sampling and health assessments. Fish are useful for monitoring aquatic ecosystems because they accumulate pesticides and other contaminants. Biomarker data on individual fish, generated at the USGS National Wetlands Research Center (Lafayette, La.), included percent white blood cells in whole blood, spleen weight to body weight ratio, liver weight to body weight ratio, condition factor, splenic macrophage aggregates, and liver microsomal 7-ethoxyresorufin-o-deethylase (EROD) activity. Fish age was estimated by comparing total lengths with values from the same species in the Southeast United States as determined from the literature. Contaminant analyses were coordinated by the U.S. Fish and Wildlife Service (USFWS) Analytical Control Facility (Laurel, Md.), where residues of organochlorine (OC) pesticides, total polychlorinated biphenyls (PCBs), polycyclic aromatic hydrocarbons (PAHs), aliphatic hydrocarbons (AHs), and trace elements were determined. The organic contaminant data were generated at the Mississippi State University Chemical Lab (Mississippi State, Miss.), and the inorganic contaminant data were generated by the Texas A&M University Geochemical and Environmental Research Group (College Station, Tex.). Statistical tests were performed to assess relationships among contaminants, fish age, fish species, and collection sites. Trends in interspecific differences among fish in concentrations of contaminants were noted. Striped mullet (hereafter mullet) frequently displayed the highest chemical concentrations. Levels of contaminants were generally higher in samples obtained from the Mississippi River than in those collected from the diversion area and were higher in mussels and clams (hereafter bivalves) than in fish. Because the Mississippi River sampling site for mullet and largemouth bass was downriver of the structure and south of New Orleans and the catfish site was upriver, the downriver data may not be directly reflective of the results from the receiving waters at the Davis Pond structure. Compared to the Caernarvon freshwater prediversion study in 1990 that assessed possible influx of contaminants with the freshwater diversion, contaminant levels in fishes and bivalves in this study were generally lower, yet three nontoxic inorganic elements in Davis Pond fish samples exhibited ranges of concentrations that were more than two times higher than did those from Caernarvon. Levels in bivalves were different between diversions but about equal in the numbers of trace elements showing high levels per location. Contaminant values were compared to those listed in various literature and agency sources, both regional and national, including the National Contaminant Biomonitoring Program (NCBP), in which the 85th percentile and above represents what is considered to be an elevated contaminant concentration and cause for concern. Generally, bivalves were at the high end of their ranges for both organic and inorganic contaminants. In this study, OCs were detectable in 67 percent of fish from the Mississippi River site, ranging from 0.15 to 1.09 μg/g wet weight (ww) or fresh weight (fw), and in 11 percent of the fish from the marsh sites, ranging from 0.06 to 0.612 μg/g ww. Bivalves from the Mississippi River had OC levels of 0.096 μg/g ww, whereas none were detectable in bivalves at the marsh sites. In this study, p,p ’-dichlorodiphenyldichloroethylene ( p,p’ -DDE) (a biodegradation product of DDT [dichlorodiphenyl trichloroethane]) and total PCBs were the most frequently detected OCs and were primarily from the Mississippi River. For total OC content, using adjusted least squares means, some significant interactions were noted between fish species and sites. PAHs were detected in aquatic animals at all sites (range of 0.017–17.534 μg/g ww), as were AHs (range of 0.423–4.549 μg/g ww); the highest levels of PAHs and AHs were found in bivalves from the Mississippi River. When analysis of variance (α = 0.05) was performed with data from aquatic animals, there were only two significant relationships between PAHs, AHs, and OCs between species, site, and age or the interaction among these variables. There was an interaction between fish species and n-decane (an AH) in that mullet and largemouth bass had significantly higher levels than did catfish ( P = 0.0175). When general linear means were used to investigate associations of inorganic contaminants among fish species, site, and age or any interactions among these variables, no significant results were noted for arsenic, cadmium, lead, beryllium, boron, molybdenum, or nickel. The range of mercury in fish in this study was 0.04–0.14 μg/g ww (0.14– 0.48 μg/g dry weight [dw]), with the most elevated levels detected in predatory largemouth bass at the sampling point farthest downstream from the structure and within the marsh area. Mercury was positively correlated with fish age ( P = 0.0152), where levels were estimated to increase 0.0253 parts per million (ppm) dw per year. In the Mississippi River, catfish showed significantly higher levels of mercury than did mullet or largemouth bass ( P = 0.00167). Among fish species, mullet displayed the highest levels in fish of aluminum, barium, manganese, and iron, all considered to have low toxicity in hydrologic systems. An interaction between fish and site was seen with aluminum ( P = 0.0031), where concentrations in mullet were significantly higher in the Mississippi River than at the other sites, as was also seen with barium ( P = 0.0009), chromium ( P = <0.0001), manganese ( P = 0.0004), strontium ( P = 0.0074), vanadium ( P = 0.0156), and zinc ( P = 0.0059). For iron ( P = 0.0.0001), mullet and largemouth bass at both the Mississippi River and Lake Salvador showed higher levels than did catfish, and these two species showed higher levels at two of the four sites. An interaction between fish and site was also seen with chromium ( P = <0.0001) in that concentrations in mullet were significantly higher in the Mississippi River than at the other sites, as was also seen with strontium ( P = 0.0074), vanadium ( P = 0.0156), and zinc ( P = 0.0059), metals for which deleterious effects have been demonstrated in other ecosystems. The NCBP program lists the 85th percentile for zinc at 34.2 μg/g fw (117.9 μg/g dw). In the Davis Pond prediversion biomonitoring study (hereafter the current study), one fish (MUL31RIVER, fish ID 8) showed values higher than that (125.4 μg/g dw or 37.54 μg/g ww), and the Mississippi River bivalve sample (MUSSRIVER) had a value of 140 μg/g dw (41.2 μg/g ww). In the current study, approximately 86 percent of the fish had measurable selenium levels, yet none reached the 85th percentile. The 85th percentile for selenium from the NCBP was 0.73 μg/g ww. Significantly higher levels of selenium were seen in mullet than in largemouth bass and catfish ( P = 0.0023). The NCBP 85th percentile for lead is 0.22 μg/g ww (0.76 μg/g dw). In the current study, the range of concentrations of lead was as much as 18.3 ppm dw (MUL31RIVER, fish ID 8), with the three most elevated values (range of 3.46–5.31 μg/g ww) coming from mullet from the Mississippi River. Biomarker data are measurable and directly reflect the condition of the animal, and measuring more than one biomarker in an individual increases confidence in health assessments. In the current study, biomarkers included macrophage aggregates (MAs), liver (hepatosomatic index [HSI]) and spleen (splenosomatic index [SSI]) weight to body weight ratios, percent white blood cells (WBCs) in whole blood, and condition factor. Few significant differences were noted with any of the biomarkers between sites, and there were no relationships between species and sites. For improved use of biomarker assessments, an increase in fish sample size would be useful for postdiversion sampling, as would comparisons of fish of the same sex and reproductive condition. During the current study, success for eagle nests in the diversion area and reference sites was similar as determined by numbers of nestlings fledged. When temperatures were below average during winter 2000, nests in both regions similarly failed. At each nest, the primary evidence of food items was small mammals. Eaglets (n = 6) generally appeared healthy, and whole blood concentrations of organic contaminants exceeded detection limits with three incidences of p,p’ -DDE (0.002–0.006 μg/L ww) and one incidence of oxychlordane (0.002 μg/L ww). The levels of p,p’ -DDE were well below those that have been inversely correlated with productivity and success rates of nesting bald eagles on a regional scale. The low values found in the whole blood samples for OC pesticides and PCBs were even lower when corrected for plasma volume. Aluminum values were 3.66 and 5.75 μg/L in two samples, zinc ranged from 5.21 to 6.77 μg/L ww in six samples, and silicon ranged from 1.7 to 4.6 μg/L in four samples. Selenium was detectable in each bird with the range at 0.332–0.566 μg/L ww, and strontium ranged from 0.0581 to 0.0975 μg/L ww. Mercury was detectable in blood samples from each bird and ranged from 0.0254 to 0.0845 μg/L ww, whereas lead was detectable in four samples and ranged from 0.0042 to 0.0136 μg/L ww. Although no detectable levels of total PCBs were found (also correlated with decreased reproductive productivity), 70 percent of the aquatic animals from the Mississippi River contained total PCBs (range 0.13–0.79 μg/L), whereas only about 7 percent of the aquatic animals sampled from the marsh area contained PCBs. Suggestions for postdiversion sampling include lowering the analytical detection limit for some metals, sampling aquatic animals over the course of a single season, obtaining a higher sample number of mature fish of one species (for example, blue catfish) within a range of total lengths for biomarker analyses, obtaining otoliths for estimating fish ages, assessing dioxins in eaglet blood, examining triazines in water, and obtaining all Mississippi River fish samples as close to the Davis Pond structure intake as possible. Because contaminants found in blood of eaglets reflect their prey species and because of the contaminant levels found in fish in the current study, eaglets may not be consuming primarily these species; therefore, obtaining juvenile nutria ( Myocastor coypus ) or turtle species for contaminant analyses might be considered, as well as collecting greater blood volume and using plasma to measure OCs and PCBs. Data obtained postdiversion will be compared with prediversion data to monitor changes.

Louisiana↗