Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Geoderma Regional”

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.

586 records · Page 6Linked to original sources
USGS2024

The spatially adaptable filter for error reduction (SAFER) process: Remote sensing-based LANDFIRE disturbance mapping updates

LANDFIRE (LF) has been producing periodic spatially explicit vegetation change maps (i.e., LF disturbance products) across the entire United States since 1999 at a 30 m spatial resolution. These disturbance products include data products produced by various fire programs, field-mapped vegetation and fuel treatment activity (i.e., events) submissions from various agencies, and disturbances detected by the U.S. Geological Survey Earth Resources Observation and Science (EROS)-based Remote Sensing of Landscape Change (RSLC) process. The RSLC process applies a bi-temporal change detection algorithm to Landsat satellite-based seasonal composites to generate the interim disturbances that are subsequently reviewed by analysts to reduce omission and commission errors before ingestion them into LF’s disturbance products. The latency of the disturbance product is contingent on timely data availability and analyst review. This work describes the development and integration of the Spatially Adaptable Filter for Error Reduction (SAFER) process and other error and latency reduction improvements to the RSLC process. SAFER is a random forest-based supervised classifier and uses predictor variables that are derived from multiple years of pre- and post-disturbance Landsat band observations. Predictor variables include reflectance, indices, and spatial contextual information. Spatial contextual information that is unique to each contiguous disturbance region is parameterized as Z scores using differential observations of the disturbed regions with its undisturbed neighbors. The SAFER process was prototyped for inclusion in the RSLC process over five regions within the conterminous United States (CONUS) and regional model performance, evaluated using 2016 data. Results show that the inclusion of the SAFER process increased the accuracies of the interim disturbance detections and thus has potential to reduce the time needed for analyst review. LF does not track the time taken by each analyst for each tile, and hence, the relative effort saved was parameterized as the percentage of 30 m pixels that are correctly classified in the SAFER outputs to the total number of pixels that are incorrectly classified in the interim disturbance and are presented. The SAFER prototype outputs showed that the relative analysts’ effort saved could be over 95%. The regional model performance evaluation showed that SAFER’s performance depended on the nature of disturbances and availability of cloud-free images relative to the time of disturbances. The accuracy estimates for CONUS were inferred by comparing the 2017 SAFER outputs to the 2017 analyst-reviewed data. As expected, the SAFER outputs had higher accuracies compared to the interim disturbances, and CONUS-wide relative effort saved was over 92%. The regional variation in the accuracies and effort saved are discussed in relation to the vegetation and disturbance type in each region. SAFER is now operationally integrated into the RSLC process, and LANDFIRE is well poised for annual updates, contingent on the availability of data.

Fire↗
USGS2025

SURF: An automated method for building nonplanar 3D fault models from earthquake hypocenters

Accurately characterizing 3D fault geometry is vital for improving our understanding of earthquake behavior and informing the development of seismic hazard models. Despite their importance, subsurface fault structures tend to be poorly constrained because of limitations in observational data. Improvements to the seismic networks and earthquake detection algorithms have increased the precision and volume of earthquake catalogs, which help illuminate detailed subsurface fault structure and provide the most direct information available about fault geometries at depth. We present a Python package to automate generating 3D fault geometries directly from hypocentral seismicity patterns. This method begins with clustering events based on their spatial density, identifying coherent patterns. Nearby clusters are then merged based on the similarity of their orientations. We fit nonplanar surfaces using support vector regression to balance surface accuracy with minimal deviations from planarity. The fault models are output as quadrilateral meshes at user‐defined resolution. In the process of generating the 3D fault surfaces, we compute the spatial density of seismicity around the surface and the planarity as quantitative metrics of the model outputs. As a proof of concept, we apply this approach to the San Andreas–Calaveras fault junction region and the 2019 Ridgecrest earthquake sequence, both in California, which contain complex subparallel faults well defined at the Earth’s surface and abundant microseismicity. These case studies demonstrate the method’s ability to model complex fault structures, including long continuous fault surfaces, crossing faults, variably dipping segments, and subparallel faults. We test the method on both standard network catalogs and double‐difference relocated catalogs. We find that our seismicity‐based fault model results align with published 3D models that incorporate additional constraints and interpretations ( Plesch et al. , 2020 ; Aagaard and Hirakawa, 2021 ). This workflow provides a low‐user‐input solution for estimating fault geometries at depth from earthquake catalogs.

California↗
USGS2026

Effects of human activity on Khumbu Glacier: Towards a sustainable Everest Base Camp

