Search USGS⌕ Search

SEARCH · Search USGS

Results for “Scientific Investigations Map”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,657 records · Page 92Linked to original sources

Analysis and mapping of post-fire hydrologic hazards for the 2002 Hayman, Coal Seam, and Missionary Ridge wildfires, Colorado

Wildfires caused extreme changes in the hydrologic, hydraulic, and geomorphologic characteristics of many Colorado drainage basins in the summer of 2002. Detailed assessments were made of the short-term effects of three wildfires on burned and adjacent unburned parts of drainage basins. These were the Hayman, Coal Seam, and Missionary Ridge wildfires. Longer term runoff characteristics that reflect post-fire drainage basin recovery expected to develop over a period of several years also were analyzed for two affected stream reaches: the South Platte River between Deckers and Trumbull, and Mitchell Creek in Glenwood Springs. The 10-, 50-, 100-, and 500-year flood-plain boundaries and water-surface profiles were computed in a detailed hydraulic study of the Deckers-to-Trumbull reach. The Hayman wildfire burned approximately 138,000 acres (216 square miles) in granitic terrain near Denver, and the predominant potential hazard in this area is flooding by sediment-laden water along the large tributaries to and the main stem of the South Platte River. The Coal Seam wildfire burned approximately 12,200 acres (19.1 square miles) near Glenwood Springs, and the Missionary Ridge wildfire burned approximately 70,500 acres (110 square miles) near Durango, both in areas underlain by marine shales where the predominant potential hazard is debris-flow inundation of low-lying areas. Hydrographs and peak discharges for pre-burn and post-burn scenarios were computed for each drainage basin and tributary subbasin by using rainfall-runoff models because streamflow data for most tributary subbasins were not available. An objective rainfall-runoff model calibration method based on nonlinear regression and referred to as the ?objective calibration method? was developed and applied to rainfall-runoff models for three burned areas. The HEC-1 rainfall-runoff model was used to simulate the pre-burn rainfall-runoff processes in response to the 100-year storm, and HEC-HMS was used for runoff hydrograph generation. Post-burn rainfall-runoff parameters were determined by adjusting the runoff-curve numbers on the basis of a weighting procedure derived from the U.S. Soil Conservation Service (now the National Resources Conservation Service) equation for precipitation excess and the effect of burn severity. This weighting procedure was determined to be more appropriate than simple area weighting because of the potentially marked effect of even small burned areas on the runoff hydrograph in individual drainage basins. Computed water-peak discharges from HEC-HMS models were increased volumetrically to account for increased sediment concentrations that are expected as a result of accelerated erosion after burning. Peak discharge estimates for potential floods in the South Platte River were increased by a factor that assumed a volumetric sediment concentration (Cv) of 20 percent. Flood hydrographs for the South Platte River and Mitchell Creek were routed down main-stem channels using watershed-routing algorithms included in the HEC-HMS rainfall-runoff model. In areas subject to debris flows in the Coal Seam and Missionary Ridge burned areas, debris-flow discharges were simulated by 100-year rainfall events, and the inflow hydrographs at tributary mouths were simulated by using the objective calibration method. Sediment concentrations (Cv) used in debris-flow simulations were varied through the event, and were initial Cv 20 percent, mean Cv approximately 31 percent, maximum Cv 48 percent, Cv 43 percent at the time of the water hydrograph peak, and Cv 20 percent for the duration of the event. The FLO-2D flood- and debris-flow routing model was used to delineate the area of unconfined debris-flow inundation on selected alluvial fan and valley floor areas. A method was developed to objectively determine the post-fire recovery period for the Hayman and Coal Seam burned areas using runoff-curve numbers (RCN) for all drainage basins for a 50-year period. A

Scientific Investigations Report↗

Water-level conditions in the confined aquifers of the New Jersey Coastal Plain, 2008

Groundwater-level altitudes in 10 confined aquifers of the New Jersey Coastal Plain were measured and evaluated to provide an overview of regional groundwater conditions during fall 2008. Water levels were measured in more than 900 wells in New Jersey, eastern Pennsylvania, and northern Delaware and potentiometric surface maps prepared for the confined Cohansey aquifer of Cape May County, the Rio Grande water-bearing zone, the Atlantic City 800-foot sand, the Piney Point, Vincentown, and the Wenonah-Mount Laurel aquifers, the Englishtown aquifer system, and the Upper, Middle, and Lower aquifers of the Potomac-Raritan-Magothy aquifer system. In 2008, the highest water-level altitudes were observed in the Vincentown aquifer (median, 78 ft) and the lowest in the Atlantic City 800-foot sand (median, -45 ft). Persistent, regionally extensive cones of depression were present within the potentiometric surfaces of the Englishtown aquifer system in east-central New Jersey, the Wenonah-Mount Laurel aquifer in east-central and southern New Jersey, the Upper, Middle, and Lower Potomac-Raritan-Magothy aquifers in southern New Jersey, and the Atlantic City 800-foot sand in the southeastern part of the State. Cones of depression in the potentiometric surfaces of the Upper Potomac-Raritan-Magothy and the Piney Point aquifers in east-central and southwestern New Jersey had broadened and deepened since 2003. Declining water levels in many of New Jersey’s confined Coastal Plain aquifers intensified during the late 1970s and early 1980s, prompting the designation of two water-supply Critical Areas by the New Jersey Department of Environmental Protection; Critical Areas 1 and 2 continued to be of concern. To address that concern, water-level changes were assessed in nearly 800 wells measured during the fall of 2003 and 2008, and potentiometric-surface difference maps for each aquifer were constructed and evaluated. In addition, water-level trends were calculated for 77 wells for the periods 2003–8 and 1998–2008 and for 73 wells for the period 1978–2008. From 2003 to 2008 small to moderate water-level changes were observed in many Coastal Plain aquifers in New Jersey, but in places, groundwater levels continued to decline substantially as a result of pumping. Groundwater levels in the Atlantic City 800-foot sand were lower in 2008 than in 2003; declines were greatest near pumping centers in eastern Atlantic County. Changes were less pronounced in Cape May County where water levels were, on average, 1 to 3 feet (ft) lower than those during the previous study (2003), except near Rio Grande where a localized cone of depression had formed as a result of increased withdrawals. Large and widespread declines occurred in the Piney Point aquifer in Cumberland County where water levels in and around the city of Bridgeton had fallen in excess of 100 ft since 2003, and by 30 ft to more than 60 ft in surrounding areas. Groundwater levels in the Wenonah-Mount Laurel aquifer and Englishtown aquifer system continued to recover in east-central New Jersey; however, groundwater levels in the Wenonah-Mount Laurel aquifer throughout the southern part of the State continued to decline. In the Upper Potomac-Raritan-Magothy aquifer, groundwater levels were substantially lower than in 2003 in parts of northern Ocean County but were stable in the area adjacent to Raritan Bay (Critical Area 1), and water levels continued to recover in southern New Jersey. In the Middle Potomac-Raritan-Magothy aquifer, water levels rose near Raritan Bay in Middlesex County; however, modest declines were recorded in interior areas of Monmouth and Ocean Counties. Groundwater levels in both the Middle and Lower Potomac-Raritan-Magothy aquifers were stable or rising within the regional cone of depression in Critical Area 2; beyond the critical area in southern New Jersey, however, water levels were slightly lower than in 2003. Analyses of long-term water-level changes indicate that from 1978 to 2008 downward trends occurred at 20 wells (27 percent), upward trends at 27 wells (37 percent), and trends at 26 wells (36 percent) were insubstantial. Sustained, long-term declines were observed most often at wells within the Atlantic City 800-foot sand and at wells in the Piney Point aquifer in southern New Jersey, in which rates of decline were as great as 1.4 feet/year. Upward water-level trends were observed frequently at wells screened in the Englishtown aquifer system and the Wenonah-Mount Laurel aquifer in Critical Area 1 in east-central New Jersey, and in the Potomac-Raritan-Magothy aquifer system in parts of Critical Area 1 and throughout most of Critical Area 2 in southern New Jersey. Annual rates of upward change were as great as 3.9 and 5.6 ft/yr in the Englishtown aquifer system and Wenonah-Mount Laurel aquifer, respectively. Among the units of the Potomac-Raritan-Magothy aquifer system, annual rates of recovery were greatest in the Lower aquifer. From 1998 to 2008, downward water-level trends were observed at 22 wells (29 percent), upward trends were observed at 21 wells (27 percent), and insubstantial trends at 34 wells (44 percent). Downward trends were detected most often at wells open to the Piney Point aquifer and the Atlantic City 800-foot sand. Upward water-level trends were most frequent in wells open to the Englishtown aquifer system in Critical Area 1 and in wells within the Potomac-Raritan-Magothy aquifer system in southern New Jersey.

New Jersey↗

Assessment of ethylene dibromide, dibromochloropropane, other volatile organic compounds, radium isotopes, radon, and inorganic compounds in groundwater and spring water from the Crouch Branch and McQueen Branch aquifers near McBee, South Carolina, 2010-2012

