USGS Science⌕ Search

SEARCH · USGS Science

Results for “Avian Conservation and Ecology”

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

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

1,656 records · Page 46Linked to original sources
USGS2025

Using public participatory geographic information systems (PPGIS) to explore uses and values for Mojave Trails National Monument, California

Many people ascribe a variety of values to public lands and waters, but some values are more difficult to assess and quantify than others. Public participatory geographic information systems (PPGIS) are tools that have been used to help quantify and map the public’s diverse values for a landscape. This work describes the first known Office of Management and Budget–approved use of PPGIS by a Department of the Interior bureau. The U.S. Geological Survey developed an internet-based application to aid in gathering PPGIS data, called Values Mapping for Planning in Regional Ecosystems (VaMPIRE). Further, this work describes the first pilot of the VaMPIRE application in coordination with the Bureau of Land Management to collect spatial data and other survey data regarding the public’s uses of and values for locations within Mojave Trails National Monument. We emailed the link to the VaMPIRE application to an interested party email list in 2024 with 207 valid emails and received 74 responses; we also received 47 responses from members of an off-roading social media group. Of the list of 16 value options, recreation was the most popular value for the monument, followed by wilderness and inspirational. Over 1,000 points were placed throughout the monument, indicating locations people use or value, with the locations spread throughout the entire monument. Additionally, most survey respondents stated their ability to receive benefits in locations they mapped would not change in response to a hypothetical scenario related to recreational facility development. This report describes exploratory results from the first use of the VaMPIRE tool in Mojave Trails National Monument and includes reflections on how the process went and considerations for future use of VaMPIRE.

California↗
USGS2024

Application of non-stationary shear-wave velocity randomization approach to predict 1D seismic site response and its variability at two downhole array recordings

Accounting for uncertainties in seismic site response is crucial to improving the performance of one-dimensional (1D) ground response analyses (GRAs) at downhole array recording sites. In addition to site effects, uncertainties in 1D-GRAs can also be contributed from the seismic source and/or path. Though often representing not more than one percent of the distance (path) from the source, site conditions are known to have an enormous influence on ground shaking. In this study, we focus on the site shear-wave velocity ( V S ) structure, which is the main ingredient for estimating the variability of site response. As such, V S can manifest aleatory uncertainties related to the effects of small-scale spatial heterogeneities within the near surface, thus V S can substantially modify ground shaking during earthquakes. We apply a novel V S randomization approach to propagate the small-scale heterogeneities of V S to estimate seismic site response within a non-stationary probabilistic framework. The randomization approach generates samples of V S profiles that are used to perform several 1D-GRAs and obtain an averaged site response and related variability. The proposed method is implemented on data recorded at two downhole array sites with different subsurface soil conditions: a soft soil site on Treasure Island (California, United States of America) and a rock outcrop site in Cadarache (South-East France). We show that synthetic surface-to-borehole transfer functions from 1D-GRAs provide an acceptable fit to the empirical transfer functions from low-motion earthquake records and succeed in reproducing most of the site-specific seismic response variability. The remaining mismatch between transfer functions is likely due to insufficient precision on the seismic bedrock and the impedance contrast. The variability in site response is discussed with emphasis on the role of V S small-scale heterogeneities, attenuation, and input motion incidence angle in ground motion variability for the site and soil conditions at both locations.

California↗
USGS2025

Fiber-optic distributed temperature sensing of hydrologic processes—Diverse deployments and new applications by the U.S. Geological Survey

Fiber-optic distributed temperature sensing instruments harness the temperature-dependent properties of glass to measure temperature continuously along optical fibers by using precise pulses of laser light. In the mid-2000s, this technology was refined for environmental monitoring purposes such as snowpack-air exchange, groundwater/surface-water exchange, and lake-water stratification. Fiber-optic distributed temperature sensing has revealed unprecedented details about preferential flow processes; however, the method is labor intensive and requires specific training, resulting in limited use by the broader water community. With the ongoing national implementation of the U.S. Geological Survey Next Generation Water Observing System, there has been renewed interest in harnessing the unique spatiotemporal monitoring capabilities of fiber-optic distributed temperature sensing. This fact sheet briefly describes this technology, highlights uses by the U.S. Geological Survey, and discusses current applications and future opportunities.

Fact Sheet↗
USGS2025