The Everest Base Camp (EBC), a critical staging area for mountaineers attempting to summit Mount Everest, has seen a significant increase in human activity over the past few decades. The increasing number of climbers, trekkers, and support staff, particularly at the campsite area, has intensified environmental impacts on the Khumbu Glacier. This study investigates the effects of human activity at EBC, focusing on three contributors to glacier melt: global climate change, local fossil fuel consumption, and human urinary discharge. Through field data collected during the spring 2023 climbing season, this study examines the energy released through the use of liquefied petroleum gas (LPG), kerosene, and petrol for cooking, heating, and electricity, as well as the heat generated by human urinary discharge. The total heat energy released from these activities at EBC in the spring of 2023 was 849,174 ± 179,774 MJ, which would be sufficient to melt 2492 ± 528 tons of glacier ice and snow. Trends in land surface temperature at EBC and adjacent environments during the 1991–2023 period were calculated from Landsat 5, 7, and 9 satellite data. The 32-year Landsat record reveals that EBC surface temperature increased by 0.28 °C year −1 , which was roughly twice the rate of warming at the surface of Khumbu Glacier adjacent to EBC and about 15% higher than the temperature increase of the debris-covered terrain immediately to the north of EBC. The findings suggest that addressing these anthropogenic influences could help to preserve the Khumbu Glacier and support the sustainability of mountaineering in the region. This paper also takes a transformative approach and explores stable and safer locations in case there is a necessity to relocate the current EBC. Two sites have been identified southwest of the existing EBC, which, unlike the current supraglacial site, are situated on stable ground.

Everest Base Camp, Khumbu Glacier↗
USGS2026

As above, so below? A framework for integrating long-term water quantity trends reveals divergent patterns in groundwater and low streamflow across the United States

Climate, land-use, and disturbance drive long-term global trends in groundwater levels and streamflow. At large scales, these trends are typically considered separately, despite the well-established concept that groundwater and surface water comprise a single resource. Joint trend assessment at national scales is challenging because it requires pairing and aggregating data from spatially disparate streamflow and groundwater monitoring sites for which no established framework exists. Here, we evaluate alternative approaches for integrating groundwater and streamflow data to enable joint trend analysis—a critical step toward understanding how water-budget components respond concurrently and interactively to environmental drivers. Mann–Kendall trends were computed for individual groundwater (annual mean depth) and streamflow (annual low of 7 d averages) sites across the U.S over 21- (2000–2020), 31- (1990–2020), and 41-year (1980–2020) periods. Regional Kendall trends were calculated using five spatially contiguous and noncontiguous regional classifications for aggregation based on subsurface (e.g. aquifer, geology) and surface (e.g. watershed, landscape) characteristics. Site-level results revealed contrasting trends, with tendencies toward increasing low flows (wetting) and increasing groundwater depths (drying). Agreement between streamflow and groundwater trends increased with regional aggregation and longer timeframes, though persistent skew toward streamflow wetting and groundwater drying remained. Results varied by region and trend period, with notable consistencies: unified drying in the West/Southwest and wetting in the Upper Midwest. Directional mismatches in long-term trends were prominent in the High Plains and Mississippi Alluvial Plain, whereas near-term mismatches were most evident in the Northwest. Aggregation by hydrologic landscape regions (HLR) yielded the greatest agreement between groundwater and streamflow trends. These findings indicate that coupled responses may represent combined influences of climate, relief, and geology, as captured by HLR, more strongly than geography or geology alone. Integrated water availability assessments may benefit from a multi-characteristic classification framework to treat groundwater and surface water as a unified resource.

Environmental Research: Water↗
USGS2024

Comparing subduction ground-motion models to observations for Cascadia

We evaluate Cascadia subduction ground-motion models (GMMs), considered for the 2023 US National Seismic Hazard Model (NSHM) update, by comparing observations to model predictions. The observations comprise regional recordings from intraslab earthquakes, including contributions from 2021 and 2022 events in southern Cascadia and global records from interface earthquakes. Since the 2018 NSHM update, new GMMs for Cascadia have been published by the Next Generation Attenuation (NGA)-Subduction Project that require independent evaluation. In the regional intraslab comparisons, we highlight a characteristic frequency dependence for Cascadia data, with short periods having lower ground motions and longer periods being comparable to other subduction zones. We evaluate differences in northern and southern Cascadia and find that the NGA-Subduction GMMs developed using southern Cascadia data perform better in this region than the model that did not consider these data. We compare ground-motion variability in Cascadia with the NGA-Subduction model predictions and find differences at short periods ( T = 0.1 s) due to the use of global versus regional data in the development of these models. Moreover, the within-event component of aleatory variability from the GMMs overpredicts the standard deviation of Cascadia recordings at very short periods ( T < 0.05 s). Using global interface earthquakes as a proxy to evaluate the Cascadia GMMs, we find long-period overprediction from a simulation-based GMM and some of the empirical GMMs. When comparing recent observations, we find a similar misfit to GMMs and the 2010 and 2022 Ferndale earthquakes. Finally, we observe different basin amplification factors arising in different subsets of the data, which indicate that differences in basin factors between empirical GMMs could arise from the data selection choices by the developers. As part of evaluating the regional basin terms, we apply basin amplification factors from the magnitude 9 Cascadia earthquake simulations to the empirical GMMs for interface earthquakes. The comparisons presented in this study indicate that the NGA-Subduction GMMs for Cascadia perform well relative to observations and older subduction GMMs.

