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 433 records · Page 24Linked to original sources
USGS2020

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↗
USGS2021

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↗
USGS2014

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↗
USGS2012

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↗
USGS1991

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↗
USGS2023

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↗
USGS2008

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↗
USGS2023

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↗
USGS2009

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↗
USGS1997

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↗
USGS2005

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↗
USGS2005

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↗
USGS2009

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↗
USGS1997

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↗
USGS2021

Analysis of Escherichia coli, total recoverable iron, and dissolved selenium concentrations, loading, and identifying data gaps for selected 303(d) listed streams, Grand Valley, western Colorado, 1980–2018

Tributaries to the Colorado River in the Grand Valley in western Colorado (segment COLCLC13b) have been placed on the State of Colorado 303(d) list as impaired for Escherichia coli (E. coli) , total recoverable iron, and dissolved selenium. The Colorado Department of Public Health and Environment Water Quality Control Division is required to develop total maximum daily loads for these constituents in these tributaries. The U.S. Geological Survey, in cooperation with the Grand Valley Drainage District and Colorado Water Conservation Board, conducted a study to (1) characterize concentrations, loads, and load reductions for E. coli , total recoverable iron, and dissolved selenium using existing data and (2) identify water-quality data gaps to inform future monitoring strategies. This study analyzed water-quality and streamflow data for 3 main-stem sites (2 sites along the Colorado River and 1 site along the Gunnison River) and 29 selected sites on tributaries to the Colorado River. Sample data were available at five sites along Adobe Creek and at six sites along Leach Creek, the two tributaries in the study area that are impaired for E. coli . All geometric mean E. coli concentrations at sites along Adobe Creek and Leach Creek exceeded the State recreational use standard of 126 colony forming units per 100 milliliters (CFU/100 mL). In Adobe Creek, E. coli concentrations in samples ranged from 45.7 to more than 2,420 CFU/100 mL (method upper reporting limit for undiluted samples), and geometric mean concentrations at sites ranged from 301 to 1,180 CFU/100 mL. The E. coli concentrations generally increased in the downstream direction in Adobe Creek; however, increases were not seen between all sites. The largest downstream increase in E. coli concentration was measured between the two most upstream sites. In Leach Creek, concentrations of E. coli in samples ranged from 25.9 to more than 2,420 CFU/100 mL, and geometric mean concentrations at sites ranged from 160 to 259 CFU/100 mL. The E. coli concentrations showed no consistent downgradient increase in Leach Creek. In fact, some of the highest E. coli concentrations were measured at the most upstream site, Leach Creek at Summer Hill Drive. Total recoverable iron concentrations and loads were evaluated at 15 tributary sites for samples collected from August 1993 to February 2018. Median total recoverable iron concentrations ranged from 211 to 4,670 micrograms per liter (µg/L). The chronic aquatic-life water-quality standard (1,000 µg/L) was exceeded in most irrigation season (April through October) samples but was rarely exceeded in nonirrigation season (November through March) samples. Concentrations were often an order of magnitude higher in samples collected during irrigation season than in samples collected during nonirrigation season. None of the sites had enough concurrent total recoverable iron and streamflow data to compute annual loads. As with E. coli , the lack of concurrent total recoverable iron and streamflow information represents a data gap, which needs to be addressed to compute annual loads. Dissolved selenium concentrations and loads were evaluated at 20 tributary sites using discrete water-quality data collected 1991–2018. Dissolved selenium concentrations were higher during nonirrigation season than during irrigation season at tributary sites. However, irrigation season dissolved selenium loads were generally higher than nonirrigation selenium loads, because streamflows were higher during irrigation season. Regression analysis was used to estimate daily dissolved selenium concentrations and loads at three main-stem sites for water years (WYs) 1980–2018 (Gunnison River near Grand Junction and Colorado River near Colorado-Utah State Line) and WYs 2002–18 (Colorado River near Cameo). A trend analysis of dissolved selenium concentrations and loads was completed for these sites from the same respective starting dates but ending in 2017. A continuing downward trend in dissolved selenium concentration was observed at all sites and across all seasonal designations of the analysis. The dissolved selenium concentration decreased by 0.12 µg/L from WY 2002 to 2017 at Colorado River near Cameo, representing an 18-percent decrease during the time period. The dissolved selenium concentration at Gunnison River near Grand Junction decreased by 4.2 µg/L from WY 1980 to 2017, representing a 56-percent decrease overall. During the same time period, dissolved selenium concentration at Colorado River near Colorado-Utah State Line decreased by 3.8 µg/L, representing a 56-percent decrease overall. A downward trend in dissolved selenium load was also observed at all sites and across all seasonal designations of the analysis. The relative contribution of dissolved selenium from the Grand Valley near Grand Junction was estimated by comparing loads at main-stem sites bracketing the study area. The two upstream sites, Colorado River near Cameo and Gunnison River near Grand Junction, contributed 60,300 cumulative pounds and 251,000 cumulative pounds, respectively, during WYs 2002–18. At the furthest downstream site, Colorado River near Colorado-Utah State Line, 490,000 cumulative pounds were estimated during the same time period, indicating that the region between Whitewater and State line contributed approximately 179,000 cumulative pounds or a mean annual load of 10,500 lb/yr. Grand Valley dissolved selenium contributions appear to be stable during WYs 2002–18.

