SEARCH · Geology Reports

Results for “Remote Sensing Applications: Society and Environment”

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

SourceResearch connectionFrom yearThrough year

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

1,659 records · Page 11Linked to original sources

USGS2025

Application of Hydrologic Simulation Program—FORTRAN (HSPF) as part of an integrated hydrologic model for the Salinas Valley, California

The U.S. Geological Survey (USGS), in cooperation with the Monterey County Water Resources Agency, conducted studies to help evaluate the surface-water and groundwater resources of the Salinas Valley study area, consisting of the entire Salinas River watershed and several smaller, adjacent coastal watersheds draining into Monterey Bay. The Salinas Valley study area is a highly productive agricultural region that depends on the coordinated use of surface water and groundwater to meet demand for irrigation and public water supply. To continue to meet these demands, a better understanding of the historical water balance and the effects of water-resource development on the long-term sustainability of water resources in the Salinas Valley study area is needed.

California↗

USGS2025

International data gaps at the Center for Engineering Strong Motion Data

The Center for Engineering Strong Motion Data (CESMD) is utilized by seismologists, engineers, and disaster management professionals in the US and has historically achieved and distributed waveforms from across the globe for significant earthquakes. The increased access to the waveforms via Web API (Application Programming Interface) offers a unique opportunity to provide the community complete datasets, sampling a variety of tectonic environments and geologic conditions, increasing the number of available ground motion records for use in ground motion models (GMMs) and improving the accuracy of earthquake engineering evaluations. The objective of this study is to programmatically identify gaps in global event data from the past decade and backfill missing data gaps at CESMD. We first compare the CESMD catalog with the Advanced National Seismic System (ANSS) Comprehensive Earthquake Catalog identifying regions and time periods where strong-motion data is limited or inadequate. To backfill datasets at CESMD for significant events, we pinpoint regions and time intervals that lack information, creating a list of events for which we’d like to obtain data. An important facet of this work is identifying the source of data and metadata across earthquake repositories around the world and integrating these data repositories into our current strong-motion data processing workflow. In parallel with these newly processed datasets, we are developing a script to produce data origination citations to include provenance and attribution information to associate with respective datasets at CESMD. We showcase our methodology for identifying and filling data gaps at CESMD using three case studies (the 2018 Anchorage Alaska earthquake sequence, seismicity associated with the 2018 Hawaiian Kilauea volcano eruption, and several earthquakes in Turkey) and then outline our strategy to apply our data gap backfilling methods on an international scale.

Conference Paper↗

USGS2026

Deep subsurface organic-rich shale supports abundant, diverse, and novel fungi

As Earth’s principal reservoir of organic carbon and microbial biomass, the deep subsurface hosts microorganisms capable of mobilizing this once-sequestered carbon. Contrary to standard assumptions of eukaryotic scarcity, this study documents abundant fungal communities, ranging from 4.2 × 10 3 to 6.8 × 10 3 fungal cells mL −1 , across a methane-producing organic-rich shale 247–556 meters below the surface. Although fungal:bacterial cell ratios ranged from 1:7028 to 1:713, application of biomass conversion factors developed for oceanic systems yielded a median fungal:bacterial biomass ratio of 1:4.7. 16S rRNA gene amplicons revealed bacterial and archaeal communities mirroring those found in well-characterized extremophilic, carbon-degrading environments, while sequencing of 18S rRNA gene and ITS rRNA spacer amplicons collectively identified a eukaryotic hotspot with 689 fungal OTUs across six phyla. The dominant fungal classes, Agaricomycetes and Dothideomycetes, are well-established degraders of recalcitrant carbon compounds at the surface, suggesting they may similarly contribute to organic matter degradation and ecosystem maintenance in the subsurface. Cultivation and isolation efforts yielded 205 fungal strains, including 13 candidate novel taxa, underscoring the deep subsurface as an underexplored eukaryotic habitat. Stable carbon isotopes indicate methane is predominantly generated via microbial conversion of the fossil carbon, while water isotopes suggest in situ geochemical conditions have been relatively stable since the Late Pleistocene, with subglacial recharge as a plausible mechanism for microbial introduction. Collectively, these findings suggest that fungi are underrecognized contributors to organic matter transformation and functional diversity in the deep biosphere, revealing a critical gap in our understanding of deep subsurface ecosystem processes.

Indiana, Michigan, Ohio↗

USGS2024

Recent applications of the USGS National Crustal Model for Seismic Hazard Studies

