SEARCH · USGS Science

Results for “Environmental Data Science”

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

SourceResearch connectionFrom yearThrough year

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

At least 1,099 records · Page 61Linked to original sources

USGS2006

Water-Quality Data for Selected National Park Units within the Southern Colorado Plateau Network, Arizona, Utah, Colorado, and New Mexico, Water Years 2005 and 2006

The National Park Service initiated a Level 1 Water-Quality Inventory program to provide water-quality data to park managers so informed natural resource management decisions could be made. Level 1 water-quality data were collected by the U.S. Geological Survey Arizona Water Science Center at 57 sites in 13 National Park units located in the Southern Colorado Plateau Inventory and Monitoring network in water years 2005 and 2006. These data describe the current water-quality at selected sites within the park units and provide information for monitoring future trends. Water samples were collected three times at each type of site including wells, springs, seeps, tinajas, rivers, a lake, and an irrigation ditch. Field measurements were taken at each site and they included pH, specific conductance, temperature, barometric pressure, dissolved oxygen, alkalinity, turbidity, and discharge rates where applicable. Water samples collected were sent to the U.S. Geological Survey National Water Quality Laboratory and analyzed for major ions, trace elements, and nutrients. The National Water Quality Laboratory also analyzed selected samples for mercury and petroleum hydrocarbons. Additional samples at selected sites were collected and analyzed for cyanide, radiochemistry, and suspended sediment by U.S. Geological Survey contract labs. Fecal-indicator bacteria (Escherichia coli) were sampled for at selected sites as another indicator of water quality. Quality control for this study was achieved through proper training of field personnel, use of standard U.S. Geological Survey field and laboratory protocols, collection of sample blanks and replicates, and a thorough review of the water-quality analyses. Measured field pH ranged from 6.0 to 8.8, within normal range for springs and rivers, at most sites. Concentrations of dissolved solids ranged from 48 to 8,680 mg/L and the majority of samples had concentrations of dissolved solids below 900 mg/L. Trace-element concentrations at most sites were at or near the laboratory reporting levels. The highest overall trace-element concentrations were found at U.S. Highway 160 Spring near Park Entrance to Mesa Verde National Park. Concentrations of uranium in samples at all sites ranged from below the detection limit to 55.7 ?g/L. Water samples from selected sites were analyzed for total petroleum hydrocarbons and concentrations of total petroleum hydrocarbons were at or above the laboratory detection limit in samples at six National Park units. Ten sites were sampled for Escherichia coli and positive counts were found at 9 out of the ten sites, the highest colony counts were found at Chinle Creek at Chinle, AZ in Canyon de Chelly National Monument. Measured concentrations of dissolved ammonia, nitrite, and nitrate were at or near laboratory reporting levels at most sites; nitrate concentrations ranged from below the reporting limit (0.047 mg/L) to 9.77 mg/L. Samples that were analyzed for mercury had concentrations below or at the laboratory reporting level. Concentrations of cyanide were less than the laboratory reporting level for all samples except two, Spruce Tree House Spring in Mesa Verde National Park and Pine Tree Canyon Tinaja in Canyon de Chelly National Monument, which had average concentrations of .042 and .011 ?g/L respectively. Gross alpha/beta radioactivity counts were below the U.S. Environmental Protection Agency maximum contaminant level except for samples from Casa Chiquita Well Middle at Chaco Culture National Historical Park which averaged 35 pCi/L. Suspended-sediment concentrations were variable and ranged from 10 to 150,000 mg/L.

Open-File Report↗

USGS1998

Occurrence of pesticides in shallow groundwater of the United States: initial results from the National Water-Quality Assessment program

The first phase of intensive data collection for the National Water-Quality Assessment (NAWQA) was completed during 1993−1995 in 20 major hydrologic basins of the United States. Groundwater land-use studies, designed to sample recently recharged groundwater (generally within 10 years) beneath specific land-use and hydrogeologic settings, are a major component of the groundwater quality as sessment for NAWQA. Pesticide results from the 41 land-use studies conducted during 1993−1995 indicate that pesticides were commonly detected in shallow groundwater, having been found at 54.4% of the 1034 sites sampled in agricultural and urban settings across the United States. Pesticide concentrations were generally low, with over 95% of the detections at concentrations less than 1 μg/L. Of the 46 pesticide compounds examined, 39 were detected. The compounds detected most frequently were atrazine (38.2%), deethylatrazine (34.2%), simazine (18.0%), metolachlor (14.6%), and prometon (13.9%). Statistically significant relations were observed between frequencies of detection and the use, mobility, and persistence of these compounds. Pesticides were commonly detected in both agricultural (56.4%; 813 sites) and urban (46.6%; 221 sites) settings. Frequent detections of pesticides in urban areas indicate that, as is the case with agricultural pesticide use in agricultural areas, urban and suburban pesticide use significantly contribute to pesticide occurrence in shallow groundwater. Although pesticides were detected in groundwater sampled in urban areas and all nine of the agricultural land-use categories examined, significant variations in occurrence were observed among these categories. Maximum contaminant levels (MCLs) established by the U.S. Environmental Protection Agency for drinking water were exceeded for only one pesticide (atrazine, 3 μg/L) at a single location. However, MCLs have been established for only 25 of the 46 pesticide compounds examined, do not cover pesticide degradates, and, at present, do not take into account additive or synergistic effects of combinations of pesticide compounds or potential effects on nearby aquatic ecosystems.

Environmental Science & Technology↗

USGS2023

Ecological effects of pinyon-juniper removal in the Western United States—A synthesis of scientific research, January 2014–March 2021

