Geology Reports⌕ Search

SEARCH · Geology Reports

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,747 records · Page 97Linked to original sources

Relations among water levels, specific conductance, and depths of bedrock fractures in four road-salt-contaminated wells in Maine, 2007–9

Data on groundwater-level, specific conductance (a surrogate for chloride), and temperature were collected continuously from 2007 through 2009 at four bedrock wells known to be affected by road salts in an effort to determine the effects of road salting and fractures in bedrock that intersect the well at a depth below the casing on the presence of chloride in groundwater. Dissolved-oxygen data collected periodically also were used to make inferences about the interaction of fractures and groundwater flow. Borehole geophysical tools were used to determine the depths of fractures in each well that were actively contributing flow to the well, under both static and pumped conditions; sample- and measurement-depths were selected to correspond to the depths of these active fractures. Samples of water from the wells, collected at depths corresponding to active bedrock fractures, were analyzed for chloride concentration and specific conductance; from these analyses, a linear relation between chloride concentration and specific conductance was established, and continuous and periodic measurements of specific conductance were assumed to represent chloride concentration of the well water at the depth of measurement. To varying degrees, specific conductance increased in at least two of the wells during winter and spring thaws; the shallowest well, which also was closest to the road receiving salt treatment during the winter, exhibited the largest changes in specific conductance during thaws. Recharge events during summer months, long after application of road salt had ceased for the year, also produced increases in specific conductance in some of the wells, indicating that chloride which had accumulated or sequestered in the overburden was transported to the wells throughout the year. Geophysical data and periodic profiles of water quality along the length of each well’s borehole indicated that the greatest changes in water quality were associated with active fractures; in one case, high concentration of dissolved oxygen at the bottom of the well indicated the presence of a highly transmissive fracture that was in good connection with a surficial feature (stream or atmosphere). Data indicated that fractures have a substantial influence on the transport of chlorides to the subsurface; that elevated specific conductance occurred throughout the year, not just when road salts were applied; and that chloride contamination, as indicated by elevated specific conductance, may persist for years.

Maine↗

Flood-inundation maps for the South Platte River at Fort Morgan, Colorado, 2018

In 2017, the U.S. Geological Survey (USGS), in cooperation with the Colorado Water Conservation Board (CWCB), studied floods in the historic record to produce a library of flood-inundation maps for the South Platte River at Fort Morgan, Colorado. Digital flood-inundation maps for a 4.5-mile (7.2-kilometers) reach of the South Platte River at Fort Morgan from Morgan County Road 16 to Morgan County Road 20.5, were created. The flood-inundation maps depict estimates of the areal extent and depth of flooding corresponding to selected water levels (stages) at the U.S. Geological Survey streamgage on the South Platte River at Fort Morgan (streamgage number 06759500). Water-surface profiles were computed for the stream reach by means of a one-dimensional, step-backwater model. The September 15, 2013, and May 20, 2017, floods were used to calibrate the model, and the June 15, 2015, and May 29, 2017, floods were used to independently validate the model. Nine pressure transducers were deployed to record the stage at nine different locations along the reach and to document the floods of May 20 and 29, 2017, at the South Platte River at Fort Morgan streamgage. The hydraulic model was then used to compute 16 water-surface profiles for flood stages at 1-foot (ft; 0.3-meter [m]) intervals referenced to the streamgage datum and ranging from 12 ft (3.66 m) or below bankfull to 27 ft (8.23 m), which is 1 ft (0.3 m) greater than the highest recorded water level (25.73 ft [7.84 m] on September 15, 2013) at the South Platte River at Fort Morgan streamgage during its period of record; the 2013 flood exceeds the major flood stage of 21.5 ft (6.55 m) by more than 4 ft (1.2 m), as defined by the National Weather Service. The simulated water-surface profiles were then combined with a geographic information system digital elevation model (derived from light detection and ranging data having a 0.37-ft [0.11-m] vertical accuracy and 3.28-ft [1.00-m] horizontal resolution) to delineate the area flooded for stages ranging from 12 to 27 ft (3.66 to 8.23 m). These flood-inundation maps, in conjunction with the real-time stage data from the USGS streamgage on the South Platte River at Fort Morgan, are intended to help guide the general public in taking individual safety precautions and provide emergency management personnel with a tool to efficiently manage emergency flood operations and post flood recovery efforts.

Colorado↗

California’s exposure to volcanic hazards

The potential for damaging earthquakes, landslides, floods, tsunamis, and wildfires is widely recognized in California. The same cannot be said for volcanic eruptions, despite the fact that they occur in the state about as frequently as the largest earthquakes on the San Andreas Fault. At least ten eruptions have taken place in the past 1,000 years, and future volcanic eruptions are inevitable. The U.S. Geological Survey’s (USGS) national volcanic threat assessment identifies eight young volcanic areas in California as moderate, high, or very high threat. Of the eight volcanic areas that exist in California, molten rock resides beneath at least seven of these—Medicine Lake volcano, Mount Shasta, Lassen Volcanic Center, Clear Lake volcanic field, the Long Valley volcanic region, Coso volcanic field, and Salton Buttes—and are therefore considered “active” volcanoes producing volcanic earthquakes, toxic gas emissions, hot springs, geothermal systems, and (or) ground movement. The USGS California Volcano Observatory in Menlo Park, California, monitors these potentially hazardous volcanoes to help communities and government authorities understand, prepare for, and respond to volcanic activity. Although volcanic activity can sometimes be forecast, eruptions, like earthquakes or tsunamis, cannot be prevented. Understanding the hazards and identifying what and who is in harm’s way is the first step in mitigating volcanic risk and building community resilience to volcanic hazards. This report, which was prepared in collaboration with the California Governor’s Office of Emergency Services and the California Geological Survey, provides a broad perspective on the state’s exposure to volcanic hazards by integrating volcanic hazard information with geospatial data on at-risk populations, infrastructure, and resources. This information is intended to prompt site- and sector-specific vulnerability analyses and preparation of hazard mitigation and response plans.

