USGS Science⌕ Search

SEARCH · USGS Science

Results for “Computer contribution”

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 487 records · Page 27Linked to original sources

Estimated effects of pumping on groundwater storage and Walker River stream efficiencies in Smith and Mason Valleys, west-central Nevada

The Walker River originates in the Sierra Nevada Mountains and flows nearly 160 miles to its terminus at Walker Lake in west-central Nevada. The river provides a source of irrigation water for tens of thousands of acres of agricultural lands in California and Nevada and is the principal source of inflow to Walker Lake. Extraction of groundwater for agricultural use became prevalent in the late 1950s and early 1960s to supplement irrigation demands not met by surface-water diversions during times of drought. There is growing concern that continued groundwater withdrawals within the Walker River Basin are likely contributing to depleted streamflow of the Walker River and the long-term depletion of groundwater storage in the basin. This report documents changes in groundwater storage-volume and trends in Walker River stream efficiency, a measure of change in flow due to gaining or losing conditions, in the two largest agricultural valleys in the Walker River Basin, Smith and Mason Valleys, for a multi-decade period. Groundwater-level maps from previous studies were used for the beginning (1970) and middle (2006) points of this study. Groundwater levels measured from 1991–95 and 2016–20 were used to construct median groundwater-level maps that represented conditions in 1995 and 2020. Valley wide groundwater-level change was calculated by comparing groundwater-level maps for the periods 1970–95, 1996–2006, and 2007–20 and by observing the overall change from 1970 to 2020. Groundwater storage-volume change was calculated using groundwater-level change and previously defined specific yield values. Between 1970 and 2020, groundwater storage-volume declined 287,600 acre-feet in Smith Valley and 269,000 acre-feet in Mason Valley. Using groundwater storage-volume decline and annual groundwater pumpage rates, a maximum groundwater pumpage rate can be computed to support management of water resources. In Smith Valley, groundwater pumping in excess of 22,300 acre-feet per year would likely result in groundwater storage decline. In Mason Valley, groundwater pumping in excess of 75,200 acre-feet per year would likely result in groundwater storage decline. Stream efficiency was calculated using continuous streamflow data and monthly diversion volumes on two reaches: (1) the West Walker River in Smith Valley, from 1948 to 2020 and (2) the Walker River in Mason Valley, from 1958 to 2020. Stream efficiency during non-irrigation season in Smith and Mason Valleys declined at a statistically significant rate of 1.1 and 0.6 percent per year, respectively. Trends in stream efficiency corresponded to occurrence of prolonged drought, deviation from average annual streamflows, and total groundwater pumpage. Long-term declines in groundwater storage-volume and stream efficiency demonstrate that the alluvial aquifer system is becoming increasingly depleted, such that the river can no longer replenish groundwater storage while simultaneously balancing groundwater and surface-water withdrawals. The introduction of supplemental groundwater pumpage was intended to offset surface-water deficits during dry years; however, pumpage occurs even in years when average or above average streamflows meet surface-water demands. Reliance on supplemental groundwater pumpage has resulted in widespread groundwater storage-volume decline and decreased stream efficiency. With each successive drought cycle, the ability of Walker River to sustain streamflows and convey water downstream has diminished. Above average wet periods have a marginal and short-lived effect on rebounding the groundwater levels outside of the river corridor. Moreover, if the trend continues, each future drought cycle may further reduce groundwater supplies and that may further decrease streamflow reliability.

Nevada↗

Suspended-sediment and nutrient loads for Waiakea and Alenaio Streams, Hilo, Hawaii, 2003-2006

Suspended sediment and nutrient samples were collected during wet-weather conditions at three sites on two ephemeral streams in the vicinity of Hilo, Hawaii during March 2004 to March 2006. Two sites were sampled on Waiakea Stream at 80- and 860-foot altitudes during March 2004 to August 2005. One site was sampled on Alenaio Stream at 10-foot altitude during November 2005 to March 2006. The sites were selected to represent different land uses and land covers in the area. Most of the drainage area above the upper Waiakea Stream site is conservation land. The drainage areas above the lower site on Waiakea Stream, and the site on Alenaio Stream, are a combination of conservation land, agriculture, rural, and urban land uses. In addition to the sampling, continuous-record streamflow sites were established at the three sampling sites, as well as an additional site on Alenaio Stream at altitude of 75 feet and 0.47 miles upstream from the sampling site. Stage was measured continuously at 15-minute intervals at these sites. Discharge, for any particular instant, or for selected periods of time, were computed based on a stage-discharge relation determined from individual discharge measurements. Continuous records of discharge were computed at the two sites on Waiakea Stream and the upper site on Aleniao Stream. Due to non-ideal hydraulic conditions within the channel of Alenaio Stream, a continuous record of discharge was not computed at the lower site on Alenaio Stream where samples were taken. Samples were analyzed for suspended sediment, and the nutrients total nitrogen, dissolved nitrite plus nitrate, and total phosphorus. Concentration data were converted to instantaneous load values: loads are the product of discharge and concentration, and are presented as tons per day for suspended sediment or pounds per day for nutrients. Daily-mean loads were computed by estimating concentrations relative to discharge using graphical constituent loading analysis techniques. Daily-mean loads were computed at the two Waiakea Stream sampling sites for the analyzed constituents, during the period October 1, 2003 to September 30, 2005. No record of daily-mean load was computed for the Alenaio Stream sampling site due to the problems with computing a discharge record. The maximum daily-mean loads for the upper site on Waiakea Stream for suspended sediment was 79 tons per day, and the maximum daily-mean loads for total nitrogen, dissolved nitrite plus nitrate, and total phosphorus were 1,350, 13, and 300 pounds per day, respectively. The maximum daily-mean loads for the lower site on Waiakea Stream for suspended sediment was 468 tons per day, and the maximum daily-mean loads for total nitrogen, nitrite plus nitrate, and total phosphorus were 913, 8.5, and 176 pounds per day, respectively. From the estimated continuous daily-mean load record, all of the maximum daily-mean loads occurred during October 2003 and September 2004, except for suspended sediment load for the lower site, which occurred on September 15, 2005. Maximum values were not all caused by a single storm event. Overall, the record of daily-mean loads showed lower loads during storm events for suspended sediments and nutrients at the downstream site of Waiakea Stream during 2004 than at the upstream site. During 2005, however, the suspended sediment loads were higher at the downstream site than the upstream site. Construction of a flood control channel between the two sites in 2005 may have contributed to the change in relative suspended-sediment loads.