Executive Summary Increasing density of pinyon ( Pinus spp. ) and juniper ( Juniperus spp. ) woodlands (hereinafter “pinyon-juniper”), as well as expansion of these woodlands into adjacent shrublands and grasslands, has altered ecosystem function and wildlife habitat across large areas of the interior western United States. Although there are many natural and human-caused drivers of woodland infilling and expansion, restoration of sagebrush ( Artemisia spp. ) habitat through removal of pinyon-juniper is considered an urgent management objective in many locations, particularly in support of sagebrush-dependent wildlife species of conservation concern. In December 2020, the Bureau of Land Management (BLM) established the Pinyon-Juniper Management Categorical Exclusion (PJCX) to expedite the regulatory process for pinyon-juniper removal projects on public lands, largely intended to benefit mule deer ( Odocoileus hemionus ) and greater sage-grouse ( Centrocercus urophasianus ) habitats. During final preparation of this report, the BLM discontinued use of the PJCX (as of November, 2022), but the pinyon-juniper tree removal techniques assessed in this report are commonly used and understanding their effects remains relevant to land use planning. To address areas of uncertainty relative to potential ecological effects of the PJCX, we conducted a review of the peer-reviewed science literature to better understand the likely responses of vegetation, environmental (for example, soils), and wildlife variables to specific tree removal techniques permitted by the PJCX. In brief, the PJCX permitted removal of trees by either manual cutting, mechanical cutting, or mastication; allowed certain methods to redistribute or remove resulting tree biomass after treatment; and prohibited broadcast burning, roadbuilding, removal of old-growth, and seeding of non-native species. Specifically, we conducted our review to address the following questions: How will PJCX removal techniques affect plant communities, soils, and abiotic resources? How do these pinyon-juniper removal techniques affect wildlife communities, including both woodland- and sagebrush-dependent species? What are the potential ecological implications of different pinyon-juniper removal treatment types and implementation strategies (for example, treatment sizes) over time? What are the most important gaps in our scientific understanding of how treatments might affect targeted ecosystems over space and time (for example, potential effects of climate change)? To answer these questions, we considered studies related to pinyon-juniper ecosystems, focusing on research that occurred over a large portion of the interior western United States that is the primary focus of the PJCX. We also focused on papers published from 2014 onward, to avoid excessive overlap with other recent reviews on pinyon-juniper management effects. Using strict criteria, including only considering research that tested responses for statistical significance, we identified 48 papers that primarily examined treatment effects on vegetation and other environmental variables (1,709 responses), and 11 papers that addressed effects on wildlife (132 responses). Responses to the PJCX-permitted treatments were summarized as either positive (that is, a significant increase), negative (that is, a significant decrease), or non-significant (that is, no significant difference). Responses were assigned to categories (for example, Native Annual Grass/Forb Abundance) and hierarchical treatment levels. We found that there were large proportions of non-significant responses among all categories combined, with roughly half or more of all responses non-significant (48 percent for wildlife, 60 percent for vegetation-environmental), comparable to other recent systematic reviews of pinyon-juniper treatment effects. However, we also found that when there were significant responses, some important trends potentially emerged. Important undesirable outcomes included far more positive than negative responses of exotic grass and forb abundance among nearly all treatment types. Cutting treatments were also more likely to decrease biocrust cover and microbial activity. Potentially beneficial outcomes included mostly positive responses among sagebrush obligate species, including more positive than negative responses for mule deer and sage-grouse. Some treatment types (for example, mastication) also resulted in more positive than negative responses for native grasses and forbs (although, non-significant responses were the majority). We also highlighted many limitations of this review, including how responses often come from few studies, and how some response-treatment category combinations lack adequate response data. Moreover, the existing research is often insufficient to address many key questions about treatment effects, largely owing to short time-scales and limited spatial extents of observations, which do not match the size of treatments being implemented by land managers, nor capture long-term, post-treatment ecological dynamics. We also identify a lack of research that addresses key interactions that could undermine restoration objectives, including potential effects of climate change and grazing on post-treatment environments. Thus, we emphasize the importance of integrating these factors into future pinyon-juniper treatment research, and we stress the need for use of monitoring programs and research studies that partake in data collection and analysis over long durations and broad spatial scales.

Open-File Report↗

USGS2021

Improved wetland soil organic carbon stocks of the conterminous U.S. through data harmonization

Wetland soil stocks are important global repositories of carbon (C) but are difficult to quantify and model due to varying sampling protocols, and geomorphic/spatio-temporal discontinuity. Merging scales of soil-survey spatial extents with wetland-specific point-based data offers an explicit, empirical and updatable improvement for regional and continental scale soil C stock assessments. Agency-collected (U.S. Department of Agriculture, U.S. Environmental Protection Agency) and community-contributed soil datasets were compared for representativeness and bias, with the goal of producing a harmonized national map of wetland soil C stocks with error quantification for wetland areas of the conterminous United States (CONUS) identified by the USGS National Landcover Change Dataset (NLCD). This allowed application of an empirical predictive model of SOC density to be applied across the entire CONUS using relational %OC distribution alone. A broken-stick quantile-regression model identified %OC with its relatively high analytical confidence as a key predictor of SOC density in soil segments; soils less than 6%OC (hereafter, mineral wetland soils, 85% of the dataset) had a strong linear relationship of %OC to SOC density (RMSE = 0.0059, ~4% mean RMSE) and soils greater than 6%OC (organic wetland soils, 15% of the dataset) had virtually no predictive relationship of %OC to SOC density (RMSE = 0.0348 g C cm-3, ~56% mean RMSE). Disaggregation by vegetation type (woody v. emergent herbaceous), or region did not alter the breakpoint significantly (6% OC) nor improve model accuracies for inland and tidal wetlands. Similarly, SOC stocks in tidal wetlands were related to %OC, but without a mappable product for disaggregation to improve accuracy by soil class, region or depth. Our layered, harmonized CONUS wetland soil maps have now revised wetland SOC stock estimates downward by 24% (9.5 vs. 12.5Pg C) with the overestimation being entirely an issue of inland, organic wetland soils, (35% lower than SSURGO-derived SOC stocks). Further, SSURGO underestimated soil carbon stocks at depth, as modeled wetland SOC stocks for organic-rich soils showed significant preservation downcore in the NWCA dataset (<3% loss between 0-30 cm and 30-100 cm depths) in contrast to mineral-rich soils (37% downcore stock loss). Future CONUS wetland soil C assessments will benefit from focused attention on improved organic wetland soil measurements, land history, and spatial representativeness.

Frontiers in Soil Science↗

USGS1985

Management of ground water and evolving hydrogeologic studies in New Jersey : a heavily urbanized and industrialized state in the northeastern United States

New Jersey is the most densely populated and one of the most industrialized states in the United States. An abundance of freshwater and proximity to major northeastern metropolitan centers has facilitated this development. Pumpage of freshwater from all aquifers in the State in 1980 was 730 million gallons per day (2.76 million cubic meters per day). Management and efficient development of the ground-water resources of the State are the responsibility of the New Jersey Department of Environmental Protection. Laws have been enacted and updated by the State legislature to manage water allocation and to control the disposal of hazardous wastes. Present resource management is guided by the New Jersey Water-Supply Master Plan of 1981. Funding for management activities is partially derived from the sale of state-approved bonds. Effective planning and regional management require accurate and up-to-date hydrologic information and analyses. The U.S. Geological Survey, in cooperation with the New Jersey Geological Survey, is conducting three intensive ground-water studies involving the collection and interpretation of hydrologic data to meet the urgent water-management needs of New Jersey. These studies are part of a long-term cooperative program and are funded through the Water-Supply Bond Act of 1981. They began in 1983 and are scheduled to be completed in 1988. The project areas are situated in the New Jersey part of the Atlantic Coastal Plain in and near Atlantic City, Camden, and South River. They range in size from 400 to 1,200 mil (1,040 to 3,120 km2). The studies are designed to define the geology, hydrology, and geochemistry of the local ground-water systems. The results of these studies will enable the State to address more effectively major problems in these areas such as declining water levels, overpumping, saltwater intrusion, and ground-water contamination resulting from the improper disposal of hazardous wastes. Specific objectives of these studies by the U.S. Geological Survey are to (1) develop an accurate and up-to-date hydrogeologic data base, (2) design and implement a data-collection program and establish a computerized information management system, (3) refine the conceptualization of the ground-water flow system, and (4) define the geochemistry of the aquifer system by conducting a water-quality appraisal. The objectives are accomplished by standard hydrogeologic methods. Information concerning hydrogeologic framework, ground-water levels, water use, hydraulic characteristics, and water quality in the study areas is compiled from all available sources. Additional data needed are collected through well inventories, surface geophysical surveys, water-quality samplings, water-level measurements, and a well-drilling program. Interpretation of the flow system is based on the use of standard analytical techniques and digital flow modeling. Calibrated flow models will provide ground-water managers with a mechanism to develop and test regional water-supply strategies. Definition of the geochemistry of the aquifer system is accomplished through a variety of methods which depend on the problems and available data in the particular study area. The approach includes statistical analysis of water-quality data, reaction-path modeling, and determination of the movement of chemical constituents using analytical and numerical modeling techniques. A combined staff of 25 to 30 professionals and technicians from the New Jersey District office of the U.S. Geological Survey is committed to the three studies. The staff has specialists in geohydrology, numerical modeling, geochemistry, geophysics, and computer science. The findings of these studies will be published in data reports, interpretive reports, instructional manuals and journal articles.

