Geology ReportsSearch

SEARCH · Geology Reports

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 415 records · Page 23Linked to original sources

Hydrology of Area 62, Northern Great Plains and Rocky Mountain Coal Provinces, New Mexico and Arizona

This report summarizes available hydrologic data for Area 62 and will aid leasing decisions, and the preparation and appraisal of environmental impact studies and mine-permit applications. Area 62 is located at the southern end of the Rocky Mountain Coal Province in parts of New Mexico and Arizona and includes approximately 9,500 square miles. Surface mining alters, at least temporarily, the environment; if the areas are unreclaimed, there can be long-term environmental consequences. The land-ownership pattern in Area 62 is complicated. The checkerboard pattern created by several types of ownership makes effective management of these lands difficult. The climate generally is semiarid with average annual precipitation ranging from 10 to 20 inches. Pinons, junipers, and grasslands cover most of the area, and much of it is used for grazing by livestock. Soils vary with landscape, differing from flood plains and hillslopes to mountain slopes. The major structural features of this area were largely developed during middle Tertiary time. The main structural features are the southern San Juan Basin and the Mogollon slope. Coal-bearing rocks are present in four Cretaceous rock units of the Mesaverde Group: the Gallup Sandstone, the Dileo Coal Member, and the Gibson Coal Member of the Crevasse Canyon Formation, and the Cleary Coal Member of the Menefee Formation. Area 62 is drained by Black Creek, the Puerco River, the Zuni River, Carrizo Wash-Largo Creek, and the Rio San Jose. Only at the headwaters of the Zuni River is the flow perennial. The streamflow-gaging station network consists of 25 stations operated for a variety of needs. Streamflow changes throughout the year with variation related directly to rainfall and snowmelt. Base flow in Area 62 is zero indicating no significant ground-water discharge. Mountainous areas contribute the highest mean annual runoff of 1.0 inch. Very few water-quality data are available for the surface-water stations. Of the nine surface-water stations that have water-quality data, only one has chemical analyses from more than 10 samples. Therefore, sufficient data to characterize the area in detail are not available. Suspended sediment data are available only for a few surface- water stations in the area. Erosion rates generally are less than 1 acre-foot per square mile per year. Greater erosion rates are found within the badland areas. Water levels are periodically measured at 21 selected wells in Area 62. These observation wells are located mostly along the Rio San Jose and northeast of Gallup, New Mexico. The recharge to ground-water aquifers generally coincide with areas of greater precipitation in the mountainous areas. Depth to water below land surface is generally less than 200 feet. Well yields of 100 gallons per minute are common in most of the area. Ground-water quality is variable both within each aquifer and between aquifers. Water quality generally is best near recharge areas. Historical and current data related to stream discharge, water quality, and suspended sediment are available from computer files in the U.S. Geological Survey's National Water Data Storage and Retrieval System (WATSTORE) and through the National Water Data Exchange (NAWDEX).

Open-File Report

Hydrology of Area 62, Northern Great Plains and Rocky Mountain Coal Provinces, New Mexico and Arizona

This report summarizes available hydrologic data for Area 62 and will aid leasing decisions, and the preparation and appraisal of environmental impact studies and mine-permit applications. Area 62 is located at the southern end of the Rocky Mountain Coal Province in parts of New Mexico and Arizona and includes approximately 9,500 square miles. Surface mining alters, at least temporarily, the environment; if the areas are unreclaimed, there can be long-term environmental consequences. The land-ownership pattern in Area 62 is complicated. The checkerboard pattern created by several types of ownership makes effective management of these lands difficult. The climate generally is semiarid with average annual precipitation ranging from 10 to 20 inches. Pinons, junipers, and grasslands cover most of the area, and much of it is used for grazing by livestock. Soils vary with landscape, differing from flood plains and hillslopes to mountain slopes. The major structural features of this area were largely developed during middle Tertiary time. The main structural features are the southern San Juan Basin and the Mogollon slope. Coal-bearing rocks are present in four Cretaceous rock units of the Mesaverde Group: the Gallup Sandstone, the Dileo Coal Member, and the Gibson Coal Member of the Crevasse Canyon Formation, and the Cleary Coal Member of the Menefee Formation. Area 62 is drained by Black Creek, the Puerco River, the Zuni River, Carrizo Wash-Largo Creek, and the Rio San Jose. Only at the headwaters of the Zuni River is the flow perennial. The streamflow-gaging station network consists of 25 stations operated for a variety of needs. Streamflow changes throughout the year with variation related directly to rainfall and snowmelt. Base flow in Area 62 is zero indicating no significant ground-water discharge. Mountainous areas contribute the highest mean annual runoff of 1.0 inch. Very few water-quality data are available for the surface-water stations. Of the nine surface-water stations that have water-quality data, only one has chemical analyses from more than 10 samples. Therefore, sufficient data to characterize the area in detail are not available. Suspended sediment data are available only for a few surface- water stations in the area. Erosion rates generally are less than 1 acre-foot per square mile per year. Greater erosion rates are found within the badland areas. Water levels are periodically measured at 21 selected wells in Area 62. These observation wells are located mostly along the Rio San Jose and northeast of Gallup, New Mexico. The recharge to ground-water aquifers generally coincide with areas of greater precipitation in the mountainous areas. Depth to water below land surface is generally less than 200 feet. Well yields of 100 gallons per minute are common in most of the area. Ground-water quality is variable both within each aquifer and between aquifers. Water quality generally is best near recharge areas. Historical and current data related to stream discharge, water quality, and suspended sediment are available from computer files in the U.S. Geological Survey's National Water Data Storage and Retrieval System (WATSTORE) and through the National Water Data Exchange (NAWDEX).