Open-File Report↗

Benthic invertebrates of fixed sites in the western Lake Michigan drainages, Wisconsin and Michigan, 1993-95

This report describes the variability in family-level benthic-invertebrate population data and the reliability of the data as a water-quality indicator for 11 fixed surface-water sites in the Western Lake Michigan Drainages study area of the National Water-Quality Assessment Program. Benthic-invertebrate-community measures were computed for the following: number of individuals, Hilsenhoff’s Family-Level Biotic Index, number and percent EPT (Ephemeroptera, Plecoptera, and Tricoptera), Margalef’s Diversity Index, and mean tolerance value. Relations between these measures and environmental setting, habitat, and of chemical water quality are examined. Benthic-invertebrate communities varied greatly among fixed sites and within individual streams among multiple-reach and multiple-year sampling. The variations between multiple reaches and years were sometimes larger than those found between different fixed sites. Factors affecting benthic invertebrates included both habitat and chemical quality. Generally, fixed-site streams with the highest diversity, greatest number of benthic invertebrates, and those at which community measures indicated the best water quality also had the best habitat and chemical quality. Variations among reaches are most likely related to differences in habitat. Variations among years are most likely related to climatic changes, which create variations in flow and/or chemical quality. The variability in the data analyzed in this study shows how benthic invertebrates are affected by differences in both habitat and water quality, making them useful indicators of stream health; however, a single benthic-invertebrate sample alone cannot be relied upon to accurately describe water quality of the streams in this study. Benthic-invertebrate data contributed valuable information on the biological health of the 11 fixed sites when used as one of several data sources for assessing water quality.

Michigan, Wisconsin↗

Synthesis of nutrient and sediment data for watersheds within the Chesapeake Bay drainage basin

Nutrient and sediment data collected by Federal and state agencies from 1972 through 1992 at 1,058 surface-water sites in nontidal parts of the Chesapeake Bay Basin were compiled into a large database. Adequate nutrient, sediment, and streamflow data were not available to compute annual loads for all sites because water-quality monitoring at many of the sites was either short term or noncontinuous or because stream-flow was not measured. Annual nutrient and sediment loads were calculated at a total of 127 sites. Annual loads of dissolved nitrate were calculated for 108 sites, but total nitrogen loads could be calculated for only 48 of these sites because ammonia plus organic nitrogen data were not available for many of these 108 sites. Annual loads of total phosphorus were calculated for 99 sites, and annual loads of suspended sediment were calculated for 33 sites. Loads could be calculated for only a very few sites in the Juniata River Basin (a tributary to the Susquehanna River), the York River Basin, the middle and lower reaches of the James River, and the nontidal parts of the eastern shore of the Bay. Geographic Information System (GIS) spatial data sets of land use, physiographic province, rock type, and watershed delineation were compiled for the entire Chesapeake Bay Basin (approximately 64,000 square miles). The nutrient- and sediment-yield were evaluated with respect to land use, physiographic province, rock type, and hydrologic characteristics. During years that the mean streamflow was about equal to the long-term mean streamflow, the Susquehanna River contributed about 50 percent of the freshwater, 66 percent of the total nitrogen, and 40 percent of the total phosphorus transported by tributaries to the Bay. Nutrient and sediment data were available for less than 18 percent of the predominantly agricultural areas underlain by siliciclastic rock and for less than 35 percent of the predominantly agricultural areas underlain by either carbonate rock or unconsolidated rock. Nutrient and sediment data were available for about 91 percent of the predominantly forested areas underlain by siliciclastic rock. Spatial and temporal gaps in the water-quality data and GIs data sets limited some data analysis. Correlations of annual yields or nutrients and sediment with respect to land use, physiographic province, and rock type indicated (1) basins with larger percentages of agricultural land had larger nutrient and sediment yields, (2) basins with larger percentages of forest land had smaller nutrient and sediment loads, (3) the largest total nitrogen yields were from agricultural basins underlain by carbonate rock, (4) yields or nutrients from urban basins were substantially less than yields from agricultural basins, and (5) basins with small amounts of agricultural and urban land had relatively small nutrient and sediment yields.

Chesapeake Bay drainage basin↗

Genesis and evolution of ferromanganese crusts from the summit of Rio Grande Rise, southwest Atlantic Ocean

The Rio Grande Rise (RGR) is a large elevation in the Atlantic Ocean and known to host potential mineral resources of ferromanganese crusts (Fe–Mn), but no investigation into their general characteristics have been made in detail. Here, we investigate the chemical and mineralogical composition, growth rates and ages of initiation, and phosphatization of relatively shallow-water (650–825 m) Fe–Mn crusts dredged from the summit of RGR by using computed tomography, X-ray diffraction, 87 Sr/ 86 Sr ratios, U–Th isotopes, and various analytical techniques to determine their chemical composition. Fe–Mn crusts from RGR have two distinct generations. The older one has an estimated age of initiation around 48–55 Ma and was extensively affected by post-depositional processes under suboxic conditions resulting in phosphatization during the Miocene (from 20 to 6.8 Ma). As a result, the older generation shows characteristics of diagenetic Fe–Mn deposits, such as low Fe/Mn ratios (mean 0.52), high Mn, Ni, and Li contents and the presence of a 10 Å phyllomanganate, combined with the highest P content among crusts (up to 7.7 wt %). The younger generation is typical of hydrogenetic crusts formed under oxic conditions, with a mean Fe/Mn ratio of 0.75 and mean Co content of 0.66 wt %, and has the highest mean contents of Bi, Nb, Ni, Te, Rh, Ru, and Pt among crusts formed elsewhere. The regeneration of nutrients from local biological productivity in the water column is the main source of metals to crusts, providing mainly metals that regenerate rapidly in the water column and are made available at relatively shallow water depths (Ni, As, V, and Cd), at the expense of metals of slower regeneration (Si and Cu). Additionally, important contributions of nutrients may derive from various water masses, especially the South Atlantic Mode Water and Antarctic Intermediate Water (AAIW). Bulk Fe–Mn crusts from the summit of RGR plateau are generally depleted in metals considered of greatest economic interest in crusts like Co, REE, Mo, Te, and Zr, but are the most enriched in the critical metals Ni and Li compared to other crusts. Further investigations are warranted on Fe–Mn crusts from deeper-water depths along the RGR plateau and surrounding areas, which would less likely be affected by phosphatization.

