USGS Science⌕ Search

USGS · sir20215039

Occurrence, fate, and transport of aerially applied herbicides to control invasive buffelgrass within Saguaro National Park Rincon Mountain District, Arizona, 2015–18

Abstract

The spread of the invasive and fire-adapted buffelgrass ( Cenchrus ciliaris L.) threatens desert ecosystems by competing for resources, increasing fuel loads, and creating wildfire connectivity. The Rincon Mountain District of Saguaro National Park addressed this natural resource threat with the use of glyphosate-based herbicides (GBHs). In 2010, the Rincon Mountain District initiated an aerial restoration plan to control dense buffelgrass patches in remote areas and implemented a trial project to evaluate the effects of aerial restoration techniques that included the helicopter application of GBHs. In 2014, more than 250 acres of buffelgrass in the Rincon Mountain District were treated with the aerial application of GBHs. This widespread aerial application of GBHs continued through 2018, but the potential transport and effects to aquatic ecosystems were unknown. In 2015–18, the U.S. Geological Survey, in cooperation with the National Park Service, studied the occurrence, distribution, fate, and transport of glyphosate in surface water and sediments derived from areas that were treated during past and current aerial herbicide applications. Three watersheds, treated with different regimens of GBHs, were sampled for glyphosate and the primary metabolite of glyphosate, aminomethylphosphonic acid (AMPA), during various hydrologic flow conditions. Water and aquatic sediment were collected from three watersheds, each in a different stage of application during the U.S. Geological Survey study. The unnamed watershed above the Loma Verde Trailhead referred to by the National Park Service as “Loma Verde canyon” had received no aerial treatment since 2014, whereas the Box Canyon watershed was aerially treated every year beginning in 2014. The Madrona Canyon watershed was first sprayed in 2016 and aerial application continued once a year though the entirety of the study. In addition, terrestrial soil samples were sampled from areas sprayed to understand dissipation rates and herbicide transport via sediments washing away during rainfall runoff. The concentrations present in water and sediment samples were compared to ecological benchmarks and characterized within the context of the environmental conditions of the park setting. Of the 48 water samples collected and analyzed for glyphosate and AMPA, 10.4 percent and 14.6 percent were detected above the laboratory minimum detection limit, respectively. Mean water concentrations, calculated using specific statistical methods for non-detects, were equal to the laboratory minimum detection limit of 0.02 microgram per liter for samples collected in all the watersheds. In aquatic sediments, glyphosate and AMPA were detected in 10.7 and 25.0 percent of the samples, whereas 89.5 and 100 percent of the terrestrial soil samples had detections for glyphosate and AMPA, respectively. Mean aquatic sediment concentrations were 1.13 and 4.42 micrograms per kilogram (μg/kg) for glyphosate and AMPA, respectively. Mean terrestrial soil concentrations were orders of magnitude greater than water and aquatic sediment with concentrations of 678 μg/kg for AMPA and 1,240 μg/kg for glyphosate. Hours after glyphosate-based herbicide was applied, the concentrations of glyphosate and AMPA were present in terrestrial soil samples near or above the laboratory maximum detection limit of 5,000 μg/kg. The Box Canyon watershed was the most intensively treated watershed in terms of total land area treated, total amount of GBH applied, and number of years treated. The frequent and large volume of treatment resulted in the highest number of detections of glyphosate and AMPA in water (3 and 7 detections, respectively) and in aquatic sediment (2 and 6 detections, respectively) samples. In comparison, the other two watersheds had two or fewer detections for glyphosate and AMPA in water and aquatic sediment. Glyphosate detected in pools was associated with increased rainfall closer in time to the last herbicide treatment. Glyphosate and AMPA concentration ratios above one, along with stable-isotope and tritium results, indicated that runoff processes were the primary transport mechanism for the two compounds when found in streams and pools rather than subsurface recharge or deeper flow paths. One pool in a small tributary of Box Canyon consistently had detections of glyphosate and AMPA in aquatic sediments, but these frequent concentrations were likely related to the intensive application upstream, near the steep terrain above the head of the channel that supplies the downstream pool. Intense flows during summer rainfall events move treated sediments into this channel where vegetation and the incised bedrock banks of the pool retained those sediments and ultimately led to frequent detections of both compounds. Isotope results in most of the pools and tinajas indicated that the water source had residence time representative of recently recharged waters, on the order of years. No water concentrations exceeded published criteria for human health or aquatic life. Median and maximum glyphosate and AMPA water concentrations were lower than those reported in other national assessments, but maximum concentrations observed in individual runoff samples were higher than median concentrations measured in the national assessments. A similar finding was observed with aquatic sediment concentrations measured in the Rincon Mountain District. Results from the study were compared and assessed in the context of other studies examining GBHs and their effects on amphibians, fish, and macroinvertebrates. This comparison was used to generalize the potential risk to aquatic species similar to those species in the Rincon Mountain District. Concentrations of published effect levels were several orders of magnitude greater than the highest concentration detected in water at the Rincon Mountain District. Most published studies evaluate acute and chronic toxicity for glyphosate and GBHs, and these criteria may not be representative of environmental conditions in the Rincon Mountain District. The classic lethal dose studies conducted in a controlled laboratory setting may not be suitable for comparison to the longer, variable, low-dose exposure conditions in the pools and tinajas in the Rincon Mountain District. However, this study determined that the fate of GBHs transported from treated areas to potential aquatic habitat was highly variable in occurrence, timing, and concentrations. This variability in glyphosate concentrations was too high, and the potential exposure was determined to be far too complex to directly compare with the results from controlled studies. This study provides the first information collected on GBHs used to control invasive buffelgrass in a remote, mountainous, and semiarid setting. The information about the transport and fate of herbicide application near aquatic habitat will help to inform managers about the broader ecosystem implications and provide useful information to other agencies implementing buffelgrass remediation strategies near aquatic habitat.

