Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Nature Scientific Reports”

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,639 records · Page 91Linked to original sources
USGS2021

Hydraulic characterization of carbonate-rock and basin-fill aquifers near Long Canyon, Goshute Valley, northeastern Nevada

Understanding groundwater flow and pumping effects near pending mining operations requires accurate subsurface hydraulic characterization. To improve conceptual models of groundwater flow and development in the complex hydrogeologic system near Long Canyon Mine, in northwestern Goshute Valley, northeastern Nevada, the U.S. Geological Survey characterized the hydraulic properties of carbonate rocks and basin-fill aquifers using an integrated analysis of steady-state and stressed aquifer conditions informed by water chemistry and aquifer-test data. Hydraulic gradients and groundwater-age data in northern Goshute Valley indicate carbonate rocks in the Pequop Mountains just west and south of the Long Canyon Mine project area constitute a more permeable and active flow system than saturated rocks in the northern Pequop Mountains, western Toano Range, and basin fill. Permeable carbonate rocks in the northern Pequop Mountains, in part, discharge to the Johnson Springs wetland complex (JSWC), where mean groundwater ages range from 500 to 2,400 years and samples all contain a small fraction of modern waters, relative to mean ages of 8,600 to more than 22,000 years for most groundwater sampled to the north and east. Recharge to the JSWC occurs from a roughly 27-square-mile area in the upgradient Pequop Mountains to the west, composed mostly of permeable carbonate rock and fractured quartzite, and bounded by low-permeability shales and marbleized and siliclastic rocks. Single-well aquifer-test analyses provided transmissivity estimates at pumped wells. Transmissivity estimates ranged from 7,000 to 400,000 feet squared per day (ft 2 /d) in carbonate rocks and from 2,000 to 80,000 ft 2 /d in basin fill near the Long Canyon Mine. Water-level drawdown from multiple-well aquifer testing and rise from unintentional leakage into the overlying basin-fill aquifer were estimated and distinguished from natural fluctuations in 93 pumping and monitoring sites using analytical water-level models. Leakage of disposed aquifer-test pumpage occurred south of the aquifer test area through an unlined irrigation ditch. Drawdown was detected at distances of as much as 3 miles (mi) from pumping wells at all but one carbonate-rock site, at basin-fill sites on the alluvial fan immediately downgradient from pumping wells, and in Big Spring and spring NS-05. Similar drawdowns in carbonate rocks within the drawdown detection area suggest all wells penetrate a highly transmissive zone (HTZ) that is bounded by low-permeability rocks. Drawdown was not detected in carbonate rocks to the west of Canyon fault, in any basin-fill sites on the valley floor east of the Hardy fault, or at volcanic sites to the north, indicating that these major fault structures and (or) permeability contrasts between hydrogeologic units impeded groundwater flow or obscured pumping signals. Alternatively, unintentional leakage might have obscured drawdown at basin-fill sites on the valley floor, where water-level rise was detected at nine sites over 3 mi. Consistent hydraulic properties were estimated by simultaneously interpreting steady-state flow during predevelopment conditions and changes in groundwater levels and springflows from the 2016 carbonate-rock aquifer test with an integrated groundwater-flow model. Hydraulic properties were distributed across carbonate rocks, basin fill, volcanic rocks, and siliciclastic rocks with a hydrogeologic framework developed from geologic mapping and hydraulic testing. Estimated transmissivity distributions spanned at least three orders of magnitude in each rock unit. In the HTZ, simulated transmissivities ranged from 10,000 to 23,000,000 ft 2 /d, with the most transmissive areas occurring around Big Spring. Comparatively low carbonate-rock transmissivities of less than 10,000 ft 2 /d were estimated in the northern Pequop Mountains and poorly defined values of less than 1,000 ft 2 /d were estimated in the western Toano Range. Transmissivities in basin fill ranged from less than 10 to 80,000 ft 2 /d and were minimally constrained by the 2016 carbonate-rock aquifer test because poorly quantified leakage affected water levels more so than pumping. The most transmissive areas were informed by single-well aquifer tests along the eastern edge of the Pequop Mountains near Long Canyon Mine and could be indicative of a hydraulic connection between basin fill and more transmissive underlying carbonate rocks. Simulated transmissivities of volcanic and low-permeability rocks mostly are less than 1,000 ft 2 /d. The estimated hydraulic-property distributions and informed interpretation of hydraulic connections among hydrogeologic units improved the characterization and representation of groundwater flow near the Long Canyon Mine.

Nevada↗
USGS2013

Pesticides in Wyoming Groundwater, 2008-10

Groundwater samples were collected from 296 wells during 1995-2006 as part of a baseline study of pesticides in Wyoming groundwater. In 2009, a previous report summarized the results of the baseline sampling and the statistical evaluation of the occurrence of pesticides in relation to selected natural and anthropogenic (human-related) characteristics. During 2008-10, the U.S. Geological Survey, in cooperation with the Wyoming Department of Agriculture, resampled a subset (52) of the 296 wells sampled during 1995-2006 baseline study in order to compare detected compounds and respective concentrations between the two sampling periods and to evaluate the detections of new compounds. The 52 wells were distributed similarly to sites used in the 1995-2006 baseline study with respect to geographic area and land use within the geographic area of interest. Because of the use of different types of reporting levels and variability in reporting-level values during both the 1995-2006 baseline study and the 2008-10 resampling study, analytical results received from the laboratory were recensored. Two levels of recensoring were used to compare pesticides—a compound-specific assessment level (CSAL) that differed by compound and a common assessment level (CAL) of 0.07 microgram per liter. The recensoring techniques and values used for both studies, with the exception of the pesticide 2,4-D methyl ester, were the same. Twenty-eight different pesticides were detected in samples from the 52 wells during the 2008-10 resampling study. Pesticide concentrations were compared with several U.S. Environmental Protection Agency drinking-water standards or health advisories for finished (treated) water established under the Safe Drinking Water Act. All detected pesticides were measured at concentrations smaller than U.S. Environmental Protection Agency drinking-water standards or health advisories where applicable (many pesticides did not have standards or advisories). One or more pesticides were detected at concentrations greater than the CAL in water from 16 of 52 wells sampled (about 31 percent) during the resampling study. Detected pesticides were classified into one of six types: herbicides, herbicide degradates, insecticides, insecticide degradates, fungicides, or fungicide degradates. At least 95 percent of detected pesticides were classified as herbicides or herbicide degradates. The number of different pesticides detected in samples from the 52 wells was similar between the 1995-2006 baseline study (30 different pesticides) and 2008-2010 resampling study (28 different pesticides). Thirteen pesticides were detected during both studies. The change in the number of pesticides detected (without regard to which pesticide was detected) in groundwater samples from each of the 52 wells was evaluated and the number of pesticides detected in groundwater did not change for most of the wells (32). Of those that did have a difference between the two studies, 17 wells had more pesticide detections in groundwater during the 1995-2006 baseline study, whereas only 3 wells had more detections during the 2008-2010 resampling study. The difference in pesticide concentrations in groundwater samples from each of the 52 wells was determined. Few changes in concentration between the 1995-2006 baseline study and the 2008-2010 resampling study were seen for most detected pesticides. Seven pesticides had a greater concentration detected in the groundwater from the same well during the baseline sampling compared to the resampling study. Concentrations of prometon, which was detected in 17 wells, were greater in the baseline study sample compared to the resampling study sample from the same well 100 percent of the time. The change in the number of pesticides detected (without regard to which pesticide was detected) in groundwater samples from each of the 52 wells with respect to land use and geographic area was calculated. All wells with land use classified as agricultural had the same or a smaller number of pesticides detected in the resampling study compared to the baseline study. All wells in the Bighorn Basin geographic area also had the same or a smaller number of pesticides detected in the resampling study compared to the baseline study.

Wyoming↗
USGS2007

Introduction to Regional Geology, Tectonics, and Metallogenesis of Northeast Asia

