Search USGS⌕ Search

SEARCH · Search USGS

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 739 records · Page 41Linked to original sources

Conceptual Model of Hydrologic and Thermal Conditions of the Eastbank Aquifer System near Rocky Reach Dam, Douglas County, Washington

The Lower and Combined Aquifers of the Eastbank Aquifer system, located in a river-terrace deposit along the Columbia River near Rocky Reach Dam, Washington, are primarily recharged by the Columbia River and provide water to the Eastbank Hatchery and the regional water system servicing the cities of Wenatchee, East Wenatchee, and parts of unincorporated Chelan and Douglas Counties. In 2006, mean annual pumpage from the aquifers by the hatchery and regional water system was about 43 and 16 cubic feet per second, respectively. Reportedly, temperatures of ground water pumped by the hatchery have been increasing, thereby making water potentially too warm for salmonid fish production. An evaluation of hourly ground-water and river temperatures from January 1991 through August 2007 indicates increasing interannual trends in temperatures in most of the Lower and Combined Aquifers from 1999 through 2006 that correspond to increasing trends in the annual mean and annual maximum river temperatures during the same period of 0.07 and 0.17?C per year, respectively. There were no trends in the annual minimum river temperatures from 1999 through 2006, and there were no trends in the annual minimum, mean, and maximum river temperatures from 1991 through 1998 and from 1991 through 2007. Increases in river temperatures from 1999 through 2006 are within the natural variability of the river temperatures. Most of the Lower and Combined Aquifers reached thermal equilibrium?defined by constant time lags between changes in river temperatures and subsequent changes in ground-water temperatures?during 1991?98. The only exceptions are the Combined Aquifer north of the well field of the regional water system, which had not reached thermal equilibrium by 2006, and the Lower Aquifer west of the well fields of the hatchery and the regional water system, which reached thermal equilibrium prior to 1991. Because most of the Lower and Combined Aquifers were in thermal equilibrium from 1999 through 2006 and seasonal pumpage patterns were relatively stable, reported trends of increasing temperatures of water pumped by the hatchery well field are most likely explained by increasing trends in river temperatures. Most of the water pumped by the hatchery well field recharges in an area west to southwest of the well field about 2 months prior to the time it is pumped from the aquifer. The northern extent of the hatchery well field may pump some colder water from a bedrock depression to the north and west of the well field. The conceptual model of hydrologic and thermal conditions is supported by analyses of historical water temperatures, water-level data collected on July 18, 2007, and dissolved-constituent and bacterial concentrations in samples collected on August 20?22, 2007.

Scientific Investigations Report↗

TOPMODEL simulations of streamflow and depth to water table in Fishing Brook Watershed, New York, 2007-09

TOPMODEL, a physically based, variable-source area rainfall-runoff model, was used to simulate streamflow and depth to water table for the period January 2007-September 2009 in the 65.6 square kilometers of Fishing Brook Watershed in northern New York. The Fishing Brook Watershed is located in the headwaters of the Hudson River and is predominantly forested with a humid, cool continental climate. The motivation for applying this model at Fishing Brook was to provide a simulation that would be effective later at this site in modeling the interaction of hydrologic processes with mercury dynamics. TOPMODEL uses a topographic wetness index computed from surface-elevation data to simulate streamflow and subsurface-saturation state, represented by the saturation deficit. Depth to water table was computed from simulated saturation-deficit values using computed soil properties. In the Fishing Brook Watershed, TOPMODEL was calibrated to the natural logarithm of streamflow at the study area outlet and depth to water table at Sixmile Wetland using a combined multiple-objective function. Runoff and depth to water table responded differently to some of the model parameters, and the combined multiple-objective function balanced the goodness-of-fit of the model realizations with respect to these parameters. Results show that TOPMODEL reasonably simulated runoff and depth to water table during the study period. The simulated runoff had a Nash-Sutcliffe efficiency of 0.738, but the model underpredicted total runoff by 14 percent. Depth to water table computed from simulated saturation-deficit values matched observed water-table depth moderately well; the root mean squared error of absolute depth to water table was 91 millimeters (mm), compared to the mean observed depth to water table of 205 mm. The correlation coefficient for temporal depth-to-water-table fluctuations was 0.624. The variability of the TOPMODEL simulations was assessed using prediction intervals grouped using the combined multiple-objective function. The calibrated TOPMODEL results for the entire study area were applied to several subwatersheds within the study area using computed hydrogeomorphic properties of the subwatersheds.

