Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Water Surface Techniques”

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

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

At least 847 records · Page 47Linked to original sources

Volcanic gas monitoring

Introduction As magma rises through the crust, decreasing pressure conditions allow volatiles to exsolve from the magma. These volatiles then migrate upward through the crust, where they can be stored at shallower levels or escape to the atmosphere. Rising magma also heats rock masses beneath volcanic centers, causing water in shallow aquifers and hydrothermal systems to boil and release additional gases and steam (see chapter F , this volume; Ingebritsen and Hurwitz, 2024). The chemistry and quantity of gases that reach the surface during periods of quiescence or volcanic unrest can reveal that gas-rich magma is ascending, crystallizing, or alternatively stalling, with important implications for volcanic hazard (for example, Sutton and others, 1992; Aiuppa and others, 2007, 2021; Werner and others, 2009, 2011, 2012; Moretti and others, 2013; de Moor and others, 2016; Lewicki and others, 2019; Edmonds and others, 2022; Kern and others, 2022; Kunrat and others, 2022). Most volcanoes in Alaska and the western United States are characterized by weak degassing, with one or more low-temperature fumaroles (typically near the local boiling temperature of water) and connect to a deeper and sometimes extensive hydrothermal system (for example, McGee and others, 2001; Symonds and others, 2003a, b). Hydrothermal systems will affect the chemistry of rising gases exsolved from deeper magma (Symonds and others, 2001), including sulfur dioxide (SO 2 ), hydrogen chloride (HCl), and water vapor (for example, Doukas and Gerlach, 1995; Gerlach and others, 1998, 2008; Symonds and others, 2001; Werner and others, 2013). As an example, depending on factors such as temperature, pressure, and oxidation state, rising SO 2 will react with groundwater to form hydrogen sulfide (H 2 S) gas, dissolved sulfate (SO 4 2− ), or elemental sulfur (Christenson, 2000; Symonds and others, 2001; Werner and others, 2008). The reaction and dissolution of SO 2 into shallow groundwater is commonly referred to as scrubbing, and can reduce the likelihood that ascending, degassing magma can be detected. Carbon dioxide, however, in addition to exsolving from magma early in the ascent process, is not easily removed by hydrothermal fluids (Lowenstern, 2001). As scrubbing and other processes take place, the SO 2 /H 2 S, CO 2 /SO 2 , and CO 2 /H 2 S ratios may change. High rates of SO 2 emission indicate that magma has moved to relatively shallow levels in the volcano and that the system has heated up enough to establish dry pathways from depth to the surface. Monitoring multiple gas species and the total output of those species is thereby useful for volcano monitoring during both periods of quiescence, to establish background degassing conditions, and during unrest, when gas geochemistry and emission rates can provide information on changing conditions, such as magma ascent. To provide context for multidisciplinary volcano forecasts, we focus on the following two key required capabilities: (1) characterizing baseline geochemistry and gas discharge from volcanoes and volcanic regions and (2) monitoring changes in gas geochemistry and discharge to inform forecasts of volcanic eruptions and their effects. Sufficient baseline data must be collected to identify and interpret anomalous degassing associated with volcanic unrest (for example, Sorey and others, 1998; Rouwet and others, 2014). Differences in volcano type, baseline degassing rates, local hydrology, and geography (for example, high versus low latitude) will result in a different baseline for each volcano. Volcanoes of any threat level that exhibit one or more degassing phenomena would ideally be monitored by techniques needed to establish baseline degassing data, with the sampling frequency of baseline data dictated by the threat level (table E1). Additional monitoring techniques become necessary during periods of unrest. In general, three of the most important techniques for gas monitoring are (1) direct sampling of fumarole, spring, and soil gases for laboratory geochemical measurements, (2) measurements of the chemical composition of the volcanic plume and emission rates of major gas species (for example, H 2 O, CO 2 , SO 2 , and H 2 S) by satellite, airborne, or ground-based techniques, and (3) measurements of diffuse emissions of CO 2 and other gases through soils. Various methods and instruments may be useful both for baseline studies and during unrest.

Scientific Investigations Report↗

Organic geochemistry applied to environmental assessments of Prince William Sound, Alaska, after the Exxon Valdez oil spill—a review

Organic geochemistry played a major role in the environmental assessments conducted following the Exxon Valdez oil spill, which occurred on March 24, 1989, and released about 258,000 bbls (41 million liters) of Alaska North Slope crude oil into Prince William Sound. Geochemical analyses of more than 15,000 sediment, tar, and biological samples and about 5000 water samples provide the largest database yet collected on oil-spill chemistry, and we review the results here. The marine environment of the Sound has a complex background of petrogenic, pyrogenic, and biogenic hydrocarbons from natural and anthropogenic sources. Geochemical evaluation of the fate and effects of the spilled oil required that this oil and its residues be distinguished from the background. A variety of molecular and isotopic techniques were employed to identify various hydrocarbon sources and to distinguish quantitatively among mixed sources in the samples. Although the specific criteria used to distinguish multiple sources in the region affected by the Exxon Valdez spill are not necessarily applicable to all spill situations, the principles that governed their selection are. Distributions of polycyclic aromatic hydrocarbons (PAH) and dibenzothiophenes distinguish Exxon Valdez oil and its weathered residues from background hydrocarbons in benthic sediments. Ratios of C 2 -dibenzothiophene C 2 -phenanthrene "> C 2 -dibenzothiopheneC 2 -phenanthrene and C 3 -dibenzothiophene C 3 -phenanthrene "> C 3 -dibenzothiopheneC 3 -phenanthrene were particularly useful. Carbon isotopes and terpane distributions distinguished Exxon Valdez residues found on shorelines from tars from other sources. Diesel and diesel soot were identified by the absence of alkylated chrysenes and a narrow distribution of n -alkanes, whereas pyrogenic products were distinguished by the dominance of 4- to 6-ring PAH over 2- to 3-ring PAH and by the dominance of non-alkylated over alkylated homologues of each PAH series. The presence of 18α(H)-oleanane in benthic sediments, coupled with its absence in Exxon Valdez oil and its residues, confirm another petrogenic source. Results of geochemical studies suggest that the petrogenic component in the background of benthic sediments is derived from oil seeps in the eastern Gulf of Alaska. In 1990 and 1991, Exxon Valdez residues, generally forming a small increment to the pre-spill background, were found to be only sporadically distributed in some shallow, near shore sediments adjacent to shorelines that had been heavily oiled in 1989. In 1994, occurrences of Exxon Valdez tars on shoreline surfaces were rare, although residues could be found buried in shoreline sediments at some isolated locations along the spill path where they were protected from wave action. Spilled oil residues collected 16 months after the spill were degraded, on average, by nearly 50%. Shoreline residues from sources other than the spill were also identified and are widespread throughout the Sound. These residues include (1) geochemically distinct tars and oils imported from California oil fields to Alaska for fuel and construction purposes prior to the discovery of the Cook Inlet and North Slope oil fields, (2) diesel and diesel soot, and (3) more highly refined products. Of the more than 2700 chemical analyses of biological samples of higher life forms (fish, birds, and mammals) about 150 (6%) indicate recognizable residues of Exxon Valdez oil, which were identified by their distribution of polycyclic aromatic hydrocarbons (PAH). Most of these samples (138) were collected in 1989 and most were associated with external surfaces or the gastrointestinal tract. Rarely do internal tissues or fluids contain recognizable fingerprints of spilled oil. This observation includes samples from marine mammals that were visibly oiled externally. Other hydrocarbon sources, including diesel and a non-petroleum artifact that occurs when concentrations of individual PAH are at or near their method detection limit, are also identified in biological samples.

Alaska↗

Federal standards and procedures for the National Watershed Boundary Dataset (WBD)

The Watershed Boundary Dataset (WBD) is a seamless, national hydrologic unit dataset developed using the guidelines and specifications outlined in this document. A hydrologic unit defines the areal extent of surface-water drainage to an outlet point on a dendritic stream network or to outlet points where the stream network is not dendritic. A hydrologic unit may represent all or only part of the total drainage area to the outlet point so that multiple hydrologic units may be required to define the entire drainage area for a given outlet. Hydrologic unit boundaries in the WBD are determined on the basis of topographic, hydrologic, and other relevant landscape characteristics without regard for administrative, political, or jurisdictional boundaries. The WBD seamlessly represents hydrologic units at six required and two optional hierarchical levels mapped at a minimum of 1:24,000-scale in the United States, except for in Hawaii, the Caribbean, and the Pacific Islands, which are at 1:25,000-scale, and in Alaska, where the data range from the minimum required 1:24,000-scale to 1:63,360-scale. Hydrologic units in the WBD provide a standardized base for water-resources organizations to locate, store, retrieve, and exchange hydrologic data; to index and inventory hydrologic data and information; to catalog water-data acquisition activities; and to use in a variety of other applications. The specifications and procedures established in this document are intended to provide guidelines and best practices to local, regional, and national partners for delineating and updating the hydrologic units of the WBD. Maintaining the WBD using consistent practices improves watershed management through efficient sharing of information and resources and by ensuring that digital geographic data can be used with other related geographic information system data. Edits and additions to the hydrologic units are reviewed against this specification prior to inclusion to the national WBD.

