USGS Science⌕ Search

SEARCH · USGS Science

Results for “Scientific Data - Nature”

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 1,063 records · Page 59Linked to original sources

Testing models of Laramide orogenic initiation by investigation of Late Cretaceous magmatic-tectonic evolution of the central Mojave sector of the California arc

The Mojave Desert region is in a critical position for assessing models of Laramide orogenesis, which is hypothesized to have initiated as one or more seamounts subducted beneath the Cretaceous continental margin. Geochronological and geochemical characteristics of Late Cretaceous magmatic products provide the opportunity to test the validity of Laramide orogenic models. Laramide-aged plutons are exposed along a transect across the Cordilleran Mesozoic magmatic system from Joshua Tree National Park in the Eastern Transverse Ranges eastward into the central Mojave Desert. A transect at latitude ∼33.5°N to 34.5°N includes: (1) the large upper-crustal Late Cretaceous Cadiz Valley batholith, (2) a thick section of Proterozoic to Jurassic host rocks, (3) Late Cretaceous stock to pluton-sized bodies at mesozonal depths, and (4) a Jurassic to Late Cretaceous midcrustal sheeted complex emplaced at ∼20 km depth that transitions into a migmatite complex truncated along the San Andreas fault. This magmatic section is structurally correlative with the Big Bear Lake intrusive suite in the San Bernardino Mountains and similar sheeted rocks recovered in the Cajon Pass Deep Scientific Drillhole. Zircon U-Pb geochronology of 12 samples via secondary ionization mass spectrometry (SIMS) (six from the Cadiz Valley batholith and six from the Cajon Pass Deep Scientific Drillhole) indicates that all Cretaceous igneous units investigated were intruded between 83 and 74 Ma, and Cajon Pass samples include a Jurassic age component. A compilation of new and published SIMS geochronological data demonstrates that voluminous magmatism in the Eastern Transverse Ranges and central Mojave Desert was continuous throughout the period suggested for the intersection and flat-slab subduction of the Shatsky Rise conjugate deep into the interior of western North America. Whole-rock major-element, trace-element, and isotope geochemistry data from samples from a suite of 106 igneous rocks represent the breadth of Late Cretaceous units in the transect. Geochemistry indicates an origin in a subduction environment and intrusion into a crust thick enough to generate residual garnet. The lack of significant deflections of compositional characteristics and isotopic ratios in igneous products through space and time argues against a delamination event prior to 74 Ma. We argue that Late Cretaceous plutonism from the Eastern Transverse Ranges to the central Mojave Desert represents subduction zone arc magmatism that persisted until ca. 74 Ma. This interpretation is inconsistent with the proposed timing of the docking of the Shatsky Rise conjugate with the margin of western North America, particularly models in which the leading edge of the Shatsky Rise was beneath Wyoming at 74 Ma. Alternatively, the timing of cessation of plutonism precedes the timing of the passage of the Hess Rise conjugate beneath western North America at ca. 70–65 Ma. The presence, geochemical composition, and age of arc products in the Eastern Transverse Ranges and central Mojave Desert region must be accounted for in any tectonic model of the transition from Sevier to Laramide orogenesis.

California↗

Rapid Method for Escherichia coli in the Cuyahoga River

This study is a continuation of a previous U.S. Geological Survey (USGS) project in cooperation with the National Park Service at Cuyahoga Valley National Park in Brecksville, Ohio. A rapid (1-hour) method for detecting Escherichia coli (E. coli) in water was tested and compared to the standard (24-hour) method for determining E. coli concentrations. Environmental data were collected to determine turbidity, rainfall, and streamflow at the time of sampling. In the previous study (2004-5), data collected were used to develop predictive models to determine recreational water quality in the river at two sites within the park. Data collected during this continued study (2006) were used to test these models. At Jaite, a centrally located site within the park, the model correctly predicted exceedances or nonexceedances of the Ohio Environmental Protection Agency maximum for recreational water quality in 80 percent of samples. At Old Portage, a site near the upstream boundary of the park, the model correctly predicted recreational water quality in 58 percent of samples. All of the data collected in 2004-6 will be used to develop more accurate models for use in future studies. Analysis and discussion of model results are scheduled to be included in an upcoming USGS Scientific Investigations Report.

Open-File Report↗

Conifer health classification for Colorado, 2008

