USGS Science⌕ Search

USGS · sir20165114

Hydrological conditions and evaluation of sustainable groundwater use in the Sierra Vista Subwatershed, Upper San Pedro Basin, southeastern Arizona

Abstract

This study assessed progress toward achieving sustainable groundwater use in the Sierra Vista Subwatershed of the Upper San Pedro Basin, Arizona, through evaluation of 14 indicators of sustainable use. Sustainable use of groundwater in the Sierra Vista Subwatershed requires, at a minimum, a stable rate of groundwater discharge to, and thus base flow in, the San Pedro River. Many of the 14 indicators are therefore related to long-term or short-term effects on base flow and provide us with a means to evaluate groundwater discharge to and base flow in the San Pedro River. The indicators were based primarily on 10 to 20 years of data monitoring in the subwatershed, ending in 2012, and included subwatershedwide indicators, riparian-system indicators, San Pedro River indicators, and springs indicators. Groundwater management actions including voluntary retirement of irrigation pumping in the subwatershed resulted in about a 5,100 acre-feet (acre-ft) reduction in net human use from 2002 to 2012. Subwatershed population increased more than 10,000 during the same period. Most of the reduction occurred during 2002–07 and included reductions in groundwater pumping and increases in managed recharge; net human use varied annually by a few hundred acre-ft during 2007–12. The groundwater budget for 2012 showed a deficit of about 5,000 acre-ft, although the total water-budget uncertainty was about 5,500 acre-ft. In the vicinity of the U.S. Army’s Fort Huachuca, regional-aquifer water levels were in steady decline beginning in at least the mid-1990s (in older wells since at least the early-1970s), as the cone of depression centered on the Sierra Vista and Fort Huachuca pumping centers continued to deepen. This was evident in the individual water levels on Fort Huachuca, as well as from the horizontal hydraulic gradients that extend from the pumping centers toward the San Pedro and Babocomari Rivers. Basin water levels in wells southeast of Sierra Vista, away from the river, were also experiencing declines, while some water levels closer to the river were rising. Near-stream vertical gradients along the San Pedro River showed no clear increasing or decreasing trends that would indicate a shift in the direction of subsurface flow between the riverbed and the alluvial aquifer, or a trend in the magnitude of groundwater/surface-water exchange. Annual streamflow permanence data showed no clear change in streamflow permanence trends in any of the river reaches, other than those related to precipitation trends. Similarly, the single-day, dry-season, wet-dry streamflow analysis of all subwatershed river reaches indicated no change in condition over the past 14 years, with the exception of the Hereford reach, which has seen a statistically significant increase in wetted length. Dry-season, alluvial-aquifer water levels in the Hereford reach also showed a statistically significant increase. These improvements are attributed to the end of irrigation pumping in the area. Although data indicate that the length of the Fairbank North wetted reach may be in decline, it is not yet statistically significant. Stable-isotope data indicated reduced groundwater discharge to the Babocomari River in the vicinity of the Babocomari River near Tombstone gaging station and to the San Pedro River near the San Pedro River at Palominas gaging station and near the Lewis Springs DCP stage recorder. The Babocomari River near Tombstone gaging station is downgradient of the major pumping centers. The change in isotopic signature at the Lewis Springs stage recorder could have been the result of alterations in groundwater/surface-water interactions there caused by beaver damming of the river. Base flow in the San Pedro River declined over the periods of record at the three San Pedro River gaging stations in the subwatershed (Palominas, Charleston, and Tombstone), as well as at the Babocomari River near Tombstone gaging station. Precipitation declined slightly from the 1990s to the 2000s, although there is no statistically significant trend in subwatershed precipitation from 1991 to 2012. The occurrence of large winter discharge events appeared to decline and that of large summer discharge events appeared to increase over this same period. Data for physical parameters, general chemistry, nutrient species, select trace elements, and suspended sediment were collected at San Pedro River at Charleston stream-gaging station. These data were summarized over time and analyzed in relation to discharge and season as a means to assess trends over the period of analysis. Federal and State of Arizona drinking-water and human-contact standards were all met and few exceedances occurred for the ecological thresholds investigated. Several constituents showed a significant trend over the period of analysis, but only concentration and flux data for total phosphate, orthophosphate, total nitrogen, suspended sediment, and sulfate were suitable to be used in a weighted regression analysis that statistically accounted for time, discharge, and season. Sulfate concentrations and flux showed a significant downward trend over the period of analysis, whereas total phosphorus and ortho-phosphate showed a relatively small magnitude upward trend relative to standards. Suspended sediment concentrations and flux both showed a significant downward trend in the 1980s, an effect attributed to reduction of cattle in the subwatershed at about this time, and (or) increased cottonwood ( Populus fremontii ) and willow ( Salix goodingii ) recruitment, and (or) the curtailment of sand and gravel mining adjacent to the San Pedro River with the designation of the San Pedro Riparian National Conservation Area in 1988. A spike in sediment flux in 2006 may be attributable to the more than 100 debris flows in the Huachuca Mountains during the summer monsoon of that year. Spring discharge along the San Pedro River generally increased at three sites proximate to the Sierra Vista treated effluent recharge facility and varied somewhat with climate at two other sites. Median annual discharge at the recharge facility peaked in 2006, and at Murray Springs and Horsethief Spring, downgradient of the recharge facility, in 2009. Sampling for trace organic compounds in flow from springs was carried out using both discrete sampling and passive sampling methods. Spring samples thus collected showed the presence of trace-organic compounds. Lewis Springs (background site) had the least number of detections, whereas Murray Springs, located directly downgradient of the City of Sierra Vista’s treated effluent recharge facility, had the greatest number of detections of all the springs. Discrete samples from the recharge facility had more than twice the detections found in discrete samples from Murray Spring and at much higher concentrations. Few similar trace-organic compounds were detected at both the springs and the treated effluent recharge facility, and the number of detections did not increase during the collection period. Limitations of the study prevented the determination of trace-organic concentration in passive samplers and also prevented linking trace organic compounds detected at the treated effluent recharge facility with compounds detected from the springs. In particular, trace organic compounds could also derive from other sources such as septic systems. Looking at the subwatershed as a whole, base flow was in decline along the entire river reach, but determination of the specific cause of the decline was beyond the scope of this report. Conditions in the area from the municipal pumping center of Sierra Vista and Fort Huachuca northeast to the river (from about the Charleston to Tombstone gaging stations) were more commonly in decline than in regions further south. Both long-term indicators, such as regional aquifer groundwater levels and horizontal gradients, and the isotope analysis indicated that groundwater discharge to the river and thus base flow may continue to decline in that area. South of Charleston, indicators were more mixed. Some indicators in the Hereford reach suggest groundwater discharge to the San Pedro River may be increasing there, whereas some indicators in the Palominas reach suggest groundwater discharge to the river there may be declining.

