Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Earth System Science Data”

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,193 records · Page 48Linked to original sources
USGS2026

U.S. Geological Survey Monitoring Milestones—Oe-151 at Woodgate, NY (433112075091501)

On July 9, 1926, monitoring well Oe-151 at Woodgate, New York (USGS ID 433112075091501) recorded its first groundwater data. Since then, the well has provided water data nearly continuously and has now reached a 100-year milestone for data collection. The well is part of the U.S. Geological Survey (USGS) Climate Response Network (CRN), which is a national network of wells selected to monitor natural groundwater conditions. Well Oe-151 is the first well in the network to reach a 100-year Centennial milestone.

New York↗
USGS1984

Techniques for estimating magnitude and frequency of floods on streams in Indiana

Equations are presented for estimating the magnitude and frequency of floods at ungaged sites on unregulated and nonurban streams in Indiana. The equations were developed by multiple-regression, analysis of basin characteristics and peak-flow statistical data from 242 gaged locations in Indiana, Ohio, and Illinois. The State of Indiana was divided into seven areas on the basis of the regression analysis. A set of equations for estimating peak discharges with recurrence intervals of 2, IO, 25, 50, and 100 years was developed for each area. Significant basin characteristics in the equations are drainage area, channel length, channel slope, mean annual precipitation, storage, precipitation intensity, and a runoff coefficient. Standard errors of estimate for the equations range from 24 to 45 percent. Methods are also presented for estimating flood magnitude and frequency at sites on gaged streams. Flood-frequency data based on observed peaks are given for 270 gaged locations. Twenty of these are on regulated streams, and six are on urban streams. Basin characteristics are also included car 245 of the gaged locations on unregulated and nonurban streams. No techniques are given for estimating flood magnitude and frequency at ungaged sites on regulated or urban streams. A rainfall-runoff model was tlsed to synthesize long-term peak data at 11 gaged locations on small streams. Flood-frequency curves developed from the long-term synthetic data were combined with curves based on short-term observed data to provide weighted estimates of flood magnitude and frequency at the rainfall-runoff stations.

Indiana↗
USGS2025

Statistical approaches for modeling correlated grade and tonnage distributions and applications for mineral resource assessments

Correlations between grade and tonnage exist in mineral resource data compiled from published reports, but they are not always addressed during quantitative assessment of undiscovered mineral resources. Failure to account for correlated grade and tonnage distributions can result in geologically unrealistic assessment results. Current software tools simulate univariate ore tonnage and multivariate resource grades of undiscovered deposits independently. As a result, analysts are forced to rely on ad-hoc solutions to minimize the correlation issues by: 1) creating subsets of data with restricted criteria; 2) truncating grade and tonnage distributions; and 3) testing model robustness using exploratory data analysis. While these methods represent pragmatic solutions, the statistical solutions presented here provide additional options to address real correlations in grade and tonnage data used for mineral resource assessments. We present a modified version of the MapMark4 package in R that introduces two alternatives for modeling grade and tonnage distributions, consisting of a multivariate solution that accounts for correlations between ore tonnage and metal grades and an empirical solution that utilizes simple random sampling with replacement to reproduce coupled grades and tonnages from the input data. We present simulations for contained ore and metal for three case studies representing tungsten skarn, komatiite-hosted nickel, and sediment-hosted carbonate amagmatic zinc-lead (Mississippi Valley-type) deposits. Employing the methods presented here yields quantitative mineral resource assessment results that more closely reflect the empirical distributions of grades and tonnages observed in nature and expands the applicability of these tools for ongoing critical mineral resource assessments.

Applied Computing and Geosciences↗
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↗
USGS2026

A case study of maximum depositional ages from terrestrial sandstones near the Cretaceous—Paleogene transition, western Williston Basin, USA

We present a new Bayesian method for deriving maximum depositional ages from detrital K-feldspar using total-fusion 40 Ar/ 39 Ar data. Individual analyses weighted by K/Ca ratio, age uncertainty, and percent radiogenic argon prioritize a result to come from accurately measured volcanic sources. Results from this method applied to sandstones from the Hell Creek region of northeastern Montana show that detrital K-feldspar maximum depositional ages align closely with detrital U-Pb zircon ages as well as tephra based chronostratigraphic constraints from both 40 Ar/ 39 Ar and U-Pb data. An age model informed by maximum depositional ages and available tephra data provides an estimate of 1.74 ± 1.04 Ma for the duration of deposition for the Hell Creek Formation. Combined age model and maximum depositional age data suggest ≤ 1.84 Ma of missing time is represented by the Hell Creek and Fox Hills formational disconformity, coinciding with the final regression of the Western Interior seaway.