Quantitative evaluations of earthquake early warning performance using “Did You Feel It?” and post-alert surveys

We examine responses to the U.S. Geological Survey’s “Did You Feel It?” (DYFI) survey and its companion earthquake early warning (EEW) questionnaire to assess the performance of the U.S. ShakeAlert EEW system directly from the alert recipients’ perspectives. ShakeAlert rapidly detects earthquakes and develops alert information, but as official alert delivery partners issue these alerts, it is thus difficult to determine how many people were alerted and when. We investigate DYFI reports for six California earthquakes that had EEW alerts and substantial responses to the DYFI EEW questionnaire. Comparisons of ShakeAlert predictions to reported intensities demonstrate that magnitude estimation accuracy is not necessarily indicative of ground-motion prediction accuracy. Perceived warning time distributions indicate that estimating maximum-expected warning times using the S-wave arrival is a reasonable assumption when discussing public EEW performance. However, we also find many reports of shorter warning times, late alerts, and missed alerts than expected based on ShakeAlert publication times, suggesting alert delivery latencies are substantial and highly variable. The novelty of our analysis is that we demonstrate that the DYFI EEW survey provides useful EEW efficacy information—independent of the specific alerting pathway—that can be used to inform our choices for conveying EEW performance.

California↗
USGS2024

Testing rate‐and‐state predictions of aftershock decay with distance

We analyze aftershocks of the 2019 M 7.1 Ridgecrest mainshock and isolated M 5–6 mainshocks in southern California to test predictions made by the rate‐and‐state friction model of Dieterich (1994) . Rate‐and‐state friction predicts that the seismicity rate after a stress step follows Omori decay, where the Omori c ‐value, which is the saturation in aftershock rate observed at small times, is larger for smaller stress steps. Put in the context of an aftershock sequence, this predicts that the Omori c ‐value will be systematically larger at greater distances from the mainshock. To our knowledge, this predicted effect has not been observed. In part this may be because the Omori c ‐value is difficult to measure because it often reflects short‐term catalog incompleteness rather than a true saturation in aftershock rate. We explore the dependence of the Omori c ‐value on the distance to the mainshock by applying the “ a ‐positive” method ( van der Elst and Page, 2023 ). This method is insensitive to short‐term aftershock incompleteness and allows resolution of the true aftershock rate deep into the mainshock coda. For aftershocks of the Ridgecrest mainshock and stacked M 5–6 mainshocks, we observe systematic differences in early aftershock rates, relative to mainshock distance, consistent with the predictions of rate‐and‐state friction. Furthermore, for the larger Ridgecrest dataset, we observe that aftershocks nearer to the mainshock start earlier, and we resolve a flattening of the Omori curve consistent with a larger Omori c ‐value for the farthest aftershocks, as predicted by Dieterich (1994) .

Seismological Research Letters↗
USGS2026

Surface-wave relocation and characterization of the 2024–2025 dyking episode along the Fentale–Dofen segment of the Ethiopian rift

From late-December 2024 to mid-March 2025, a 50-km-long dyke intrusion triggered over 300 earthquakes (magnitude 4 to 5.9) between Fentale and Dofen volcanoes along the Northern Main Ethiopian Rift. Dyke intrusions periodically occur along the Fentale–Dofen magmatic segment and are an expression of ongoing rift extension. Preliminary analyses using interferometric synthetic aperture radar revealed extensive ground deformation (up to 60 cm), which closely matched the temporal and spatial evolution of surface manifestations and earthquake locations from global catalogues. While global catalogues are critical for real-time monitoring, the precision of locations in remote and or sparsely instrumented regions can be low. In this investigation, we present surface-wave relocation results of the dyking episode that began near Fentale volcano in December 2024. We estimate relative locations using differential traveltimes measured from regional-to-teleseismic distance surface-wave observations of earthquakes reported by the U.S. Geological Survey. Relative relocations reduce the initial region of diffuse seismicity to a 50-km-long narrow band bounding the strike of surface manifestations and the zone of maximum surface deformation. We demonstrate the precision of surface-wave relocations over incremental time periods, capturing the progression of dyking from seismic onset through seismic migration and caldera subsidence. Results showcase the utility of surface-wave relocations in the characterization of dyking episodes and provide complementary insights into the current understanding of the Fentale–Dofen volcanic plumbing system.