British Columbia, California, Oregon, Washington↗
USGS2025

Wetter winters, drier summers: Quantifying the change in hydrological response around the Puget Sound area using the wflow_sbm hydrological model and CMIP6 projections

Climate change is expected to impact hydrological regimes worldwide, including the Pacific Northwest of the United States. This study investigates how climate change will affect river discharge in the Puget Sound region of the State of Washington, with a focus on King and Pierce Counties. We simulated river discharge under historical and future conditions using the physically based, spatially distributed wflow_sbm hydrological model, which was calibrated and validated against U.S. Geological Survey discharge records. Future forcing was based on an ensemble of six high-resolution CMIP6 climate models, which were bias corrected using empirical quantile mapping. The results indicate a decrease in summer discharges (5–10%) and an increase in winter discharges (5–10%) across the study region. The high discharges (90th percentile) are projected to increase in winter, and the low discharges are projected to decrease in summer, due to shifts in precipitation regimes, snowpack hydrology, and evapotranspiration. However, variability between individual CMIP6 models often exceeds the magnitude of ensemble mean changes, underscoring substantial uncertainty in climate projections and the importance of including multiple climate models in climate change analysis. Furthermore, model consensus increased with elevation, which could be the result of the higher elevation areas being driven by less diverse hydrological processes. These findings highlight potential challenges for regional water management, ecosystem health, and flood risk mitigation in the Puget Sound region under future climate conditions.

EarthArXiv↗
USGS2024

Effective site coefficients for the 2024 International Building Code (IBC)

The U.S. National Seismic Hazard Models (NSHMs), developed by the U.S. Geological Survey (USGS), have long been the scientific foundation of seismic design guidelines and have been used to compute design ground motions for construction of new buildings and retrofit of existing buildings in the United States and its territories. The 2018 NSHM is adopted by the 2024 International Building Code (IBC). Prior to the 2018 NSHM update, hazard calculations were required only at one reference site condition defined by V S30 =760 m/s (where V S30 is the time-averaged shear wave velocity from the surface to a depth of 30 m) and three periods (peak ground acceleration, PGA, and pseudo spectral accelerations at periods of 0.2 and 1 s, S s and S 1 ). Site coefficients, F PGA , F a , and F v , were then defined by the Building Seismic Safety Council (BSSC) Provisions Update Committee (PUC) in the site-specific procedures of National Earthquake Hazard Reduction Program (NEHRP) Recommended Seismic Provisions to calculate ground motions for other site classes with different V S30 values at the given periods. The design ground motions at other periods were then estimated using a generic spectral shape that was also defined by the BSSC PUC in NEHRP provisions. In recent years, the engineering community has realized there were deficiencies with the F a and F v site coefficients and the generic spectral shape. To avoid potentially dangerous underestimations of design ground motions for long period structures on soft site conditions in high seismicity regions, the BSSC PUC recommended the use of multi-period response spectra (MPRS) in 2017. As a result, the USGS produced multi-period and multi-V S30 response spectral values in the 2018 NSHM for calculations of design ground motions and the site coefficients F a and F v were eliminated from the 2020 NEHRP Provisions. As these site coefficients were widely used inside and outside of the United States, in this study we back-calculate the “effective” site coefficients F a,eff , and F v,eff by comparing MPRS for various site classes with the MPRS for the reference site condition, and discuss the changes that are observed in the 2024 IBC compared to its previous version in 2021. The effective site coefficients are presented for test site locations and their dependence on various factors including period, ground motion intensity, and regional models are discussed. Ratio maps between the new effective site coefficients and the old ones are then presented for soft site classes and for short and long periods. For soft site classes at short periods, the new effective site coefficients are lower than the old site coefficients for high seismicity regions and higher for low seismicity regions. As it was expected, for soft site classes at long periods and high seismicity regions, the new effective site coefficients are much larger than the old site coefficients without imposing the 50% increase as a penalty that was suggested in the 2021 IBC, whereas they could be much smaller if the 50% increase would have been imposed particularly around New Madrid and Charleston high seismicity regions. For low seismicity regions, the long period effective site coefficients can be smaller or larger by 20% compared to the 2021 IBC coefficients.