New Jersey↗

USGS2019

Accuracy of methods for reporting inorganic element concentrations and radioactivity in oil and gas wastewaters from the Appalachian Basin, U.S. based on an inter-laboratory comparison.

Accurate and precise analyses of oil and gas (O&G) wastewaters and solids ( e.g. , sediments and sludge) are important for the regulatory monitoring of O&G development and tracing potential O&G contamination in the environment. In this study, 15 laboratories participated in an inter-laboratory comparison on the chemical characterization of three O&G wastewaters from the Appalachian Basin and four solids impacted by O&G development, with the goal of evaluating the quality of data and the accuracy of measurements for various analytes of concern. Using a variety of different methods, analytes in the wastewaters with high concentrations ( i.e. , >5 mg L −1 ) were easily detectable with relatively high accuracy, often within ±10% of the most probable value (MPV). In contrast, often less than 7 of the 15 labs were able to report detectable trace metal(loid) concentrations ( i.e. , Cr, Ni, Cu, Zn, As, and Pb) with accuracies of approximately ±40%. Despite most labs using inductively coupled plasma mass spectrometry (ICP-MS) with low instrument detection capabilities for trace metal analyses, large dilution factors during sample preparation and low trace metal concentrations in the wastewaters limited the number of quantifiable determinations and likely influenced analytical accuracy. In contrast, all the labs measuring Ra in the wastewaters were able to report detectable concentrations using a variety of methods including gamma spectroscopy and wet chemical approaches following Environmental Protection Agency (EPA) standard methods. However, the reported radium activities were often greater than ±30% different to the MPV possibly due to calibration inconsistencies among labs, radon leakage, or failing to correct for self-attenuation. Reported radium activities in solid materials had less variability (±20% from MPV) but accuracy could likely be improved by using certified radium standards and accounting for self-attenuation that results from matrix interferences or a density difference between the calibration standard and the unknown sample. This inter-laboratory comparison illustrates that numerous methods can be used to measure major cation, minor cation, and anion concentrations in O&G wastewaters with relatively high accuracy while trace metal(loid) and radioactivity analyses in liquids may often be over ±20% different from the MPV.

Appalachian Basin↗

USGS2001

Constraining controls on carbonate sequences with high-resolution chronostratigraphy: Upper Miocene, Cabo de Gata region, SE Spain

A high-resolution chronostratigraphy has been developed for Miocene shallow-water carbonate strata in the Cabo de Gata region of SE Spain for evaluation of local, regional and global factors that controlled platform architecture prior to and during the Messinian salinity crisis. Paleomagnetic data were collected from strata at three localities. Mean natural remanent magnetization (NRM) ranges between 1.53 ?? 10-8 and 5.2 ?? 10-3 Am2/kg. Incremental thermal and alternating field demagnetization isolated the characteristic remanent magnetization (ChRM). Rock magnetic studies show that the dominant magnetic mineral is magnetite, but mixtures of magnetite and hematite occur. A composite chronostratigraphy was derived from five stratigraphic sections. Regional stratigraphic data, biostratigraphic data, and an 40Ar/39Ar date of 8.5 ?? 0.1 Ma, for an interbedded volcanic flow, place the strata in geomagnetic polarity Chrons C4r to C3r. Sequence-stratigraphic and diagenetic evidence indicate a major unconformity at the base of depositional sequence (DS)3 that contains a prograding reef complex, suggesting that approximately 250 000 yr of record (Subchrons C3Br.2r to 3Br.1r) are missing near the Messinian-Tortonian boundary. Correlation to the GPTS shows that the studied strata represent five third- to fourth-order DSs. Basal units are temperate to subtropical ramps (DS1A, DS1B, DS2); these are overlain by subtropical to tropical reefal platforms (DS3), which are capped by subtropical to tropical cyclic carbonates (Terminal Carbonate Complex, TCC). Correlation of the Cabo de Gata record to the Melilla area of Morocco, and the Sorbas basin of Spain indicate that early - Late Tortonian ramp strata from these areas are partially time-equivalent. Similar strata are extensively developed in the Western Mediterranean and likely were influenced by a cool climate or influx of nutrients during an overall rise in global sea-level. After ramp deposition, a sequence boundary (SB3) in Cabo de Gata correlates with a sequence boundary in Morocco and a published third-order eustatic fall suggesting at least a partial eustatic control for the sequence boundary. Coral reefs began to develop earlier in Cabo de Gata than at Melilla or Sorbas, arguing for local factors affecting this major environmental transition. Later Messinian reefs (DS3) from all areas are time-equivalent, suggesting a regional or global control on their formation. Some Halimeda-rich horizons in the Western Mediterranean are not time-equivalent event strata as hypothesized by others. Correlation of the relative sea-level curve for the fringing reef complex (DS3) with a published eustatic curve suggests at least a partial third-order global eustatic control for the highstand part of the sequence. Downstepping DS3 reefs and initial subaerial exposure of earlier DS3 reef strata approximately correlate with initiation of a series of subaerial unconformities in the South Pacific. The longer-term relative fall in sea-level during DS3 downstepping reef progradation does not correlate with a published third-order eustatic fall. Eustatic sea-level fluctuations may have been associated with initiation of the Mediterranean Messinian salinity crisis, but the longer-term fall may have been linked to tectonic uplift in the Mediterranean region. Widespread distribution of 'TCC-style' cycles of approximately the same age suggests a regional (Western Mediterranean) or global control on sea-level change responsible for TCC cycles. In addition, four subaerial exposure-capped TCC cycles may correlate with similar subaerial unconformities in the South Pacific, suggesting at least a partial eustatic control on TCC cyclicity. The high rates of relative sea-level change needed to generate a minimum of 25-30 m sea-level changes associated with each cycle are consistent with glacio-eustacy along with rapid evaporitic drawdown in the Mediterranean. ?? 2001 Elsevier Science B.V. All rights reserved.

Palaeogeography, Palaeoclimatology, Palaeoecology↗

USGS2010

Isotopic composition and origin of indigenous natural perchlorate and co-occurring nitrate in the southwestern United States