Geophysical Journal International↗
USGS2025

A global database of soil microbial phospholipid fatty acids and enzyme activities

Soil microbes drive ecosystem function and play a critical role in how ecosystems respond to global change. Research surrounding soil microbial communities has rapidly increased in recent decades, and substantial data relating to phospholipid fatty acids (PLFAs) and potential enzyme activity have been collected and analysed. However, studies have mostly been restricted to local and regional scales, and their accuracy and usefulness are limited by the extent of accessible data. Here we aim to improve data availability by collating a global database of soil PLFA and potential enzyme activity measurements from 12,258 georeferenced samples located across all continents, 5.1% of which have not previously been published. The database contains data relating to 113 PLFAs and 26 enzyme activities, and includes metadata such as sampling date, sample depth, and soil pH, total carbon, and total nitrogen. This database will help researchers in conducting both global- and local-scale studies to better understand soil microbial biomass and function.

Scientific Data↗
USGS2024

Remote sensing evapotranspiration in ensemble-based framework to enhance cascade routing and re-infiltration concept in integrated hydrological model applied to support decision making

Integrated hydrological models (IHMs) help characterize the complexity of surface–groundwater interactions. The cascade routing and re-infiltration (CRR) concept, recently applied to a MODFLOW 6 IHM, improved conceptualization and simulation of overland flow processes. The CRR controls the transfer of rejected infiltration and groundwater exfiltration from upslope areas to adjacent downslope areas where that water can be evaporated, re-infiltrated back to subsurface, or discharged to streams as direct runoff. The partitioning between these three components is controlled by uncertain parameters that must be estimated. Thus, by quantifying and reducing those uncertainties, next to uncertainties of the other model parameters (e.g. hydraulic and storage parameters), the reliability of the CRR is improved and the IHM is better suited for decision support modelling, the two key objectives of this work. To this end, the remotely sensed MODIS-ET product was incorporated into the calibration process for complementing traditional hydraulic head and streamflow observations. A total of approximately 150,000 observations guided the calibration of a 13-year MODFLOW 6 IHM simulation of the Sardon catchment (Spain) with daily stress periods. The model input uncertainty was represented by grid-cell-scale parameterization, yielding approximately 500,000 unknown input parameters to be conditioned. The calibration was carried out through an iterative ensemble smoother. Incorporating the MODIS-ET data improved the CRR implementation, and reduced uncertainties associated with other model parameters. Additionally, it significantly reduced the uncertainty associated with net recharge, a critical flux for water management that cannot be directly measured and rather is commonly estimated by IHM simulations.

Sardon catchment↗
USGS2026

Regional characterization of coal resources in the Gulf Coast Coal Region—In support of coal as a potential source for critical minerals and rare earth elements

There is increasing interest in extracting critical minerals (CM), including rare earth elements (REE), from coals in the United States to address the overreliance on imported REE. The Gulf Coast and the Fort Union regions are the two major lignite-bearing regions within the United States. Recent REE and CM studies of the lignite in these regions have indicated that the coals may be a viable source material for REE and CM extraction. To evaluate in-place coal as a potential source of REE and CM, the coal resources need to be quantified. This study presents the results of a regional analysis of Gulf Coast Coal Region lignite and bituminous coal resources that might be available as potential sources of REE and CM. The resource analysis used kriging methods to develop isopleth maps of cumulative coal thickness throughout the region using data from drill holes and measured mine and outcrop sections for a total of 31,181 data points. The estimated total remaining coal resource in the Gulf Coast Coal Region is about 83 billion metric tons in the upper 90 meters (approximately 300 feet) of the subsurface. The current estimate of remaining coal resources includes approximately 1.88 billion metric tons of coal that have been extracted in the region since 1983. Texas accounted for 39 percent (32 billion metric tons) of the remaining total resource, followed by Mississippi (24 percent, 20 billion metric tons ), Louisiana (14 percent, 12 billion metric tons), Tennessee (10 percent, 8.5 billion metric tons), and Arkansas (6 percent, 5.1 billion metric tons). Alabama, Georgia, Kentucky, and Missouri each accounted for less than 5 percent of the total resource; Georgia had the smallest resource estimated at 7 million metric tons. Here we report the first known statewide lignite resource estimates for Georgia, Kentucky (820 million metric tons), and Missouri (1,800 million metric tons). A comparison of the results of this study with those of previous Gulf Coast Coal Region and Fort Union Coal Region resource studies is difficult because each study used different data sources, assessment methodologies, overburden depths, and qualifying coal thicknesses.