Public-supply wells near the rural town of McBee, in southwestern Chesterfield County, South Carolina, have provided potable water to more than 35,000 residents throughout Chesterfield County since the early 1990s. Groundwater samples collected between 2002 and 2008 in the McBee area by South Carolina Department of Health and Environmental Control (DHEC) officials indicated that groundwater from two public-supply wells was characterized by the anthropogenic compounds ethylene dibromide (EDB) and dibromochloropropane (DBCP) at concentrations that exceeded their respective maximum contaminant levels (MCLs) established by the U.S. Environmental Protection Agency’s (EPA) National Primary Drinking Water Regulations (NPDWR). Groundwater samples from all public-supply wells in the McBee area were characterized by the naturally occurring isotopes of radium-226 and radium-228 at concentrations that approached, and in one well exceeded, the MCL for the combined isotopes. The local water utility installed granulated activated carbon filtration units at the two EDB- and DBCP-contaminated wells and has, since 2011, shut down these two wells. Groundwater pumped by the remaining public-supply wells is currently (2014) centrally treated at a water-filtration plant. To assess the occurrence, distribution, and potential sources of the anthropogenic and naturally occurring compounds detected in groundwater in the McBee area, samples of groundwater and spring water were collected from public-supply, domestic-supply, agricultural-supply, and monitoring wells and springs, respectively, between 2010 and 2012 by the U.S. Geological Survey. The water samples were analyzed for concentrations of EDB, DBCP, other volatile organic compounds (VOCs), radium-226 and radium-228, radon, and inorganic compounds. All wells sampled were screened in the shallow Crouch Branch aquifer, the deeper McQueen Branch aquifer, or, for most public-supply wells, both aquifers. In areas where no wells existed or wells could not be installed, passive samplers that adsorb EDB, DBCP, and various VOCs, were installed in the shallow subsurface. A representative groundwater flow pathway to each public supply well and selected other wells was determined by using a calibrated three-dimensional groundwater-flow model of the Atlantic Coastal Plain in Chesterfield County and particle-tracking analysis. The aerial extent of the groundwater flow pathway to public-supply wells was mapped by using chlorofluorocarbon-concentration based, apparent-age dates of the groundwater. The water-quality data collected between 2010 and 2012, in conjunction with groundwater flow pathways and historical aerial photographs of land uses near McBee, indicate an area where EDB-, DBCP-, 1,2-dichloropropane-, 1,3-dichloropropane-, and carbon disulfide-contaminated groundwater exists in the Crouch Branch aquifer in the Cedar Creek Basin and north of McBee and is most likely related to the past use of these compounds between the early 1900s and the 1980s as soil fumigants in predominately agricultural areas north of McBee. The highest EDB concentration detected (18.6 micrograms per liter) during the 3-year study was in a groundwater sample from an agricultural-supply well located north of McBee. Other VOCs, such as dichloromethane and 1,1,2-trichloroethane, also were detected in groundwater samples from this EDB-contaminated agricultural-supply well but are from unknown source(s). The fact that the agricultural area north of McBee is located in a recharge area for the Crouch Branch aquifer most likely facilitated the groundwater contamination in this area. DBCP-contaminated groundwater detected in three public-supply wells south of McBee in the deeper McQueen Branch aquifer appears to be related to past soil fumigation practices that used DBCP in agricultural areas located south of McBee. One of the three DBCP-contaminated public-supply wells also contained EDB, most likely present in groundwater due to the release of leaded gasolines that contained EDB as a fuel additive between the 1940s and 1970s. A gasoline-source of EDB, rather than a soil-fumigation source, is supported by the co-detection in groundwater from the well of 1,2-dichloroethane, a lead scavenger compound also added to leaded gasoline. Groundwater pumped from two public-supply wells located within and to the east of the McBee town limits and one domestic-supply well east of McBee was characterized by the detection of 1,1-dichloroethane, trichloroethylene, 1,1-dichloroethylene, and perchloroethylene. Groundwater flow pathways determined for these wells indicate that the potential source(s) of these compounds detected in one public-supply well and the domestic-supply well may be located within the McBee town limits, and that the potential source(s) of these compounds detected in the public-supply well to the east of McBee may be located in an area north of McBee formerly used for agriculture, but used for industry since at least the 1970s. Radium isotopes (defined in this study as the sum of radium-226 and radium-228 concentrations) and radon were detected in all wells sampled in the McBee area between 2010 and 2012. Wells characterized by radium isotope concentrations in groundwater that exceeded the MCL of 5.0 picocuries per liter were also characterized by specific conductance values greater than 30 microsiemens per centimeter and clustered north of McBee in a predominately agricultural area, and in agricultural and urban areas located within and east of McBee. The elevated specific conductance values measured in groundwater from these wells most likely are due to recharge by water mineralized by fertilizer application in agricultural areas, or due to the recharge by water mineralized by septic-tank drain-field effluent near urban areas. Radon was detected in groundwater from all wells sampled, and radon concentrations in groundwater from three monitoring wells exceeded the proposed MCL of 300 picocuries per liter. Concentrations of uranium in groundwater in the McBee area increased with increased groundwater-sample depth, most likely due to the proximity of the sample-collection location to basement rock that contains uranium-bearing minerals.

South Carolina↗

Bathymetric and Velocimetric Survey and Assessment of Habitat for Pallid Sturgeon on the Mississippi River in the Vicinity of the Proposed Interstate 70 Bridge at St. Louis, Missouri

A bathymetric and velocimetry survey was conducted on the Mississippi River in the vicinity of a proposed new bridge for Interstate 70 at St. Louis, Missouri. A multibeam echo sounder mapping system and an acoustic Doppler current profiler were used to obtain channel-bed elevations and vertically averaged and near-bed velocities for a 3,545-foot (1,080-meter) long reach of the Mississippi River approximately 1,935 feet (590 meters) wide from the Illinois to Missouri banks. Data from the 2009 survey were used to determine the conditions of the benthic habitat in the vicinity of the proposed Interstate 70 bridge. The channel-bed elevations ranged from approximately 346 feet (105.46 meters) to 370 feet (112.78 meters) above the North American Vertical Datum of 1988 in a majority of the channel except for the channel banks. Large dune features up to 12.5 feet (3.81 meters) high were present in the middle of the channel, and numerous smaller dunes and many ripples as smaller features were superimposed on the larger dunes. However, it is uncertain if the large dune features present in mid-channel are long-term features or an artifact of the seasonal flooding on the Mississippi River. A substantial scour depression was present on the right descending bank (Missouri side) near the downstream end of the study area, as well as other smaller scour holes near the instream barge mooring structures on the Missouri bank. The vertically averaged velocities acquired with the acoustic Doppler current profiler ranged from approximately 2 feet per second (0.61 meters per second) along the channel margins to approximately 7.0 feet per second (2.13 meters per second) in the main channel, with an average velocity of 5.5 feet per second (1.68 meters per second) in mid-channel. The orientation of the vertically averaged velocity vectors showed flow crossing from the Illinois bank to the Missouri bank from upstream to downstream in the study area, which was confirmed by the orientation of the large dune features in mid-channel and a shift in the channel thalweg from the Illinois bank to the Missouri bank. The near-bottom velocities acquired with the acoustic Doppler current profiler ranged from 0.3 to 7.0 feet per second (0.09 to 2.13 meters per second), and the effects of the large dune features were apparent in the more random scattering of the velocity vectors, the low velocities downstream from the dunes, and higher velocities near the crests of the dunes. Despite the considerable physical complexity of this site because of the arrangement of large sand dunes in the middle of the channel, existing studies do not document persistent use of these deep, fast, main-channel habitats by pallid sturgeon. Narrow channel-margin areas on both banks having relatively low velocity, high depth slope, and high velocity gradients are similar to adult migration habitats as documented on the Missouri River downstream from Kansas City, Missouri. Although the reach generally lacks features associated with sturgeon habitat selection on the Middle Mississippi River, the barge mooring areas on the right descending bank have topographic complexity and contain large woody debris and small patches of probable gravel-cobble substrate that may have positive habitat value for sturgeon or other species. Furthermore, telemetry studies have documented sturgeon migrating upstream and downstream through this reach as adults, and they probably drift downstream through this reach as free-embryo larvae. Successful upstream migration may depend on availability of areas with hydraulic complexity and relatively low velocities, as presently exist on the margins of the site. Additionally, complexity at the channel margin may provide areas where larvae settle out from drifting in the main current or may act to slow bulk drift rates. Construction of bridge piers close to the banks will likely alter hydraulics and sediment transport on the channel margins and may result in substanti

Scientific Investigations Report↗

Sand resources, regional geology, and coastal processes of the Chandeleur Islands Coastal System: An evaluation of the Breton National Wildlife Refuge

Breton National Wildlife Refuge, the Chandeleur Islands chain in Louisiana, provides habitat and nesting areas for wildlife and is an initial barrier protecting New Orleans from storms. The U.S. Geological Survey (USGS) in partnership with the University of New Orleans Pontchartrain Institute for Environmental Sciences undertook an intensive study that included (1) an analysis of island change based on historical maps and remotely sensed shoreline and topographic data; (2) a series of lidar surveys at 3- to 4-month intervals after Hurricane Katrina to determine barrier island recovery potential; (3) a discussion of sea level rise and effects on the islands; (4) an analysis of sea floor evolution and sediment dynamics in the refuge over the past 150 years; (5) an assessment of the local sediment transport and sediment resource availability based on the bathymetric and subbottom data; (6) a carefully selected core collection effort to groundtruth the geophysical data and more fully characterize the sediments composing the islands and surrounds; (7) an additional survey of the St. Bernard Shoals to assess their potential as a sand resource; and (8) a modeling study to numerically simulate the potential response of the islands to the low-intensity, intermediate, and extreme events likely to affect the refuge over the next 50 years. Results indicate that the islands have become fragmented and greatly diminished in subaerial extent over time: the southern islands retreating landward as they reorganize into subaerial features, the northern islands remaining in place. Breton Island, because maintenance of the Mississippi River-Gulf Outlet (MRGO) outer bar channel requires dredging, is deprived of sand sufficient to sustain itself. Regional sediment transport trends indicate that large storms are extremely effective in transporting sand and controlling the shoreline development and barrier island geometry. Sand is transported north and south from a divergent zone near Monkey Bayou at the southern end of the Chandeleur Islands. Numerical simulation of waves and sediment transport supports the geophysical results and indicates that vast areas of the lower shoreface are affected and are undergoing erosion during storm events, that there is little or no fair weather mechanism to rework material into the littoral system, and that as a result, there is a net loss of sediment from the system. Lidar surveys revealed that the island chain immediately after Hurricane Katrina lost about 84 percent of its area and about 92 percent of its prestorm volume. Marsh platforms that supported the islands’ sand prior to the storm were reduced in width by more than one-half. Repeated lidar surveys document that in places the shoreline has retreated about 100 m under the relatively low-energy waves since Hurricanes Katrina and Rita; however, this retreat is nonuniform. Recent high-resolution geophysical surveys of the sea floor and subsurface within 5–6 km of the Chandeleur Islands during 2006 and 2007 show that, in addition to the sand that is rebuilding portions of the island chain, a large volume of sand is contained in Hewes Point, in an extensive subtidal spit platform that has formed at the northern end of the Chandeleur Islands. Hewes Point appears to be the depositional terminus of the alongshore transport system. In the southern Chandeleurs, sand is being deposited in a broad tabular deposit near Breton Island called the southern offshore sand sheet. These two depocenters account for approximately 70 percent of the estimated sediment volume located in potential borrow sites. An additional large potential source of sand for restoration lies in the St. Bernard Shoals, which are estimated to contain approximately 200 × 10 6 m 3 of sand. Successful restoration planning for the Breton National Wildlife Refuge should mimic the natural processes of early stages of barrier island evolution including lateral transport to the flanks of the island chain from a centralized sand source that will ultimately enhance the ability of the islands to naturally build backbarrier marsh, dunes, and a continuous sandy shoreline. Barrier island sediment nourishment should be executed with the understanding that gulf shoreline erosion is inevitable but that island area can be maintained and enhanced during retreat (thus significantly prolonging the life of the island chain) with strategic sand placement.