Conference Paper↗
USGS2026

An improved empirical model for predicting postfire debris-flow volume in the western United States

Reliable estimates of debris-flow volume can be used to help predict the magnitude of debris-flow hazards following wildfire in the western United States. In this study, we compiled and used a database of 227 postfire debris-flow volumes that were collected across the western United States to develop a multiple linear regression model for predicting postfire debris-flow volume. We explored 36 predictor variables related to rainfall, terrain, and fire characteristics, and selected the model with the combination of variables that yielded the most accurate predictions of debris-flow volume. We evaluated model performance against the entire volume database, as well as against four subsets of volume data from southern California, the Intermountain West, the Southwest, and regions with limited volume data, such as northern California and Washington. We also compared model performance against 3 existing postfire debris-flow volume models that were developed for use in southern California, the Intermountain West, and the Southwest. We demonstrate that the new volume model performs as well as the regional models in the regions for which they were developed and outperforms existing models when applied to volumes from data-limited regions in the western United States. These results indicate that the debris-flow volume model introduced in this study can be used to improve postfire hazard assessments across the western United States, especially outside of southern California.

Arizona, California, Colorado, New Mexico, Utah, W↗
USGS2024

Global geologic map of Europa

Discovered by Galileo Galilei more than 400 years ago and imaged in detail by the Voyager 2 Galileo spacecraft, Jupiter’s icy moon Europa has been a source of intrigue. A range of science investigations indicate that it contains the key ingredients for habitability, notably energy, chemistry, and liquid water. Europa’s surface is geologically complex and, based on the dearth of impact craters, interpreted to be as young as ~60 Ma. The array of geologic features that characterize the surface include extensive ridged plains, regions of broad disruption termed chaotic terrain, long, quasi-linear ridges that span thousands of kilometers, and bands as much as 60 kilometers wide and that extend 100s of kilometers. These features, along with other geophysical measurements, indicate the presence of a global briny liquid water ocean beneath the ice shell. It was not until the arrival of the Galileo spacecraft in 1995 that the true nature and level of complexity of the surface was revealed. Although image data returned by Galileo provided insight into the structure of a variety of regions, the entire satellite has yet to be observed at a regional scale (less than 250 meters per pixel) and the detailed geologic nature of much of its surface remains a mystery. Establishing the global context of the distribution and timing of Europan geologic units forms a basis to understand regional and local scale processes, serves as a tool for the planning of future missions, and most of all is essential to gaining insight into the potential habitability of this icy world.

Scientific Investigations Map↗
USGS2024

Separating signals in elevation data improves supervised machine learning predictions for hydrothermal favorability

A recent study identified topography (land surface elevation above sea level) as an important input dataset (feature) for predicting the location of hydrothermal systems in the Great Basin in Nevada. Yet, topography is generally a result of more than one geological process and may consequently contain multiple distinct signals. For example, the geologic evolution of the Great Basin has produced both crustal thickening (i.e., regional-scale trends in elevation) and thinning via Basin and Range extensional faulting (i.e., valley-scale topographic relief). We postulate that these geologic processes may affect the occurrence of hydrothermal systems differently. Therefore, we separate the regional trend from the valley-scale signal in the Great Basin, and then use them separately to evaluate the importance of each as predictors for hydrothermal favorability. Our prior work applying supervised machine learning (ML) using the data from the Nevada Machine Learning Project demonstrated that employing a training strategy that randomly selects negative training sites produces better performing models for predicting hydrothermal favorability than a training strategy that uses expert-selected negatives. The models created using both training strategies exhibited a west-east geographic trend in the predictions for the favorability of hydrothermal resources. These models generally predicted higher favorability in western Nevada and lower favorability in eastern Nevada. This west-east trend in predicted favorability correlates with elevation across the Great Basin, which trends higher from west to east. By separating the original elevation feature into distinct features for elevation trend (i.e., regional-scale topography) and detrended elevation (i.e., valley-scale or local relative topography), we find that models using the separated topographic signals consistently outperform competing models that use the original elevation feature. Although western Nevada still exhibits higher favorability than eastern Nevada, using separated signals for regional elevation and local structure reduces the west-east prediction trend in the region and emphasizes structures associated with hydrothermal upflow. This work emphasizes how carefully engineering features to represent geological conditions relevant to hydrothermal systems allows ML algorithms to detect important patterns for predicting hydrothermal resource favorability and leads to better model performance.

