USGS Science⌕ Search

SEARCH · USGS Science

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

Occurrence of Selected Nutrients, Trace Elements, and Organic Compounds in Streambed Sediment in the Lower Chena River Watershed near Fairbanks, Alaska, 2002-03

In 2002-03, the U.S. Geological Survey collected samples of streambed sediment at 18 sites in the lower Chena River watershed for analysis of selected nutrients, traces elements, and organic compounds. The purpose of the project was to provide Federal, State, and local agencies as well as neighborhood committees, with information for consideration in plans to improve environmental conditions in the watershed. The exploratory sampling program included analysis of streambed sediment from the Chena River and Chena Slough, a tributary to the Chena River. Results were compared to streambed-sediment guidelines for the protection of aquatic life and to 2001-02 sediment data from Noyes Slough, a side channel of the lower Chena River. The median total phosphorus concentration in Chena Slough sediment samples, 680 milligrams per kilogram (mg/kg), was two orders of magnitude greater than median total phosphorus concentration in Chena River sediment samples of 5.2 mg/kg. Median concentrations of chloride and sulfate also were greater in Chena Slough samples. Low concentrations of nitrate were detected in most of the Chena Slough samples; nitrate concentrations were below method reporting limits or not detected in Chena River sediment samples. Streambed-sediment samples were analyzed for 24 trace elements. Arsenic, nickel, and zinc were the only trace elements detected in concentrations that exceeded probable-effect levels for the protection of aquatic life. Concentrations of arsenic in Chena Slough samples ranged from 11 to 70 mg/kg and concentrations in most of the samples exceeded the probable-effect guideline for arsenic of 17 mg/kg. Arsenic concentrations in samples from the Chena River ranged from 9 to 12 mg/kg. The background level for arsenic in the lower Chena River watershed is naturally elevated because of significant concentrations of arsenic in local bedrock and ground water. Sources of elevated concentrations of zinc in one sample, and of nickel in two samples, are unknown. With the exception of elevated arsenic levels in samples from Chena Slough, the occurrence and concentration of trace elements in the streambed sediments of Chena Slough and Chena River were similar to those in Noyes Slough sediment. Sediment samples were analyzed for 78 semivolatile organic compounds and 32 organochlorine pesticides and polychlorinated biphenyls (PCBs). Low concentrations of dimethylnaphthalene and p-Cresol were detected in most Chena Slough and Chena River sediment samples. The number of semivolatile organic compounds detected ranged from 5 to 21 in most Chena Slough sediment samples. In contrast, three or fewer semivolatile organic compounds were detected in Chena River sediment samples, most likely because chemical-matrix interference resulted in elevated reporting limits for organochlorine compounds in the Chena River samples. Low concentrations of fluoranthene, pyrene, and phenanthrene were detected in Chena Slough sediment. Relatively low concentrations of DDT or its degradation products, DDD and DDE, were detected in all Chena Slough samples. Concentrations of total DDT (DDT+DDD+DDE) in two Chena Slough sediment samples exceeded the effectsrange median aquatic-life criteria of 46.1 micrograms per kilogram (ug/kg). DDT concentrations in Chena River streambed-sediment samples were less than 20 ug/kg. Low concentrations of PCB were detected in two Chena Slough streambed-sediment samples. None of the concentrations of the polychlorinated biphenyls or semivolatile organic compounds for which the samples were analyzed exceeded available guidelines for the protection of aquatic life. With the exception of elevated total DDT in two Chena Slough samples, the occurrence and concentration of organochlorine compounds in Chena Slough and Chena River sediment were similar to those in samples collected from Noyes Slough in 2001-02.

Scientific Investigations Report↗

Hydrology-driven chemical loads transported by the Green River to the Lower Duwamish Waterway near Seattle, Washington, 2013–17

The sediments in the Lower Duwamish Waterway Superfund site in Seattle, Washington, are contaminated with chemicals including metals such as arsenic, polychlorinated biphenyls (PCBs), carcinogenic polycyclic aromatic hydrocarbons (cPAHs), and dioxins/furans from decades of intense anthropogenic activities. The U.S. Geological Survey, in cooperation with the Washington State Department of Ecology, collected new data from 2013 to 2017 to estimate sediment and chemical loads transported by the Green/Duwamish River to the Lower Duwamish Waterway Superfund site (the final 8-kilometer reach of the river) in support of sediment remediation within the site. Chemical loads were calculated as the product of river suspended-sediment loads and suspended sediment-bound chemical concentrations measured at river kilometer 16.7. Using four different approaches, annual suspended sediment-bound chemical load estimates transported by the river to the Lower Duwamish Waterway were in the range of 1,120–1,470 kilograms arsenic, 2,810–8,200 grams (g) toxic equivalent cPAHs, 205–407 milligrams toxic equivalent dioxins/furans, and 340–1,180 g PCBs. Storm events contributed a disproportionately large amount of the load of anthropogenic organic compounds such as cPAHs (54 percent), dioxins/furans (44 percent), and PCBs (52 percent) as compared to overall time (17 percent). Chemical concentrations and load estimates often were underestimated using results from unfiltered water analysis only, especially in samples with high suspended-sediment concentrations and for hydrophobic organic chemicals such as cPAHs that prefer to sorb to particulates and are at low concentrations near or below the analytical limits of water methods. For metals and PCBs, the dissolved concentration was relatively low and consistent between sampling events, whereas the suspended sediment-bound chemical concentrations contributed most of the chemical concentration in the water column during periods of high river suspended-sediment concentrations. However, the dissolved fraction, on average, contributed more than one-third of the estimated total chemical load in the river system for arsenic and PCBs, even given the hydrophobic nature of the chemicals. These results suggest that the sum of the chemical concentrations measured on two separate fractions—the particulate fraction and the dissolved fraction—more fully represents the total chemical concentration as compared to analysis of an unfiltered water sample, especially in samples with high suspended-sediment concentrations. Most of the suspended-sediment load (97 percent) and sediment-bound chemical load (92–94 percent) occurred during the wet winter half of the year from October 15 to April 14. However, the highest sediment-bound chemical concentrations often occurred during short intense storms or “first flush” autumn runoff events during the dry summer half of the year from April 15 to October 14. Because of the highly variable and dynamic river system characteristics (including precipitation, discharge, sediment concentration, and tidal fluctuations), it is critical to characterize the occurrence, frequency, concentrations, and loads during extreme conditions (for example, when the river is affected by storm-derived runoff) rather than time-averaged conditions. These short extreme events have a high potential for acute effects on ecological and human health, and may have a great influence on the effectiveness of the sediment remediation activities that are underway in the Lower Duwamish Waterway.

Washington↗

Predicting lake trophic state by relating Secchi-disk transparency measurements to Landsat-satellite imagery for Michigan inland lakes, 2003-05 and 2007-08