Louisiana↗

Geohydrology, geochemistry, and groundwater simulation (1992-2011) and analysis of potential water-supply management options, 2010-60, of the Langford Basin, California

Groundwater withdrawals began in 1992 from the Langford Basin within the Fort Irwin National Training Center (NTC), California. From April 1992 to December 2010, approximately 12,300 acre-feet of water (averaging about 650 acre-feet per year) has been withdrawn from the basin and transported to the adjacent Irwin Basin. Since withdrawals began, water levels in the basin have declined by as much as 40 feet, and the quality of the groundwater withdrawn from the basin has deteriorated. The U.S. Geological Survey collected geohydrologic data from Langford Basin during 1992–2011 to determine the quantity and quality of groundwater available in the basin. Geophysical surveys, including gravity, seismic refraction, and time-domain electromagnetic induction surveys, were conducted to determine the depth and shape of the basin, to delineate depths to the Quaternary-Tertiary interface, and to map the depth to the water table and changes in water quality. Data were collected from existing wells and test holes, as well as 11 monitor wells that were installed at 5 sites as part of this study. Water-quality samples collected from wells in the basin were used to determine the groundwater chemistry within the basin and to delineate potential sources of poor-quality groundwater. Analysis of stable isotopes of oxygen and hydrogen in groundwater indicates that present-day precipitation is not a major source of recharge to the basin. Tritium and carbon-14 data indicate that most of the basin was recharged prior to 1952, and the groundwater in the basin has an apparent age of 12,500 to 30,000 years. Recharge to the basin, estimated to be less than 50 acre-feet per year, has not been sufficient to replenish the water that is being withdrawn from the basin. A numerical groundwater-flow model was developed for the Langford Basin to better understand the aquifer system used by the Fort Irwin NTC as part of its water supply, and to provide a tool to help manage groundwater resources at the NTC. Measured groundwater-level declines since the initiation of withdrawals (1992–2011) were used to calibrate the groundwater-flow model. The simulated recharge was about 46 acre-feet per year, including approximately 6 acre-feet per year of natural recharge derived from precipitation runoff and as much as 40 acre-feet per year of underflow from the Irwin Basin. Between April 1992 and December 2010, an average of about 650 acre-feet per year of water was withdrawn from the Langford Basin. Groundwater withdrawals in excess of natural recharge resulted in a net loss of 11,670 acre-feet of groundwater storage within the basin for the simulation period. The Fort Irwin NTC is considering various groundwater-management options to address the limited water resources in the Langford Basin. The calibrated Langford Basin groundwater-flow model was used to evaluate the hydrologic effects of four groundwater-withdrawal scenarios being considered by the Fort Irwin NTC over the next 50 years (January 2011 through December 2060). Continuation of the 2010 withdrawal rate in the three existing production wells will result in 70 feet of additional drawdown in the central part of the basin. Redistributing the 2010 withdrawal rate equally to the three existing wells and two proposed new wells in the northern and southern parts of the basin would result in about 10 feet less drawdown in the central part of the basin but about 100 feet of additional drawdown in the new well in the northern part of the basin and about 50 feet of additional drawdown in the new well in the southern part of the basin. Reducing the withdrawals from the three existing production wells in the central part of the basin from about 45,000 acre-feet to about 32,720 acre-feet would result in about 40 feet of additional drawdown in the central basin near the pumping wells, about 25 feet less than if withdrawals were not reduced. The combination of reducing and redistributing the cumulative withdrawals to the three existing and two proposed new wells results in about 40 feet of additional drawdown in the central and southern parts of the basin and about 70 feet in the northern part of the basin. These results show that reducing and redistributing the groundwater withdrawals would maintain the upper aquifer at greater than 50 percent of its predevelopment saturated thickness throughout the groundwater basin. The scenarios simulated for this study demonstrate how the calibrated model can be utilized to evaluate the hydrologic effects of different water-management strategies.

California↗

Integrated monitoring of hydrogeomorphic, vegetative, and edaphic conditions in riparian ecosystems of Great Basin National Park, Nevada

In semiarid regions such as the Great Basin, riparian areas function as oases of cooler and more stable microclimates, greater relative humidity, greater structural complexity, and a steady flow of water and nutrients relative to upland areas. These qualities make riparian areaʼs attractive not only to resident and migratory wildlife, but also to visitors in recreation areas such as Great Basin National Park in the Snake Range, east-central Nevada. To expand upon the system of ten permanent plots sampled in 1992 (Smith et al. 1994) and 2001 (Beever et al. in press ), we established a collection of 31 cross-sectional transects of 50-m width across the mainstems of Strawberry, Lehman, Baker, and Snake creeks. Our aims in this research were threefold: a) map riparian vegetative communities in greater detail than had been done by past efforts; b) provide a monitoring baseline of hydrogeomorphology; structure, composition, and function of upland- and riparianassociated vegetation; and edaphic properties potentially sensitive to management; and c) test whether instream conditions or physiographic variables predicted vegetation patterns across the four target streams. In each of the four watersheds, we performed walking transects from the lower-elevation boundary of the park along creek mainstems to a point well above the point at which vehicle access stopped. In these transects, we ranked, by cover, the riparian and upland woody species on each side of the creek, in 0.32-km segments. These walking transects also facilitated selection of a suite of cross-sectional transects that might serve as an early-warning signal of change for natural (e.g., aggradative) and anthropogenic changes (e.g., due to visitor impacts or climate change). At each cross-sectional transect, we used several methods: a) measurement of the number, approximate volume, and total length of instream logs greater than 10 cm in diameter that were within 5 m up- or downstream of the transect; b) counts of pebbles by size class, following Wolman (1954); c) line-point intercepts, which provided various measures of percent cover; d) gap-intercept transects, following Herrick et al. (in press), to measure susceptibility of uplands to erosion by wind or water; e) 1-m 2 quadrats, to obtain frequency of woody species; f) nested-frequency plots, to measure frequency of all plant species in quadrats of varying size; g) a field-based soil aggregate stability test following Herrick et al. (2001); and h) an impact penetrometer, to measure penetration resistance of soil horizons. We used species-accumulation curves to assess the ability of our methods to detect the majority of plant species at sites, using the most species-rich and species-poor sites as illustrations. We compared characteristics of hydrogeomorphic valley types (designated by Frissell and Liss 1993), vegetation types, and creeks individually and, using multivariate analyses for the first two ʻtypes,ʼ simultaneously. For the latter, using both the nested-frequency and 1-m 2 frequency data, we first used nonmetric multidimensional scaling (NMS) to assess relationships of plant communities among sites. Secondly, we used multi-response permutation procedures (MRPP) to test whether plant-community differences existed among either hydrogeomophic valley types or vegetation types. To increase the value of these comparisons for management, we used indicator species analyses to quantify the indicator value of each individual plant species for separating groups. In contrast to the more incised riparian channels of central Nevada, we observed knickzones, downcutting, and incision only rarely and usually with limited extent in the walking surveys. Downcutting occurred most frequently and extensively in Strawberry and Snake creeks, due in part to their more erodible soils. According to a hydrogeomorphologist with extensive experience in Great Basin riparian systems, the sediment-delivery and hydrologic systems appeared relatively undisturbed in most reaches, with respect to grazing animals and other types of anthropogenic alteration. Site elevation of the 31 transects ranged from 1,950-2,987 m, and stream slope (i.e., gradient) was relatively steep (mean = 9.3%, range 3-16%). Strawberry Creek averaged the lowest maximum water depth, and correspondingly had greatest width/depth ratios. Baker Creek sites averaged the smallest amount of tree-canopy gaps, whereas Snake Creek sites on average had the largest proportion of gaps in understory vegetation. Sites in terrace-bound valley types averaged the lowest slope in the channel as well as the least cover of trees, litter, and vegetation overall, whereas alluviated, boulder-bed canyon sites averaged the greatest widths of the active channel. Sites in Lehman Creek averaged nearly twice as much coarse woody debris as sites from any other creek, whereas Baker Creek sites averaged greatest tree cover (mean = 67%, range 40 – 96%) and species richness (mean = 17.3 species). Multivariate ordinations suggested that sites in leveed outwash valleys and alluvial-fan-influenced valleys had the greatest inter-site heterogeneity in plant composition, whereas sites in incised moraine-filled valleys appeared most homogeneous. Differences among homogeneity of sites within vegetation types were less pronounced, but sites dominated by either aspen and Woodsʼ rose or narrow-leaved cottonwood had the most similar plant communities among sites of the same vegetation type. A number of species were faithful indicators of various valley and vegetation types, using either set of plant-frequency data. We estimate that all 31 sites could be subsequently re-sampled in 14-18 field days by individuals possessing familiarity of the riparian flora of the southern Snake Range. As with any research, monitoring-focused investigations must balance the concerns for number of ecosystem attributes measured, extensiveness in time and space of sampling periods and locations, and the time and cost of sampling.

Nevada↗

Linkage of the Soil and Water Assessment Tool and the Texas Water Availability Model to simulate the effects of brush management on monthly storage of Canyon Lake, south-central Texas, 1995-2010