Alabama, Arkansas, Florida, Georgia, Illinois, Lou↗
USGS2025

The influence of pre-existing structures on geothermal springs: Inferences from potential field mapping in Surprise Valley, CA and other sites In the northwestern Great Basin

Surprise Valley, located in the northwestern Great Basin, is an asymmetric extensional basin that marks a major tectonic transition between the relatively un-extended volcanic Modoc Plateau to the west, and the Basin and Range to the east that has undergone 10-15% extension. In addition, it sits just north of the Walker Lane which accommodates up to 20% of dextral slip associated with Pacific-North American plate interactions. Thermal springs issue from eight areas within Surprise Valley. Most of these occur within the basin and are not situated on the main basin forming range-front faults. As a result, efforts to resolve the structural setting of the valley’s hydrothermal system have relied on geophysics to characterize basin structure and geology. Extensive efforts to map the basin with ground and airborne magnetics have revealed a >35 km-long linear, intra-basin magnetic high, interpreted as a buried dike swarm. Geothermal springs on the eastern side of the valley, including Seifert hot springs, Leonards hot springs, and Surprise Valley hot springs (SVHS), are all situated along the magnetic high and occur at local breaks and bends in the anomaly, suggesting that fracture permeability is enhanced along the feature and particularly at these discontinuities. Recent studies, including drilling over the anomaly near SVHS that likely intersected dike material, as well as subsequent mapping and sampling of dikes outcropping along the anomaly on the playa surface south of SVHS, confirm (as previously inferred) that mafic intrusives are the principal source of the anomaly. Similar interpretations made in two other valleys (in southern Oregon and northwestern Nevada), where inferred intra-basin dikes appear to be spatially correlated with hot springs or prospective geothermal resource areas, suggest that the impact of pre-existing basement structure on hydrothermal activity may pertain more generally to other hydrothermal settings throughout the Great Basin. If so, efforts to map basement may enhance understanding structural controls on some geothermal systems. Furthermore, similarities across these disparate sites suggest that magmatism may play a much larger role in accommodating extension and influencing basin evolution across the western Great Basin than previously recognized.

Conference Paper↗
USGS2024

Post-glacial stratigraphy and late Holocene record of great Cascadia earthquakes in Ozette Lake, Washington, USA

Ozette Lake is an ~100-m-deep coastal lake located along the outer coast of the Olympic Peninsula (Washington, USA); it is situated above the locked portion of the northern Cascadia megathrust but also relatively isolated from active crustal faults and intraslab earthquakes. Here we present a suite of geophysical and geological evidence for earthquake-triggered mass transport deposits (MTDs) and related turbidite deposition in Ozette Lake since ca. 14 ka. Comprehensive high-resolution bathymetry data, seismic reflection profiles, and sediment cores are used to characterize the post-glacial stratigraphic framework and examine paleoseismic evidence in the lacustrine sediments. Stacked sequences of MTDs along the steep eastern flanks of the lake appear to grade basin-ward from thick, chaotic, blocky masses to thin, parallel-bedded turbidite beds. The discrete turbidite event layers are separated by fine-grained (silt and clay) lake sedimentation. The event layers are observed throughout the lake, but the physical characteristics of the deposits vary considerably depending on proximity to primary depocenters, steep slopes, and subaqueous deltas. A total of 30–34 event deposits are observed in the post-glacial record. Radiometric dating was used to reconstruct a detailed sedimentation history over the last ~5.5 k.y., develop an age model, and estimate the recurrence (365–405 yr) for the most recent 12 event layers. Based on sedimentological characteristics, temporal overlap with other regional paleoseismic chronologies, and recurrence estimates, at least 10 of the dated event layers appear to be sourced from slope failures triggered by intense shaking during megathrust ruptures; the recurrence interval for these 10 events is 440–560 yr. Thus, Ozette Lake contains one of the longest and most robust geological records of repeated shaking along the northern Cascadia subduction zone.

Washington↗
USGS2026