Inland lakes are an important economic and environmental resource for Michigan. The U.S. Geological Survey and the Michigan Department of Natural Resources and Environment have been cooperatively monitoring the quality of selected lakes in Michigan through the Lake Water Quality Assessment program. Sampling for this program began in 2001; by 2010, 730 of Michigan’s 11,000 inland lakes are expected to have been sampled once. Volunteers coordinated by the Michigan Department of Natural Resources and Environment began sampling lakes in 1974 and continue to sample (in 2010) approximately 250 inland lakes each year through the Michigan Cooperative Lakes Monitoring Program. Despite these sampling efforts, it still is impossible to physically collect measurements for all Michigan inland lakes; however, Landsat-satellite imagery has been used successfully in Minnesota, Wisconsin, Michigan, and elsewhere to predict the trophic state of unsampled inland lakes greater than 20 acres by producing regression equations relating in-place Secchi-disk measurements to Landsat bands. This study tested three alternatives to methods previously used in Michigan to improve results for predicted statewide Trophic State Index (TSI) computed from Secchi-disk transparency (TSI (SDT)). The alternative methods were used on 14 Landsat-satellite scenes with statewide TSI (SDT) for two time periods (2003– 05 and 2007–08). Specifically, the methods were (1) satellitedata processing techniques to remove areas affected by clouds, cloud shadows, haze, shoreline, and dense vegetation for inland lakes greater than 20 acres in Michigan; (2) comparison of the previous method for producing a single open-water predicted TSI (SDT) value (which was based on an area of interest (AOI) and lake-average approach) to an alternative Gethist method for identifying open-water areas in inland lakes (which follows the initial satellite-data processing and targets the darkest pixels, representing the deepest water, before regression equations are created); and (3) checking to see whether the predicted TSI (SDT) values compared well between two regression equations, one previously used in Michigan and an alternative equation from the hydrologic literature. The combination of improved satellite-data processing techniques and the Gethist method to identify open-water areas in inland lakes during 2003–05 and 2007–08 provided a stronger relation and statistical significance between predicted TSI (SDT) and measured TSI than did the AOI lake-average method; differences in results for the two methods were significant at the 99-percent confidence level. With regard to the comparison of the regression equations, there were no statistically significant differences at the 95-percent confidence level between results from the two equations. The previously used equation, in combination with the Gethist method, yielded coefficient of determination (R 2 ) values of 0.71 and 0.77 for the periods 2003–05 and 2007–08, respectively. The alternative equation, in combination with the Gethist method, yielded R 2 values of 0.74 and 0.75 for 2003–05 and 2007–08, respectively. Predicted TSI (SDT) and measured TSI (SDT) values for lakes used in the regression equations compared well, with R 2 values of 0.95 and 0.96 for predicted TSI (SDT) for 2003–05 and 2007–08, respectively. The R 2 values for statewide predicted TSI (SDT) for all inland lakes with available open-water areas for 2003–05 and 2007–08 were 0.91 and 0.93, respectively. Although the two equations predicted similar trophic-state classes, the alternative equation is planned to be used for future prediction of TSI (SDT) values for Michigan inland lakes, to promote consistency in comparing predicted values between States and for potential use in trend analysis.

Michigan↗

Questa baseline and pre-mining ground-water quality investigation. 19. Leaching characteristics of composited materials from mine waste-rock piles and naturally altered areas near Questa, New Mexico

The goal of this study is to compare and contrast the leachability of metals and the acidity from individual mine waste-rock piles and natural erosional scars in the study area near Questa, New Mexico. Surficial multi-increment (composite) samples less than 2 millimeters in diameter from five waste-rock piles, nine erosional-scar areas, a less-altered site, and a tailings slurry-pipe sample were analyzed for bulk chemistry and mineralogy and subjected to two back-to-back leaching procedures. The first leaching procedure, the U.S. Geological Survey Field Leach Test (FLT), is a short-duration leach (5-minute shaking and 10-minute settling) and is intended to leach readily soluble materials. The FLT was immediately followed by an 18-hour, end-over-end rotation leaching procedure. Comparison of results from the back-to-back leaching procedures can provide information about reactions that may take place upon migration of leachates through changing geochemical conditions (for example, pH changes), both within the waste-rock and scar materials and away from the source materials. For the scar leachates, the concentrations of leachable metals varied substantially between the scar areas sampled. The scar leachates have low pH (pH 3.2-4.1). Under these low-pH conditions, cationic metals are solubilized and mobile, but anionic species, such as molybdenum, are less soluble and less mobile. Generally, metal concentrations in the waste-rock leachates did not exceed the upper range of those metal concentrations in the erosional-scar leachates. One exception is molybdenum, which is notably higher in the waste-rock leachates compared with the scar leachates. Most of the waste-rock leachates were at least mildly acidic (pH 3.0-6.2). The pH values in the waste-rock leachates span a large pH range that includes some pH-dependent solubility and metal-attenuation reactions. An increase in pH with leaching time and agitation indicates that there is pH-buffering capacity in some of the waste-rock piles. As pH increased in the waste-pile leachates, concentrations of several metals decreased with increasing time and agitation. Similar pH-dependent reactions may take place upon migration of the leachates in the waste-rock piles. Bulk chemistry, mineralogy, and leachate sulfur-isotope data indicate that the Capulin and Sugar Shack West waste-rock piles are compositionally different from the younger Sugar Shack South, Sugar Shack Middle, and Old Sulphur Gulch piles. The Capulin and Sugar Shack West piles have the lowest-pH leachates (pH 3.0-4.1) of the waste-pile samples, and the source material for the Capulin and Sugar Shack West piles appears to be similar to the source material for the erosional-scar areas. Calcite dissolution, in addition to gypsum dissolution, appears to produce the calcium and sulfate concentrations in leachates from the Sugar Shack South, Sugar Shack Middle, and Old Sulphur Gulch piles.

Scientific Investigations Report↗

Stratiform chromite deposit model

A new descriptive stratiform chromite deposit model was prepared which will provide a framework for understanding the characteristics of stratiform chromite deposits worldwide. Previous stratiform chromite deposit models developed by the U.S. Geological Survey (USGS) have been referred to as Bushveld chromium, because the Bushveld Complex in South Africa is the only stratified, mafic-ultramafic intrusion presently mined for chromite and is the most intensely researched. As part of the on-going effort by the USGS Mineral Resources Program to update existing deposit models for the upcoming national mineral resource assessment, this revised stratiform chromite deposit model includes new data on the geological, mineralogical, geophysical, and geochemical attributes of stratiform chromite deposits worldwide. This model will be a valuable tool in future chromite resource and environmental assessments and supplement previously published models used for mineral resource evaluation. Stratiform chromite deposits are found throughout the world, but the chromitite seams of the Bushveld Complex, South Africa, are the largest and most intensely researched. The chromite ore is located primarily in massive chromitite seams and, less abundantly, in disseminated chromite-bearing layers, both of which occur in the ultramafic section of large, layered mafic-ultramafic stratiform complexes. These mafic-ultramafic intrusions mainly formed in stable cratonic settings or during rift-related events during the Archean or early Proterozoic, although exceptions exist. The chromitite seams are cyclic in nature as well as laterally contiguous throughout the entire intrusion. Gangue minerals include olivine, pyroxenes (orthopyroxene and clinopyroxene), plagioclase, sulfides (pyrite, chalcopyrite, pyrrhotite, pentlandite, bornite), platinum group metals (mainly laurite, cooperite, braggite), and alteration minerals. A few deposits also contain rutile and ilmenite. The alteration phases include serpentine, chlorite, talc, magnetite, kaemmererite, uvarovite, hornblende, and carbonate minerals, such as calcite and dolomite. Stratiform chromite deposits are primarily hosted by peridotites, harzburgites, dunites, pyroxenites, troctolites, and anorthosites. Although metamorphism may have altered the ultramafic regions of layered intrusions postdeposition, only igneous processes are responsible for formation. From a diagnostic standpoint and for assessment purposes, they have no temporal or spatial relation to sedimentary rocks. The exact mechanisms responsible for the development of stratiform chromite deposits and the large, layered mafic-ultramafic intrusions where they are found are highly debated. The leading argument postulates that a parent magma mixed with a more primitive magma during magma chamber recharge. The partially differentiated magma could then be forced into the chromite stability field, resulting in the massive chromitite layers found in stratiform complexes. Contamination of the parent magma by localized assimilation of felsic country rock at the roof of the magma chamber has also been proposed as a mechanism of formation. Others suggest that changes in pressure or oxygen fugacity may be responsible for the occurrence of massive chromitite seams in layered mafic, ultramafic intrusions. The massive chromitite layers contain high levels of chromium and strong associations with platinum group elements. Anomalously high magnesium concentrations as well as low sodium, potassium, and phosphorus concentrations are also important geochemical features of stratiform chromite deposits. The presence of orthopyroxenite in many of the deposits suggests high silica and high magnesium concentrations in the parent magma. Most environmental concerns associated with the mining and processing of chromite ore focus on the solubility of chromium and its oxidation state. Although trivalent chromium (Cr 3+ ) is an essential micronutrient for humans, hexavalent chromium (Cr 6+ ) is highly toxic. Chromium-bearing solid phases that occur in the chromite ore-processing residue, for example, can effect the geochemical behavior and oxidation state of chromium in the environment.