Explore related subjects

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

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Nicholas V. Paretti, Kimberly R. Beisner, Bruce Gungle, Michael T. Meyer, Bethany K. Kunz, Edyth Hermosillo, Jay R. Cederberg, Justine P. Mayo. 2021-09-01. Occurrence, fate, and transport of aerially applied herbicides to control invasive buffelgrass within Saguaro National Park Rincon Mountain District, Arizona, 2015–18. https://doi.org/10.3133/sir20215039

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

KEEP EXPLORING

Related USGS reports

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↗

Streamflow assessment of Birch Creek, southern boundary of the Blackfeet Indian Reservation, Pondera County, northwestern Montana—Synoptic and continuous streamflow measurements through 2022

The U.S. Geological Survey, in cooperation with the Blackfeet Water Department, performed this study to assist water-resource managers in understanding and quantifying irrigation diversions, instream flows, and the spatial and temporal streamflow dynamics of Birch Creek in northwestern Montana. Birch Creek forms the southernmost boundary of the Blackfeet Indian Reservation and is highly valued culturally as an aquatic life and wildlife resource and as a principal source of water for irrigation. The study area included Birch Creek from Swift Dam to 53.3 miles downstream just above its confluence with the Two Medicine River. The study period was calendar year 2022, but selected historical information was also included in the study. Birch Creek streamflow originates through releases from Swift Reservoir, with the largest releases typically occurring during the irrigation season (May–October) to primarily meet the needs of downstream irrigators. This study included the compilation and collection of (1) continuous streamflow measured at streamgages installed for this study, (2) instantaneous streamflow measurements made during four synoptic events, and (3) reported daily values of Swift Reservoir outflow and irrigation diversions. Continuous streamflow data were available from two U.S. Geological Survey streamgages: Birch Creek below Heart Butte Road, near Heart Butte, Montana (upper streamgage; station 06094900), at river mile 12.9 and Birch Creek at Bullhead Bridge near Valier, Mont. (lower streamgage; station 06098120), at river mile 53.3. The upper streamgage (station 06094900) was primarily affected by the regulated releases from Swift Dam with daily mean streamflow typically slightly higher than the reported Swift Reservoir outflows. Daily mean streamflow ranged from 9.87 cubic feet per second (ft 3 /s) (November 18, 2022) to 640 ft 3 /s (June 22, 2022) during September 16, 2021–December 31, 2022. The mean annual streamflow was 122 ft 3 /s in water year 2022 (a water year is the 12-month period from October 1 to September 30 and is designated by the year in which it ends), about 13 percent greater than the mean annual reported Swift Reservoir outflows. The lower streamgage (station 06098120) was downstream from all major tributaries and diversions within the Birch Creek drainage basin. Daily mean streamflow ranged from 6.02 ft 3 /s (September 8, 2022) to 95.0 ft 3 /s (June 5, 2021) during August 26, 2020–December 31, 2022. Mean annual streamflow was 27.8 ft 3 /s for water year 2022. Flows at the lower streamgage (station 06098120) were typically lower than at the upper streamgage (station 06094900) from late April to early October, whereas flows from early October to late spring were similar to or greater than those at the upper streamgage. The upper streamgage had a total streamflow volume of 87,960 acre-feet (acre-ft) for water year 2022, whereas the lower streamgage had a total streamflow of 20,357 acre-ft. A net volume of about 67,603 acre-ft was removed from Birch Creek through diversions, losses to groundwater, and evapotranspiration along that reach during water year 2022. Most of the changes in main-stem streamflow could be accounted for by the tributary inflows or diversion outflows. The most dynamic streamflow gains and losses occurred along the upper few miles of Birch Creek. Streamflow was consistently lost to the groundwater system immediately below Swift Dam and then returned to levels that slightly exceeded the outflow at the dam near river mile 11.0. The overall net gains for the entire 53.3-mile study reach ranged from 15.5 to 42.6 ft 3 /s among the synoptic events. On several occasions, streamflow fell below minimum instream flow conditions specified in the Blackfeet Compact (water rights compact among the Blackfeet Tribe, State of Montana, and United States), which are subject to a 15-year deferral period (until June 12, 2033). Notably, a section of Birch Creek immediately downstream from Swift Dam was completely dry during the October 18, 2022, synoptic event.

Montana↗