Techniques and Methods↗

Petrography, geochemistry, and depositional setting of the San Pedro and Santo Tomas coal zones: Anomalous algae-rich coals in the middle part of the Claiborne Group (Eocene) of Webb County, Texas

Two coal zones, the San Pedro and the overlying Santo Tomas, are present for nearly 35 km in outcrop, surface and underground mines, and shallow drill holes along the strike of the middle part of the Claiborne Group (Eocene) in Webb County, Texas. A sandstone-dominated interval of 25 to 35 m separates the two coal zones, which range up to 3 m in thickness. Each coal zone contains carbonaceous shales, thin (<0.75 m) impure coal beds, and thin (<0.85 m) but commercially significant nonbanded coal beds. The nonbanded coals are different from other Tertiary coals of the Gulf of Mexico Coastal Plain: unlike lignites that are typical of the older Wilcox Group (Paleocene-Eocene) and younger Jackson Group (Eocene), nonbanded coals of the Claiborne Group have high vitrinite-reflectance values (0.53 R max ) and high calorific yields (average 6670 kcal/kg or 12,000 Btu, dry basis). The coals are weakly agglomerating (free-swelling index is 1.5–2.0) and have an apparent rank of high-volatile bituminous. The coal-bearing portion of the middle Claiborne Group in the Rio Grande area represents a fining-upward transition from sandstone-dominated, marine-influenced, lower delta plain depositional environments to more inland, mudstone-rich, predominantly freshwater deltaic settings. Discontinuities within the San Pedro coal zone are attributed mainly to the influence of contemporaneous deposition of distributary mouth-bar sand bodies. The less variable nature of the Santo Tomas coal zone reflects its origin in the upper part of an interlobe basin that received only minor clastic influx. Petrographic attributes of the nonbanded coals indicate that they formed subaqueously in fresh to possibly brackish waters. A highly degraded groundmass composed of eugelinite is the main petrographic component (approximately 71%, mineral-matter-free basis). An enriched liptinite fraction (approximately 23%) probably accounts for unusually high calorific values. There is negligible inertinite. Petrographic study of polished blocks indicates that approximately 10 percent of the nonbanded coal from both coal zones is composed of green algae fructifications, which also occur in clastic rocks of the coal-bearing interval. Such algal material cannot be identified or quantified by conventional coal petrographic techniques that utilize particle pellets or by palynological analyses that include acid preparation.

International Journal of Coal Geology↗

Simulation of a long-term aquifer test conducted near the Rio Grande, Albuquerque, New Mexico

A long-term aquifer test was conducted near the Rio Grande in Albuquerque during January and February 1995 using 22 wells and piezometers at nine sites, with the City of Albuquerque Griegos 1 production well as the pumped well. Griegos 1 discharge averaged about 2,330 gallons per minute for 54.4 days. A three-dimensional finite-difference ground-water-flow model was used to estimate aquifer properties in the vicinity of the Griegos well field and the amount of infiltration induced into the aquifer system from the Rio Grande and riverside drains as a result of pumping during the test. The model was initially calibrated by trial-and-error adjustments of the aquifer properties. The model was recalibrated using a nonlinear least-squares regression technique. The aquifer system in the area includes the middle Tertiary to Quaternary Santa Fe Group and post-Santa Fe Group valley- and basin-fill deposits of the Albuquerque Basin. The Rio Grande and adjacent riverside drains are in hydraulic connection with the aquifer system. The hydraulic-conductivity values of the upper part of the Santa Fe Group resulting from the model calibrated by trial and error varied by zone in the model and ranged from 12 to 33 feet per day. The hydraulic conductivity of the inner-valley alluvium was 45 feet per day. The vertical to horizontal anisotropy ratio was 1:140. Specific storage was 4 x 10-6 per foot of aquifer thickness, and specific yield was 0.15 (dimensionless). The sum of squared errors between the observed and simulated drawdowns was 130 feet squared. Not all aquifer properties could be estimated using nonlinear regression because of model insensitivity to some aquifer properties at observation locations. Hydraulic conductivity of the inner-valley alluvium, middle part of the Santa Fe Group, and riverbed and riverside-drain bed and specific yield had low sensitivity values and therefore could not be estimated. Of the properties estimated, hydraulic conductivity of the upper part of the Santa Fe Group was estimated to be 12 feet per day, the vertical to horizontal anisotropy ratio was estimated to be 1:82, and specific storage was estimated to be 1.2 x 10-6 per foot of aquifer thickness. The overall sum of squared errors between the observed and simulated drawdowns was 87 feet squared, a significant improvement over the model calibrated by trial and error. At the end of aquifer-test pumping, induced infiltration from the Rio Grande and riverside drains was simulated to be 13 percent of the total amount of water pumped. The remainder was water removed from aquifer storage. After pumping stopped, induced infiltration continued to replenish aquifer storage. Simulations estimated that 5 years after pumping began (about 4.85 years after pumping stopped), 58 to 72 percent of the total amount of water pumped was replenished by induced infiltration from the Rio Grande surface-water system.

Water-Resources Investigations Report↗

Action plan for restoration of coral reef coastal protection services: Case study example and workbook

This report was prepared by the U.S. Environmental Protection Agency (USEPA), Office of Research and Development, as part of the Air, Climate and Energy (ACE) research program, with support from Tetra Tech, Inc., and in collaboration with the National Oceanic and Atmospheric Administration, the U.S. Geological Survey, and The Nature Conservancy. The ACE research program provides scientific information and tools to support USEPA’s commitment to clean air, clean water and sustainable natural resources, even as environmental conditions change. A key component of this is the development of sound science to support adaptation. Adaptation involves preparing for and adjusting to the effects of climate change and its interactions with other global and local stressors. Because these effects are diverse, interactive, and difficult to predict, adapting management of natural resources in this context can be very challenging. Coral reefs—which provide valued ecosystem services such as fisheries, coastal protection, and tourism—are threatened by the effects of increased sea surface temperatures, sea level rise, and intensifying storms. These large-scale stressors are interacting with local stressors such as pollution, overfishing, and recreational misuse to drive ongoing and accelerating declines in coral reef ecosystems. Thus, there is a rising urgency to design and implement climate change adaptation measures that will enable reef resilience in the face of these changes. This includes accounting for, and adjusting to, the combined effects of climate change and local stressors in coral reef protection and restoration efforts. The action plan, example case study, and workbook found in this report demonstrate a structured process for integrating climate-smart design considerations into restoration planning using A Manager’s Guide to Coral Reef Restoration Planning and Design. The focus is a hypothetical coral reef restoration project that has a goal of recovering nature-based coastal protection services using restoration interventions. The intent is to provide readers with a completed example of how to use the Guide workbook to inform a draft action plan, centering on the topic of coastal protection as a burgeoning area of interest in coral reef science and management communities. The information in this hypothetical case study is not intended for direct use; rather, it provides a starting point for more detailed planning that would occur in specific places. And while a full review of the current literature on reef restoration methods is outside the scope of this report, readers are encouraged to use the examples herein as well as in the Guide as a jumping-off point for exploring the rapidly growing body of information on methods, techniques, successes, failures, monitoring challenges and future directions of coral reef restoration in a changing world. The workbook, together with the action plan, can serve as a valuable record of the planning thought process as well as a living document for adaptive management, to be updated through time as improved information becomes available.

EPA Report↗

Hydrologic budget of the Beaverdam Creek basin, Maryland