Scientific Investigations Report↗

Organic wastewater compounds, pharmaceuticals, and coliphage in ground water receiving discharge from onsite wastewater treatment systems near La Pine, Oregon: Occurrence and implications for transport

The occurrence of organic wastewater compounds (components of 'personal care products' and other common household chemicals), pharmaceuticals (human prescription and nonprescription medical drugs), and coliphage (viruses that infect coliform bacteria, and found in high concentrations in municipal wastewater) in onsite wastewater (septic tank effluent) and in a shallow, unconfined, sandy aquifer that serves as the primary source of drinking water for most residents near La Pine, Oregon, was documented. Samples from two types of observation networks provided basic occurrence data for onsite wastewater and downgradient ground water. One observation network was a group of 28 traditional and innovative (advanced treatment) onsite wastewater treatment systems and associated downgradient drainfield monitoring wells, referred to as the 'innovative systems network'. The drainfield monitoring wells were located adjacent to or under onsite wastewater treatment system drainfield lines. Another observation network, termed the 'transect network', consisted of 31 wells distributed among three transects of temporary, stainless-steel-screened, direct-push monitoring wells installed along three plumes of onsite wastewater. The transect network, by virtue of its design, also provided a basis for increased understanding of the transport of analytes in natural systems. Coliphage were frequently detected in onsite wastewater. Coliphage concentrations in onsite wastewater were highly variable, ranging from less than 1 to 3,000,000 plaque forming units per 100 milliliters. Coliphage were occasionally detected (eight occurrences) at low concentrations in samples from wells located downgradient from onsite wastewater treatment system drainfield lines. However, coliphage concentrations were below method detection limits in replicate or repeat samples collected from the eight sites. The consistent absence of coliphage detections in the replicate or repeat samples is interpreted to indicate that the detections reported for ground-water samples represented low-level field or laboratory contamination, and it would appear that coliphage were effectively attenuated to less than 1 PFU/100 mL over distances of several feet of transport in the La Pine aquifer and (or) overlying unsaturated zone. Organic wastewater compounds were frequently detected in onsite wastewater. Of the 63 organic wastewater compounds in the analytical schedule, 45 were detected in the 21 samples of onsite wastewater. Concentrations of organic wastewater compounds reached a maximum of 1,300 ug/L (p-cresol). Caffeine was detected at concentrations as high as 320 ug/L. Fourteen of the 45 compounds were detected in more than 90 percent of onsite wastewater samples. Fewer (nine) organic wastewater compounds were detected in ground water, despite the presence of nitrate and chloride likely from onsite wastewater sources. The nine organic wastewater compounds that were detected in ground-water samples were acetyl-hexamethyl-tetrahydro-naphthalene (AHTN), caffeine, cholesterol, hexahydrohexamethyl-cyclopentabenzopyran, N,N-diethyl-meta-toluamide (DEET), tetrachloroethene, tris (2-chloroethyl) phosphate, tris (dichloroisopropyl) phosphate, and tributyl phosphate. Frequent detection of household-chemical type organic wastewater compounds in onsite wastewater provides evidence that some of these organic wastewater compounds may be useful indicators of human waste effluent dispersal in some hydrologic environments. The occurrence of organic wastewater compounds in ground water downgradient from onsite wastewater treatment systems demonstrates that a subgroup of organic wastewater compounds is transported in the La Pine aquifer. The consistently low concentrations (generally less than 1 ug/L) of organic wastewater compounds in water samples collected from wells located no more than 19 feet from drainfield lines indicates that the reactivity (sorption, degradation) of this suite of organic waste.

Oregon↗

Characteristics and Classification of Least Altered Streamflows in Massachusetts

Streamflow records from 85 streamflow-gaging stations at which streamflows were considered to be least altered were used to characterize natural streamflows within southern New England. Period-of-record streamflow data were used to determine annual hydrographs of median monthly flows. The shapes and magnitudes of annual hydrographs of median monthly flows, normalized by drainage area, differed among stations in different geographic areas of southern New England. These differences were gradational across southern New England and were attributed to differences in basin and climate characteristics. Period-of-record streamflow data were also used to analyze the statistical properties of daily streamflows at 61 stations across southern New England by using L-moment ratios. An L-moment ratio diagram of L-skewness and L-kurtosis showed a continuous gradation in these properties between stations and indicated differences between base-flow dominated and runoff-dominated rivers. Streamflow records from a concurrent period (1960-2004) for 61 stations were used in a multivariate statistical analysis to develop a hydrologic classification of rivers in southern New England. Missing records from 46 of these stations were extended by using a Maintenance of Variation Extension technique. The concurrent-period streamflows were used in the Indicators of Hydrologic Alteration and Hydrologic Index Tool programs to determine 224 hydrologic indices for the 61 stations. Principal-components analysis (PCA) was used to reduce the number of hydrologic indices to 20 that provided nonredundant information. The PCA also indicated that the major patterns of variability in the dataset are related to differences in flow variability and low-flow magnitude among the stations. Hierarchical cluster analysis was used to classify stations into groups with similar hydrologic properties. The cluster analysis classified rivers in southern New England into two broad groups: (1) base-flow dominated rivers, whose statistical properties indicated less flow variability and high magnitudes of low flow, and (2) runoff-dominated rivers, whose statistical properties indicated greater flow variability and lower magnitudes of low flow. A four-cluster classification further classified the runoff-dominated streams into three groups that varied in gradient, elevation, and differences in winter streamflow conditions: high-gradient runoff-dominated rivers, northern runoff-dominated rivers, and southern runoff-dominated rivers. A nine-cluster division indicated that basin size also becomes a distinguishing factor among basins at finer levels of classification. Smaller basins (less than 10 square miles) were classified into different groups than larger basins. A comparison of station classifications indicated that a classification based on multiple hydrologic indices that represent different aspects of the flow regime did not result in the same classification of stations as a classification based on a single type of statistic such as a monthly median. River basins identified by the cluster analysis as having similar hydrologic properties tended to have similar basin and climate characteristics and to be in close proximity to one another. Stations were not classified in the same cluster on the basis of geographic location alone; as a result, boundaries cannot be drawn between geographic regions with similar streamflow characteristics. Rivers with different basin and climate characteristics were classified in different clusters, even if they were in adjacent basins or upstream and downstream within the same basin.

