Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Nature, Scientific Data”

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 721 records · Page 40Linked to original sources

Completion summary for borehole TAN-2312 at Test Area North, Idaho National Laboratory, Idaho

In 2017, the U.S. Geological Survey, in cooperation with the U.S. Department of Energy, drilled and constructed borehole TAN-2312 for stratigraphic framework analyses and long-term groundwater monitoring of the eastern Snake River Plain aquifer at the Idaho National Laboratory in southeast Idaho. The location of borehole TAN-2312 was selected because it was downgradient from TAN and believed to be the outer extent of waste plumes originating from the TAN facility. Borehole TAN-2312 initially was cored to collect continuous geologic data, and then re-drilled to complete construction as a monitor well. The final construction for borehole TAN-2312 required 16- and 10-inch (in.) diameter carbon-steel well casing to 37 and 228 feet below land surface (ft BLS), respectively, and 9.9-in. diameter open-hole completion below the casing to 522 ft BLS. Depth to water is measured near 244 ft BLS. Following construction and data collection, a temporary submersible pump and water-level access line were placed near 340 ft BLS to allow for aquifer testing, for collecting periodic water samples, and for measuring water levels. Borehole TAN-2312 was cored continuously, starting at the first basalt contact (about 37 ft BLS) to a depth of 568 ft BLS. Not including surface sediment (0–37 ft), recovery of basalt and sediment core at borehole TAN-2312 was about 93 percent; however, core recovery from 170 to 568 ft BLS was 100 percent. Based on visual inspection of core and geophysical data, basalt examined from 37 to 568 ft BLS consists of about 32 basalt flows that range from approximately 3 to 87 ft in thickness and 4 sediment layers with a combined thickness of approximately 76 ft. About 2 ft of total sediment was described for the saturated zone, observed from 244 to 568 ft BLS, near 296 and 481 ft BLS. Sediment described for the saturated zone were composed of fine-grained sand and silt with a lesser amount of clay. Basalt texture for borehole TAN-2312 generally was described as aphanitic, phaneritic, and porphyritic. Basalt flows varied from highly fractured to dense with high to low vesiculation. Geophysical and borehole video logs were collected after core drilling and after final construction at borehole TAN-2312. Geophysical logs were examined synergistically with available core material to suggest zones where groundwater flow was anticipated. Natural gamma log measurements were used to assess sediment layer thickness and location. Neutron and gamma-gamma source logs were used to identify fractured areas for aquifer testing. Acoustic televiewer logs, fluid logs, and electromagnetic flow meter results were used to identify fractures and assess groundwater movement when compared against neutron measurements. Furthermore, gyroscopic deviation measurements were used to measure horizontal and vertical displacement for borehole TAN-2312. After construction of borehole TAN-2312, a single-well aquifer test was completed September 27, 2017, to provide estimates of transmissivity and hydraulic conductivity. Estimates for transmissivity and hydraulic conductivity were 1.51×102 feet squared per day and 0.23 feet per day, respectively. During the 220-minute aquifer test, well TAN-2312 had about 23 ft of measured drawdown at sustained pumping rate of 27.2 gallons per minute. The transmissivity and hydraulic conductivity estimates for well TAN-2312 were lower than the values determined from previous aquifer tests in other wells near Test Area North. Water samples were analyzed for cations, anions, metals, nutrients, volatile organic compounds, stable isotopes, and radionuclides. Water samples for most of the inorganic constituents showed concentrations near background levels for eastern regional groundwater. Water samples for stable isotopes of oxygen, hydrogen, and sulfur indicated some possible influence of irrigation on the water quality. The volatile organic compound data indicated that this well had some minor influence by wastewater disposal practices at Test Area North.

Idaho↗

Estimated reductions in phosphorus loads from removal of leaf litter in the Lake Champlain drainage area, Vermont

Excess nutrient loading and other factors are driving eutrophication and other negative effects on water-quality conditions in Lake Champlain and other receiving waters in Vermont. Two common best management practices were evaluated to determine how these practices can be optimized by targeting maintenance and operation to align better with seasonally driven needs, specifically to help municipalities remove a greater proportion of seasonal leaves and organic debris, reduce nutrient loading, and achieve water-quality goals. To characterize solid materials typically removed by the municipal BMPs of catch-basin (CB) cleaning and street cleaning (SC), subsamples of CB and SC materials were collected each month from nine participating municipalities in central and northwestern Vermont between September 2017 and November 2018. Monthly and seasonal composites of CB and SC samples were created from the subsamples of available materials from all municipalities. Samples were analyzed for concentrations of total organic carbon, total Kjeldahl nitrogen, and total phosphorus (P), and separated into three particle-size fractions. Distribution of particle-size fractions was similar between CB and SC as both practices generally collect the coarser fraction of solid materials (greater than 125 micrometers in diameter). In the fall, however, the range of the coarser fraction of materials increased. This is attributed to the ability of SC to collect leaves and other light organic materials that commonly pass through a CB system designed to trap heavier materials. Total organic carbon, total Kjeldahl nitrogen, and total P concentrations were highest in the catch-basin samples in the fall of 2017, and concentrations in the SC samples were highest in the fall of 2018. The collection of fewer samples in 2017 may account for some of the variability between fall 2017 and fall 2018 results. A subset of SC samples collected from piles representing specific street-cleaning routes in September and November 2018 were also analyzed. Materials collected in November were dominated by leaves, and the concentrations of the analyzed species of carbon, nitrogen, and phosphorus in some samples were more than double those in samples collected on the same street-cleaning routes in September. The Vermont Department of Environmental Conservation and the University of Vermont developed estimates of load-reduction credits for CB and SC practices based on a policy developed by the Wisconsin Department of Natural Resources that determined the potential for credits associated with leaf-removal activities. This process also considered BMPs that were initiated during the U.S. Environmental Protection Agency’s Lake Champlain Basin Total Maximum Daily Load monitoring period (2000 to 2009) and adapted the Wisconsin Department of Natural Resources policies to apply to existing SC routes in the cooperating Vermont municipalities that possessed at least 17 percent tree cover. This exercise demonstrated that applying the Wisconsin Department of Natural Resources policy to existing street-cleaning routes possessing 17 percent or more tree cover would result in reductions in total P loads up to 65 percent of mandated target reductions, and about a 25 percent reduction on average. Continuous simulations of stormwater runoff volume, and of loads of suspended sediments and total P, also were created for Englesby Brook Basin, an urbanized basin in Burlington and South Burlington that drains to Lake Champlain. Although the basin is more developed than the average of the nine cooperating municipalities, streamflow and P loading data collected by the U.S. Geological Survey were available to evaluate model performance. Simulations based on a year of average climatic conditions projected potential small reductions in total P of 0.08 to 0.10 percent as a result of CB cleaning and SC practices. Simulated weekly SC practices, however, reduced street-solid loads by as much as 7 percent. When the proportion of total P seen in fall SC materials collected in Vermont was applied to these simulated street-solid loads, estimated reductions of total P were about 29 percent. The combination of analytical results, estimated load-reduction credits, and simulated reductions indicate that targeted increases of SC activities to reduce leaf loading in the fall have the potential to reduce loading to receiving waters and could help regulated communities meet their water-quality goals.