California↗

Geochemical assessment of groundwater in the Big Chino subbasin, Arizona, 2011–18

A geochemical characterization of groundwater in the Big Chino subbasin of Arizona was conducted by the U.S. Geological Survey, in cooperation with the City of Prescott, the Town of Prescott Valley, and the Salt River Project, to understand groundwater evolution through the study area and the source of water to springs along the gaining reach of the Verde River just downstream from its confluence with Granite Creek. Samples were collected between 2011 and 2018 in groundwater wells completed in basin-fill and carbonate aquifers and at selected springs, including two discrete springs discharging along the aforementioned stretch of the Verde River. Five newly installed monitoring wells completed in the carbonate aquifer were sampled in 2018. Water-quality results obtained from these samples include the first known geochemical data for carbonate groundwater beneath the basin-fill in the Big Chino subbasin downgradient from Walnut Creek near Paulden, Arizona, as well as other parts of the study area without previous data. Groundwater samples were collected and analyzed for major ions, arsenic, nutrients, stable isotopes of oxygen and hydrogen (δ 18 O and δ 2 H), strontium isotopes ( 87 Sr/ 86 Sr), carbon-14, isotopes of carbon (δ 13 C), and noble gases. Significant differences in groundwater geochemistry between the basin-fill and carbonate aquifers were driven primarily by higher pH, tritium, and δ 18 O and δ 2 H in the basin-fill aquifer samples and higher specific conductance and higher concentrations of calcium, sodium, bicarbonate, fluoride, and arsenic in the carbonate aquifer samples. All but one sample from the carbonate aquifer and two samples from the basin-fill aquifer exceeded the U.S. Environmental Protection Agency (EPA) drinking water standard for arsenic of 10 micrograms per liter. One basin-fill aquifer sample exceeded the EPA drinking water standard for fluoride of 4 milligrams per liter, and one carbonate aquifer sample exceeded the EPA secondary drinking water standard for fluoride of 2 milligrams per liter. A component of modern groundwater recharged following aboveground nuclear testing beginning in the mid-1950s is present in some basin-fill and spring groundwater from this study. Groundwater that can be dated using radiocarbon decay is also present in the study area, with four groundwater samples indicating possible recharge during the Pleistocene with groundwater ages ranging from approximately 34,600 to 13,300 years before present. Other groundwater sampled during this study that can dated using radiocarbon decay ranged in age from about 7,500 to 1,100 years before present, indicating possible recharge during the Holocene. The gaining reach of the Verde River downstream from the confluence with Granite Creek shows areal changes in temperature, pH, and specific conductance, indicating multiple zones of groundwater input. Surface-water samples for analyses of δ 18 O and δ 2 H have been collected at the Verde River near Paulden, Ariz. streamgage (09503700) during discharge measurements since 2009, and a trend analysis of the δ 18 O and δ 2 H data indicated no significant trend exists for the 10-year period of record. Additional groundwater samples from the carbonate aquifer beneath the basin-fill upgradient and downgradient from Walnut Creek would provide valuable information to understand groundwater evolution along the Big Chino subbasin.

Arizona↗

Flood-inundation maps for the Cuyahoga River at Jaite, Ohio, 2024

Digital flood-inundation maps for a nearly 6-mile reach of the Cuyahoga River at Jaite, Ohio, were created by the U.S. Geological Survey (USGS) in cooperation with the Northeast Ohio Regional Sewer District Board of Trustees. The maps depict estimates of the extent and depth of flooding corresponding to selected water levels (stages) at USGS streamgage 04206425 on the Cuyahoga River at Jaite, Ohio. Water-surface profiles were computed for the stream reach by using a one-dimensional steady-state step-backwater model. The hydraulic model was calibrated to the current USGS streamgage data and then used to compute 15 water-surface profiles for flood stages at 1-foot intervals referenced to the streamgage datum and ranging from 6 to 20 feet, which correspond to below “action stage” to “major flood stage” as reported by the National Weather Service. The simulated water-surface profiles were then used with a geographic information system digital elevation model derived from light detection and ranging data to delineate the areas flooded at each stage. These maps, along with current stage data from the USGS streamgage and forecasted high-flow stages from the National Weather Service, can provide emergency management personnel and residents with information that is critical for flood response activities such as evacuations and road closures, as well as for postflood recovery efforts.

Ohio↗

Scientist’s guide to developing explanatory statistical models using causal analysis principles