Scientific Investigations Report↗

Geophysical bed sediment characterization of the Androscoggin River from the former Chlor-Alkali Facility Superfund Site, Berlin, New Hampshire, to the state border with Maine, August 2009

The former Chlor-Alkali Facility in Berlin, New Hampshire, was listed on the U.S. Environmental Protection Agency National Priorities List in 2005 as a Superfund site. The Chlor-Alkali Facility lies on the east bank of the Androscoggin River. Elemental mercury currently discharges from that bank into the Androscoggin River. The nature, extent, and the speciation of mercury and the production of methyl mercury contamination in the adjacent Androscoggin River is the subject of continuing investigations. The U.S. Geological Survey, in cooperation with Region I of the U.S. Environmental Protection Agency, used geophysical methods to determine the distribution, thickness, and physical properties of sediments in the Androscoggin River channel at a small area of an upstream reference reach and downstream from the site to the New Hampshire–Maine State border. Separate reaches of the Androscoggin River in the study area were surveyed with surface geophysical methods including ground-penetrating radar and step-frequency electromagnetics. Results were processed to assess sediment characteristics including grain size, electrical conductivity, and pore-water specific conductance. Specific conductance measured during surface- and pore-water sampling was used to help interpret the results of the geophysical surveys. The electrical resistivity of sediment samples was measured in the laboratory with intact pore water for comparison with survey results. In some instances, anthropogenic features and land uses, such as roads and power lines affected the detection of riverbed properties using geophysical methods; when this occurred, the data were removed. Through combining results, detailed riverbed sediment characterizations were made. Results from ground-penetrating radar surveys were used to image and measure the depth to the riverbed, depth to buried riverbeds, riverbed thickness and to interpret material-type variations in terms of relative grain size. Fifty two percent of the riverbed in the study area was covered with gravel and finer sediments. The electrically resistive river water and sediment in this study area were conducive to the penetration of the ground-penetrating radar and step-frequency electromagnetic signals and allowed for effective sediment characterization by geophysical methods. The reach between the former Chlor-Alkali Facility and the Riverside Dam, had small areas of fine sediment (estimated 11 percent of riverbed area), found on the upstream left bank and the downstream right bank, with an electromagnetic conductivity (31.4 millisiemens per meter (mS/m) maximum) that was higher than the upstream reference reach. The greatest electromagnetic conductivity (195 mS/m), pore-water specific conductance (324 mS/m) and lab measured sediment conductivity of (76.8 mS/m, measured with a direct-current resistivity test box) in the study were measured approximately 1 mile (mi) downstream of the site from a sandbar on the left bank. Reaches adjacent to and within 2 mi downstream from the site had elevated electromagnetic conductivity despite having lower estimated percentages of riverbed area covered in sediment (11, 25, and 61 percent, respectively) than the reference reach (97). Typically finer grained sediment with similar mineralogy will be more conductive. The Shelburne Reservoir is approximately 8 mi downstream from the site had the second greatest pore-water specific conductance measured, 45.8 mS/m. Many of the locations with the largest step-frequency electromagnetic values have not been sampled for pore water and sediment.

New Hampshire↗

Development and analysis of regional curves for streams in the non-urban valley and ridge physiographic province, Maryland, Virginia, and West Virginia

Regression relations for bankfull stream characteristics based on drainage area (often called 'regional curves') are used in natural stream channel design to verify field determinations of bankfull discharge and stream channel characteristics. Bankfull stream characteristics were assessed for stream reaches at 41 streamflow-gaging stations in the Valley and Ridge Physiographic Province in Maryland, Virginia, and West Virginia. Data collected included bankfull cross-sectional geometry, flood plain geometry, and longitudinal profile data. In addition, particle-size distributions of streambed material were determined and data on basin characteristics were compiled for each reach. Regional curves were developed for bankfull cross-sectional area, width, and discharge with R2 values of 0.95, 0.89, 0.87, and 0.91, respectively. Examination of the regional curves residuals indicates that there is more variability in bankfull cross-sectional area, width, and discharge for smaller streams than for larger streams. In contrast, there is more variability for bankfull mean depth for larger streams than for smaller streams. Geographic analysis of regional curve residuals indicated that there were no further subdivisions within the Valley and Ridge Physiographic Province in the three-state study area for which individual sets of regional curves should be developed. In addition, two separate sets of regional curves were developed with data from the 41 sites to examine potential differences in the relations between the southern (n = 9) and central (n = 32) sections of the province. There were differences in slope and intercept between the two bankfull discharge test relations and a difference in intercept for the width test relations at the 95-percent confidence level. However, the results of this analysis were inconclusive and therefore one set of regional curves for the study area is presented in this report. The regional curves were compared to regression models developed from similar data collected in the Pennsylvania and Maryland portions of the province. No statistical difference in the slope or intercept of regression lines of the three data sets was detected for any of the four bankfull parameters at the 95-percent confidence level. Basin characteristics such as percentage of basin forested (percent forested) and percentage of basin underlain by carbonate bedrock (percent carbonate) were analyzed to evaluate variability among regression points. Multivariate regression relations including explanatory terms for percent carbonate and drainage area produced higher R2 values than the regional curves for bankfull cross-sectional area (R2 = 0.95), bankfull width (R2 = 0.92), and bankfull discharge (R2 = 0.93). There was no improvement for the bankfull mean depth relation from adding the additional term. Inclusion of the other basin characteristics in multivariate relations did not improve the regression models. Regression models developed for the 1.5-year discharge for all streamflow-gaging stations with peak discharge data throughout Virginia (n = 486) and throughout the Valley and Ridge Physiographic Province in Virginia (n = 147) were compared to the regional curve relating bankfull discharge to drainage area. A similar trend in decreasing variability with increasing drainage area was observed for the 1.5-year discharge for all stations in Virginia . This indicates that the change in variability observed in the discharge regional curve likely would exist with a larger data set. There was no statistical difference at the 95-percent confidence level between regression relations for the southern section of the province (n = 40) and the central section (n = 107). This finding supports maintaining only one set of regional curves for the study area. Not all of the variability in the regional curves is explained by drainage area alone. Causes of the remaining variability likely vary among study sites. Users of the regional curves developed in this study are cautioned that, because of inherent variability in hydrologic data, a thorough analysis of the fluvial system and the drainage basin is necessary when evaluating the bankfull characteristics of a stream.

Maryland, Virginia, West Virginia↗

Water resources of Sweetwater County, Wyoming