Vermont↗

Watershed Boundary Dataset for Mississippi

The U.S. Geological Survey, in cooperation with the Mississippi Department of Environmental Quality, U.S. Department of Agriculture-Natural Resources Conservation Service, Mississippi Department of Transportation, U.S. Department of Agriculture-Forest Service, and the Mississippi Automated Resource Information System developed a 1:24,000-scale Watershed Boundary Dataset for Mississippi including watershed and subwatershed boundaries, codes, names, and areas. The Watershed Boundary Dataset for Mississippi provides a standard geographical framework for water-resources and selected land-resources planning. The original 8-digit subbasins (Hydrologic Unit Codes) were further subdivided into 10-digit watersheds (62.5 to 391 square miles (mi2)) and 12-digit subwatersheds (15.6 to 62.5 mi2) - the exceptions being the Delta part of Mississippi and the Mississippi River inside levees, which were subdivided into 10-digit watersheds only. Also, large water bodies in the Mississippi Sound along the coast were not delineated as small as a typical 12-digit subwatershed. All of the data - including watershed and subwatershed boundaries, subdivision codes and names, and drainage-area data - are stored in a Geographic Information System database, which are available at: http://ms.water.usgs.gov/. This map shows information on drainage and hydrography in the form of U.S. Geological Survey hydrologic unit boundaries for water-resource 2-digit regions, 4-digit subregions, 6-digit basins (formerly called accounting units), 8-digit subbasins (formerly called cataloging units), 10-digit watershed, and 12-digit subwatersheds in Mississippi. A description of the project study area, methods used in the development of watershed and subwatershed boundaries for Mississippi, and results are presented in Wilson and others (2008). The data presented in this map and by Wilson and others (2008) supersede the data presented for Mississippi by Seaber and others (1987) and U.S. Geological Survey (1977).

Scientific Investigations Map↗

Anatomical and genetic variation of western Oxyloma (Pulmonata: Succineidae) concerning the endangered Kanab ambersnail ( Oxyloma haydeni kanabense ) in Arizona and Utah

The land snail genus Oxyloma (Pulmonata: Succineidae) includes the Federally endangered Kanab ambersnail (Oxyloma haydeni kanabense Pilsbry), which is known at the time of this study from only two locations in the United States: Three Lakes, Utah, and Vaseys Paradise, Arizona, on the Colorado River in Grand Canyon National Park. Since 1994, the Kanab ambersnail has received much attention because its presence at Vaseys Paradise has implications for the ecosystem-wide management of the Colorado River. This attention is primarily because an experimental high-flow release of water from Glen Canyon Dam in 1996 destroyed or degraded Kanab ambersnail habitat at Vaseys Paradise. This experimental high flow was designed to replicate natural flow regimes throughout the Grand Canyon river corridor. However, as a result of the habitat destruction at Vaseys Paradise, in 1996, the U.S. Fish and Wildlife Service ruled that no further experimental high-discharge floods could be carried out until additional Kanab ambersnail populations were discovered or established. This mandate created a situation where the management of a single endangered species conflicted directly with the management of an entire ecosystem. Although since 1996, the U.S. Fish and Wildlife Service has permitted the use of flows as high as stage heights equivalent to 44,000 cubic feet per second, higher flows were requested by various Grand Canyon stakeholders and scientists but were not possible owing to low storage of Lake Powell. Adding to the controversy about Oxyloma and the Kanab ambersnail were previous anatomical and genetic analyses of the genus, which showed that genetic characteristics of specimens did not correspond with their identifications based on traditional taxonomic criteria, raising questions about the validity of the taxonomy of Oxyloma and the protected status of Kanab ambersnails. Specifically, a previous study suggested that the endangered Kanab ambersnail population at Three Lakes was more closely related to other, non-endangered ambersnail populations across the Southwest. In contrast, the Kanab ambersnail population at Vaseys Paradise appeared to be genetically distinct from all other ambersnail populations studied. Management options for the ambersnail population at Vaseys Paradise, at the time of this study, conflict with ecosystem-wide measures proposed to benefit other natural resources in the Grand Canyon. The U.S. Fish and Wildlife Service will not revise the 1995 Kanab Ambersnail Recovery Plan until further genetic and anatomical analyses provide more fine-scale taxonomic resolution of the identity of Oxyloma populations on the Colorado Plateau and elsewhere in the American Southwest. Likewise, interagency cooperators cannot revise down-listing criteria for the Kanab ambersnail until substantial evidence is provided identifying distinct Oxyloma taxa or a larger group of conspecifics that reasonably could be managed as one species. Therefore, given the current controversy about the taxonomy of Oxyloma and the endangered Kanab ambersnail, new detailed analyses were completed of morphological and genetic variation from many Oxyloma specimens collected at 12 western North American locations. These new data have allowed us to evaluate many issues related to Kanab ambersnail taxonomy. Using this dataset, the study of shells and anatomy indicates that the holotype of Oxyloma haydeni kanabense plausibly can be regarded as a member of the same species as the populations of Oxyloma analyzed in this study. Additionally, the presence of gene flow among all populations is evidence that they are members of the same species. Almost all the observed genetic diversity can be accounted for by short-distance or long-distance dispersal events between populations in this study. Our major taxonomic conclusion is that all samples collected for this study were drawn from populations of the same species.

Arizona;Utah↗

Landscape and climatic influences on actual evapotranspiration and available water using the Operational Simplified Surface Energy Balance (SSEBop) Model in eastern Bernalillo County, New Mexico, 2015