Conference Paper↗
USGS2026

Moment magnitude for small earthquakes in the Delaware basin of west Texas and southeast New Mexico, USA

The Delaware Basin region of west Texas and southeast New Mexico has become one of the most prolific regions of seismic activity in the continental United States due to widespread hydraulic fracturing and wastewater disposal injection. In response to the increased number of earthquakes in this region, rapid and accurate characterization of earthquake sources is necessary to understand the evolution of seismic activity and level of seismic hazard associated with these earthquakes. This study re-evaluates earthquake magnitudes, estimating moment magnitude (MW) for small earthquakes in the Delaware Basin using 1) moment-rate spectra derived from S-wave coda envelopes, and 2) a relative magnitude method that relies exclusively on the ratio of waveform amplitudes between highly correlated waveform pairs. The coda-envelope method produces accurate M W estimates for small earthquakes ( M 1.5 – 3) that are consistent with independent, waveform modeled moment magnitudes for events with M W > 3 . Using the relative amplitudes method to extend these M W magnitudes to many other events, we successfully provide relative moment magnitude ( M W,rel ) values for 81% of the Texas Seismological Network catalog in the Delaware Basin region, and 45% of the USGS Induced Seismicity Project’s catalog of events in southeast New Mexico. The adoption and integration of the calibrated M W,rel method with current magnitude estimation methods offers valuable insights into the relationships between local and moment magnitude and will contribute to improved characterization of widespread induced seismicity.

New Mexico, Texas↗
USGS2024

The dominance and growth of shallow groundwater resources in continuous permafrost environments

Water is a limited resource in Arctic watersheds with continuous permafrost because freezing conditions in winter and the impermeability of permafrost limit storage and connectivity between surface water and deep groundwater. However, groundwater can still be an important source of surface water in such settings, feeding springs and large aufeis fields that are abundant in cold regions and generating runoff when precipitation is rare. Whether groundwater is sourced from suprapermafrost taliks or deeper regional aquifers will impact water availability as the Arctic continues to warm and thaw. Previous research is ambiguous about the role of deep groundwater, leading to uncertainty regarding Arctic water availability and changing water resources. We analyzed chemistry and residence times of spring, stream, and river waters in the continuous permafrost zone of Alaska, spanning the mountains to the coastal plain. Water chemistry and age tracers show that surface waters are predominately sourced from recent precipitation and have short (<50 y) subsurface residence times. Remote sensing indicates trends in the areal extent of aufeis over the last 37 y, and correlations between aufeis extent and previous year summer temperature. Together, these data indicate that surface waters in continuous permafrost regions may be impacted by short flow paths and shallow suprapermafrost aquifers that are highly sensitive to climatic and hydrologic change over annual timescales. Despite the lack of connection to regional aquifers, continued warming and permafrost thaw may promote deepening of the shallow subsurface aquifers and creation of shallow taliks, providing some resilience to Arctic freshwater ecosystems.

Proceedings of the National Academy of Sciences↗
USGS2024

Simulation of groundwater flow in the Long Island, New York regional aquifer system for pumping and recharge conditions from 1900 to 2019