Recent discussions of model selection and multimodel inference highlight a general challenge for researchers, which is how to clearly convey the explanatory content of a hypothesized model or set of competing models. The advice from statisticians for scientists employing multimodel inference is to develop a well‐thought‐out set of candidate models for comparison, though precise instructions for how to do that are typically not given. A coherent body of knowledge, which falls under the general term causal analysis, now exists for examining the explanatory scientific content of candidate models. Much of the literature on causal analysis has been recently developed and we suspect may not be familiar to many ecologists. This body of knowledge comprises a set of graphical tools and axiomatic principles to support scientists in their endeavors to create “well‐formed hypotheses”, as statisticians are asking them to do. Causal analysis is complementary to methods such as structural equation modeling, which provides the means for evaluation of proposed hypotheses against data. In this paper, we summarize and illustrate a set of principles that can guide scientists in their quest to develop explanatory hypotheses for evaluation. The principles presented in this paper have the capacity to close the communication gap between statisticians, who urge scientists to develop well‐thought‐out coherent models, and scientists, who would like some practical advice for exactly how to do that.

Ecology↗

TriNet "ShakeMaps": Rapid generation of peak ground motion and intensity maps for earthquakes in southern California

Rapid (3-5 minutes) generation of maps of instrumental ground-motion and shaking intensity is accomplished through advances in real-time seismographic data acquisition combined with newly developed relationships between recorded ground-motion parameters and expected shaking intensity values. Estimation of shaking over the entire regional extent of southern California is obtained by the spatial interpolation of the measured ground motions with geologically based frequency and amplitude-dependent site corrections. Production of the maps is automatic, triggered by any significant earthquake in southern California. Maps are now made available within several minutes of the earthquake for public and scientific consumption via the World Wide Web; they will be made available with dedicated communications for emergency response agencies and critical users.

California↗

Data report: Stable isotopic measurements of sedimentary organic matter and N. pachyderma (s.) from site 1166, Prydz Bay continental shelf

We report the results of downhole stable isotopic ( ẟ 13 C org [organic carbon] and ẟ 15 N) and elemental measurements (total organic carbon [TOC], total nitrogen [TN], and carbon/nitrogen [C/N]) of sedimentary organic matter (SOM) along with stable isotopic measurements ( ẟ 18 O and ẟ 13 C) of left-coiling Neogloboquadrina pachyderma planktonic foraminifers from Ocean Drilling Program Site 1166. TOC and TN measurements indicate a large change from organic-rich preglacial sediments with primary organic matter to organic-poor early glacial and glacial sediments, with mainly recycled organic matter. Results of the stable isotopic measurements of SOM show a range of values that are typical of both marine and terrestrial organic matter, probably reflecting a mixture of the two. However, C/N values are mostly high (>15), suggesting greater input and/or preservation of terrestrial organic matter. Foraminifers are only present in glacial/glaciomarine sediments of latest Pliocene to Pleistocene age at Site 1166 (lithostratigraphic Unit I). The majority of this unit has ẟ 13 C org and TOC values that are similar to those of glacial sediments recovered at Site 1167 (lithostratigraphic Unit II) on the slope and may have the same source(s). Although the low resolution of the N. pachyderma (s.) ẟ 18 O and ẟ 13 C data set precludes any specific paleoclimatic interpretation, downcore variations in foraminifer ẟ 18 O and ẟ 13 C values of 0.5 to 1 amplitude may indicate glacial-interglacial changes in ice volume/temperature in the Prydz Bay region.

Proceedings of the Ocean Drilling Program: Scienti↗

Precambrian crystalline basement map of Idaho: An interpretation of aeromagnetic anomalies

Idaho lies within the northern sector of the U.S. Cordillera astride the boundary between the Proterozoic continent (Laurentia) to the east and the Permian to Jurassic accreted terranes to the west. The continental basement is mostly covered by relatively undeformed Mesoproterozoic metasedimentary rocks and intruded or covered by Phanerozoic igneous rocks; accordingly, knowledge of the basement geology is poorly constrained. Incremental knowledge gained since the pioneering studies by W. Lindgren, C.P. Ross, A.L. Anderson, A. Hietanen, and others during the early- and mid-1900's has greatly advanced our understanding of the general geology of Idaho. However, knowledge of the basement geology remains relatively poor, partly because of the remoteness of much of the region plus the lack of a stimulus to decipher the complex assemblage of high-grade gneisses and migmatite of central Idaho. The availability of an updated aeromagnetic anomaly map of Idaho (North American Magnetic Anomaly Group, 2002) provides a means to determine the regional Precambrian geologic framework of the State. The combined geologic and aeromagnetic data permit identification of previously unrecognized crystalline basement terranes, assigned to Archean and Paleoproterozoic ages, and the delineation of major shear zones, which are expressed in the aeromagnetic data as linear negative anomalies (Finn and Sims, 2004). Limited geochronologic data on exposed crystalline basement aided by isotopic studies of zircon inheritance, particularly Bickford and others (1981) and Mueller and others (1995), provide much of the geologic background for our interpretation of the basement geology. In northwestern United States, inhomogeneities in the basement inherited from Precambrian tectogenesis controlled many large-scale tectonic features that developed during the Phanerozoic. Two basement structures, in particular, provided zones of weakness that were repeatedly rejuvenated: (1) northeast-trending ductile shear zones developed on the northwest margin of the Archean Wyoming province during the Paleoproterozoic Trans-Montana orogeny (Sims and others, 2004), and (2) northwest-trending intra-continental faults of the Mesoproterozoic Trans-Rocky Mountain strike-slip fault system (Sims, unpub. data, 2003). In this report, geologic ages are reported in millions of years (Ma) and generalized ages are given in billions of years (Ga). The subdivision of Precambrian rocks used herein is the time classification recommended by the International Union of Geological Sciences (Plumb, 1991).