Purpose This introduction presents an overview of the regional geology, tectonics, and metallogenesis of Northeast Asia. The major purposes are to provide a relatively short summary of these features for readers who are unfamiliar with Northeast Asia; a general scientific introduction for the succeeding chapters of this volume; and an overview of the methodology of metallogenic and tectonic analysis employed for Northeast Asia. The introduction also describes how a high-quality metallogenic and tectonic analysis, including synthesis of an associated metallogenic-tectonic model will greatly benefit refinement of mineral deposit models and deposit genesis; improvement of assessments of undiscovered mineral resources as part of quantitative mineral resource assessment studies; land-use and mineral exploration planning; improvement of interpretations of the origins of host rocks, mineral deposits, and metallogenic belts; and suggestions for new research. The compilation, synthesis, description, and interpretation of metallogenesis and tectonics of major regions, such as Northeast Asia (Eastern Russia, Mongolia, northern China, South Korea, and Japan) and the Circum-North Pacific (Russian Far East, Alaska, and Canadian Cordillera) requires a complex methodology. The methodology includes: (1) definitions of key terms; (2) compilation of a regional geologic base map that can be interpreted according to modern tectonic concepts and definitions; (3) compilation of a mineral deposit database that enables the determination of mineral deposit models, and relations of deposits to host rocks and tectonic origins; (4) synthesis of a series of mineral deposit models that characterize the known mineral deposits and inferred undiscovered deposits of the region; (5) compilation of a series of maps of metallogenic belts constructed on the regional geologic base map; and (6) formulation of a unified metallogenic and tectonic model. The summary of regional geology and metallogenesis in this introduction is based on publications of the major international collaborative studies of the metallogenesis and tectonics of Northeast Asia that were led by the U.S. Geological Survey. These studies have produced two broad types of publications. One type is a series of regional geologic, mineral deposit, and metallogenic belt maps and companion descriptions for the regions. Examples of major publications of this type are Obolenskiy and others (2003a, b, 2004), Parfenov and others (2003, 2004a, b), Nokleberg and others (2004), Rodionov and others (2004), and Naumova and others (2006). The other type is a suite of metallogenic and tectonic analyses of these same regions. Examples of major publications of this type are Rodionov and others (2004), Nokleberg and others (2000, 2004, 2005), and Naumova and others (2006). The Northeast Asia project area consists of eastern Russia (most of Siberia and most of the Russian Far East), Mongolia, Northern China, South Korea, Japan, and adjacent offshore areas. This area is approximately bounded by 30 to 82? N. latitude and 75 to 144? E. longitude. The major participating agencies are the Russian Academy of Sciences, Academy of Sciences of the Sakha Republic (Yakutia), VNIIOkeangeologia and Ministry of Natural Resources of the Russian Federation, Mongolian Academy of Sciences, Mongolian University of Science and Technology, Mongolian National University, Jilin University, Changchun, China, the China Geological Survey, the Korea Institute of Geosciences and Mineral Resources, the Geological Survey of Japan/AIST, University of Texas Arlington, and the U.S. Geological Survey. The Northeast Asia project extends and build on data and interpretations from a previous project on the Major Mineral Deposits, Metallogenesis, and Tectonics of the Russian Far East, Alaska, and the Canadian Cordillera that was conducted by the USGS, the Russian Academy of Sciences, the Alaska Division of Geological and Geophysical Surveys,

Open-File Report↗
USGS2005

Hydrogeology and ground-water/surface water interactions in the Des Moines River valley, southwestern Minnesota, 1997-2001

Increased water demand in and around Windom led the U.S. Geological Survey, in cooperation with the Minnesota Department of Natural Resources, local water suppliers, and Cottonwood County, to study the hydrology of aquifers in the Des Moines River Valley near Windom. The study area is the watershed of a 30-kilometer (19-mile) reach of the Des Moines River upstream from Windom. Based on stratigraphic analysis, two hydrologically and genetically separate surficial aquifers underlie the study area. The Windom aquifer has a saturated thickness of 34 meters (111 feet), and the Des Moines aquifer has a saturated thickness of 33 meters (108 ft). The surficial aquifers are relatively isolated from deeper aquifers by till, but some leakage probably occurs. Recharge to the aquifers is from areal recharge, from Cottonwood Lake, and from edge recharge. Pumping at the Windom well field induces substantial amounts of Cottonwood Lake water into the aquifer. During this study, the water level in a well located between two Red Rock wells and the river was lower than the river level during two periods. During those periods, water in the Des Moines River had the potential to recharge the aquifer. Discharge from the aquifers is primarily to municipal wells, the Des Moines River, and other surface waters. Most of the ground-water samples collected in the study area consisted of calcium-magnesium bicarbonate waters. Corn and soybean herbicides and their degradates were detected at low concentrations in 14 of 27 ground-water samples and in all 3 river samples. Metolachlor ethane sulfonic acid was the most commonly detected compound and also was detected at the highest concentrations. Nutrient concentrations in ground-water samples were skewed low with high outliers, and nutrient concentrations in river samples generally were less than analytical reporting limits. Nearly all recharge to the aquifer in the ground-water simulation was from edge recharge (80 percent). Calibrated net areal recharge ranged from 17 to 30 percent of the average annual precipitation. Isotopic composition of ground water and Cottonwood Lake water indicated about one-half of the water withdrawn from the Windom aquifer is from Cottonwood Lake. Scenarios tested with the calibrated model involved increased ground-water withdrawals and changes in recharge to simulate drier or wetter weather conditions. Doubling the withdrawals from all wells in the model had a small effect except in the Windom well-field area. Maximum head declines in the Red Rock well field and the Jeffers city well were less than 40 centimeters (15 inches). In the Windom well field, the maximum head decline was 11 meters (36 feet). The Windom well field does not induce recharge from the Des Moines River. The addition of a new well that pumped 2,000 cubic meters per day (0.44 million gallons per day) in the Augusta Lake Valley area caused a 0.83-meter-deep (2.72-foot-deep) cone of depression that extended to the valley walls. The drought scenario and the high-precipitation scenario resulted in head changes in the northern part of the Augusta Lake Valley area, in the southwestern part of the Red Rock area, and near the valley edges. Long-term withdrawals of water for public supplies may cause a net decrease in ground-water discharge to surface water. Water that does not evaporate, or that is not exported, is discharged to the Des Moines River but with changed water quality. Because ground-water and surface-water qualities in the study area are similar, the ground-water discharge probably has little effect on river water quality.

Minnesota↗
USGS2012

Sediment transport to and from small impoundments in northeast Kansas, March 2009 through September 2011

The U.S. Geological Survey, in cooperation with the Kansas Water Office, investigated sediment transport to and from three small impoundments (average surface area of 0.1 to 0.8 square miles) in northeast Kansas during March 2009 through September 2011. Streamgages and continuous turbidity sensors were operated upstream and downstream from Atchison County, Banner Creek, and Centralia Lakes to study the effect of varied watershed characteristics and agricultural practices on sediment transport in small watersheds in northeast Kansas. Atchison County Lake is located in a predominantly agricultural basin of row crops, with wide riparian buffers along streams, a substantial amount of tile drainage, and numerous small impoundments (less than 0.05 square miles; hereafter referred to as “ponds”). Banner Creek Lake is a predominantly grassland basin with numerous small ponds located in the watershed, and wide riparian buffers along streams. Centralia Lake is a predominantly agricultural basin of row crops with few ponds, few riparian buffers along streams, and minimal tile drainage. Upstream from Atchison County, Banner Creek, and Centralia Lakes 24, 38, and 32 percent, respectively, of the total load was transported during less than 0.1 percent (approximately 0.9 days) of the time. Despite less streamflow in 2011, larger sediment loads during that year indicate that not all storm events transport the same amount of sediment; larger, extreme storms during the spring may transport much larger sediment loads in small Kansas watersheds. Annual sediment yields were 360, 400, and 970 tons per square mile per year at Atchison County, Banner, and Centralia Lake watersheds, respectively, which were less than estimated yields for this area of Kansas (between 2,000 and 5,000 tons per square mile per year). Although Centralia and Atchison County Lakes had similar percentages of agricultural land use, mean annual sediment yields upstream from Centralia Lake were about 2.7 times those at Atchison County or Banner Creek Lakes. These data indicate larger yields of sediment from watersheds with row crops and those with fewer small ponds, and smaller yields in watersheds which are primarily grassland, or agricultural with substantial tile drainage and riparian buffers along streams. These results also indicated that a cultivated watershed can produce yields similar to those observed under the assumed reference (or natural) condition. Selected small ponds were studied in the Atchison County Lake watershed to characterize the role of small ponds in sediment trapping. Studied ponds trapped about 8 percent of the sediment upstream from the sediment-sampling site. When these results were extrapolated to the other ponds in the watershed, differences in the extent of these ponds was not the primary factor affecting differences in yields among the three watersheds. However, the selected small ponds were both 45 years old at the time of this study, and have reduced capacity because of being filled in with sediments. Additionally, trapping efficiency of these small ponds decreased over five observed storms, indicating that processes that suspended or resuspended sediments in these shallow ponds, such as wind and waves, affected their trapping efficiencies. While small ponds trapped sediments in small storms, they could be a source of sediment in larger or more closely spaced storm events. Channel slope was similar at all three watersheds, 0.40, 0.46, and 0.31 percent at Atchison County, Banner Creek, and Centralia Lake watersheds, respectively. Other factors, such as increased bank and stream erosion, differences in tile drainage, extent of grassland, or riparian buffers, could be the predominant factors affecting sediment yields from these basins. These results show that reference-like sediment yields may be observed in heavily agricultural watersheds through a combination of field-scale management activities and stream channel protection. When computing loads using published erosion rates obtained by single-point survey methodology, streambank contributions from the main stem of Banner Creek are three times more than the sediment load observed by this study at the sediment sampling site at Banner Creek, 2.6 times more than the sediment load observed by this study at the sediment sampling site at Clear Creek (upstream from Atchison County Lake), and are 22 percent of the load observed by this study at the sediment sampling site at Black Vermillion River above Centralia Lake. Comparisons of study sites to similarly sized urban and urbanizing watersheds in Johnson County, Kansas indicated that sediment yields from the Centralia Lake watershed were similar to those in construction-affected watersheds, while much smaller sediment yields in the Atchison County and Banner Creek watersheds were comparable to stable, heavily urbanized watersheds. Comparisons of study sites to larger watersheds upstream from Tuttle Creek Lake indicate the Black Vermillion River watershed continues to have high sediment yields despite 98 percent of sediment from the Centralia watershed (a headwater of the Black Vermillion River) being trapped in Centralia Lake. Estimated trapping efficiencies for the larger watershed lakes indicated that Banner Creek and Centralia Lakes trapped 98 percent of incoming sediment, whereas Atchison County Lake trapped 72 percent of incoming sediment during the 3-year study period.