The U.S. Geological Survey, in cooperation with the Bernalillo County Public Works Division, conducted a 1-year study in 2015 to assess the spatial and temporal distribution of evapotranspiration (ET) and available water within the East Mountain area in Bernalillo County, New Mexico. ET and available water vary spatiotemporally because of complex interactions among environmental factors, including vegetation characteristics, soil characteristics, topography, and climate. Precipitation data from the Parameter-Elevation Regressions on Independent Slopes Model (PRISM) ( P ) were used in conjunction with actual ET ( ETa ) data from the Operational Simplified Surface Energy Balance (SSEBop) model to estimate available water ( P – ETa ) at 100-meter (m) resolution in the study area. Maps, descriptive statistics, boxplots, regression analyses (continuous data), and multiple comparison tests (categorical data) were used to characterize P , ETa , and available water and their relations to topographic, soil, and vegetation datasets in the East Mountain area. Five categories of the natural land-cover type (evergreen forest, shrub, herbaceous, deciduous forest, and mixed forest) and four categories of developed land-cover type specific to residential intensity (developed open, developed low, developed medium, and developed high) were analyzed individually and in interaction with multiple elevation, tree canopy, and soil texture classes. Annual mean P in 2015 in the East Mountain area was 608 millimeters (mm), and annual mean ETa was 543 mm (89 percent of annual P in 2015), indicating that in 2015, a spatial mean of about 65 mm of water was available for runoff, soil moisture replenishment, or groundwater recharge. Monthly ETa was greatest in July and smallest in January. The intervening months did not show smooth temporal or consistent spatial changes from month to month. Months with lower ETa (January to March, October to December) also tended to have greater available water, indicating that soil moisture (water supply) and potential ET (water demand) may have been out of phase. Regression analyses showed that monthly ETa data had the highest correlation with annual ETa among the atmospheric, topographic, soil, or vegetation datasets, particularly during the early and late growing season (March, April, May, and September). In contrast, monthly P was highly variable and not as highly correlated with annual ETa . Among landscape variables, correlations with annual ETa were highest for tree canopy cover (coefficient of determination [R 2 ] = 0.46). Correlations between ETa and other landscape variables were lower (R 2 = 0.06–0.19): available soil water in the top 100 centimeters, soil bulk density of layer 1, slope, sand content of soil layer 1, soil depth, available soil water in the top 25 centimeters, leaf area index, aspect eastness, and elevation. Evergreen forest areas had the highest annual median ETa , followed by mixed forest, open residential areas, and deciduous forest. Available water typically was higher in landcover types with lower ETa : herbaceous cover, followed by deciduous forest, high-intensity developed areas, and shrub. Deciduous forest had the second highest median available water, despite having the fourth highest ETa , because deciduous forest had greater P than most other areas. Annual median ETa typically was greatest in the second highest elevation band (2,401–2,800 m above the North American Vertical Datum of 1988 [NAVD 88]), and lower in the highest elevation band (2,801–3,254 m above NAVD 88), despite having greater P , likely because of decreased tree canopy cover or a shift from evergreen to deciduous trees at the highest elevations. Annual median ETa increased with tree canopy cover, regardless of landcover type. ETa correlation was higher with tree canopy than with leaf area index or normalized difference vegetation index. This result indicates that it is important to include the thermal band (from satellite multispectral data) in vegetation indices used to describe ETa , perhaps to account for the influence of energy limitation or water limitation on ET. Of all natural landcover types, finer soils had the most available water, whereas coarser soils had the least available water. Relations of soil type with P – ETa were different than with ETa , indicating ET and available water have a complex response to differences in soil type. Further modeling would be useful in determining soils’ infiltration, storage, conductivity, and plant-water availability relations to individual storms for each position in the landscape, as well as the corresponding effects of these processes on ET and available water. The best multivariate linear model for annual ETa had an R 2 value of 0.62. Monthly ETa models had R 2 values between 0.16 and 0.65. Models usually, but not always, performed best during the growing season. These results indicate that even the best multivariate linear models cannot explain a notable amount of the variability in ET. The monthly ETa models with the highest correlations (August and September) followed a July having almost twice the mean precipitation for July (1981–2010), which indicates that a soil-moisture variable is needed to more accurately model monthly ETa . Further study is needed to better characterize this system, the variables that affect ET and available water, and the partitioning of available water into runoff, soil moisture storage, and groundwater recharge.

New Mexico↗

Quaternary sediment thickness and bedrock topography of the glaciated United States east of the Rocky Mountains

Beginning roughly 2.6 million years ago, global climate entered a cooling phase known as the Pleistocene Epoch. As snow in northern latitudes compacted into ice several kilometers thick, it flowed as glaciers southward across the North American continent. These glaciers extended across the northern United States, dramatically altering the landscape they covered. East of the Rocky Mountains, the ice coalesced into continental glaciers (called the Laurentide Ice Sheet) that at times blanketed much of the north-central and northeastern United States. To the west of the Laurentide Ice Sheet, glaciers formed in the mountains of western Canada and the United States and coalesced into the Cordilleran ice sheet; this relatively smaller ice mass extended into the conterminous United States in the northernmost areas of western Montana, Idaho, and Washington. Throughout the Pleistocene, landscape alteration occurred by (1) glacial erosion of the rocks and sediments; (2) redeposition of the eroded earth materials in a form substantially different from their source rocks, in terms of texture and overall character; and (3) disruption of preexisting drainage patterns by the newly deposited sediments. In many cases, pre-glacial drainage systems (including, for example, the Mississippi River) were rerouted because their older drainage courses became blocked with glacial sediment. The continental glaciers advanced and retreated many times across those areas. During each ice advance, or glaciation, erosion and deposition occurred, and the landscape was again altered. Through successive glaciations, the landscape and the bedrock surface gradually came to resemble their present configurations. As continental ice sheets receded and the Pleistocene ended, erosion and deposition of sediment (for example in stream valleys) continued to shape the landscape up to the present day (albeit to a lesser extent than during glaciation). The interval of time since the last recession of the glaciers is called the Holocene and, together with the Pleistocene, constitutes the Quaternary Period of geologic time; this publication characterizes the three-dimensional geometry of the Quaternary sediments and the bedrock surface that lies beneath. The pre-glacial landscape was underlain mostly by weathered bedrock generally similar in nature to that found in many areas of the non-glaciated United States. Glacial erosion and redeposition of earth materials produced a young, mineral-rich soil that formed the basis for the highly productive agricultural economy in the U.S. midcontinent. Extensive buried sands and gravels within the glacial deposits also provided a stimulus to other economic sectors by serving as high-quality aquifers supplying groundwater to the region’s industry and cities. An understanding of the three-dimensional distribution of these glacial sediments has direct utility for addressing various societal issues including groundwater quality and supply, and landscape and soil response to earthquake-induced shaking. The Quaternary sediment thickness map and bedrock topographic map shown here provide a regional overview and are intended to supplement the more detailed work on which they are based. Detailed mapping is particularly useful in populated areas for site-specific planning. In contrast, regional maps such as these serve to place local, detailed mapping in context; to permit the extrapolation of data into unmapped areas; and to depict large-scale regional geologic features and patterns that are beyond the scope of local, detailed mapping. They also can enhance the reader’s general understanding of the region’s landscape and geologic history and provide a source of information for regional decision making that could benefit by improved predictability of bedrock depth beneath the unconsolidated Quaternary sediments. To enable these maps to be analyzed in conjunction with other types of information, this publication also includes the map data in GIS compatible format.