A hydrologic budget is a statement accounting for the water gains and losses for selected periods in an area. Weekly measurements of precipitation streamflow, surface-water storage, ground-water stage, and soil resistivity were made during a 2year period, April 1, 1950, to March 28, 1952, in the Beaverdam Creek basin, Wicomico County, Md. The hydrologic measurements are summarized in two budgets, a total budget and a ground-water budget, and in supporting tables and graphs. The results of the investigation have some potentially significant applications because they describe a method for determining the annual replenishment of the water supply of a basin and the ways of water disposal under natural conditions. The information helps to determine the 'safe' yield of water in diversion from natural to artificial discharge. The drainage basin of Beaverdam Creek was selected because it appeared to have fewer hydrologic variables than are generally found. However, the methods may prove applicable in many places under a variety of conditions. The measurements are expressed in inches of water over the area of the basin. The equation of the hydrologic cycle is the budget balance: P= R+E+ASW+ delta SW + delta SM + delta GW where P is precipitation; R is runoff; ET is evapotranspiration; delta SW is change in surface-water storage; delta SM is change in soil moisture; and delta GW is change in ground-water storage. In this report 'change' is the final quantity minus the initial quantity and thus is synonymous with 'increase.' Further, ,delta GW= delta H .x Yg, in which delta H is the change in ground-water stage and Yg is the gravity yield, or the specific yield of the sediments as measured during the short periods of declining ground-water levels characteristic of the area. The complex sum of the revised equation P ? R - delta SW ? ET - delta SM, which is equal to delta H. x Yg, has been named the 'infiltration residual'; it is equivalent to ground-water recharge. Two unmeasured, but not entirely unknown, quantities, evapotranspiration, (ET) and gravity yield, (Yg), are included in the equation. They are derived statistically by a method of convergent approximations, one of the contributions of this investigation. On the basis of laboratory analysis, well-field tests, and general information on rates of drainage from saturated sediments, a gravity yield of 14 percent was assumed as a first approximation. The equation was then solved, by weeks, for evapotranspiration, ET. The evapotranspiration losses were plotted against the calendar week. Using the time of year as a control, a smooth curve was fitted to the evapotranspiration data, and modified values of ET were read from the curve. These were used to compute weekly values of the infiltration residual which were plotted against ground-water stage. The slope of the line of best fit gave a closer approximation of gravity yield, Yg. The process was repeated. The approximations converged, so that a fourth and final approximation resulted in a close grouping of all the points along a line whose slope indicated a Yg of 11.0 percent, and a slightly asymmetric bell-shaped curve of total evapotranspiration by weeks was obtained that is considered representative of this area. Check calculations of gravity yield were made during periods of low evapotranspiration and high infiltration, which substantiate the computed average of 11.0 percent. Refinements in the method of deriving the ground-water budget were introduced to supplement the techniques developed by Meinzer and Stearns in the study of the Pomperaug River basin in Connecticut in 1913 and 1916. The hydrologic equation for the ground-water cycle may be written Gr=D + delta H. x Yg + ETg, in which Gr is ground-water recharge (infiltration); D is ground-water drainage; delta H is the change in mean ground-water stage (final stage minus initial stage); Yg is gravity yield (taken as 11.0 percent in computations here); an

Water Supply Paper↗

An integrated hydrologic model to support the Central Platte Natural Resources District Groundwater Management Plan, central Nebraska

The groundwater and surface-water supply of the Central Platte Natural Resources District supports a large agricultural economy from the High Plains aquifer and Platte River, respectively. This study provided the Central Platte Natural Resources District with an advanced numerical modeling tool to assist with the update of their Groundwater Management Plan. An integrated hydrologic model, called the Central Platte Integrated Hydrologic Model, was constructed using the MODFLOW-One-Water Hydrologic Model code with the Newton solver. This code integrates climate, landscape, surface water, and groundwater-flow processes in a fully coupled approach. Model framework included 163 rows; 327 columns; 2,640 feet cell sides; and 3 vertical layers. A predevelopment model simulated steady-state hydrologic conditions prior to April 30, 1895, and a development period model discretized into 610 stress periods simulated transient hydrologic conditions from May 1, 1895, to December 31, 2016, using 170 biannual stress periods from 1895 to 1980, and monthly stress periods from May 1, 1980, to December 31, 2016. Calibration of the Central Platte Integrated Hydrologic Model involved two phases: a manual adjustment of parameters, followed by the automated calibration completed using BeoPEST that was facilitated by the employment of the singular value decomposition-assist features of PEST that specified 50 super parameters assembled from the 435 adjustable parameters and Tikhonov regularization. The average absolute groundwater-level residuals for model layers one, two, and three were 6.1, 12.4, and 7.4 feet, respectively. Calibrated horizontal hydraulic conductivity was about 70, 32, and 35 feet per day for layers 1, 2, and 3, respectively. The largest development period inflow to groundwater was recharge from deep percolation past the root zone, averaging 1,122,257 acre-feet per year (2.7 inches per year), and the largest outflow was to irrigation wells, averaging 693,171 acre-feet per year (10.2 inches per year for the Central Platte Natural Resources District). Other substantial groundwater outflows included evapotranspiration and base flow. For the total development period, there was a net change in storage of −122,393 acre-feet per year (−0.3 inch per year). The calibrated Central Platte Integrated Hydrologic Model was used to simulate eight different potential future climate and irrigation pumping conditions from January 1, 2017, to December 31, 2049. Simulated future groundwater levels within the Central Platte Natural Resources District varied significantly between scenarios and locally, from 13.8 feet below to 7.6 feet above baseline 1982 groundwater levels. Most areas exhibited groundwater-level declines for the drought scenarios and rises for the alternate irrigation scenarios. Changes in scenario groundwater levels correlated with the relations between farm net recharge and irrigation pumping. Linear “first order second moment” techniques indicated that the uncertainty in projected groundwater altitudes was reduced by 15.33 feet through model calibration.

Nebraska↗

Groundwater quality and geochemistry of West Virginia’s southern coal fields