Sweetwater County is located in the southwestern part of Wyoming and is the largest county in the State. A study to quantify the availability and describe the chemical quality of surface-water and ground-water resources in Sweetwater County was conducted by the U.S. Geological Survey in cooperation with the Wyoming State Engineer’s Office. Most of the county has an arid climate. For this reason a large amount of the flow in perennial streams within the county is derived from outside the county. Likewise, much of the ground-water recharge to aquifers within the county is from flows into the county, and occurs slowly. Surface-water data were not collected as part of the study. Evaluations of streamflow and stream-water quality were limited to analyses of historical data and descriptions of previous investigations. Forty-six new ground-water-quality samples were collected as part of the study and the results from an additional 782 historical ground-water-quality samples were reviewed. Available hydrogeologic characteristics for various aquifers throughout the county also are described. Flow characteristics of streams in Sweetwater County vary substantially depending on regional and local basin characteristics and anthropogenic factors. Because precipitation amounts in the county are small, most streams in the county are ephemeral, flowing only as a result of regional or local rainfall or snowmelt runoff. Flows in perennial streams in the county generally are a result of snowmelt runoff in the mountainous headwater areas to the north, west, and south of the county. Flow characteristics of most perennial streams are altered substantially by diversions and regulation. Water-quality characteristics of selected streams in and near Sweetwater County during water years 1974 through 1983 were variable. Concentrations of dissolved constituents, suspended sediment, and bacteria generally were smallest at sites on the Green River because of resistant geologic units, increased vegetative cover, large diluting streamflows, and large reservoirs. Concentrations of dissolved constituents, suspended sediment, and bacteria generally were largest at sites in the Big Sandy River and Bitter Creek Basins. Some nutrient concentrations and bacteria counts exceeded various State and Federal water-quality criteria. Historical and recent anthropogenic activities contributed to natural sources of many dissolved constituents and suspended sediment. Both water-table and artesian conditions occur in aquifers within the county. Shallow ground water is available throughout the county, although much of it is only marginally suitable or is unsuitable for domestic and irrigation uses mainly because of high total dissolved solids (TDS) concentrations. Suitable ground water for livestock use can be found in most areas of the county. Ground-water quality tends to deteriorate with increasing distance from recharge areas and with increasing depth below land surface. Ground water from depths of greater than a few thousand feet tends to have TDS concentrations that make it moderately saline to briny. In some areas even shallow ground water has moderately saline TDS concentrations. Specific constituents in parts of some aquifers in the county occur in relatively high concentrations when compared to U.S. Environmental Protection Agency drinking-water standards. Relatively high concentrations of sulfate, fluoride, boron, iron, and manganese were found in several aquifers. Many ground-water samples from the Battle Spring aquifer in the Great Divide Structural Basin had high radionuclide concentrations. The estimated mean daily water use in Sweetwater County in 2000 was 170.73 million gallons per day. Irrigation was the largest single use of water in the county with an estimated mean use of more than 92 million gallons per day. Surface water irrigation accounted for nearly 90 percent of the total irrigation water used in 2000. Although ground water is used to a much lesser extent than surface water, in many areas of the county it is the only source of water available.

Wyoming↗

Potential effects of alterations to the hydrologic system on the distribution of salinity in the Biscayne aquifer in Broward County, Florida

To address concerns about the effects of water-resource management practices and rising sea level on saltwater intrusion, the U.S. Geological Survey in cooperation with the Broward County Environmental Planning and Community Resilience Division, initiated a study to examine causes of saltwater intrusion and predict the effects of future alterations to the hydrologic system on salinity distribution in eastern Broward County, Florida. A three-dimensional, variable-density solute-transport model was calibrated to conditions from 1970 to 2012, the period for which data are most complete and reliable, and was used to simulate historical conditions from 1950 to 2012. These types of models are typically difficult to calibrate by matching to observed groundwater salinities because of spatial variability in aquifer properties that are unknown, and natural and anthropogenic processes that are complex and unknown; therefore, the primary goal was to reproduce major trends and locally generalized distributions of salinity in the Biscayne aquifer. The methods used in this study are relatively new, and results will provide transferable techniques for protecting groundwater resources and maximizing groundwater availability in coastal areas. The model was used to (1) evaluate the sensitivity of the salinity distribution in groundwater to sea-level rise and groundwater pumping, and (2) simulate the potential effects of increases in pumping, variable rates of sea-level rise, movement of a salinity control structure, and use of drainage recharge wells on the future distribution of salinity in the aquifer. Results from the simulation of historical conditions indicate that the model generally represents the observed greater westward extent of elevated salinity in the central part of the intruded area relative to the northern and southernmost parts of the intruded area. Results of sensitivity testing indicate that the extent of elevated salinity is most sensitive to pumping in areas where the source of saltwater is largely offshore, from the Atlantic Ocean, and is most sensitive to sea-level rise in areas where the source of salinity is downward leakage of brackish water from canals. Simulations of future scenarios indicate that increases in pumping near the existing interface may cause the interface to advance and decreases in pumping may cause it to retreat. Climatic effects, such as periods of prolonged drought or high precipitation, may augment or counteract long-term effects of changes in pumping on aquifer salinity at well fields. With increasing rates of sea-level rise, the freshwater-saltwater interface advances progressively inland, and flow-averaged salinities at well fields near the existing interface increase commensurately. Hypothetical southeastward (downstream) re-positioning of the existing G–54 salinity-control structure may prevent the interface from moving northwestward along and near the North New River canal, but beneficial effects are localized. Implementation of freshwater recharge wells in the city of Hallandale Beach may also have only a localized freshening effect in the aquifer and little appreciable effect on the freshwater-saltwater interface or on concentrations of salinity at well fields. Model accuracy and use are limited by uncertainty in the physical properties and boundary conditions of the system, uncertainty in historical and future conditions, and generalizations made in the mathematical relationships used to describe the physical processes of groundwater flow and transport. Because of these limitations, model results should be considered in relative rather than absolute terms. Nonetheless, model results do provide useful information on the relative scale of response of the system to changes in pumping distribution, sea-level rise, and mitigation activities.

Florida↗

Bathymetric contour maps, surface area and capacity tables, and bathymetric change maps for selected water-supply lakes in northwestern Missouri, 2019 and 2020

Bathymetric data were collected at 12 water-supply lakes in northwestern Missouri by the U.S. Geological Survey, in cooperation with the Missouri Department of Natural Resources and in collaboration with various local agencies, as part of a multiyear effort to establish or update the surface area and capacity tables for the surveyed lakes. Ten of the lakes were surveyed from July to September 2019, one of the original 10 was resurveyed in March 2020, and two lakes of high interest near Maryville were surveyed in June 2020. Six of the lakes had been surveyed by the U.S. Geological Survey before, and the recent surveys were compared to the earlier surveys to document the changes in the bathymetric surface and capacity of the lake and to produce a bathymetric change map. Bathymetric data were collected using a high-resolution multibeam mapping system mounted on a boat. Supplemental depth data were collected in shallow areas with an acoustic Doppler current profiler on a remote-controlled boat. At Hamilton Reservoir, a Global Navigation Satellite System survey receiver was used to collect additional bathymetric data at several points across four transects and around the perimeter of a substantial shallow area filled with aquatic vegetation upstream from a low-clearance bridge on the northern arm. Data points from the various sources were exported at a gridded data resolution appropriate to each lake. Data outside the multibeam echosounder survey extent and greater than the surveyed water-surface elevation generally were obtained from data collected using aerial light detection and ranging point cloud data, 1/9 arc-second National Elevation Dataset data based on aerial light detection and ranging data, or both. A linear enforcement technique was used to add points to the dataset in areas of sparse data (the upper ends of coves where the water was shallow or aquatic vegetation precluded data acquisition) based on surrounding multibeam and upland data values. The various point datasets were used to produce a three-dimensional triangulated irregular network surface of the lake-bottom elevations for each lake. A surface area and capacity table was produced from the three-dimensional surface showing surface area and capacity at specified lake water-surface elevations. Various quality-assurance tests were conducted to ensure quality data were collected with the multibeam, including beam angle checks and patch tests. Additional quality-assurance tests were conducted on the gridded bathymetric data from the survey, the bathymetric surface created from the gridded data, and the contours created from the bathymetric survey. If data from a previous bathymetric survey existed at a given lake, a bathymetric change map was generated from the elevation difference between the previous survey and the 2019 bathymetric survey data points. After applying any vertical elevation changes to the previous survey data to ensure a match to the 2019 survey datum, coincident points between the surveys were found, and a bathymetric change map was generated using the coincident point data. A decrease in capacity was observed at all the lakes for which a previous survey existed. The decrease in capacity at the primary spillway or intake elevation ranged from 0.8 percent at Lake Viking to 21.4 percent at Middle Fork Grand River Reservoir. The mean bathymetric change ranged from 0.33 foot at Willow Brook Lake to 1.18 feet at Middle Fork Grand River Reservoir. The computed sedimentation rate generally ranged from 0.54 to 4.19 acre-feet per year at Maysville Lake and Middle Fork Grand River Reservoir, respectively; however, Lake Viking had the largest sedimentation rate of 14.9 acre-feet per year, despite having the smallest decrease in capacity at the spillway elevation of only 0.8 percent and a mean bathymetric change of only 0.4 foot. Evidence of dredging was observed in the bathymetric surface for Lake Viking. Some changes observed in some bathymetric change maps are hypothesized to result from the difference in data collection equipment and techniques between the previous and present bathymetric surveys. Certain erosional features around the perimeter of certain lakes may be the result of wave action during low-water years.