Montana↗
USGS2026

High-resolution landscape-scale maps of nutrients, phytoplankton, and related water-quality constituents in Suisun Bay and the Sacramento–San Joaquin Delta, California, during 2018—Results from high-resolution underway surveys

We examined the abundance and distribution of nutrients and phytoplankton in the tidal aquatic environments of the Sacramento–San Joaquin Delta (Delta) and Suisun, Grizzly, and Honker Bays in California by completing three spatial surveys that used continuous underway high frequency (1 hertz) measurements and sampling onboard a high-speed boat. Surveys were conducted in May, July, and October 2018. Water was sampled continuously during surveys to simultaneously collect information about the concentration and spatial distribution of all major nutrient forms together with information about the major classes of phytoplankton and associated standard field measurements of water quality, such as dissolved oxygen, dissolved organic matter, pH, salinity, specific conductance, turbidity, and water temperature. The results show a greater than 50-fold variation in nutrient and phytoplankton concentrations across space and time, providing evidence of the dynamic environmental processes that shape the ways nutrients interact with and affect Delta aquatic habitats.

Data Report↗
USGS2024

Are field observations of surface rupture useful? An example from the 2023 Mw 7.8 Pazarcık, Turkey (Türkiye), earthquake

Field investigations have long been an important component of the scientific response to surface‐faulting earthquakes. However, in light of advances in remote data and models, the question arises whether field‐based observations of surface rupture remain useful for understanding rupture processes and seismic hazards. We approach this question using a field‐based study of the central 2023 M w 7.8 Pazarcık, Turkey (Türkiye), earthquake rupture, at the intersection of the east Anatolian fault (EAF) and Narlı fault. Our field observations include the surface rupture expression and extent of the central EAF and northernmost Narlı fault in generally forested and steep terrain and 68 measurements of left‐lateral surface displacement. These data improve our understanding of the Pazarcık rupture complexity, resolve the surface geometry of the Narlı fault–EAF connection, and exhibit a clear (>2 m) change in surface displacement across this intersection zone that confirms remote‐based coseismic slip models. Our study shows that focusing field efforts in areas of obscured or low‐resolution remote data can yield essential data for refining rupture extent, documenting perishable on‐fault displacement, and improving postearthquake situational awareness. A comparison of similarly large‐magnitude continental surface‐rupturing earthquakes indicates that displacement uncertainties relate to a complex set of factors, including measurement methods, rupture complexity, and displacement magnitude. Our study validates the need for postearthquake field observations, which, when driven by clear motivating questions and knowledge of methodological strengths and limitations, provide high‐resolution rupture data that complement remote‐based models.

Pazarcık↗
USGS2025

The 3D National Topography Model Call for Action—Part 2: The Next Generation 3D Elevation Program

The three-dimensional (3D) National Topography Model initiative to integrate elevation and hydrography data includes the next generation of hydrography data from the 3D Hydrography Program and the next generation of elevation data from the 3D Elevation Program (3DEP). The first-ever collection of light detection and ranging (lidar) data for the nation (IfSAR for Alaska) provides a critical baseline reference, and the addition of multiple repeat elevation mapping projects as part of the next generation of 3DEP would substantially expand analysis capabilities. As the U.S. Geological Survey (USGS) is closing in on our goal of complete coverage with 3DEP data available or in progress for 98.3 percent of the Nation at the end of fiscal year 2024, the USGS is already transitioning to the next generation of 3DEP. Based on the 3D Nation Study results and input from a broad range of stakeholders, the USGS National Geospatial Program has finalized a new design for 3DEP that provides increased lidar quality levels and refresh rates. The new program is designed with more flexibility to meet changing user needs and take advantage of improvements in mapping technologies. The program will aim to expand the level of interagency coordination for topobathymetric lidar acquisition for inland rivers. The next generation of 3DEP will also aim to emphasize research, including advancing program design, products, and services and engaging and leveraging the evolving 3D industry. Research goals also include becoming more flexible in meeting user needs and taking advantage of evolving remote-sensing technologies. The program also plans to move from focusing on producing standard products to producing a concept of a 3D Nation Ecosystem with a variety of inputs, products, and services.

Circular↗
USGS2026