Scientific Investigations Map↗

Implementation of local grid refinement (LGR) for the Lake Michigan Basin regional groundwater-flow model

The U.S. Geological Survey is evaluating water availability and use within the Great Lakes Basin. This is a pilot effort to develop new techniques and methods to aid in the assessment of water availability. As part of the pilot program, a regional groundwater-flow model for the Lake Michigan Basin was developed using SEAWAT-2000. The regional model was used as a framework for assessing local-scale water availability through grid-refinement techniques. Two grid-refinement techniques, telescopic mesh refinement and local grid refinement, were used to illustrate the capability of the regional model to evaluate local-scale problems. An intermediate model was developed in central Michigan spanning an area of 454 square miles (mi2) using telescopic mesh refinement. Within the intermediate model, a smaller local model covering an area of 21.7 mi2 was developed and simulated using local grid refinement. Recharge was distributed in space and time using a daily output from a modified Thornthwaite-Mather soil-water-balance method. The soil-water-balance method derived recharge estimates from temperature and precipitation data output from an atmosphere-ocean coupled general-circulation model. The particular atmosphere-ocean coupled general-circulation model used, simulated climate change caused by high global greenhouse-gas emissions to the atmosphere. The surface-water network simulated in the regional model was refined and simulated using a streamflow-routing package for MODFLOW. The refined models were used to demonstrate streamflow depletion and potential climate change using five scenarios. The streamflow-depletion scenarios include (1) natural conditions (no pumping), (2) a pumping well near a stream; the well is screened in surficial glacial deposits, (3) a pumping well near a stream; the well is screened in deeper glacial deposits, and (4) a pumping well near a stream; the well is open to a deep bedrock aquifer. Results indicated that a range of 59 to 50 percent of the water pumped originated from the stream for the shallow glacial and deep bedrock pumping scenarios, respectively. The difference in streamflow reduction between the shallow and deep pumping scenarios was compensated for in the deep well by deriving more water from regional sources. The climate-change scenario only simulated natural conditions from 1991-2044, so there was no pumping stress simulated. Streamflows were calculated for the simulated period and indicated that recharge over the period generally increased from the start of the simulation until approximately 2017, and decreased from then to the end of the simulation. Streamflow was highly correlated with recharge so that the lowest streamflows occurred in the later stress periods of the model when recharge was lowest.

Scientific Investigations Report↗

Quality of stormwater runoff discharged from Massachusetts highways, 2005-07

The U.S. Geological Survey (USGS), in cooperation with U.S. Department of Transportation Federal Highway Administration and the Massachusetts Department of Transportation, conducted a field study from September 2005 through September 2007 to characterize the quality of highway runoff for a wide range of constituents. The highways studied had annual average daily traffic (AADT) volumes from about 3,000 to more than 190,000 vehicles per day. Highway-monitoring stations were installed at 12 locations in Massachusetts on 8 highways. The 12 monitoring stations were subdivided into 4 primary, 4 secondary, and 4 test stations. Each site contained a 100-percent impervious drainage area that included two or more catch basins sharing a common outflow pipe. Paired primary and secondary stations were located within a few miles of each other on a limited-access section of the same highway. Most of the data were collected at the primary and secondary stations, which were located on four principal highways (Route 119, Route 2, Interstate 495, and Interstate 95). The secondary stations were operated simultaneously with the primary stations for at least a year. Data from the four test stations (Route 8, Interstate 195, Interstate 190, and Interstate 93) were used to determine the transferability of the data collected from the principal highways to other highways characterized by different construction techniques, land use, and geography. Automatic-monitoring techniques were used to collect composite samples of highway runoff and make continuous measurements of several physical characteristics. Flowweighted samples of highway runoff were collected automatically during approximately 140 rain and mixed rain, sleet, and snowstorms. These samples were analyzed for physical characteristics and concentrations of 6 dissolved major ions, total nutrients, 8 total-recoverable metals, suspended sediment, and 85 semivolatile organic compounds (SVOCs), which include priority polyaromatic hydrocarbons (PAHs), phthalate esters, and other anthropogenic or naturally occurring organic compounds. The distribution of particle size of suspended sediment also was determined for composite samples of highway runoff. Samples of highway runoff were collected year round and under various dry antecedent conditions throughout the 2-year sampling period. In addition to samples of highway runoff, supplemental samples also were collected of sediment in highway runoff, background soils, berm materials, maintenance sands, deicing compounds, and vegetation matter. These additional samples were collected near or on the highways to support data analysis. There were few statistically significant differences between populations of constituent concentrations in samples from the primary and secondary stations on the same principal highways (Mann-Whitney test, 95-percent confidence level). Similarly, there were few statistically significant differences between populations of constituent concentrations for the four principal highways (data from the paired primary and secondary stations for each principal highway) and populations for test stations with similar AADT volumes. Exceptions to this include several total-recoverable metals for stations on Route 2 and Interstate 195 (highways with moderate AADT volumes), and for stations on Interstate 95 and Interstate 93 (highways with high AADT volumes). Supplemental data collected during this study indicate that many of these differences may be explained by the quantity, as well as the quality, of the sediment in samples of highway runoff. Nonparametric statistical methods also were used to test for differences between populations of sample constituent concentrations among the four principal highways that differed mainly in traffic volume. These results indicate that there were few statistically significant differences (Mann-Whitney test, 95-percent confidence level) for populations of concentrations of most total-recoverable metals

Scientific Investigations Report↗

Thermal mapping of Hawaiian volcanoes with ASTER satellite data