The U.S. Geological Survey (USGS), in cooperation with the Texas State Soil and Water Conservation Board, developed and applied an approach to create a linkage between the published upper Guadalupe River Soil Water Assessment Tool (SWAT) brush-management (ashe juniper [ Juniperus ashei ]) model and the full authorization version Guadalupe River Water Availability Model (WAM). The SWAT model was published by the USGS, and the Guadalupe River WAM is available from the Texas Commission on Environmental Quality. The upper Guadalupe River watershed is a substantial component of the Guadalupe River WAM. This report serves in part as documentation of a proof of concept on the feasibility of linking these two water-resources planning models for the purpose of simulating possible increases in water storage in Canyon Lake as a result of different brush-management scenarios. The SWAT-WAM linkage for the upper Guadalupe River is documented with a principal objective to evaluate the distributional characteristics of the monthly water storage of Canyon Lake during selected drought conditions. Focus is on the relative evaluation of select scenarios of large-scale or “extensive” brush management within the upper Guadalupe River watershed. There are six SWAT simulations for the upper Guadalupe River watershed that include a baseline (0-percent management of treatable ashe juniper, the baseline scenario from a previous study in which no percentage of ashe juniper is numerically replaced with grassland) along with five scenarios (extensions of SWAT simulations from a previous study) of 20-, 40-, 60-, 80-, and 100-percent random (numerical) replacement of treatable ashe juniper with grasslands throughout the upper Guadalupe River watershed in south-central Texas. SWAT is a process-based, semidistributed, water-balance model designed to predict the effects of landscape management decisions on water yields. A watershed is subdivided into subbasins, and each subbasin is associated with a single reach on the stream network. In general a WAM, such as the Guadalupe River WAM, provides analysis of generalized water rights in a river and reservoir framework. A WAM accommodates hydrology and water usage through several input files containing water rights, watershed parameters, and naturalized streamflow time series. A WAM is generalized for application to rivers and reservoir systems, and input datasets are uniquely developed for a river basin of concern. The extractions of SWAT output for the five extensive brush-management and baseline scenarios were offset by –21 years and, in general, the results were then mapped to the WAM input-flow file. The offset of –21 years was chosen arbitrarily for technical reasons and means that the period of monthly record 1995–2010 of the upper Guadalupe River SWAT became the synthetic period of monthly record 1974–89, hereinafter 1974–89 (synthetic) period, of the Guadalupe River WAM. The relative (between scenario to baseline) effects of extensive brush-management scenarios by using the SWAT-WAM linkage were evaluated, and two critical intermediate results were total inflow to Canyon Lake from 1995 to 2010 and the monthly storage of Canyon Lake from 1974 to 1989 (synthetic). The first quartile or lower 25th percentile of monthly storage of Canyon Lake for the baseline scenario is 381,000 acre-feet (acre-ft) for the hereinafter 1974–89 (synthetic) period. This lower quartile was chosen for analysis for two critical purposes. First, Canyon Lake is managed with a conservation pool of about 386,200 acre-ft capacity (as recognized by the WAM) and is at or near conservation capacity about 50 percent or more of the time; further, there is intrinsic data censoring that occurs for the monthly storage distribution because Canyon Lake is at or near conservation pool elevation the majority of the time. This intrinsic censoring has the effect of creating a bounded distribution with a left or low-volume tail. Statistical assessment of the brush-management scenarios beginning with the 381,000 acre-ft censoring threshold provides readily interpretable results. Second, the quantification of brush management during periods lacking abundant rainfall, which were defined in this study as months for which Canyon Lake storage was below the 25th percentile for the simulation period, are of substantial interest to water-resource managers and stakeholders in the context of water-supply enhancement. A statistical assessment of the SWAT-WAM linkage for the low-volume tail of the distribution of monthly storage of Canyon Lake is the focus of analysis and interpretation. Drought periods for the analysis are defined as the months (consecutive or not) during which Canyon Lake is below the 25th percentile of storage (381,000 acre-ft) for the baseline scenario. Such months are referred to as being within the “Drought Quartile.” The Drought Quartile is a conceptual and heuristically determined waypoint for the analysis and is not related to any administrative definition of drought by stakeholders or policy makers. The five scenarios and the baseline scenario simulated in the upper Guadalupe River SWAT were all passed through the Guadalupe River WAM by the SWAT-WAM linkage described in this report. A comparison of the mean increase per month in reservoir storage for Canyon Lake conditioned for the Drought Quartile was made. For each of the five brush-management and baseline scenarios, the months with storage below 381,000 acre-ft were extracted. The mean monthly storages during the Drought Quartile were computed for each of the five scenarios and the baseline scenario. The mean of the baseline scenario was 376,458 acre-ft and subsequently was subtracted from the mean monthly storage during the Drought Quartile for each of the five scenarios. The mean monthly offset storages of Canyon Lake during the Drought Quartile were 110 acre-ft (20 percent); 448 acre-ft (40 percent); 754 acre-ft (60 percent); 1,080 acre-ft (80 percent); and 1,090 acre-ft (100 percent). A particular mean was interpreted as follows: the value of 754 acre-ft for the 60-percent brush-management scenario implies that, on average, this scenario indicates an additional 754 acre-ft per month of storage in Canyon Lake relative to the baseline during the Drought Quartile. All of the five scenarios resulted in an increase on average to water supply relative to the baseline scenario during the Drought Quartile through the SWAT-WAM linkage.

Texas↗

Storage capacity and sedimentation characteristics of the San Antonio Reservoir, California, 2018

The San Antonio Reservoir is a large water storage facility in Alameda County, California, and is a major component of the Hetch Hetchy Regional Water System (RWS). The RWS is a water-supply system owned and operated by the San Francisco Public Utilities Commission (SFPUC) and provides water for about 2.7 million people in the San Francisco, Santa Clara, Alameda, and San Mateo Counties. The San Antonio Reservoir is one of two RWS reservoirs in Alameda County and the third largest of the RWS reservoirs in the San Francisco Bay Area. The reservoir was formed by the James H. Turner Dam, which was completed in 1965. At the time of construction, the reservoir was estimated to have 50,500 acre-feet (acre-ft) of storage capacity. That early estimate was based on a 1963 pre-construction topographic map, which was drawn from aerial photographs. The capacity of the reservoir was later surveyed in 1994 and 2000. These two later surveys did not include the upper 18 feet (ft) of the reservoir, which represents roughly 30 percent of the overall storage volume. To determine the storage capacity and provide updated stage-capacity curves up to the spillway, the U.S. Geological Survey, in cooperation with the SFPUC, surveyed the bathymetry and shoreline of the reservoir in April 2018. The bathymetric survey was performed by making depth soundings using a boat-mounted, multibeam echosounder. At the time of the survey, the water level was between 13 and 14 ft below the spillway elevation. To measure capacity between the water line up to the spillway elevation, topography along most of the shoreline was surveyed from the boat using a terrestrial Light Detection and Ranging (LiDAR) scanner and in other areas by using ground-survey techniques. Location during bathymetric and topographic data collection was determined using a Global Navigation Satellite System-Real Time Network system. Vertical profiles of sound speed were collected periodically. The sound-speed profiles were used to spatially and temporally adjust the sound-speed calculations used to determine depth from the soundings. Approximately 125 kilometers (78 miles) of transects with a total of about 560 million depth soundings and topographic LiDAR points were collected (about 160 per square meter). In addition, approximately 500 topographic survey points were collected in shallow, wadable areas and on land near the upper reservoir area using a Global Navigation Satellite System receiver attached to a fixed length survey rod. Depth soundings, terrestrial LiDAR points, topographic survey points, and a digitized shoreline were merged and interpolated to generate a digital elevation model (DEM) of the reservoir. Gridded elevation data extracted from the DEM were then tabulated to determine total reservoir capacity and create reservoir stage-surface area and stage-storage capacity tables. Results of the reservoir capacity analysis indicated that the reservoir has 53,266 (plus or minus 140) acre-ft of storage capacity, which is an increase of 2,766 acre-ft (or 5.5 percent) greater than the original 1965 estimate; the increase is likely due to improved survey methods. Also, at the time of this 2018 survey, Intake #1 (the lowest intake) was not in operation. Intake #1 is estimated to be buried approximately 10 ft below the bed, whereas Intake #2 is about 20 ft above the bed. There are five intakes at different elevation levels; however, when consecutive lower intakes become inoperable due to sedimentation, the live storage capacity (capacity available for use) is reduced. At the time of this survey, the remaining live storage (above Intake #2) was approximately 52,363 acre-ft. The 2018 stage-capacity curve was compared to the original 1965 stage-capacity curve. Although overall, the changes indicate an increase in storage capacity, the change in volume at 372.7 ft North American Vertical Datum of 1988 (370 ft National Geodetic Vertical Datum of 1929, NGVD 29) shows a decrease of 733 acre-ft (the elevation of 370 ft NGVD 29 was used because it is the lowest elevation available for the 1965 stage-capacity curves). This finding agrees with the observed accumulation of sediment over Intake #1. That volume was converted to an annual sediment yield of 0.35 acre-ft per square mile (or 165 cubic meters per square kilometer), which is of the same order of magnitude as that found in other watersheds for the Coast Ranges in California. A decrease of 733 acre-ft between 1965 and 2018 thus represents a loss of 1.5 percent of the overall storage capacity in the reservoir. The updated stage-surface area and stage-capacity tables provided in this report and online ( https://doi.org/10.5066/P9KC9DU8 ) can be used by the SFPUC to improve reservoir operations and serve as an accurate baseline to monitor bathymetric changes in the future.

California↗

Assessment of channel morphology, hydraulics, and bedload transport along the Siletz River, western Oregon