Coal mining has been the dominant industry and land use in West Virginia’s southern coal fields since the mid-1800s. Mortality rates for a variety of serious chronic conditions, such as diabetes, heart disease, and some forms of cancer in Appalachian coal mining regions, are higher than in areas lacking substantial coal mining activity within the Appalachian Region or elsewhere in the United States. Causes of the increased mortality and morbidity are not clear, but poor diet, high rates of smoking, socioeconomic factors, and the quality of groundwater used by area residents are all possible contributing factors. This study was conducted by the U.S. Geological Survey in cooperation with the West Virginia Department of Health and Human Resources and the West Virginia Department of Environmental Protection, with grant support from the Centers for Disease Control and Prevention (CDC) to assess the quality of groundwater in southern West Virginia. The data from this assessment of groundwater quality may be used by the CDC and other agencies to potentially investigate the role or lack thereof of groundwater quality with respect to mortality and morbidity rates in the region. The study was conducted in a region where a high density of current or past coal mining combined with a lack of advanced sewage treatment could affect concentrations of commonly occurring constituents plus contaminants, including nitrate, trace metals, major ions, indicator bacteria, radon, hydrogen sulfide, and dissolved hydrocarbons. Because rural residential wells and mine outfalls are considered private sources of water in the region, and are therefore unregulated and unmonitored, water-quality data are sparse. To fill the data gap and assess the groundwater quality in the region, water-quality samples were collected from 60 sites in a 10-county area. The 60 sites sampled included 46 rural residential homeowner wells and 14 mine outfall discharges used for residential supply. For this study, all samples were collected prior to any filtration or other treatments, typically at the pressure tank, and are indicative of total and dissolved constituents in the untreated water. Generally, data for the 60 sites indicate that most waters sampled do not exceed thresholds for most U.S. Environmental Protection Agency (EPA) drinking-water standards and U.S. Geological Survey (USGS) drinking-water screening criteria. However, there were several notable exceptions. Turbidity exceeded the 5-Nephelometric Turbidity Unit (NTU) EPA treatment technique (TT) drinking-water standard in 14 of 60 (23 percent) sites sampled and exceeded the 1-NTU TT standard in 51 of 60 (85 percent) sites sampled. Turbidity is common in many wells in southern West Virginia and may be attributed to iron oxyhydroxide precipitates, sediment carried into the aquifers from the shallow soil zone due to improperly constructed or cased wells or transported to the aquifer in shallow stress-relief fracture zones or through permeable bedding-plane partings. For the sites sampled, 31 of 60 (52 percent) had pH values at, above, or below the upper and lower range of the EPA Secondary Maximum Contaminant Level (SMCL, 6.5–8.5 standard units). Of those 31 sites, 28 (90 percent) were indicative of acidic corrosive water and 3 (10 percent) were indicative of alkaline water. The Langelier Saturation Index (LSI), which is a measure of the corrosivity of the water, was computed for all sites sampled for the study. Eighty-two percent of the sites sampled had waters that were classified as corrosive, based on a LSI less than −0.5. Corrosive water has the potential to leach lead, copper, and other metals from lead, copper, galvanized, or lead-tin soldered connections in water lines. The chloride to sulfate mass ratio also was assessed with the alkalinity to indicate the potential to promote galvanic corrosion (PPGC) of water lines and plumbing fixtures. Only one of the sites (1.7 percent) classified as a corrosive water site, had a PPGC considered high; the remaining sites were classified as having either a moderate (53.3 percent) or low (45 percent) PPGC. Therefore, the type of plumbing systems sampled for this study may be affected by corrosive water, but the potential for leaching trace metals and other constituents from residential plumbing systems containing older galvanized pipes or lead-tin soldered copper pipes is moderate to low. The indicator bacteria total coliform and Escherichia coli ( E. coli ) also were detected in groundwater samples to varying degrees. Total coliforms, which are a broad class of indicator bacteria, are common in groundwater in southern West Virginia and were detected in 39 of the 60 sites (65 percent) sampled. The presence of total coliform bacteria is a potential indicator of surface contamination, due to improperly constructed or cased wells, or infiltration of soil or other surface contaminants into the aquifer or well bore. E. coli bacteria, however, are much more indicative of fecal contamination of groundwater from either human or animal sources, and 14 of the 60 (23 percent) sites sampled had detections of E. coli . Although only a few strains of E. coli are known pathogens, their presence in groundwater may be an indicator of other related pathogens such as viruses and should be regarded as a serious potential issue. Water treatment such as chlorination, ozonation, or ultraviolet light may be appropriate to kill potential pathogenic bacteria or viruses in the source water. Manganese and iron were prevalent contaminants in the groundwater samples collected for this study, with 30 of 60 (50 percent) sites analyzed for manganese and 25 of 60 (42 percent) sites analyzed for iron exceeding the proposed 50- and 300-micrograms per liter (µg/L) SMCL drinking-water standards, respectively, for aesthetic criteria such as taste, odor, or staining of plumbing fixtures. Fourteen of the 60 sites sampled (23 percent) had concentrations of manganese that exceeded the 300-µg/L USGS health-based screening level, and 1 site exceeded the 1,600-µg/L EPA drinking-water equivalent level, which is based on a lifetime exposure level. Sodium is another common constituent in groundwater within the study area. Sodium has an EPA health-based value (HBV) of 20 milligrams per liter (mg/L) for individuals who are on a sodium-restricted diet for blood pressure or other health reasons. Sodium concentrations exceeded the 20-mg/L EPA HBV in 27 of 60 (45 percent) samples. Radon, a naturally occurring carcinogenic radioactive gas known to cause lung cancer, was detected at concentrations at or exceeding the proposed 300-picocuries per liter (pCi/L) EPA Maximum Contaminant Level (MCL) in 12 of the 60 (20 percent) sites sampled. Sites with radon gas concentrations exceeding the 300-pCi/L proposed MCL have the potential for airborne concentrations of radon to exceed the 4-pCi/L indoor air standard. Inhalation of radon can cause lung cancer, and the 4-pCi/L indoor air standard is based on an inhalation standard. Therefore, homeowners whose wells have radon gas concentrations exceeding 300 pCi/L may be advised to have their indoor air tested to determine if indoor air concentrations exceed the 4-pCi/L indoor air standard established by the EPA. Various factors were analyzed statistically and graphically to determine whether they have an influence on groundwater quality within the study area, including topographic setting, well depth, type of mining (surface or underground), type of site (well or mine outfall), and geologic formation. Only geologic formation and the type of site sampled had strong statistical correlations with one or more of the constituents of concern for this study. The overall chemistry of outfalls (mine outfalls) and wells was significantly different, with a much higher dissolved oxygen content in outfalls than in wells. The dissolved oxygen content is the primary component driving the oxidation and reduction of minerals, and the precipitation of minerals that are saturated or super saturated with respect to various cations and anions. Median dissolved oxygen concentrations for the outfalls sampled was 8.75 mg/L, and only 0.4 mg/L for the wells sampled. Median concentrations of sulfate and selenium were much higher in waters from the outfalls sampled, with median concentrations of 73.75 mg/L and 2.35 µg/L, respectively, compared to the wells sampled, which had median concentrations of 18.3 mg/L and less than (<) the 0.05-µg/L method detection limit, respectively. The maximum selenium concentration was for a well, with a concentration of 16.6 µg/L. The geochemical processes that control sulfate and selenium concentrations in groundwater are similar and are the result of the oxidation of sulfide minerals such as pyrite and ferroselite. Iron and manganese concentrations were elevated in most of the wells sampled, with median concentrations of 269.5 and 124.5 µg/L, respectively, but were rarely detected in the outfalls sampled, with median concentrations of < 4.0 and < 0.4 µg/L, respectively. The difference in iron and manganese between wells and outfalls is indicative of the role of dissolved oxygen on processes controlling groundwater chemistry in the region. Three principal geologic formations were assessed for the study, and the overall chemistry for the Pocahontas, New River, and Kanawha Formations varied substantially with respect to several constituents. Concentrations of calcium, magnesium, and total dissolved solids were highest for sites sampled in the Pocahontas Formation, with median concentrations of 41.9, 18.6, and 312 mg/L, respectively. For constituents that are commonly associated with mining activity, the highest concentrations were for sites sampled in the New River Formation, with median concentrations of iron and manganese of 2,450 µg/L and 482 µg/L, respectively, and a median pH of 6.35 standard units. Concentrations of barium also were elevated in samples collected from sites in the New River Formation, with a median barium concentration of 184 µg/L. The source of the barium is not fully known but may be associated with commingling of shallow groundwater with deeper brines or dissolution of the mineral barite. The highest median sulfate concentrations were from sites sampled in the Pocahontas Formation, with a median concentration of 64.0 mg/L. Of the 12 sites at or exceeding the 300-pCi/L proposed drinking-water standard for radon, 8 (67 percent of MCL exceedances) were for sites deriving water from the Kanawha Formation, 3 (25 percent of MCL exceedances) were for sites deriving water from the New River Formation, and only 1 site was for water from the Pocahontas Formation (8 percent of proposed MCL exceedances). Dissolved hydrocarbons, including methane, ethane, propane, propene, n - and i -butane, 1-butene, n - and i -pentane, pentane, 2- and 3-ethyl pentane, hexane, and benzene were analyzed in samples collected from 59 of the 60 sites to assess the potential occurrence and sources of these trace gases in groundwater within the study area. Results of the analysis indicate that most of the gas is of shallow biogenic origin, possibly associated with coal-bed methane, but a subset of samples has a gas signature and a chloride to bromide ratio indicative of potential mixing with deeper thermogenic gases. Only 2 of the 59 (3.3 percent) sites sampled had concentrations of methane gas, which is a highly combustible and explosive gas, exceeding the 10 milligrams per kilogram level of concern established by the U.S. Office of Surface Mining Reclamation and Enforcement. Principal components analysis was used to assess the primary geochemical processes occurring in the aquifers sampled. The first principal component had significant positive loadings for bromide, chloride, silica, ammonia, barium, iron, manganese, and arsenic, and significant negative loadings for dissolved oxygen, potassium, nitrate, and uranium, and reflects reduction and oxidation (redox) processes occurring in deeper anoxic groundwater or shallow oxic groundwater. The strong positive loadings for iron, manganese, barium, and arsenic are correlated with reducing conditions often found deeper in the aquifer. More oxic water is correlated with oxidation of nitrogen species to nitrate and environmental mobilization of uranium and sulfate in shallow wells and mine outfalls.

West Virginia↗

Estimation of recharge rates to the sand and gravel aquifer using environmental tritium, Nantucket Island, Massachusetts

Estimation of the average annual rate of ground-water recharge to sand and gravel aquifers using elevated tritium concentrations in ground water is an alternative to traditional steady-state and water-balance recharge-rate methods. The concept of the tritium tracer method is that the average annual rate of ground-water recharge over a period of time can be calculated from the depth of the peak tritium concentration in the aquifer. Assuming that ground-water flow is vertically downward and that aquifer properties are reasonably homogeneous, and knowing the date of maximum tritium concentration in precipitation and the current depth to the tritium peak from the water table, the average recharge rate can be calculated. The method, which is a direct-measurement technique, was applied at two sites on Nantucket Island, Massachusetts. At site 1, the average annual recharge rate between 1964 and 1983 was 26.1 inches per year, or 68 percent of the average annual precipitation, and the estimated uncertainty is ?15 percent. At site 2, the multilevel water samplers were not constructed deep enough to determine the peak concentration of tritium in ground water. The tritium profile at site 2 resembles the upper part of the tritium profile at site 1 and indicates that the average recharge rate was at least 16 .7 inches per year, or at least 44 percent of the average annual precipitation. The Nantucket tritium recharge rates clearly are higher than rates determined elsewhere in southeastern Massachusetts using the tritium, water-table-fluctuation, and water-balance (Thornthwaite) methods, regardless of the method or the area. Because the recharge potential on Nantucket is so high (runoff is only 2 percent of the total water balance), the tritium recharge rates probably represent the effective upper limit for ground-water recharge in this region. The recharge-rate values used by Guswa and LeBlanc (1985) and LeBlanc (1984) in their ground-water-flow computer models of Cape Cod are 20 to 30 percent lower than this upper limit. The accuracy of the tritium method is dependent on two key factors: the accuracy of the effective-porosity data, and the sampling interval used at the site. For some sites, the need for recharge-rate data may require a determination as statistically accurate as that which can be provided by the tritium method. However, the tritium method is more costly and more time consuming than the other methods because numerous wells must be drilled and installed and because many water samples must be analyzed for tritium, to a very small level of analytical detection. For many sites, a less accurate, less expensive, and faster method of recharge-rate determination might be more satisfactory. The factor that most seriously limits the usefulness of the tritium tracer method is the current depth of the tritium peak. Water with peak concentrations of tritium entered the ground more than 20 years ago, and, according to the Nantucket data, that water now is more than 100 feet below the land surface. This suggests that the tracer method will work only in sand and gravel aquifers that are exceedingly thick by New England standards. Conversely, the results suggest that the method may work in areas where saturated thicknesses are less than 100 feet and the rate of vertical ground-water movement is relatively slow, such as in till and in silt- and clay-rich sand and gravel deposits.