Thermal mapping of volcanoes is important to determine baseline thermal behavior in order to judge future thermal activity that may precede an eruption. We used cloud-free kinetic temperature images from the ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer) sensor obtained between 2000 and 2010 to produce thermal maps for all five subaerial volcanoes in Hawai‘i that have had eruptions in the Holocene (Kīlauea, Mauna Loa, Hualālai, Mauna Kea, and Haleakalā). We stacked the images to provide time-averaged thermal maps, as well as to analyze temperature trends through time. Thermal areas are conspicuous at the summits and rift zones of Kīlauea and Mauna Loa, and the summit calderas of these volcanoes contain obvious arcuate, concentric linear thermal areas that probably result from channeling of rising gas along buried, historical intracaldera scarps. The only significant change in thermal activity noted in the study period is the opening of the Halema‘uma‘u vent at Kīlauea's summit in 2008. Several small thermal anomalies are coincident with pit craters on Hualālai. We suspect that these simply result from the sheltered nature of the depression, but closer inspection is warranted to determine if genuine thermal activity exists in the craters. Thermal areas were not detected on Haleakalā or Mauna Kea. The main limitation of the study is the large pixel size (90 m) of the ASTER images, which reduces our ability to detect subtle changes or to identify small, low-temperature thermal activity. This study, therefore, is meant to characterize the broad, large-scale thermal features on these volcanoes. Future work should study these thermal areas with thermal cameras and thermocouples, which have a greater ability to detect small, low-temperature thermal features.

Hawai'i↗

Potential sources of analytical bias and error in selected trace element data-quality analyses

Potential sources of analytical bias and error associated with laboratory analyses for selected trace elements where concentrations were greater in filtered samples than in paired unfiltered samples were evaluated by U.S. Geological Survey (USGS) Water Quality Specialists in collaboration with the USGS National Water Quality Laboratory (NWQL) and the Branch of Quality Systems (BQS). Causes for trace-element concentrations in filtered samples to exceed those in associated unfiltered samples have been attributed to variability in analytical measurements, analytical bias, sample contamination either in the field or laboratory, and (or) sample-matrix chemistry. These issues have not only been attributed to data generated by the USGS NWQL but have been observed in data generated by other laboratories. This study continues the evaluation of potential analytical bias and error resulting from matrix chemistry and instrument variability by evaluating the performance of seven selected trace elements in paired filtered and unfiltered surface-water and groundwater samples collected from 23 sampling sites of varying chemistries from six States, matrix spike recoveries, and standard reference materials. Filtered and unfiltered samples have been routinely analyzed on separate inductively coupled plasma-mass spectrometry instruments. Unfiltered samples are treated with hydrochloric acid (HCl) during an in-bottle digestion procedure; filtered samples are not routinely treated with HCl as part of the laboratory analytical procedure. To evaluate the influence of HCl on different sample matrices, an aliquot of the filtered samples was treated with HCl. The addition of HCl did little to differentiate the analytical results between filtered samples treated with HCl from those samples left untreated; however, there was a small, but noticeable, decrease in the number of instances where a particular trace-element concentration was greater in a filtered sample than in the associated unfiltered sample for all trace elements except selenium. Accounting for the small dilution effect (2 percent) from the addition of HCl, as required for the in-bottle digestion procedure for unfiltered samples, may be one step toward decreasing the number of instances where trace-element concentrations are greater in filtered samples than in paired unfiltered samples. The laboratory analyses of arsenic, cadmium, lead, and zinc did not appear to be influenced by instrument biases. These trace elements showed similar results on both instruments used to analyze filtered and unfiltered samples. The results for aluminum and molybdenum tended to be higher on the instrument designated to analyze unfiltered samples; the results for selenium tended to be lower. The matrices used to prepare calibration standards were different for the two instruments. The instrument designated for the analysis of unfiltered samples was calibrated using standards prepared in a nitric:hydrochloric acid (HNO 3 :HCl) matrix. The instrument designated for the analysis of filtered samples was calibrated using standards prepared in a matrix acidified only with HNO 3 . Matrix chemistry may have influenced the responses of aluminum, molybdenum, and selenium on the two instruments. The best analytical practice is to calibrate instruments using calibration standards prepared in matrices that reasonably match those of the samples being analyzed. Filtered and unfiltered samples were spiked over a range of trace-element concentrations from less than 1 to 58 times ambient concentrations. The greater the magnitude of the trace-element spike concentration relative to the ambient concentration, the greater the likelihood spike recoveries will be within data control guidelines (80–120 percent). Greater variability in spike recoveries occurred when trace elements were spiked at concentrations less than 10 times the ambient concentration. Spike recoveries that were considerably lower than 90 percent often were associated with spiked concentrations substantially lower than what was present in the ambient sample. Because the main purpose of spiking natural water samples with known quantities of a particular analyte is to assess possible matrix effects on analytical results, the results of this study stress the importance of spiking samples at concentrations that are reasonably close to what is expected but sufficiently high to exceed analytical variability. Generally, differences in spike recovery results between paired filtered and unfiltered samples were minimal when samples were analyzed on the same instrument. Analytical results for trace-element concentrations in ambient filtered and unfiltered samples greater than 10 and 40 μg/L, respectively, were within the data-quality objective for precision of ±25 percent. Ambient trace-element concentrations in filtered samples greater than the long-term method detection limits but less than 10 μg/L failed to meet the data-quality objective for precision for at least one trace element in about 54 percent of the samples. Similarly, trace-element concentrations in unfiltered samples greater than the long-term method detection limits but less than 40 μg/L failed to meet this data-quality objective for at least one trace-element analysis in about 58 percent of the samples. Although, aluminum and zinc were particularly problematic, limited re-analyses of filtered and unfiltered samples appeared to improve otherwise failed analytical precision. The evaluation of analytical bias using standard reference materials indicate a slight low bias for results for arsenic, cadmium, selenium, and zinc. Aluminum and molybdenum show signs of high bias. There was no observed bias, as determined using the standard reference materials, during the analysis of lead.

Scientific Investigations Report↗

Maintaining momentum for collaborative working groups in a post-pandemic world

Scientific progress depends in part on our ability to synthesize heterogeneous data and ideas into new models and paradigms. In environmental sciences, such synthesis has been particularly effective when conducted by collaborative working groups: diverse groups of researchers and practitioners brought together for a concentrated period of collaboration on key questions. Such work is often done at synthesis centres: organizations that promote, fund, organize and host working groups and other collaborative research and training activities 1 . However, because of the COVID-19 pandemic, synthesis centres have had to rapidly adapt to supporting fully virtual working groups; the eight centres we direct supported 68 virtual working groups in the past year. Based on this experience, we conclude — contrary to a recent editorial on conferences published in this journal 2 — that virtual gatherings, while providing a bridge during the pandemic, cannot replace immersive, in-person collaborations. In-person working group meetings involve productive and varied interactions for many hours over consecutive days. We have found that virtual sessions lose effectiveness after a few hours, as participants become fatigued from staring at a screen or juggling local demands. While virtual meetings can work for short, well-delineated tasks, they are less suited for unstructured and free-flowing discussions — and thus struggle to create the social cohesion and trust known to fuel creative breakthroughs during week-long in-person meetings 3 .