The U.S. Geological Survey has developed a transient, groundwater-flow model that simulates hydrologic conditions in the Long Island aquifer system as part of an ongoing (since 2016) multiyear, cooperative investigation with the New York State Department of Environmental Conservation. The goals of this investigation are to assist stakeholders and resource managers to evaluate the response of the hydrologic system to changes in future hydraulic stresses. Responses in the hydrologic system include changes in water levels in the hydrogeologic units; discharge to streams, coastal waters, and subsurface infrastructure; and the extent of saline groundwater in the aquifers. Hydraulic stresses include future water-supply management and changes in land use and infrastructure. The numerical model synthesizes a diverse set of physiographic, geologic, climatic, land-use, and historical population, water use, and infrastructure data to physically represent the Long Island aquifer system from land surface to bedrock and to simulate annual hydrologic conditions between 1900 and 2019. A three-dimensional hydrogeologic framework was developed from existing and recently collected borehole geologic and geophysical data collected as part of a companion drilling program. Water-transmitting properties of the principal aquifer sediments were defined in three dimensions from new and existing lithologic logs. The distribution of recharge from precipitation was estimated from landscape characteristics and climate data. Anthropogenic recharge from wastewater, leaky infrastructure, and storm runoff were estimated from population, infrastructure, and pumping data. Water-use data, including well locations, depths, and pumping rates, were obtained from historical sources and records and used to estimate pumping stresses continuously in time and space, at an annual average time scale. The data were incorporated into a three-dimensional numerical model using the U.S. Geological Survey finite difference modeling code MODFLOW 6; the model encompassed all of Long Island and surrounding surface waters and simulated historical hydrologic conditions from 1900 to 2019. The calibration process involved trial and error adjustments using prior knowledge to improve general fit to observations followed by an inverse calibration to update and optimize input parameters, using an iterative ensemble smoother algorithm implemented in PEST++ version 5.0. This resulted in a model that generally was in good agreement with observed, dynamically varying hydrologic conditions from 1900 to 2019. The calibrated model was used to develop two base-case models for scenario testing of future, hypothetical conditions where one represented average-annual conditions, and one represented average-seasonal conditions from 2010 to 2019. The model representing average-annual conditions was modified further to represent an alternate sea-level position of 6 feet above the North American Vertical Datum of 1988, and the model representing average-seasonal conditions was modified to represent the average seasonal effects of a 5-year drought imposed upon current hydrologic conditions. Recharge is the sole source of water to the aquifer system; groundwater discharges to coastal water and streams and is withdrawn by pumped wells. Model-estimated annual recharge ranged from about 11 inches in 1965 to 41 inches in 1983. On average, from 2010 to 2019, about 23 percent of water was pumped from wells, and about 47 and 27 percent discharged to coastal waters and streams, respectively; the remaining 4 percent was water that moved into storage in the aquifer matrix. Water levels on Long Island vary naturally during time in response to changes in recharge; the amount of variation is largest in the interior of the island, in areas with highest water table altitudes near groundwater divides and lowest near streams and the coastal waters. The total range of water table altitudes on Long Island between 1900 and 2019 ranged from near 0 to more than 70 feet in western parts of Long Island. The largest range in altitudes is in New York City and is associated with areas of large historical withdrawals between the 1920s and the late 1980s. Water table altitudes generally varied by less than 10 feet in eastern Suffolk County, where the aquifer is under more natural conditions. Saltwater intrusion is of great concern on Long Island, particularly in western Long Island where both the unconfined and confined parts of the aquifer system have been intruded in response to large-scale groundwater withdrawals; however, the volume of freshwater in the islandwide aquifer system only has changed by about 5 percent between 1900 and 2019. The decadal change in the freshwater volume was largest during the early and mid-20th century, corresponding to the largest historical pumping, but that volume change did not exceed 1 percent. The negligible change in freshwater volume suggests that saltwater intrusion as of 2019 was limited at an islandwide scale but continues to occur in local areas of Queens and Nassau Counties, adversely affecting current water supplies and limiting future water supplies for affected communities. The regional groundwater model developed for this investigation is a tool that can be used to help determine the viability of current and future water supplies at a regional scale and can be used to support development of additional models at finer scale to support more focused assessments of groundwater sustainability.

New York↗
USGS2025

Solution-collapse breccia pipe uranium deposits of the southern Colorado Plateau, northwestern Arizona, USA

Some of the highest-grade uranium deposits in the United States occur in breccia pipes that formed by solution and collapse of sedimentary strata, which occur in the southern portion of the Colorado Plateau in northwestern Arizona. The host breccia pipes are up to 1200 m in vertical extent, average about 90 m in diameter, and can cross-cut strata from their base in the Mississippian Redwall Limestone to as stratigraphically high on some plateaus as the Triassic Chinle Formation. These uranium-base metal deposits are up to 600 m thick and formed within the breccia pipes where they transect the Permian Coconino Sandstone, Hermit Formation, and the Esplanade Sandstone. Of the hundreds of breccia pipes identified across this region, only a small percentage are known to contain mineralization. The main uranium ore mineral is uraninite that is intergrown with at least 20 base-metal sulfide minerals, which contribute Fe, Cu, Co, As, Pb, Zn, Ni, and Ag to the deposits. This study considered regional stratigraphy, sulfur isotope systematics, mineralogy, in situ dating, and compilation and analysis of previous work on the deposits. A comprehensive deposit model has not been published for these deposits. This analysis identified new additions to update the deposit model for these unusual, possibly unique deposits. Proposed modifications to the model include: (1) the source, mechanisms, timing of the base-metal sulfide mineral assemblages, and (2) the source, mechanism, and timing of the uranium mineralization. Sulfide and uranium deposition are shown to be separate mineralization events. The study proposes the possible role of gypsum as a source of sulfur for the sulfide minerals in the deposits. Groundwaters carrying uranium encountered the preexisting sulfides in breccia pipes, reducing the uranyl ions, and precipitating U oxide (as uraninite). Analysis of the regional stratigraphy recognized that numerous beds of gypsum are in the strata that lie only tens of meters above the breccia pipe deposits. In the breccia pipe region, if these stratigraphic units (Toroweap and Kaibab Formations) do not contain gypsum layers then the underlying pipes are not mineralized; where these Permian gypsum layers do occur, breccia pipes can host mineralization. This new understanding should be useful in identifying the prospective region for mineralized pipes.