Perchlorate (ClO 4 − ) has been detected widely in groundwater and soils of the southwestern United States. Much of this ClO 4 − appears to be natural, and it may have accumulated largely through wet and dry atmospheric deposition. This study evaluates the isotopic composition of natural ClO 4 − indigenous to the southwestern U.S. Stable isotope ratios were measured in ClO 4 − (δ 18 O, Δ 17 O, δ 37 Cl) and associated NO 3 − (δ 18 O, Δ 17 O, δ 15 N) in groundwater from the southern High Plains (SHP) of Texas and New Mexico and the Middle Rio Grande Basin (MRGB) in New Mexico, from unsaturated subsoil in the SHP, and from NO 3 − -rich surface caliche deposits near Death Valley, California. The data indicate natural ClO 4 − in the southwestern U.S. has a wide range of isotopic compositions that are distinct from those reported previously for natural ClO 4 − from the Atacama Desert of Chile as well as all known synthetic ClO 4 − . ClO 4 − in Death Valley caliche has a range of high Δ 17 O values (+8.6 to +18.4 ‰), overlapping and extending the Atacama range, indicating at least partial atmospheric formation via reaction with ozone (O 3 ). However, the Death Valley δ 37 Cl values (−3.1 to −0.8 ‰) and δ 18 O values (+2.9 to +26.1‰) are higher than those of Atacama ClO 4 − . In contrast, ClO 4 − from western Texas and New Mexico has much lower Δ 17 O (+0.3 to +1.3‰), with relatively high δ 37 Cl (+3.4 to +5.1 ‰) and δ 18 O (+0.5 to +4.8 ‰), indicating either that this material was not primarily generated with O 3 as a reactant or that the ClO 4 − was affected by postdepositional O isotope exchange. High Δ 17 O values in ClO 4 − (Atacama and Death Valley) are associated with high Δ 17 O values in NO 3 − , indicating that both compounds preserve characteristics of O 3 -related atmospheric production in hyper-arid settings, whereas both compounds have low Δ 17 O values in less arid settings. Although Δ 17 O variations in terrestrial NO 3 − can be attributed to mixing of atmospheric (high Δ 17 O) and biogenic (low Δ 17 O) NO 3 − , variations in Δ 17 O of terrestrial ClO 4 − are not readily explained in the same way. This study provides important new constraints for identifying natural sources of ClO 4 − in different environments by multicomponent isotopic characteristics, while presenting the possibilities of divergent ClO 4 − formation mechanisms and(or) ClO 4 − isotopic exchange in biologically active environments.

Environmental Science & Technology↗

USGS2007

Geochemical and geophysical examination of submarine groundwater discharge and associated nutrient loading estimates into Lynch Cove, Hood Canal, WA

Geochemical tracer data (i.e., 222 Rn and four naturally occurring Ra isotopes), electromagnetic (EM) seepage meter results, and high-resolution, stationary electrical resistivity images were used to examine the bi-directional (i.e., submarine groundwater discharge and recharge) exchange of a coastal aquifer with seawater. Our study site for these experiments was Lynch Cove, the terminus of Hood Canal, WA, where fjord-like conditions dramatically limit water column circulation that can lead to recurring summer-time hypoxic events. In such a system a precise nutrient budget may be particularly sensitive to groundwater-derived nutrient loading. Shore-perpendicular time-series subsurface resistivity profiles show clear, decimeter-scale tidal modulation of the coastal aquifer in response to large, regional hydraulic gradients, hydrologically transmissive glacial terrain, and large (4-5 m) tidal amplitudes. A 5-day 222 Rn time-series shows a strong inverse covariance between 222 Rn activities (0.5−29 dpm L -1 ) and water level fluctuations, and provides compelling evidence for tidally modulated exchange of groundwater across the sediment/water interface. Mean Rn-derived submarine groundwater discharge (SGD) rates of 85 ± 84 cm d -1 agree closely in the timing and magnitude with EM seepage meter results that showed discharge during low tide and recharge during high tide events. To evaluate the importance of fresh versus saline SGD, Rn-derived SGD rates (as a proxy of total SGD) were compared to excess 226Ra-derived SGD rates (as a proxy for the saline contribution of SGD). The calculated SGD rates, which include a significant (>80%) component of recycled seawater, are used to estimate associated nutrient (NH 4+ , Si, PO 4 3- , NO 3 + NO 2 , TDN) loads to Lynch Cove. The dissolved inorganic nitrogen (DIN = NH 4 + NO 2 + NO 3 ) SGD loading estimate of 5.9 × 10 4 mol d -1 is 1−2 orders of magnitude larger than similar estimates derived from atmospheric deposition and surface water runoff, respectively.

Washington↗

USGS2021

Capacity assessment for Earth Monitoring, Analysis, and Prediction (EarthMAP) and future integrated monitoring and predictive science at the U.S. Geological Survey

Executive Summary Managers of our Nation’s resources face unprecedented challenges driven by the convergence of increasing, competing societal demands and a changing climate that affects the stability, vulnerability, and predictability of those resources. To help meet these challenges, the scientific community must take advantage of all available technologies, data, and integrative Earth systems modeling capacity to better inform resource and risk management decisions. This is the overarching goal of the U.S. Geological Survey (USGS) Earth Monitoring, Analysis, and Prediction (EarthMAP) vision: “By 2030, the USGS will deliver well integrated observations and predictions of the future state of natural systems—water, ecosystems, energy, minerals, hazards—at regional and national scales, working primarily with federal, state, and academic partners to develop and operate the capability” (U.S. Geological Survey, 2021). Providing more integrated Earth systems science and actionable information to decision makers, stakeholders, and the public requires a better understanding of the depth and distribution of existing capacity (capabilities, tools, and techniques) across the Bureau. Identifying existing capacity is also a critical first step toward gap analysis and targeted investments to increase capacity over time. The USGS formed a Capacity Assessment Team (CAT) and charged it with (1) conducting a Request for Information (RFI) to identify existing USGS expertise and activities supportive of integrated and predictive science to inform decision making, (2) developing a strategy and proof-of-concept for a continuously updated capacity assessment capability, and (3) identifying lessons learned to inform development of best practices for future capacity assessment efforts. The RFI took the form of a survey, with content guided by the science and technology needs identified in a USGS report titled “Grand Challenges for Integrated U.S. Geological Survey Science—A Workshop Report” (Jenni and others, 2017). The 44-question survey provided respondents the ability to rate their level of experience with a suite of priority disciplines, analysis and modeling approaches, technologies, and stakeholder engagement strategies and to enter optional narrative text for supporting context. An introductory portion focused on general science capacity assessment, followed by three sections targeting capabilities related to the foundational components of EarthMAP: (1) data and information integration, (2) integrated predictive science, and (3) actionable information. The survey results provided a high-level snapshot of USGS capacity in the targeted areas. Respondents (1,035 individuals) represented approximately 13 percent of the USGS across all mission areas and regions. Seventy-four percent of the respondents held a science-focused position title and the remainder had position titles in information technology, computer science, management, administrative, or other (contractors, volunteers, emeritus, and unknown). To provide greater insight into respondent capabilities and activities, information from the U.S. Department of the Interior and USGS enterprise information systems were used to further characterize topical expertise and organizational associations of survey respondents. To address the ongoing need to assess the Bureau’s capacity to address integrated predictive science priorities, the CAT developed a software-based proof-of-concept called the Integrated Science Assessment Information Database (iSAID) for assembling various information sources together toward making the full extent of USGS capabilities and scientific assets available for routine capacity assessment. This proof-of-concept is intended to serve as a catalyst for further development. The process of implementing the EarthMAP capacity assessment survey, analyzing survey responses, and developing the proof-of-concept resulted in lessons learned, findings, and recommendations. Example scenarios throughout the report demonstrate how capacity assessment data can inform science planning. Three overarching findings and recommendations are: (1) Finding: Capacity is limited in some critical disciplines, skills, and technology applications, but “sufficient” depends on the question and the need relative to availability at a given point in time. Recommendation: Develop an on-demand capacity assessment framework that enables rapid identification and evaluation of existing and available expertise to support decision needs as they arise. (2) Finding: Institutional barriers and lack of awareness constrain the ability of USGS staff to adopt new technologies, collaborate across administrative boundaries, and deliver actionable information to stakeholders in a timely manner. However, these barriers are not universally experienced. Recommendation: Pursue more targeted inquiries to clarify which institutional barriers are obstructing the adoption of new technologies and approaches or the sharing of expertise and equipment across organizational and regional boundaries. These inquiries should inform USGS leadership, mission areas, and regions whether policies can be revised or whether a lack of understanding is creating perceived obstacles. Highlight cases when staff have successfully adopted new technologies and approaches to advance EarthMAP priorities and provide actionable information in a timely manner to spread awareness of how perceived obstacles can be navigated and overcome when appropriate. (3) Finding: Examples of people and projects integrating across disciplines and scales and applying advanced approaches to meet complex stakeholder needs exist. Such examples provide transfer value across the spectrum from approach to decision making. Many projects, already underway, appear to meet elements of the EarthMAP vision, and the USGS has people who can provide leadership in multiple types of specific integrated science efforts. Recommendation: Use these findings as a starting point for near-term strategic planning for integrated science. Highlight, incentivize, and build on existing interdisciplinary predictive science and information delivery activities across the USGS to advance toward further realization of an EarthMAP capacity. The CAT efforts to develop and assess existing USGS capacity to advance the EarthMAP vision revealed a fundamental challenge for not only this effort but any effort to assess existing capacity: A considerable amount of thought, time, and effort is required to survey and assess capabilities and tools available to support a given need, yet best results are still likely to provide an incomplete assessment. To better meet the frequent need to assess capabilities, tools, products, and projects that address an expressed strategic priority, the CAT proposes the concept of an on-demand capacity assessment framework supported by a software package that dynamically pulls and integrates information from existing USGS information systems and public domain registries. Although existing USGS enterprise information systems currently lack the structure, cross-system consistency, interoperability, and stability to support a continuously updated capacity assessment capability, we identify reasonable near-term steps to improve the utility of information gathered on expertise and project capacity and to improve the consistency and completeness of information and the ability of USGS systems to share that information. The ability to search and characterize this information will make future assessments of capacity faster, more complete, more efficient, and more targeted. This approach would grow the Bureau’s capacity knowledge over time, iteratively improving the ability to access, leverage, and synthesize existing capabilities and assets as well as identify and fill critical gaps. The greatest promise for developing integrated science could lie in linking across existing projects and expertise to create a multi-project capacity for addressing large, complex environmental issues.