Nature Ecology and Evolution↗

Hydrologic assessment of three drainage basins in the Pinelands of southern New Jersey, 2004-06

The New Jersey Pinelands is an ecologically diverse area in the southern New Jersey Coastal Plain, most of which overlies the Kirkwood-Cohansey aquifer system. The demand for groundwater from this aquifer system is increasing as local development increases. Because any increase in groundwater withdrawals has the potential to affect streamflows and wetland water levels, and ultimately threaten the ecological health and diversity of the Pinelands ecosystem, the U.S. Geological Survey, in cooperation with the New Jersey Pinelands Commission, began a multi-phase hydrologic investigation in 2004 to characterize the hydrologic system supporting the aquatic and wetland communities of the New Jersey Pinelands area (Pinelands). The current investigation of the hydrology of three representative drainage basins in the Pinelands (Albertson Brook, McDonalds Branch, and Morses Mill Stream basins) included a compilation of existing data; collection of water-level and streamflow data; mapping of the water-table altitude and depth to the water table; and analyses of water-level and streamflow variability, subsurface gradients and flow patterns, and water budgets. During 2004-06, a hydrologic database of existing and new data from wells and stream sites was compiled. Methods of data collection and analysis were defined, and data networks consisting of 471 wells and 106 surface-water sites were established. Hydrographs from 26 water-level-monitoring wells and four streamflow-gaging stations were analyzed to show the response of water levels and streamflow to precipitation and recharge with respect to the locations of these wells and streams within each basin. Water-level hydrographs show varying hydraulic gradients and flow potentials, and indicate that responses to recharge events vary with well depth and proximity to recharge and discharge areas. Results of the investigation provide a detailed characterization of hydrologic conditions, processes, and relations among the components of the hydrologic cycle in the Pinelands. In the Pinelands, recharge replenishes the aquifer system and contributes to groundwater flow, most of which moves to wetlands and surface water where natural discharge occurs. Some groundwater flow is intercepted by supply wells. Recharge rates generally are highest during the non-growing season and are inversely related to evapotranspiration. Analysis of subsurface hydraulic gradients, water-table fluctuations, and streamflow variability indicates a strong linkage between groundwater and wetlands, lakes and streams. Gradient analysis indicates that most wetlands are in groundwater discharge areas, but some wetlands are in groundwater recharge areas. The depth to the water table ranges from zero at surface-water features up to about 10 meters in topographically high areas. Depth to water fluctuates seasonally, and the magnitude of these fluctuations generally increases with distance from surface water. Variations in the permeability of the soils and sediments of the aquifer system strongly affect patterns of water movement through the subsurface and the interaction of groundwater with wetlands, lakes and streams. Mean annual streamflow during 2004-06 ranged from 83 to 106 percent of the long-term mean annual discharge, indicating that the data-collection period can be considered representative of average conditions. Measurements of groundwater levels, stream stage, and stream discharge and locations of start-of-flow are illustrated in basin-wide maps of water-table altitude, depth to the water table, and stream base flow during the period. Water-level data collected along 15 hydrologic transects that span the range of environments from uplands through wetlands to surface water were used to determine hydraulic gradients, potential flow directions, and areas of recharge and discharge. These data provide information about the localized interactions of groundwater with wetlands and surface water. Wetlands were categorized with r

New Jersey↗

Hydrologic budget of the Harney Basin groundwater system, southeastern Oregon

Groundwater-level declines and limited quantitative knowledge of the groundwater-flow system in the Harney Basin prompted a cooperative study between the U.S. Geological Survey and the Oregon Water Resources Department to evaluate the groundwater-flow system and budget. This report provides a hydrologic budget of the Harney Basin groundwater system that includes separate groundwater budgets for upland and lowland areas to avoid double counting water that recharges in the uplands, discharges to streams and springs in the uplands, flows downstream to the lowlands, and recharges the lowland groundwater system. Lowlands generally represent the conterminous valleys within the center of the basin, including floodplains of the major streams and uplands represent all other areas in the basin. The upland groundwater budget is minimally affected by groundwater development and generally represents the budget of the natural system. In upland areas during 1982–2016, mean-annual recharge totaled 288,000 acre-feet (acre-ft) and mean-annual discharge totaled 239,000 acre-ft, resulting in a net recharge of 49,000 acre-ft. Upland groundwater recharge occurs as infiltration of precipitation and snowmelt and was estimated using the USGS Soil-Water-Balance model calibrated to estimates of runoff, evapotranspiration (ET), base flow, and snow-water equivalent. Groundwater discharge to streams is the predominant discharge mechanism in upland areas and was estimated as 225,000 acre-feet per year (acre-ft/yr) during 1982–2016 using hydrograph separation and summer low-flow estimates in streamgaged watersheds and a linear relation between estimated streamflow and base flow in ungaged watersheds. The remaining upland discharge occurs through springs (14,000 acre-ft/yr) that either emerge downgradient of locations where groundwater discharge to streams was estimated or are routed to irrigated areas. Spring discharge was estimated as a compilation of current and historical measurements. The net upland recharge, which is 17 percent of total upland recharge, ultimately recharges lowland areas as groundwater flow from uplands to lowlands. The lowland groundwater budget for the Harney Basin represents a combination of natural conditions and human activity as more than 99 percent of groundwater development has occurred either inside or within 2 miles of the lowland boundary. In lowland areas during 1982–2016, mean annual groundwater recharge totaled 173,000 acre-ft and groundwater discharge totaled 283,000 acre-ft, indicating discharge exceeded recharge by more than 60 percent. Excluding groundwater pumping, the lowland groundwater budget is more in balance with a mean annual recharge of 165,000 acre-ft and a mean annual discharge of 131,000 acre-ft during 1982–2016. The 23-percent difference between non-pumping recharge and discharge mostly represents the cumulative uncertainty in the estimates of the various groundwater budget components but also likely includes a small reduction in natural groundwater discharge captured by pumping. Lowland groundwater is predominantly recharged by infiltration of surface water (116,000 acre-ft/yr) through streams, floodwater, and irrigation, with a lesser amount as groundwater inflow from uplands and minimal recharge beneath Malheur and Harney Lakes. Recharge from streams and floodwater (natural and irrigation) was estimated using a balance of measured and estimated surface-water inflow to and outflow from lowland areas including streamflow, springflow, and ET where a portion of surface-water inflow to lowland areas is comprised of upland discharge to streams and springs. Groundwater ET (119,000 acre-ft/yr) is the predominant natural discharge mechanism in lowland areas and was estimated as the mean from two remote-sensing based approaches incorporating groundwater ET measurements from other similar basins and 23 years (1987–2015) of Landsat imagery. Discharge of lowland groundwater into Malheur and Harney Lakes is about 700 acre-ft/yr and is represented in groundwater ET estimates. The remaining natural groundwater discharge from lowland areas issues from Sodhouse Spring (8,900 acre-ft/yr) and as groundwater flow to the Malheur River Basin through Virginia Valley (3,100 acre-ft/yr). The relatively large amount of groundwater discharged to springs in Warm Springs Valley (25,000 acre-ft/yr) is accounted for in groundwater ET estimates. Natural groundwater discharge in lowland areas of the Harney Basin has remained relatively constant during the last 80 years based on comparisons with estimates north of Malheur Lake and west of Harney Lake published in the 1930s. Annual net amount of groundwater pumped (pumpage) from the Harney Basin during 2017–18 averaged 144,000 acre-ft. The net value is the difference between pumpage (about 152,000 acre-ft/yr) and reinfiltration of groundwater pumped for irrigation and non-irrigation purposes (about 8,000 acre-ft/yr). Net pumpage was estimated in concurrent studies that compiled groundwater-use data and coupled reported groundwater pumpage data from wells with remote-sensing-based ET estimates from groundwater-irrigated fields. Total pumpage for irrigation has increased from about 54,000 acre-ft/yr during 1991–92 to 145,000 acre-ft/yr during 2017–18. Presently, pumpage is greatest in the lowland region north of Malheur Lake (81,000 acre-ft/yr), with lesser amounts to the north and northwest of Harney Lake (41,000 acre-ft/yr) and to the south and east of Malheur Lake (22,000 acre-ft/yr). During this study, mean annual lowland groundwater discharge (including pumpage) exceeded mean annual recharge, indicating that the lowland hydrologic budget is out of balance. Net groundwater pumpage during 2017–18 is similar to groundwater discharge from all other sources in the lowlands and is four times the imbalance between non-pumping lowland recharge and discharge (34,000 acre-ft/yr). Declining groundwater levels at depth across many parts of the Harney Basin lowlands indicate that pumpage is depleting aquifer storage and is likely capturing a small amount of natural groundwater discharge to springs and ET in some lowland areas. If pumping continues, aquifer storage depletion will continue until the capture rate of natural discharge to springs and ET is equal to the pumping rate. If groundwater development occurs in upland areas and reduces either the streamflow or groundwater inflow to lowland areas, the deficit in the lowland water budget will increase.