Massachusetts↗

Plankton communities and summertime declines in algal abundance associated with low dissolved oxygen in the Tualatin River, Oregon

Phytoplankton populations in the Tualatin River in northwestern Oregon are an important component of the dissolved oxygen (DO) budget of the river and are critical for maintaining DO levels in summer. During the low-flow summer period, sufficient nutrients and a long residence time typically combine with ample sunshine and warm water to fuel blooms of cryptophyte algae, diatoms, green and blue-green algae in the low-gradient, slow-moving reservoir reach of the lower river. Algae in the Tualatin River generally drift with the water rather than attach to the river bottom as a result of moderate water depths, slightly elevated turbidity caused by suspended colloidal material, and dominance of silty substrates. Growth of algae occurs as if on a “conveyor belt” of streamflow, a dynamic system that is continually refreshed with inflowing water. Transit through the system can take as long as 2 weeks during the summer low-flow period. Photosynthetic production of DO during algal blooms is important in offsetting oxygen consumption at the sediment-water interface caused by the decomposition of organic matter from primarily terrestrial sources, and the absence of photosynthesis can lead to low DO concentrations that can harm aquatic life. The periods with the lowest DO concentrations in recent years (since 2003) typically occur in August following a decline in algal abundance and activity, when DO concentrations often decrease to less than State standards for extended periods (nearly 80 days). Since 2003, algal populations have tended to be smaller and algal blooms have terminated earlier compared to conditions in the 1990s, leading to more frequent declines in DO to levels that do not meet State standards. This study was developed to document the current abundance and species composition of phytoplankton in the Tualatin River, identify the possible causes of the general decline in algae, and evaluate hypotheses to explain why algal blooms diminish in midsummer. Plankton and water-quality sample data from 2006 to 2008 were combined with parts of a larger discrete-sample and continuous water-quality monitoring dataset and examined to identify patterns in water-quality and algal conditions since 1991, with a particular emphasis on 2003–08. Longitudinal plankton surveys were conducted in 2006–08 at six sites between river miles (RM) 24.5 and 3.4 at 2- to 3-week intervals, or 5–6 per season, and in-situ bioassay experiments were conducted in 2008 to examine the potential effects of wastewater treatment facility (WWTF) effluent and phosphorus additions on phytoplankton biomass and algal photosynthesis. Phytoplankton and zooplankton community composition, streamflow, and water-quality data were analyzed using multivariate statistical techniques to gain insights into plankton dynamics to determine what factors might be most tied to the abundance and characteristics of the phytoplankton assemblages, and identify possible causes of their declines. The connection between low-DO events and algal declines was clearly evident, as bloom crashes were nearly always followed by periods of low DO. Algal blooms occurred each year during 2006–08, producing maximum chlorophyll-a (Chl-a) values in June or July generally in the range of 50–80 micrograms per liter (µg/L). Bloom crashes and absence of sufficient algal photosynthesis in mid- to late-summer contributed to minimum DO concentrations that were less than the State standard of 6.5 milligrams per liter (mg/L) based on the 30-day mean daily concentration, for 62–74 days each year. At times, the absolute minimum State standard (4 mg/L DO) also was not met. To learn more about why low-DO events occurred, specific algal declines during 2003–08 were scrutinized to determine their likely causal factors. From this information, a series of hypotheses were formulated and evaluated in terms of their ability to explain recent declines in algal populations in the river in late summer. Meteorological, streamflow, turbidity, water temperature, and conductance conditions in the Tualatin River during the 2006–08 summer seasons were not atypical. Natural flow comprised the majority (70–80 percent) of flow each year during spring, but then reduced to 38–40 percent during midsummer when WWTF effluent—which contributed as much as 36 percent—and flow augmentation releases comprised a greater fraction of the flow. Summer 2008 was unusual, however, in the prolonged influence from the Wapato Lake agricultural area near Gaston in the upper part of the basin. The previous winter flooding and levee breach at Wapato Lake caused a much greater area of inundation. As a result, drainage from this area continued into July, much later than normal. A subsequent algal bloom in Wapato Lake then seeded the upper Tualatin River, and this drainage had a profound effect on the downstream plankton community. A large blue-green algae bloom developed—the largest in recent memory—prompting a public health advisory for recreational contact for about two weeks. Algal growths and surface blooms are a common feature of the Tualatin River. Most of the dominant algae have growth forms and morphologies that are well suited for planktonic life, employing spines and gas vacuoles to resist settling, forming colonies, and producing mucilage (or toxins) to resist zooplankton grazing. In 2006–08, 143 algal taxa were identified in 117 main-stem samples; diatoms and green algae were more diverse than blue-green, golden, and cryptophyte algae, although these later groups sometimes dominated the overall volumetric abundance (biovolume). The most frequently occurring taxa, occurring in 97–99 percent of samples, were flagellated cryptophytes Cryptomonas erosa and Rhodomonas minuta. Other important algal taxa included centric diatoms Stephanodiscus, Cyclotella, and Melosira species and colonial green algae Scenedesmus and Actinastrum. These taxa comprised the majority of the algal biovolume during much of the growing season. A general seasonal trend in the phytoplankton assemblages was observed, with dominance by filamentous centric diatoms Stephanodiscus and Melosira in spring and early summer, and flagellated cryptophytes and green algae, particularly Chlamydomonas sp., in late-summer; or, in 2008, dominance by blue-green algae Anabaena flos-aquae and Aphanizomenon flos-aquae during the Wapato Lake bloom event. There were 99 zooplankton taxa identified from the Tualatin River in 2006–08, composed primarily of cladocerans, copepods, and rotifers. A seasonal increase in zooplankton abundance was observed in early summer just as or shortly after the phytoplankton population began to increase, with populations growing to 15,000−120,000 organisms per cubic meter in the lower river. Zooplankton abundance showed a predictable and distinct longitudinal downstream increase, particularly downstream of Highway 99W (RM 11.6). Although grazing rates were not measured, the data suggest that, at times, zooplankton grazing may affect algal abundance and species composition in the Tualatin River, with diatoms becoming relatively less abundant and flagellated cryptophytes and green algae relatively more abundant during periods when zooplankton densities were highest. Multivariate statistical analyses identified soluble reactive phosphorus (SRP), natural flow, flow augmentation, and WWTF effluent as important factors influencing Tualatin River phytoplankton populations, with zooplankton density (particularly rotifers and copepods), specific conductance, chloride, and water temperature also having an important influence. Although SRP was highly correlated with the plankton communities, that correlation was likely the result of high or low algal activity (uptake) as SRP concentrations were often reduced to low levels during blooms. While previous studies have already established that phosphorus, among other factors such as flow, places a theoretical cap on the size of the phytoplankton population in the river, sometimes algal declines occur when SRP concentrations are apparently sufficient. To identify alternative causal factors, additional analyses were performed without SRP to focus on other water-quality parameters, zooplankton density, and flow factors. Considering data for all 3 years and including just those samples from the lower Tualatin River not affected by the 2008 Wapato Lake drainage event, three factors (percentage of reservoir flow augmentation, total natural flow, and rotifer density) best explained variations in the phytoplankton assemblages. Analyses focusing on the possible causes of algal declines included the above multivariate analyses, scrutiny of 10 specific instances of declines in algal populations during 2003–08 including several bloom–crash sequences, and analyses of historic routine watershed monitoring data from Clean Water Services. Six factors were hypothesized to be important in causing bloom crashes or impeding blooms from rebounding in August: (1) light limitation from cloudy weather, (2) a reduction in the plankton inocula or “seed” entering the lower river from upstream sources, (3) increased summer streamflows, (4) changes in the dominant sources of flow as the percentage of flow augmentation and WWTF discharges have increased, (5) zooplankton grazing, and (6) low concentrations of bioavailable phosphorus (<0.015 milligram per liter). All of these hypotheses are supported in some fashion by the available data and statistical analyses. Zooplankton grazing, short-term declines in photosynthesis from cloudy weather, total flow as it affects residence time, and the dominant source of flow are primary factors responsible for the low-DO events caused by declines in algae in the lower Tualatin River during late summer. Cloudy weather and increased turbidity are known to inhibit algal growth in the Tualatin River, and slight increases in turbidity in recent years may be a problem. Upstream sources of algae are critical in determining the characteristics and size of downstream populations, as illustrated by the Wapato Lake bloom in 2008, but more data are needed from upstream to fully define the importance of this connection. The sources of flow, through their differential contribution of plankton inocula (quality and amount), were, at times, important factors affecting phytoplankton populations. While SRP concentrations were often most highly correlated with phytoplankton species community, the bioavailability of phosphorus is still somewhat unknown and there are several sources to consider. Preliminary bioassay tests suggested that while treated wastewater effluent may stimulate algae at 30 percent concentrations, negative effects (or decreased stimulation) on Chl-a and DO production may occur at concentrations of 50 percent. Targeted data collection and future research will be needed to further understand the importance of these factors on Tualatin River phytoplankton. While the data and analysis completed for this report provide insights into future research and monitoring that would be useful to continue, additional monitoring of turbidity, Chl-a, and plankton abundance and species composition in the upper part of the basin would enhance our understanding of plankton dynamics and factors affecting phytoplankton abundance in the lower river. Assessment of the key upstream sources of algal inocula via surveys of the major flow sources as well as tributaries and wetlands would provide useful information for the management of river water quality. Other studies that could prove useful for developing management strategies include targeted experiments to evaluate the bioavailability of phosphorus from a variety of sources. New research on phytoplankton–zooplankton interactions, and studies of planktivorous fish, might also provide insight about food web dynamics and potential “top-down” effects of fish predation on the plankton communities. In addition, further development of neural-network or other water-quality models would help to evaluate management strategies and provide forecasts of water-quality conditions. Finally, periodic future reassessments of the available data with the multivariate statistical tools used in this study would be helpful to assess whether and how plankton communities are changing, and to continue to shed light on the importance of factors shaping the plankton. Although certain types and sizes of algal blooms are undesirable, minimum phytoplankton populations are an important part of aquatic food webs and are needed to maintain healthy levels of DO in the river. By understanding the sources, characteristics, causal factors, and responses of the plankton communities, management strategies can be developed to improve DO conditions in the lower Tualatin River during the important summer low-flow period.