The U.S. Geological Survey is developing the National Crustal Model (NCM) for seismic hazard studies to facilitate modeling site, path, and source components of seismic hazard across the conterminous United States. The NCM is composed of a 1km grid of geophysical profiles, extending from the Earth’s surface into the upper mantle. It is constructed from a threedimensional (3D) geologic framework and geophysical rules that use (1) a petrologic and mineral physics database; (2) a 3D temperature model; and (3) a calibrated rock type- and age-dependent porosity model. Parameters needed to estimate site response for existing ground motion models (GMMs), including the time-averaged velocity in the upper 30 meters (VS30), the depths to 1.0 and 2.5 km/s shear-wave velocity (Z1.0 and Z2.5), and sediment thickness, can be computed from the NCM. As GMMs continue to improve in the future, other metrics could also be extracted or derived from the NCM, such as fundamental period, site attenuation (ko), a fully frequency-dependent site response function, or 3D geophysical volumes for wavefield simulations. Application of the NCM may also benefit other aspects of seismic hazard analysis, including better accounting for path-dependent attenuation and geometric spreading, more accurate estimation of earthquake source properties such as hypocentral location and stress drop, and calculation of crustal strength profiles that inform estimates of the base of seismicity.

conterminous United States↗

USGS2024

Reference 1D seismic velocity models for volcano monitoring and imaging: Methods, models, and applications

Seismic velocity models of the crust are an integral part of earthquake monitoring systems at volcanoes. 1D models that vary only in depth are typically used for real‐time hypocenter determination and serve as critical reference models for detailed 3D imaging studies and geomechanical modeling. Such models are usually computed using seismic tomographic methods that rely on P ‐ and S ‐wave arrival‐time picks from numerous earthquakes recorded at receivers around the volcano. Traditional linearized tomographic methods that jointly invert for source locations, velocity structure, and station corrections depend critically on having reasonable starting values for the unknown parameters, are susceptible to local misfit minima and divergence, and often do not provide adequate uncertainty information. These issues are often exacerbated by sparse seismic networks, inadequate distributions of seismicity, and/or poor data quality common at volcanoes. In contrast, modern probabilistic global search methods avoid these issues only at the cost of increased computation time. In this article, we review both approaches and present example applications and comparisons at several volcanoes in the United States, including Mount Hood (Oregon), Mount St. Helens (Washington), the Island of Hawai’i, and Mount Cleveland (Alaska). We provide guidance on the proper usage of these methods as relevant to challenges specific to volcano monitoring and imaging. Finally, we survey‐published 1D P ‐wave velocity models from around the world and use them to derive a generic stratovolcano velocity model, which serves as a useful reference model for comparison and when local velocity information is sparse.

Seismological Research Letters↗

USGS2025

Site-specific, extended ShakeMaps for earthquake engineering applications

The U.S. Geological Survey (USGS) routinely produces ShakeMaps of shaking intensity across the globe. Due to practical constraints, the number of response spectral periods was limited to three standard periods (0.3, 1.0, and 3.0 sec). We have recently developed the tools that are necessary to expand this functionality to include 22 periods (matching the current U.S. National Seismic Hazard Model periods) as well as the orientation-independent components (e.g., “RotD50”). We refer to ShakeMap products that include these extensions as “extended ShakeMaps.” The added level of complexity motivated us to also develop a user-friendly tool called the “ShakeMap Sampling Tool” (SST) that gives all the estimated shaking metrics for a specific location (or list of locations). Additionally, we develop a web application where users can input locations of interest and view/download the SST results. We further familiarize users with the concept of “Composite ShakeMaps.” For earthquakes sequences such as a mainshock and larger foreshocks and aftershocks, this provides a map of the maximum value of each shaking metric, which is useful for overall loss estimates, the full extent of ground failure triggering potential, and a better portrayal of the repeated shaking levels at a given point for a series of earthquakes. Such a site-specific shaking history facilitates earthquake forensics at building or infrastructure sites for which damage may be of concern, as described in the Disproportionate Damage Earthquake trigger specified in the IEBC (2018, Section 405.2.2) and in developing ATC-145 guidelines (Guidelines for Post-Earthquake Assessment, Repair, and Retrofit of Buildings). The composite ShakeMap can be combined with the SST for a variety of earthquake-hazard applications, such as systematically inferring triggering shaking estimates at specific sites of geotechnical interest for landsliding, liquefaction, and lateral-spreading hazards.

Kahramanmaraş↗

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↗