Water-Resources Investigations Report

Nitrogen concentrations and loads for the Connecticut River at Middle Haddam, Connecticut, computed with the use of autosampling and continuous measurements of water quality for water years 2009 to 2014

The daily and annual loads of nitrate plus nitrite and total nitrogen for the Connecticut River at Middle Haddam, Connecticut, were determined for water years 2009 to 2014. The analysis was done with a combination of methods, which included a predefined rating curve method for nitrate plus nitrite and total nitrogen for water years 2009 to 2011 and a custom rating curve method that included sensor measurements of nitrate plus nitrite nitrogen concentration and turbidity along with mean daily flow to determine total nitrogen loads for water years 2011 to 2014. Instantaneous concentrations of total nitrogen were estimated through the use of a regression model based on sensor measurements at 15-minute intervals of nitrate plus nitrite nitrogen and turbidity for water years 2011 to 2014. Annual total nitrogen loads at the Connecticut River at Middle Haddam ranged from 12,900 to 19,200 metric tons, of which about 42 to 49 percent was in the form of nitrate plus nitrite. The mean 95-percent prediction intervals on daily total nitrogen load estimates were smaller from the custom model, which used sensor data, than those calculated by the predefined model. Annual total nitrogen load estimates at the Connecticut River at Middle Haddam were compared with the upstream load estimates at the Connecticut River at Thompsonville, Conn. Annual gains in total nitrogen loads between the two stations ranged from 3,430 to 6,660 metric tons. These increases between the two stations were attributed to the effects of increased urbanization and to combined annual discharges of 1,540 to 2,090 metric tons of nitrogen from 24 wastewater treatment facilities in the drainage area between the two stations. The contribution of total nitrogen from wastewater discharge between the two stations had declined substantially before the beginning of this study and accounted for from 31 to 52 percent of the gain in nitrogen load between the Thompsonville and Middle Haddam sites.

Connecticut

The use of NOAA AVHRR data for assessment of the urban heat sland effect

A vegetation index and a radiative surface temperature were derived from satellite data acquired at approximately 1330 LST for each of 37 cities and for their respective nearby rural regions from 28 June through 8 August 1991. Urban–rural differences for the vegetation index and the surface temperatures were computed and then compared to observed urban–rural differences in minimum air temperatures. The purpose of these comparisons was to evaluate the use of satellite data to assess the influence of the urban environment on observed minimum air temperatures (the urban heat island effect). The temporal consistency of the data, from daily data to weekly, biweekly, and monthly intervals, was also evaluated. The satellite-derived normalized difference (ND) vegetation-index data, sampled over urban and rural regions composed of a variety of land surface environments, were linearly related to the difference in observed urban and rural minimum temperatures. The relationship between the ND index and observed differences in minimum temperature was improved when analyses were restricted by elevation differences between the sample locations and when biweekly or monthly intervals were utilized. The difference in the ND index between urban and rural regions appears to be an indicator of the difference in surface properties (evaporation and heat storage capacity) between the two environments that are responsible for differences in urban and rural minimum temperatures. The urban and rural differences in the ND index explain a greater amount of the variation observed in minimum temperature differences than past analyses that utilized urban population data. The use of satellite data may contribute to a globally consistent method for analysis of urban heat island bias.

Journal of Applied Meteorology and Climatology

Partly cloudy with a chance of migration: Weather, radars, and aeroecology

Aeroecology is an emerging scientific discipline that integrates atmospheric science, Earth science, geography, ecology, computer science, computational biology, and engineering to further the understanding of biological patterns and processes. The unifying concept underlying this new transdisciplinary field of study is a focus on the planetary boundary layer and lower free atmosphere (i.e., the aerosphere), and the diversity of airborne organisms that inhabit and depend on the aerosphere for their existence. Here, we focus on the role of radars and radar networks in aeroecological studies. Radar systems scanning the atmosphere are primarily used to monitor weather conditions and track the location and movements of aircraft. However, radar echoes regularly contain signals from other sources, such as airborne birds, bats, and arthropods. We briefly discuss how radar observations can be and have been used to study a variety of airborne organisms and examine some of the many potential benefits likely to arise from radar aeroecology for meteorological and biological research over a wide range of spatial and temporal scales. Radar systems are becoming increasingly sophisticated with the advent of innovative signal processing and dual-polarimetric capabilities. These capabilities should be better harnessed to promote both meteorological and aeroecological research and to explore the interface between these two broad disciplines. We strongly encourage close collaboration among meteorologists, radar scientists, biologists, and others toward developing radar products that will contribute to a better understanding of airborne fauna.

Bulletin of the American Meteorological Society

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