Idaho↗

Flood-inundation maps for the Saluda River from Old Easley Bridge Road to Saluda Lake Dam near Greenville, South Carolina

Digital flood-inundation maps for a 3.95-mile reach of the Saluda River from approximately 815 feet downstream from Old Easley Bridge Road to approximately 150 feet downstream from Saluda Lake Dam near Greenville, South Carolina, were developed by the U.S. Geological Survey (USGS). The inundation maps, which can be accessed through the USGS Flood Inundation Mapping Science Web site at http://water.usgs.gov/osw/flood_inundation/, depict estimates of the areal extent and depth of flooding corresponding to selected water levels (stages) at the USGS streamgage at Saluda River near Greenville, South Carolina (station 02162500). Current conditions at the USGS streamgage may be obtained through the National Water Information System Web site at http://waterdata.usgs.gov/sc/nwis/uv/?site_no=02162500&PARAmeter_cd=00065,00060,00062. The National Weather Service (NWS) forecasts flood hydrographs at many places that are often collocated with USGS streamgages. Forecasted peak-stage information is available on the Internet at the NWS Advanced Hydrologic Prediction Service (AHPS) flood-warning system Web site (http://water.weather.gov/ahps/) and may be used in conjunction with the maps developed in this study to show predicted areas of flood inundation.In this study, flood profiles were computed for the stream reach by means of a one-dimensional step-backwater model. The model was calibrated using the most current stage-streamflow relations at USGS streamgage station 02162500, Saluda River near Greenville, South Carolina. The hydraulic model was then used to determine water-surface profiles for flood stages at 1.0-foot intervals referenced to the streamgage datum and ranging from approximately bankfull to 2 feet higher than the highest recorded water level at the streamgage. The simulated water-surface profiles were then exported to a geographic information system, ArcGIS, and combined with a digital elevation model (derived from Light Detection and Ranging [LiDAR] data with a 0.6-foot vertical Root Mean Square Error [RMSE] and a 3.0-foot horizontal RMSE), using HEC-GeoRAS tools in order to delineate the area flooded at each water level. The availability of these maps, along with real-time stage data from the USGS streamgage station 02162500 and forecasted stream stages from the NWS, can provide emergency management personnel and residents with information that is critical during flood-response and flood-recovery activities, such as evacuations, road closures, and disaster declarations.

South Carolina↗

Flood-inundation maps for the Saddle River from Upper Saddle River Borough to Saddle River Borough, New Jersey, 2013

Digital flood-inundation maps for a 4.1-mile reach of the Saddle River from 0.6 miles downstream from the New Jersey-New York State boundary in Upper Saddle River Borough to 0.2 miles downstream from the East Allendale Road bridge in Saddle River Borough, New Jersey, were created by the U.S. Geological Survey (USGS) in cooperation with the New Jersey Department of Environmental Protection (NJDEP). The inundation maps, which can be accessed through the USGS Flood Inundation Mapping Science Web site at http://water.usgs.gov/osw/flood_inundation/, depict estimates of the areal extent and depth of flooding corresponding to select water levels (stages) at the USGS streamgage 01390450, Saddle River at Upper Saddle River, New Jersey. Current conditions for estimating near real-time areas of inundation using USGS streamgage information may be obtained on the Internet at http://waterdata.usgs.gov/nwis/uv?site_no=01390450. The National Weather Service (NWS) forecasts flood hydrographs at many places that are often collocated with USGS streamgages. NWS-forecasted peak-stage information may be used in conjunction with the maps developed in this study to show predicted areas of flood inundation. In this study, flood profiles were computed for the stream reach by means of a one-dimensional step-backwater model. The model was calibrated by using the most current stage-discharge relations (in effect March 2013) at USGS streamgage 01390450, Saddle River at Upper Saddle River, New Jersey, and documented high-water marks from recent floods. The hydraulic model was then used to determine eight water-surface profiles for flood stages at 0.5-foot (ft) intervals referenced to the streamgage datum, North American Vertical Datum of 1988 (NAVD 88), and ranging from bankfull, 0.5 ft below NWS Action Stage, to the upper extent of the stage-discharge rating which is approximately 1 ft higher than the highest recorded water level at the streamgage. Action Stage is the stage which when reached by a rising stream the NWS or a partner needs to take some type of mitigation action in preparation for possible significant hydrologic activity. The simulated water-surface profiles were then combined with a geographic information system 3-meter (9.84 ft) digital elevation model (derived from Light Detection and Ranging (LiDAR) data) in order to delineate the area flooded at each water level. The availability of these maps along with real-time streamflow data and information regarding current stage from USGS streamgages and forecasted stream stages from the NWS provide emergency management personnel and residents with information that is critical for flood response activities, such as evacuations and road closures, as well as for post-flood recovery efforts.

New Jersey↗

Geologic map of the north half of the Lake Walcott 30'×60' quadrangle, Idaho