Utilization of multiple geochronology techniques to constrain the age of laterization and mineralization of the world-class Mount Weld rare earth element deposit, Western Australia

Pervasive chemical weathering on stable cratons may form thick regoliths and elemental enrichment, but constraining the age of regolith formation is challenging. In this study we utilize multiple geochronological techniques on different minerals from the world-class Mount Weld rare earth element (REE) deposit, formed by lateritic weathering of a carbonatite, to constrain the age of formation and provide insight into landscape evolution. The oldest dates, ca. 100 to 38 Ma, are from Lu-Hf dating of churchite [HREE(PO 4 )·2(H 2 O)], a heavy REE phosphate mineral. Growth bands on individual minerals show a younging outwards. 40 Ar/ 39 Ar geochronology of cryptomelane [K(Mn 4+ ,Mn 2+ )₈O₁₆] yielded dates from ca. 40 to 27 Ma. Similarly, (U-Th)/He geochronology of goethite [FeO(OH)] yielded dates ranging from ca. 45 to 19 Ma. Integrating results into regional constraints, suggests 1) churchite formed by mineral saturation in a karst-like setting below the water table from ca. 100 to 40 Ma, 2) with minor uplift and erosion, cryptomelane and goethite formed at or near the water table between ca. 45 and 19 Ma, 3) after ca. 15 to 10 Ma chemical weathering within the profile had ended. Other studies document that the region experienced minimal uplift and a wet, warm climate from ca. 100 Ma to 15 Ma. These conditions and the high carbonate content of the carbonatite promote extensive chemical weathering, a deep weathering profile, and the preservation of the weathered section. This study highlights the use of multiple geochronological techniques utilizing different minerals to provide insight into how laterites form and to constrain the timing and history of the formation of this important mineral deposit.

Western Australia↗
USGS2025

Multi-scale geophysical mapping of the brine and bedrock surfaces along the Dolores River, Paradox Valley, Colorado, December 2023

Total dissolved solids derived from salt dome–sourced brine in the underlying alluvial aquifer substantially increase with distance in the reach of the Dolores River that passes through Paradox Valley in southwestern Colorado. The area has been the site of salinity control operations since the 1990s to reduce salt loading to the downstream Colorado River. Previous airborne and ground/water-based electromagnetic (EM) geophysical data have successfully mapped the top of the brine surface, albeit with relatively coarse near-surface resolution and limited spatial coverage. This present December 2023 study used ground-based high-resolution EM and passive seismic (horizontal-to-vertical spectral ratio, HVSR) tools to map in detail the depth and thickness of the brine zone in the alluvial aquifer (top of the brine down to bedrock contact) in areas immediately surrounding the Dolores River where previous airborne EM (AEM) results indicated brine within 10 m of land surface. Results indicate the deepest bedrock is generally associated with the shallowest brine and local depressions in the collapse breccia (caprock to the Paradox Formation salt) may facilitate vertical migration of brine into the alluvial aquifer. Additionally, the ground-based EM mapping corroborated general patterns in depth to brine that were observed in previous AEM results while also revealing additional detail, including suspected focused brine discharge zones to the Dolores River. A river-based survey wherein EM data, channel depth and river water–specific conductance information were collected augmented these findings and indicated specific areas on both the western and eastern sides of the river where focused brine discharge may occur. This study comprises a large-scale, ground- and water-based geophysical mapping effort, including hundreds of HVSR soundings and 100s of kilometres of EM data, which were successfully translated into digital brine and bedrock surfaces that could be incorporated into groundwater modelling efforts, future well siting or other decision-making.

Colorado↗
USGS2026

Preliminary bedrock geologic map of the Port Henry quadrangle, Essex County, New York, and Addison County, Vermont