USGS2025

Modeling byproduct and coproduct mine production and mineral substitution using multidimensional supply curves: Application to the Cu-Co-Ni system and beyond

Rapid demand growth is expected for many metals used in the energy transition. Many of these metals are byproducts of other commodities. Byproduct production’s price response is tied to host mineral economics, complicating its supply dynamics. Moreover, many of these metals are used in applications where the material properties desired are difficult to substitute; effectively, limiting how quickly demand can adapt to changes in commodity price. Previous work has demonstrated the interconnectivity of jointly produced mineral commodities from the supply side, where the copper–cobalt–nickel system was used and demand was assumed independent across commodities. Studies to understand byproduct-coproduct market interconnectivity on the demand side are limited, while studies on the interconnectivity of supply and demand simultaneously are even more so. We propose a modification to the multicommodity supply curve method to enable inter-commodity effects on demand simultaneous with supply. In batteries, high cobaltprices may push consumers to transition to high-nickel chemistries, causing the nickel demand surface to decrease with nickel price but increase with cobaltprice, creating a two-dimensional demand surface. Below cross-price elasticities of 0.05, inter-commodity effects were found to be negligible, potentially permitting exclusion of these effects for many commodities. This additional demand curve complexity introduces potential computation challenges alongside the capacity to model many interrelated commodity systems such as rare earth elements, ferroalloys, country-oriented subsidies or restrictions, and bifurcated sustainable metals markets. By presenting the work done on multicommodity supply surfaces to date and potential new directions, this work aims to catalyze the next round of innovative approaches to modeling jointly produced commodities.

Conference Paper↗

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↗

USGS2025

Broadband stochastic simulation of earthquake ground motions with multiple strong phases with an application to the 2023 Kahramanmaraş, Turkey (Türkiye), earthquake

Stochastic ground motion simulation models are often less accurate at lower frequencies than at higher frequencies when fitting recorded data unless supplemented by a deterministic forward directivity velocity pulse model. Moreover, time-modulated stochastic models, which adjust ground motion amplitudes over time, typically use functions that fail to capture multiple strong-motion phases. The February 2023 Turkey (Türkiye) earthquake exhibited diverse recordings, including near-fault and far-field motions with pulse-like and non-pulse-like characteristics, along with single and multiple strong-motion phases. To better represent such a diverse set of recordings, this study enhances a fully non-stationary site-based stochastic model without combining it with a deterministic model. Improvements include a new band-pass filter with upper- and lower-frequency limits, which refines the representation of the low-frequency content. Moreover, a time-modulating function that can represent energy arrival in multiple strong phases is introduced. The reference model’s parameters are identified by fitting to the energy content, zero-level crossings, and cumulative counts of positive-minima and negative-maxima of a target accelerogram. This fitting procedure is modified to address the increased number of parameters. These improvements broaden the reference model’s applicability while preserving its simplicity, a key aspect appealing to engineering practitioners. The improved model’s applicability is demonstrated by simulating a dataset from the February 2023 Türkiye earthquake, and the accuracy is tested using a pulse-like Next Generation Attenuation Relationships for Western United States dataset. Validations are performed based on total energy, zero-level crossings, Fourier amplitude spectrum, elastic response spectra, and peak ground motion parameters. Validations are performed schematically in the time and frequency domains and quantitatively using goodness-of-fit scores, various validation-metrics errors, and inter-period correlations. Overall, the improved stochastic model can effectively simulate a set of diverse ground motion recordings, including near-fault pulse-like records, records with multiple strong phases, and far-field motions across a broad frequency range.

Kahramanmaraş↗

USGS2025

Towards mobile wind measurements using joust configured ultrasonic anemometer for applications in gas flux quantification

Small uncrewed aerial systems (sUASs) can be used to quantify emissions of greenhouse and other gases, providing flexibility in quantifying these emissions from a multitude of sources, including oil and gas infrastructure, volcano plumes, wildfire emissions, and natural sources. However, sUAS-based emission estimates are sensitive to the accuracy of wind speed and direction measurements. In this study, we examined how filtering and correcting sUAS-based wind measurements affects data accuracy by comparing data from a miniature ultrasonic anemometer mounted on a sUAS in a joust configuration to highly accurate wind data taken from a nearby eddy covariance flux tower (aka the Tower). These corrections had a small effect on wind speed error, but reduced wind direction errors from 50° to >120° to 20–30°. A concurrent experiment examining the amount of error due to the sUAS and the Tower not being co-located showed that the impact of this separation was 0.16–0.21 ms − 1 "> ms − 1 , a small influence on wind speed errors. Lower wind speed errors were correlated with lower turbulence intensity and higher relative wind speeds. There were also some loose trends in diminished wind direction errors at higher relative wind speeds. Therefore, to improve the quality of sUAS-based wind measurements, our study suggested that flight planning consider optimizing conditions that can lower turbulence intensity and maximize relative wind speeds as well as include post-flight corrections.