Arizona↗
USGS2026

Unprecedented burning in tropical peatlands during the 20th century compared to the previous two millennia

Tropical peatland wildfire incidence has risen in recent decades, driven by drainage for land use and intensified by severe droughts with global climate change. These disturbances have altered vegetation structure, disrupted ecosystem functioning, and increased carbon emissions, particularly in Southeast Asia. However, the long-term history and characteristics of wildfires in tropical peatlands remain largely unknown. Here, we compiled fifty-eight macro-charcoal records from peatlands across the tropics, ranging from lowland forested to montane peatlands, to assess millennia-scale changes and controlling factors of tropical peatland burning. We divided the datasets into four main sub-regions: Neotropical, Afrotropical, Indomalayan and Australasian ecoregions to explore regional variability. Tropical peatlands had high burning levels between 0 and 850 ce , followed by a relatively low and stable period until a marked increase during the 20th century. The general trend in tropical peatland burning follows changes in global temperature, and climate variables that control the length and severity of drought events have a notable influence on peat burning before 1900 ce . During the 20th century, regional differences were observed, with declining fire trends in the Neotropical and Afrotropical regions and increasing fire trends in the Indomalayan and Australasian regions. This difference is likely attributable to human activities, and such intervention is also evident in palm swamps and hardwood swamps under similar wet, weakly seasonal climates. With the increase in anthropogenic pressures on peatlands and greater climate variability, future wildfires in peatlands are likely to become more frequent and widespread across all tropical ecoregions. Conservation and sustainable land-use practices could be used to mitigate and control peatland burning and protect these carbon-rich sinks.

Global Change Biology↗
USGS2026

Integrating detrital magnetite geochemistry and (U-Th)/He chronometry as a sediment provenance tool in geologic and metallogenic terranes

Magnetite is ubiquitous in porphyry Cu systems and in sediment sourcing both barren and mineralized regions, with potential as an indicator mineral in concealed and coarsely-mapped terranes. We develop and test a workflow for integrated geochemistry and (U–Th)/He (He) dating for inferring detrital magnetite (DMt) provenance in these settings. The ca. 70 Ma Taurus porphyry Cu–Mo(–Au) district in eastern interior Alaska serves as a test case. DMt from streams draining porphyry-related mineralization was characterized by geochemistry and mineral inclusion and microstructure observations, complemented by similar data for potential porphyry and host rock sources. Principal component analysis and clustering of DMt geochemical data resolve multiple populations in our samples geochemically and texturally compatible with derivation from metamorphic, porphyry-related hydrothermal, and igneous sources. Hydrothermal magnetite comprises ∼16–50% of DMt nearest porphyry mineralization but diminishes to ∼4% ∼15 km downstream. Subsampled grains within populations yield ∼160–110 Ma, ∼70 Ma, ∼55 Ma, and ∼20 Ma magnetite He date modes. Combined with provenance, He dates capture Early Cretaceous regional exhumation of metamorphic host rock and Late Cretaceous porphyry Cu mineralization. DMt grains showing partial hematite replacement yield ca. 55–20 Ma dates regardless of source, overlapping regional warm/wet climatic intervals. We interpret Cenozoic dates to reflect exhumation to near-surface oxidizing conditions and(or) supergene weathering. Magnetite is thus a promising target phase for (1) tracking the spatiotemporal distribution of porphyry systems, and (2) linking the formation and exhumation of these systems to a regional geologic history, both in Alaska and globally.

Alaska, Yukon↗
USGS2026

High-resolution transboundary vegetation community maps of the Sonoran and Mojave Desert ecoregion to support critical landscape conservation planning and habitat management needs