Kansas↗
USGS2011

Effects of brush management on the hydrologic budget and water quality in and adjacent to Honey Creek State Natural Area, Comal County, Texas, 2001-10

The U.S. Geological Survey, in cooperation with the U.S. Department of Agriculture Natural Resources Conservation Service, the Edwards Region Grazing Lands Conservation Initiative, the Texas State Soil and Water Conservation Board, the San Antonio River Authority, the Edwards Aquifer Authority, Texas Parks and Wildlife, the Guadalupe Blanco River Authority, and the San Antonio Water System, evaluated the hydrologic effects of ashe juniper ( Juniperus ashei ) removal as a brush management conservation practice in and adjacent to the Honey Creek State Natural Area in Comal County, Tex. By removing the ashe juniper and allowing native grasses to reestablish in the area as a brush management conservation practice, the hydrology in the watershed might change. Using a simplified mass balance approach of the hydrologic cycle, the incoming rainfall was distributed to surface water runoff, evapotranspiration, or groundwater recharge. After hydrologic data were collected in adjacent watersheds for 3 years, brush management occurred on the treatment watershed while the reference watershed was left in its original condition. Hydrologic data were collected for another 6 years. Hydrologic data include rainfall, streamflow, evapotranspiration, and water quality. Groundwater recharge was not directly measured but potential groundwater recharge was calculated using a simplified mass balance approach. The resulting hydrologic datasets were examined for differences between the watersheds and between pre- and post-treatment periods to assess the effects of brush management. The streamflow to rainfall relation (expressed as event unit runoff to event rainfall relation) did not change between the watersheds during pre- and post-treatment periods. The daily evapotranspiration rates at the reference watershed and treatment watershed sites exhibited a seasonal cycle during the pre- and post-treatment periods, with intra- and interannual variability. Statistical analyses indicate the mean difference in daily evapotranspiration rates between the two watershed sites is greater during the post-treatment than the pre-treatment period. Average annual rainfall, streamflow, evapotranspiration, and potential groundwater-recharge conditions were incorporated into a single hydrologic budget (expressed as a percentage of the average annual rainfall) applied to each watershed before and after treatment to evaluate the effects of brush management. During the post-treatment period, the percent average annual unit runoff in the reference watershed was similar to that in the treatment watershed, however, the difference in percentages of average annual evapotranspiration and potential groundwater recharge were more appreciable between the reference and treatment watersheds than during the pre-treatment period. Using graphical comparisons, no notable differences in major ion or nutrient concentrations were found between samples collected at the reference watershed (site 1C) and treatment watershed (site 2C) during pre- and post-treatment periods. Suspended-sediment loads were calculated from samples collected at sites 1C and 2T. The relation between suspended-sediment loads and streamflow calculated from samples collected from sites 1C and 2T did not exhibit a statistically significant difference during the pre-treatment period, whereas during the post-treatment period, relation between suspended-sediment loads and streamflow did exhibit a statistically significant difference. The suspended-sediment load to streamflow relations indicate that for the same streamflow, the suspended-sediment loads calculated from site 2T were generally less than suspended-sediment loads calculated from site 1C during the post-treatment period.

Texas↗
USGS2013

Sources and characteristics of organic matter in the Clackamas River, Oregon, related to the formation of disinfection by-products in treated drinking water