Significant Findings Chinook salmon ( Oncorhynchus tshawytscha ) and Pacific lamprey ( Entosphenus tridentatus ) are native, anadromous fish species in the Siletz River Basin, western Oregon, that face many threats to their survival in freshwater and the ocean. The Confederated Tribes of Siletz Indians of Oregon seek to mitigate freshwater threats to Chinook salmon and Pacific lamprey, where possible, with habitat conservation and restoration efforts. This study was conducted to assist the Confederated Tribes of Siletz Indians of Oregon in documenting and understanding the hydrogeomorphic processes shaping present-day habitat conditions and assessing future habitat implications for Chinook salmon and Pacific lamprey along the main-stem Siletz River. As such, this study focused on understanding geomorphic processes and patterns of channel change, including lateral and vertical adjustments in channel position and changes in bed-material sediment (sands, gravels, and cobbles that mantle the channel bed), which collectively determine overall patterns of channel morphology and fluvial habitats. Objective One was to evaluate lateral changes in channel position, vertical changes in bed elevation, and longitudinal patterns in bed-material particle size along the Siletz River using detailed channel maps developed from aerial photographs collected from 1939 to 2016, long-term records of stage and discharge collected by the U.S. Geological Survey (USGS) near the City of Siletz, and sediment particle size data. Objective Two was to assess hydraulic conditions using one- and two-dimensional hydraulic models and transport capacity of bed-material sediment using bedload transport models and sediment particle size data for a range of discharge conditions. Objective Three was to identify potential burrowing habitat for lamprey larvae (PBH) along the Siletz River network and provide insights in local factors influencing PBH along the main-stem Siletz River. The overall findings are synthesized to describe habitat implications for Chinook salmon and Pacific lamprey under present-day and future conditions. Results of Objective One, an evaluation of changes in channel position and bed elevations and longitudinal patterns in bed-material particle size along the Siletz River, include the following From 1939 to 2016, channel planform and the mapped area of gravel bars did not change considerably along the 97.2-kilometers (km) of Siletz River between Elk Creek and Millport Slough, except for in short sections generally less than 1 km long. This inherent lateral channel stability results from the resistant bedrock and terraces that bound most of the channel and limit lateral changes in channel position. Intermittent sections along the study reach where the active channel widened at channel bends displayed noticeable planform changes and increases in mapped bar area. From water year (WY; a 12-month period from October 1 through September 30 and named for the year in which it ends) 1906 to 2021, changes in the stage-discharge relation interpreted as rising and lowering channel bed levels were observed at the USGS streamgage 14305500 (Siletz River at Siletz, OR) in response to floods (such as high flows in February 1996 and January 2002 that exceeded 0.1 and 0.667 annual exceedance probabilities [AEP] events, respectively). However, the rating curve representing the stage-discharge relation did not change in response to high-magnitude floods between 2007 and 2021. Along the approximately 54-km of the Siletz River between Moonshine Park and the Bulls Bag area, surficial particle distributions varied considerably between sampling sites in response to changes in channel width and gradient, sediment inputs from tributaries, and basin geology. Despite this variability, median particle sizes tended to decrease in size in the downstream direction over the 54 km. Results of Objective Two, an evaluation of hydraulic and bedload transport conditions along the Siletz River, include the following The most substantial increases in maximum and mean water velocity and bedload sediment transport capacity occurred at events between the 0.995 and 0.50 AEPs. Events of these magnitudes occur approximately every 1–2 years. Smaller events (0.995 AEP) are generally contained by the banks of the main channel, whereas larger events (0.50 AEP) generally spill over the banks and inundate high-elevation bar and low-elevation floodplain surfaces. Multiple smaller floods within a WY that exceed the 0.995 and 0.50 AEPs (such as in WY 1996) can transport as much or more bed-material sediment than a single, higher magnitude event (such as the maximum event recorded in WY 2000 with an estimated AEP of around 0.002). Bedload transport capacity generally exceeds sediment supply (greater than 2 millimeters [mm]) for most of the study area from Wildcat Creek to the City of Siletz as evidenced by substantial in-channel bedrock. Despite overall conditions of limited sediment supply, transport capacity still varies considerably within and between years with discharge magnitude and spatially in relation to local hydraulics imposed by bedrock, channel morphology, and human infrastructure (such as bridges). Results of Objective Three, an analysis of PBH for lamprey larvae, include the following About 28 percent of the Siletz River network meets the mean annual suspended sediment loads and channel slope criteria for PBH. Along the main-stem Siletz River, in-channel bedrock outcrops and high transport capacity are expected to further constrain PBH. Together, these results suggest that most of the Siletz River between Wildcat Creek and the City of Siletz has had only modest vertical and lateral change between the 1930s and 2010s because of the bedrock in and along the main channel and the river’s relatively high transport capacity relative to bed-material sediment supply. However, localized sections of the Siletz River where the active channel widens, particularly at channel bends, exhibited some change in channel planform and the locations and area of gravel bars. In the future, moderate increases in autumn-winter discharge may not result in substantial changes in coarse gravel bars along the Siletz River but may result in selective transport of finer bed-material sediment (gravel, sands, and silts) that provide spawning habitats for Chinook salmon and Pacific lamprey and burrowing habitats for lamprey larvae. Assuming no substantial changes in bed-material sediment supply, increased bedload transport capacity may cause frequent entrainment of lamprey larvae that are burrowed in coarse sand deposits, suspension and downstream transport of salmon eggs incubating in gravels, and reductions in the areas of spawning gravels for Chinook salmon and Pacific lamprey. Exact implications of current and future discharge conditions for these species along the Siletz River depends on many factors, including sediment supply, local hydraulics, and the timing of flood events relative to fish life stages.

Oregon↗

Depth of groundwater used for drinking-water supplies in the United States

Groundwater supplies 35 percent of drinking water in the United States. Mapping the quantity and quality of groundwater at the depths used for potable supplies requires an understanding of locational variation in the characteristics of drinking-water wells (depth and open interval). Typical depths of domestic- and public-drinking-water supply wells vary by and within aquifer across the United States. The depths to the top and bottom of the zones from which drinking water is withdrawn are important predictor variables in regional- and national-scale statistical water models, but spatially extensive maps of the depths to drinking-water-supply sources are not consistently available in modeled regions. Therefore, it was necessary to generate a set of grids representing surfaces of the approximate common depth and length of open intervals in the wells from which water is withdrawn for domestic- and public-drinking-water supply (withdrawal zones) within the conterminous United States. Well data (about 7.6 million records) were compiled from several sources, including the U.S. Geological Survey’s National Water Information System (600,922 records), the U.S. Environmental Protection Agency’s Safe Drinking Water Information System dataset (66,540 records, primarily public-supply wells), a groundwater ambient monitoring dataset (31,448 records, primarily domestic-supply wells), individual State data (6,096,503 records), a national brackish aquifer study (96,885 records), and a glacial aquifer study (729,564 records). Fifty-seven principal aquifers and 65 secondary hydrogeologic regions have been designated in the conterminous United States. The principal aquifers and secondary hydrogeologic regions vary in depth, thickness, lithology, and transmissivity characteristics. Some principal aquifers underlie secondary hydrogeologic regions, and may in turn be overlain by glacial sediment or basin and valley fill aquifers, which may also be used as drinking-water sources. The principal aquifer and secondary hydrogeologic region polygons were merged with overlying sediment polygons, where present, including glacial sediment, coarse glacial sediment, and stream valley alluvium (alluvium) polygons, to generate unique hydrogeologic settings across the conterminous United States. A total of 288 distinct hydrogeologic settings resulted from the merging of principal aquifer, secondary hydrogeologic region, glacial sediment, coarse glacial sediment, and alluvium polygons. Each well was assigned to a hydrogeologic setting on the basis of location. Hydrogeologic setting well groupings were used to guide calculations of the median value for well depth and depth to and length of open intervals across the hydrogeologic setting. Where well data were sparse or missing, wells from hydrogeologic settings with similar well construction properties, geology, physiography, and topography were grouped and used to calculate the moving median depth (if less than five wells in a 100-kilometer [62.1-mile] radius) and to estimate open interval length (if not available within hydrogeologic setting). Grids were generated to represent what might be considered as the “typical” or “median” domestic- and public-supply well in an area. The well properties are defined with moving median grids of top depth, bottom depth, and length of open interval at a 1-square-kilometer (0.38-square-mile) grid cell scale. Median depths and open intervals of domestic- and public-supply wells varied by lithology of the hydrogeologic setting and overlying sediment. Overall, the median depths were 142 feet (43.3 meters) for all domestic-supply wells and 202 feet (61.6 meters) for all public-supply wells. The median open intervals were 21 feet (6.4 meters) for domestic-supply wells and 49 feet (14.9 meters) for public-supply wells. The shallowest median bottom open interval depths for domestic-supply wells were in the secondary hydrogeologic regions with coarse glacial sediment, which suggests that the wells are most commonly completed in the permeable coarse glacial sediment and not in the underlying secondary hydrogeologic region. Public-supply wells were completed at relatively shallow median depths when drilled in permeable sediment that overlie secondary hydrogeologic regions. When public-supply wells were completed in principal aquifers, the median depths were typically greater than wells completed in secondary hydrogeologic regions. Well data used in this study were limited to those available from national or State digital databases. Several quality-assurance checks were performed during data compilation, but a comprehensive quality assurance inspection for each of the data sources was outside the scope of this study. Grids defining typical open intervals in domestic- and public-supply wells are presented. Although there are many places where multiple aquifers are stacked, these results correspond primarily to the aquifer with the highest documented number of wells for each use.

Scientific Investigations Report↗

Hydrogeology of sand and gravel aquifers in the Owasco Inlet watershed, Cayuga and Tompkins Counties, New York