Introduction The bedrock geology of the 7.5-minute Port Henry quadrangle consists of deformed and metamorphosed Mesoproterozoic gneisses of the Adirondack Highlands unconformably overlain by weakly deformed lower Paleozoic sedimentary rocks of the Champlain Valley. The Mesoproterozoic rocks occur on the eastern edge of the Adirondack Highlands and represent an extension of the Grenville Province of Laurentia. Mesoproterozoic paragneiss, marble, and amphibolite hosted the emplacement of an anorthosite-mangerite-charnockite-granite (AMCG) suite, now exposed mostly as orthogneiss, at approximately 1.18–1.15 Ga (giga-annum). In the Port Henry quadrangle, the AMCG metaigneous rocks (Yhg, Ygb, Yanw) intruded older, mostly metasedimentary rocks of the Grenville Complex during the middle to late Shawinigan orogeny (~1,160–1,150 Ma [mega-annum]). All rocks were subsequently metamorphosed to upper amphibolite to granulite facies conditions during the 1,080–1,050 Ma Ottawan orogeny. New mapping reveals four periods of deformation: (1) D1 produced rarely preserved isoclinal folds in the paragneiss and marble and predates AMCG magmatism. (2) Subsequent D2 deformation produced the dominant gneissic fabric preserved in the rock, recumbent folding, and deformed all the Proterozoic units in the map area. Syn- to late-D2 felsic magmatism resulted in the regionally extensive Lyon Mountain Granite Gneiss, which hosts numerous magnetite ore bodies. (3) Mylonitic extensional shear zones and core complex formation marked the beginning of D3 deformation. Protracted D3 deformation resulted in F3 upright folding, dome and basin formation, pegmatite intrusion, reactivation of the S2 foliation, partial melting, metamorphism, metasomatism, iron-ore remobilization, and intrusion of magnetite-bearing pegmatite both as layer-parallel sills and crosscutting dikes. (4) D4 created northeast- and northwest-trending local high-grade ductile shear zones and boudinage, northwest-trending regional kilometer (km)-wide ductile shear zones, and crosscutting granitic pegmatite dikes. The development of the late-stage regional shear zones (D4) was likely due to the continuation of extensional doming and uplift from upper amphibolite facies conditions at the end of the Ottawan orogeny. The majority of iron-ore deposits in the Port Henry and adjacent Witherbee quadrangles are in the hanging wall of these extensional shear zones. In the Port Henry quadrangle, the km-wide Cheney Mountain shear zone is the result of D4 deformation. Kilometer-scale lineaments readily observed in lidar data are Ediacaran mafic dikes and Phanerozoic brittle faults. The Paleozoic rocks are part of the Early Cambrian to Late Ordovician carbonate bank on the ancient margin of Laurentia. The approximately 1-km-thick Cambrian to Ordovician stratigraphy records a transition from synrift clastics to passive-margin peritidal carbonate buildups to gradually deeper-water subtidal- to shelf-carbonates during foreland basin development associated with the Taconic orogeny. The Paleozoic rocks are weakly folded and block faulted. Large areas of the Champlain Valley are covered by undifferentiated glacial deposits, some of which contain mapped landslides. The map also shows waste rock piles and tailings from historical mining operations. This study was undertaken to improve our understanding of the bedrock geology in the Adirondack Highlands, establish a modern framework for 1:24,000-scale bedrock geologic mapping in the Adirondacks, provide a context for historical iron mines in the eastern Adirondacks, and update the stratigraphy of the Champlain Valley in New York and Vermont. This Open-File Report includes a bedrock geologic map; a description of map units; a correlation of map units; and a geographic information system database that includes bedrock geologic units, faults, outcrops, and structural geologic information.

New York, Vermont↗
USGS2026

Determining Volcanic Risk in Auckland (DEVORA) Research Programme—A transdisciplinary approach to address the challenge of distributed volcanism in an urban environment