This study characterized the amount and quality of organic matter in the Clackamas River, Oregon, to gain an understanding of sources that contribute to the formation of chlorinated and brominated disinfection by-products (DBPs), focusing on regulated DBPs in treated drinking water from two direct-filtration treatment plants that together serve approximately 100,000 customers. The central hypothesis guiding this study was that natural organic matter leaching out of the forested watershed, in-stream growth of benthic algae, and phytoplankton blooms in the reservoirs contribute different and varying proportions of organic carbon to the river. Differences in the amount and composition of carbon derived from each source affects the types and concentrations of DBP precursors entering the treatment plants and, as a result, yield varying DBP concentrations and species in finished water. The two classes of DBPs analyzed in this study-trihalomethanes (THMs) and haloacetic acids (HAAs)-form from precursors within the dissolved and particulate pools of organic matter present in source water. The five principal objectives of the study were to (1) describe the seasonal quantity and character of organic matter in the Clackamas River; (2) relate the amount and composition of organic matter to the formation of DBPs; (3) evaluate sources of DBP precursors in the watershed; (4) assess the use of optical measurements, including in-situ fluorescence, for estimating dissolved organic carbon (DOC) concentrations and DBP formation; and (5) assess the removal of DBP precursors during treatment by conducting treatability "jar-test" experiments at one of the treatment plants. Data collection consisted of (1) monthly sampling of source and finished water at two drinking-water treatment plants; (2) event-based sampling in the mainstem, tributaries, and North Fork Reservoir; and (3) in-situ continuous monitoring of fluorescent dissolved organic matter (FDOM), turbidity, chlorophyll- a , and other constituents to continuously track source-water conditions in near real-time. Treatability tests were conducted during the four event-based surveys to determine the effectiveness of coagulant and powdered activated carbon (PAC) on the removal of DBP precursors. Sample analyses included DOC, total particulate carbon (TPC), total and dissolved nutrients, absorbance and fluorescence spectroscopy, and, for regulated DBPs, concentrations of THMs and HAAs in finished water and laboratory-based THM and HAA formation potentials (THMFP and HAAFP, respectively) for source water and selected locations throughout the watershed. The results of this study may not be typical given the record and near record amounts of precipitation that occurred during spring that produced streamflow much higher than average in 2010-11. Although there were algal blooms, lower concentrations of chlorophyll- a were observed in the water column during the study period compared to historical data. Concentrations of DBPs in finished (treated) water averaged 0.024 milligrams per liter (mg/L) for THMs and 0.022 mg/L for HAAs; maximum values were about 0.040 mg/L for both classes of DBPs. Although DBP concentrations were somewhat higher within the distribution system, none of the samples collected for this study or for the quarterly compliance monitoring by the water utilities exceeded levels permissible under existing U.S. Environmental Protection Agency (USEPA) regulations: 0.080 mg/L for THMs and 0.060 mg/L for HAAs. DOC concentrations were generally low in the Clackamas River, typically about 1.0-1.5 mg/L. Concentrations in the mainstem occasionally increased to nearly 2.5 mg/L during storms; DOC concentrations in tributaries were sometimes much higher (up to 7.8 mg/L). The continuous in-situ FDOM measurements indicated sharp rises in DOC concentrations in the mainstem following rainfall events; concentrations were relatively stable during summer base flow. Even though the first autumn storm mobilized appreciable quantities of carbon, higher concentrations of DBPs in finished water were observed 3-weeks later, after the ground was saturated from additional rainfall. The majority of the DOC in the lower Clackamas River appears to originate from the upper basin, suggesting terrestrial carbon was commonly the dominant source. Lower-basin tributaries typically contained the highest concentrations of DOC and DBP precursors and contributed substantially to the overall loads in the mainstem during storms. During low-flow periods, tributaries were not major sources of DOC or DBP precursors to the Clackamas River. Although the dissolved fraction of organic carbon contributed the majority of DBP precursors, at times the particulate fraction (inorganic sediment and organic particles including detritus and algal material) contributed a substantial fraction of DBP precursors. Considering just the main-stem sites, on average, 10 percent of THMFP and 32 percent of HAAFP were attributed to particulate carbon. This finding suggests water-treatment methods that remove particles prior to chlorination would reduce finished-water DBP concentrations to some degree. Overall, concentrations of THM and HAA precursors were closely linked to DOC concentrations; laboratory DBP formation potentials (DBPFPs) clearly showed that THMFP and HAAFP were greatest in the downstream tributaries that contained elevated carbon concentrations. However, carbon-normalized "specific" formation potentials for THMs and HAAs (STHMFP and SHAAFP, respectively) revealed changes in carbon character over time that affected the two types of DBP classes differently. HAA precursors were elevated in waters containing aromatic-rich soil-derived material arising from forested areas. In contrast, THM precursors were associated with carbon having a lower aromatic content; highest STHMFP occurred in autumn 2011 in the mainstem from North Fork Reservoir downstream to LO DWTP. This pattern suggests the potential for a link between THM precursors and algal-derived carbon. The highest STHMFP value was measured within North Fork Reservoir, indicating reservoir derived carbon may be important for this class of DBPs. Weak correlations between STHMFP and SHAAFP emphasize that precursor sources for these types of DBPs may be different. This highlights not only that different locations within the watershed produce carbon with different reactivity (specific DBPFP), but also that different management approaches for each class of DBP precursors could be required for control. Treatability tests conducted on source water during four basin-wide surveys demonstrated that an average of about 40 percent of DOC can be removed by coagulation. While the decrease in THMFP following coagulation was similar to DOC, the decrease in HAAFP was much greater (approximately 70 percent), indicating coagulation is particularly effective at removing HAA precursors'likely because of the aromatic nature of the carbon associated with HAA precursors. Several findings from this study have direct implications for managing drinking-water resources and for providing useful information that may help improve treatment-plant operations. For example, the use of in-situ fluorometers that measure FDOM provided an excellent proxy for DOC concentration in this system and revealed short-term, rapid changes in DOC concentration during storm events. In addition, the strong correlation between FDOM values measured in-situ and HAA5 concentrations in finished water may permit estimation of continuous HAA concentrations, as was done here. As part of this study, multiple in-situ FDOM sensors were deployed continuously and in real-time to characterize the composition of dissolved organic matter. Although the initial results were promising, additional research and engineering developments will be needed to demonstrate the full utility of these sensors for this purpose. In conclusion, although DBPFPs were strongly correlated to DOC concentration, some DBPs formed from particulate carbon, including terrestrial leaf material and algal material such as planktonic species of blue-green algae and sloughed filaments, stalks, and cells of benthic algae. Different precursor sources in the watershed were evident from the data, suggesting specific actions may be available to address some of these sources. In-situ measurements of FDOM proved to be an excellent proxy for DOC concentration as well as HAA formation during treatment, which suggests further development and refinement of these sensors have the potential to provide real-time information about complex watershed processes to operators at the drinking-water treatment plants. Follow-up studies could examine the relative roles that terrestrial and algal sources have on the DBP precursor pool to better understand how watershed-management activities may be affecting the transport of these compounds to Clackamas River drinking-water intakes. Given the low concentrations of algae in the water column during this study, additional surveys during more typical river conditions could provide a more complete understanding of how algae contribute DBP precursors. Further development of FDOM-sensor technology can improve our understanding of carbon dynamics in the river and how concentrations may be trending over time. This study was conducted in collaboration with Clackamas River Water and the City of Lake Oswego water utilities. Other research partners included Oregon Health and Science University in Hillsboro, Oregon, Alexin Laboratory in Tigard, Oregon, U.S. Geological Survey National Research Program Laboratory in Denver, Colorado, and the U.S. Geological Survey Water Science Centers in Portland, Oregon, and Sacramento, California. This project was supported with funding from Clackamas River Water, City of Lake Oswego, the U.S. Geological Survey, and the Water Research Foundation.

Oregon↗
USGS2020

Effects of legacy sediment removal and effects on nutrients and sediment in Big Spring Run, Lancaster County, Pennsylvania, 2009–15

Big Spring Run is a 1.68-square mile watershed underlain by mostly carbonate rock in a mixed land-use setting (part agricultural and part developed) in Lancaster County, Pennsylvania. Big Spring Run is a subwatershed of Mill Creek, a tributary to the Conestoga River. These watersheds are known contributors of nutrient and sediment loads to the Chesapeake Bay and several stream reaches are on the Pennsylvania impaired waters list. Big Spring Run is listed as impaired and was selected by the Pennsylvania Department of Environmental Protection to evaluate a novel best management practice to restore natural aquatic ecosystems by removing legacy sediment. The study was designed to quantify sediment and nutrient contributions in pre- and postrestoration periods (water years 2009–11 and 2012–15, respectively) using an intensive monitoring approach at three surface-water sites within the watershed. Instrumentation at each site continuously measured (15-minute intervals) streamflow, water temperature, and turbidity. Water-quality samples were collected routinely (generally monthly and during selected storms); sampling frequency varied by site and constituent at the three monitoring sites. Effects of legacy sediment removal and restoration on nutrient concentrations varied in surface water samples depending on the form (particulate, dissolved, organic, inorganic). For example, total phosphorus concentrations at the downstream site decreased from a median of 0.19 milligram per liter (mg/L) to 0.04 mg/L, pre- and postrestoration periods, respectively. Concentrations of orthophosphate, the dissolved form of phosphorus, were not significantly different pre- to postrestoration at the downstream site. Similarly, nitrate concentrations, the dominant form of nitrogen in Big Spring Run surface-water samples (92.3 percent of total nitrogen) were not significantly different in the pre- compared to the postrestoration periods. Legacy sediment removal and restoration had significant effects on suspended-sediment concentrations and loads. Median suspended-sediment concentrations at the downstream site decreased from 556 mg/L prerestoration to 74 mg/L postrestoration even though streamflow hydrographs during the two periods were similar. In the postrestoration period, the mean annual suspended-sediment load conveyed to the restoration area from the upstream sites was 839 tons, whereas mean annual suspended-sediment load at the downstream site was reduced to 242 tons. Streamflow during storms transports a large proportion of the suspended-sediment load; there were a total of 320 storms over the study period. In Big Spring Run, a single storm event can transport more than 25 percent of the annual suspended-sediment load. The greatest single-storm contribution to suspended-sediment load was 38 percent in water year 2015 at the downstream site. Although streamflow magnitudes during storms varied greatly over the study period, median streamflow was 17.5 cubic feet per second and median duration was about 3 hours and 24 minutes. Results observed for this study using the newly proposed best management practice were compared with other best management practices intended to reduce sediment. For example, during a previous study, statistically significant reductions in suspended-sediment concentration were observed when streambank fencing was implemented in an adjacent watershed; however, suspended-sediment reductions were an order of magnitude less than the reductions observed in the current study. Median suspended-sediment concentration at the downstream site was reduced by 482 mg/L in the current study compared to only 30 to 46 mg/L as a result of streambank fencing.

Pennsylvania↗
USGS2024

Peak streamflow trends in South Dakota and their relation to changes in climate, water years 1921–2020