This study is a continuation of a series of hydrogeologic appraisals that have been conducted since 1980, as part of a cooperative, long-term, detailed aquifer mapping program by the U.S. Geological Survey and the New York State Department of Environmental Conservation. These appraisals provide a foundation for wellhead protection programs, water-resource management and planning decisions, and groundwater remediation in upstate New York. The Owasco Inlet watershed drains north directly to Owasco Lake, one of New York’s Finger Lakes. The watershed is similar in form to other watersheds of the Finger Lakes region of New York State and is characterized by broad, smooth uplands punctuated by the 18-mile-long Owasco Inlet valley and the secondary valleys of Decker, Dresserville, and Hemlock Creeks. All these streams occupy “through valleys”—the valleys are continuous into adjacent watersheds such that watershed divides are high points within the valley bottoms, rather than in uplands. Most glacial deposits in the watershed were deposited during the Valley Heads Readvance and subsequent retreat (with at least one minor readvance). Estimates are that the Valley Heads Readvance likely peaked in the Cayuga basin about 17,000 (calendar) years ago based on dates from western New York and the Mohawk Valley in east-central New York. It was the last major ice advance of the Pleistocene Epoch in this region. This readvance covered the entire Owasco watershed. Valley Heads ice in the Owasco watershed was part of the Cayuga ice lobe, which flowed parallel to the Owasco Inlet valley as far south as Moravia, but which became progressively more eastward south of Locke, effectively raking across the valley from the west. The Owasco Inlet watershed drains north, like other Finger Lakes watersheds, and during deglaciation, northward meltwater drainage was largely blocked by the ice, which resulted in development of a series of regional proglacial lakes in which fine lacustrine sediments predominated. Lacustrine sediments constitute much of the Owasco Inlet valley fill, but it has been pointed out that there are no obvious outlet channels exiting the Owasco Inlet watershed, and it has been proposed that there was northward drainage of meltwater into the ice within the larger, neighboring Cayuga Trough. The lacustrine deposits in the Owasco Inlet valley are commonly underlain by thinner coarse-grained stratified deposits (subaqueous fans and eskers) of variable sorting and permeability. An exception occurs at Groton, N.Y., where retreating ice paused long enough for coarse-grained sediments to fill much of the valley. Smaller, higher elevation valleys have a wider variety of glacial valley-fill deposits, ranging from fine lacustrine sediments to sand and gravel to till. Groundwater is the sole source of water supply in the area; glacial sand and gravel aquifers are the primary water source in the Owasco Inlet and Decker Creek valleys, and fractured bedrock aquifers typically supply domestic wells in the remaining valleys and upland areas. Municipal supplies tap a variety of aquifer types in the Owasco Inlet valley and nearby uplands. The hamlet of Locke taps the extensive confined aquifer beneath fine-grained lacustrine deposits in the valley. The Village of Moravia taps a semiconfined aquifer that overlies proglacial lake deposits and is partly confined by overlying recent lake deposits and alluvium. Withdrawals from this aquifer may also induce water from the Owasco Inlet into the aquifer. The Village of Groton draws from two aquifers: (1) a thin upland unconfined sand and gravel aquifer tapped by an infiltration gallery and (2) the local unconfined or semiconfined aquifer in the Owasco Inlet valley. Another aquifer with potential for municipal supply is the unconfined aquifer in the Decker Creek valley near Wilson Corners. The confined aquifer in the lower Dresserville Creek valley may have some water-resource potential, but it is largely untested. Unconfined aquifers are the most susceptible to contamination from activities at land surface directly above the aquifer because precipitation and subsequent recharge can transport contaminants directly to the water table. Adjacent upland areas can also contribute contaminants. Confined aquifers are less susceptible to contamination from overlying land surface areas because confining units largely prevent downward movement of water. Recharge occurs elsewhere at unconfined upvalley locations and along valley walls where alluvial fans, ice-contact deposits, or stream incision into the valley wall may provide pathways for downward movement of groundwater.

New York↗

Simulation of groundwater flow and analysis of the effects of water-management options in the North Platte Natural Resources District, Nebraska

The North Platte Natural Resources District (NPNRD) has been actively collecting data and studying groundwater resources because of concerns about the future availability of the highly inter-connected surface-water and groundwater resources. This report, prepared by the U.S. Geological Survey in cooperation with the North Platte Natural Resources District, describes a groundwater-flow model of the North Platte River valley from Bridgeport, Nebraska, extending west to 6 miles into Wyoming. The model was built to improve the understanding of the interaction of surface-water and groundwater resources, and as an optimization tool, the model is able to analyze the effects of water-management options on the simulated stream base flow of the North Platte River. The groundwater system and related sources and sinks of water were simulated using a newton formulation of the U.S. Geological Survey modular three-dimensional groundwater model, referred to as MODFLOW–NWT, which provided an improved ability to solve nonlinear unconfined aquifer simulations with wetting and drying of cells. Using previously published aquifer-base-altitude contours in conjunction with newer test-hole and geophysical data, a new base-of-aquifer altitude map was generated because of the strong effect of the aquifer-base topography on groundwater-flow direction and magnitude. The largest inflow to groundwater is recharge originating from water leaking from canals, which is much larger than recharge originating from infiltration of precipitation. The largest component of groundwater discharge from the study area is to the North Platte River and its tributaries, with smaller amounts of discharge to evapotranspiration and groundwater withdrawals for irrigation. Recharge from infiltration of precipitation was estimated with a daily soil-water-balance model. Annual recharge from canal seepage was estimated using available records from the Bureau of Reclamation and then modified with canal-seepage potentials estimated using geophysical data. Groundwater withdrawals were estimated using land-cover data, precipitation data, and published crop water-use data. For fields irrigated with surface water and groundwater, surface-water deliveries were subtracted from the estimated net irrigation requirement, and groundwater withdrawal was assumed to be equal to any demand unmet by surface water. The groundwater-flow model was calibrated to measured groundwater levels and stream base flows estimated using the base-flow index method. The model was calibrated through automated adjustments using statistical techniques through parameter estimation using the parameter estimation suite of software (PEST). PEST was used to adjust 273 parameters, grouped as hydraulic conductivity of the aquifer, spatial multipliers to recharge, temporal multipliers to recharge, and two specific recharge parameters. Base flow of the North Platte River at Bridgeport, Nebraska, streamgage near the eastern, downstream end of the model was one of the primary calibration targets. Simulated base flow reasonably matched estimated base flow for this streamgage during 1950–2008, with an average difference of 15 percent. Overall, 1950–2008 simulated base flow followed the trend of the estimated base flow reasonably well, in cases with generally increasing or decreasing base flow from the start of the simulation to the end. Simulated base flow also matched estimated base flow reasonably well for most of the North Platte River tributaries with estimated base flow. Average simulated groundwater budgets during 1989–2008 were nearly three times larger for irrigation seasons than for non-irrigation seasons. The calibrated groundwater-flow model was used with the Groundwater-Management Process for the 2005 version of the U.S. Geological Survey modular three-dimensional groundwater model, MODFLOW–2005, to provide a tool for the NPNRD to better understand how water-management decisions could affect stream base flows of the North Platte River at Bridgeport, Nebr., streamgage in a future period from 2008 to 2019 under varying climatic conditions. The simulation-optimization model was constructed to analyze the maximum increase in simulated stream base flow that could be obtained with the minimum amount of reductions in groundwater withdrawals for irrigation. A second analysis extended the first to analyze the simulated base-flow benefit of groundwater withdrawals along with application of intentional recharge, that is, water from canals being released into rangeland areas with sandy soils. With optimized groundwater withdrawals and intentional recharge, the maximum simulated stream base flow was 15–23 cubic feet per second (ft 3 /s) greater than with no management at all, or 10–15 ft 3 /s larger than with managed groundwater withdrawals only. These results indicate not only the amount that simulated stream base flow can be increased by these management options, but also the locations where the management options provide the most or least benefit to the simulated stream base flow. For the analyses in this report, simulated base flow was best optimized by reductions in groundwater withdrawals north of the North Platte River and in the western half of the area. Intentional recharge sites selected by the optimization had a complex distribution but were more likely to be closer to the North Platte River or its tributaries. Future users of the simulation-optimization model will be able to modify the input files as to type, location, and timing of constraints, decision variables of groundwater withdrawals by zone, and other variables to explore other feasible management scenarios that may yield different increases in simulated future base flow of the North Platte River.

Nebraska↗

Hydrogeologic characterization, groundwater chemistry, and vulnerability assessment, Ute Mountain Ute Reservation, Colorado and Utah

The U.S. Geological Survey, in cooperation with the Ute Mountain Ute Tribe (UMUT), initiated a study in 2016 to increase understanding of the hydrogeology and chemistry of groundwater within select areas of the Ute Mountain Ute Reservation (UMUR) in Colorado and Utah, identify vulnerabilities to the system and other natural resources, and outline information needs to aid in the understanding and protection of groundwater resources. The results presented for this study can be used to support the UMUT’s goal of protecting their vital groundwater resources on the UMUR. Hydrogeologic conditions were characterized for the surficial aquifer contained in Quaternary-age unconsolidated surficial deposits and the Dakota aquifer contained in the Cretaceous-age Dakota Sandstone. In the surficial aquifer, median depth to water ranges from about 5.4 to 17.2 feet below land surface in the Farm and Ranch Enterprise area and 11 to 34 feet below land surface in the Towaoc area, and the water table slopes generally southwest or south. A map of depth to the top of the Dakota Sandstone was constructed from existing well data. Depths range from zero in outcrop areas to more than 3,000 feet below land surface on mesas in the southeastern part of the UMUR. Groundwater-chemistry data were collected by the UMUT from 13 springs and 31 wells from 1996 through 2017. Specific conductance was much lower for samples from springs than from wells; median values were 512 and 6,024 microsiemens per centimeter at 25 degrees Celsius, respectively. Spring samples were well oxygenated. A few well samples were anoxic (dissolved oxygen concentrations less than 0.5 milligrams per liter [mg/L]), indicating reducing conditions in the aquifer. About 75 percent of spring samples had fresh water (total dissolved solids concentrations less than 1,000 mg/L), and about 85 percent of well samples had brackish or highly saline water (total dissolved solids concentrations greater than 1,000 mg/L). Water type for springs on the Ute Mountains was calcium bicarbonate. Lower-altitude springs had a calcium-sulfate water type. Most well samples had sodium as the dominant cation, and sulfate, bicarbonate, and chloride as the dominant anions. Fluoride concentrations in about 45 percent of well samples were greater than an agricultural-use standard of 2 mg/L. Nitrate plus nitrite concentrations in most spring and well samples were less than about 1.6 mg/L per liter. Concentrations in samples from wells in the irrigated agricultural area were elevated; the maximum concentration was 78.5 mg/L. About one-half of the trace-element samples had concentrations that were less than laboratory reporting limits. Only aluminum, arsenic, and selenium in spring samples, and boron and selenium in well samples, were detected at concentrations greater than surface-water standards or water-quality standards for agricultural use of groundwater. Only three organic compounds, the pesticides alachlor and atrazine and the volatile organic compound di(2-ethylhexyl) phthalate, were detected in well samples. The Escherichia coli bacteria was detected in 47 and 23 percent of samples from wells and springs, respectively. The E. coli detections included samples from three culturally significant springs, which did not meet the UMUT cultural-use standard of total absence of E. coli. Tritium and carbon-14 were the primary environmental tracers used for interpreting groundwater ages for Lopez 2 Spring and five wells (AP–1, 5000 Block, Cottonwood Spring, Goodknight, and SE Toe). Water from the AP–1 well contained a mixture of pre- and post-1950s recharge. Tritium and carbon-14 recharge ages for Lopez 2 Spring (post-1950s in age), Goodknight and SE Toe wells (pre-1950s in age), and Cottonwood Spring well (primarily pre-1950s in age) are supported by helium-4 data. The helium-4 data for the 5000 Block well are inconsistent with the tritium and carbon-14 age of pre-1950s recharge because of interference caused by high methane concentrations in the water. Springs and surficial deposits are more vulnerable to contamination from anthropogenic chemicals than deeper bedrock wells. Bedrock aquifers are vulnerable in areas where the geologic formations containing the aquifers are exposed at the land surface. Groundwater in deep bedrock aquifers is likely thousands of years old and is not currently affected by present-day land uses. Both shallow and deep groundwater are vulnerable to naturally occurring salts and minerals, such as of total dissolved solids, major ions, nitrate, and trace elements. Effects of a changing climate on water resources and other ecological characteristics of the UMUR could include changes in evapotranspiration, a decrease in snowpack, decreased aquifer recharge and flow of springs, a decrease in soil moisture, and increased occurrence of wildfires and forest mortality. Of particular interest for the UMUT are possible effects of a changing climate on medicinal and culturally important plants and springs Several information needs were identified during this study that would aid in the understanding and protection of groundwater resources on the UMUR. These include well-completion information for bedrock wells, the collection of environmental tracer data at additional wells, the addition of methane and hydrocarbon analysis to well sampling plans, and the resampling of springs and wells that were last sampled in 2002 or earlier.