The Determining Volcanic Risk in Auckland (DEVORA) Research Programme was launched in 2008 to address the challenges associated with monogenetic volcanism in an urban setting and to enhance volcanic risk management in Tāmaki Makaurau Auckland in Aotearoa New Zealand. It is a multi-agency, increasingly transdisciplinary (defined here as research that transcends traditional disciplinary boundaries by integrating diverse types of knowledge, perspectives, and methods from academic and non-academic participants to create novel solutions to complex problems), and collaborative research program jointly led by Waipapa Taumata Rau University of Auckland and Earth Sciences New Zealand (ESNZ; formerly GNS Science), with core funding from Natural Hazards Commission Toka Tū Ake (NHC; formerly the Earthquake Commission, EQC) and Te Kaunihera o Tāmaki Makaurau Auckland Council (AC). The primary research focus of DEVORA is to investigate the geologic history, volcanic hazards, and risk posed by the basaltic intraplate Auckland Volcanic Field. Disruption from ash fall and gas from other Aotearoa New Zealand volcanoes is also considered. DEVORA’s work to explore exposure and vulnerability in Tāmaki Makaurau Auckland is also useful for assessing risks from other non-volcanic natural hazards, such as seismic and tsunami hazards. The greater Tāmaki Makaurau Auckland region has an ethnically and socio-economically diverse population of approximately 1.7 million, representing about one-third of the Aotearoa New Zealand population, and hosts critical infrastructure of national significance. The size and nature of the populace, consequential economic base, and important infrastructure within Tāmaki Makaurau Auckland mean that the effects of a volcanic eruption would be felt nationally, including through the disruption of air travel to Aotearoa New Zealand. The hazards from such an eruption could potentially affect hundreds of thousands of people, businesses, and lifelines (critical infrastructure). A considerable challenge for emergency and risk managers is the monogenetic nature of the volcanic field. It is not known where or when the next eruption will occur, how much warning we may get before an eruption, nor how an eruption and its effects might unfold. In this contribution, we highlight the concept and collaborative intent of the DEVORA Programme and show how it has evolved over the 16 years since its inception. We describe how DEVORA has unified more than 100 researchers (including more than 50 graduate students) and numerous stakeholders to address key issues facing Tāmaki Makaurau Auckland and describe how research findings are being implemented into policy and communicated to stakeholder agencies and the public. We also illustrate the broader influence of the DEVORA Programme and provide some learnings that might benefit others embarking on similar integrated projects, especially those focused on distributed volcanism in and near populated areas.

Auckland↗
USGS2026

Vegetation cover and composition in environments surrounding uranium mines in the Grand Canyon ecosystem, Northern Arizona

Mining uranium from breccia-pipe deposits in the greater Grand Canyon region has occurred since the mid-1900s. However, possible ecosystem contamination with harmful levels of radionuclides may have occurred due to mining activities in the 21st century. In response, a 20-year Federal moratorium on new mining claims in the Grand Canyon watershed was initiated in 2012, to allow time to evaluate the potential effects of uranium exploration and mining on human health, wildlife, and water resources. This moratorium, nor the 2023 designation of the “Baaj Nwaavjo I’tah Kukveni–Ancestral Footprints of the Grand Canyon National Monument,” precludes operation or development of mining claims predating 2012. Vegetation is a core ecosystem component that may be affected by uranium mining (for instance, through uptake and storage of radionuclides from the air or soil) or may act as a vector of exposure to wildlife, livestock, and humans (for instance, via their consumption of contaminated plant tissues). To provide baseline information about the plant communities associated with uranium mines in the Grand Canyon region, the U.S. Geological Survey surveyed an approximately 200-meter-wide buffer surrounding four breccia-pipe deposits, each in a unique stage of mine development, and at one reference area (a livestock water tank) that underwent ground disturbance but contains no mineral deposits. We sectioned the buffer zones into 0.65–4.52 hectare plots, within which we (1) inventoried all plant species, (2) measured percent cover of plant species, plant functional groups, and ground surface types (dark cyanobacteria, lichen, moss, bedrock, rock, embedded litter, duff, plant bases, and bare soil) using line-point intercept, and (3) measured length and frequency of gaps between perennial plant canopies using canopy gap intercept. We found that plant composition at the mines and the reference area differed from one another but were all characteristic of expected regional vegetation patterns. We provide this data summary as potential baseline information for future research and management efforts.

Arizona↗
USGS2026

Small earthquake moment magnitude and implications for frequency–magnitude scaling of injection induced earthquakes of the Raton Basin

Accurate estimation of earthquake source parameters—such as moment magnitudes, corner frequencies, and stress drops—is essential for improving seismic hazard assessments and understanding earthquake physics. In this study, moment magnitudes ( M W ) are calculated for 31,581 earthquakes associated with wastewater injection in the Raton Basin (located along the border between northern New Mexico and southern Colorado) between 2016 and 2024 using radiative transfer theory to fit coda decay envelopes. Our results show that it is feasible to estimate moment magnitudes down to M W ~1 with coda envelopes from a small local monitoring network. Significant differences were found between M W and local magnitudes ( M L ) for small earthquakes ( M < 3.0). A linear relationship was optimized to convert M L to M W : M W = 0.7 M L + 0.96 and M W = 0.73 M L + 0.99 (for the events reported by the U.S. Geological Survey), which can be applied in future studies of Raton Basin seismicity. We find that b -values calculated employing different methods and using M L are approximately 1.0, while those using M W range from 1.2 to 1.4. A larger estimate of the b -value could influence interpretations of the statistical behavior of earthquakes associated with injection and consequently seismic hazard assessments based on a magnitude–frequency distribution. The potential differences between local versus moment magnitude-based earthquake statistics should be considered in other seismically active regions.