Colorado has undergone substantial changes in forests due to urbanization, wildfires, insect-caused tree mortality, and other human and environmental factors. The U.S. Geological Survey Rocky Mountain Geographic Science Center evaluated and developed a methodology for applying remotely-sensed imagery for assessing conifer health in Colorado. Two classes were identified for the purposes of this study: healthy and unhealthy (for example, an area the size of a 30- x 30-m pixel with 20 percent or greater visibly dead trees was defined as ?unhealthy?). Medium-resolution Landsat 5 Thematic Mapper imagery were collected. The normalized, reflectance-converted, cloud-filled Landsat scenes were merged to form a statewide image mosaic, and a Normalized Difference Vegetation Index (NDVI) and Renormalized Difference Infrared Index (RDII) were derived. A supervised maximum likelihood classification was done using the Landsat multispectral bands, the NDVI, the RDII, and 30-m U.S. Geological Survey National Elevation Dataset (NED). The classification was constrained to pixels identified in the updated landcover dataset as coniferous or mixed coniferous/deciduous vegetation. The statewide results were merged with a separate health assessment of Grand County, Colo., produced in late 2008. Sampling and validation was done by collecting field data and high-resolution imagery. The 86 percent overall classification accuracy attained in this study suggests that the data and methods used successfully characterized conifer conditions within Colorado. Although forest conditions for Lodgepole Pine (Pinus contorta) are easily characterized, classification uncertainty exists between healthy/unhealthy Ponderosa Pine (Pinus ponderosa), Pi?on (Pinus edulis), and Juniper (Juniperus sp.) vegetation. Some underestimation of conifer mortality in Summit County is likely, where recent (2008) cloud-free imagery was unavailable. These classification uncertainties are primarily due to the spatial and temporal resolution of Landsat, and of the NLCD derived from this sensor. It is believed that high- to moderate-resolution multispectral imagery, coupled with field data, could significantly reduce the uncertainty rates. The USGS produced a four-county follow-up conifer health assessment using high-resolution RapidEye remotely sensed imagery and field data collected in 2009.

Scientific Investigations Map↗

Estimation of reservoir storage capacity using multibeam sonar and terrestrial lidar, Randy Poynter Lake, Rockdale County, Georgia, 2012

The U.S. Geological Survey, in cooperation with the Rockdale County Department of Water Resources, conducted a bathymetric and topographic survey of Randy Poynter Lake in northern Georgia in 2012. The Randy Poynter Lake watershed drains surface area from Rockdale, Gwinnett, and Walton Counties. The reservoir serves as the water supply for the Conyers-Rockdale Big Haynes Impoundment Authority. The Randy Poynter reservoir was surveyed to prepare a current bathymetric map and determine storage capacities at specified water-surface elevations. Topographic and bathymetric data were collected using a marine-based mobile mapping unit to estimate storage capacity. The marine-based mobile mapping unit operates with several components: multibeam echosounder, singlebeam echosounder, light detection and ranging system, navigation and motion-sensing system, and data acquisition computer. All data were processed and combined to develop a triangulated irregular network, a reservoir capacity table, and a bathymetric contour map.

Georgia↗

Improved Automated Identification and Mapping of Iron Sulfate Minerals, Other Mineral Groups, and Vegetation using Landsat 8 Operational Land Imager Data, San Juan Mountains, Colorado, and Four Corners Region

Multispectral remote sensing data acquired by the Landsat 8 Operational Land Imager (OLI) sensor were analyzed using a new, automated technique to generate a map of exposed mineral and vegetation groups in the western San Juan Mountains, Colo., and the Four Corners Region of the United States. Band ratio results were combined into displayed mineral and vegetation groups using Boolean algebra. New analysis logic has been implemented to exploit the coastal aerosol band in Landsat 8 OLI data and identify concentrations of iron sulfate minerals. These may indicate the presence of near-surface pyrite, which can be a potential nonpoint source of acid rock drainage. Zoned occurrences of iron sulfate minerals in areas surrounding and down gradient of known sources of pyrite have been mapped using this technique. They show high correlation with occurrences of jarosite-bearing mineral assemblages, as mapped using airborne imaging spectrometer data and supporting field verification surveys. Mapping the occurrence of iron sulfate minerals produced by the weathering and oxidation of pyrite in exposed hydrothermally altered rocks can also provide an important indicator of the genetic environment of alteration and the associated mineral deposit type. The automated analysis methodology is being employed to rapidly and cost-effectively generate maps of large regions of the United States in support of U.S. Geological Survey mineral resource and mineral-environmental assessments. This map, which includes the ERDAS IMAGINE thematic raster format in the data release, has been attributed by pixel value with mineral and vegetation group classification data, which can be queried in most image processing and GIS software packages.

Colorado, Utah↗

Development of an approach for integrating components of the U.S. Geological Survey Biomonitoring of Environmental Status and Trends (BEST) and National Stream Quantity Accounting Network (NASQAN) programs for large U.S. rivers

A national-scale framework for monitoring environmental contaminants in fish and effects of contaminant exposure on fish in large U.S. rivers has been proposed by the Biomonitoring of Environmental Status and Trends (BEST) Program of the U.S. Geological Survey (USGS). The framework shares many features and objectives with the USGS National Stream Quantity Accounting Network (NASQAN) Program, which monitors water quality in large U.S. river basins–those with drainage areas of 250,000 to 1,200,000 square miles at their most downstream stations. Because the two programs appear to be complementary, this study was initiated in 2001 to investigate alternative techniques for summarizing and integrating the water-quality data with the fish-contaminant and fish-health data, and to provide recommendations to the BEST program for future integrated studies.

Scientific Investigations Report↗

Summary of significant results from studies of triazine herbicides and their degradation products in surface water, ground water, and precipitation in the midwestern United States during the 1990s