We produced a 30-m resolution binational land cover map of Bird Conservation Region 33 (BCR 33) for the U.S. North American Bird Conservation Initiative. The region covers large portions of the Sonoran and Mojave Deserts. The map can support the U.S. Fish and Wildlife Service (FWS) Migratory Bird Program’s recovery planning efforts and constitutes the first known binational land cover dataset spanning sections of the United States–Mexico border and using a consistent classification system for both countries. The mapped region includes 152 distinct land cover classes, covering a total area of 38,421,453 ha (148,345 mi 2 ), of which 13,148,345 ha (52,706 mi 2 ) are located in Mexico and 24,770,640 ha (95,639 mi 2 ) in the United States. We primarily used Landsat 8 (OLI) imagery, supplemented by limited ground surveys from two field campaigns, drone-based aerial data, and existing vegetation classification frameworks from both countries. The classification applied a data-fusion approach integrating 30-m Landsat 8 imagery, decadal phenology metrics from vegetation indices, and a random forest model trained mainly with datasets from a comprehensive national mapping project from the U.S. Geological Survey (USGS) GAP Analysis Project (GAP) and federal wildland fire agencies’ Landscape Fire and Resource Management Planning Tools (LANDFIRE) (GAP/LANDFIRE) [United States side] and the National Institute of Statistics and Geography (INEGI) [Mexico side] as well as land cover maps and opportunistic open-access and field observations. Mapping of the full BCR 33 region was carried out in two phases: 1) Phase I, the prototype map, covered a smaller portion of the transboundary area and identified 31 land cover classes, and 2) Phase II, the full BCR 33 map (refer to Figure 1), which resulted in 152 land cover classes. Using a Random Forest classifier, we achieved an overall prediction accuracy of 92% for the Phase I map and 87% for the Phase II full region map. This slight decrease can be attributed to working on a larger, more complex area with a greater number of land cover classes. No formal validation was conducted, aside from using a subset of the collected field observations and training data to assess model performance during and after training. The training sites were further verified using Google Earth (Google, 2026) imagery. Two undergraduate students who worked for over a year visually inspected imagery and open access public images to confirm each training site during model training using in-house developed, online, visual tools. A portion of this field training data was reserved for model validation, and the corresponding results are to be presented in later sections. The project developed an end-to-end, medium- and fine-resolution remote sensing–based data fusion mapping approach. This effort produced a map (Nagler et al., 2025) and the online tools to support a dynamic, live, online map for visualizing the transboundary vegetation communities in BCR 33. The toolset is currently hosted by the University of Arizona (UofA) Vegetation Index and Phenology (VIP) Lab to support FWS partners (https://vip.arizona.edu/viplab_data_explorer?LCM_BCR33). The online map is designed to allow rapid updates using new training, validation, or correction data, making it dynamic and maintainable. The approach we took established a framework for rapid updating and correction of land cover maps, as the model can be quickly retrained with new field observations, updated training data, or other sources. This enables dynamic mapping and change detection of the region’s vegetation. This framework is an advance in data fusion and crowdsourced mapping of complex, vulnerable regions, providing support to regional stakeholders and the wider user community. This transboundary map can inform the protection, conservation, and restoration of vegetation, habitat, and ecosystems, particularly for threatened and endangered species across the two nations using consistent and harmonized binational mapping systems. Beyond supporting land management decisions and stakeholders in the transboundary desert ecoregions, this BCR 33 mapping effort establishes a foundation for future rapid, low-cost, cross-border land cover mapping that can benefit and advance ecosystem management.

Arizona, Baja California, California, Nevada, Sina↗
USGS2026

New measurements indicate that natural geologic methane emissions from microseepage in the Michigan Basin are likely negligible

The magnitude of natural geologic methane (CH 4 ) emissions to the atmosphere (including emissions of fossil CH 4 from offshore and onshore gas and oil seeps, diffuse microseepage, mud volcanoes, volcanic vents, and geothermal areas) is highly uncertain. The largest component of geologic emissions is thought to be microseepage, which is the diffuse flux of CH 4 from soils across large areas of productive hydrocarbon basins. The accuracy of existing bottom-up estimates of microseepage is limited by low spatial coverage of published microseepage measurements. We present the first soil—atmosphere CH 4 flux measurements from Michigan Basin, USA. Results from 335 measurements taken during summer and winter seasons across a large portion of the basin suggest microseepage is nonexistent in the sampled region. Even areas with predictive features for microseepage (e.g., underlying mature, organically rich source rocks, proven gas accumulations, faults, and lineaments) yield null or negative fluxes, suggesting that CH 4 emissions from microseepage are negligible throughout our study region. A Monte Carlo method was used to place an upper limit on the regional-mean microseepage flux, in which synthetic patchy microseepage distributions were generated and tested against our observations to assess the impact of possible emission hot spots that were missed by our sampling strategy. Our analysis finds it is very unlikely that regional-mean emissions are as high as assumed in a previous global microseepage study. The observed lack of seepage may be explained by groundwater flow, active methanotrophy, glacial sediments, and bedded salt deposits, which could inhibit vertical gas migration and release to atmosphere.

Michigan↗