Missouri↗

Sand resources, regional geology, and coastal processes of the Chandeleur Islands Coastal System: An evaluation of the Breton National Wildlife Refuge

Breton National Wildlife Refuge, the Chandeleur Islands chain in Louisiana, provides habitat and nesting areas for wildlife and is an initial barrier protecting New Orleans from storms. The U.S. Geological Survey (USGS) in partnership with the University of New Orleans Pontchartrain Institute for Environmental Sciences undertook an intensive study that included (1) an analysis of island change based on historical maps and remotely sensed shoreline and topographic data; (2) a series of lidar surveys at 3- to 4-month intervals after Hurricane Katrina to determine barrier island recovery potential; (3) a discussion of sea level rise and effects on the islands; (4) an analysis of sea floor evolution and sediment dynamics in the refuge over the past 150 years; (5) an assessment of the local sediment transport and sediment resource availability based on the bathymetric and subbottom data; (6) a carefully selected core collection effort to groundtruth the geophysical data and more fully characterize the sediments composing the islands and surrounds; (7) an additional survey of the St. Bernard Shoals to assess their potential as a sand resource; and (8) a modeling study to numerically simulate the potential response of the islands to the low-intensity, intermediate, and extreme events likely to affect the refuge over the next 50 years. Results indicate that the islands have become fragmented and greatly diminished in subaerial extent over time: the southern islands retreating landward as they reorganize into subaerial features, the northern islands remaining in place. Breton Island, because maintenance of the Mississippi River-Gulf Outlet (MRGO) outer bar channel requires dredging, is deprived of sand sufficient to sustain itself. Regional sediment transport trends indicate that large storms are extremely effective in transporting sand and controlling the shoreline development and barrier island geometry. Sand is transported north and south from a divergent zone near Monkey Bayou at the southern end of the Chandeleur Islands. Numerical simulation of waves and sediment transport supports the geophysical results and indicates that vast areas of the lower shoreface are affected and are undergoing erosion during storm events, that there is little or no fair weather mechanism to rework material into the littoral system, and that as a result, there is a net loss of sediment from the system. Lidar surveys revealed that the island chain immediately after Hurricane Katrina lost about 84 percent of its area and about 92 percent of its prestorm volume. Marsh platforms that supported the islands’ sand prior to the storm were reduced in width by more than one-half. Repeated lidar surveys document that in places the shoreline has retreated about 100 m under the relatively low-energy waves since Hurricanes Katrina and Rita; however, this retreat is nonuniform. Recent high-resolution geophysical surveys of the sea floor and subsurface within 5–6 km of the Chandeleur Islands during 2006 and 2007 show that, in addition to the sand that is rebuilding portions of the island chain, a large volume of sand is contained in Hewes Point, in an extensive subtidal spit platform that has formed at the northern end of the Chandeleur Islands. Hewes Point appears to be the depositional terminus of the alongshore transport system. In the southern Chandeleurs, sand is being deposited in a broad tabular deposit near Breton Island called the southern offshore sand sheet. These two depocenters account for approximately 70 percent of the estimated sediment volume located in potential borrow sites. An additional large potential source of sand for restoration lies in the St. Bernard Shoals, which are estimated to contain approximately 200 × 10 6 m 3 of sand. Successful restoration planning for the Breton National Wildlife Refuge should mimic the natural processes of early stages of barrier island evolution including lateral transport to the flanks of the island chain from a centralized sand source that will ultimately enhance the ability of the islands to naturally build backbarrier marsh, dunes, and a continuous sandy shoreline. Barrier island sediment nourishment should be executed with the understanding that gulf shoreline erosion is inevitable but that island area can be maintained and enhanced during retreat (thus significantly prolonging the life of the island chain) with strategic sand placement.

Louisiana↗

Bathymetric contour maps, surface area and capacity tables, and bathymetric change maps for selected water-supply lakes in north-central and west-central Missouri, 2020

Bathymetric data were collected at 10 water-supply lakes in north-central and west-central Missouri by the U.S. Geological Survey (USGS) in cooperation with the Missouri Department of Natural Resources and in collaboration with various local agencies, as part of a multiyear effort to establish or update the surface area and capacity tables for the surveyed lakes. The lakes were surveyed in June and July 2020. Seven of the lakes had been surveyed by the USGS between 2002 and 2007, and the recent surveys were compared to the earlier surveys to document changes in the bathymetric surface and capacity of the lake and produce a bathymetric change map. Bathymetric data were collected using a high-resolution multibeam mapping system mounted on a boat. Supplemental depth data at two of the lakes were collected in shallow areas with an acoustic Doppler current profiler on a remote-controlled boat. Data points from the various sources were exported at a gridded data resolution appropriate to each lake, either 0.82 foot or 1.64 feet. Data outside the multibeam survey extent and greater than the surveyed water-surface elevation generally were obtained from data collected using aerial light detection and ranging (lidar) point cloud data, except at Holden City Lake. A linear enforcement technique was used to add points to the dataset in areas of sparse data (the upper ends of coves where the water was shallow or aquatic vegetation precluded data acquisition) based on surrounding multibeam and upland data values. The various point datasets were used to produce a three-dimensional triangulated irregular network surface of lake-bottom elevations for each lake. A surface area and capacity table was produced from the three-dimensional surface for each lake showing surface area and capacity at specified lake water-surface elevations. Various quality-assurance tests were conducted to ensure quality data were collected with the multibeam, including beam angle checks and patch tests. Additional quality-assurance tests were conducted on the gridded bathymetric data from the survey, the bathymetric surface created from the gridded data, and the contours created from the bathymetric survey. If data from a previous bathymetric survey existed at a given lake, a bathymetric change map was generated from the elevation difference between the previous survey and the 2020 bathymetric survey data points. After reconciling any vertical datum disagreement between the previous survey and the 2020 survey, coincident points between the surveys were identified, and a bathymetric change map was generated using the coincident point data. A decrease in capacity was observed at nearly all the lakes for which a previous survey existed, and the mean bathymetric change between the surveys was positive at all the lakes. The decrease in capacity at the primary spillway elevation ranged from –0.4 percent at Edwin A Pape Lake to 9.4 percent at upper Higginsville Reservoir. The mean bathymetric change ranged from 0.03 foot at Garden City New Lake to 1.75 feet at Harrisonville City Lake, which corresponds to a time-averaged mean bathymetric change ranging from 0.002 foot per year at Garden City New Lake to 0.132 foot per year at Harrisonville City Lake. The computed volumetric sedimentation rate generally ranged from 0.04 to 4.91 acre-feet per year at Garden City New Lake and Holden City Lake, respectively; however, Harrisonville City Lake had a substantially larger volumetric sedimentation rate of 42.7 acre-feet per year, corresponding to the substantial mean bathymetric change of 1.75 feet and combined with the relatively shorter interval between surveys. Harrisonville City Lake also had the second-largest decrease in capacity at the spillway elevation of 5.9 percent. As with the 2019 surveys, some changes observed in the bathymetric change maps likely result from the difference in data collection equipment and techniques between the surveys. Certain apparent erosional features around the perimeter of certain lakes may be the result of wave action or compaction of sediments exposed to air during low-water years, or may indicate an unidentified but systemic error in the older singlebeam echosounder survey data.