Colorado↗
USGS2006

Suspended-sediment yields and stream-channel processes on Judy's Branch watershed in the St. Louis Metro East region in Illinois

Judy's Branch watershed, a small basin (8.64 square miles) in the St. Louis Metro East region in Illinois, was selected as a pilot site to determine suspended-sediment yields and stream-channel processes in the bluffs and American Bottoms (expansive low-lying valley floor in the region). Suspended-sediment and stream-chan-nel data collected and analyzed for Judy's Branch watershed are presented in this report to establish a baseline of data for water-resource managers to evaluate future stream rehabilitation and manage-ment alternatives. The sediment yield analysis determines the amount of sediment being delivered from the watershed and two subwatersheds: an urban tributary and an undeveloped headwater (pri-marily agricultural). The analysis of the subwater-sheds is used to compare the effects of urbanization on sediment yield to the river. The stream-channel contribution to sediment yield was determined by evaluation of the stream-channel processes operat-ing on the streambed and banks of Judy's Branch watershed. Bank stability was related to hydrologic events, bank stratigraphy, and channel geometry through model development and simulation. The average suspended-sediment yield from two upland subwatersheds (drainage areas of 0.23 and 0.40 sq.mi. was 1,163 tons per square mile per year (tons/sq.mi.-year) between July 2000 and June 2004. The suspended-sediment yield at the Route 157 station was 2,523 tons/sq.mi.-year, near the outlet of Judy's Branch watershed (drainage area = 8.33 sq.mi.). This is approximately 1,360 tons/sq.mi.-year greater than the average at the upland stations for the same time period. This result is unexpected in that, generally, the suspended-sediment yield decreases as the watershed area increases because of sediment stored in the channel and flood plain. The difference indicates a possible increase in yield from a source, such as bank retreat, and supports the concept that land-use changes increase stream-flows that may in turn result in higher rates of bank retreat. Utilizing both bank-rod data and resurveyed cross-section data, it was determined that approxi-mately half of the suspended- sediment yield at Route 157 during July 2000-June 2004 came from bank retreat. Given that bank retreat can be a substantial portion of the sediment yield, understanding bank stability processes is important. Bank stability can be assessed mathematically by computing the factor of safety, which is defined by the ratio of the shear strength (resisting force) along the failure surface and the shear stress (driving gravitational force). Once the factor of safety falls below one, the bank theoretically becomes unstable. Bank-stability conditions were related to hydrologic events, bank type, and channel geometry through model develop-ment and simulation. The most common type of bank in the watershed consists of cohesive alluvial soil deposits overlying a stiff glacial till. A stabil-ity chart for different bank types was developed using a bank-stability analysis. Banks steeper than 70 degrees and higher than from 10 to 11.5 feet (depending on bank type) become at risk for mass failure in the watershed under conditions that pro-mote saturation of the bank and a sudden drop in the river level.

Scientific Investigations Report↗
USGS2017

Observations and 3D hydrodynamics-based modeling of decadal-scale shoreline change along the Outer Banks, North Carolina