Scientific Investigations Report↗

Navigability potential of Washington rivers and streams determined with hydraulic geometry and a geographic information system

Using discharge and channel geometry measurements from U.S. Geological Survey streamflow-gaging stations and data from a geographic information system, regression relations were derived to predict river depth, top width, and bottom width as a function of mean annual discharge for rivers in the State of Washington. A new technique also was proposed to determine bottom width in channels, a parameter that has received relatively little attention in the geomorphology literature. These regression equations, when combined with estimates of mean annual discharge available in the National Hydrography Dataset, enabled the prediction of hydraulic geometry for any stream or river in the State of Washington. Predictions of hydraulic geometry can then be compared to thresholds established by the Washington State Department of Natural Resources to determine navigability potential of rivers. Rivers with a mean annual discharge of 1,660 cubic feet per second or greater are 'probably navigable' and rivers with a mean annual discharge of 360 cubic feet per second or less are 'probably not navigable'. Variance in the dataset, however, leads to a relatively wide range of prediction intervals. For example, although the predicted hydraulic depth at a mean annual discharge of 1,660 cubic feet per second is 3.5 feet, 90-percent prediction intervals indicate that the actual hydraulic depth may range from 1.8 to 7.0 feet. This methodology does not determine navigability - a legal concept determined by federal common law - instead, this methodology is a tool for predicting channel depth, top width, and bottom width for rivers and streams in Washington.

Washington↗

Volcanic aquifers of Hawai‘i—Hydrogeology, water budgets, and conceptual models

Hawai‘i’s aquifers have limited capacity to store fresh groundwater because each island is small and surrounded by saltwater. Saltwater also underlies much of the fresh groundwater. Fresh groundwater resources are, therefore, particularly vulnerable to human activity, short-term climate cycles, and long-term climate change. Availability of fresh groundwater for human use is constrained by the degree to which the impacts of withdrawal—such as lowering of the water table, saltwater intrusion, and reduction in the natural discharge to springs, streams, wetlands, and submarine seeps—are deemed acceptable. This report describes the hydrogeologic framework, groundwater budgets (inflows and outflows), conceptual models of groundwater occurrence and movement, and the factors limiting groundwater availability for the largest and most populated of the Hawaiian Islands—Kaua‘i, O‘ahu, Maui, and Hawai‘i Island. The bulk of each of Hawai‘i’s islands is built of many thin lava flows erupted from shield volcanoes; the great piles of lava flows form highly permeable aquifers. In some areas, low-permeability dikes cutting across the lava flows, or low-permeability ash and soil horizons interlayered with the lava flows, can substantially alter groundwater flow. On some islands, sedimentary rocks form thick semiconfining coastal-plain deposits, locally known as caprock, that impede natural groundwater discharge to the ocean. In some regions, thick lava flows that ponded in preexisting depressions form aquifers that are much less permeable than aquifers formed by thin lava flows. Fresh groundwater inflow to Hawai‘i’s aquifers comes from recharge. For predevelopment conditions (1870), estimates of groundwater recharge from this study are 871, 675, 1,279, and 5,291 million gallons per day (Mgal/d) for Kaua‘i, O‘ahu, Maui, and Hawai‘i Island, respectively. Estimates of recharge for recent conditions (2010 land cover and 1978–2007 rainfall for Kaua‘i, O‘ahu, and Maui; 2008 land cover and 1916–1983 rainfall for Hawai‘i Island) are 875, 660, 1,308, and 6,595 Mgal/d for Kaua‘i, O‘ahu, Maui, and Hawai‘i Island, respectively. Recent recharge values differ from predevelopment recharge values by only a few percent for all islands except Hawai‘i Island, where changes in forest cover affected recharge. Spatial distribution of recharge mimics the orographic rainfall pattern—recharge is high on windward slopes and mountain peaks below the top of the trade-wind inversion. Human activity such as irrigation also contributes to recharge in some areas. Outflows from Hawai‘i’s aquifers include withdrawals from wells and natural groundwater discharge to springs, streams, wetlands, and submarine seeps. Under predevelopment conditions, groundwater withdrawal is assumed to be negligible and natural groundwater discharge probably was equal, or close, to recharge. Under recent conditions (2000–2010), groundwater withdrawal averaged 19, 209, 104, and 103 Mgal/d on Kaua‘i, O‘ahu, Maui, and Hawai‘i Island, respectively. If recent withdrawal and recharge rates are maintained until steady state is achieved, natural groundwater discharge will be reduced by an amount equal to the withdrawal rate. Total recent withdrawal for the four islands is only about 5 percent of total recharge, but about half of the withdrawal comes from O‘ahu, whereas O‘ahu receives only 7 percent of the total recharge. Effects of high withdrawals on O‘ahu cannot be mitigated by the lower withdrawals on other islands because no freshwater flows between islands. Even within an island, high withdrawals from one area cannot be completely mitigated by recharge in another area. Water-level, saltwater/freshwater-transition-zone, spring, and stream base-flow data indicate an overall reduction in storage for most areas where groundwater has been developed. Groundwater occurrence and movement in Hawai‘i’s volcanic aquifers can be described in terms of four conceptual models: (1) fresh groundwater lenses in high-permeability lava-flow aquifers, (2) aquifers with groundwater impounded by dikes, (3) thickly saturated low-permeability aquifers, and (4) perched aquifers. In Hawai‘i, most fresh groundwater withdrawn for human use comes from freshwater lenses in the dike-free high-permeability lava-flow aquifers where the principal limiting factor to groundwater availability is saltwater intrusion, but impacts of reduced natural groundwater discharge may also limit availability. Dike-impounded groundwater is common near the center of Hawaiian shield volcanoes, where water moves and is stored in permeable lava flows between the dikes; groundwater availability in these aquifers is primarily limited by storage depletion and reduction of flow to adjacent aquifers and natural groundwater discharge. Thickly saturated low-permeability aquifers have been identified on Kaua‘i and Maui; groundwater availability is primarily limited by streamflow depletion and water-table decline. Perched groundwater is postulated to exist in some areas of Hawai‘i, but store much less water than other modes ofgroundwater occurrence. Limits on groundwater availability in perched aquifers include the potential of reducing inflow to other groundwater settings and reducing natural discharge and stream seepage. Some groundwater bodies in Hawai‘i are enigmatic; consequences of groundwater development in these bodies and their relation to groundwater availability are not completely understood.