Rio Grande Rise↗

Discharge data collection and analysis and implications for surface-water/groundwater interactions in the lower Las Vegas Wash, Clark County, Nevada, 2016–18

The lower Las Vegas Wash represents the terminal surface drainage for the Las Vegas Valley in southern Nevada. In 1997, high concentrations of perchlorate were found in seeps contributing to discharge in this area and traced to an industrial byproduct from manufacturing operations in the mid-1900s at the nearby Basic Magnesium, Incorporated, plant. The discovery prompted a water-resources investigation by the Nevada Department of Environmental Protection (NDEP) to develop an understanding of the nearby groundwater flow system and the dynamics associated with surface-water flow in the Wash. In 2016, the U.S. Geological Survey was tasked with evaluating surface-water discharge in the lower Las Vegas Wash near locations where perchlorate concentrations from the groundwater system had been detected. Results of this study will assist NDEP with identifying areas of groundwater and surface-water interaction and help guide future cleanup and monitoring efforts. Streamflow discharge is evaluated along a 4-mile section of the lower Las Vegas Wash (referred to as the Wash) and used to describe surface-water and groundwater interactions between the Wash channel and bank sediments. Continuous discharge data were collected during a 2-year period (2016–18) at 5 gaging stations along the Wash. Additionally, multiple discrete measurements between gaging stations were collected during 4 synoptic sampling events between 2016 and 2018. A diurnal discharge pattern, controlled by upstream treated wastewater releases, provided high- and low-discharge markers that are used to compute downstream time-lags of peak and minimum flows. Computed time-lags are used to establish travel times between measurement sites, and difference in upstream and time-lagged downstream hydrographs are used to compute increases (gain) or decreases (loss) in discharge between gaging stations or between gaging stations and discrete measurements. Tributary surface-water inflows to the lower Las Vegas Wash from wastewater discharge, remediation efforts, and periodic flooding from rainfall runoff are included in computing differences in discharge. Differences between discharge data from delineated reaches are used to define locations of daily, monthly, and yearly streamflow gains from or losses to adjacent bank sediments. Construction of additional channel-stabilization weirs have occurred since the completion of this study and the associated change to streamflow dynamics may limit study results to the period analyzed; however, methods and processes described in this report can be used in future evaluations.

Nevada↗

Assessing streamflow sensitivity to variations in glacier mass balance