Oregon↗

Radiometric calibration of a non-imaging airborne spectrometer to measure the Greenland ice sheet surface

Methods to radiometrically calibrate a non-imaging airborne visible-to-shortwave infrared (VSWIR) spectrometer to measure the Greenland ice sheet surface are presented. Airborne VSWIR measurement performance for bright Greenland ice and dark bare rock/soil targets is compared against the MODerate resolution atmospheric TRANsmission (MODTRAN ® ) radiative transfer code (version 6.0), and a coincident Landsat 8 Operational Land Imager (OLI) acquisition on 29 July 2015 during an in-flight radiometric calibration experiment. Airborne remote sensing flights were carried out in northwestern Greenland in preparation for the Ice, Cloud, and land Elevation Satellite 2 (ICESat-2) laser altimeter mission. A total of nine science flights were conducted over the Greenland ice sheet, sea ice, and open-ocean water. The campaign's primary purpose was to correlate green laser pulse penetration into snow and ice with spectroscopic-derived surface properties. An experimental airborne instrument configuration that included a nadir-viewing (looking downward at the surface) non-imaging Analytical Spectral Devices (ASD) Inc. spectrometer that measured upwelling VSWIR (0.35 to 2.5 µ m) spectral radiance ( W m - 2 sr - 1 &#xB5; m - 1 "> W m − 2 sr − 1 µ m − 1 ) in the two-color Slope Imaging Multi-polarization Photon-Counting Lidar's (SIMPL) ground instantaneous field of view, and a zenith-viewing (looking upward at the sky) ASD spectrometer that measured VSWIR spectral irradiance (W m −2 nm −1 ) was flown. National Institute of Standards and Technology (NIST) traceable radiometric calibration procedures for laboratory, in-flight, and field environments are described in detail to achieve a targeted VSWIR measurement requirement of within 5 % to support calibration/validation efforts and remote sensing algorithm development. Our MODTRAN predictions for the 29 July flight line over dark and bright targets indicate that the airborne nadir-viewing spectrometer spectral radiance measurement uncertainty was between 0.6 % and 4.7 % for VSWIR wavelengths (0.4 to 2.0 µ m) with atmospheric transmittance greater than 80 %. MODTRAN predictions for Landsat 8 OLI relative spectral response functions suggest that OLI is measuring 6 % to 16 % more top-of-atmosphere (TOA) spectral radiance from the Greenland ice sheet surface than was predicted using apparent reflectance spectra from the nadir-viewing spectrometer. While more investigation is required to convert airborne VSWIR spectral radiance into atmospherically corrected airborne surface reflectance, it is expected that airborne science flight data products will contribute to spectroscopic determination of Greenland ice sheet surface optical properties to improve understanding of their potential influence on ICESat-2 measurements.

Atmospheric Measurement Techniques↗

The Connecticut Streamflow and Sustainable Water Use Estimator—A decision-support tool to estimate water availability at ungaged stream locations in Connecticut

Freshwater streams in Connecticut are subject to many competing demands, including public water supply; agricultural, commercial, and industrial water use; and ecosystem and habitat needs. In recent years, drought has further stressed Connecticut’s water resources. To sustainably allocate and manage water resources among these competing uses, Federal, State, and local water-resource managers require data and modeling tools to estimate the water availability at a variety of temporal and spatial scales for planning purposes. The Connecticut Streamflow and Sustainable Water Use Estimator (CT SSWUE), developed by the U.S. Geological Survey in cooperation with the Connecticut Department of Energy and Environmental Protection, is a decision-support tool for estimating daily unaltered streamflow and sustainable water use at ungaged sites in Connecticut. The CT SSWUE estimates unaltered daily mean streamflow and water-use-adjusted streamflow for the period from October 1, 1960, to September 30, 2015, and the monthly sustainable net withdrawal at ungaged sites in Connecticut. Unaltered streamflow is the estimated daily mean streamflow in a drainage basin in the absence of any water withdrawals or wastewater discharges and with minimal human development. Sustainable net withdrawal is the maximum net withdrawal (withdrawal minus wastewater discharges) that can be drawn from a basin without critically depleting the water available through natural streamflow patterns. Sustainable net withdrawal is defined for this study as the difference between the unaltered daily mean streamflow and a user-defined target minimum streamflow. Weighted least squares and Tobit regression techniques were used to develop equations for estimating streamflow at ungaged sites at 19 streamflow quantiles with exceedance probabilities ranging from 0.005 to 99.995 percent. Regressions were based on streamflow quantiles and basin characteristics from 36 reference streamgages in and around Connecticut. Four basin characteristics—drainage area, mean of the soil permeability, mean of the average annual precipitation, and ratio of the length of streams that overlay sand and gravel deposits to the total length of streams in the basin—are used as explanatory variables in the equations. At an ungaged site, interpolation between the streamflow quantiles estimated from the regression equations produces a continuous flow-duration curve. A time series of daily mean streamflow at an ungaged site is then estimated by assuming that for each day, the streamflow quantile occurs on the same date at both a reference streamgage and the ungaged site. In a remove-one cross validation, estimated unaltered daily mean streamflow agreed well with observed values at reference streamgages, with a few exceptions. Nash Sutcliffe efficiency ranged from −0.43 to 0.97 with a median value of 0.88. The normalized root-mean-square error ranged from 16.6 to 120.4 percent with a median value of 34.5 percent. An empirical method for estimating 95-percent prediction intervals for unaltered daily and monthly mean streamflow was developed and tested by using the cross-validation data. Prediction intervals for unaltered daily mean streamflow at the cross-validation reference streamgages performed well in most cases. Gaged streamflow values from the cross-validation data fell within the prediction intervals a median 96.6 percent of the time for daily mean time series and 93.9 percent of the time for monthly mean time series. The CT SSWUE computes water-use-adjusted streamflow using spatially referenced water-use information provided by the Connecticut Department of Energy and Environmental Protection. Available water-use information included permitted and registered water withdrawals and permitted wastewater discharges during 1998 to 2015 for the Thames River Basin and central coastal drainage basins. Water-use information was incorporated into the U.S. Geological Survey StreamStats web application for Connecticut and can be used for computing water-use-adjusted streamflow and sustainable net withdrawal at selected points of interest. Altered daily streamflow is computed by applying average daily withdrawals and wastewater discharges to the water balance equation. Average daily surface water withdrawals and wastewater discharges are applied directly to the daily water balance equation. Time-lagged alterations on streamflow from groundwater withdrawals or wastewater discharges are estimated by using a response-coefficient method developed from results of previously published, calibrated groundwater models.