Streamflow, water quality, and constituent loads and yields, Scituate Reservoir drainage area, Rhode Island, water year 2023

As part of a long-term cooperative program to monitor water quality within the Scituate Reservoir drainage area, the U.S. Geological Survey, in cooperation with Providence Water (formerly the Providence Water Supply Board), collected streamflow and water-quality data in tributaries to the Scituate Reservoir, Rhode Island. Streamflow and concentrations of chloride and sodium estimated from records of specific conductance for 16 tributaries were used to calculate loads of chloride and sodium during water year 2023 (October 1, 2022, through September 30, 2023). Water-quality samples were collected by Providence Water at 38 sampling stations on tributaries to the Scituate Reservoir during water year 2023. These water-quality data are summarized by using values of central tendency and are used, in combination with measured (or estimated) streamflows, to calculate loads and yields of selected water-quality constituents for water year 2023. Annual mean streamflows for monitoring stations in this study ranged from 0.47 to 34.4 cubic feet per second during water year 2023. At the 16 continuous-record streamgages, tributaries transported about 3,000 metric tons of chloride and 1,900 metric tons of sodium to the Scituate Reservoir; annual chloride yields for the tributaries ranged from 19 to 200 metric tons per square mile, and annual sodium yields ranged from 13 to 120 metric tons per square mile. At the stations where water-quality samples were collected by Providence Water, the medians of the median daily loads were 65,000 million colony forming units per day for coliform bacteria, 870 million colony forming units per day for Escherichia coli , 200 kilograms per day for chloride, 8.2 grams per day as nitrogen for nitrite, 520 grams per day as nitrogen for nitrate, and 300 grams per day as orthophosphate for phosphate. The medians of the median yields were 25,000 million colony forming units per day per square mile for coliform bacteria, 600 million colony forming units per day per square mile for Escherichia coli , 65 kilograms per day per square mile for chloride, 4.6 grams per day per square mile as nitrogen for nitrite, 180 grams per day per square mile as nitrogen for nitrate, and 180 grams per day per square mile as orthophosphate for phosphate.

Rhode Island↗
USGS2025

High-pass corner frequency selection and review tool for use in ground-motion processing

Raw seismological waveform data contain noise from the instrument’s surroundings and the instrument itself that can dominate recordings at low and high frequencies. To use these data in ground‐motion modeling, the effects of noise on the signals must be reduced and the signals’ usable frequency range identified. We present automated procedures to efficiently reduce low‐frequency noise that are implemented in the software package gmprocess. These procedures check for, and as needed remove, low‐frequency artifacts in the displacement record using polynomial fits, which can be used in combination with existing signal‐to‐noise ratio (SNR)‐based corner‐frequency selection procedures. The automated selections are then efficiently verified and refined using a graphical user interface (GUI) that plots relevant ground‐motion time series and spectra and tracks modifications to signal processing parameters. We demonstrate these procedures using recordings from the 2020 M 5.1 Sparta, North Carolina, and the 2013 M 4.7 southern Ontario earthquakes. Data processed with the SNR‐only and polynomial criteria for these events contain displacement artifacts in 37% and 23% of processed traces, respectively. Records with remaining artifacts are corrected manually using the GUI. These processing steps illustrate the workflow for efficient data processing with quality control.

Seismological Research Letters↗
USGS1997

Application of acoustical methods for estimating water flow and constituent loads in Perdido Bay, Florida

Water flow and quality data were collected from December 1994 to September 1995 to evaluate variations in discharge, water quality, and chemical fluxes (loads) through Perdido Bay, Florida. Data were collected at a cross section parallel to the U.S. Highway 98 bridge. Discharges measured with an acoustic Doppler current profiler (ADCP) and computed from stage-area and velocity ratings varied roughly between + or - 10,000 cubic feet per second during a typical tidal cycle. Large reversals in flow direction occurred rapidly (less than 1 hour), and complete reversals (resulting in near peak net-upstream or downstream discharges) occurred within a few hours of slack water. Observations of simultaneous upstream and downstream flow (bidirectional flow) were quite common in the ADCP measurements, with opposing directions of flow occurring predominantly in vertical layers. Continuous (every 15 minutes) discharge data were computed for the period from August 18, 1995, to September 28, 1995, and filtered daily mean discharge values were computed for the period from August 19 to September 26, 1995. Data were not computed prior to August 18, 1995, either because of missing data or because the velocity rating was poorly defined (because of insufficient data) for the period prior to landfall of hurricane Erin (August 3, 1995). The results of the study indicate that acoustical techniques can yield useful estimates of continuous (instantaneous) discharge in Perdido Bay. Useful estimates of average daily net flow rates can also be obtained, but the accuracy of these estimates will be limited by small rating shifts that introduce bias into the instantaneous values that are used to compute the net flows. Instantaneous loads of total nitrogen ranged from -180 to 220 grams per second for the samples collected during the study, and instantaneous loads of total phosphorous ranged from -10 to 11 grams per second (negative loads indicate net upstream transport). The chloride concentrations from the water samples collected from Perdido Bay indicated a significant amount of mixing of saltwater and freshwater. Mixing effects could greatly reduce the accuracy of estimates of net loads of nutrients or other substances. The study results indicate that acoustical techniques can yield acceptable estimates of instantaneous loads in Perdido Bay. However, estimates of net loads should be interpreted with great caution and may have unacceptably large errors, especially when saltwater and freshwater concentrations differ greatly.