Open-File Report↗

USGS2025

When the wild things are: Defining mammalian diel activity and plasticity

Circadian rhythms are a mechanism by which species adapt to environmental variability and fundamental to understanding species behavior. However, we lack data and a standardized framework to accurately assess and compare temporal activity for species during rapid ecological change. Through a global network representing 38 countries, we leveraged 8.9 million mammalian observations to create a library of 14,587 standardized diel activity estimates for 445 species. We found that less than half the species’ estimates were in agreement with diel classifications from the reference literature and that species commonly used more than one diel classification. Species diel activity was highly plastic when exposed to anthropogenic change. Furthermore, body size and distributional extent were strongly associated with whether a species is diurnal or nocturnal. Our findings provide essential knowledge of species behavior in an era of rapid global change and suggest the need for a new, quantitative framework that defines diel activity logically and consistently while capturing species plasticity.

Science Advances↗

USGS2026

Small creatures, big challenge: A review of current knowledge and assessment of research needs to inform water withdrawal management for the Duck River's freshwater mussel fauna

The Duck River is a remarkable national resource, boasting outstanding biodiversity (Ahlstedt et al. 2017). It contains 151 fish species, over 60 mussel species, and 22 snail species, many of which are federally listed (Etnier and Starnes 1993; Johnson et al. 2013; Ahlstedt et al. 2017; Womble and Rosenberger 2025a). The Duck River watershed is at a pivotal time in its management and ecological history due to significant growth in human population and economic activity. The consequent increase in surface water demand presents challenges for sustainable water resource management. The Tennessee Department of Environment and Conservation (TDEC) has permitted or proposed to permit eight water utilities to withdraw over 77 million gallons per day (mgd) from the Duck River, representing > 40% increase from previously permitted amounts (USFWS 2025). On November 20, 2024, Tennessee Governor Bill Lee signed Executive Order No 108 to ensure long-term stability of the Duck River watershed through balancing economic growth, water resource management, and environmental conservation. Some primary components of this plan are the creation of a Comprehensive Regional Drought Management Plan, the development of a Habitat Conservation Plan, and the establishment of the Duck River Watershed Planning Partnership (DRWPP). Although the planning partnership is not focused on a sole aspect of the Duck River, but the watershed as a whole, the mussel assemblages of the Duck River have been a particular focus due to their importance in providing ecosystem services, sensitivity to water quality, conservation status, and sensitivity to periodic drought conditions in this system. With multiple state and federal agencies involved in managing the Duck River’s unique water resources, a focused review on how future water management could affect the freshwater mussel communities in the Duck River will contribute to this process. Our objectives of this review are to: 1. Review the effects of changing flow regimes on mussels, with emphasis on the effect of droughts/low flows, highlighting key processes/factors affecting mussels. 2. Highlight critical data gaps that could allow future, evidence-based water management decisions in this system. 3. Propose research topics to address critical data gaps and enable evidence-based water management decisions that consider the needs of freshwater mussels. Freshwater mussels are a critical component of riverine ecosystems, providing ecosystem services including water filtration, nutrient cycling, and habitat provision (Vaughn et al. 2004; Haag 2012) that benefit human communities. Environmental flows, also referred to as ecological or e-flows, refer to the quantity, timing, and quality of water flows necessary to maintain the ecological health of rivers, wetlands, and other aquatic ecosystems (Acreman and Dunbar 2004; Acreman 2016). For freshwater mussels, environmental flows include maintaining a flow regime and associated water quality regimes that are favorable for both juveniles and adults while also seasonally supporting host fish and the survival of settling glochidia (i.e., larval mussels) in suitable environments (Gates et al. 2015). Freshwater mussels are adapted to the natural flow regime of their native streams, including seasonal variability in flows. As sedentary filter feeders, mussels require adequate water flow to provide oxygen and food (Allen and Vaughn 2009); however, they also require periods of low flow to increase access to fish hosts and to ensure juveniles are not displaced (Gates et al 2015). Mussels are sensitive to water quality conditions, including pollutants, sedimentation, and artificial alterations of temperature and pH (Strayer and Malcom 2012). Mussel community characteristics, including presence, abundance, and species richness, generally increase with stream size (e.g., stream order); however, responses to changes in flow conditions are often species-specific, and the magnitude and timing of extreme flow conditions have more predictive power for mussel community response than general changes in flow (Lopez and Vaughn 2021). Complex hydraulic variables, such as relative shear stress and Reynolds number, can be strong predictors of mussel community responses because they combine simple hydraulic variables (e.g., velocity, depth) with stream channel characteristics, such as substrate composition and dynamics. Droughts directly and indirectly affect mussels by influencing their physiology, behavior, reproduction, and ecosystem roles. Mass mussel die-offs during droughts and low water conditions are widely documented (Vaughn et al. 2015; Sousa et al. 2018; Dubose et al. 2019). Droughts can also lead to a shift in mussel community composition, with drought-tolerant species becoming more dominant while more sensitive species experience declines (Galbraith et al, 2010; Lopez et al. 2022; Vaughn and Atkinson 2025). Droughts affect mussels through direct mortality from stranding as well as physiological stress from altered water temperature and quality, reduced habitat and food availability and quality, increased predation and competition, disruption of reproductive cycles, and increased susceptibility to disease and pathogens. Although the effects of droughts, water withdrawals, and low flow conditions on mussels are an increasing area of study, there is a general paucity of data-driven consensus regarding optimal mitigation approaches and strategies (Cushway et al. 2024). Long-term management decisions regarding water withdrawals, particularly during extended or multi-year drought conditions, may consider factors beyond direct mortality such as the sublethal effects of drought and low flow conditions on freshwater mussels for a comprehensive understanding of the effects of prolonged low flows on freshwater mussels and aquatic communities.