Connecticut↗

Ground-water hydrology and simulation of ground-water flow at Operable Unit 3 and surrounding region, U.S. Naval Air Station, Jacksonville, Florida

The Naval Air Station, Jacksonville (herein referred to as the Station), occupies 3,800 acres adjacent to the St. Johns River in Duval County, Florida. Operable Unit 3 (OU3) occupies 134 acres on the eastern side of the Station and has been used for industrial and commercial purposes since World War II. Ground water contaminated by chlorinated organic compounds has been detected in the surficial aquifer at OU3. The U.S. Navy and U.S. Geological Survey (USGS) conducted a cooperative hydrologic study to evaluate the potential for ground water discharge to the neighboring St. Johns River. A ground-water flow model, previously developed for the area, was recalibrated for use in this study. At the Station, the surficial aquifer is exposed at land surface and forms the uppermost permeable unit. The aquifer ranges in thickness from 30 to 100 feet and consists of unconsolidated silty sands interbedded with local beds of clay. The low-permeability clays of the Hawthorn Group form the base of the aquifer. The USGS previously conducted a ground-water investigation at the Station that included the development and calibration of a 1-layer regional ground-water flow model. For this investigation, the regional model was recalibrated using additional data collected after the original calibration. The recalibrated model was then used to establish the boundaries for a smaller subregional model roughly centered on OU3. Within the subregional model, the surficial aquifer is composed of distinct upper and intermediate layers. The upper layer extends from land surface to a depth of approximately 15 feet below sea level; the intermediate layer extends from the upper layer down to the top of the Hawthorn Group. In the northern and central parts of OU3, the upper and intermediate layers are separated by a low-permeability clay layer. Horizontal hydraulic conductivities in the upper layer, determined from aquifer tests, range from 0.19 to 3.8 feet per day. The horizontal hydraulic conductivity in the intermediate layer, determined from one aquifer test, is 20 feet per day. An extensive stormwater drainage system is present at OU3 and the surrounding area. Some of the stormwater drains have been documented to be draining ground water from the upper layer of the surficial aquifer, whereas other drains are only suspected to be draining ground water. The subregional model contained 78 rows and 148 columns of square model cells that were 100 feet on each side. Vertically, the surficial aquifer was divided into two layers; layer 1 represented the upper layer and layer 2 represented the intermediate layer. Steady-state ground-water flow conditions were assumed. The model was calibrated to head data collected on October 29 and 30, 1996. After calibration, the model matched all 67 measured heads to within the calibration criterion of 1 foot; and 48 of 67 simulated heads (72 percent) were within 0.5 foot. Model simulated recharge rates ranged from 0.4 inch per year in areas that were largely paved to 13.0 inches per year in irrigated areas. Simulated hydraulic conductivities in the upper layer at OU3 ranged from 0.5 foot per day in the north to 1.0 foot per day in the south. Simulated vertical leakance between the upper and intermediate layers ranged from 1.0x10-6 per day in an area with low-permeability clays to 4.3x10-2 per day in an area that had been dredged. Simulated transmissivities in the intermediate layer ranged from 25 feet squared per day in an area of low-permeability channel-fill deposits to a high of 1,200 feet squared per day in areas covering most of OU3. Simulated riverbed conductances ranged from 4 to 60 feet squared per day and simulated bottom conductances of leaking stormwater drains ranged from 5 to 20 feet squared per day. The direction and velocity of ground-water flow was determined using particle-tracking techniques. Ground-water flow in the upper layer was generally eastward toward the St. Johns River. However, leaking stormwater drains locally modified the flow system to create small areas with flow that was diverted to the drains. The flow velocities in the upper layer at OU3 were slow, averaging about 2 feet per year. The slow velocities were primarily the result of the low horizontal hydraulic conductivity and, secondarily, the result of the low recharge rate. The simulated rate at which ground water leaked into the stormwater drains was low, averaging about 0.0011 cubic feet per second per 100 feet of stormwater drainage conduit. Ground-water flow in the intermediate layer moved eastward toward and discharged into the St. Johns River. Flow velocities were significantly higher in this layer than in the upper layer. The velocity was about 35 and 12 feet per year in the northern and southern parts of OU3, respectively.

Florida↗

Geographic relations of landslide distribution and assessment of landslide hazards in the Blanco, Cibuco, and Coamo basins, Puerto Rico

Landslide occurrence is common in mountainous areas of Puerto Rico where mean annual rainfall and the frequency of intense storms are high and hillslopes are steep. Each year, landslides cause extensive damage to property and occasionally result in loss of life. Landslide maps developed from 1:20,000 scale aerial photographs in combination with a computerized geographic information system were used to evaluate the landslide potential in the Blanco, Cibuco, and Coamo Basins of Puerto Rico. These basins, ranging in surface area from 276 to 350 square kilometers, are described in this report. The basins represent a broad range of the climatologic, geographic, and geologic conditions that occur in Puerto Rico. In addition, a variety of landslide types were documented. Rainfall-triggered debris flows, shallow soil slips, and slumps were most abundant. The most important temporal control on landslide occurrence in Puerto Rico is storm rainfall. Forty-one storms triggered widespread landsliding about 1 to 2 times per year during the last three decades. These storms were frequently of 1 to 2 days duration in which, on average, several hundred millimeters of rainfall triggered tens to hundreds of landslides in the central mountains. Most of these storms were tropical disturbances that occurred during the hurricane season of June through November. Land use and the topographic characteristics of hillslope angle, elevation, and aspect are the most important spatial controls governing landslide frequency. Hillslopes in the study area that have been anthropogenically modified, exceed 12 degrees in gradient and about 350 meters in elevation, and face the east-northeast are most prone to landsliding. Bedrock geology and soil order seem less important in the determination of landslide frequency, at least when considered at a generalized level. A rainfall accumulation-duration relation for the triggering of numerous landslides throughout the central mountains, and a set of simplified matrices representing geographic conditions in the three river basins were developed and are described in this report. These two elements provide a basis for the estimation of the temporal and spatial controls on landslide occurrence in Puerto Rico. Finally, this approach is an example of a relatively inexpensive technique for landslide hazard analysis that may be applicable to other settings.

Water-Resources Investigations Report↗

Methods for estimating flow-duration and annual mean-flow statistics for ungaged streams in Oklahoma

Flow statistics can be used to provide decision makers with surface-water information needed for activities such as water-supply permitting, flow regulation, and other water rights issues. Flow statistics could be needed at any location along a stream. Most often, streamflow statistics are needed at ungaged sites, where no flow data are available to compute the statistics. Methods are presented in this report for estimating flow-duration and annual mean-flow statistics for ungaged streams in Oklahoma. Flow statistics included the (1) annual (period of record), (2) seasonal (summer-autumn and winter-spring), and (3) 12 monthly duration statistics, including the 20th, 50th, 80th, 90th, and 95th percentile flow exceedances, and the annual mean-flow (mean of daily flows for the period of record). Flow statistics were calculated from daily streamflow information collected from 235 streamflow-gaging stations throughout Oklahoma and areas in adjacent states. A drainage-area ratio method is the preferred method for estimating flow statistics at an ungaged location that is on a stream near a gage. The method generally is reliable only if the drainage-area ratio of the two sites is between 0.5 and 1.5. Regression equations that relate flow statistics to drainage-basin characteristics were developed for the purpose of estimating selected flow-duration and annual mean-flow statistics for ungaged streams that are not near gaging stations on the same stream. Regression equations were developed from flow statistics and drainage-basin characteristics for 113 unregulated gaging stations. Separate regression equations were developed by using U.S. Geological Survey streamflow-gaging stations in regions with similar drainage-basin characteristics. These equations can increase the accuracy of regression equations used for estimating flow-duration and annual mean-flow statistics at ungaged stream locations in Oklahoma. Streamflow-gaging stations were grouped by selected drainage-basin characteristics by using a k-means cluster analysis. Three regions were identified for Oklahoma on the basis of the clustering of gaging stations and a manual delineation of distinguishable hydrologic and geologic boundaries: Region 1 (western Oklahoma excluding the Oklahoma and Texas Panhandles), Region 2 (north- and south-central Oklahoma), and Region 3 (eastern and central Oklahoma). A total of 228 regression equations (225 flow-duration regressions and three annual mean-flow regressions) were developed using ordinary least-squares and left-censored (Tobit) multiple-regression techniques. These equations can be used to estimate 75 flow-duration statistics and annual mean-flow for ungaged streams in the three regions. Drainage-basin characteristics that were statistically significant independent variables in the regression analyses were (1) contributing drainage area; (2) station elevation; (3) mean drainage-basin elevation; (4) channel slope; (5) percentage of forested canopy; (6) mean drainage-basin hillslope; (7) soil permeability; and (8) mean annual, seasonal, and monthly precipitation. The accuracy of flow-duration regression equations generally decreased from high-flow exceedance (low-exceedance probability) to low-flow exceedance (high-exceedance probability) . This decrease may have happened because a greater uncertainty exists for low-flow estimates and low-flow is largely affected by localized geology that was not quantified by the drainage-basin characteristics selected. The standard errors of estimate of regression equations for Region 1 (western Oklahoma) were substantially larger than those standard errors for other regions, especially for low-flow exceedances. These errors may be a result of greater variability in low flow because of increased irrigation activities in this region. Regression equations may not be reliable for sites where the drainage-basin characteristics are outside the range of values of independent vari