Florida↗
USGS2026

Quantitative mineral resource assessment of lithium pegmatite deposits in the southern Appalachian orogen

The first quantitative mineral resource assessment for undiscovered lithium pegmatite deposits in the southern Appalachian region of the United States was conducted. Permissive tracts for lithium pegmatite deposits were delineated by integrating lithological, tectonic, geochemical, geophysical and mineral occurrence data. Lithium pegmatite prospectivity of the tracts was ranked with simplified mappable criteria, including proximity to Paleozoic felsic intrusions and major lithotectonic structures, stream sediment geochemical anomalies, and pegmatite occurrence data. The geospatial data and permissive tracts were used to estimate the number of undiscovered lithium pegmatite deposits. These estimates were integrated into probabilistic simulations along with a new global lithium pegmatite grade and tonnage dataset to quantify potential contained undiscovered lithium resources. An economic filter was applied to convert the probabilistic estimates of contained lithium into recoverable material. The identified lithium pegmatite resources for the Carolina Lithium and Kings Mountain deposits, North Carolina, contain 1589 thousand tons (kt) of Li 2 O. The median contained undiscovered resource for the southern Appalachian orogen was estimated to be 2240 kt Li 2 O. At 90% confidence, the region contains at least 130 kt Li 2 O, and 10,700 kt at 10% confidence. After applying economic filters, the median recoverable contained resource was 1430 kt Li 2 O, corresponding to approximately 201 years of current lithium imports for consumption in the United States. North and South Carolina are likely to contain most of these resources. Coarse data resolution and intra-state variations in the geological data contribute to uncertainty of undiscovered lithium pegmatite resources. Continued efforts to harmonize disparate geospatial datasets with updated or new information can improve the accuracy and precision of estimated undiscovered lithium pegmatite resources in the study area and at broader scales.

Alabama, Georgia, Maryland, North Carolina, South ↗
USGS2025

Multi-scale predictors of Northern Long-eared Bat (Myotis septentrionalis) occupancy in the United States

Historically, Myotis septentrionalis (Northern Long eared Bat) was among the most common forest-interior species in North America. Largely due to high mortality from white-nose syndrome, this species has experienced severe population declines across its range. To create an updated species distribution map representing summer occupancy probabilities from 2017 to 2022, we integrated stationary acoustic data with live-capture data from the database of the North American Bat Monitoring Program into a multi-scale, multi-method occupancy modeling framework. Our results provide data-driven predictions with quantified uncertainty for summer occupancy probabilities for Northern Long-eared Bats at 2 spatial scales across the range of the species, while also accounting for inherent observation biases (e.g., imperfect detection).

Journal of North American Bat Research↗
USGS2024

Introduction and methods of analysis for peak streamflow trends and their relation to changes in climate in Illinois, Iowa, Michigan, Minnesota, Missouri, Montana, North Dakota, South Dakota, and Wisconsin