Nonpoint-source contamination of water resources from triazine herbicides has been a major water-quality issue during the 1990s in the United States. To address this issue, studies of surface water, ground water, and precipitation have been carried out by the U.S. Geological Survey in the Midwestern United States. Reconnaissance studies of 147 streams were conducted to determine the geographic and seasonal distribution of atrazine, cyanazine, propazine, and simazine. These studies showed that high concentrations of herbicides were flushed from cropland and transported through the stream system as pulses in response to spring and summer rainfall. The studies also revealed the persistence of herbicides and their degradation products in streams. An investigation of 76 reservoirs showed that the occurrence and temporal distribution of herbicides and their degradation products in reservoir outflow could be related to reservoir and drainage-basin characteristics, water and land use, herbicide use, and climate. Significant findings showed that concentrations of atrazine and its degradation products remained elevated all summer and into the fall and that recently applied atrazine mixed with atrazine applied the previous year as water moved through a reservoir. Reconnaissance studies of 303 ground-water wells were completed to determine hydrogeological and seasonal occurrence, concentration, and distribution of herbicides and their degradation products. Samples collected from across the Midwestern United States consistently revealed that triazine herbicide degradation products commonly were found more frequently than their parent herbicide and that ground-water age could be an important factor in explaining variations in herbicide contamination. A final study investigated precipitation in the Midwestern United States, northeast to the Atlantic Ocean, and northward to the Canadian border. It found that the highest herbicide concentrations in precipitation occurred following herbicide application to cropland. Atrazine was detected most often, followed by deethylatrazine, cyanazine, and deisoproplyatrazine. Mass deposition of herbicides by precipitation was greatest in areas where herbicide use was intense and decreased with distance from the Midwest. Findings of the 1990s studies include an improved understanding of the occurrence, persistence, chemistry, and transport of triazine herbicides and their degradation products in the hydrologic environment. A significant increase in knowledge of triazine herbicides and development and improvement of analytical methods were accomplished in the past decade. The results produced are not only significant for the present (2005) but provide an important data set for future use.

Midwest region↗

Simulation of Water Levels and Salinity in the Rivers and Tidal Marshes in the Vicinity of the Savannah National Wildlife Refuge, Coastal South Carolina and Georgia