Scientific Investigations Report↗

Monitoring marine eruptions

Introduction Submarine volcanoes produce much of the same seismicity and eruptive activity as subaerial volcanoes and can pose hazards to society. Although they can be monitored with similar techniques and methods as described in other chapters of this volume, their submerged location brings unique challenges. This chapter addresses these challenges and provides recommendations for monitoring volcanoes fully or partly in marine environments to meet the capabilities described in other chapters of this volume. The United States and its territories host dozens of submarine volcanoes with most (around 60) in the Commonwealth of the Northern Mariana Islands. Approximately 20 of the Northern Mariana Islands submarine volcanoes are known to be hydrothermally active, and 10 have confirmed eruptions since the 1950s (for example, Baker and others, 2008; Tepp and others, 2019a). Nine of those volcanoes were considered by the National Volcanic Threat Assessment (Ewert and others, 2018) to have a combination of eruptive type and summit depth that poses a higher risk of hazardous eruptions, although only one was listed as a moderate (level 3) threat. Other notable submarine volcanoes of interest to the United States that have historically erupted are Axial Seamount off the Washington State coast, Kamaʻehuakanaloa in Hawaiʻi, and Vailuluʻu seamount in American Samoa. All of these, however, have a low risk of hazards because of their depth (greater than 600 meters below sea level) and eruptive type and so are not included in the National Volcanic Threat Assessment. In addition to submarine volcanoes, the submerged flanks of island volcanoes can also be a source of hazardous submarine eruptions—for example, the 1877 eruption of Mauna Loa, Hawai‘i, in Kealakekua Bay (Wanless and others, 2006). The most notable submarine eruption in recent times was the 2022 eruption of Hunga Tonga–Hunga Haʻapai in Tonga, which was one of the largest eruptions on Earth in the past 100 years. It created a massive volcanic plume, unprecedented shock waves, and far-reaching tsunami (Lynett and others, 2022). Other recent submarine eruptions in the Pacific Ocean Basin have produced subaerial plumes that reached aircraft heights (Carey and others, 2014) and large pumice rafts that can affect marine traffic and harbors (for example, Jutzeler and others, 2014; Kornei, 2019). These examples illustrate the potential hazards of major submarine eruptions. Yet, submarine volcanoes are largely unmonitored, and many eruptions occur that are unnoticed or only identified hours or days afterward. Within U.S. territory, submarine volcanoes in the Northern Mariana Islands have been known to produce eruptive activity that can affect society. Reports from fishermen and other marine vessels in the Northern Mariana Islands have noted underwater explosions, sea-surface discoloration, and bubbling water, all of which are known to be signs of submarine volcanic activity. South Sarigan seamount, located about 160 kilometers (km) north of Saipan, erupted in 2010 from greater than 150 meters below the sea surface, resulting in a gas and ash plume that reached more than 11.9 km into the atmosphere (for example, Searcy, 2013; Embley and others, 2014), high enough to affect international air traffic. Precursory and co-eruptive seismicity was detected on the regional Northern Mariana Islands seismic network (Searcy, 2013) and on global monitoring instruments (Green and others, 2013). Monitoring of submarine volcanoes is best accomplished with marine-based instrumentation, which is also useful for monitoring small island volcanoes that may not have the land area necessary for comprehensive subaerial monitoring. The primary marine-based instrumentation used for submarine volcanoes includes ocean-bottom pressure sensors to assess sea-floor deformation, ocean-bottom seismometers (OBSs) to detect seismicity, and both moored and ocean-bottom hydrophones to detect submarine explosions. Other sensors offer important monitoring data, such as turbidity, temperature, and chemistry of hydrothermal emissions. Marine-based instruments are typically deployed in campaign-style networks with no real-time telemetry owing to cost considerations and technical limitations. However, when necessary, marine instruments can be operated in real time using cables to transmit data to land-based facilities; other technologies for this purpose are in use or in development, such as acoustic transmission from the instrument to a moored buoy (Matsumoto and others, 2016) and a winch-based system with a satellite antenna that is part of the instrument mooring (Matsumoto and others, 2019). Emerging technologies for marine-based monitoring may be considered as part of a future monitoring plan. These technologies include ocean gliders and floats with on-board hydrophones that have been used to record earthquakes and submarine eruptions (for example, Matsumoto and others, 2013; Sukhovich and others, 2015) and fiber-optic cables that have been used as strainmeters to detect earthquakes (for example, Marra and others, 2018; Lindsey and others, 2019). Land-based instruments and satellites can also provide some capability for monitoring submarine volcanoes, but they provide more limited observations than marine-based instrumentation.

Scientific Investigations Report↗

Water-quality conditions during low flow in the lower Youghiogheny River basin, Pennsylvania, October 5-7, 1998

In October 1998, a chemical synoptic survey was conducted by the U.S. Geological Survey, in cooperation with the U.S. Department of Energy, National Energy Technology Laboratory, in the Lower Youghiogheny River Basin in Pennsylvania to give a snapshot of present (1998) water quality during low-flow conditions. Water samples from 38 sites—12 mainstem sites, 22 tributaries, and 4 mine discharges that discharge directly to the Youghiogheny River—were used to identify sources of contaminants from mining operations. Specific conductance, water temperature, pH, and dissolved oxygen were measured in the field at each site and concentrations of major ions and trace elements were measured in the laboratory.z Unaccounted for gains and losses in streamflow were measured during the study. Unaccounted for losses in streamflow might be attributed to water loss through streambed fractures. Extensive mine tunnels are present in the basin and loss of water to these tunnels seems likely. Unaccounted for gains in streamflow may be from unmeasured tributaries or surface seeps, but most of the gains are suspected to come from artesian flow through fractures in the streambed from underground mine pools. Influent flows of rust-colored water were noted in some river sections. The pH values for all the samples collected during this survey were above 5.8, and most (33 of 38 samples) were above 7.0. Samples from the four mine-discharge sites also had pH values between 6.3 and 6.7. The lowest pH (5.8) was in a tributary, Galley Run. All 38 sampling sites had net alkalinity. The alkalinity load in the Youghiogheny River increased between Connellsville and McKeesport from 35 to 79 tons per day. Above Smithton, the measured alkalinity load in the Lower Youghiogheny River agreed well with the estimated alkalinity load. Below Smithton, measured alkalinity loads in the Lower Youghiogheny River are greater than calculated loads, resulting in unaccounted for gains in alkalinity. These gains are believed to be from seeps in the streambed. Approximately one-third of the load of total alkalinity in the Youghiogheny River at McKeesport is attributed to Sewickley Creek, which contributes 14 tons per day. Sulfate concentrations in the Youghiogheny River steadily increase from 33 milligrams per liter at Connellsville to 77 milligrams per liter near McKeesport. The measured concentrations of sulfate exceeded Pennsylvania water-quality standards at four tributary sites (Galley Run, Hickman Run, Sewickley Creek, and Gillespie Run) and all four mine-discharge sites but not at any main-stem sites. A large increase in sulfate load between West Newton and Sutersville can be attributed almost entirely to the contribution from Sewickley Creek (49 tons per day). Approximately 25 percent of the load measured between Connellsville and McKeesport is unaccounted for. These gains are believed to be from seeps in the streambed from underground mine pools. Similar patterns also were observed for loads of sodium, calcium, and magnesium. Unmeasured inputs from mine drainage are believed to be the source of these loads. Elevated concentrations (above background levels) of chemicals associated with drainage from coal-mining operations were measured in samples from tributaries, especially from Galley Run, Gillespie Run, and Sewickley Creek, and from the mine-discharge sites. The synoptic survey conducted for this study was successful in identifying generalized reaches of the Youghiogheny River where unaccounted for loads of constituents associated with mining activities are entering the river. However, the survey was not able to pinpoint the location of these loads. Remote-sensing techniques, such as thermal infrared imaging by the National Energy Technology Laboratory, could be useful for determining the precise locations of these inputs.

Water-Resources Investigations Report↗