The geologic map of the northern half of the Lake Walcott 30ʹ×60ʹ quadrangle shows the volcanic geology of the southern part of the Craters of the Moon lava field, the complex geologic features of the Holocene Kings Bowl and Wapi lava fields, and the southern part of the Great Rift volcanic rift zone. The long extent and distribution of skylights in lava-tube systems of the Horse Butte and Wapi Park lava fields are depicted on this map. 40 Ar/ 39 Ar and K/Ar age determinations give detail to the Holocene, late Pleistocene, and late middle Pleistocene volcanic lava fields in this quadrangle. Most of the younger basalt eruptions (less than 150 thousand years [ka]) have occurred along the Great Rift volcanic rift zone, but two of the younger lava fields are located in the western part of the quadrangle. Kimama Butte, a shield volcano, is 87±11 ka, and Shale Butte is dated at 11±6 ka. Paleomagnetic studies have shown that the Horse Butte-Inferno Chasm eruptive fissure system has at least five paleomagnetic-correlative lava fields, the Claasen vent complex consists of at least seven correlative lava fields, and the Streifling-Flat Top vent complex includes at least four correlative lava fields. This map provides geologic, geochronologic, and paleomagnetic data for Holocene lava fields along the southern part of the Great Rift, and for late Pleistocene and late middle Pleistocene lava fields in the central and western parts of the quadrangle. These data can contribute to wise management and preservation of the Craters of the Moon National Monument and for broad-scale understanding of the basaltic-volcanic evolution of the eastern Snake River Plain.

Idaho↗

A graphical method for estimation of barometric efficiency from continuous data: Concepts and application to a site in the Piedmont, Air Force Plant 6, Marietta, Georgia

A graphical method that uses continuous water-level and barometric-pressure data was developed to estimate barometric efficiency. A plot of nearly continuous water level (on the y-axis), as a function of nearly continuous barometric pressure (on the x-axis), will plot as a line curved into a series of connected elliptical loops. Each loop represents a barometric-pressure fluctuation. The negative of the slope of the major axis of an elliptical loop will be the ratio of water-level change to barometric-pressure change, which is the sum of the barometric efficiency plus the error. The negative of the slope of the preferred orientation of many elliptical loops is an estimate of the barometric efficiency. The slope of the preferred orientation of many elliptical loops is approximately the median of the slopes of the major axes of the elliptical loops. If water-level change that is not caused by barometric-pressure change does not correlate with barometric-pressure change, the probability that the error will be greater than zero will be the same as the probability that it will be less than zero. As a result, the negative of the median of the slopes for many loops will be close to the barometric efficiency. The graphical method provided a rapid assessment of whether a well was affected by barometric-pressure change and also provided a rapid estimate of barometric efficiency. The graphical method was used to assess which wells at Air Force Plant 6, Marietta, Georgia, had water levels affected by barometric-pressure changes during a 2003 constant-discharge aquifer test. The graphical method was also used to estimate barometric efficiency. Barometric-efficiency estimates from the graphical method were compared to those of four other methods: average of ratios, median of ratios, Clark, and slope. The two methods (the graphical and median-of-ratios methods) that used the median values of water-level change divided by barometric-pressure change appeared to be most resistant to error caused by barometric-pressure-independent water-level change. The graphical method was particularly resistant to large amounts of barometric-pressure-independent water-level change, having an average and standard deviation of error for control wells that was less than one-quarter that of the other four methods. When using the graphical method, it is advisable that more than one person select the slope or that the same person fits the same data several times to minimize the effect of subjectivity. Also, a long study period should be used (at least 60 days) to ensure that loops affected by large amounts of barometric-pressure-independent water-level change do not significantly contribute to error in the barometric-efficiency estimate.

Georgia↗

Simulation of Hydrodynamics and Water Quality in Pueblo Reservoir, Southeastern Colorado, for 1985 through 1987 and 1999 through 2002