Colorado↗

Borehole geophysical investigation of a formerly used defense site, Machiasport, Maine, 2003-2006

The U.S. Geological Survey, in cooperation with the U.S. Army Corps of Engineers, collected borehole geophysical logs in 18 boreholes and interpreted the data along with logs from 19 additional boreholes as part of an ongoing, collaborative investigation at three environmental restoration sites in Machiasport, Maine. These sites, located on hilltops overlooking the seacoast, formerly were used for military defense. At each of the sites, chlorinated solvents, used as part of defense-site operations, have contaminated the fractured-rock aquifer. Borehole geophysical techniques and hydraulic methods were used to characterize bedrock lithology, fractures, and hydraulic properties. In addition, each geophysical method was evaluated for effectiveness for site characterization and for potential application for further aquifer characterization and (or) evaluation of remediation efforts. Results of borehole geophysical logging indicate the subsurface is highly fractured, metavolcanic, intrusive, metasedimentary bedrock. Selected geophysical logs were cross-plotted to assess correlations between rock properties. These plots included combinations of gamma, acoustic reflectivity, electromagnetic induction conductivity, normal resistivity, and single-point resistance. The combined use of acoustic televiewer (ATV) imaging and natural gamma logs proved to be effective for delineating rock types. Each of the rock units in the study area could be mapped in the boreholes, on the basis of the gamma and ATV reflectivity signatures. The gamma and mean ATV reflectivity data were used along with the other geophysical logs for an integrated interpretation, yielding a determination of quartz monzonite, rhyolite, metasedimentary units, or diabase/gabbro rock types. The interpretation of rock types on the basis of the geophysical logs compared well to drilling logs and geologic mapping. These results may be helpful for refining the geologic framework at depth. A stereoplot of all fractures intersecting the boreholes indicates numerous fractures, a high proportion of steeply dipping fractures, and considerable variation in fracture orientation. Low-dip-angle fractures associated with unloading and exfoliation are also present, especially at a depth of less than 100 feet below the top of casing. These sub-horizontal fractures help to connect the steeply dipping fractures, making this a highly connected fracture network. The high variability in the fracture orientations also increases the connectivity of the fracture network. A preliminary comparison of all fracture data from all the boreholes suggests fracturing decreases with depth. Because all the boreholes were not drilled to the same depth, however, there is a clear sampling bias. Hence, the deepest boreholes are analyzed separately for fracture density. For the deepest boreholes in the study, the intensity of fracturing does not decline significantly with depth. It is possible the fractures observed in these boreholes become progressively tighter or closed with depth, but this is difficult to verify with the borehole methods used in this investigation. The fact that there are more sealed fractures at depth (observed in optical televiewer logs in some of the boreholes) may indicate less opening of the sealed fractures, less water moving through the rock, and less weathering of the fracture infilling minerals. Although the fracture orientation remained fairly constant with depth, differences in the fracture patterns for the three restoration sites indicate the orientation of fractures varies across the study area. The fractures in boreholes on Miller Mountain predominantly strike northwest-southeast, and to a lesser degree they strike northeast. The fractures on or near the summit of Howard Mountain strike predominantly east-west and dip north and south, and the fractures near the Transmitter Site strike northeast-southwest and dip northwest and southeast. The fracture populations for the boreholes on or near the summit of Howard Mountain show more variation than at the other two sites. This variation may be related to the proximity of the fault, which is northeast of the summit of Howard Mountain. In a side-by-side comparison of stereoplots from selected boreholes, there was no clear correspondence between fracture orientation and proximity to the fault. There is, however, a difference in the total populations of fractures for the boreholes on or near the summit of Howard Mountain and the boreholes near the Transmitter Site. Further to the southwest and further away from the fault, the fractures at the Transmitter Site predominantly strike northeast-southwest and northwest-southeast.Heat-pulse flowmeter (HPFM) logging was used to identify transmissive fractures and to estimate the hydraulic properties along the boreholes. Ambient downflow was measured in 13 boreholes and ambient upflow was measured in 9 boreholes. In nine other bedrock boreholes, the HPFM did not detect measurable vertical flow. The observed direction of vertical flow in the boreholes generally was consistent with the conceptual flow model of downward movement in recharge locations and upward flow in discharge locations or at breaks in the slope of land surface. Under low-rate pumping or injection rates [0.25 to 1 gallon per minute (gal/min)], one to three inflow zones were identified in each borehole. Two limitations of HPFM methods are (1) the HPFM can only identify zones within 1.5 to 2 orders of magnitude of the most transmissive zone in each borehole, and (2) the HPFM cannot detect flow rates less than 0.010 + or - 0.005 gal/min, which corresponds to a transmissivity of about 1 foot squared per day (ft2/d). Consequently, the HPFM is considered an effective tool for identifying the most transmissive fractures in a borehole, down to its detection level. Transmissivities below that cut-off must be measured with another method, such as packer testing or fluid-replacement logging. Where sufficient water-level and flowmeter data were available, HPFM results were numerically modeled. For each borehole model, the fracture location and measured flow rates were specified, and the head and transmissivity of each fracture zone were adjusted until a model fit was achieved with the interpreted ambient and stressed flow profiles. The transmissivities calculated by this method are similar to the results of an open-hole slug test; with the added information from the flowmeter, however, the head and transmissivity of discrete zones also can be determined. The discrete-interval transmissivities ranged from 0.16 to 330 ft2/d. The flowmeter-derived open-hole transmissivity, which is the combined total of each of the transmissive zones, ranged from 1 to 511 ft2/d. The whole-well open-hole transmissivity values determined with HPFM methods were compared to the results of open-hole hydraulic tests. Despite the fact that the flowmeter-derived transmissivities consistently were lower than the estimates derived from open-hole hydraulic tests alone, the correlation was very strong (with a coefficient of determination, R2, of 0.9866), indicating the HPFM method provides a reasonable estimate of transmissivities for the most transmissive fractures in the borehole. Geologic framework, fracture characterization, and estimates of hydraulic properties were interpreted together to characterize the fracture network. The data and interpretation presented in this report should provide information useful for site investigators as the conceptual site groundwater flow model is refined. Collectively, the results and the conceptual site model are important for evaluating remediation options and planning or implementing the design of a well field and borehole completions that will be adequate for monitoring flow, remediation efforts, groundwater levels, and (or) water quality. Similar kinds of borehole geophysical logging (specifically the borehole imaging, gamma, fluid logs, and HPFM) should be conducted in any newly installed boreholes and integrated with interpretations of any nearby boreholes. If boreholes are installed close to existing or other new boreholes, cross-hole flowmeter surveys may be appropriate and may help characterize the aquifer properties and connections between the boreholes.

Maine↗

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

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

New Mexico↗

Assessing spatial variability of nutrients, phytoplankton, and related water-quality constituents in the California Sacramento–San Joaquin Delta at the landscape scale—Comparison of four (2018, 2020, 2021, 2022) spring high-resolution mapping surveys