The mountains ringing the Gulf of Alaska (GOA) receive upwards of 4–8 m yr −1 of precipitation (Simpson et al.2005; Weingartner et al. 2005; O’Neel 2012), much of which runs off into productive coastal waters. The alpine landscape is heavily glacierized, and storage and turnover of water by glaciers substantially influences the regional surface water balance (Neal et al. 2010). In turn, the land-to-ocean flux of freshwater impacts the biogeochemistry, physical oceanography, freshwater and marine ecology of the downstream components of the GOA ecosystem (e.g., Royer et al. 2001; Hood and Scott 2008). In this way, the links between terrestrial and ocean ecosystems along the GOA have widespread impacts on regional socioeconomic issues including water and hydropower resources, fish populations, and sea level change (Dorava and Milner 2000; Royer and Grosch 2006; Cherry et al. 2010; Gardner et al. 2013). Moreover, predicting future changes in physical, chemical and biological processes in near-shore ecosystems along the GOA hinges, in part, on developing a robust understanding of water storage and transfer by glaciers through streams to the ocean. Glacierized basins (i.e. presently ice covered as opposed to glaciated, or historically ice covered) are very efficient producers of runoff, yielding 2–10 times greater runoff than similarly sized, non-glacierized basins (Mayo 1984). The unique energy balance that characterizes these basins (Jansson et al. 2003; Hock 2005) results in substantial alterations to streamflow, even when fractional ice coverage is very small (Stahl and Moore 2006). Consistent and precise treatment of glacier runoff is key to accurate assessment of hydrologic, ecological and socioeconomic impacts, but previously used definitions for glacier runoff are variable. They include: 1) meltwater produced as negative annual mass balance (e.g., Fountain and Tangborn 1985); 2) storage changes in the monthly water budget, where solid precipitation is balanced by melt and evaporation (Huss 2011, concept #2); 3) meltwater derived from melting ice only (irrespective of melting snow or mass balance) (Nolin et al. 2010; Huss 2011, concept #1); 4) all meltwater derived from the glacier surface (Cogley et al. 2011, meltwater runoff); 5) total runoff from the glacier surface (meltwater runoff plus rain on the glacier) (Neal et al. 2010). Total glacier runoff (Definitions 4 and 5 above) includes a contribution from annual mass balance, i.e. the sum of accumulation and ablation through a mass balance year (Definition 1), or what has historically been referred to as the “net” balance (Cogley et al. 2011). Indeed, annual balance has been shown to be an important driver of streamflow trends in glacierized basins, with periods of persistent negative annual balance resulting in statistically significant increases in streamflow (e.g., Pellicciotti et al. 2010). However, in maritime climates, anomalies in glacier runoff can be disconnected from annual balance because of the high variability in winter precipitation. For example, positive anomalies in winter accumulation can result in elevated levels of glacier runoff in times of positive annual mass balance (Thayyen and Gergan 2010). Quantifying the impacts of changing glacier geometries (annual balance) on glacier runoff is essential for predicting future changes in streamflow in glacierized basins. However, determining the role that this component plays in total glacier runoff (Definition 5) requires consistent measurements of seasonal (or shorter period) mass balances, measurements of precipitation at multiple locations within a basin, and streamflow measurements in close proximity to a glacier’s terminus. Practical and logistical challenges associated with assembling such data sets typically preclude such partitioning. As a result, most analyses of the relationship between annual mass balance and streamflow rely on some component of model output to compute glacier runoff (e.g. Huss et al. 2008; Kaser et al. 2010). Ultimately, developing an understanding of how total glacier runoff will change in the future is critical for predicting downstream ecological impacts associated with changes in riverine fluxes of water, sediment, and solutes (e.g., metals and nutrients) to near-shore coastal ecosystems. The purpose of this paper is to evaluate relationships among seasonal and annual glacier mass balances, glacier runoff and streamflow in two glacierized basins in different climate settings. We use long-term glacier mass balance and streamflow datasets from the United States Geological Survey (USGS) Alaska Benchmark Glacier Program to compare and contrast glacier-streamflow interactions in a maritime climate (Wolverine Glacier) with those in a continental climate (Gulkana Glacier). Our overall goal is to improve our understanding of how glacier mass balance processes impact streamflow, ultimately improving our conceptual understanding of the future evolution of glacier runoff in continental and maritime climates.

Alaska↗

Directional connectivity in hydrology and ecology

Quantifying hydrologic and ecological connectivity has contributed to understanding transport and dispersal processes and assessing ecosystem degradation or restoration potential. However, there has been little synthesis across disciplines. The growing field of ecohydrology and recent recognition that loss of hydrologic connectivity is leading to a global decline in biodiversity underscore the need for a unified connectivity concept. One outstanding need is a way to quantify directional connectivity that is consistent, robust to variations in sampling, and transferable across scales or environmental settings. Understanding connectivity in a particular direction (e.g., streamwise, along or across gradient, between sources and sinks, along cardinal directions) provides critical information for predicting contaminant transport, planning conservation corridor design, and understanding how landscapes or hydroscapes respond to directional forces like wind or water flow. Here we synthesize progress on quantifying connectivity and develop a new strategy for evaluating directional connectivity that benefits from use of graph theory in ecology and percolation theory in hydrology. The directional connectivity index (DCI) is a graph-theory based, multiscale metric that is generalizable to a range of different structural and functional connectivity applications. It exhibits minimal sensitivity to image rotation or resolution within a given range and responds intuitively to progressive, unidirectional change. Further, it is linearly related to the integral connectivity scale length—a metric common in hydrology that correlates well with actual fluxes—but is less computationally challenging and more readily comparable across different landscapes. Connectivity-orientation curves (i.e., directional connectivity computed over a range of headings) provide a quantitative, information-dense representation of environmental structure that can be used for comparison or detection of subtle differences in the physical-biological feedbacks driving pattern formation. Case-study application of the DCI to the Everglades in south Florida revealed that loss of directional hydrologic connectivity occurs more rapidly and is a more sensitive indicator of declining ecosystem function than other metrics (e.g., habitat area) used previously. Here and elsewhere, directional connectivity can provide insight into landscape drivers and processes, act as an early-warning indicator of environmental degradation, and serve as a planning tool or performance measure for conservation and restoration efforts.

Ecological Applications↗

Geohydrology and evaluation of water-resource potential of the Upper Floridan aquifer in the Albany area, southwestern Georgia

In the Albany area of southwestern Georgia, the Upper Floridan aquifer lies entirely within the Dougherty Plain district of the Coastal Plain physiographic province, and consists of the Ocala Limestone of late Eocene age. The aquifer is divided throughout most of the study area into an upper and a lower lithologic unit, which creates an upper and a lower water-bearing zone. The lower water-bearing zone consists of alternating layers of sandy limestone and medium-brown, recrystallized dolomitic limestone, and ranges in thickness from about 50 to 100 feet. It is highly fractured, and exhibits well-developed permeability by solution features that are responsible for transmitting most of the ground water in the aquifer. Transmissivity of the lower water-bearing zone ranges from about 90,000 to 178,000 feet squared per day. The upper water-bearing zone is a finely crystallized-to-oolitic, locally dolomitic limestone having an average thickness of about 60 feet. Transmissivities in the upper water-bearing zone are considerably less than those in the lower water-bearing zone. The Upper Floridan aquifer is overlain by about 20 to 120 feet of undifferentiated overburden consisting of fine-to-coarse quartz sand and noncalcareous clay. A clay zone about 10 to 30 feet thick may be continuous throughout the southwestern part of the Albany area, and where present, causes confinement of the Upper Floridan aquifer and creates perched ground water after periods of heavy rainfall. The Upper Floridan aquifer is confined below by the Lisbon Formation, a mostly dolomitic limestone that contains trace amounts of glauconite. The Lisbon Formation is at least 50 feet thick in the study area, and acts as an impermeable base to the Upper Floridan aquifer. The quality of ground-water in the Upper Floridan aquifer is suitable for most uses; wells generally yield water of the hard, calcium-bicarbonate type that generally meets the U.S. Environmental Protection Agency's Primary or Secondary Drinking Water Regulations. The water-resource potential of the Upper Floridan aquifer was evaluated by compiling results of test drilling and aquifer testing in the study area, and by conducting computer simulations of the ground-water-flow system under the seasonal-low conditions of November 1985, and under conditions of pumping within a 12square-mile area located southwest of Albany. Results of test drilling, aquifer testing, and water-quality analyses indicate that, in the area southwest of Albany, geohydrologic conditions in the Upper Floridan aquifer, undifferentiated overburden, and Lisbon Formation were favorable for the aquifer to provide a large quantity of water without having adverse effects on the ground-water system. The confinement of the Upper Floridan aquifer by the undifferentiated overburden and the rural setting of the area of potential development decreases the likelihood that chemical constituents will enter the aquifer during development of the ground-water resources. Computer simulations of ground-water flow in the Upper Floridan aquifer, incorporating conditions for regional flow across model boundaries, leakage from rivers and other surface-water features, and vertical leakage from the undifferentiated overburden, were conducted by using a finite-element model for groundwater flow in two dimensions. Comparison of computed and measured water levels in the Upper Floridan aquifer for November 1985 at 74 locations indicated that computed water levels generally were within 5 feet of the measured values, which is the accuracy to which measured water levels were known. Water-level altitudes ranged from about 260 feet to 130 feet above sea level in the study area during calibration. Aquifer discharge to the Flint River downstream from the Lake Worth dam was computed by the calibrated model to be about 1 billion gallons per day; about 300 million gallons per day greater than was measured for similar low-flow conditions. The excess computed discharge was attributed partially to stream withdrawals for industrial use, non-reported use, and channel evaporation, but mostly to increased gradients and increased flow from the aquifer to the river than existed during calibration. Results from the calibrated finite-element model indicate that ground-water flow is dominated by inflow from regional-flow components to the west, north, and east of the study area, and by outflow to the Flint River downstream from the Lake Worth dam. Simulation results indicated that directions of ground-water flow were not changed appreciably by pumping at the November 1985 rates. However, vertical leakage from the undifferentiated overburden caused local deviations in the regional flow pattern. A sensitivity analysis that was performed on 18 hydrologic factors affecting the flow system in the Upper Floridan aquifer showed that computed water levels changed the most (were the most sensitive) in response to changes in hydraulic conductivity of the aquifer, vertical leakage coefficient and water level in the undifferentiated overburden, and stage of the Flint River downstream from the Lake Worth dam. Computed water levels were least sensitive to changes in well pumpage, flow across the northern boundary and from Lake Worth, the boundary coefficient for the Flint River downstream from the Lake Worth dam, and flow from Cooleewahee Creek. Simulations of six pumping scenarios in the area of potential development southwest of Albany showed that the Upper Floridan aquifer is capable of providing at least 72 million gallons per day from five locations (14.4 million gallons per day each) within this area without causing adverse affects on the flow system. The 72million-gallon-per-day scenario yielded a maximum drawdown of about 9.4 feet, which placed the water level in the Upper Floridan aquifer about 50 feet above the top of the lower water-bearing zone. Hence, the likelihood of aquifer dewatering, well interference, or sinkhole development from pumping as much as 72 million gallons per day from within the area of potential development is small. All pumping scenarios showed that about 81 percent of the ground-water pumpage was derived from regional flow that would have discharged to the Flint River downstream from the Lake Worth dam. The dominant ground-water-flow direction toward the Flint River was not changed and no induced recharge from the Flint River entered the potential-development area. Induced recharge from the undifferentiated overburden contributed to about 1.5 percent of the total volume pumped during the simulations.

Georgia↗

Bathymetric and velocimetric surveys at highway bridges crossing the Missouri and Mississippi Rivers near St. Louis, Missouri, August 3–10, 2020

Bathymetric and velocimetric data were collected by the U.S. Geological Survey, in cooperation with the Missouri Department of Transportation, near 15 bridges at 10 highway crossings of the Missouri and Mississippi Rivers near Washington, Louisiana, and St. Louis, Missouri, on August 3–10, 2020. A multibeam echosounder mapping system was used to obtain channel-bed elevations for river reaches about 1,640 to 1,970 feet longitudinally and generally extending laterally across the active channel from bank to bank during moderate flood-flow conditions. These surveys provided channel geometry and hydraulic conditions at the time of the surveys and provided characteristics of scour holes that may be useful in developing predictive guidelines or equations for computing potential scour depth. These data also may be useful to the Missouri Department of Transportation as a low to moderate flood-flow assessment of the bridges for stability and integrity issues with respect to bridge scour during floods. Bathymetric data were collected around every in-channel pier. Scour holes were present at most piers for which bathymetry could be obtained, except those on banks or surrounded by riprap. All the bridge sites in this study were previously surveyed and documented in previous studies, including the two new bridge structures at Louisiana and Washington (structures A8141 and A8504, sites 22 and 32, respectively). Comparisons between bathymetric surfaces from the previous surveys and those of the current (2020) study do not indicate any consistent correlation in channel-bed elevations with streamflow conditions. The comparisons of the 2020 surveys to two previous surveys at the new bridge structure A8141 at Washington (site 22) resulted in net erosion of the channel bed in both comparisons, despite the 2020 streamflow being less than either previous survey. Alternatively, there was a net gain of sediment at new bridge structure A8504 at Louisiana (site 32) between 2014 and 2020, which was the most substantial increase in the surveys detailed in this report; substantially less flow in 2020 than in 2014 or changes to the channel and spur dikes near the bridge may have contributed to the observed sediment gain. Pier size, nose shape, and skew to approach flow had a substantial effect on the size of the scour hole observed at a given pier. Larger and deeper scour holes were present at piers with wide or blunt noses caused by exposed footings, seal courses, or caissons. When a pier was skewed to primary approach flow, the scour hole was generally deeper and larger than at a similar pier without skew; however, the shape of the scour hole near skewed piers in this study generally was longer and deeper on the leeward side, contrary to the general shape of scour holes for skewed piers. However, this phenomenon has been observed historically at these sites, and likely is exacerbated by debris rafts or other turbulence-inducing features near the atypical scour holes. A substantial scour hole was observed near pier 11 of structure A6500 (site 33), which was deeper than in the 2016 survey. The scour holes observed at pier 17 of structure L0561 (site 25) and piers 3 and 4 of structure A1500 (site 34) also were slightly deeper and wider in 2020 than in 2016. At new bridge structures A8141 at Washington (site 22) and A8504 at Louisiana (site 32), the smaller cross-sectional area and configuration of the piers of the new bridges resulted in substantially less scour than with the wider old piers.

Illinois, Missouri↗

Simulated Effects of Year 2030 Water-Use and Land-Use Changes on Streamflow near the Interstate-495 Corridor, Assabet and Upper Charles River Basins, Eastern Massachusetts

Continued population growth and land development for commercial, industrial, and residential uses have created concerns regarding the future supply of potable water and the quantity of ground water discharging to streams in the area of Interstate 495 in eastern Massachusetts. Two ground-water models developed in 2002-2004 for the Assabet and Upper Charles River Basins were used to simulate water supply and land-use scenarios relevant for the entire Interstate-495 corridor. Future population growth, water demands, and commercial and residential growth were projected for year 2030 by the Metropolitan Area Planning Council. To assess the effects of future development on subbasin streamflows, seven scenarios were simulated by using existing computer-based ground-water-flow models with the data projected for year 2030. The scenarios incorporate three categories of projected 2030 water- and land-use data: (1) 2030 water use, (2) 2030 land use, and (3) a combination of 2030 water use and 2030 land use. Hydrologic, land-use, and water-use data from 1997 through 2001 for the Assabet River Basin study and 1989 through 1998 for the Upper Charles River Basin study were used to represent current conditions - referred to as 'basecase' conditions - in each basin to which each 2030 scenario was compared. The effects of projected 2030 land- and water-use change on streamflows in the Assabet River Basin depended upon the time of year, the hydrologic position of the subbasin in the larger basin, and the relative areas of new commercial and residential development projected for a subbasin. Effects of water use and land use on streamflow were evaluated by comparing average monthly nonstorm streamflow (base flow) for March and September simulated by using the models. The greatest decreases in streamflow (up to 76 percent in one subbasin), compared to the basecase, occurred in September, when streamflows are naturally at their lowest level. By contrast, simulated March streamflows decreased less than 6.5 percent from basecase streamflows in all subbasins for all scenarios. The simulations showed similar effects in the Upper Charles River Basin, but increased water use contributed to decreased simulated streamflow in most subbasins. Simulated changes in March streamflows for 2030 in the Upper Charles River Basin were within +- 6 percent of the basecase for all scenarios and subbasins. Percentage decreases in simulated September streamflows for 2030 were greater than in March but less than the September decreases that resulted for some subbasins in the Assabet River Basin. Only two subbasins of the Upper Charles River Basin had projected decreases greater than 5 percent. In the Mill River subbasin, the decrease was 11 percent, and in the Mine Brook subbasin, 6.6 percent. Changes in water use and wastewater return flow generally were found to have the greatest effect in the summer months when streamflow and aquifer recharge rates are low and water use is high. September increases in main-stem streamflow of both basins were due mainly to increased discharge of treated effluent from wastewater-treatment facilities on the main-stem rivers. In the Assabet River Basin, wastewater-treatment-facility discharge became a smaller proportion of total streamflow with distance downstream. In contrast, wastewater-treatment facility discharge in the Upper Charles River Basin became a greater proportion of streamflow with distance downstream. The effects of sewer-line extension and low-impact development on streamflows in two different subbasins of the Assabet River Basin also were simulated. The result of extending sewer lines with a corresponding decrease in septic-system return flow caused September streamflows to decrease as much as 15 percent in the Fort Pond Brook subbasin. The effect of low-impact development was simulated in the Hop Brook subbasin in areas projected for commercial development. In this simulation, the greater the area where low-i

Scientific Investigations Report↗

Documentation of linear regression models for computing water-quality constituent concentrations using continuous real-time water-quality data for the North Fork Ninnescah River and Cheney Reservoir, Kansas, 2014–21

Cheney Reservoir, in south-central Kansas, was constructed to provide a reliable municipal water supply for the city of Wichita, Kansas, and to provide downstream flood control, wildlife habitat, and recreation. Cheney Reservoir will continue to be important for municipal water supply use as needs increase with ongoing population growth and urban development. Advanced notification of changing water-quality conditions near water-treatment facility intakes and in source waters allows water-treatment facilities and resource planning officials to proactively monitor changing conditions. The U.S. Geological Survey (USGS), in cooperation with the City of Wichita, collected water-quality data at the North Fork Ninnescah River above Cheney Reservoir (USGS station 07144780) and Cheney Reservoir near Cheney, Kans. (USGS station 07144790), monitoring sites to update and develop regression models relating continuous water-quality constituents, streamflow, reservoir storage, and seasonal components to discretely sampled water-quality constituent concentrations of interest. Linear regression analysis was used to update and develop models for alkalinity, major ions, nutrients (nitrogen and phosphorus species), total and dissolved organic carbon, total suspended solids, suspended sediment, fecal indicator bacteria, and atrazine at the North Fork Ninnescah River site and total dissolved solids, major ions, hardness as calcium carbonate, nutrients (nitrogen and phosphorus species), chlorophyll a , and suspended sediment at the Cheney Reservoir site. New and updated models for both sites are applicable to the period of YSI EXO water-quality monitor and sensor deployment (November 14, 2015, through September 30, 2021, at the North Fork Ninnescah River site; October 1, 2014, through September 30, 2021, at the Cheney Reservoir site). Models and resulting water-quality information included in this report can be used in real time, potentially as guidance for water-treatment processes, and can be used to characterize changes in water-quality conditions over time in Cheney Reservoir and its contributing drainage basin provided that the deployed equipment, sensors, and location do not change.

Kansas↗

Coastal groundwater dynamics off Santa Barbara, California: combining geochemical tracers, electromagnetic seepmeters, and electrical resistivity

This paper presents repeat field measurements of 222 Rn and 223,224,226,228 Ra, electromagnetic seepage meter-derived advective fluxes, and multi-electrode, stationary and continuous marine resistivity surveys collected between November 2005 and April 2007 to study coastal groundwater dynamics within a marine beach in Santa Barbara, California. The study provides insight into magnitude and dynamics of submarine groundwater discharge (SGD) and associated nutrient loadings into near-shore coastal waters, where the predominant SGD drivers can be both spatially and temporally separated. Rn-222 and 223,224,226,228 Ra were utilized to quantify the total and saline contribution, respectively, of SGD. The two short-lived 224,223 Ra isotopes provided an estimate of apparent near-shore water mass age, as well as an estimate of the Ra-derived eddy diffusion coefficient, K h ( 224 Ra = 2.86 ± 0.7 m 2 s −1 ; 223 Ra = 1.32 ± 0.5 m 2 s −1 ). Because 222 Rn ( t ½ = 3.8 day) and 224 Ra ( t ½ = 3.66 day) have comparable half-lives and production terms, they were used in concert to examine respective water column removal rates. Electromagnetic seepage meters recorded the physical, bi-directional exchange across the sediment/water interface, which ranged from −6.7 to 14.5 cm day −1 , depending on the sampling period and position relative to the low tide line. Multi-day time-series 222 Rn measurements in the near-shore water column yielded total (saline + fresh) SGD rates that ranged from 3.1 ± 2.6 to 9.2 ± 0.8 cm day −1 , depending on the sampling season. Offshore 226 Ra ( t ½ = 1600 year) and 222 Rn gradients were used with the calculated K h values to determine seabed flux estimates (dpm m −2 day −1 ), which were then converted into SGD rates (7.1 and 7.9 cm day −1 , respectively). Lastly, SGD rates were used to calculate associated nutrient loads for the near-shore coastal waters off Santa Barbara. Depending on both the season and the SGD method utilized, the following SGD-derived nutrient inputs were computed (mol per day per meter of shoreline): NH 4 + = 0.06–0.29 mol day −1 m −1 ; SiO 4 = 0.22–0.29 mol day −1 m −1 ; PO 4 3− = 0.04–0.17 mol day −1 m −1 ; [NO 2 − + NO 3 − ] = 0–0.52 mol day −1 m −1 ; dissolved inorganic nitrogen (DIN) = 0.01–0.17 mol day −1 m −1 , and dissolved organic nitrogen (DON) = 0.08–0.09 mol day −1 m −1 . Compared to the ephemeral nature of fluvial and marine inputs into this region, such SGD-derived loadings can provide a sustained source of select nutrients to the coastal waters off Santa Barbara, California that should be accounted for in mass balance estimates.

California↗

Characterization of stormwater runoff from the Naval Air Station and Naval Wepons Industrial Reserve Plant, Dallas, Texas, 1994-96

The characterization of stormwater runoff from the Naval Air Station (NAS) and the Naval Weapons Industrial Reserve Plant (NWIRP), Dallas, Texas, is necessary to determine if runoff from the facilities is contributing to off-site contamination of surface waters, A network of five fixed sites and four grab sites was established to collect stormwater-runoff samples from a substantial part of the drainage area of each facility. Fixed sites were instrumented to measure and store precipitation, stage, discharge, and runoff-volume data and to collect flow-weighted composite samples during a storm. Grab and composite samples were collected for six storms at each of the five fixed sites from October 1994 to March 1996. The grab samples were analyzed for about 100 properties and constituents including specific conductance, pH, water temperature, bacteria, trace elements, oil and grease, total phenols, and volatile organic compounds. The composite samples were analyzed for about 220 properties and constituents including specific conductance, pH, chemical oxygen demand, biochemical oxygen demand, major ions, suspended and dissolved solids, nutrients, trace elements, total organic carbon, volatile organic compounds, semivolatile organic compounds, and organochlorine and organophosphorus pesticides. Grab samples were collected for two storms (September 18,1995, and October 2,1995) at each of the four grab sites. The grab samples were analyzed for about 80 constituents including specific conductance, pH, water temperature, trace elements, and volatile organic compounds. Composite samples were collected for two of the six storms sampled at the fixed sites and analyzed for aquatic toxicity. Fathead minnow growth and survival toxicity tests and water flea reproduction and survival toxicity tests were done. Median event-mean concentrations computed for 12 selected constituents in samples from NAS and NWIRP fixed sites were compared to median event-mean concentrations for residential, commercial, industrial, and highway land uses within the Dallas-Fort Worth area computed from data collected for the National Pollutant Discharge Elimination System program. NAS and NWIRP median event-mean concentrations also were compared to those for residential and commercial land uses from the Nationwide Urban Runoff Program.

Open-File Report↗

Volcano surveillance by ACR silver fox

Recent growth in the business of unmanned air vehicles (UAVs) both in the US and abroad has improved their overall capability, resulting in a reduction in cost, greater reliability and adoption into areas where they had previously not been considered. Uses in coastal and border patrol, forestry and agriculture have recently been evaluated in an effort to expand the observed area and reduce surveillance and reconnaissance costs for information gathering. The scientific community has both contributed and benefited greatly in this development. A larger suite of light-weight miniaturized sensors now exists for a range of applications which in turn has led to an increase in the gathering of information from these autonomous vehicles. In October 2004 the first eruption of Mount St Helens since 1986 caused tremendous interest amoUg people worldwide. Volcanologists at the U.S. Geological Survey rapidly ramped up the level of monitoring using a variety of ground-based sensors deployed in the crater and on the flanks of the volcano using manned helicopters. In order to develop additional unmanned sensing methods that can be used in potentially hazardous and low visibility conditions, a UAV experiment was conducted during the ongoing eruption early in November. The Silver Fox UAV was flown over and inside the crater to perform routine observation and data gathering, thereby demonstrating a technology that could reduce physical risk to scientists and other field operatives. It was demonstrated that UAVs can be flown autonomously at an active volcano and can deliver real time data to a remote location. Although still relatively limited in extent, these initial flights provided information on volcanic activity and thermal conditions within the crater and at the new (2004) lava dome. The flights demonstrated that readily available visual and infrared video sensors mounted in a small and relatively low-cost aerial platform can provide useful data on volcanic phenomena. This was made possible by utilizing GPS and computer-controlled flight direction and stabilization to acquire and track target areas within the Mount St. Helens crater. It was also determined that additional light-weight sensor development will be needed to enable autonomous measurements of volcanic gasses and imaging in poor-weather conditions. Copyright ?? 2005 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Conference Paper↗

Development of a comprehensive watershed model applied to study stream yield under drought conditions

We developed a model code to simulate a watershed's hydrology and the hydraulic response of an interconnected stream-aquifer system, and applied the model code to the Lower Republican River Basin in Kansas. The model code links two well-known computer programs: MODFLOW (modular 3-D flow model), which simulates ground water flow and stream-aquifer interaction; and SWAT (soil water assessment tool), a soil water budget simulator for an agricultural watershed. SWAT represents a basin as a collection of subbasins in terms of soil, land use, and weather data, and simulates each subbasin on a daily basis to determine runoff, percolation, evaporation, irrigation, pond seepage, and crop growth. Because SWAT applies a lumped hydrologic model to each sub-basin, spatial heterogeneities with respect to factors such as soil type and land use are not resolved geographically, but can instead be represented statistically. For the Republican River Basin model, each combination of six soil types and three land uses, referred to as a hydrologic response unit (HRU), was simulated with a separate execution of SWAT. A spatially weighted average was then taken over these results for each hydrologic flux and time step by a separate program, SWBAVG. We wrote a package for MODFLOW to associate each subbasin with a subset of aquifer grid cells and stream reaches, and to distribute the hydrologic fluxes given for each subbasin by SWAT and SWBAVG over MODFLOW's stream-aquifer grid to represent tributary flow, surface and ground water diversions, ground water recharge, and evapotranspiration from ground water. The Lower Republican River Basin model was calibrated with respect to measured ground water levels, streamflow, and reported irrigation water use. The model was used to examine the relative contributions of stream yield components and the impact on stream yield and base flow of administrative measures to restrict irrigation water use during droughts. Model results indicate that tributary flow is the dominant component of stream yield and that reduction of irrigation water use produces a corresponding increase in base flow and stream yield. However, the increase in stream yield resulting from reduced water use does not appear to be of sufficient magnitude to restore minimum desirable streamflows.

Groundwater↗

Comparison of groundwater flow in Southern California coastal aquifers

Development of the coastal aquifer systems of Southern California has resulted in overdraft, changes in streamflow, seawater intrusion, land subsidence, increased vertical flow between aquifers, and a redirection of regional flow toward pumping centers. These water-management challenges can be more effectively addressed by incorporating new understanding of the geologic, hydrologic, and geochemical setting of these aquifers. Groundwater and surface-water flow are controlled, in part, by the geologic setting. The physiographic province and related tectonic fabric control the relation between the direction of geomorphic features and the flow of water. Geologic structures such as faults and folding also control the direction of flow and connectivity of groundwater flow. The layering of sediments and their structural association can also influence pathways of groundwater flow and seawater intrusion. Submarine canyons control the shortest potential flow paths that can result in seawater intrusion. The location and extent of offshore outcrops can also affect the flow of groundwater and the potential for seawater intrusion and land subsidence in coastal aquifer systems. As coastal aquifer systems are developed, the source and movement of ground-water and surface-water resources change. In particular, groundwater flow is affected by the relative contributions of different types of inflows and outflows, such as pump-age from multi-aquifer wells within basal or upper coarse-grained units, streamflow infiltration, and artificial recharge. These natural and anthropogenic inflows and outflows represent the supply and demand components of the water budgets of ground-water within coastal watersheds. They are all significantly controlled by climate variability related to major climate cycles, such as the El Niño–Southern Oscillation and the Pacific Decadal Oscillation. The combination of natural forcings and anthropogenic stresses redirects the flow of groundwater and either mitigates or exacerbates the potential adverse effects of resource development, such as declining water levels, sea-water intrusion, land subsidence, and mixing of different waters. Streamflow also has been affected by development of coastal aquifer systems and related conjunctive use. Saline water is the largest water-quality problem in Southern California coastal aquifer systems. Seawater intrusion is a significant source of saline water, but saline water is also known to come from other sources and processes. Seawater intrusion is typically restricted to the coarse-grained units at the base of fining-upward sequences of terrestrial deposits, and at the top of coarsening upward sequences of marine deposits. This results in layered and narrow intrusion fronts. Maintaining the sustainability of Southern California coastal aquifers requires joint management of surface water and groundwater (conjunctive use). This requires new data collection and analyses (including research drilling, modern geohydrologic investigations, and development of detailed computer groundwater models that simulate the supply and demand components separately), implementation of new facilities (including spreading and injection facilities for artificial recharge), and establishment of new institutions and policies that help to sustain the water resources and better manage regional development.

California↗

Comparison of two methods for estimating discharge and nutrient loads from tidally affected reaches of the Myakka and Peace Rivers, west-central Florida

The Myakka and Peace River Basins constitute more than 60 percent of the total inflow area and contribute more than half the total tributary inflow to the Charlotte Harbor estuarine system. Water discharge and nutrient enrichment have been identified as significant concerns in the estuary, and consequently, it is important to accurately estimate the magnitude of discharges and nutrient loads transported by inflows from both rivers. Two methods for estimating discharge and nutrient loads from tidally affected reaches of the Myakka and Peace Rivers were compared. The first method was a tidal-estimation method, in which discharge and nutrient loads were estimated based on stage, water-velocity, discharge, and water-quality data collected near the mouths of the rivers. The second method was a traditional basin-ratio method in which discharge and nutrient loads at the mouths were estimated from discharge and loads measured at upstream stations. Stage and water-velocity data were collected near the river mouths by submersible instruments, deployed in situ, and discharge measurements were made with an acoustic Doppler current profiler. The data collected near the mouths of the Myakka River and Peace River were filtered, using a low-pass filter, to remove daily mixed-tide effects with periods less than about 2 days. The filtered data from near the river mouths were used to calculate daily mean discharge and nutrient loads. These tidal-estimation-method values were then compared to the basin-ratio-method values. Four separate 30-day periods of differing streamflow conditions were chosen for monitoring and comparison. Discharge and nutrient load estimates computed from the tidal-estimation and basin-ratio methods were most similar during high-flow periods. However, during high flow, the values computed from the tidal-estimation method for the Myakka and Peace Rivers were consistently lower than the values computed from the basin-ratio method. There were substantial differences between discharges and nutrient loads computed from the tidal-estimation and basin-ratio methods during low-flow periods. Furthermore, the differences between the methods were not consistent. Discharges and nutrient loads computed from the tidal-estimation method for the Myakka River were higher than those computed from the basin-ratio method, whereas discharges and nutrients loads computed by the tidal-estimation method for the Peace River were not only lower than those computed from the basin-ratio method, but they actually reflected a negative, or upstream, net movement. Short-term tidal measurement results should be used with caution, because antecedent conditions can influence the discharge and nutrient loads. Continuous tidal data collected over a 1- or 2-year period would be necessary to more accurately estimate the tidally affected discharge and nutrient loads for the Myakka and Peace River Basins.

Florida↗