Relations between rainfall–runoff-induced erosion and aeolian deposition at archaeological sites in a semi-arid dam-controlled river corridor

Brian D. Collins; David Bedford; Skye C. Corbett; Helen C. Fairley; Collin Cronkite-Ratcliff

doi:10.1002/esp.3874

Home/Archive/USGS

USGS · 70173713

Relations between rainfall–runoff-induced erosion and aeolian deposition at archaeological sites in a semi-arid dam-controlled river corridor

Brian D. Collins·David Bedford·Skye C. Corbett·Helen C. Fairley·Collin Cronkite-Ratcliff

View original source ↗Export metadata ↓

Abstract

Process dynamics in fluvial-based dryland environments are highly complex with fluvial, aeolian, and alluvial processes all contributing to landscape change. When anthropogenic activities such as dam-building affect fluvial processes, the complexity in local response can be further increased by flood- and sediment-limiting flows. Understanding these complexities is key to predicting landscape behavior in drylands and has important scientific and management implications, including for studies related to paleoclimatology, landscape ecology evolution, and archaeological site context and preservation. Here we use multi-temporal LiDAR surveys, local weather data, and geomorphological observations to identify trends in site change throughout the 446-km-long semi-arid Colorado River corridor in Grand Canyon, Arizona, USA, where archaeological site degradation related to the effects of upstream dam operation is a concern. Using several site case studies, we show the range of landscape responses that might be expected from concomitant occurrence of dam-controlled fluvial sand bar deposition, aeolian sand transport, and rainfall-induced erosion. Empirical rainfall-erosion threshold analyses coupled with a numerical rainfall–runoff–soil erosion model indicate that infiltration-excess overland flow and gullying govern large-scale (centimeter- to decimeter-scale) landscape changes, but that aeolian deposition can in some cases mitigate gully erosion. Whereas threshold analyses identify the normalized rainfall intensity (defined as the ratio of rainfall intensity to hydraulic conductivity) as the primary factor governing hydrologic-driven erosion, assessment of false positives and false negatives in the dataset highlight topographic slope as the next most important parameter governing site response. Analysis of 4+ years of high resolution (four-minute) weather data and 75+ years of low resolution (daily) climate records indicates that dryland erosion is dependent on short-term, storm-driven rainfall intensity rather than cumulative rainfall, and that erosion can occur outside of wet seasons and even wet years. These results can apply to other similar semi-arid landscapes where process complexity may not be fully understood.

Explore related subjects

Earth Surface Processes and Landforms Article

Source-reported bounding extent: 35.715298012125295° to 36.910372213522535° latitude; -113.9886474609375° to -111.5277099609375° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

Brian D. Collins, David Bedford, Skye C. Corbett, Helen C. Fairley, Collin Cronkite-Ratcliff. 2016-01-14. Relations between rainfall–runoff-induced erosion and aeolian deposition at archaeological sites in a semi-arid dam-controlled river corridor. https://doi.org/10.1002/esp.3874

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

USGS2026

Deep learning models for mapping surficial geology in selected physiographic regions of New York

Surficial geologic mapping is required for decisions including infrastructure and water-resource development and management. Deep learning, a type of machine learning that uses training data to self-learn and perform tasks, is explored as a tool to help increase efficiency of the labour- and time-intensive mapping process. Deep learning models were trained, and their potential to aid in mapping surficial geology was explored for two physiographic regions in New York: the high-relief Allegheny Plateau and the low-relief Erie-Ontario Lowlands. Key to the development of deep learning models is the availability of highly detailed surficial geologic maps in each of these regions from which the models can learn. Through exploring different groupings of surficial deposits and their corresponding spatial data, deep learning models were iteratively developed to reproduce published mapping for greater than 79% of the training areas, and with similar accuracy in test areas within the same physiographic region. These models were trained using only two inputs: high-resolution lidar data and previously published surficial geologic maps. This straightforward approach is intended to make the methods and models created easily reproducible using widely accessible datasets. This research describes the strengths and limitations of lidar-derived surficial geologic models and how broadly grouping surficial types, characteristics of different physiographic regions, and testing models in physiographic regions outside of their training areas affect model creation and performance. In addition, the models created could be used as a tool for mapping surficial geology in similar physiographic regions, in addition to establishing a framework for creating similar models elsewhere.

New York↗

USGS2026

Regional models for postfire debris-flow likelihood and rainfall thresholds across the western United States

The U.S. Geological Survey (USGS) uses an empirical model developed with logistic regression (the ‘M1’ model) to rapidly assess debris-flow likelihood and to identify quantitative rainfall thresholds for debris flows after wildfire in the western United States. The M1 model was calibrated to a debris-flow inventory from southern California (United States) and has been applied throughout the western United States. Limited spatial coverage in the calibration dataset has motivated evaluation of M1 model accuracy outside the calibration region (e.g., the Sierra Nevada or the eastern Cascade Range, United States). Previous test cases showed that M1 overpredicts debris-flow likelihood and underpredicts rainfall thresholds for some locations (e.g., Arizona, northern California, Colorado, New Mexico, United States). We sought to improve the regional applicability of a debris-flow likelihood model by expanding the debris-flow inventory used for calibration, testing multiple potential models and generating an updated model framework. The updated inventory includes 3788 observations from 67 burned areas paired with short duration rainfall ratios. The updated model framework consists of a modified model structure and sets of coefficients calibrated separately to the entire updated inventory and to subsets of the inventory that intersect three Environmental Protection Agency (EPA) Level 2 ecoregions (Mediterranean California, Upper Gila Mountains and Western Cordillera). Comparisons of predictions from the updated models with observed rainfall and debris-flow activity show that the updated models outperform the M1 model by ~15%–60% and improve the uniformity of predictive performance across the western United States. The updated models also reduce false positive rates relative to M1 and generate rainfall thresholds that are better aligned with relative differences in regional climatology and debris-flow activity.

Arizona, California, Colorado, Idaho, Montana, Nev↗

USGS2026

Numerical modeling of Late Pleistocene to Holocene earthquake-induced progressive rock slope damage in the paraglacial Serpentine valley, Prince William Sound, Alaska

Landslides in deglaciating fjords pose a potential tsunami threat to nearby communities; however, processes contributing to long-term progressive rock damage and landslide conditioning remain poorly constrained in paraglacial settings. Here, we analyse the role of earthquake-induced rock mass damage over late Pleistocene to Holocene time scales, as a conditioning factor for modern landslides, using distinct element numerical modelling to assess spatial and temporal patterns of fracture propagation influenced by varying glacier thickness. Conceptualised numerical models were parameterised by in situ rock mass, glacial and topographic conditions in Serpentine valley, located in Prince William Sound, Alaska, where several large landslides are actively developing along the western valley wall. Results show that although rigid glacier buttressing reduces co-seismic rock mass damage, it does not suppress it completely, and damage occurs both above and below the glacier surface elevation. We demonstrate that topography, and especially steep slopes with topographic convexity, as well as preexisting damage of joint networks and faults inherited from tectonic and exhumation induced stresses, exert primary control on the location of new co-seismic damage. Simulations representing a simplified deglaciation sequence over the past 25 ka, with a series of 10 evenly spaced earthquakes, generated rock mass damage patterns that qualitatively match in situ landslide structural and kinematic observations at one instability in Serpentine valley. Our conceptual study helps clarify the role of repeated seismicity over glacial timescales as a long-term conditioning process for paraglacial rock slope failure, highlighting spatial and temporal patterns of progressive damage accumulation, with outcomes relevant for modern landslide hazard assessment.

Alaska↗