Alaska↗

USGS2026

Assessment of extreme subsurface hydrologic conditions captured during atmospheric river storms in the San Francisco Bay area (California, USA) with applications to shallow landslide early warning

An increase in soil pore water pressure is the typical trigger for the majority of landslides caused by rainfall. Atmospheric river storms, common to the west coast of North America during the winter season, can deliver landslide triggering rainfall resulting in severe impacts to coastal communities. Using a network of hydrologic monitoring stations situated within landslide-prone terrain in the San Francisco Bay area of California (USA), we assess the meteorologic conditions and resulting hydrologic and landslide response resulting from eight consecutive storm events that caused thousands of shallow landslides during the winter of 2022–2023. We find disparate hydrological responses and resultant degrees of observed landsliding ranging from < 1 landslide/km2 to 18 landslides/km2 at the monitoring sites that reflect the interplay and differences between rainfall delivery, subsurface hydrological characteristics, and geotechnical properties at each site. Antecedent soil moisture from both early season rainfall and the first storm in the sequence played a critical role in setting up some hillslopes for failure. Subsequent storms then generated elevated pore water pressures for several hours with associated landsliding. However, we find that the occurrence of widespread landsliding required not only sufficient pore water pressure magnitude in susceptible hillslopes, but also full and prolonged effective soil saturation throughout hillslope profiles. Landslides may still occur at lower values and durations of effective saturation but are likely to be less extensive regionally. We present these findings within the context of research directions and improvements to landslide early warning systems first suggested by researchers 40 years ago.

California↗

USGS2025

Deterministic physics-based earthquake sequence simulators match empirical ground-motion models and enable extrapolation to data poor regimes: Application to multifault multimechanism ruptures

We use the deterministic earthquake simulator RSQSim to generate complex sequences of ruptures on fault systems used for hazard assessment. We show that the source motions combined with a wave propagation code create surface ground motions that fall within the range of epistemic uncertainties for the Next Generation Attenuation‐West2 set of empirical models. We show the model is well calibrated where there are good data constraints, and has good correspondence in regions with fewer data constraints. We show magnitude, distance, and mechanism dependence all arising naturally from the same underlying friction. The deterministic physics‐based approach provides an opportunity for better understanding the physical origins of ground motions. For example, we find that reduced stress drops in shallow layers relative to constant stress drop with depth lead to peak ground velocities in the near field that better match empirical models. The simulators may also provide better extrapolations into regimes that are poorly empirically constrained by data because physics, rather than surface shaking data parameterizations, is underlying the extrapolations. Having shown the model is credible, we apply it to a problem where observations are lacking. We examine the case of crustal faults above a shallow subduction interface seen to break coseismically in simulations of the New Zealand fault system. These types of events were left out of consideration in the most recent New Zealand national seismic hazard model due to the modeling complexity and lack of observational data to constrain ground‐motion models (GMMs). Here, we show that in the model, by breaking up the coseismic crustal and interface rupturing fault motions into two separate subevents, and then recombining the resulting ground‐motion measures in a square‐root‐of‐sum‐of‐squares incoherent manner, we reproduce well the ground‐motion measures from the full event rupture. This provides a new method for extrapolating GMMs to more complex multifault ruptures.

Seismological Research Letters↗

USGS2026

Structural evolution of iron coordination in proteins across Earth’s oxygenation history