Peak-flow (flood) frequency analysis is essential to water-resources management applications, including the design of critical infrastructure such as bridges and culverts, and floodplain mapping. Federal guidelines for performing peak-flow flood frequency analyses are presented in a U.S. Geological Survey Techniques and Methods Report known as Bulletin 17C. A basic assumption within Bulletin 17C, which documents the guidelines for determining annual peak streamflow frequency, is that, for basins without major hydrologic alterations (for example, regulation, diversion, and urbanization), statistical properties of the distribution of annual peak streamflows are stationary; that is, the mean, variance, and skew are constant through time. Nonstationarity is a statistical property of a peak-flow series such that the long-term (on the order of decades) distributional properties change one or more times either gradually or abruptly through time. Individual nonstationarities may be attributed to one source such as flow regulation, land-use change, or climate but are often the result of a combination of sources, making detection and attribution of nonstationarities challenging. In response to a growing concern regarding nonstationarity in peak streamflows in the region, the U.S. Geological Survey, in cooperation with the Departments of Transportation of Illinois, Iowa, Michigan, Minnesota, Missouri, South Dakota, and Wisconsin; the Montana Department of Natural Resources and Conservation; and the North Dakota Department of Water Resources, assessed the potential nonstationarity in peak streamflows in the north-central United States. This chapter characterizes the effects of natural hydroclimatic shifts and potential climate change on annual peak streamflows in the State of South Dakota. Annual peak and daily streamflow as well as model-simulated gridded climatic data were examined for temporal monotonic trends, change points, and other statistical properties indicative of changing climatic and environmental conditions. Changes in annual peak and daily flows were evaluated among 13, 35, and 81 qualifying U.S. Geological Survey streamgages for the 75-, 50-, and 30-year trend periods through water year 2020 (the period from October 1, 2019, to September 30, 2020) in South Dakota, respectively. No qualifying streamgages were in the 100-year trend period in the State. Statistical tests for autocorrelation (independent and identically distributed assumption), monotonic trends, and change points in the median and scale are analyzed to evaluate potential stationarity violations (nonstationarity) for performing at-site peak-flow flood-frequency analysis. The trends are reported using a likelihood approach as an alternative to simply reporting significant trends with an arbitrary p -value cutoff point. A distinct east-west spatial pattern of likely upward and downward monotonic trends and change points, respectively, was detected in 75- and 50-year trend periods, but an inconsistent spatial pattern was detected in the 30-year trend period. Additionally, change points in the median annual peak streamflows were detected in the late 1970s and early 1980s in the western part of the State, but in the east, the change point was more commonly detected in 1992–93. A similar east-west spatial pattern of likely upward and downward trends was detected in the annual peak-flow timing, the day of the year of the annal peak streamflow. In the western part of the State, the annual peak streamflows are arriving earlier, but in the east, the annual peak streamflows are arriving later. A peaks-over-threshold (POT) analysis where, on average, there are two events per year (POT2) and four events per year (POT4) was also used to evaluate changes in the frequency (count) of daily streamflows exceeding the threshold. Similar to detected changes in the annual peak streamflow, an east-west likely upward or downward change corresponding to an increase or decrease, respectively, in the frequency of daily streamflow greater than a POT2 and POT4 threshold was detected. A monthly water-balance model was used to evaluate hydroclimatic variation in annual and seasonal precipitation, snowfall, potential evapotranspiration, and soil moisture storage for all qualifying streamgages in the 75-, 50-, and 30-year trend periods. Detected trends in the annual hydroclimatic metrics for the 75- and 50-year trend periods indicate a spatially consistent statewide increase in precipitation, decrease in snowfall, increase in potential evapotranspiration, and increase in soil moisture storage. Furthermore, detected trends in seasonal precipitation in the 75- and 50-year trend periods highlight a pronounced change in precipitation in winter and later into the summer season, especially in the 50-year trend period in the eastern part of the State. Statewide increases in seasonal soil moisture storage were also detected, highlighting year-round increasing flood magnitudes, particularly in the eastern part of the State. Based on the results of these stationarity tests for the qualifying streamgages in South Dakota among the 75-, 50-, and 30-year trend periods, consistent temporal and spatial patterns of nonstationarity were detected among the 75- and 50-year trend periods. Furthermore, when nonstationarity is detected in daily streamflow, increased streamflow and volume (increasing frequency in POT), as well as potentially bridge scour, may have implications on culvert and highway design in the eastern part of South Dakota. Thus, when performing at-site peak-flow flood-frequency analyses in South Dakota, potential nonstationarities and alternative approaches are important considerations.

South Dakota↗
USGS2015

Estimating natural recharge in San Gorgonio Pass watersheds, California, 1913–2012

A daily precipitation-runoff model was developed to estimate spatially and temporally distributed recharge for groundwater basins in the San Gorgonio Pass area, southern California. The recharge estimates are needed to define transient boundary conditions for a groundwater-flow model being developed to evaluate the effects of pumping and climate on the long-term availability of groundwater. The area defined for estimating recharge is referred to as the San Gorgonio Pass watershed model (SGPWM) and includes three watersheds: San Timoteo Creek, Potrero Creek, and San Gorgonio River. The SGPWM was developed by using the U.S. Geological Survey INFILtration version 3.0 (INFILv3) model code used in previous studies of recharge in the southern California region, including the San Gorgonio Pass area. The SGPWM uses a 150-meter gridded discretization of the area of interest in order to account for spatial variability in climate and watershed characteristics. The high degree of spatial variability in climate and watershed characteristics in the San Gorgonio Pass area is caused, in part, by the high relief and rugged topography of the area. Daily climate data developed from a network of monitoring sites and published average monthly precipitation maps were used to develop the climate inputs for the SGPWM. Geographic Information System (GIS) data defining land surface altitude, vegetation, soils, surficial geology, and land cover were used to define input parameters representing the physical characteristics of the land surface, root zone, and shallow subsurface underlying the root zone. Model parameterization was based on a previous INFILv3 model developed for an area including the upper parts of the San Timoteo Creek and Potrero Creek drainages and the western part of the San Gorgonio River watershed. The previous INFILv3 model was calibrated by using available streamflow records from the model area. The SGPWM uses an updated INFILv3 version to represent shallow groundwater flow better beneath the root zone that contributes to lateral, downslope seepage rather than deep recharge. The SGPWM calibration was tested by using available streamflow records in the San Gorgonio Pass region. The SGPWM was used to simulate a 100-year water budget, including recharge and runoff, for water years 1913 through 2012. Results indicated that most recharge came from episodic infiltration of surface-water runoff in the larger stream channels. Results also indicated periods of great variability in recharge and runoff in response to variability in precipitation. More recharge was simulated for the area of the groundwater basin underlying the more permeable alluvial fill of the valley floor compared to recharge in the neighboring upland areas of the less permeable mountain blocks. The greater recharge was in response to the episodic streamflow that discharged from the mountain block areas and quickly infiltrated the permeable alluvial fill of the groundwater basin. Although precipitation at the higher altitudes of the mountain block was more than double precipitation at the lower altitudes of the valley floor, recharge for inter-channel areas of the mountain block was limited by the lower permeability bedrock underlying the thin soil cover, and most of the recharge in the mountain block was limited to the main stream channels underlain by alluvial fill.

California↗
USGS2020

Simulation of discharge, water-surface elevations, and water temperatures for the St. Louis River estuary, Minnesota-Wisconsin, 2016–17

