Comparison of soil thickness in a zero-order basin in the Oregon Coast Range using a soil probe and electrical resistivity tomography

Michael S. Morse; Ning Lu; Jonathan W. Godt; André Revil; Jeffrey A. Coe

doi:10.1061/(asce)gt.1943-5606.0000717

Home/Archive/USGS

USGS · 70045123

Comparison of soil thickness in a zero-order basin in the Oregon Coast Range using a soil probe and electrical resistivity tomography

Michael S. Morse·Ning Lu·Jonathan W. Godt·André Revil·Jeffrey A. Coe

View original source ↗Export metadata ↓

Abstract

Accurate estimation of the soil thickness distribution in steepland drainage basins is essential for understanding ecosystem and subsurface response to infiltration. One important aspect of this characterization is assessing the heavy and antecedent rainfall conditions that lead to shallow landsliding. In this paper, we investigate the direct current (DC) resistivity method as a tool for quickly estimating soil thickness over a steep (33–40°) zero-order basin in the Oregon Coast Range, a landslide prone region. Point measurements throughout the basin showed bedrock depths between 0.55 and 3.2 m. Resistivity of soil and bedrock samples collected from the site was measured for degrees of saturation between 40 and 92%. Resistivity of the soil was typically higher than that of the bedrock for degrees of saturation lower than 70%. Results from the laboratory measurements and point-depth measurements were used in a numerical model to evaluate the resistivity contrast at the soil-bedrock interface. A decreasing-with-depth resistivity contrast was apparent at the interface in the modeling results. At the field site, three transects were surveyed where coincident ground truth measurements of bedrock depth were available, to test the accuracy of the method. The same decreasing-with-depth resistivity trend that was apparent in the model was also present in the survey data. The resistivity contour of between 1,000 and 2,000 Ωm that marked the top of the contrast was our interpreted bedrock depth in the survey data. Kriged depth-to-bedrock maps were created from both the field-measured ground truth obtained with a soil probe and interpreted depths from the resistivity tomography, and these were compared for accuracy graphically. Depths were interpolated as far as 16.5 m laterally from the resistivity survey lines with root mean squared error (RMSE) = 27 cm between the measured and interpreted depth at those locations. Using several transects and analysis of the subsurface material properties, the direct current (DC) resistivity method is shown to be able to delineate bedrock depth trends within the drainage basin.

Explore related subjects

Journal of Geotechnical and Geoenvironmental Engineering Article

Source-reported bounding extent: 42° to 46.29° latitude; -124.61° to -116.46° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

Michael S. Morse, Ning Lu, Jonathan W. Godt, André Revil, Jeffrey A. Coe. 2012. Comparison of soil thickness in a zero-order basin in the Oregon Coast Range using a soil probe and electrical resistivity tomography. https://doi.org/10.1061/(asce)gt.1943-5606.0000717

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

KEEP EXPLORING

Related USGS reports

USGS2026

Discrete element investigation of the influence of shallow soil density on the manifestations of strike-slip surface fault rupture

This study investigates the influence of soil relative density on strike-slip surface fault rupture manifestations using three-dimensional numerical simulations performed with the discrete element method (DEM). The simulations capture the formation of distinctive fault strands within complex flower structures using tens of millions of grains. The tendency for dense soils to localize shear manifests as multiple localized shear bands within wide zones of deformation in strike-slip fault rupture, whereas diffuse shear deformation develops within narrow zones in loose soils. The spatial extents of soil deformation are consistent between simulations having similar relative densities but with different quantities of grains. However, individual shears are more distinguishable in assemblages of finer grains than in assemblages of coarser grains. The simulations show the progressive development of new shears within the bounds of previously developed shears. Shear activity transitions inward as fault activity diminishes along the outermost shears and continues along newly developed shears until a vertically dipping throughgoing shear structure develops that accommodates most of the fault displacement thereon. The throughgoing fault develops at smaller fault displacements in looser soils because the first shear rupture propagates closer to the vertical direction and does not undergo as much inward translation of shear activity, as is observed in denser soils. In all simulations, ground surface uplift develops between nonintersecting active shears, and ground surface subsidence tends to develop where new shears intersect previous shears. The surface traces in these simulations are shown to be consistent with analog models and case histories of surface fault rupture occurring in different shallow subsurface materials. Although computationally costly, these modeling results are valuable for providing a strong numerical supplement to traditional analog models used to represent the mechanics of strike-slip zones in soil, and they provide quantifiable stresses and large-strain deformations throughout the model domain.