Tennessee↗

USGS2022

Pharmaceutical pollution of the world’s rivers

Environmental exposure to active pharmaceutical ingredients (APIs) can have negative effects on the health of ecosystems and humans. While numerous studies have monitored APIs in rivers, these employ different analytical methods, measure different APIs, and have ignored many of the countries of the world. This makes it difficult to quantify the scale of the problem from a global perspective. Furthermore, comparison of the existing data, generated for different studies/regions/continents, is challenging due to the vast differences between the analytical methodologies employed. Here, we present a global-scale study of API pollution in 258 of the world’s rivers, representing the environmental influence of 471.4 million people across 137 geographic regions. Samples were obtained from 1,052 locations in 104 countries (representing all continents and 36 countries not previously studied for API contamination) and analyzed for 61 APIs. Highest cumulative API concentrations were observed in sub-Saharan Africa, south Asia, and South America. The most contaminated sites were in low- to middle-income countries and were associated with areas with poor wastewater and waste management infrastructure and pharmaceutical manufacturing. The most frequently detected APIs were carbamazepine, metformin, and caffeine (a compound also arising from lifestyle use), which were detected at over half of the sites monitored. Concentrations of at least one API at 25.7% of the sampling sites were greater than concentrations considered safe for aquatic organisms, or which are of concern in terms of selection for antimicrobial resistance. Therefore, pharmaceutical pollution poses a global threat to environmental and human health, as well as to delivery of the United Nations Sustainable Development Goals.

Proceedings of the National Academy of Sciences↗

USGS2013

Borehole-explosion and air-gun data acquired in the 2011 Salton Seismic Imaging Project (SSIP), southern California: description of the survey

The Imperial and Coachella Valleys are being formed by active plate-tectonic processes. From the Imperial Valley southward into the Gulf of California, plate motions are rifting the continent apart. In the Coachella Valley, the plates are sliding past one another along the San Andreas and related faults (fig. 1). These processes build the stunning landscapes of the region, but also produce damaging earthquakes. Rupture of the southern section of the San Andreas Fault (SAF), from the Coachella Valley to the Mojave Desert, is believed to be the greatest natural hazard that California will face in the near future. With an estimated magnitude between 7.2 and 8.1, such an event would result in violent shaking, loss of life, and disruption of infrastructure (freeways, aqueducts, power, petroleum, and communication lines) that might bring much of southern California to a standstill. As part of the nation’s efforts to avert a catastrophe of this magnitude, a number of projects have been undertaken to more fully understand and mitigate the effects of such an event. The Salton Seismic Imaging Project (SSIP), funded jointly by the National Science Foundation (NSF) and the U.S. Geological Survey (USGS), seeks to understand, through seismic imaging, the structure of the Earth surrounding the SAF, including the sedimentary basins on which cities are built. The principal investigators (PIs) of this collaborative project represent the USGS, Virginia Polytechnic Institute and State University (Virginia Tech), California Institute of Technology (Caltech), Scripps Institution of Oceanography (Scripps), University of Nevada, Reno (UNR), and Stanford University. SSIP will create images of underground structure and sediments in the Imperial and Coachella Valleys and adjacent mountain ranges to investigate the earthquake hazards posed to cities in this area. Importantly, the images will help determine the underground geometry of the SAF, how deep the sediments are, and how fast earthquake energy can travel through the sediments. All of these factors determine how hard the earth will shake during a major earthquake. If we can improve on our understanding of how and where earthquakes will occur, and how strong their resultant shaking will be, then buildings can be designed or retrofitted accordingly in order to resist damage and collapse, and emergency plans can be adequately prepared. In addition, SSIP will investigate the processes of rifting and magmatism in the Salton Trough in order to better understand this important plate-boundary region. The Salton Trough is a unique rift in that subsidence is accompanied by huge influxes of infilling sediment from the Colorado River. Volcanism that accompanies the subsidence here is muted by these influxes of sediment. The Salton Trough, in the central part of the Imperial Valley, is apparently made up of entirely new crust: young sediment in the upper crust and basaltic intrusive rocks in the mid-to-lower crust (Fuis and others, 1984). Similar to the ultrasound and computed tomography (CT) scans performed by the medical industry, seismic imaging is a collection of techniques that enable scientists to obtain a picture of what is underground. The petroleum industry routinely uses these techniques to search for oil and gas at relatively shallow depths; however, the scope of this project demanded that we image as much as 30 km into the Earth’s crust. This project generated and recorded seismic waves, similar to sound waves, which move downward into the Earth and are bent (refracted) or echoed (reflected) back to the surface. SSIP acquired data in a series of intersecting lines that cover key areas of the Salton Trough. The sources of sound waves were detonations (shots) in deep boreholes, designed to create energy equivalent to magnitude 1–2 earthquakes. The study region routinely experiences earthquakes of these magnitudes, but earthquakes are not located in such a way as to permit us to create the detailed images we need for earthquake hazard assessment. Air gun bursts, generated in the Salton Sea along extensions of our onshore seismic lines, also were utilized as sound-wave sources. Temporary deployments of portable land seismometers, as well as ocean-bottom seismometers (OBSs) on the floor of the Salton Sea, recorded the energy from the land shots and air gun bursts. SSIP is similar to the Los Angeles Regional Seismic Experiments of 1994 and 1999 (LARSE I and II, respectively; Murphy and others, 1996; Fuis and others, 2001). The LARSE surveys demonstrated that the USGS and collaborators can safely and effectively conduct seismic imaging surveys in urban and nonurban areas, on lands owned and/or managed by many different types of agencies and entities. Information was produced that could not have been obtained any other way, and this information was key to changing the leading ideas about earthquake hazards at that time in the Los Angeles region. These surveys produced no significant environmental impact or damage to structures, and they did not trigger earthquakes.

California↗

USGS2025

Late Quaternary environmental change in eastern Beringia