Oregon↗

Methods for estimating the magnitude and frequency of peak streamflows for unregulated streams in Oklahoma developed by using streamflow data through 2017

The U.S. Geological Survey (USGS), in cooperation with the Oklahoma Department of Transportation, updated peak-streamflow regression equations for estimating flows with annual exceedance probabilities from 50 to 0.2 percent for the State of Oklahoma. These regression equations incorporate basin characteristics to estimate peak-streamflow magnitude and frequency throughout the State by use of a generalized least-squares regression analysis. The most statistically significant independent variables required to estimate peak-streamflow magnitude and frequency for unregulated streams in Oklahoma are contributing drainage area, mean-annual precipitation, and main-channel slope. The regression equations are applicable for stream basins with drainage areas less than 2,510 square miles that are not affected by regulation. The standard model error ranged from 31.28 to 49.32 percent for the different annual exceedance probabilities that were computed. Annual-maximum peak flows observed at 212 USGS streamgages through water year 2017 were used for the regression analysis, excluding the Oklahoma Panhandle region. The USGS StreamStats web application was used to obtain the independent variables required for the peak-streamflow regression equations. Limitations on the use of the regression equations and the reliability of regression estimates for natural unregulated streams are described. Log-Pearson Type III analysis information, basin and climate characteristics, and the peak-streamflow frequency estimates for the 212 streamgages in and near Oklahoma are provided in this report. This report contains descriptions of the methods that can be used to estimate peak streamflows at ungaged sites by using estimates from streamgages on unregulated streams. For ungaged sites on urban streams and streams regulated by small floodwater-retarding structures, an adjustment of the statewide regression equations for natural unregulated streams can be used to estimate peak-streamflow magnitude and frequency.

Arkansas, Kansas, Missouri, Oklahoma, Texas↗

Assessment of natural gas pipeline construction on stream temperature and turbidity in southwestern Virginia, 2017–25

Despite the extensive natural gas pipeline network in the United States that intersects streams and other sensitive habitats, few case studies use a comparative upstream-to-downstream approach to evaluate potential short- and long-term effects of pipeline stream crossings from pre-construction through post-restoration. In 2017, the U.S. Geological Survey, in cooperation with the Virginia Department of Environmental Quality, deployed real-time continuous stream monitoring stations upstream and downstream from six proposed Mountain Valley Pipeline stream crossings in southwestern Virginia. Water temperature and turbidity data collected at the upstream and downstream sites were compared across three periods—before stream crossing construction, during stream crossing construction, and after stream crossing construction—to determine potential influences from the pipeline stream crossing. Additionally, the monitoring network was used to notify regulators of potentially anomalous conditions throughout the entire monitoring period. The results of this study indicate that pipeline stream crossing did not affect long-term or short-term upstream-to-downstream water temperature conditions or long-term upstream-to-downstream turbidity conditions in any of the six monitored streams. Some short-term anomalously elevated turbidity conditions were observed and attributable to pipeline stream crossing; however, the magnitudes and durations were not sufficient to alter the long-term turbidity regimes of the streams in which they were observed. The application of the monitoring network as a real-time alert system successfully alerted regulators to potentially anomalous conditions.

Virginia↗

The U.S. Geological Survey Geologic Collections Management System (GCMS)—A master catalog and collections management plan for U.S. Geological Survey geologic samples and sample collections