Explore related subjects

90° N90° S · 180° W ← longitude → 180° E
Source-reported bounding extent: 31.335° to 31.8° latitude; -110.59° to -109.86328125° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Bruce Gungle, James B. Callegary, Nicholas V. Paretti, Jeffrey R. Kennedy, Christopher J. Eastoe, Dale S. Turner, Jesse E. Dickinson, Lainie R. Levick, Zachary P. Sugg. 2016-08-18. Hydrological conditions and evaluation of sustainable groundwater use in the Sierra Vista Subwatershed, Upper San Pedro Basin, southeastern Arizona. https://doi.org/10.3133/sir20165114

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

Baseline water-quality inventory and summary for the Paterson Great Falls National Historical Park, New Jersey

The Paterson Great Falls National Historical Park (PAGR), established in 2011 in the City of Paterson, New Jersey, contains many cultural, historical, and natural resources. Because of the PAGR’s young age and limited resources, a long record of natural resource data is not available for use to guide future park management and planning decisions. A detailed natural resource inventory was identified as an immediate need in the “Foundation Document—Paterson Great Falls National Historical Park,” given the history of extensive industrial activity throughout the PAGR and the Passaic River Basin. The U.S. Geological Survey, in cooperation with the National Park Service and the Urban Waters Federal Partnership for the Lower Passaic River, inventoried water-quality monitoring efforts and compiled historical water-quality data from secondary data sources to define baseline conditions in the PAGR and the surrounding study area. This report documents the methods used to harmonize water-quality data aggregated from multiple sources for the PAGR study area. Procedures implemented during data harmonization include compiling data from multiple sources, assessing the dataset for completeness, addressing metadata concerns, and creating a consistent and quality-assured dataset for future analysis. The harmonized dataset contains 14 aggregated surface-water sites, 1,103 individual sampling days, and 32,000 water-quality records throughout the period of record defined for this study (1952–2021). Summary statistics and comparison to water-quality criteria are provided for focus parameters that can support PAGR management, including field properties, suspended material, total dissolved solids, specific conductance, chloride, sulfate, Escherichia coli, nutrients, contaminants of concern related to fish consumption advisories, and litter and floatable debris. A spatial data gap was identified for the recreational area near the base of the Great Falls and immediately downstream from a known combined sewer overflow discharge location. Analytic data gaps were identified for continuous streamflow and water-quality parameters, floatable debris characterization, and harmful algal blooms identification. Continuous collection of contaminants of concern was identified as a temporal data gap, as these parameters were not consistently collected within the study area. The water-quality data synthesis and summary presented in this report can assist PAGR staff in understanding the complex baseline water resources within the study area and can serve as a baseline to design future water-quality monitoring efforts and evaluate changes in water quality over time. The harmonized dataset, as well as future New Jersey Department of Environmental Protection and New Jersey Harbor Discharge Group monitoring data and assessments, can also be used to examine the status and trends in surface-water quality, support compliance with the Clean Water Act (33 U.S.C. 1251 et seq.), and further investigate changes to park water resources that can affect ecological and human health.