Colorado, New Mexico↗
USGS2000

Monthly variability and possible sources of nitrate in ground water beneath mixed agricultural land use, Suwannee and Lafayette Counties, Florida

In an area of mixed agricultural land use in Suwannee and Lafayette Counties of northern Florida, water samples were collected monthly from 14 wells tapping the Upper Floridan aquifer during July 1998 through June 1999 to assess hydrologic and land-use factors affecting the variability in nitrate concentrations in ground water. Unusually high amounts of rainfall in September and October 1998 (43.5 centimeters total for both months) resulted in an increase in water levels in all wells in October 1998. This was followed by unusually low amounts of rainfall during November 1998 through May 1999, when rainfall was 40.7 centimeters below 30-year mean monthly values. The presence of karst features (sinkholes, springs, solution conduits) and the highly permeable sands that overlie the Upper Floridan aquifer provide for rapid movement of water containing elevated nitrate concentrations to the aquifer. Nitrate was the dominant form of nitrogen in ground water collected at all sites and nitrate concentrations ranged from less than 0.02 to 22 milligrams per liter (mg/L), as nitrogen. Water samples from most wells showed substantial monthly or seasonal fluctuations in nitrate concentrations. Generally, water samples from wells with nitrate concentrations higher than 10 mg/L showed the greatest amount of monthly fluctuation. For example, water samples from six of eight wells had monthly nitrate concentrations that varied by at least 5 mg/L during the study period. Water from most wells with lower nitrate concentrations (less than 6 mg/L) also showed large monthly fluctuations. For instance, nitrate concentrations in water from four sites showed monthly variations of more than 50 percent. Large fluctuations in nitrate concentrations likely result from seasonal agricultural practices (fertilizer application and animal waste spreading) at a particular site. For example, an increase in nitrate concentrations observed in water samples from seven sites in February or March 1999 most likely results from application of synthetic fertilizers during the late winter months. Lower nitrate concentrations were detected in water samples from five of eight wells sampled during high-flow conditions for the Suwannee River in March 1998 compared to low-flow conditions in November 1998. Evidence for reduction of nitrate due to denitrification reactions was observed at one site (AC-1), as indicated by elevated concentrations of nitrogen gas and a corresponding increase in nitrogen isotope (d15N-NO3) values with a decrease in nitrate concentrations. Denitrification is unlikely at other sites based on the presence of dissolved oxygen concentrations greater than 2 mg/L in ground water and no observed trend between nitrate concentrations and values d15N-NO3 values. Nitrate was the dominant nitrogen species in most monthly rainfall samples; however, ammonium concentrations were similar or greater than nitrate during November and December 1998. During February through May 1999, both nitrate and ammonium concentrations were substantially higher in monthly rainfall samples collected at the study area compared to mean monthly concentrations at the Bradford Forest site located east of the study area, which is part of the National Atmospheric Deposition Program/National Trends Network. Also, higher nitrogen deposition rates in the study area compared to those at Bradford Forest could indicate that substantial amounts of ammonia are volatilized from fertilizers and animal wastes, released to the atmosphere, and incorporated as nitrate and ammonium in rainfall deposited in the middle Suwannee River Basin. Ground-water samples from most sites had d15N-NO3 values that indicated a mixture of inorganic and organic sources of nitrogen, which corresponded to multiple land uses where both synthetic fertilizers and manure are used on fields near these sites. Distinct d15N-NO3 signatures, however, were observed at some sites. For example, water samples from areas of row-crop farming as the dominant land use had d15N-NO3 values less than 4 per mil, indicating an inorganic nitrogen source such as synthetic fertilizer. In contrast, d15N-NO3 values greater than 9 per mil were found in water samples from three sites where manure from dairy and or poultry operations is spread on fields throughout the year.

Florida↗