Flood-frequency analysis, also called peak-flow frequency or flood-flow frequency analysis, is essential to water resources management applications including critical structure design and floodplain mapping. Federal guidelines for doing flood-frequency analyses are presented in a U.S. Geological Survey Techniques and Methods Report known as Bulletin 17C. A basic assumption within Bulletin 17C is that for drainage basins without major hydrologic alterations, statistical properties of the distribution of annual peak streamflows (peak flows) are stationary; that is, the mean, variance, and skew are constant. The stationarity assumption has been widely accepted within the flood-frequency community; however, a better understanding of long-term climatic persistence and concerns about potential climate change and land-use change has caused a reexamination of the stationarity assumption. Flood-frequency analyses that do not incorporate observed trends and abrupt changes may result in a poor representation of the true flood risk. Bulletin 17C does not offer guidance for incorporating nonstationarities when estimating floods, and it describes a need for studies that incorporate changing climate or basin characteristics. In response to this need and a history of concern regarding nonstationarity peak flows in the region, this study was done by the U.S. Geological Survey, in cooperation with the Departments of Transportation of Illinois, Iowa, Michigan, Minnesota, Missouri, South Dakota, and Wisconsin; the Montana Department of Natural Resources and Conservation; and the North Dakota Department of Water Resources, to assess potential nonstationarity in peak flows in the north-central United States. This chapter summarizes the methods used to detect hydroclimatic changes in peak-flow data in the study region. A wide range of analyses and statistical approaches are applied to document the primary mechanisms controlling floods and characterize temporal changes in hydroclimatic variables and peak flow. Four periods were selected for analysis of peak flow, daily streamflow, and climate data. The periods are (1) a 100-year period, 1921–2020; (2) a 75-year period, 1946–2020; (3) a 50-year period, 1971–2020; and (4) a 30-year period, 1991–2020. The climate data consist of monthly time series estimates of temperature, precipitation, potential evapotranspiration, actual evapotranspiration, snowfall, soil moisture storage, snow water equivalent, and runoff on a 3.1-mile by 3.1-mile grid for the conterminous United States. Statistical and graphical analyses were used to investigate potential changes in hydrology and climate. The starting point for these analyses was the initial data analysis of peak flow described in Bulletin 17C, which includes plotting the peak flow and checking for autocorrelation, monotonic trends, and changes points. Analyses were added to examine additional features in the data. To examine potential causal drivers of changes, the climate data were analyzed graphically and statistically. Results are provided in a U.S. Geological Survey data release. The study limitations are documented for users of the results.

Iowa, Illinois, Michigan, Missouri, Minnesota, Mon↗
USGS2024

Peak streamflow trends and their relation to changes in climate in Illinois, Iowa, Michigan, Minnesota, Missouri, Montana, North Dakota, South Dakota, and Wisconsin

Flood-frequency analysis, also called peak-flow frequency or flood-flow frequency analysis, is essential to water resources management applications including critical structure design and floodplain mapping. Federal guidelines for doing flood-frequency analyses are presented in a U.S. Geological Survey Techniques and Methods report known as Bulletin 17C. A basic assumption within Bulletin 17C is that for drainage basins without major hydrologic alterations, statistical properties of the distribution of annual peak streamflows (peak flows) are stationary; that is, the mean, variance, and skew are constant. The stationarity assumption has been widely accepted within the flood-frequency community; however, a better understanding of long-term climatic persistence and concerns about potential climate change and land-use change has caused a reexamination of the stationarity assumption. This work is part of that reexamination. The stationarity assumption is a concern because flood-frequency analyses that do not incorporate observed trends and abrupt changes may result in a poor representation of the true flood risk. Bulletin 17C does not offer guidance for incorporating nonstationarities when estimating floods, and it describes a need for studies that incorporate changing climate or basin characteristics. In response to this need and a history of concern regarding nonstationarity peak flows in the region, this study was done to assess potential nonstationarity in peak flows in the north-central United States. This report summarizes the methods used to detect hydroclimatic changes in peak-flow data in the study region. Four periods were selected for analysis of peak flow, daily streamflow, and climate data. The periods are (1) a 100-year period, 1921–2020; (2) a 75-year period, 1946–2020; (3) a 50-year period, 1971–2020; and (4) a 30-year period, 1991–2020. The starting point for these analyses was the initial data analysis of peak flow described in Bulletin 17C, which includes plotting the peak flow and checking for autocorrelation, monotonic trends, and changes points. Analyses were added to examine additional features in the data. Results are provided in a U.S. Geological Survey data release. The study limitations are documented for users of the results.

Scientific Investigations Report↗
USGS2024

Stream nitrate dynamics driven primarily by discharge and watershed physical and soil characteristics at intensively monitored sites: Insights from deep learning