Long-term decadal-scale shoreline change is an important parameter for quantifying the stability of coastal systems. The decadal-scale coastal change is controlled by processes that occur on short time scales (such as storms) and long-term processes (such as prevailing waves). The ability to predict decadal-scale shoreline change is not well established and the fundamental physical processes controlling this change are not well understood. Here we investigate the processes that create large-scale long-term shoreline change along the Outer Banks of North Carolina, an uninterrupted 60 km stretch of coastline, using both observations and a numerical modeling approach. Shoreline positions for a 24-yr period were derived from aerial photographs of the Outer Banks. Analysis of the shoreline position data showed that, although variable, the shoreline eroded an average of 1.5 m/yr throughout this period. The modeling approach uses a three-dimensional hydrodynamics-based numerical model coupled to a spectral wave model and simulates the full 24-yr time period on a spatial grid running on a short (second scale) time-step to compute the sediment transport patterns. The observations and the model results show similar magnitudes (O(10 5 m 3 /yr)) and patterns of alongshore sediment fluxes. Both the observed and the modeled alongshore sediment transport rates have more rapid changes at the north of our section due to continuously curving coastline, and possible effects of alongshore variations in shelf bathymetry. The southern section with a relatively uniform orientation, on the other hand, has less rapid transport rate changes. Alongshore gradients of the modeled sediment fluxes are translated into shoreline change rates that have agreement in some locations but vary in others. Differences between observations and model results are potentially influenced by geologic framework processes not included in the model. Both the observations and the model results show higher rates of erosion (∼−1 m/yr) averaged over the northern half of the section as compared to the southern half where the observed and modeled averaged net shoreline changes are smaller (<0.1 m/yr). The model indicates accretion in some shallow embayments, whereas observations indicate erosion in these locations. Further analysis identifies that the magnitude of net alongshore sediment transport is strongly dominated by events associated with high wave energy. However, both big- and small- wave events cause shoreline change of the same order of magnitude because it is the gradients in transport, not the magnitude, that are controlling shoreline change. Results also indicate that alongshore momentum is not a simple balance between wave breaking and bottom stress, but also includes processes of horizontal vortex force, horizontal advection and pressure gradient that contribute to long-term alongshore sediment transport. As a comparison to a more simple approach, an empirical formulation for alongshore sediment transport is used. The empirical estimates capture the effect of the breaking term in the hydrodynamics-based model, however, other processes that are accounted for in the hydrodynamics-based model improve the agreement with the observed alongshore sediment transport.

North Carolina↗
USGS2025

Methods for estimating selected low-flow statistics at gaged and ungaged stream sites in Massachusetts

The U.S. Geological Survey, in cooperation with the Massachusetts Department of Conservation and Recreation, Office of Water Resources, computed selected at-site streamflow statistics at U.S. Geological Survey streamgages in and near Massachusetts and developed regional regression equations for estimating selected streamflows at ungaged stream sites in Massachusetts. Two sets of regional regression equations were developed: (1) the “mainland” equations, for mainland Massachusetts excluding the area covered by the second set, and (2) the “southeastern” equations, for the Plymouth-Carver-Kingston-Duxbury aquifer area in southeastern Massachusetts and for Cape Cod. The regression equations and at-site statistics may be used by Federal, State, and local water managers in addressing water-resources issues relevant in Massachusetts. Regional regression analyses for the mainland equations were developed to estimate the following 27 streamflow statistics: 99-, 98-, 95-, 90-, 85-, 80-, 75-, 70-, 60-, and 50-percent flow durations; monthly June, July, August, and September 90- and 50-percent flow durations; February, June, and August median of the monthly means; harmonic mean; and medians of the following annual low-flow frequency statistics: 7-day; 7-day, 2-year; 7-day, 10-year; 30-day, 2-year; and 30-day, 10-year. The analyses used 81 streamgages with minimal to no regulations in and near Massachusetts. The regression analyses determined that four basin characteristics—drainage area, combined hydrologic soils A and B, streamflow variability index, and annual mean temperature—were the only significant explanatory variables for the different mainland equations. Regional regression equations were developed for the Plymouth-Carver-Kingston-Duxbury aquifer area in southeastern Massachusetts and Cape Cod, because surface-water drainage areas and groundwater contributing areas do not always coincide in this area of the State. The regression analyses to estimate 10 flow durations from the 99th to 50th percentiles used 18 streamflow sites with some occasional minor regulations—because there are few unregulated streams in southeastern Massachusetts. The analyses determined that groundwater contributing area and storage (combined water bodies and wetlands) were the only significant explanatory variables in the southeastern equations.

Massachusetts↗