Eastern Beringia (Alaska and western Yukon) is an extensive, high-latitude region of North America that remained largely unglaciated throughout the Quaternary. Consequently, its sedimentary deposits preserve long-term environmental records that have intrigued scientists for nearly a century. Recent advances in palaeoecological proxies and dating methods have proved critical in addressing long-standing questions about regional late Quaternary environmental change. At the same time, they have led to new and sometimes controversial hypotheses. This review covers recent discoveries and unresolved questions focused on the period 57,000–10,000 calendar years before C.E. 1950 (cal yr BP). The middle Wisconsin interstadial (57,000–30,000 cal yr BP) was a period of relative warmth in eastern Beringia, compared with the late Wisconsin (30,000–14,000 cal yr BP). Early in the interstadial occasional Picea woodland was present amongst widespread shrub tundra. Palaeoecological, sedimentary and isotopic data indicate that climate was cooler and drier than the Holocene, with high rates of aeolian activity. Megafauna typically associated with the ‘mammoth steppe’ ecosystem (woolly mammoth [ Mammuthus primigenius ], horse [ Equus ] and steppe-bison [ Bison priscus ]) were present in some abundance. The transition towards late Wisconsin cold-stage conditions (35,000–30,000 cal yr BP) coincided with the establishment of the Bering Land Bridge and featured expansion of spatially varied, herbaceous vegetation, sometimes associated with deep active layers. Sedimentary DNA ( seda DNA) and macrofossil evidence show vegetation was not a prairie-like grassland, and the term “steppe-tundra” is a better descriptor. Permafrost pore-ice isotopic (δ 18 O) records suggest a step change in one or more climate drivers ca. 30,000 cal yr BP, by which time steppe-tundra was established across eastern Beringia. It remains uncertain whether Picea survived cold-stage conditions within isolated refugia, or whether it recolonized from south of the Laurentide-Cordilleran ice sheets. Genetic data suggest that Picea probably survived in situ ; however, there is no definitive fossil evidence to support this. The end-Pleistocene transition from steppe-tundra to shrub tundra began ca. 15,000 cal yr BP and took place within decades at local scales. The expansion of woody taxa coincided with rising sea levels, reduced sea-ice extent and an abrupt shift in atmospheric circulation that enhanced precipitation. During this time, Earth's orbital configuration caused high early-summer temperatures and strong seasonality, creating growing conditions very different from today. The vegetation consisted of Salix and Betula shrub tundra with open areas of herbs and graminoids. During the deglacial warming trend, the Younger Dryas oscillation (12,800–11,700 cal yr BP) was variably expressed. It is generally evident in records affected by adjacent oceans but can be absent at sites in continental areas. These past conditions and paleoenvironmental changes have implications for contemporary issues: hypotheses about Pleistocene mammalian extinction; sensitivity of eastern Beringia to major oceanic reorganizations and high-frequency climate variability; the nature of woody plant expansion with climate warming; grazing, hydroclimate and fire as controls over ecosystems; the efficacy of “Pleistocene rewilding” for carbon capture.

Beringia↗

USGS2021

Synoptic analysis of a decade of daily measurements of SO2 emission in the troposphere from volcanoes of the global ground-based Network for Observation of Volcanic and Atmospheric Change

Volcanic plumes are common and far-reaching manifestations of volcanic activity during and between eruptions. Observations of the rate of emission and composition of volcanic plumes are essential to recognize and, in some cases, predict the state of volcanic activity. Measurements of the size and location of the plumes are important to assess the impact of the emission from sporadic or localized events to persistent or widespread processes of climatic and environmental importance. These observations provide information on volatile budgets on Earth, chemical evolution of magmas, and atmospheric circulation and dynamics. Space-based observations during the last decades have given us a global view of Earth's volcanic emission, particularly of sulfur dioxide ( SO 2 ). Although none of the satellite missions were intended to be used for measurement of volcanic gas emission, specially adapted algorithms have produced time-averaged global emission budgets. These have confirmed that tropospheric plumes, produced from persistent degassing of weak sources, dominate the total emission of volcanic SO 2 . Although space-based observations have provided this global insight into some aspects of Earth's volcanism, it still has important limitations. The magnitude and short-term variability of lower-atmosphere emissions, historically less accessible from space, remain largely uncertain. Operational monitoring of volcanic plumes, at scales relevant for adequate surveillance, has been facilitated through the use of ground-based scanning differential optical absorption spectrometer (ScanDOAS) instruments since the beginning of this century, largely due to the coordinated effort of the Network for Observation of Volcanic and Atmospheric Change (NOVAC). In this study, we present a compilation of results of homogenized post-analysis of measurements of SO 2 flux and plume parameters obtained during the period March 2005 to January 2017 of 32 volcanoes in NOVAC. This inventory opens a window into the short-term emission patterns of a diverse set of volcanoes in terms of magma composition, geographical location, magnitude of emission, and style of eruptive activity. We find that passive volcanic degassing is by no means a stationary process in time and that large sub-daily variability is observed in the flux of volcanic gases, which has implications for emission budgets produced using short-term, sporadic observations. The use of a standard evaluation method allows for intercomparison between different volcanoes and between ground- and space-based measurements of the same volcanoes. The emission of several weakly degassing volcanoes, undetected by satellites, is presented for the first time. We also compare our results with those reported in the literature, providing ranges of variability in emission not accessible in the past. The open-access data repository introduced in this article will enable further exploitation of this unique dataset, with a focus on volcanological research, risk assessment, satellite-sensor validation, and improved quantification of the prevalent tropospheric component of global volcanic emission.

Earth System Science Data↗

USGS2019

Case studies in groundwater contaminant fate and transport

A case study of groundwater contamination is a detailed study of a single site contaminated with a chemical or mixture that is known to be a problem at many sites. The goal of case studies is to provide insights into the physical, chemical, and biological processes controlling migration, natural attenuation, or remediation of common groundwater contaminants. Ideally, processes occurring at a case study site are representative of other sites so that knowledge gained from these intensive studies can be applied at thousands of sites where fewer data are available. Several characteristics of case studies contribute to their value. First, they may have tens to hundreds of monitoring wells, compared to fewer than ten wells at some contaminated sites. Second, some case studies continue for many years or even decades, providing insights into temporal progression of slow processes. Third, analytical methods prohibitively expensive for routine use or under development may be tested at case study sites. Finally, the ongoing characterization typical of case study sites builds a foundation of knowledge that facilitates sophisticated experimental design and testing of new methods. This article is divided into sections based on the contaminant type because the chemical and biological processes required for remediation vary for each contaminant. Most importantly, some contaminants can be biodegraded whereas metals and radionuclides cannot be destroyed but can be immobilized or rendered less toxic. The emphasis is on case studies of natural processes that control the fate and transport of contaminants in groundwater rather than on active remediation methods. The principles learned from these studies may form the basis for design of remedial strategies. The organic contaminants are divided into: petroleum hydrocarbons, fuel oxygenates, coal tar and wastes from manufactured gas plants, and chlorinated solvents. The inorganic contaminants covered are metals and radionuclides, arsenic, and nitrate. Case studies of mixed waste plumes from landfills are also described. Experimental sites where contaminants have been introduced into an aquifer as an emplaced source or a controlled release may not meet the above definition of case studies, but some are included because the overall goal is to impart lessons learned from detailed field studies. It is impossible to cover all case studies in this short format. Conversely, focusing on one or two does not convey the breadth of research results in entire range of case studies. Instead, the strategy is to describe the evolution of knowledge for each contaminant class while providing citations of relevant case studies. Much of the progress in understanding of the fate of contaminants in groundwater is based on laboratory studies; thus whenever possible, papers that included both field and laboratory results have been included among the citations. Two topics of growing importance have not been covered. These are the fate of pharmaceuticals in groundwater and discharge of contaminant plumes to surface water. These topics merit coverage in the future as knowledge grows and case studies increase in number.

Environmental Science↗

USGS2020

Population and habitat analyses for greater sage-grouse (Centrocercus urophasianus) in the bi-state distinct population segment—2018 update