The St. Louis River estuary is a large freshwater estuary, next to Duluth, Minnesota, that encompasses the headwaters of Lake Superior. The St. Louis River estuary is one of the most complex and compromised near-shore systems in the upper Great Lakes with a long history of environmental contamination caused by logging, mining, paper mills, and other heavy industrial activities. Presently (2020), a widely available, science-based assessment tool capable of evaluating ecosystem-level responses to remediation and restoration projects has not existed for the estuary. To address this need, the U.S. Geological Survey (USGS) built a predictive, mechanistic, three-dimensional hydrodynamic model for the estuary using the Environmental Fluid Dynamics Code framework. In the current version, the model can simulate continuous discharge, water-surface elevations, water temperature, and flow velocity, although the modular framework allows for future additions of water-quality modeling. The model was calibrated using data collected from April 2016 through November 2016 and validated with data collected from April 2017 through November 2017. The four types of data used to evaluate model performance were water-surface elevations, discharge, water temperature, and flow velocities. Streamflow and temperature boundary condition data included a mixture of USGS streamgage data, Minnesota Department of Natural Resources gage data, and estimates derived from the gage data. The model was able to simulate the water-surface elevations with generally good agreement between the simulated and measured values for both years at the daily time step. Specifically, the model was able to demonstrate excellent agreement with the measured data with Nash-Sutcliffe efficiency coefficients greater than 0.8 for all three locations; however, the model was unable to produce hourly water-surface elevations with such accuracy for 2016–17. Discharge was more dynamic than the water-surface elevations, both for the measured and simulated data. Generally, most of the discharge ranged from −650 to 1,200 cubic meters per second, but the constantly changing flux exiting the estuary into Lake Superior (positive flows) and entering the estuary from Lake Superior (negative flows) occurred throughout the year. Even upstream at the St. Louis River at Oliver, Wisconsin, gage (USGS station 0402403250), the effect of flows into the estuary from Lake Superior did occur, demonstrating the strong effect of the Lake Superior seiche on flows for the estuary. From a performance standpoint, the model was able to simulate discharge with generally good agreement in both years, although the 2017 validation was better than the 2016 calibration period. For the daily Nash-Sutcliffe efficiency coefficients, the simulated values were 0.98, 0.62, 0.49, and 0.71 for the Oliver gage; the Superior Bay entry channel at Superior, Wisc., (USGS station 464226092005600); the Superior Bay Duluth Ship Canal at Duluth, Minn., (USGS station 464646092052900); and total entries (combination of the Superior entry and Duluth entry), respectively. For the hourly evaluation criteria, the model performed poorly, with Nash-Sutcliffe efficiency coefficients less than 0 for the two entries into Lake Superior; therefore, as a predictor of discharge at the hourly scale, the model performed worse than using the measured data average. Similar to discharge, the model was a good predictor of flow velocity at the daily time scale but had difficulty matching the measured data at the hourly scale. For discharge and flow velocity, matching at subdaily time steps for a system as complicated as the St. Louis River estuary is considered difficult because the match is highly sensitive to coordinating the exact measurement location to the simulated value. The final calibration target was water temperature, calibrated for the Oliver gage and the Duluth entry. For calibration purposes, the Duluth entry was the more important water temperature target because the Oliver gage was more of an internal check on the model. The Nash-Sutcliffe efficiency coefficients for the Duluth entry were high; hourly Nash-Sutcliffe efficiency coefficients at the Duluth entry were either at or greater than 0.7 for both years, and daily values were 0.84 and 0.82 for 2016 and 2017, respectively.

Minnesota, Wisconsin↗
USGS2009

Occurrence of volatile organic compounds in selected urban streams in the United States, 1995-2003

As part of the U.S. Geological Survey's (USGS) National Water-Quality Assessment (NAWQA) Program, urban indicator sites were monitored to (1) characterize the stream quality from drainage basins with predominantly residential and commercial land use, and (2) determine which selected natural and anthropogenic factors affect stream quality. A total of 869 water samples were collected from 37 urban streams during 1995-2003 and were analyzed for 87 volatile organic compounds (VOCs). The occurrence of VOCs in urban streams is described in this report for (1) all samples as a single dataset, (2) all samples grouped by streamflow pentiles, and (3) all samples grouped by warmer (April through September) and cooler (October through March) months by the detection frequency and (or) concentration of (a) any VOC, (b) VOC groups, and (c) individual compounds. An assessment level of 0.02 microgram per liter (ug/L) was used to compute the detection frequencies and concentrations of VOCs. Concentrations of VOCs were compared to (1) U.S. Environmental Protection Agency's (USEPA) drinking-water Maximum Contaminant Levels (MCLs) or Drinking Water Advisories, (2) Health-Based Screening Levels (HBSLs) developed by the USGS in collaboration with the USEPA and other agencies, and (3) USEPA and Canadian aquatic-life criteria. One or more VOCs were detected in 97.1 percent of 869 samples, and one or more VOCs were detected frequently (greater than 80 percent) at all sites. The median total VOC concentration for all samples was 0.57 ug/L, and total VOC concentrations in a single sample ranged from not detected to 698 ug/L. About 85 percent of the samples contained two or more VOCs, and about one-half contained five or more VOCs. The gasoline hydrocarbons were the most frequently occurring VOC group followed by solvents, trihalomethanes (THMs), gasoline oxygenates, organic synthesis compounds, fumigants, and refrigerants. Concentration ranges for most VOC groups were distributed over at least two orders of magnitude. Fifty-seven of the 87 VOCs analyzed were detected in at least one sample at an assessment level of 0.02 ug/L. More than one-half of the 30 VOCs not detected in samples were organic synthesis compounds. Fifteen compounds had detection frequencies greater than or equal to 10 percent. With the exception of toluene and chloroform, the median concentration of each VOC for all samples was less than the assessment level. Furthermore, the median concentrations of detections for the 15 most frequently occurring VOCs ranged from 0.03 to 3.9 ug/L, and typically were less than or equal to 0.10 ug/L. The 869 samples from the 37 sites were stratified into five streamflow pentiles (less than 20, 20-less than 40, 40-less than 60, 60-less than 80, and greater than or equal to 80 percent of estimated long-term streamflow statistics) for comparison of the occurrence of VOCs. The detection frequency of one or more VOCs by streamflow pentile varied only slightly from 96.7 to 97.7 percent. The median total VOC concentrations in samples for the five streamflow pentiles ranged from 0.39 to 1.0 ug/L. Two or more VOCs were present in more than 80 percent of samples in each of the five pentiles. The gasoline hydrocarbons, solvents, THMs, and gasoline oxygenates occurred frequently (greater than 30 percent) in all streamflow pentiles, in contrast to the organic synthesis compounds, fumigants, and refrigerants that occurred less frequently in urban streams under all streamflow conditions. The median total VOC concentrations for gasoline hydrocarbons, solvents, gasoline oxygenates, and organic synthesis compounds generally increased as streamflow increased. In contrast, the median total VOC concentrations for THMs and fumigants generally decreased as streamflow increased. The median total VOC concentrations for refrigerants showed no pattern as streamflow increased. Because differences between VOC occurrence and streamflow pentiles were small for most compariso

Scientific Investigations Report↗
USGS2010

Effects of Changes in Irrigation Practices and Aquifer Development on Groundwater Discharge to the Jobos Bay National Estuarine Research Reserve near Salinas, Puerto Rico

Since 1990, about 75 acres of black mangroves have died in the Jobos Bay National Estuarine Research Reserve near Salinas, Puerto Rico. Although many factors can contribute to the mortality of mangroves, changes in irrigation practices, rainfall, and water use resulted in as much as 25 feet of drawdown in the potentiometric surface of the aquifer in the vicinity of the reserve between 1986 and 2002. To clarify the issue, the U.S. Geological Survey, in cooperation with the Puerto Rico Department of Natural and Environmental Resources, conducted a study to ascertain how aquifer development and changes in irrigation practices have affected groundwater levels and groundwater flow to the Mar Negro area of the reserve. Changes in groundwater flow to the mangrove swamp and bay from 1986 to 2004 were estimated in this study by developing and calibrating a numerical groundwater flow model. The transient simulations indicate that prior to 1994, high irrigation return flows more than offset the effect of reduced groundwater withdrawals. In this case, the simulated discharge to the coast in the modeled area was 19 million gallons per day. From 1994 through 2004, furrow irrigation was completely replaced by micro-drip irrigation, thus eliminating return flows and the simulated average coastal discharge was 7 million gallons per day, a reduction of 63 percent. The simulated average groundwater discharge to the coastal mangrove swamps in the reserve from 1986 to 1993 was 2 million gallons per day, compared to an average simulated discharge of 0.2 million gallons per day from 1994 to 2004. The average annual rainfall for each of these periods was 38 inches. The groundwater discharge to the coastal mangrove swamps in the Jobos Bay National Estuarine Research Reserve was estimated at about 0.5 million gallons per day for 2003-2004 because of higher than average annual rainfall during these 2 years. The groundwater flow model was used to test five alternatives for increasing groundwater discharge to the coastal mangrove swamps to approximately 1.4 million gallons per day: (1) artificially recharging the aquifer with injection wells or (2) by increasing irrigation return flow by going back to furrow irrigation; (3) termination of groundwater withdrawals near the mangroves; (4) reduction of groundwater withdrawals at irrigation wells by 50 percent; and (5) a combination of alternatives 2 and 4 increasing irrigation return flows and decreasing irrigation withdrawals. Each alternative assumed average climatic conditions and groundwater withdrawals at 2004 rates. Alternative 1 required 1.5 million gallons per day of injected water. Alternative 2 required flooding 958 acres with a rate of 1.84 million gallons per day if no crops are grown. Alternative 3 required the termination of 2.44 million gallons per day of withdrawals to achieve 1.34 million gallons per day of discharge to the mangroves. Alternative 4 did not achieve the objective with only 0.80 million gallons per day simulated discharge to the mangroves, while requiring a 1.26 million gallon per day reduction in groundwater withdrawals. Alternative 5 required flooding fields with additional 1.13 million gallons of day and the same reduction in groundwater withdrawals, but did achieve the objective of about 1.4 million gallons per day discharge to the mangroves. Alternative 1, incorporating injection wells near the reserve required the least amount of water to raise groundwater levels and maintain discharge of 1.4 million gallons per day through the mangroves.