Executive Summary This report summarizes results from boat-based, high-resolution water-quality mapping surveys completed before, during, and after upgrades to the EchoWater Resource Recovery Facility (EchoWater Facility), the regional wastewater facility for the City of Sacramento and surrounding areas, near Elk Grove, California. Surveys were completed in the tidal aquatic environments of the Sacramento–San Joaquin Delta (Delta) in spring (May or June) 2018, 2020, 2021, and 2022. In each survey, a suite of in situ sensors were used to continuously (one measurement per second) measure water-quality conditions, nutrients, phytoplankton abundance, and species composition. In addition to in situ data collection, discrete water samples were collected about every 2 miles while underway for determination of phosphate, ammonium, and nitrate concentration. The boat stopped at about 30 locations to collect discrete samples for a suite of additional analytes, including phytoplankton enumeration. The four surveys represent snapshots in time across different phases of the EchoWater Facility Biological Nutrient Reduction (BNR) upgrade. The May 2018 survey represents conditions before the upgrade. The second survey (June 2020) represents conditions after implementation of the Nitrifying Sidestream Treatment. The third survey (May 2021) was completed immediately after the completion of the BNR upgrade and represents a transitional period, and the final survey (May 2022) represents post-upgrade conditions. Relevant hydrologic and climatic context such as water-year type, X 2 position (the distance from the Golden Gate Bridge to the point upstream where bottom salinity is 2 parts per thousand; Jassby and others, 1995), water export to import ratio, and management actions like the Delta Cross Channel gate operations are presented for each survey so they may be considered in comparisons among surveys. Differences in water-quality parameters, like turbidity, temperature, salinity, pH, and dissolved oxygen (DO) improve understanding of nutrient cycling and phytoplankton dynamics. Because the Delta is a complex system, we divided the study area into hydrologic zones to better examine general trends and obtain a broadscale view of differences among the 4 study years. Results are presented for each survey and parameter using box plots to compare the different hydrologic zones. We also present each parameter using contour maps by survey to display gradients across the system. The most evident change to water quality in the Delta across surveys is related to the EchoWater Facility BNR upgrade, which included nitrification and denitrification processes. Through this upgrade, effluent ammonium (NH 4 + ) concentrations were reduced by more than 95 percent (from about 2,000 micromolars [μM] to below the reporting limit of 35 μM), and nitrate (NO 3 − ) concentrations increased from near zero to about 500 μM; therefore, the concentration of dissolved inorganic nitrogen (DIN; the sum of NH 4 + and NO 3 − ) in the effluent was reduced by about 75 percent between May 2018 and May 2022. The BNR upgrade resulted in a reduction in NH 4 + concentrations in aquatic habitats immediately below the facility, designated as the “north Delta tidal transition zone” (Bergamaschi and others, 2024), from about 30 μM pre-upgrade to near zero during the 2022 spring survey, whereas effluent NO 3 − increased from median concentrations of about 7 μM to about 15 μM. Because of the reduced effluent nitrogen loads and variability in Sacramento River nitrogen loads from upstream sources, DIN concentrations in the north Delta tidal transition zone decreased from a median of 53.3 μM in 2018 to 35.3 μM in 2020, 20.7 μM in 2021, and 11.3 μM in 2022 during the spring surveys. The changes in DIN concentration and form observed in the north Delta tidal transition zone after the EchoWater Facility upgrade extended downstream but were rapidly altered by hydrologic mixing, biogeochemical processes, and other nutrient source inputs. Most of the Delta indicated near-zero concentrations of NH 4 + 1 year after the completion of the EchoWater Facility upgrades represented by the 2022 survey. Exceptions to this finding were observed in the San Joaquin River near Stockton and in Suisun Bay, indicating there are NH 4 + inputs to these locations from other sources (for example, Stockton Regional Wastewater Control Facility and Central Contra Costs Sanitary District wastewater treatment plants or agricultural and urban runoff). Although there was an increase in NO 3 − concentrations in the north Delta tidal transition zone after the upgrade, increases in NO 3 − in other zones were not apparent, presumably because nitrification of effluent derived ammonium was no longer a source of NO 3 − . Concentrations of DIN in many Delta zones were lower in 2022 compared to 2018 and 2020, with concentrations near or below what is considered potentially nitrogen limiting conditions for phytoplankton growth in the North Delta tidal transition zone and the Cache Slough complex channel system. Unrelated to the EchoWater Facility upgrade, NO 3 − and therefore DIN concentrations increased in the San Joaquin River near Stockton and in adjacent water bodies by survey date (likely associated with increasing drought conditions). The Mokelumne River had low DIN concentrations, except in 2018 when the Delta Cross Channel was open, which allowed nutrient-rich Sacramento River water to flow into this section of the river. Data from these surveys also support the hypothesis that nutrient drawdown during phytoplankton blooms may create localized nitrogen limiting conditions. The BNR upgrade resulted in lower effluent phosphate (PO 4 3− ) concentrations, which lowered PO 4 3− concentrations in some zones of the Delta during the four spring surveys; however, PO 4 3− concentrations throughout the Delta remained above 0.3 μM, indicating that primary productivity was not limited by phosphorous availability. DIN and PO 4 3− decreased after the upgrade in many areas of the Delta, and the DIN to dissolved inorganic phosphorus (DIN:DIP) ratio remained similar to pre-upgrade conditions and was often below the Redfield Ratio of 16, indicating nitrogen is more likely to limit phytoplankton growth than phosphorous. Inputs of dissolved organic carbon (DOC) from the EchoWater Facility are a minor source of this constituent to the Delta, so the upgrade had little to no effect on DOC concentrations across the Delta. Because phytoplankton abundance and species composition in the Delta are shaped by multiple factors other than nutrients (for example, light availability, temperature, salinity, and predation), it is important to consider these factors (as well as long-term monitoring) in addition to the EchoWater Facility upgrade. Although phytoplankton populations were low across much of the Delta during the spring surveys, several localized phytoplankton blooms (defined here as greater than 15 micrograms per liter [μg/L] of chlorophyll) provide insight into conditions that may favor the growth of beneficial and harmful species.

Callifornia↗

Hydrologic assessment of the shallow groundwater flow system beneath the Shinnecock Nation tribal lands, Suffolk County, New York

Defining the distribution and flow of shallow groundwater beneath the Shinnecock Nation tribal lands in Suffolk County, New York, is a crucial first step in identifying sources of potential contamination to the surficial aquifer and coastal ecosystems. The surficial or water table aquifer beneath the tribal lands is the primary source of potable water supply for at least 6 percent of the households on the tribal lands. Oyster fisheries and other marine ecosystems are critical to the livelihood of many residents living on the tribal lands, but are susceptible to contamination from groundwater entering the embayment from the surficial aquifer. Contamination of the surficial aquifer from flooding during intense coastal storms, nutrient loading from fertilizers, and septic effluent have been identified as potential sources of human and ecological health concerns on tribal lands. The U.S. Geological Survey (USGS) facilitated the installation of 17 water table wells on and adjacent to the tribal lands during March 2014. These wells were combined with other existing wells to create a 32-well water table monitoring network that was used to assess local hydrologic conditions. Survey-grade, global-navigation-satellite systems provided centimeter-level accuracy for positioning wellhead surveys. Water levels were measured by the USGS during May (spring) and November (fall) 2014 to evaluate seasonal effects on the water table. Water level measurements were made at high and low tide during May 2014 to identify potential effects on the water table caused by changes in tidal stage (tidal flux) in Shinnecock Bay. Water level contour maps indicate that the surficial aquifer is recharged by precipitation and upgradient groundwater flow that moves from the recharge zone located generally beneath Sunrise Highway, to the discharge zone beneath the tribal lands, and eventually discharges into the embayment, tidal creeks, and estuaries that bound the tribal lands to the east, south, and west. Water levels in many of the wells in the network fluctuated in response to precipitation, upgradient groundwater flow, and tidal flux in Shinnecock Bay. Water level altitudes ranged from 6.66 to 0.47 feet (ft) above the North American Vertical Datum of 1988 during the spring measurement period, and from 5.25 to -0.24 ft (NAVD 88) during fall 2014. Historically, annual and seasonal precipitation seem to indicate long-term water level trends in an index well located in the town of Southampton, correlates with changes in storage in the upper glacial aquifer, but does not necessarily indicate water level extremes in the shallow groundwater system. To place the study period in perspective, calendar year 2014 was the 32d wettest year on record, with precipitation for the year totaling 48.1 inches, a 2.6-percent increase from the annual average (46.9 inches per year), based on 81 years of complete record at the National Oceanographic and Atmospheric Administration, National Weather Service cooperative meteorological station at Bridgehampton, New York. Estimated recharge to the water table beneath the tribal lands from precipitation for 2014 is 25.4 inches. Tidal flux caused water levels in wells to fluctuate from 0.30 to -0.24 ft during May 2014. Water levels in wells located north of Old Fort Pond and beneath the southernmost extent of the tribal lands were most influenced by tidal flux. During June 2014, hydrographs indicate that tidal flux influenced water levels by 0.48 ft in a well located near the southernmost extent of the tribal lands approximately 0.3 miles north of Shinnecock Bay, and was zero at a well located approximately 0.5 miles south of Montauk Highway, and 0.4 miles west of Heady Creek, near the geographic center of the tribal lands. Tidal-influence delay time (time interval between peak high-tide stage and corresponding peak high-water level) ranged from 1.75 hours at the well located near the southernmost extent of the tribal lands, to more than 4 hours at a well located north of Old Fort Pond, near the northwestern part of the tribal lands. Estimated hydraulic-conductivity values derived from the results of specific-capacity tests that were completed at nine observation wells during March 2015 were used to calculate average linear velocity. Average linear velocity along conceptualized flow-path segments of the upper glacial aquifer located beneath the tribal lands was estimated using an assumed effective porosity value, and hydraulic-conductivity and hydraulic-head values that were interpolated from measured values. Groundwater travel times were estimated by dividing the length of the flow-path segment by the average linear velocity along the flow-path segment. Total estimated groundwater travel time along a conceptualized flow path, beginning near Sunrise Highway and terminating at Shinnecock Bay, is approximately 45 years using a porosity value of 30 percent. A surficial-silty unit was identified from approximately 0 to 10 ft below land surface at multiple locations beneath the tribal lands. The lithology of the surficial unit was verified by interpreted gamma log results obtained from select wells, and auger-rig drill cuttings from an observation well located near the geographic center of the tribal lands. The altitude of the unit varies with topography and was delineated along a cross section line that trends north-south along the approximate centerline (spine) of the tribal lands. The altitude of the hydrogeologic contact between the upper glacial and the Magothy aquifers generally decreases from northwest to southeast, occurs at a depth ranging from about 150 to 200 ft beneath the tribal lands, and was identified at two locations north of the tribal lands, near Sunrise Highway and Sebonac Road. Results of electrical geophysical surveys indicate that the depth to the freshwater/saltwater interface decreases from north to south with decreasing water level altitude, and the Magothy and upper glacial aquifers contain saltwater at varying depths along the north-south trending section. Results of the surveys also indicate that the Magothy aquifer beneath the tribal lands contains brackish and salty water and is not considered a source of potable water supply. In general, depth to the interface increases with increasing geographic distance from the coastline. Low water table altitudes can result in increased saltwater encroachment into the surficial aquifer beneath the tribal lands. This upward movement and shallow depth of the freshwater/saltwater interface can jeopardize water quality in wells that supply water for domestic use.

New York↗