The Savannah Harbor is one of the busiest ports on the East Coast of the United States and is located downstream from the Savannah National Wildlife Refuge, which is one of the Nation?s largest freshwater tidal marshes. The Georgia Ports Authority and the U.S. Army Corps of Engineers funded hydrodynamic and ecological studies to evaluate the potential effects of a proposed deepening of Savannah Harbor as part of the Environmental Impact Statement. These studies included a three-dimensional (3D) model of the Savannah River estuary system, which was developed to simulate changes in water levels and salinity in the system in response to geometry changes as a result of the deepening of Savannah Harbor, and a marsh-succession model that predicts plant distribution in the tidal marshes in response to changes in the water-level and salinity conditions in the marsh. Beginning in May 2001, the U.S. Geological Survey entered into cooperative agreements with the Georgia Ports Authority to develop empirical models to simulate the water level and salinity of the rivers and tidal marshes in the vicinity of the Savannah National Wildlife Refuge and to link the 3D hydrodynamic river-estuary model and the marsh-succession model. For the development of these models, many different databases were created that describe the complexity and behaviors of the estuary. The U.S. Geological Survey has maintained a network of continuous streamflow, water-level, and specific-conductance (field measurement to compute salinity) river gages in the study area since the 1980s and a network of water-level and salinity marsh gages in the study area since 1999. The Georgia Ports Authority collected water-level and salinity data during summer 1997 and 1999 and collected continuous water-level and salinity data in the marsh and connecting tidal creeks from 1999 to 2002. Most of the databases comprise time series that differ by variable type, periods of record, measurement frequency, location, and reliability. Understanding freshwater inflows, tidal water levels, and specific conductance in the rivers and marshes is critical to enhancing the predictive capabilities of a successful marsh succession model. Data-mining techniques, including artificial neural network (ANN) models, were applied to address various needs of the ecology study and to integrate the riverine predictions from the 3D model to the marsh-succession model. ANN models were developed to simulate riverine water levels and specific conductance in the vicinity of the tidal marshes for the full range of historical conditions using data from the river gaging networks. ANN models were also developed to simulate the marsh water levels and pore-water salinities using data from the marsh gaging networks. Using the marsh ANN models, the continuous marsh network was hindcasted to be concurrent with the long-term riverine network. The hindcasted data allow ecologists to compute hydrologic parameters?such as hydroperiods and exposure frequency?to help analyze historical vegetation data. To integrate the 3D hydrodynamic model, the marsh-succession model, and various time-series databases, a decision support system (DSS) was developed to support the various needs of regulatory and scientific stakeholders. The DSS required the development of a spreadsheet application that integrates the database, 3D hydrodynamic model output, and ANN riverine and marsh models into a single package that is easy to use and can be readily disseminated. The DSS allows users to evaluate water-level and salinity response for different hydrologic conditions. Savannah River streamflows can be controlled by the user as constant flow, a percentage of historical flows, a percentile daily flow hydrograph, or as a user-specified hydrograph. The DSS can also use output from the 3D model at stream gages near the Savannah National Wildlife Refuge to simulate the effects in the tidal marshes. The DSS is distributed with a two-dimensional (

Georgia, South Carolina↗

Velocities and Attenuations of Gas Hydrate-Bearing Sediments

Monopole and dipole logging data at the Mallik 5L-38, Mackenzie Delta, Canada, provide a challenge for sonic velocity and attenuation models used to remotely estimate pore-space gas hydrate content. Velocity and attenuation are linked, with velocity dispersion causing increased attenuation. Sonic waveforms for Mallik 5L-38, however, show no velocity dispersion in gas hydrate-bearing layers, yet are highly attenuated. Attenuation models applied to Mallik 5L-38 data are shown to be inconsistent with the observed velocity measurements, and therefore are suspect in their ability to predict gas hydrate content. A model explicitly linking velocity and attenuation data is presented, accurately predicting gas hydrate content from velocity data alone while demonstrating that the attenuation mechanisms at the Mallik 5L-38 site have not yet been identified.

Scientific Investigations Report↗

Water quality and ecological condition of urban streams in Independence, Missouri, June 2005 through December 2008

To identify the sources of selected constituents in urban streams and better understand processes affecting water quality and their effects on the ecological condition of urban streams and the Little Blue River in Independence, Missouri the U.S. Geological Survey in cooperation with the City of Independence Water Pollution Control Department initiated a study in June 2005 to characterize water quality and evaluate the ecological condition of streams within Independence. Base-flow and stormflow samples collected from five sites within Independence, from June 2005 to December 2008, were used to characterize the physical, chemical, and biologic effects of storm runoff on the water quality in Independence streams and the Little Blue River. The streams draining Independence-Rock Creek, Sugar Creek, Mill Creek, Fire Prairie Creek, and the Little Blue River-drain to the north and the Missouri River. Two small predominantly urban streams, Crackerneck Creek [12.9-square kilometer (km2) basin] and Spring Branch Creek (25.4-km2 basin), were monitored that enter into the Little Blue River between upstream and downstream monitoring sites. The Little Blue River above the upstream site is regulated by several reservoirs, but streamflow is largely uncontrolled. The Little Blue River Basin encompasses 585 km2 with about 168 km2 or 29 percent of the basin lying within the city limits of Independence. Water-quality samples also were collected for Rock Creek (24.1-km2 basin) that drains the western part of Independence. Data collection included streamflow, physical properties, dissolved oxygen, chloride, metals, nutrients, common organic micro-constituents, and fecal indicator bacteria. Benthic macroinvertebrate community surveys and habitat assessments were conducted to establish a baseline for evaluating the ecological condition and health of streams within Independence. Additional dry-weather screenings during base flow of all streams draining Independence were conducted to identify point-source discharges and other sources of potential contamination. Regression models were used to estimate continuous and annual flow-weighted concentrations, loadings, and yields for chloride, total nitrogen, total phosphorus, suspended sediment, and Escherichia coli bacteria densities. Base-flow and stormflow water-quality samples were collected at five sites within Independence. Base-flow samples for Rock Creek and two tributary streams to the Little Blue River exceeded recommended U.S. Environmental Protection Agency standards for the protection of aquatic life for total nitrogen and total phosphorus in about 90 percent of samples, whereas samples collected at two Little Blue River sites exceeded both the total nitrogen and total phosphorus standards less often, about 30 percent of the time. Dry-weather screening identified a relatively small number (14.0 percent of all analyses) of potential point-source discharges for total chlorine, phenols, and anionic surfactants. Stormflow had larger median measured concentrations of total common organic micro-constituents than base flow. The four categories of common organic micro-constituents with the most total detections in stormflow were pesticides (100 percent), polyaromatic hydrocarbons and combustion by-products (99 percent), plastics (93 percent), and stimulants (91 percent). Most detections of common organic micro-constituents were less than 2 micrograms per liter. Median instantaneous Escherichia coli densities for stormflow samples showed a 21 percent increase measured at the downstream site on the Little Blue River from the sampled upstream site. Using microbial source-tracking methods, less than 30 percent of Escherichia coli bacteria in samples were identified as having human sources. Base-flow and stormflow data were used to develop regression equations with streamflow and continuous water-quality data to estimate daily concentrations, loads, and yields of various water-quality contaminants.

Scientific Investigations Report↗

Effects of urban best management practices on streamflow and phosphorus and suspended-sediment transport on Englesby Brook in Burlington, Vermont, 2000-2010

An assessment of the effectiveness of several urban best management practice structures, including a wet extended detention facility and a shallow marsh wetland (together the "wet extended detention ponds"), was made using data collected from 2000 through 2010 at Englesby Brook in Burlington, Vermont. The purpose of the best management practices was to reduce high streamflows and phosphorus and suspended-sediment loads and concentrations and to increase low streamflows. Englesby Brook was monitored for streamflow, phosphorus, and suspended-sediment concentrations at a streamgage downstream of the best management practice structures for 5 years before the wet extended detention ponds were constructed in 2005 and for 4 years (phosphorus and suspended-sediment concentrations) or 5 years (streamflow) after they were constructed. The period after construction of the best management practice structures was wetter and had higher discharges than the period before construction. Despite the wetter conditions, streamflow duration curves provided evidence that the streamflow regime appeared to have shifted so that the percentages of low streamflows have increased and those of high streamflows may have slightly decreased. Two other hydrologic measures showed improvements in the years following construction of the best management practices: the percentage of annual discharge transported during the 3 days with highest discharges and the number of days with zero streamflow have both decreased. Evidence was mixed for the effectiveness of the best management practices in reducing phosphorus and suspended-sediment concentrations and loads. Annual phosphorus and suspended-sediment loads, monthly loads, low-streamflow concentrations, storm-averaged streamflow-adjusted concentrations, and total storm loads either did not change significantly or increased in the period after construction. These results likely were because of the wetter conditions in the period after construction. For example, monthly loads assessed using analysis of covariance, which compensated for the effects of streamflow on loads, suggested no difference in phosphorus or suspended-sediment loads between the two periods, whereas the comparison of monthly loads without factoring in streamflow showed an increase. This result could be viewed as evidence that the ponds may have mitigated the effect of greater discharges in the period after construction by preventing a corresponding increase in loads. In another analysis used to adjust for the difference in discharge between the two comparison periods, annual and monthly load results were grouped into dry and wet years. Large (50 percent) reductions in annual loads were observed when data from dry (or wet) years before construction were compared with data from dry (or wet) years after construction. When paired monthly loads of each constituent were grouped into dry and wet years, approximately the same number of months had increases as did decreases with the magnitudes of the decreases generally larger than the magnitudes of the increases. These differences in magnitude explain the decrease in annual loads for dry and wet years. The close association of phosphorus with suspended-sediment data suggested that most of the phosphorus was in the particulate form and was controlled by suspended-sediment dynamics.

Vermont↗

Design and testing of a process-based groundwater vulnerability assessment (P-GWAVA) system for predicting concentrations of agrichemicals in groundwater across the United States

Efforts to assess the likelihood of groundwater contamination from surface-derived compounds have spanned more than three decades. Relatively few of these assessments, however, have involved the use of process-based simulations of contaminant transport and fate in the subsurface, or compared the predictions from such models with measured data—especially over regional to national scales. To address this need, a process-based groundwater vulnerability assessment (P-GWAVA) system was constructed to use transport-and-fate simulations to predict the concentration of any surface-derived compound at a specified depth in the vadose zone anywhere in the conterminous United States. The system was then used to simulate the concentrations of selected agrichemicals in the vadose zone beneath agricultural areas in multiple locations across the conterminous United States. The simulated concentrations were compared with measured concentrations of the compounds detected in shallow groundwater (that is, groundwater drawn from within a depth of 6.3 ± 0.5 meters [mean ± 95 percent confidence interval] below the water table) in more than 1,400 locations across the United States. The results from these comparisons were used to select the simulation approaches that led to the closest agreement between the simulated and the measured concentrations. The P-GWAVA system uses computer simulations that account for a broader range of the hydrologic, physical, biological and chemical phenomena known to control the transport and fate of solutes in the subsurface than has been accounted for by any other vulnerability assessment over regional to national scales. Such phenomena include preferential transport and the influences of temperature, soil properties, and depth on the partitioning, transport, and transformation of pesticides in the subsurface. Published methods and detailed soil property data are used to estimate a wide range of model input parameters for each site, including surface albedo, surface crust permeability, soil water content, Brooks-Corey parameters, saturated hydraulic conductivity, macroporosity and sizes of microbial populations, as well as solute partition coefficients, reaction rates, and meso-micropore diffusion rates. To ensure geographic consistency among the predictions, the only site-specific input data that are used are those that are available for all of the 48 conterminous states.

Scientific Investigations Report↗

Continuous monitoring of sediment and nutrients in the Illinois River at Florence, Illinois, 2012-13

The Illinois River is the largest river in Illinois and is the primary contributing watershed for nitrogen, phosphorus, and suspended-sediment loading to the upper Mississippi River from Illinois. In addition to streamflow, the following water-quality constituents were monitored at the Illinois River at Florence, Illinois (U.S. Geological Survey station number 05586300), during May 2012–October 2013: phosphate, nitrate, turbidity, temperature, specific conductance, pH, and dissolved oxygen. The objectives of this monitoring were to (1) determine performance capabilities of the in-situ instruments; (2) collect continuous data that would provide an improved understanding of constituent characteristics during normal, low-, and high-flow periods and during different climatic and land-use seasons; (3) evaluate the ability to use continuous turbidity as a surrogate constituent to determine suspended-sediment concentrations; and (4) evaluate the ability to develop a regression model for total phosphorus using phosphate, turbidity, and other measured parameters. Reliable data collection was achieved, following some initial periods of instrument and data-communication difficulties. The resulting regression models for suspended sediment had coefficient of determination (R 2 ) values of about 0.9. Nitrate plus nitrite loads computed using continuous data were found to be approximately 8 percent larger than loads computed using traditional discrete-sampling based models. A regression model for total phosphorus was developed by using historic orthophosphate data (important during periods of low flow and low concentrations) and historic suspended-sediment data (important during periods of high flow and higher concentrations). The R 2 of the total phosphorus regression model using orthophosphorus and suspended sediment was 0.8. Data collection and refinement of the regression models is ongoing.

Illinois↗

Historical water-quality data from the Harlem River, New York

Data specific to the Harlem River, New York, have been summarized and are presented in this report. The data illustrate improvements in the quality of water for the past 65 years and emphasize the importance of a continuous water-quality record for establishing trends in environmental conditions. Although there is a paucity of sediment-quality data, the New York City Department of Environmental Protection (NYCDEP) Bureau of Wastewater Treatment has maintained a water-quality monitoring network in the Harlem River (and throughout the harbor of New York City) to which 61 combined sewer outfalls discharge effluent. In cooperation with the NYCDEP, the U.S. Geological Survey evaluated water-quality data collected by the NYCDEP dating back to 1945, which indicate trends in water quality and reveal improvement following the 1972 passage of the Clean Water Act. These improvements are indicated by the steady increase in median dissolved oxygen concentrations and an overall decrease in fecal indicator bacteria concentrations starting in the late 1970s. Further, the magnitude of the highest fecal indicator bacteria concentrations (that is, the 90th percentile) in samples collected from the Harlem River have decreased significantly over the past four decades. Other parameters of water quality used to gauge the health of a water body include total suspended solids and nutrient (inorganic forms of nitrogen and phosphorus) concentrations—mean concentrations for these indicators have also decreased in the past decades. The limited sediment data available for one sample in the Harlem River indicate concentrations of copper, zinc, and lead are above sediment-quality thresholds set by the New York State Department of Environmental Conservation. However, more data are needed to better understand the changes in both sediment and water quality in the Harlem River, both as the tide cycles and during precipitation events. As a partner in the Urban Waters Federal Partnership, the U.S. Geological Survey has worked to address the chronic water-quality concerns of the Harlem River by compiling relevant data and studies, which is an important component for understanding and rectifying water-quality problems within a watershed.

New York↗

Performance evaluation testing of wells in the gradient control system at a federally operated Confined Disposal Facility using single well aquifer tests, East Chicago, Indiana

The U.S. Geological Survey (USGS) performed tests to evaluate the hydrologic connection between the open interval of the well and the surrounding Calumet aquifer in response to fouling of extraction well pumps onsite. Two rounds of air slug testing were performed on seven monitoring wells and step drawdown and subsequent recovery tests on three extraction wells on a U.S. Army Corps of Engineers Confined Disposal Facility (CDF) in East Chicago, Indiana. The wells were tested in 2014 and again in 2015. The extraction and monitoring wells are part of the gradient control system that establishes an inward gradient around the perimeter of the facility. The testing established a set of protocols that site personnel can use to evaluate onsite well integrity and develop a maintenance procedure to evaluate future well performance. The results of the slug test analysis data indicate that the hydraulic connection of the well screen to the surrounding aquifer material in monitoring wells on the CDF and the reliability of hydraulic conductivity estimates of the surrounding geologic media could be increased by implementing well development maintenance. Repeated air slug tests showed increasing hydraulic conductivity until, in the case of the monitoring wells located outside of the groundwater cutoff wall (MW–4B, MW–11B, MW–14B), the difference in hydraulic conductivity from test to test decreased, indicating the results were approaching the optimal hydraulic connection between the aquifer and the well screen. Hydraulic conductivity values derived from successive tests in monitoring well D40, approximately 0.25 mile south of the CDF, were substantially higher than those derived from wells on the CDF property. Also, values did not vary from test to test like those measured in monitoring wells located on the CDF property, which indicated that a process may be affecting the connectivity of the wells on the CDF property to the Calumet aquifer. Derived hydraulic conductivity values from the initial air slug test during the 2015 testing period for MW–11A and MW–14A are an order of magnitude less than those derived from the final test during the 2014 testing period indicating the development of a low conductivity skin between the final test of the 2014 testing period and the beginning of the 2015 testing period that created a decrease in the connection of the monitoring well screen to the surrounding aquifer material. Repeated step drawdown and recovery testing of the extraction wells tested during this study provided results that indicate a slight increase in the development of a skin and a decrease in the connectivity of the extraction wells with the Calumet aquifer. Hydraulic conductivity values obtained from the test results were relatively similar in EW–4B and EW–14A but were substantially lower for EW–11C. This difference may be due to the presence of finer grained silt deposits in the area surrounding well nest 11. Skin factors calculated during the step drawdown and recovery analysis were lowest in EW–11C and relatively similar in EW–4B and EW–14A. Calculated skin factors increased slightly in the analysis of data collected in 2015 from that collected in 2014. Comparisons of the specific-capacity values calculated from well development data collected following extraction well installation to those calculated during the single well aquifer tests at EW–4B, EW–14A and EW–11C indicate that the productivity of extraction wells on the CDF property has diminished since 2008. Values calculated for monitoring wells MW–4A, MW–11A, and MW–14A were used to evaluate the decrease in air slug derived hydraulic conductivity for monitoring wells within the groundwater cutoff wall between testing in 2014 and 2015. Results from testing by this study indicate that implementation of an air slug testing regimen of the monitoring wells that control the gradient control system at the CDF throughout the course of a year may help sustain the connectivity between the monitoring wells and the surrounding aquifer and provide data to evaluate the need for different types of well development approaches to address chemical or biological fouling issues. Repeated step drawdown and recovery testing of the extraction wells tested during this study provided results that indicate a slight increase in the development of a skin and a decrease in the connectivity of the extraction wells with the Calumet aquifer. Implementation of a specific capacity testing regimen can provide data to record and track well condition through time for individual extraction wells. Results from aquifer testing by this study indicate that specific capacity test results, when paired with recovery testing, provide useful data to measure the development of any low conductivity wellbore skin through the skin factors derived for the individual extraction wells. An initial annual schedule of specific capacity and recovery tests would provide sufficient data to identify substantial short-term changes in the operating condition of the extraction wells.

Indiana↗

Simulated groundwater flow paths, travel time, and advective transport of nitrogen in the Kirkwood-Cohansey aquifer system, Barnegat Bay–Little Egg Harbor Watershed, New Jersey

Elevated concentrations of nitrogen in groundwater that discharges to surface-water bodies can degrade surface-water quality and habitats in the New Jersey Coastal Plain. An analysis of groundwater flow in the Kirkwood-Cohansey aquifer system and deeper confined aquifers that underlie the Barnegat Bay–Little Egg Harbor (BB-LEH) watershed and estuary was conducted by using groundwater-flow simulation, in conjunction with a particle-tracking routine, to provide estimates of groundwater flow paths and travel times to streams and the BB-LEH estuary. Water-quality data from the Ambient Groundwater Quality Monitoring Network, a long-term monitoring network of wells distributed throughout New Jersey, were used to estimate the initial nitrogen concentration in recharge for five different land-use classes—agricultural cropland or pasture, agricultural orchard or vineyard, urban non-residential, urban residential, and undeveloped. Land use at the point of recharge within the watershed was determined using a geographic information system (GIS). Flow path starting locations were plotted on land-use maps for 1930, 1973, 1986, 1997, and 2002. Information on the land use at the time and location of recharge, time of travel to the discharge location, and the point of discharge were determined for each simulated flow path. Particle-tracking analysis provided the link from the point of recharge, along the particle flow path, to the point of discharge, and the particle travel time. The travel time of each simulated particle established the recharge year. Land use during the year of recharge was used to define the nitrogen concentration associated with each flow path. The recharge-weighted average nitrogen concentration for all flow paths that discharge to the Toms River upstream from streamflow-gaging station 01408500 or to the BB-LEH estuary was calculated. Groundwater input into the Barnegat Bay–Little Egg Harbor estuary from two main sources— indirect discharge from base flow to streams that eventually flow into the bay and groundwater discharge directly into the estuary and adjoining coastal wetlands— is summarized by quantity, travel time, and estimated nitrogen concentration. Simulated average groundwater discharge to streams in the watershed that flow into the BB-LEH estuary is approximately 400 million gallons per day. Particle-tracking results indicate that the travel time of 56 percent of this discharge is less than 7 years. Fourteen percent of the groundwater discharge to the streams in the BB-LEH watershed has a travel time of less than 7 years and originates in urban land. Analysis of flow-path simulations indicate that approximately 13 percent of the total groundwater flow through the study area discharges directly to the estuary and adjoining coastal wetlands (approximately 64 million gallons per day). The travel time of 19 percent of this discharge is less than 7 years. Ten percent of this discharge (1 percent of the total groundwater flow through the study area) originates in urban areas and has a travel time of less than 7 years. Groundwater that discharges to the streams that flow into the BB-LEH, in general, has shorter travel times, and a higher percentage of it originates in urban areas than does direct groundwater discharge to the Barnegat Bay–Little Egg Harbor estuary. The simulated average nitrogen concentration in groundwater that discharges to the Toms River, upstream from streamflow-gaging station 01408500 was computed and compared to summary concentrations determined from analysis of multiple surface-water samples. The nitrogen concentration in groundwater that discharges directly to the estuary and adjoining coastal wetlands is a current data gap. The particle tracking methodology used in this study provides an estimate of this concentration."

New Jersey↗

Magnitude and frequency of floods in Alabama, 2015

To improve flood-frequency estimates at rural streams in Alabama, annual exceedance probability flows at gaged locations and regional regression equations used to estimate annual exceedance probability flows at ungaged locations were developed by using current geospatial data, new analytical methods, and annual peak-flow data through September 2015 at 242 streamgages in Alabama and surrounding States. The regional regression equations were derived from statistical analyses of annual peak-flow data and basin characteristics for a subset of 217 streamgages. Four flood regions were identified based on residuals from the regional regression analyses and contain sites with similar basin characteristics. A separate set of equations was derived for estimating flood frequency and magnitude for small rural streams using a subset of 40 small basin streamgages. A large river analysis was also completed for 14 selected large-river streamgages in Alabama. Annual exceedance probability flows presented in this report reflect additional streamflow data collected since the previous study of flood magnitude and frequency in Alabama, which included streamflow through September 2003.

Alabama↗

Groundwater/surface-water interactions in the Partridge River Basin and evaluation of hypothetical future mine pits, Minnesota

The Partridge River Basin (PRB) covers 156 square miles in northeastern Minnesota with headwaters in the Mesabi Iron Range. The basin is characterized by extensive wetlands, lakes, and streams in poorly drained and often thin glacial material overlying Proterozoic bedrock. To better understand the interaction between these extensive surface water features and the groundwater system, a three-dimensional, steady-state, groundwater-flow model of the PRB was developed by the U.S. Geological Survey in cooperation with the Great Lakes Indian Fish & Wildlife Commission using the finite-difference computer code MODFLOW-NWT. The model simulates steady-state base flow in streams and groundwater interactions using the streamflow routing (SFR2) package. Existing mining features including tailings basins, stockpiles, pumped mine pits, and flooded mine pits were simulated using either high hydraulic conductivity zones or the drain (DRN) package. The unsaturated zone flow (UZF) package was used to better represent the groundwater system in areas with a high water table and for wetlands often associated with such areas. UZF typically is used to represent unsaturated zone processes but also can simulate the rejection of recharge and groundwater discharge to the land surface when the water table is near land surface. The steady-state model used data from the 2011 to 2013 period when 2011 high-resolution land surface (light detecting and ranging [lidar]) data were available that reflected land-surface and water elevations from mining activity in the basin. The parameter-estimation software suite PEST_HP was used to obtain a best fit of the modeled to measured groundwater levels, streamflow, pit inflow rates, and mapped peat deposits. The PEST calibration used the target residuals from two models with the same model parameters and targets from two separate periods: (1) a 1995–2015 calibration model, which provided a larger number of calibration targets, and (2) a 2011–2013 mining conditions model, which included calibration targets that reflected conditions consistent with the modeled mine-workings topography. Calibration of the PRB model resulted in ranges of glacial horizontal hydraulic conductivity parameters that generally agreed with literature values and other models of the region. Horizontal hydraulic conductivity of the bedrock was higher in the upper bedrock layers where numerous and continuous fractures have been observed and lower in the deeper bedrock layers. Average basin-wide calibrated infiltration was 5.3 inches per year. An average of 4.6 inches per year of infiltration crosses the water table and becomes recharge and 0.7 inch per year is rejected by UZF due to saturated conditions at the land surface. Simulated groundwater runoff (the sum of rejected recharge and groundwater seepage to the land surface) can either be routed to streams or removed from the model as evapotranspiration. The calibrated model indicates relatively shallow groundwater-flow paths dominating and approximately 50 percent of the stream base flow coming from groundwater runoff. The 2011–2013 mining conditions model was then used to develop five model scenarios simulating the response of the groundwater and surface-water system to potential hydrologic stress. The purpose of these mine pit scenarios is to present a possible workflow to quantify a model’s uncertainty for a given model forecast and serve as a possible guide for initial data collection that may improve a future model’s ability to make such a forecast. The scenarios included one scenario with the currently existing Peter Mitchell pit at final buildout and flooded to an elevation of 1,500 feet, and four scenarios with a hypothetical, new mine pit plus the flooded Peter Mitchell at final buildout. The five model scenarios were used to forecast streamflow at six locations in the PRB, pit inflow rates for the new mine pits and the flooded Peter Mitchell pit, and the average depth to water in 12 wetlands. A linear uncertainty analysis was performed using information from the PEST calibration and tools in the PyEMU python package to assess model uncertainty propagation to the model forecasts. Streamflows generally were reduced with future mining and the greatest streamflow reductions occurred from the flooded Peter Mitchell Pit, probably due to its large size. Average depth to groundwater in wetlands was most affected the closer the wetland was to a new mine pit. Linear uncertainty methods were also used to evaluate data worth, which is the ability for potential new groundwater elevation observations to reduce the uncertainty in scenario forecasts. Data worth was performed for a grid of new hydraulic head observations. Overall, areas with nonnegligible data worth generally corresponded to wetland areas with no groundwater seepage to land surface from UZF. These model behaviors indicated that the land-surface boundary condition simulated by the UZF package was pinning the groundwater elevations to the land surface in areas with groundwater seepage (33 percent of the 2011–2013 base conditions model) such that the sensitivity to new observations in these areas was minimal. Therefore, representing wetlands as boundary conditions minimized the usefulness of data worth calculations because wetland areas were present over a large part of the model domain. Probabilistic capture zones were estimated for each of the mines in the model scenarios. A capture zone represents the area contributing recharge to a model feature, like a well or a mine pit, and can be calculated by forward tracking particles from the water table. By using Monte Carlo techniques, it is possible to generate estimated capture zones that include the probability of recharge capture given the uncertainty present in the model. Monte Carlo techniques use randomly generated model parameter sets sampled from a plausible parameter range to create many possible realizations. The resulting capture zone arrays were calculated by tallying the total number of realizations in which a particle from a model cell was captured by the feature. Probabilities from the Monte Carlo runs ranged from 1 (captured in 100 percent of the runs) near the pits to 0 (captured in 0 percent of the runs) at the edges of the capture zone. Capture zones were not always spatially continuous; for example, the capture zone for the proposed mine pits south of the flooded Peter Mitchell pit was discontinuous with capture surrounding the proposed mine pit and north of the flooded Peter Mitchell pit. This northern section represents deeper groundwater flow paths that originate in the topographic high, move under the flooded pit, and discharge into the proposed pit. This pattern of capture indicates the possibility of some deeper flow through the upper fractured bedrock when the shallow groundwater flow system is modified. These results underscore that future site-specific applications of the base condition model require the input of site-specific data and recalibration to focus on the site of interest.

Minnesota↗