**Updated guidance is available in USGS Instructional Memorandum CSS 2019-01 .** Abstract The U.S. Geological Survey (USGS) is widely recognized in the earth science community as possessing extensive collections of earth materials collected by research personnel over the course of its history. In 2006, a Geologic Collections Inventory was conducted within the USGS Geology Discipline to determine the extent and nature of its sample collections, and in 2008, a working group was convened by the USGS National Geologic and Geophysical Data Preservation Program to examine ways in which these collections could be coordinated, cataloged, and made available to researchers both inside and outside the USGS. The charge to this working group was to evaluate the proposition of creating a Geologic Collections Management System (GCMS), a centralized database that would (1) identify all existing USGS geologic collections, regardless of size, (2) create a virtual link among the collections, and (3) provide a way for scientists and other researchers to obtain access to the samples and data in which they are interested. Additionally, the group was instructed to develop criteria for evaluating current collections and to establish an operating plan and set of standard practices for handling, identifying, and managing future sample collections. Policies and procedures promoted by the GCMS would be based on extant best practices established by the National Science Foundation and the Smithsonian Institution. The resulting report—USGS Circular 1410, “The U.S. Geological Survey Geologic Collections Management System (GCMS): A Master Catalog and Collections Management Plan for U.S. Geological Survey Geologic Samples and Sample Collections”—has been developed for sample repositories to be a guide to establishing common practices in the collection, retention, and disposal of geologic research materials throughout the USGS. While constructing this report, the GCMS’s potential customers and their needs were considered. Two critical definitions have been clarified: a repository is a facility for the long-term management of geologic collections, and a collection is a set of specimens that have been brought together on the basis of some common characteristic. Required sample metadata for newly collected samples, as well as for older collections, were also stipulated and are listed and explained in this report. Several basic policies are also recommended by the GCMS in order to standardize operations among the physical repositories where collections are currently housed. The GCMS Collection Management Plan provides a set of protocols and templates for the management of scientific collections, including access, storage, transfer, and disposal of physical geologic samples and data. This plan is flexible to allow each repository to adapt the practices best suited to its collections. The general consideration for implementation of the GCMS is that all active USGS geologic sample repositories will form the core of GCMS and that participating science centers will develop procedures based on proposed GCMS methodologies. The GCMS is a collective resource for the entire USGS community and the users who discover the geologic materials kept in these repositories and seek to access them.

Circular↗

Analysis of salinity intrusion in the Waccamaw River and Atlantic Intracoastal Waterway near Myrtle Beach, South Carolina, 1995-2002

Six reservoirs in North Carolina discharge into the Pee Dee River, which flows 160 miles through South Carolina to the coastal communities near Myrtle Beach, South Carolina. During the Southeast's record-breaking drought from 1998 to 2003, salinity intrusions inundated a coastal municipal freshwater intake, limiting water supplies. To evaluate the effects of regulated flows of the Pee Dee River on salinity intrusion in the Waccamaw River and Atlantic Intracoastal Waterway, the South Carolina Department of Natural Resources and a consortium of stakeholders entered into a cooperative agreement with the U.S. Geological Survey to apply data-mining techniques to the long-term time series to analyze and simulate salinity dynamics near the freshwater intakes along the Grand Strand of South Carolina. Salinity intrusion in tidal rivers results from the interaction of three principal forces—streamflow, mean tidal water levels, and tidal range. To analyze, model, and simulate hydrodynamic behaviors at critical coastal gages, data-mining techniques were applied to over 20 years of hourly streamflow, coastal water-quality, and water-level data. Artificial neural network models were trained to learn the variable interactions that cause salinity intrusions. Streamflow data from the 18,300-square-mile basin were input to the model as time-delayed variables and accumulated tributary inflows. Tidal inputs to the models were obtained by decomposing tidal water-level data into a "periodic" signal of tidal range and a "chaotic" signal of mean water levels. The artificial neural network models were able to convincingly reproduce historical behaviors and generate alternative scenarios of interest. To make the models directly available to all stakeholders along the Pee Dee and Waccamaw Rivers and Atlantic Intracoastal Waterway, an easy-to-use decision support system (DSS) was developed as a spreadsheet application that integrates the historical database, artificial neural network models, model controls, streaming graphics, and model output. An additional feature is a built-in optimizer that dynamically calculates the amount of flow needed to suppress salinity intrusions as tidal ranges and water levels vary over days and months. This DSS greatly reduced the number of long-term simulations needed for stakeholders to determine the minimum flow required to adequately protect the freshwater intakes.

South Carolina↗

Age and quality of ground water and sources of nitrogen in the surficial aquifers in Pumpkin Creek Valley, western Nebraska, 2000

Ground water is the source of drinking water for the residents of Pumpkin Creek Valley, western Nebraska. In this largely agricultural area, shallow aquifers potentially are susceptible to nitrate contamination. During the last 10 years, ground-water levels in the North Platte Natural Resources District have declined and contamination has become a major problem for the district. In 2000, the U.S. Geological Survey and the North Platte Natural Resources District began a cooperative study to determine the age and quality of the ground water and the sources of nitrogen in the aquifers in Pumpkin Creek Valley. Water samples were collected from 8 surface-water sites, 2 springs, and 88 ground-water sites during May, July, and August 2000. These samples were analyzed for physical properties, nutrients or nitrate, and hydrogen and oxygen isotopes. In addition, a subset of samples was analyzed for any combination of chlorofluorocarbons, tritium, tritium/helium, sulfur-hexafluoride, carbon-14, and nitrogen-15. The apparent age of ground water in the alluvial aquifer typically varied from about 1980 to modern, whereas ground water in the fractured Brule Formation had a median value in the 1970s. The Brule Formation typically contained ground water that ranged from the 1940s to the 1990s, but low-yield wells had apparent ages of 5,000 to 10,000 years before present. Data for oxygen-18 and deuterium indicated that lake-water samples showed the greatest effects from evaporation. Ground-water data showed no substantial evaporative effects and some ground water became isotopically heavier as the water moved downgradient. In addition, the physical and chemical ground-water data indicate that Pumpkin Creek is a gaining stream because little, if any, of its water is lost to the ground-water system. The water-quality type changed from a sodium calcium bicarbonate type near Pumpkin Creek's headwaters to a calcium sodium bicarbonate type near its mouth. Nitrate concentrations were largest in the alluvial system (median = 5 mg/L) and smallest in the surface-water system (median = 1 mg/L). Most nitrate concentrations exceeding the U.S. Environmental Protection Agency maximum contaminant level for drinking water of 10 mg/L as nitrogen were adjacent to irrigated fields and in areas where alluvial sediments are less than 50 ft thick. Sources of nitrogen in the ground water of the study area included naturally occurring nitrogen, commercial fertilizer, and animal waste. Based on nitrate concentration and delta nitrogen-15, the nitrogen in 65 percent of the water samples appears to have originated from a mixture of commercial fertilizers and animal waste. Some of the smallest nitrate concentrations in the ground-water samples contained some of the largest delta nitrogen-15 values (greater than 10 per mil), which suggests animal waste as the likely source. Commercial fertilizers were the likely source of most of the nitrogen in water samples with nitrate concentrations that exceeded 10 mg/L. The source of the nitrogen in water samples with nitrate concentrations exceeding 10 mg/L, but with delta nitrogen-15 values close to 10 per mil, could not be determined.

Nebraska↗