Missouri↗

Sediment loads in the Red River of the North and selected tributaries near Fargo, North Dakota, 2010--2011

Natural-resource agencies are concerned about possible geomorphic effects of a proposed diversion project to reduce the flood risk in the Fargo-Moorhead metropolitan area. The U.S. Geological Survey in cooperation with the U.S. Army Corps of Engineers collected data in the spring of 2010 and 2011, and from June to November 2011, during rainfall-runoff events and base-flow conditions to provide information on sediment transport. The data were used to examine sediment concentrations, loads, and particle-size distributions at nine selected sites in the Red River and its tributaries near the Fargo-Moorhead metropolitan area. Suspended-sediment concentration varied among sites in 2010 and 2011. The least suspended-sediment concentrations were measured at the Red River (site 1) and the Buffalo River (site 9), and the greatest concentrations were measured at the two Sheyenne River sites (sites 3 and 4). Estimated daily suspended-sediment loads were highly variable in 2010 and 2011 in the Red River and its tributaries, with the greatest loads occurring in the spring and the smallest loads occurring in the winter. For the Red River, daily suspended-sediment loads ranged from 26 to 3,500 tons per day at site 1 and from 30 to 9,010 tons per day at site 2. For the Sheyenne River, daily loads ranged from less than 10 to 10,200 tons per day at site 3 and from less than 10 to 4,530 tons per day at site 4. The mean daily load was 191 tons per day in 2010 and 377 tons per day in 2011 for the Maple River, and 610 tons per day in 2011 for the Wild Rice River (annual loads were not computed for 2010). For the three sites that were only sampled in 2011 (sites 7, 8 and 9), the mean daily suspended-sediment loads ranged from 40 tons per day at the Lower Branch Rush River (site 8) to 118 tons per day at the Buffalo River (site 9). For sites that had estimated loads in 2010 and 2011 (sites 1–5), estimated annual (March–November) suspended-sediment loads were greater in 2011 compared to 2010. In 2010, annual loads ranged from 68,650 tons per year at the Maple River (site 5) to 249,040 tons per year at the Sheyenne River (site 3). In 2011, when all nine sites were sampled, annual loads ranged from 8,716 tons per year at the Lower Branch Rush River (site 8) to 552,832 tons per year at the Sheyenne River (site 3). With the exception of the Sheyenne River (site 4), the greatest monthly loads occurred in March for 2010, with as little as 27 percent (site 1) and as much as 42 percent (site 3) of the annual load occurring in March. For 2011, the greatest monthly loads occurred in April, ranging from 33 percent (site 1) to 63 percent (site 7) of the 2011 annual load. A relatively small amount of sediment was transported past the nine sites as bedload in 2010 and 2011. For most of the samples collected at the nine sites, the bedload composed less than 1 percent of the calculated daily total sediment load.

North Dakota↗

Effects of surface applications of biosolids on soil, crops, ground water, and streambed sediment near Deer Trail, Colorado, 1999-2003

The U.S. Geological Survey, in cooperation with Metro Wastewater Reclamation District and North Kiowa Bijou Groundwater Management District, studied natural geochemical effects and the effects of biosolids applications to the Metro Wastewater Reclamation District properties near Deer Trail, Colorado, during 1999 through 2003 because of public concern about potential contamination of soil, crops, ground water, and surface water from biosolids applications. Parameters analyzed for each monitoring component included arsenic, cadmium, copper, lead, mercury, molybdenum, nickel, selenium, and zinc (the nine trace elements regulated by Colorado for biosolids), gross alpha and gross beta radioactivity, and plutonium, as well as other parameters. Concentrations of the nine regulated trace elements in biosolids were relatively uniform and did not exceed applicable regulatory standards. All plutonium concentrations in biosolids were below the minimum detectable level and were near zero. The most soluble elements in biosolids were arsenic, molybdenum, nickel, phosphorus, and selenium. Elevated concentrations of bismuth, mercury, phosphorus, and silver would be the most likely inorganic biosolids signature to indicate that soil or streambed sediment has been affected by biosolids. Molybdenum and tungsten, and to a lesser degree antimony, cadmium, cobalt, copper, mercury, nickel, phosphorus, and selenium, would be the most likely inorganic 'biosolids signature' to indicate ground water or surface water has been affected by biosolids. Soil data indicate that biosolids have had no measurable effect on the concentration of the constituents monitored. Arsenic concentrations in soil of both Arapahoe and Elbert County monitoring sites (like soil from all parts of Colorado) exceed the Colorado soil remediation objectives and soil cleanup standards, which were determined by back-calculating a soil concentration equivalent to a one-in-a-million cumulative cancer risk. Lead concentrations in soil slightly exceed the U.S. Environmental Protection Agency toxicity-derived ecological soil-screening levels for avian wildlife. Plutonium concentration in the soil was near zero. Wheat-grain data were insufficient to determine any measurable effects from biosolids. Comparison with similar data from other parts of North America where biosolids were not applied indicates similar concentrations. However, the Deer Trail study area had higher nickel concentrations in wheat from both the biosolids-applied fields and the control fields. Plutonium content of the wheat was near zero. Ground-water levels generally declined at most wells during 1999 through 2003. Ground-water quality did not correlate with ground-water levels. Vertical ground-water gradients during 1999 through 2003 indicate that bedrock ground-water resources downgradient from the biosolids-applied areas are not likely to be contaminated by biosolids applications unless the gradients change as a result of pumping. Ground-water quality throughout the study area varied over time at each site and from site to site at the same time, but plutonium concentrations in the ground water always were near zero. Inorganic concentrations at well D6 were relatively high compared to other ground-water sites studied. Ground-water pH and concentrations of fluoride, nitrite, aluminum, arsenic, barium, chromium, cobalt, copper, lead, mercury, nickel, silver, zinc, and plutonium in the ground water of the study area met Colorado standards. Concentrations of chloride, sulfate, nitrate, boron, iron, manganese, and selenium exceeded Colorado ground-water standards at one or more wells. Nitrate concentrations at well D6 significantly (alpha = 0.05) exceeded the Colorado regulatory standard. Concentrations of arsenic, cadmium, chromium, lead, mercury, nickel, and zinc in ground water had no significant (alpha = 0.05) upward trends. During 1999-2003, concentrations of nitrate, copper, molybdenum, and selenium