Hawai'i↗

Characterization of dissolved and particulate natural organic matter (NOM) in Neversink Reservoir, New York

Natural organic matter (NOM) was isolated from the water of the Neversink Reservoir, part of the New York City water supply, located in the Catskill Mountains of New York. The NOM was fractionated into the following nine different fractions by the isolation procedure: (1) coarse particulates, (2) fine-particulate organics, (3) solvent-extractable organics, (4) hydrophobic neutrals (HPON fraction), (5) dissolved colloids, (6) bases, (7) hydrophobic acids (HPOA), (8) transphilic acids + neutrals (TPI-A+N), and (9) hydrophilic acids + neutrals (HPI-A+N). Each of these fractions, with exception of the first and the third which were too small for the complete series of analyses, was characterized by elemental, carbohydrate, and amino acid analyses, and by nuclear magnetic resonance and infrared spectrometry. The data obtained from these analyses indicate (1) that the fine-particulate organics and colloids are mainly composed of peptidoglycans, and lipopolysaccharides derived from algal, bacterial, and fungal cell walls, (2) that the HPO-N fraction most likely consists of a mixture of alicyclic terpenes and carbohydrates, (3) that the HPOA fraction consists mainly of lignin components conjugated to carbohydrates, (4) that the TPI-A+N and the HPI-A+N fractions most likely represent complex mixtures of relatively low molecular weight carboxylic acids derived from terpenes, carbohydrates, and peptides, and (5) that the base fraction is composed of free amino acids, browning reaction products, and peptide fragments.

Scientific Investigations Report↗

Monitoring plant tissue nitrogen isotopes to assess nearshore inputs of nitrogen to Lake Crescent, Olympic National Park, Washington