Journal of Geotechnical and Geoenvironmental Engin↗

USGS2024

Framework for mapping liquefaction hazard–Targeted design ground motions

Liquefaction-induced ground failure poses substantial challenges to geotechnical earthquake engineering design. Current approaches for designing against liquefaction hazards, as specified in most seismic provisions, focus on estimating a liquefaction factor of safety ( 𝐹⁢𝑆𝐿 ) and typically characterize earthquake loading using design parameters based on probabilistic or deterministic ground motion levels. Because 𝐹⁢𝑆𝐿 is estimated deterministically, this basis of design neglects considerable uncertainties for estimating liquefaction triggering and its consequences and results in a lack of liquefaction-specific design criteria, particularly as structural design has advanced toward risk-targeted performance objectives. This study presents a framework for developing liquefaction-targeted design criteria based on a minimum acceptable return period of liquefaction, informed by probabilistic liquefaction hazard analysis (PLHA). PLHA quantifies annualized rates of liquefaction by considering contributions from (1) the full ground-motion probability space, and (2) uncertainties in liquefaction triggering using probabilistic models. PLHA is used in this study to characterize the current, effective return periods of 𝐹⁢𝑆𝐿 ( 𝑇𝑅,𝐹⁢𝑆 ) obtained from conventional liquefaction hazard analysis (CLHA) using uniform-hazard ground motions. 𝑇𝑅,𝐹⁢𝑆 is evaluated in a parametric study of nearly 100 sites throughout the conterminous United States. The results indicate large geographic variations in acceptable liquefaction hazard levels, with implied 𝑇𝑅,𝐹⁢𝑆 ranging between approximately 1,000 to 3,000 years. To address these inconsistencies without the computational demands of full PLHA, a framework is proposed for developing a liquefaction-targeted design peak ground acceleration, 𝑃⁢𝐺⁢𝐴𝐿 , for use in liquefaction models that result in consistent liquefaction design levels across all geographic locations. The mapped 𝑃⁢𝐺⁢𝐴𝐿 is shown to be somewhat sensitive to site-specific properties, and adjustment factors are developed and presented. The proposed 𝑃⁢𝐺⁢𝐴𝐿 mapping procedure produces 𝐹⁢𝑆𝐿 estimates that are consistent with those obtained from full PLHA at a target 𝑇𝑅,𝐹⁢𝑆 , providing a promising roadmap to incorporating PLHA concepts into current liquefaction design methods.

Journal of Geotechnical and Geoenvironmental Engin↗

USGS2019

Effects of infiltration characteristics on the spatial-temporal evolution of stability of an interstate highway embankment

Infiltration-induced landslides are among the most common natural disasters threatening modern civilization, but conventional methods for studying the triggering mechanisms and predicting the occurrence of these slides are limited by incomplete consideration of underlying physical processes and the lack of precision inherent in limit-equilibrium analyses. To address this problem the spatial-temporal evolution of failure is investigated in a seasonally unstable section of interstate highway embankment, known as the Straight Creek landslide, Colorado. The study includes multi-year site investigation, monitoring, and numerical simulation using a rigorous hydromechanical framework along with a field of local factor of safety method. The sensitivity of episodic landslide reactivation to infiltration characteristics is evaluated. Results indicate that annual cumulative snowmelt infiltration, which typically accounts for approximately 75% of total annual cumulative infiltration and occurs over a short period in the spring, has the most substantial impact on slide activation. The rate of snowmelt infiltration varies independently of annual cumulative snowmelt infiltration and cumulative infiltration in the previous year, but still affects antecedent soil moisture conditions at the onset of snowmelt infiltration and therefore also the level of slide activation. These findings are used to establish specific thresholds for exacerbated slide movement using annual snowpack accumulation, forecasted snowmelt rate, and the previous year’s snowmelt, an approach which may be applied for predicting movement at this and other recurring or potential slide sites.

Colorado↗