New Jersey↗

Slip tendency and dilation tendency of Quaternary faults in the Great Basin region, Arizona, California, Idaho, Nevada, Oregon, Utah, and Wyoming

This report presents the results of slip tendency and dilation tendency analysis of Quaternary faults in the Great Basin region, in the states of Arizona, California, Idaho, Nevada, Oregon, Utah, and Wyoming. This effort was undertaken as part of two studies: the U.S. Department of Energy Geothermal Technologies Office-funded Innovative Geothermal Exploration through Novel Investigations of Undiscovered Systems (INGENIOUS) project, which is a multidisciplinary, multi-institution effort to develop new methodologies and best practices to accelerate the discovery of new, commercially viable hydrothermal resources; and the U.S. Geological Survey Geothermal Resources Investigations Project. The goal of this work is to identify faults and fault segments that are appropriately oriented to be stress-loaded for slip or to dilate under the ambient stress conditions. Both conditions may make such faults more likely to host permeability and fluid flow relative to other faults in the region. Slip tendency and dilation tendency were calculated for all fault segments in the study area from the U.S. Geological Survey Quaternary Fault and Fold Database. Results demonstrate that faults with high slip tendency and dilation tendency are common in the Great Basin region. The areas containing many high slip tendency and dilation tendency fault segments are associated with high strain rates. These areas also host a relatively high proportion of the known hydrothermal systems. This indicates that faults with high slip tendency and dilation tendency in areas with high strain-rate are potentially most prospective for hosting hydrothermal systems. The database published in coordination with this report contains the slip tendency and dilation tendency results data, including the fault geometry and stress data used to calculate these quantities, for both the Great Basin study area and the smaller INGENIOUS study area.

Arizona, California, Idaho, Nevada, Oregon, Utah, ↗

A spatial analysis of the groundwater emergence flood hazard in Long Island, New York, and coastal areas surrounding Long Island Sound in New York, Connecticut, and Rhode Island

Long Island, New York, and nearby coastal areas surrounding Long Island Sound are densely populated and, like other coastal areas, are susceptible to flooding from several potential sources, including stormwater from precipitation events, tidal flooding and storm surge, and groundwater emergence flooding. Groundwater emergence fooding refers to the intersection of a rising water table with land surface or critical infrastructure. Many studies of flood drivers either neglect or only briefly discuss how shallow groundwater conditions may contribute to or exacerbate flood conditions. As a part of a comprehensive study of compound flood hazards in the coastal areas of Long Island and Long Island Sound, a spatial analysis was completed by the U.S. Geological Survey, in cooperation with the Environmental Protection Agency’s Long Island Sound Partnership, formerly called the Long Island Sound Study, using available regional datasets to characterize the potential hazard for groundwater emergence flooding. The approximately 3,100 square mile study area was subdivided into 11,407 grid cells that were 900 meters by 900 meters (approximately 3,000 feet by 3,000 feet), for the purpose of integrating the spatial datasets to calculate and map the groundwater emergence flood hazard. The depth to the water table, hydrologic soil groups, and National Land Cover Database were harmonized to the common grid. A groundwater emergence flood hazard score was calculated for each grid cell for current average conditions following a set of rules accounting for the depth to the water table and the percentage of area within each cell with slow infiltrating soils. A scenario with a higher sea-level position was also calculated for the Long Island part of the study area. The calculated groundwater emergence flood hazard score was reviewed in conjunction with the National Land Cover Database to identify developed areas and associated infrastructure that may be at risk of potential groundwater emergence flooding. Study results indicate that the groundwater emergence flood hazard is highest in coastal areas and near surface water where the water table is close to ground surface. Inland areas away from surface waterbodies are not likely to be exposed to groundwater emergence flooding. For Long Island, in the scenario with higher sea-level position, a greater groundwater emergence flood hazard is calculated for some locations closer to the coast and where land is submerged. Away from the coast and surface-water drainage, the groundwater emergence flood hazard is similar between the current average sea-level condition and the higher sea-level position scenario.

Connecticut, New York, Rhode Island↗