Mats of filamentous-periphytic algae present in some nearshore areas of Lake Crescent, Olympic National Park, Washington, may indicate early stages of eutrophication from nutrient enrichment of an otherwise highly oligotrophic lake. Natural abundance ratios of stable isotopes of nitrogen (δ 15 N) measured in plant tissue growing in nearshore areas of the lake indicate that the major source of nitrogen used by these primary producing plants is derived mainly from atmospherically fixed nitrogen in an undeveloped forested ecosystem. Exceptions to this pattern occurred in the Barnes Point area where elevated δ 15 N ratios indicate that effluent from septic systems also contribute nitrogen to filamentous-periphytic algae growing in the littoral zone of that area. Near the Lyre River outlet of Lake Crescent, the δ 15 N of filamentous-periphytic algae growing in close proximity to the spawning areas of a unique species of trout show little evidence of elevated δ 15 N indicating that nitrogen from on-site septic systems is not a substantial source of nitrogen for these plants. The δ 15 N data corroborate estimates that nitrogen input to Lake Crescent from septic sources is comparatively small relative to input from motor vehicle exhaust and vegetative sources in undeveloped forests, including litterfall, pollen, and symbiotic nitrogen fixation. The seasonal timing of blooms of filamentous-periphytic algal near the lake shoreline is also consistent with nitrogen exported from stands of red alder trees ( Alnus rubra ). Isotope biomonitoring of filamentous-periphytic algae may be an effective approach to monitoring the littoral zone for nutrient input to Lake Crescent from septic sources.

Washington↗

Hydrogeologic and geochemical characterization of groundwater resources in Rush Valley, Tooele County, Utah

The water resources of Rush Valley were assessed during 2008–2010 with an emphasis on refining the understanding of the groundwater-flow system and updating the groundwater budget. Surface-water resources within Rush Valley are limited and are generally used for agriculture. Groundwater is the principal water source for most other uses including supplementing irrigation. Most groundwater withdrawal in Rush Valley is from the unconsolidated basin-fill aquifer where conditions are generally unconfined near the mountain front and confined at lower altitudes near the valley center. Productive aquifers also occur in fractured bedrock along the valley margins and beneath the basin-fill deposits in some areas. Drillers’ logs and geophysical gravity data were compiled and used to delineate seven hydrogeologic units important to basin-wide groundwater movement. The principal basin-fill aquifer includes the unconsolidated Quaternary-age alluvial and lacustrine deposits of (1) the upper basin-fill aquifer unit (UBFAU) and the consolidated and semiconsolidated Tertiary-age lacustrine and alluvial deposits of (2) the lower basin-fill aquifer unit (LBFAU). Bedrock hydrogeologic units include (3) the Tertiary-age volcanic unit (VU), (4) the Pennsylvanian- to Permian-age upper carbonate aquifer unit (UCAU), (5) the upper Mississippian- to lower Pennsylvanian-age upper siliciclastic confining unit (USCU), (6) the Middle Cambrian- to Mississippian-age lower carbonate aquifer unit (LCAU), and (7) the Precambrian- to Lower Cambrian-age noncarbonate confining unit (NCCU). Most productive bedrock wells in the Rush Valley groundwater basin are in the UCAU. Average annual recharge to the Rush Valley groundwater basin is estimated to be about 39,000 acre-feet. Nearly all recharge occurs as direct infiltration of snowmelt and rainfall within the mountains with smaller amounts occurring as infiltration of streamflow and unconsumed irrigation water at or near the mountain front. Groundwater generally flows from the higher altitude recharge areas toward two distinct valley-bottom discharge areas: one in the vicinity of Rush Lake in northern Rush Valley and the other located west and north of Vernon. Average annual discharge from the Rush Valley groundwater basin is estimated to be about 43,000 acre-feet. Most discharge occurs as evapotranspiration in the valley lowlands, as discharge to springs and streams, and as withdrawal from wells. Subsurface discharge outflow to Tooele and Cedar Valleys makes up only a small fraction of natural discharge. Groundwater samples were collected from 25 sites (24 wells and one spring) for geochemical analysis. Dissolved-solids concentrations in water from these sites ranged from 181 to 1,590 milligrams per liter. Samples from seven wells contained arsenic concentrations that exceed the Environmental Protection Agency Maximum Contaminant Level of 10 micrograms per liter. The highest arsenic levels are found north of Vernon and in southeastern Rush Valley. Stable-isotope ratios of oxygen and deuterium, along with dissolved-gas recharge temperatures, indicate that nearly all modern groundwater is meteoric and derived from the infiltration of high altitude precipitation in the mountains. These data are consistent with recharge estimates made using a Basin Characterization Model of net infiltration that shows nearly all recharge occurring as infiltration of precipitation and snowmelt within the mountains surrounding Rush Valley. Tritium concentrations between 0.4 and 10 tritium units indicate the presence of modern (less than 60 years old) groundwater at 7 of the 25 sample sites. Apparent 3 H/ 3 He ages, calculated for six of these sites, range from 3 to 35 years. Adjusted minimum radiocarbon ages of premodern water samples range from about 1,600 to 42,000 years with samples from 11 of 13 sites being more than 11,000 years. These data help to identify areas where modern groundwater is circulating through the hydrologic system on time scales of decades or less and indicate that large parts of the principal basin-fill and the bedrock aquifers are much less active and receive little to no modern recharge.