Scientific Investigations Report↗
USGS2007

Simulation of Surface-Water Conditions in the Nontidal Passaic River Basin, New Jersey

The Passaic River Basin, the third largest drainage basin in New Jersey, encompasses 950 mi2 (square miles) in the highly urbanized area outside New York City, with a population of 2 million. Water quality in the basin is affected by many natural and anthropogenic factors. Nutrient loading to the Wanaque Reservoir in the northern part of the basin is of particular concern and is caused partly by the diversion of water at two downstream intakes that is transferred back upstream to refill the reservoir. The larger of these diversions, Wanaque South intake, is on the lower Pompton River near Two Bridges, New Jersey. To support the development of a Total Maximum Daily Load (TMDL) for nutrients in the nontidal part of the basin (805 mi2), a water-quality transport model was needed. The U.S. Geological Survey, in cooperation with the New Jersey Department of Environmental Protection and New Jersey EcoComplex, developed a flow-routing model to provide the hydraulic inputs to the water-quality model. The Diffusion Analogy Flow model (DAFLOW) described herein was designed for integration with the Water Quality Analysis Simulation Program (WASP) watershed water-quality model. The flow routing model was used to simulate flow in 108 miles of the Passaic River and major tributaries. Flow data from U.S. Geological Survey streamflow-gaging stations represent most of the model's upstream boundaries. Other model inputs include estimated flows for ungaged tributaries and unchanneled drainage along the mainstem, and reported flows for major point-source discharges and diversions. The former flows were calibrated using the drainage-area ratio method. The simulation extended over a 4+ year period representing a range in flow conditions. Simulated channel cross-sectional geometry in the DAFLOW model was calibrated using several different approaches by adjusting area and top width parameters. The model also was calibrated to observed flows for water year 2001 (low flow) at five mainstem gaging stations and one station at which flow was estimated. The model's target range was medium to low flows--the range of typical intake operations. Simulated flow mass balance, hydrographs (flood-wave speed, attenuation, and spread), flow-duration curves, and velocity and depth values were compared to observed counterparts. Mass balance and hydrograph fit were evaluated quantitatively. Simulation results generally were within the accuracy of the flow data at the measurement stations. The model was validated to observed flows for water years 2000 (average flow), 2002 (extreme low flow), and 2003 (high flow). Results for 19 of 20 comparisons indicate average mass-balance and model-fit errors of 6.6 and 15.7 percent, respectively, indicating that the model reasonably represents the time variation of streamflow in the nontidal Passaic River Basin. An algorithm (subroutine) also was developed for DAFLOW to simulate the hydraulic mixing that occurs near the Wanaque South intake upstream from the confluence of the Pompton and Passaic Rivers. The intake draws water from multiple sources, including effluent from a nearby wastewater-treatment plant, all of which have different phosphorus loads. The algorithm determines the proportion of flow from each source and operates within a narrow flow range. The equations used in the algorithm are based on the theory of diffusion and lateral mixing in rivers. Parameters used in the equations were estimated from limited available local flow and water-quality data. As expected, simulation results for water years 2000, 2001, and 2003 indicate that most of the water drawn to the intake comes from the Pompton River; however, during many short periods of low flow and high diversion, particularly in water year 2002, entrainment of the other flow sources compensated for the insufficient flow in the Pompton River. As additional verification of the flow model used in the water-quality model, a Branched Lagrangian Transport Model (B

Scientific Investigations Report↗
USGS2008

Quantifying Ground-Water and Surface-Water Discharge from Evapotranspiration Processes in 12 Hydrographic Areas of the Colorado Regional Ground-Water Flow System, Nevada, Utah, and Arizona

Rapid population growth in southern Nevada has increased the demand for additional water supplies from rural areas of northern Clark and southern Lincoln counties to meet projected water-supply needs. Springs and rivers in these undeveloped areas sustain fragile riparian habitat and may be susceptible to ground-water withdrawals. Most natural ground-water and surface-water discharge from these basins occurs by evapotranspiration (ET) along narrow riparian corridors that encompassed about 45,000 acres or about 1 percent of the study area. This report presents estimates of ground- and surface-water discharge from ET across 3.5 million acres in 12 hydrographic areas of the Colorado Regional Ground-Water Flow System. Ground-and surface-water discharge from ET were determined by identifying areas of ground- and surface-water ET, delineating areas of similar vegetation and soil conditions (ET units), and computing ET rates for each of these ET units. Eight ET units were identified using spectral-reflectance characteristics determined from 2003 satellite imagery, high-resolution aerial photography, and land classification cover. These ET units are dense meadowland vegetation (200 acres), dense woodland vegetation (7,200 acres), moderate woodland vegetation (6,100 acres), dense shrubland vegetation (5,800 acres), moderate shrubland vegetation (22,600 acres), agricultural fields (3,100 acres), non-phreatophytic areas (3,400,000 acres), and open water (300 acres). ET from diffuse ground-water and channelized surface-water is expressed as ETgs and is equal to the difference between total annual ET and precipitation. Total annual ET rates were calculated by the Bowen ratio and eddy covariance methods using micrometeorological data collected from four sites and estimated at 3.9 ft at a dense woodland site (February 2003 to March 2005), 3.6 ft at a moderate woodland site (July 2003 to October 2006), 2.8 ft at a dense shrubland site (June 2005 to October 2006), and 1.5 ft at a moderate shrubland site (April 2006 to October 2006). Annual ETgs rates were 3.4 ft for dense woodland vegetation, 3.2 ft for moderate woodland vegetation, 2.2 ft for dense shrubland vegetation, and 1.0 ft for moderate shrubland vegetation. Published annual rates of ETgs were used for the other ET units found in the study area. These rates were 3.4 ft for dense meadowland vegetation, 5.2 ft for agricultural fields, and 4.9 ft for open water. For the non-phreatophytic ET unit, ETgs was assumed to be zero. Estimated ground- and surface-water discharge from ET was calculated by multiplying the ETgs by the ET-unit acreage and equaled 24,480 acre-ft for dense woodland vegetation, 19,520 acre-ft for moderate woodland vegetation, 12,760 acre-ft for dense shrubland vegetation, 22,600 acre-ft for moderate shrubland vegetation, 680 acre-ft for dense meadowland vegetation, 16,120 acre-ft for agricultural fields, 1,440 acre-ft for open water, and 0 acre-ft for the non-phreatophytic ET unit. Estimated ground-water and surface-water discharge from ET from each hydrographic area was calculated by summing the total annual ETgs rate for ET units found within each hydrographic area and equaled 1,952 acre-ft for the Black Mountains Area, 6,080 acre-ft for California Wash, 4,090 acre-ft for the Muddy River Springs Area, 11,510 acre-ft for Lower Moapa Valley, 51,960 acre-ft for the Virgin River Valley, 16,168 acre-ft for Lower Meadow Valley Wash, 5,840 acre-ft for Clover Valley, and 0 acre-ft for Coyote Spring Valley, Kane Springs Valley, Tule Desert, Hidden Valley (North), and Garnet Valley. The annual discharge from ETgs for the study area totals about 98,000 acre-ft.

Scientific Investigations Report↗
USGS2009

Water Quality and Hydrology of Whitefish (Bardon) Lake, Douglas County, Wisconsin, With Special Emphasis on Responses of an Oligotrophic Seepage Lake to Changes in Phosphorus Loading and Water Level

Whitefish Lake, which is officially named Bardon Lake, is an oligotrophic, soft-water seepage lake in northwestern Wisconsin, and classified by the Wisconsin Department of Natural Resources as an Outstanding Resource Water. Ongoing monitoring of the lake demonstrated that its water quality began to degrade (increased phosphorus and chlorophyll a concentrations) around 2002 following a period of high water level. To provide a better understanding of what caused the degradation in water quality, and provide information to better understand the lake and protect it from future degradation, the U.S. Geological Survey did a detailed study from 2004 to 2008. The goals of the study were to describe the past and present water quality of the lake, quantify water and phosphorus budgets for the lake, simulate the potential effects of changes in phosphorus inputs on the lake's water quality, analyze changes in the water level in the lake since 1900, and relate the importance of changes in climate and changes in anthropogenic (human-induced) factors in the watershed to the water quality of the lake. Since 1998, total phosphorus concentrations increased from near the 0.005-milligrams per liter (mg/L) detection limit to about 0.010 mg/L in 2006, and then decreased slightly in 2007-08. During this time, chlorophyll a concentrations and Secchi depths remained relatively stable at about 1.5 micrograms per liter (ug/L) and 26 feet, respectively. Whitefish Lake is typically classified as oligotrophic. Because the productivity in Whitefish Lake is limited by phosphorus, phosphorus budgets were constructed for the lake. Because it was believed that much of its phosphorus comes from the atmosphere, phosphorus deposition was measured in this study. Phosphorus input from the atmosphere was greater than computed based on previously reported wetfall phosphorus concentrations. The concentrations and deposition rates can be used to estimate atmospheric loading in future lake studies. The average annual load of phosphorus to the lake was 232 pounds: 56 percent from precipitation, 27 percent from groundwater, and 16 percent from septic systems. During a series of dry years (low water levels) and wet years (high water levels), the inputs of water and phosphorus ranged by only 10-13 percent. Results from the Canfield and Bachmann eutrophication model and Carlson trophic-state-index equations demonstrated that the lake directly responds to changes in external phosphorus loading, with percent change in chlorophyll a being similar to the percent change in loading and the change in total phosphorus and Secchi depth being slightly smaller. Therefore, changes in phosphorus loading should affect the water quality of the lake. Specific scenarios that simulated the effects of anthropogenic (human-induced) and climatic (water level) changes demonstrated that: surface-water inflow (runoff) based on current development has little effect on pelagic water quality, changes in the inputs from septic systems and development in the watershed could have a large effect on water quality, and decreases in water and phosphorus loading during periods of low water level had little effect on water quality. Sustained high water levels, resulting from several wet years with relatively high water and phosphorus input, however, could cause a small degradation in water quality. Although high water levels may be associated with a degradation in water quality, it appears that anthropogenic changes in the watershed may be more important in affecting the future water quality of the lake. Fluctuations in water levels since 1998 are representative of what has occurred since 1900, with fluctuations of about 3 feet occurring about every 15 years. Based on total phosphorus concentrations inferred from sediment core analysis, there has been little long-term change in water quality and there has been a slight deterioration in water quality following most periods of high water levels. There

Scientific Investigations Report↗
USGS2007

Simulation of Multiscale Ground-Water Flow in Part of the Northeastern San Joaquin Valley, California

The transport and fate of agricultural chemicals in a variety of environmental settings is being evaluated as part of the U.S. Geological Survey (USGS) National Water-Quality Assessment Program. One of the locations being evaluated is a 2,700-km2 (square kilometer) regional study area in the northeastern San Joaquin Valley surrounding the city of Modesto, an area dominated by irrigated agriculture in a semi-arid climate. Ground water is a key source of water for irrigation and public supply, and exploitation of this resource has altered the natural flow system. The aquifer system is predominantly alluvial, and an unconfined to semiconfined aquifer overlies a confined aquifer in the southwestern part of the study area; these aquifers are separated by the lacustrine Corcoran Clay. A regional-scale 16-layer steady-state model of ground-water flow in the aquifer system in the regional study area was developed to provide boundary conditions for an embedded 110-layer steady-state local-scale model of part of the aquifer system overlying the Corcoran Clay along the Merced River. The purpose of the local-scale model was to develop a better understanding of the aquifer system and to provide a basis for simulation of reactive transport of agricultural chemicals. The heterogeneity of aquifer materials was explicitly incorporated into the regional and local models using information from geologic and drillers? logs of boreholes. Aquifer materials were differentiated in the regional model by the percentage of coarse-grained sediments in a cell, and in the local model by four hydrofacies (sand, silty sand, silt, and clay). The calibrated horizontal hydraulic conductivity values of the coarse-grained materials in the zone above the Corcoran Clay in the regional model and of the sand hydrofacies used in the local model were about equal (30?80 m/d [meter per day]), and the vertical hydraulic conductivity values in the same zone of the regional model (median of 0.012 m/d), which is dominated by the finer-grained materials, were about an order of magnitude less than that for the clay hydrofacies in the local model. Data used for calibrating both models included long-term hourly water-level measurements in 20 short-screened wells installed by the USGS in the Modesto and Merced River areas. Additional calibration data for the regional model included water-level measurements in 11 wells upslope and 17 wells downslope from these areas. The root mean square error was 2.3 m (meter) for all wells in the regional model and 0.8 m for only the USGS wells; the associated average errors were 0.9 m and 0.3 m, respectively. The root mean square error for the 12 USGS wells along a transect in the local model area was 0.08 m; the average error was 0.0 m. Particle tracking was used with the local model to estimate the concentration of an environmental tracer, sulfur hexafluoride, in 10 USGS transect wells near the Merced River that were sampled for this constituent. Measured and estimated concentrations in the mid-depth and deepest wells, which would be most sensitive to errors in hydraulic conductivity estimates, were consistent. The combined results of particle tracking and sulfur hexafluoride analysis suggest that most water sampled from the transect wells was recharged less that 25 years ago.

Scientific Investigations Report↗
USGS2009

Surface-Water and Groundwater Interactions along the Withlacoochee River, West-Central Florida

A study of the Withlacoochee River watershed in west-central Florida was conducted from October 2003 to March 2007 to gain a better understanding of the hydrology and surface-water and groundwater interactions along the river. The Withlacoochee River originates in the Green Swamp area in north-central Polk County and flows northerly through seven counties, emptying into the Gulf of Mexico. This study includes only the part of the watershed located between the headwaters in the Green Swamp and the U.S. Geological Survey gaging station near Holder, Florida. The Withlacoochee River within the study area is about 108 miles long and drains about 1,820 square miles. The Withlacoochee River watershed is underlain by thick sequences of carbonate rock that are covered by thin surficial deposits of unconsolidated sand and sandy clay. The clay layer is breached in many places because of the karst nature of the underlying limestone, and the degree of confinement between the Upper Florida aquifer and the surficial aquifer is highly variable throughout the watershed. The potential for movement of water from the surface or shallow deposits to deeper deposits, or from deeper deposits to the shallow deposits, exists throughout the Withlacoochee River watershed. Water levels were higher in deeper Upper Floridan aquifer wells than in shallow Upper Floridan aquifer wells or surficial aquifer wells at 11 of 19 paired or nested well sites, indicating potential for discharge to the surface-water system. Water levels were higher in shallow Upper Floridan aquifer or surficial aquifer wells than in deeper Upper Floridan aquifer wells at five other sites, indicating potential for recharge to the deeper Upper Floridan aquifer. Water levels in the surficial aquifer and Upper Floridan aquifer wells at the remaining three sites were virtually the same, indicating little or no confinement at the sites. Potentiometric-surface maps of the Upper Floridan aquifer indicate the pattern of groundwater flow in the aquifer did not vary greatly from season to season during the study. Potentiometric contours indicate groundwater discharge to the river in the vicinity of Dade City and Lake Panasoffkee. During dry periods, groundwater from the underlying Upper Floridan aquifer contributed to the flow in the river. During wet periods, streamflow had additional contributions from runoff and input from tributaries. Groundwater has a greater effect on streamflow downstream from the Dade City station than upstream from the Dade City station because confinement between surficial deposits and the Upper Floridan aquifer is greater in the Green Swamp area than in downstream areas. Estimates of streamflow gains and losses were made along the Withlacoochee River during base-flow conditions in May 2004, April 2005, and April 2006. Base flow was higher in April 2005 than in May 2004 and April 2006. Consistent net seepage gains were identified in 16 of 20 subreaches analyzed during all seepage runs. The direction of exchange was variable in the remaining four subreaches. Low specific conductance, pH, and calcium concentrations in water from the Withlacoochee River near the headwater area indicated a surface-water system not directly connected to the Upper Floridan aquifer. Downstream from the Dade City station, higher specific conductance, pH, and calcium concentrations in the river water indicated an increasing influence of groundwater, and were similar to groundwater during low-flow conditions. Strontium isotope ratios indicate groundwater originates from shallow parts of the Upper Floridan aquifer in the upper reaches of the river, and from increasingly deeper parts of the aquifer in the downstream direction. Mean annual base-flow estimates also indicate increasing groundwater discharge to the river in the downstream direction. Mean annual base flow estimated using standard hydrograph separation method assumptions ranged from about 4.7 to 5.1 inches per year

Scientific Investigations Report↗