Executive Summary The Bi-State Distinct Population Segment (Bi-State DPS) of greater sage-grouse ( Centrocercus urophasianus , hereinafter “sage-grouse”) represents a genetically distinct and geographically isolated population that straddles the border between Nevada and California. The primary threat to these sage-grouse populations is the expansion of single-leaf pinyon ( Pinus monophylla ) and Utah juniper ( Juniperus osteosperma ) into sagebrush ecosystems, which fragments and reduces population connectivity and survival. Other important threats include low water availability during brood-rearing, particularly during drought, and increased predation by common ravens ( Corvus corax ), a generalist predator often associated with anthropogenic resource subsidies. Although the Bi-State DPS occurs at high elevations relative to sage-grouse range-wide, changes in historical wildfire cycles and the conversion of native shrubs to invasive annual grasslands still threaten these populations. The Bi-State DPS has undergone multiple federal status assessments and associated litigation. For example, in October of 2013, the Bi-State DPS was proposed for listing as threatened under the Endangered Species Act of 1973 by the U.S. Fish and Wildlife Service (USFWS), then withdrawn in April 2015. The withdrawal decision was challenged, and in May 2018, a Federal district court ordered the withdrawal decision to be vacated, and USFWS was required to re-open the October 2013 listing evaluation. In response, the U.S. Geological Survey (USGS), with State and Federal collaborators, embarked on a multipronged analysis to provide current and best available science regarding population status of sage-grouse within the Bi-State DPS. Using data from a long-term monitoring program, we carried out four analytical study objectives, and here, we provide preliminary results of these analyses. First, we used integrated population modeling (IPM) to predict annual population abundance and annual finite rate of population change for the Bi-State DPS, as a whole, and for each subpopulation between 1995 and 2018. Because sage-grouse exhibit population cycles (periodic increases and decreases in abundance across approximately 6- to 10-year wavelengths), we estimated trends across three nested temporal scales that represent one (11 years), two (18 years), and three (24 years) complete population cycles. These estimates of relatively long-term averaged population change account for temporal (that is, interannual) variation. Our model predicted population abundance for the Bi-State DPS during 2018 at 3,305 individuals (2,247–4,683), with the majority occupying Bodie Hills and Long Valley. The model also predicted cyclic dynamics in abundance through time with evidence of 24-year population growth and slight trends of decline over the past 18 years. Specifically, across the Bi-State DPS as a whole, we estimated annual average at 0.99, 0.99, and 1.02 over the one, two, and three population cycles, which equated to a 10.5 percent, 16.6 percent decrease, and 60.0 percent increase in abundance over the 11-, 18-, and 24-year cycles. Estimated abundance in 2018 had not reached numbers lower than those predicted during 1995. However, we found spatial variation in population trends across the three cycles. Bodie Hills subpopulation comprised the greatest (1,521) and exhibited average annual greater than 1.0 across all periods resulting in average annual increases of 7 percent. This relatively large subpopulation has grown 5 times larger than what was predicted in 1995 while experiencing cyclical dynamics within that period. Conversely, other smaller subpopulations within the Bi-State DPS exhibited average annual equal to or less than 1.0 resulting in estimated 10-year risks of extirpation ranging from 2.0 to 76.1 percent. In general, evidence of decline among smaller subpopulations was greatest for the most recent period (2008–18) compared to a period that encompassed three full population cycles (24-year). This difference coincides with an intense period of drought that began in 2012. For comparative purposes as part of this first objective, we conducted a similar analysis for populations of sage-grouse within Nevada and California but outside the Bi-State DPS. We developed a region-wide and distance-weighted IPM using lek count from Nevada Department of Wildlife (NDOW) and California Department of Fish and Wildlife (CDFW) databases and with telemetry data collected by USGS across 12 sage-grouse subpopulations. Our models predicted similar patterns in population cycling outside the Bi-State DPS but with much stronger evidence of long-term declines across 24 years. Specifically, median averaged across each year of the 11-, 18-, and 24-year periods resulted in average annual values of 0.94, 0.97, and 0.99, respectively. These values equate to 41.0 percent, 38.5 percent, and 21.3 percent declines over the corresponding periods. Second, we used lek count data in a state-space modeling framework to compare trends in population abundance across different spatial scales (that is, leks versus Bi-State DPS). This hierarchical framework allowed us to disentangle declines associated with climate conditions as opposed to other local level factors that might signal the need for management intervention. Specifically, we identified 7 leks that were both declining and recently decoupled from larger spatial scale trends, typically governed by climatic conditions (referred to as soft or hard signals). The goal of this analysis was to provide an early warning system that might have implications for conservation actions at local scales. Third, we developed phenological (spring, summer–fall, and winter) and reproductive life stage (nesting, early brood-rearing, and late-brood rearing) based resource selection functions using various environmental covariates. We report rankings of variable importance for each season and life stage, developed maps of habitat selection indices (HSI), binned categories representing low, moderate, and high classes of quality (where any category greater than or equal to low indicated selected habitat) for each phenological season and life stage, and produced composite maps by selected phenological and reproductive stage to estimate annual habitat. Fourth, we used for each lek within the Bi-State DPS to carry out a spatial analysis that quantified substantial changes in the distribution of occupied habitat across long- (24-year) and short- (11-year) term periods. Owing to differences among available datasets, the long-term analysis primarily reflected spatial shifts among subpopulations comprising the majority of the Bi-State DPS (that is, Bodie Hills and Long Valley) while the short-term analysis also quantified changes among subpopulations along the periphery. Over long and short-term periods, the overall distribution of occupied habitat (as measured by 99 percent utilization distributions intersecting any quantified habitat) was reduced by 20,573 ha and 55,492 ha, respectively. Occupied core areas (as measured by 50 percent utilization distributions intersecting any quantified habitat) over long-term periods were solely located in Bodie Hills and Long Valley. Although nearly all subpopulations experienced contractions in occupied overall and core distribution, Bodie Hills experienced spatial expansion that occurred with concomitant spatial contraction at Long Valley over both periods. Subpopulations at the northern (Pine Nuts), central (Sagehen) and southern (White-Mountains) extents of the Bi-State DPS also experienced spatial contraction over the short-term period. These findings, coupled with those of population trends, indicate long-term patterns in redistribution of sage-grouse from Long Valley and peripheral subpopulations to Bodie Hills. That is, sage-grouse subpopulations at the periphery are declining while the largest population at the core is increasing, which could have meaningful impacts on overall metapopulation persistence. We provide evidence for loss of occupied habitat (reduced distribution) given local extirpation of subpopulations. Fifth, we calculated percentages of selected phenological, life stage, and annual habitat that each subpopulation contributed to the Bi-State DPS. We then intersected these maps with a composite estimate of occupied habitat from the fourth objective and calculated percentages of selected habitat likely occupied by sage-grouse that each subpopulation contributed to the Bi-State DPS. These values provide evidence for loss of occupied habitat and subsequent reductions in spatial distribution given reductions in abundance and, in some cases, extirpation of leks within subpopulations. Lastly, we carried out an initial in-depth analysis of selection for irrigated pastures and wet meadows during the brood-rearing stage for the Long Valley subpopulation. We chose this subpopulation because it represents a population core, representing 26.5 percent of total sage-grouse within the Bi-State DPS, and has exhibited long-term declines in abundance and distribution. This subpopulation is highly sensitive to precipitation and other factors that influence water availability. Models predicted higher use of the interior portions of irrigated pastures and wet meadows during late brood-rearing period, which represented a potentially risky use of habitat that was exacerbated during periods of low moisture (for example, drought, reduced water delivery, or both). Sage-grouse typically used edges of riparian areas and pastures, largely because the interior of these mesic areas consisted of considerably less overhead concealment cover (for example, shrubs) that likely resulted in a higher risk of mortality. We found that a lack of water delivery to pastures in the form of overwinter precipitation or diversion ditches increased the movements of sage-grouse to the interior of pastures. Although further investigation of water delivery impacts on chick survival are warrented, our initial findings regarding resource selection may explain recent declines in population growth at Long Valley.

California, Nevada↗