Pueblo Reservoir is west of Pueblo, Colorado, and is an important water resource for southeastern Colorado. The reservoir provides irrigation, municipal, and industrial water to various entities throughout the region. In anticipation of increased population growth, the cities of Colorado Springs, Fountain, Security, and Pueblo West have proposed building a pipeline that would be capable of conveying 78 million gallons of raw water per day (240 acre-feet) from Pueblo Reservoir. The U.S. Geological Survey, in cooperation with Colorado Springs Utilities and the Bureau of Reclamation, developed, calibrated, and verified a hydrodynamic and water-quality model of Pueblo Reservoir to describe the hydrologic, chemical, and biological processes in Pueblo Reservoir that can be used to assess environmental effects in the reservoir. Hydrodynamics and water-quality characteristics in Pueblo Reservoir were simulated using a laterally averaged, two-dimensional model that was calibrated using data collected from October 1985 through September 1987. The Pueblo Reservoir model was calibrated based on vertical profiles of water temperature and dissolved-oxygen concentration, and water-quality constituent concentrations collected in the epilimnion and hypolimnion at four sites in the reservoir. The calibrated model was verified with data from October 1999 through September 2002, which included a relatively wet year (water year 2000), an average year (water year 2001), and a dry year (water year 2002). Simulated water temperatures compared well to measured water temperatures in Pueblo Reservoir from October 1985 through September 1987. Spatially, simulated water temperatures compared better to measured water temperatures in the downstream part of the reservoir than in the upstream part of the reservoir. Differences between simulated and measured water temperatures also varied through time. Simulated water temperatures were slightly less than measured water temperatures from March to May 1986 and 1987, and slightly greater than measured data in August and September 1987. Relative to the calibration period, simulated water temperatures during the verification period did not compare as well to measured water temperatures. In general, simulated dissolved-oxygen concentrations for the calibration period compared well to measured concentrations in Pueblo Reservoir. Spatially, simulated concentrations deviated more from the measured values at the downstream part of the reservoir than at other locations in the reservoir. Overall, the absolute mean error ranged from 1.05 (site 1B) to 1.42 milligrams per liter (site 7B), and the root mean square error ranged from 1.12 (site 1B) to 1.67 milligrams per liter (site 7B). Simulated dissolved oxygen in the verification period compared better to the measured concentrations than in the calibration period. The absolute mean error ranged from 0.91 (site 5C) to 1.28 milligrams per liter (site 7B), and the root mean square error ranged from 1.03 (site 5C) to 1.46 milligrams per liter (site 7B). Simulated total dissolved solids generally were less than measured total dissolved-solids concentrations in Pueblo Reservoir from October 1985 through September 1987. The largest differences between simulated and measured total dissolved solids were observed at the most downstream sites in Pueblo Reservoir during the second year of the calibration period. Total dissolved-solids data were not available from reservoir sites during the verification period, so in-reservoir specific-conductance data were compared to simulated total dissolved solids. Simulated total dissolved solids followed the same patterns through time as the measured specific conductance data during the verification period. Simulated total nitrogen concentrations compared relatively well to measured concentrations in the Pueblo Reservoir model. The absolute mean error ranged from 0.21 (site 1B) to 0.27 milligram per liter as nitrogen (sites 3B and 7

Scientific Investigations Report↗

Development and Evaluation of Live-Bed Pier- and Contraction-Scour Envelope Curves in the Coastal Plain and Piedmont Provinces of South Carolina

The U.S. Geological Survey, in cooperation with the South Carolina Department of Transportation, used ground-penetrating radar to collect measurements of live-bed pier scour and contraction scour at 78 bridges in the Piedmont and Coastal Plain Physiographic Provinces of South Carolina. The 151 measurements of live-bed pier-scour depth ranged from 1.7 to 16.9 feet, and the 89 measurements of live-bed contraction-scour depth ranged from 0 to 17.1 feet. Using hydraulic data estimated with a one-dimensional flow model, predicted live-bed scour depths were computed with scour equations from the Hydraulic Engineering Circular 18 and compared with measured scour. This comparison indicated that predicted pier-scour depths generally exceeded the measured pier-scour depths, and at times predicted pier-scour depths were excessive (overpredictions were as large as 23.1 feet). For live-bed contraction-scour depths, predicted scour was sometimes excessive (overpredictions were as large as 14.3 feet), but often observed contraction scour was underpredicted. For live-bed pier scour, trends in laboratory and field data were compared and found to be similar. The strongest explanatory variable was pier width, and an envelope curve for assessing the upper bound of live-bed pier scour was developed using pier width as the primary explanatory variable. Relations in the live-bed contraction-scour data also were investigated, and several envelope curves were developed using the geometric-contraction ratio as the primary explanatory variable. The envelope curves developed with the field data have limitations, but the envelope curves can be used as supplementary tools for assessing the potential for live-bed pier and contraction scour in South Carolina. Data from this study were compiled into a database that includes photographs, measured scour depths, predicted scour depths, limited basin characteristics, limited soil data, and modeled hydraulic data. The South Carolina database can be used in the comparison of sites with similar characteristics to evaluate the potential for scour. In addition, the database can be used to evaluate the performance of various analytical methods for predicting live-bed pier and contraction scour.

South Carolina↗

Hydraulic and water-quality data collection for the investigation of Great Lakes tributaries for Asian carp spawning and egg-transport suitability

If the invasive Asian carps (bighead carp Hypophthalmichthys nobilis and silver carp Hypophthalmichthys molitrix ) migrate to the Great Lakes, in spite of the efforts to stop their advancement, these species will require the fast-flowing water of the Great Lakes tributaries for spawning and recruitment in order to establish a growing population. Two Lake Michigan tributaries (the Milwaukee and St. Joseph Rivers) and two Lake Erie tributaries (the Maumee and Sandusky Rivers) were investigated to determine if these tributaries possess the hydraulic and water-quality characteristics to allow successful spawning of Asian carps. To examine this issue, standard U.S. Geological Survey sampling protocols and instrumentation for discharge and water-quality measurements were used, together with differential global positioning system data for georeferencing. Non-standard data-processing techniques, combined with detailed laboratory analysis of Asian carp egg characteristics, allowed an assessment of the transport capabilities of each of these four tributaries. This assessment is based solely on analysis of observed data and did not utilize the collected data for detailed transport modeling. All four tributaries exhibited potential settling zones for Asian carp eggs both within the estuaries and river mouths and within the lower 100 kilometers (km) of the river. Dams played a leading role in defining these settling zones, with the exception of dams on the Sandusky River. The impoundments created by many of the larger dams on these rivers acted to sufficiently decelerate the flows and allowed the shear velocity to drop below the settling velocity for Asian carp eggs, which would allow the eggs to fall out of suspension and settle on the bottom where it is thought the eggs would perish. While three rivers exhibited these settling zones upstream of the larger dams, not all settling zones are likely to have such effects on egg transport. The Milwaukee River exhibited only a short settling zone upstream of the Grafton Dam, whereas the St. Joseph and Maumee Rivers both had extensive settling zones (>5 km) behind major dams. These longer settling zones are likely to capture more eggs than shorter settling reaches. All four rivers exhibited settling zones at their river mouths, with the Lake Erie tributaries having much larger settling zones extending more than 10 km up the tributaries. While hydraulic data from all four rivers indicated settling of eggs is possible in some locations, all four rivers also exhibited sufficient temperatures, water-quality characteristics, turbulence, and transport times outside of settling zones for successful suspension and development of Asian carp eggs to the hatching stage before the threat of settlement. These observed data indicate that these four Great Lakes tributaries have sufficient hydraulic and water-quality characteristics to support successful spawning and recruitment of Asian carps. The data indicate that with the right temperature and flow conditions, river reaches as short as 25 km may allow Asian carp eggs sufficient time to develop to hatching. Additionally, examining the relation between critical shear velocity and mean velocity, egg settling appears to take place at mean velocities in the range of 15–25 centimeters per second, a much lower value than is generally cited in the literature. A first-order estimate of the minimum transport velocity for Asian carp eggs in a river can be obtained by using mean flow depth and river substrate data, and curves were constructed to show this relation. These findings would expand the number of possible tributaries suitable for Asian carp spawning and contribute to the understanding of how hydraulic and water-quality information can be used to screen additional rivers in the future.

Great Lakes↗

Sediment distribution and hydrologic conditions of the Potomac aquifer in Virginia and parts of Maryland and North Carolina

Sediments of the heavily used Potomac aquifer broadly contrast across major structural features of the Atlantic Coastal Plain Physiographic Province in eastern Virginia and adjacent parts of Maryland and North Carolina. Thicknesses and relative dominance of the highly interbedded fluvial sediments vary regionally. Vertical intervals in boreholes of coarse-grained sediment commonly targeted for completion of water-supply wells are thickest and most widespread across the central and southern parts of the Virginia Coastal Plain. Designated as the Norfolk arch depositional subarea, the entire sediment thickness here functions hydraulically as a single interconnected aquifer. By contrast, coarse-grained sediment intervals are thinner and less widespread across the northern part of the Virginia Coastal Plain and into southern Maryland, designated as the Salisbury embayment depositional subarea. Fine-grained intervals that are generally avoided for completion of water-supply wells are increasingly thick and widespread northward. Fine-grained intervals collectively as thick as several hundred feet comprise two continuous confining units that hydraulically separate three vertically spaced subaquifers. The subaquifers are continuous northward but merge southward into the single undivided Potomac aquifer. Lastly, far southeastern Virginia and northeastern North Carolina are designated as the Albemarle embayment depositional subarea, where both coarse- and fine-grained intervals are of only moderate thickness. The entire sediment thickness functions hydraulically as a single interconnected aquifer. A substantial hydrologic separation from overlying aquifers is imposed by the upper Cenomanian confining unit. Potomac aquifer sediments were deposited by a fluvial depositional complex spanning the Virginia Coastal Plain approximately 100 to 145 million years ago. Westward, persistently uplifted granite and gneiss source rocks sustained a supply of coarse-grained sand and gravel. Immature, high-gradient braided streams deposited longitudinal bars and channel fills across the Norfolk arch subarea. By contrast, across the Salisbury and Albemarle embayment subareas, mature, medium- to low-gradient meandering streams deposited medium- to coarse-grained channel fills and point bars segregated from fine-grained overbank deposits. The Virginia depositional complex merged northward across the Salisbury embayment subarea with another complex in Maryland. Here, additional sediments were received from schist source rocks that underwent three cycles of initial uplift and rapid erosion followed by crustal stability and erosional leveling. Because of the predominance of coarse-grained sediments, transmissivity, hydraulic conductivity, and regional velocities of lateral flow through the Potomac aquifer are greatest across the Norfolk arch depositional subarea, but decrease progressively northward with increasingly fine-grained sediments. Confining units hydraulically separate the Potomac aquifer from overlying aquifers, as indicated by large vertical hydraulic gradients. By contrast, most of the Potomac aquifer internally functions hydraulically as a single interconnected aquifer, as indicated by uniformly small vertical gradients. Most fine-grained sediments within the aquifer do not hydraulically separate overlying and underlying coarse-grained sediments. Across the Salisbury embayment depositional subarea, however, hydraulic separation among the vertically spaced subaquifers is imposed by the intervening confining units. The Potomac aquifer is the largest and most heavily used source of groundwater in the Virginia Coastal Plain. Water-level declines as great as 200 feet create the potential for saltwater intrusion. Conventional stratigraphic correlation has been generally ineffective at accurately characterizing complexly distributed fluvial sediments that compose the Potomac aquifer. Consequently, the aquifer’s internal hydraulic connectivity and overall hydrologic function have not been well understood. Water-supply planning and development efforts have been hampered, and interpretations of regulatory criteria for allowable water-level declines have been ambiguous. An investigation undertaken during 2010–11 by the U.S. Geological Survey, in cooperation with the Virginia Department of Environmental Quality, provides a comprehensive regional description of the spatial distribution of Potomac aquifer sediments and their relation to hydrologic conditions. Altitudes and thicknesses of 2,725 vertical sediment intervals represent the spatial distribution of Potomac aquifer sediments in the Virginia Coastal Plain and adjacent parts of Maryland and North Carolina. Sediment intervals are designated as either dominantly coarse or fine grained and were determined by interpretation of geophysical logs and ancillary information from 456 boreholes. Sediment-interval and borehole summary statistical data indicate regional trends in sediment lithology and stratigraphic continuity, upon which three structurally based and hydrologically distinct sediment depositional subareas are designated. Broad patterns of sediment deposition over time are inferred from published sediment pollen-age data. Discrepancies in previously drawn hydrostratigraphic relations between southeastern Virginia and northeastern North Carolina are partly resolved based on borehole geophysical logs and a recently documented geologic map and corehole. A conceptual model theorizes the depositional history of the sediments and geologically accounts for their distribution. Documented pumping tests of the Potomac aquifer at 197 locations produced 336 values of transmissivity and 127 values of storativity. Based on effective aquifer thicknesses, 296 values of sediment hydraulic conductivity and 113 values of sediment specific storage are calculated. Vertical hydraulic gradients are calculated from 9,479 pairs of water levels measured between November 17, 1953, and October 4, 2011, in 129 closely spaced pairs of wells. Borehole sediment-interval and related data provide a means to achieve high yielding production wells in the Potomac aquifer by site-specific targeting of drilling operations toward water-bearing coarse-grained sand and gravel. Advance knowledge of the potential of different parts of the aquifer also aids in planning optimal groundwater-development areas. Depositional subareas further provide a possible context for resource management. Current (2013) regulatory limits on water-level declines are relative to top surfaces of subdivided upper, middle, and lower Potomac aquifers across the entire Virginia Coastal Plain, but have the potential to exceed the same limit relative to a single undivided Potomac aquifer. By contrast, designation of the sediments as a single aquifer in the Norfolk arch and Albemarle embayment subareas—and as a series of vertically spaced subaquifers and intervening confining units in the Salisbury embayment subarea—best reflects understanding of the Potomac aquifer and can avoid the potential for excessive water-level declines. Simulation modeling to evaluate effects of groundwater withdrawals could be designed similarly, including vertical discretization and (or) zonation of the Potomac aquifer based on depositional subareas and a geostatistical distribution of aquifer properties derived from borehole sediment-interval data. Further resource-management information needs extend beyond the developed part of the Potomac aquifer, particularly across the Northern Neck and Middle Peninsula where only the shallowest part of the aquifer is known, and include structural aspects such as faults, basement bedrock, and the Chesapeake Bay impact crater.

Maryland, North Carolina, Virginia↗

Hydrographic survey of Chaktomuk, the confluence of the Mekong, Tonlé Sap, and Bassac Rivers near Phnom Penh, Cambodia, 2012

The U.S. Geological Survey, in cooperation with the U.S. Department of State, Mekong River Commission, Phnom Penh Autonomous Port, and the Cambodian Ministry of Water Resources and Meteorology, completed a hydrographic survey of Chaktomuk, which is the confluence of the Mekong, Tonlé Sap (also spelled Tônlé Sab), and Bassac Rivers near Phnom Penh, Cambodia. The hydrographic survey used a high-resolution multibeam echosounder mapping system to map the riverbed during April 21–May 2, 2012. The multibeam echosounder mapping system was made up of several components: A RESON Seabat™ 7125 multibeam echosounder, an inertial measurement unit and navigation unit, data collection computers, and a Real-Time Kinematic (RTK) Global Navigation Satellite System (GNSS) base station. The survey area was divided into six survey subreaches and each subreach was surveyed within 3 days along survey lines oriented parallel to the flow direction. Complete coverage of the riverbed was the operational objective; however, to obtain broad spatial coverage, gaps between parallel swaths were permitted, especially in wide, shallow areas where multibeam swath widths were narrow. The survey was referenced to two existing bench marks with known geographic coordinates by establishing a GNSS base station on the bench marks each day and using real-time corrections from the base station to correct boat navigation data. The World Geodetic System of 1984 (WGS 84) ellipsoid was used during data collection to reference height, and data were adjusted to the local datum, Ha Tien 1960, during postprocessing. The quality of hydrographic surveys was described by an uncertainty estimate called total propagated uncertainty (TPU). Calculations of TPU were completed for the hydrographic survey data resulting in the maximum TPU of 0.33 meters. The mean and median TPUs were 0.18 meters, and 99.9 percent of TPU values were less than 0.25 meters. Detailed hydrographic maps of Mekong, Tonlé Sap, and Bassac Rivers showing the riverbed elevations surveyed April 21–May 2, 2012, referenced to Ha Tien 1960 were produced. The surveyed area included a 2-km stretch of the Mekong River between the confluence with the Tonlé Sap and Bassac Rivers, and extended 4 km upstream and 3.6 km downstream from the 2,000-m confluence stretch of the Mekong River. In addition, 0.7 km of the Bassac River downstream and 3.5 km of the Tonlé Sap River (from the confluence to Chroy Changvar Bridge) upstream from their confluence with the Mekong River were surveyed. Riverbed features (such as dunes, shoals, and the effects of sediment mining, which were observed during data collection) are visible on the hydrographic maps. All surveys were completed at low annual water levels as referenced to nearby Mekong River Commission streamflow-gaging stations. Riverbed elevations surveyed ranged from 24.08 m below to 1.54 m above Ha Tien 1960.

Phnom Penh↗