Scientific Investigations Report↗

Geohydrology, geochemistry, and groundwater simulation (1992-2011) and analysis of potential water-supply management options, 2010-60, of the Langford Basin, California

Groundwater withdrawals began in 1992 from the Langford Basin within the Fort Irwin National Training Center (NTC), California. From April 1992 to December 2010, approximately 12,300 acre-feet of water (averaging about 650 acre-feet per year) has been withdrawn from the basin and transported to the adjacent Irwin Basin. Since withdrawals began, water levels in the basin have declined by as much as 40 feet, and the quality of the groundwater withdrawn from the basin has deteriorated. The U.S. Geological Survey collected geohydrologic data from Langford Basin during 1992–2011 to determine the quantity and quality of groundwater available in the basin. Geophysical surveys, including gravity, seismic refraction, and time-domain electromagnetic induction surveys, were conducted to determine the depth and shape of the basin, to delineate depths to the Quaternary-Tertiary interface, and to map the depth to the water table and changes in water quality. Data were collected from existing wells and test holes, as well as 11 monitor wells that were installed at 5 sites as part of this study. Water-quality samples collected from wells in the basin were used to determine the groundwater chemistry within the basin and to delineate potential sources of poor-quality groundwater. Analysis of stable isotopes of oxygen and hydrogen in groundwater indicates that present-day precipitation is not a major source of recharge to the basin. Tritium and carbon-14 data indicate that most of the basin was recharged prior to 1952, and the groundwater in the basin has an apparent age of 12,500 to 30,000 years. Recharge to the basin, estimated to be less than 50 acre-feet per year, has not been sufficient to replenish the water that is being withdrawn from the basin. A numerical groundwater-flow model was developed for the Langford Basin to better understand the aquifer system used by the Fort Irwin NTC as part of its water supply, and to provide a tool to help manage groundwater resources at the NTC. Measured groundwater-level declines since the initiation of withdrawals (1992–2011) were used to calibrate the groundwater-flow model. The simulated recharge was about 46 acre-feet per year, including approximately 6 acre-feet per year of natural recharge derived from precipitation runoff and as much as 40 acre-feet per year of underflow from the Irwin Basin. Between April 1992 and December 2010, an average of about 650 acre-feet per year of water was withdrawn from the Langford Basin. Groundwater withdrawals in excess of natural recharge resulted in a net loss of 11,670 acre-feet of groundwater storage within the basin for the simulation period. The Fort Irwin NTC is considering various groundwater-management options to address the limited water resources in the Langford Basin. The calibrated Langford Basin groundwater-flow model was used to evaluate the hydrologic effects of four groundwater-withdrawal scenarios being considered by the Fort Irwin NTC over the next 50 years (January 2011 through December 2060). Continuation of the 2010 withdrawal rate in the three existing production wells will result in 70 feet of additional drawdown in the central part of the basin. Redistributing the 2010 withdrawal rate equally to the three existing wells and two proposed new wells in the northern and southern parts of the basin would result in about 10 feet less drawdown in the central part of the basin but about 100 feet of additional drawdown in the new well in the northern part of the basin and about 50 feet of additional drawdown in the new well in the southern part of the basin. Reducing the withdrawals from the three existing production wells in the central part of the basin from about 45,000 acre-feet to about 32,720 acre-feet would result in about 40 feet of additional drawdown in the central basin near the pumping wells, about 25 feet less than if withdrawals were not reduced. The combination of reducing and redistributing the cumulative withdrawals to the three existing and two proposed new wells results in about 40 feet of additional drawdown in the central and southern parts of the basin and about 70 feet in the northern part of the basin. These results show that reducing and redistributing the groundwater withdrawals would maintain the upper aquifer at greater than 50 percent of its predevelopment saturated thickness throughout the groundwater basin. The scenarios simulated for this study demonstrate how the calibrated model can be utilized to evaluate the hydrologic effects of different water-management strategies.

California↗

Hydrogeologic conceptual model of groundwater occurrence and brine discharge to the Dolores River in the Paradox Valley, Montrose County, Colorado

Salinity, or total dissolved solids (TDS), of the Colorado River is a major concern in the southwestern United States where the river provides water to about 40 million people for municipal and industrial use and is used to irrigate about 5.5 million acres of land. Much of the salinity in the Colorado River Basin is derived from natural interactions of surface water and groundwater with various geologic materials (rocks, soils, and alluvial deposits). The Dolores River in southwest Colorado is a major tributary of the Colorado River that historically accounts for about 6 percent of the salinity load to the Upper Colorado River Basin with the Paradox Valley being the primary source of salinity to the Dolores River. The Paradox Valley, one of several salt-anticline valleys in the region, is a fault-bounded topographic basin aligned with and exposing an underlying salt-anticline core. Salt deposits in the Pennsylvanian Paradox Formation of the Hermosa Group form an elongated salt diapir oriented northwest to southeast that is up to 12,000 feet (ft) thick beneath the present valley floor. Surface erosion, groundwater circulation, and weathering during Tertiary and Quaternary valley formation contributed to development of a cap rock, collapse features, breccia, and brine at the top of the exposed salt diapir. Today (2023), brine occurring in the brecciated cap rock and underlying salt deposits is in hydraulic connection with an overlying freshwater alluvial aquifer, and depending on seasonal river stage and hydrologic conditions, the brine discharges to the Dolores River causing the observed increase in salinity as the river crosses the Paradox Valley. To reduce salinity concentrations in the Dolores River, the Bureau of Reclamation (Reclamation) operates the Paradox Valley Unit (PVU). The PVU project consists of nine shallow brine pumping wells near the Dolores River and one deep disposal well where the brine is injected for disposal. When operational, the PVU pumping wells extract brine from the base of the alluvial aquifer that is piped and injected into a deep disposal well about 3 miles southwest of the PVU. The PVU became fully operational July 1, 1996, and by 2015, operation of the PVU had reduced salinity concentrations in the Dolores River by as much as 70 percent compared to pre-PVU conditions. In response to a 4.5 magnitude earthquake, injection operations, and thus PVU pumping, were ceased from March 2019 to June 2022. A trial period of PVU operation began in June 2022 with a reduced injection rate, and thus PVU pumping rate, of about two-thirds capacity to gather additional information and guide future operational decisions. In cooperation with Reclamation, the U.S. Geological Survey (USGS) developed this report to present the current (2023) understanding of groundwater and brine occurrence and discharge to the Dolores River in the Paradox Valley. Results from the compilation of spatial datasets, groundwater sampling and age dating, and aquifer tests are presented to provide improved understanding of the Paradox Valley hydrogeology, to supply datasets for a numerical groundwater-flow and brine-transport model, and to support future operations of the PVU. The hydrogeologic data provided herein, along with the most recent loading analysis for the Dolores River in the Paradox Valley, and a previous conceptual model for brine discharge to the river are used to present a conceptual understanding of groundwater occurrence in the Paradox Valley.

Colorado↗