We developed a suite of models using deep learning to make hindcast predictions of the 7‐day average backward‐looking nitrate concentration at 46 predominantly agricultural sites across the midwestern and eastern United States. The models used daily observations of discharge and meteorological variables and watershed attributes describing anthropogenic modification to hydrology, nitrogen application, climate, groundwater, land use, watershed physiographic attributes, and soils. Across all sites, discharge and watershed soil and physiographic attributes showed a strong influence on model performance. Analysis of drivers across sites revealed considerable regional differences related to controlling processes such as groundwater contributions. We tested several ways to pool data across sites to develop accurate models and make the most effective use of available data. Single‐site models, in which models are trained and tested at a single location, showed generally strong predictive performance (median Kling‐Gupta Efficiency = 0.66), and accuracy at poorly performing sites could be improved by grouping sites with similar characteristics. Developing a single model for all sites reduced performance at several locations with distinct characteristics, suggesting that there is a threshold of dissimilarity beyond which more data does not improve the model. While many deep learning studies have shown that national or even global models can outperform local models, it is not clear that this is true for water quality constituents. This study demonstrates how data can be combined effectively, using deep learning to develop accurate and interpretable models of instream nitrate at sites where varying processes are responsible for changes in nitrate concentration.

Water Resources Research↗
USGS2025

Multi-lithofacies alluvial characterization via airborne electromagnetic-borehole fusion using ordinary interval kriging and geologic constraints

Airborne electromagnetic (AEM) survey provides extensive spatial coverage and detailed resolution at the near surface and can be used to develop hydrogeological models. However, utilization of AEM data is not straightforward because AEM resistivity is an indirect measurement for inferring sediment types. This study develops an ordinary interval kriging (OIK) algorithm and a resistivity-to-multi-lithofacies (R2ML) data fusion workflow for multi-lithofacies alluvial characterization. OIK utilizes irregular interval data to construct three-dimensional (3D) resistivity fields from one-dimensional inverted AEM resistivity models. The R2ML workflow maps the resistivity field generated from OIK into a multi-facies lithological model, incorporating geologic constraints derived from well logs and geological observations. The numerical and real-world cases demonstrate that OIK is computationally efficient, accounts for 3D anisotropy, and minimizes the smoothing effect, thereby preserving resistivity contrasts and reducing interpolation uncertainty. The methodology is applied to lithologic characterization of the Mississippi River Valley alluvial aquifer (MRVA) in the Shellmound area, Mississippi, U.S. A frequency-domain AEM survey was conducted to support groundwater studies for the managed aquifer recharge (MAR) to the MRVA. The resulting lithological model, including four types of lithofacies—clay, very fine sand, fine-medium sands, and graveliferous sands, illustrates the geomorphological processes of the MRVA and implies potential MAR. The alignment between the lithological model and existing geological and hydrogeological investigations demonstrates that OIK and R2ML workflow effectively capture the subsurface architecture of the MRVA. The methods have broad applicability for characterizing alluvial aquifers through AEM-borehole data fusion, supporting sustainable groundwater management.

Mississippi↗
USGS2026

Is satellite-derived bathymetry vertical accuracy dependent on satellite mission and processing method?

This research focusses on three satellite-derived bathymetry methods and optical satellite instruments: (1) a stereo photogrammetry bathymetry module (SaTSeaD) developed for the NASA Ames stereo pipeline open-source software (version 3.6.0) using stereo WorldView data; (2) physics-based radiative transfer equations (PBSDB) using Landsat data; and (3) a modified composite band-ratio method for Sentinel-2 (SatBathy) with an initial simplified calibration, followed by a more rigorous linear regression against in situ bathymetry data. All methods were tested in three different areas with different geological and environmental conditions, Cabo Rojo, Puerto Rico; Key West, Florida; and Cocos Lagoon and Achang Flat Reef Preserve, Guam. It is demonstrated that all satellite derived bathymetry (SDB) methods have increased accuracy when the results are aligned with higher-accuracy ICESat-2 ATL24 track bathymetry data using the iterative closest point (ICP). SDB vertical accuracy depends more on location characteristics than the method or optical satellite instrument used. All error metrics considered (mean absolute error, median absolute deviation, and root mean square error) can be less than 5% of the maximum bathymetry depth penetration for at least one method, although not necessarily for the same method for all sites. The SDB error distribution tends to be bimodal irrespective of method, satellite instrument, alignment, site, or maximum bathymetry depth, leading to the potential ineffectiveness of traditional error metrics, such as the root mean square error. However, our analysis demonstrates that performing detrending where possible can achieve an error distribution as close to normality as possible for which error metrics are more diagnostic.

Florida↗