Protein metal-binding sites support essential biological functions shaped by protein fold, subunit interactions, and cofactor chemistry. Because these sites encode both biochemical function and environmental constraint, they offer a route to connect protein evolution with changes in Earth’s surface environment through time. Of particular interest is iron (Fe), the most widely used metal in biology and a cofactor central to both anaerobic and aerobic metabolism. Here, we systematically compare the immediate chemical environments of functionally essential Fe-binding sites in three-dimensional protein structures to test whether Fe coordination spheres differ across oxygen contexts. Using a curated dataset of experimentally determined structures, we identify a clear shift in the local chemistry of Fe-binding environments from anaerobic to aerobic proteins. Aerobic Fe sites are significantly more hydrophilic than anaerobic ones, and amino-acid composition analyses show reduced cysteine use in aerobic Fe-binding neighborhoods. These patterns suggest that as Earth’s surface environments became more oxygenated, proteins retained Fe as a core redox metal while reconfiguring local coordination chemistry in ways less vulnerable to oxidative damage. More broadly, this study introduces and applies the Coordination Sphere Analysis and Comparison (CSAC) workflow, an open and archived Python workflow for extracting local metal-binding environments from structure datasets, providing a framework for linking metalloprotein structure to evolutionary and geobiological transitions across Earth history.

Discover Life↗

USGS2026

Future water constraints on United States lithium mining under climate change

Lithium is necessary for low-carbon technologies that combat climate change, but lithium extraction is water-intensive. Changes in temperature and precipitation arising from climate change are altering water distribution, which could further strain supplies for new mines and industry, farms, and households. Here we explored how climate change, water use, and mining siting could impact lithium mining in the United States. We analyzed whether there would be sufficient water available to support the single existing and 22 proposed U.S. lithium mines at mid-century under four socioeconomic-climate scenarios and five climate models. Though dependent on socioeconomic-climate scenario, climate model, and lithium deposit type, available water supply in most subbasins would likely be unable to support new mines’ water demands, or even non-mining water demands from other sectors. Water scarcity could hinder the ability of the United States to produce enough lithium to meet domestic demand thereby necessitating higher imports.

conterminous United States↗

USGS2025

Memory and jamming in fault zone sediments

Many subsurface processes involve transitions in granular material states, from arrested to creeping to flowing. Experiments and frameworks for idealized systems reveal that granular fabrics develop during shearing, co-evolve with applied stress, and govern such transitions. We use microtomography to test whether fabrics at two San Andreas fault sites reflect slip history and whether idealized frameworks extend to nature. Near-surface sediments within the fault zone transition between deformation patterns over the seismic cycle, including bulk/localized grain re-arrangements, individual grain fracturing, and localized zones of fracturing. Aseismic and co-seismic shearing produce distinct preferred grain orientations. Co-seismic fabrics can be preserved after centuries of aseismic strain, aseismic fabrics may be overprinted, and grain size and coordination number influence the fabrics. Idealized frameworks, namely anisotropic critical state theory, frictional jamming, and material memory, can explain our observations, and fault zone sediments likely undergo cycles of memory creation and erasure that influence rigidity spatiotemporally.

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

Mapping potential sensitivity to hydrogeomorphic change in the UMRS riverscape

In 2020 the U.S. Geological Survey (USGS), as part of the Upper Mississippi River Restoration (UMRR) Program, began a new project to characterize potential hydrogeomorphic change associated with hydrogeomorphic units (HGUs) and their catenae (units linked by their association with sediment sources and flow origins). The goal of the project was to develop a geographic information system (GIS) database of HGUs for the Upper Mississippi River System (UMRS) available to both scientists and river managers working on UMRR studies and HREP planning and design studies. The characterization was based on a hydrogeomorphic change hierarchical classification developed previously for the UMRS. The products were generated with automated techniques in a GIS using systemic datasets. Landforms were mapped from the 2015 UMRS topobathymetric dataset with geomorphon (shape-based) tools tailored for the large riverscape, valley bottom environments in the UMRS. A clustering analysis was applied to the resulting landforms to identify HGUs associated with zones of perennial low flows, bankfull flows, and overbank floodplains. Catenae were assembled based on the proximity of the units to the main channel, tributary mouths, and side channels from previously published aquatic areas (USACE, 2018) coupled with least-cost flowpath linkages between potential sediment origins and planform change units developed by Rogala, Fitzpatrick, and Henderson (2020). These GIS-based analyses were successful at identifying a range of HGUs using an automated technique with available data across the entire riverscape, with emphasis on those that have the potential for hydrogeomorphic change. Most of the resulting features are depositional, as expected in a large river system. However, this is the first attempt of linking tributary inputs, side channel erosion and levee breaches with their depositional counterparts. The approach was successfully piloted in Pools 8 and 10 in the Upper Impounded Reach and Pool 14 in the Lower Impounded Reach, with next steps for application in reaches of the unimpounded section and Illinois River. This report emphasizes results from Pool 10, which was the focus of most of our attention during the pilot phase.

Illinois, Indiana, Iowa, Minnesota, Missouri, Sout↗