Utah↗

Effects of habitat characteristics and water quality on macroinvertebrate communities along the Neversink River in southeastern New York, 1991-2001

The Neversink River, in the Catskill Mountains of southeastern New York State, feeds the Neversink Reservoir, which diverts 85 percent of the river’s flow to New York City. Acidification of several headwater reaches has affected macroinvertebrate assemblages throughout the river system above the reservoir, and the alteration of flow conditions below the reservoir dam has affected macroinvertebrate assemblages for at least 10 kilometers downstream from the reservoir. In 1999, the U.S. Geological Survey, in cooperation with The Nature Conservancy, compiled data from 30 stream reaches to quantify the effects of acidification and of the reservoir on the structure and function of macroinvertebrate assemblages throughout the Neversink River. Acidic headwater reaches supported greater numbers of acid-tolerant chironomid taxa and fewer numbers of acid-sensitive Ephemeroptera and Trichoptera than neutral reaches, and fewer scraper individuals and more shredder individuals. The 14 reaches below the reservoir, with sharply decreased flows and altered flow patterns compared to reaches above the reservoir, supported more Chironomidae and fewer Ephemeroptera and Trichoptera than the upper reaches; they also had greater numbers of shredder individuals and fewer scraper and filterer individuals than reaches above the reservoir. Water-quality variables such as pH and aluminum concentration appear to have affected macroinvertebrate assemblages more strongly in the headwaters than below the reservoir, whereas physical-habitat variables such as mean channel width and water temperature have affected these assemblages more strongly downstream from the reservoir than in the headwaters. The water-quality changes due to acidification, combined with the decreased flows and lowered water temperatures below the reservoir, have disrupted downstream continuum of macroinvertebrate communities that would normally be observed from the headwaters to the mouth. The information presented herein provides a basis for further evaluation of the Neversink and similar river systems, and for assessment of the effectiveness of future conservation efforts.

New York↗

Hydrology of and Current Monitoring Issues for the Chicago Area Waterway System, Northeastern Illinois

The Chicago Area Waterway System (CAWS) consists of a combination of natural and manmade channels that form an interconnected navigable waterway of approximately 90-plus miles in the metropolitan Chicago area of northeastern Illinois. The CAWS serves the area as the primary drainage feature, a waterway transportation corridor, and recreational waterbody. The CAWS was constructed by the Metropolitan Water Reclamation District of Greater Chicago (MWRDGC). Completion of the Chicago Sanitary and Ship Canal (initial portion of the CAWS) in 1900 breached a low drainage divide and resulted in a diversion of water from the Lake Michigan Basin. A U.S. Supreme Court decree (Consent Decree 388 U.S. 426 [1967] Modified 449 U.S. 48 [1980]) limits the annual diversion from Lake Michigan. While the State of Illinois is responsible for the diversion, the MWRDGC regulates and maintains water level and water quality within the CAWS by using several waterway control structures. The operation and control of water levels in the CAWS results in a very complex hydraulic setting characterized by highly unsteady flows. The complexity leads to unique gaging requirements and monitoring issues. This report provides a general discussion of the complex hydraulic setting within the CAWS and quantifies this information with examples of data collected at a range of flow conditions from U.S. Geological Survey streamflow gaging stations and other locations within the CAWS. Monitoring to address longstanding issues of waterway operation, as well as current (2014) emerging issues such as wastewater disinfection and the threat from aquatic invasive species, is included in the discussion.

Illinois↗

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

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

Idaho↗

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

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

Vermont↗

Watershed Boundary Dataset for Mississippi

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

Scientific Investigations Map↗

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

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

Arizona;Utah↗

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

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

New Mexico↗

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

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

Scientific Investigations Map↗

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

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

Scientific Investigations Report↗

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

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

Scientific Investigations Report↗

Thermal mapping of Hawaiian volcanoes with ASTER satellite data

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

Hawai'i↗