Geology Reports⌕ Search

Geology topics

William E. Dietrich

Publications and source records attributed to William E. Dietrich.

13 recordsLinked to original sources
USGS2026

Deep critical zone controls on shallow landslides

The deep critical zone (CZ) has long been recognized for its importance in influencing shallow landslides but was not considered feasible to include in slope stability models at the watershed scale. In this study, we demonstrate that simple approximations of the CZ in a fully coupled hydrologic and soil slope stability model can effectively capture the location, timing, and likely size of shallow landslides. To achieve this, we use coupled, process-based models that incorporate the effects of 1) deep CZ structures, 2) three-dimensional transient hydrology, and 3) multidimensional slope stability, calibrated with data from an intensively monitored field site. Our results show that the hydrologically active deep CZ guides groundwater flow, influencing where it drains from or exfiltrates to the soil mantle, producing distinct patterns of soil saturation and seepage forces at the soil-bedrock boundary. Deep conductive weathered critical zone drains the soil mantle, reducing the likelihood of destabilizing pore pressures, while the downslope thinning of the CZ forces groundwater to the surface. This creates localized instability and a tendency for similar-sized landslides across the landscape. In contrast, the absence of conductive weathered bedrock results in more widespread destabilizing pore pressures, leading to larger landslides and the likelihood of landslides earlier in a storm than in landscapes underlain by a deep CZ. Our findings suggest that first-order variations of deep CZ can provide physical explanations for variations observed in the susceptibility, magnitude, and timing of shallow landslides, and that CZ structure may be inferred from patterns and timing of landsliding.

Proceedings of the National Academy of Sciences↗
USGS2024

Geological context and significance of the clay-sulfate transition region in Mount Sharp, Gale crater, Mars: An integrated assessment based on orbiter and rover data

On Mars, phyllosilicate (“clay”) minerals are often associated with older terrains, and sulfate minerals are associated with younger terrains, and this dichotomy is taken as evidence that Mars’ surface dried up over time. Therefore, in situ investigation of the Mount Sharp strata in Gale crater, which record a shift from dominantly clay-bearing to sulfate-bearing minerals, as seen in visible−near-infrared orbital reflectance spectra, is a key science objective for the Mars Science Laboratory (MSL) Curiosity rover mission. Here, we present regional (orbiter-based) and in situ (rover-based) evidence for a low-angle erosional unconformity that separates the lacustrine and marginal lacustrine deposits of the Carolyn Shoemaker formation from the dominantly eolian deposits of the lower Mirador formation within the orbitally defined clay-sulfate transition region. The up-section record of wetter (Carolyn Shoemaker formation) to drier (lower Mirador formation) depositional conditions is accompanied by distinct changes in diagenesis. Clay minerals occur preferentially within the Carolyn Shoemaker formation and are absent within the lower members of the Mirador formation. At and above the proposed unconformity, strata are characterized by an increase in diagenetic nodules enriched in X-ray amorphous Mg-sulfate. Early clay formation in the Carolyn Shoemaker formation may have created a hydraulic barrier such that later migrating magnesium- and sulfur-rich fluids accumulated preferentially within the lower members of the Mirador formation. The proposed unconformity may have also acted as a fluid conduit to further promote Mg-sulfate nodule formation at the Carolyn Shoemaker−Mirador formation boundary. These results confirm an association of the clay-sulfate transition with the drying of depositional environments, but they also suggest that at least some orbital sulfate signatures within the region are not time-congruent with the environmental signals extracted from primary sedimentology. Our findings highlight that complex interactions among primary depositional environment, erosion, and diagenesis contribute to the transition in clay-sulfate orbital signatures observed in the stratigraphy of Mount Sharp.

GSA Bulletin↗
USGS2022

The Curiosity Rover’s exploration of Glen Torridon, Gale crater, Mars: An overview of the campaign and scientific results

The Mars Science Laboratory rover, Curiosity , explored the clay mineral-bearing Glen Torridon region for one martian year between January 2019 and January 2021, including a short campaign onto the Greenheugh pediment. The Glen Torridon campaign sought to characterize the geology of the area, seek evidence of habitable environments, and document the onset of a potentially global climatic transition during the Hesperian era. Curiosity roved 5 km in total throughout Glen Torridon, from the Vera Rubin ridge to the northern margin of the Greenheugh pediment. Curiosity acquired samples from 11 drill holes during this campaign and conducted the first martian thermochemolytic-based organics detection experiment with the Sample Analysis at Mars instrument suite. The lowest elevations within Glen Torridon represent a continuation of lacustrine Murray formation deposits, but overlying widespread cross bedded sandstones indicate an interval of more energetic fluvial environments and prompted the definition of a new stratigraphic formation in the Mount Sharp group called the Carolyn Shoemaker formation. Glen Torridon hosts abundant phyllosilicates yet remains compositionally and mineralogically comparable to the rest of the Mount Sharp group. Glen Torridon samples have a great diversity and abundance of sulfur-bearing organic molecules, which are consistent with the presence of ancient refractory organic matter. The Glen Torridon region experienced heterogeneous diagenesis, with the most striking alteration occurring just below the Siccar Point unconformity at the Greenheugh pediment. Results from the pediment campaign show that the capping sandstone formed within the Stimson Hesperian aeolian sand sea that experienced seasonal variations in wind direction.

JGR Planets↗
USGS2022

Geology and stratigraphic correlation of the Murray and Carolyn Shoemaker formations across the Glen Torridon region, Gale crater, Mars

The Glen Torridon (GT) region within Gale crater, Mars, occurs in contact with the southern side of Vera Rubin ridge (VRR), a well-defined geomorphic feature that is comparatively resistant to erosion. Prior to detailed ground-based investigation of GT, its geologic relationship with VRR was unknown. Distinct lithologic subunits within the Jura member (Murray formation), which forms the upper part of VRR, made it possible to be also identified within GT. This indicates that the strata pass across the geomorphic divide between regions. Furthermore, the cross-bedded lower part of the overlying Knockfarril Hill member (Carolyn Shoemaker formation) also occurs within both VRR and GT. Correlation of both units demonstrates that the strata form a continuous stratigraphic succession regardless of large-scale geomorphic expression. The lithologic change from mudstone (Jura member) to cross-bedded sandstone (Knockfarril Hill member) heralds a significant shift in paleoenvironment from lacustrine to fluvial. The upper part of the Knockfarril Hill member consists of interbedded mudstone and sandstone that transitions to the overlying finely laminated mudstone of the Glasgow member, and a return to lacustrine deposition. In GT, the Stimson formation unconformably overlies the Glasgow member, where it demarks the southern boundary of GT. Contacts for each stratigraphic unit were defined and transferred to a high-resolution image base to make a geologic map and cross sections perpendicular to the NE strike. Stratal dips cannot exceed 2° NW to retain the positions of stratigraphic units in the locations they are exposed throughout GT.

Journal of Geophysical Research: Planets↗
USGS2020

Evidence for a diagenetic origin of Vera Rubin Ridge, Gale Crater, Mars: Summary and synthesis of Curiosity's exploration campaign

This paper provides an overview of the Curiosity rover's exploration at Vera Rubin ridge (VRR) and summarizes the science results. VRR is a distinct geomorphic feature on lower Aeolis Mons (informally known as Mount Sharp) that was identified in orbital data based on its distinct texture, topographic expression, and association with a hematite spectral signature. Curiosity conducted extensive remote sensing observations, acquired data on dozens of contact science targets, and drilled three outcrop samples from the ridge, as well as one outcrop sample immediately below the ridge. Our observations indicate that strata composing VRR were deposited in a predominantly lacustrine setting and are part of the Murray formation. The rocks within the ridge are chemically in family with underlying Murray formation strata. Red hematite is dispersed throughout much of the VRR bedrock, and this is the source of the orbital spectral detection. Gray hematite is also present in isolated, gray‐colored patches concentrated toward the upper elevations of VRR, and these gray patches also contain small, dark Fe‐rich nodules. We propose that VRR formed when diagenetic event(s) preferentially hardened rocks, which were subsequently eroded into a ridge by wind. Diagenesis also led to enhanced crystallization and/or cementation that deepened the ferric‐related spectral absorptions on the ridge, which helped make them readily distinguishable from orbit. Results add to existing evidence of protracted aqueous environments at Gale crater and give new insight into how diagenesis shaped Mars' rock record.

Journal of Geophysical Research: Planets↗
USGS2017

Expanding the role of reactive transport models in critical zone processes

Models test our understanding of processes and can reach beyond the spatial and temporal scales of measurements. Multi-component Reactive Transport Models (RTMs), initially developed more than three decades ago, have been used extensively to explore the interactions of geothermal, hydrologic, geochemical, and geobiological processes in subsurface systems. Driven by extensive data sets now available from intensive measurement efforts, there is a pressing need to couple RTMs with other community models to explore non-linear interactions among the atmosphere, hydrosphere, biosphere, and geosphere. Here we briefly review the history of RTM development, summarize the current state of RTM approaches, and identify new research directions, opportunities, and infrastructure needs to broaden the use of RTMs. In particular, we envision the expanded use of RTMs in advancing process understanding in the Critical Zone, the veneer of the Earth that extends from the top of vegetation to the bottom of groundwater. We argue that, although parsimonious models are essential at larger scales, process-based models offer tools to explore the highly nonlinear coupling that characterizes natural systems. We present seven testable hypotheses that emphasize the unique capabilities of process-based RTMs for (1) elucidating chemical weathering and its physical and biogeochemical drivers; (2) understanding the interactions among roots, micro-organisms, carbon, water, and minerals in the rhizosphere; (3) assessing the effects of heterogeneity across spatial and temporal scales; and (4) integrating the vast quantity of novel data, including “omics” data (genomics, transcriptomics, proteomics, metabolomics), elemental concentration and speciation data, and isotope data into our understanding of complex earth surface systems. With strong support from data-driven sciences, we are now in an exciting era where integration of RTM framework into other community models will facilitate process understanding across disciplines and across scales.

Earth-Science Reviews↗
USGS2015

Coevolution of bed surface patchiness and channel morphology: 1. Mechanisms of forced patch formation

Riverbeds frequently display a spatial structure where the sediment mixture composing the channel bed has been sorted into discrete patches of similar grain size. Even though patches are a fundamental feature in gravel bed rivers, we have little understanding of how patches form, evolve, and interact. Here we present a two-dimensional morphodynamic model that is used to examine in greater detail the mechanisms responsible for the development of forced bed surface patches and the coevolution of bed morphology and bed surface patchiness. The model computes the depth-averaged channel hydrodynamics, mixed-grain-size sediment transport, and bed evolution by coupling the river morphodynamic model Flow and Sediment Transport with Morphological Evolution of Channels (FaSTMECH) with a transport relation for gravel mixtures and the mixed-grain-size Exner equation using the active layer assumption. To test the model, we use it to simulate a flume experiment in which the bed developed a sequence of alternate bars and temporally and spatially persistent forced patches with a general pattern of coarse bar tops and fine pools. Cross-stream sediment flux causes sediment to be exported off of bars and imported into pools at a rate that balances downstream gradients in the streamwise sediment transport rate, allowing quasi-steady bar-pool topography to persist. The relative importance of lateral gravitational forces on the cross-stream component of sediment transport is a primary control on the amplitude of the bars. Because boundary shear stress declines as flow shoals over the bars, the lateral sediment transport is increasingly size selective and leads to the development of coarse bar tops and fine pools.

Journal of Geophysical Research F: Earth Surface↗
USGS2015

Coevolution of bed surface patchiness and channel morphology: 2. Numerical experiments

In gravel bed rivers, bed topography and the bed surface grain size distribution evolve simultaneously, but it is not clear how feedbacks between topography and grain sorting affect channel morphology. In this, the second of a pair of papers examining interactions between bed topography and bed surface sorting in gravel bed rivers, we use a two-dimensional morphodynamic model to perform numerical experiments designed to explore the coevolution of both free and forced bars and bed surface patches. Model runs were carried out on a computational grid simulating a 200 m long, 2.75 m wide, straight, rectangular channel, with an initially flat bed at a slope of 0.0137. Over five numerical experiments, we varied (a) whether an obstruction was present, (b) whether the sediment was a gravel mixture or a single size, and (c) whether the bed surface grain size feeds back on the hydraulic roughness field. Experiments with channel obstructions developed a train of alternate bars that became stationary and were connected to the obstruction. Freely migrating alternate bars formed in the experiments without channel obstructions. Simulations incorporating roughness feedbacks between the bed surface and flow field produced flatter, broader, and longer bars than simulations using constant roughness or uniform sediment. Our findings suggest that patches are not simply a by-product of bed topography, but they interact with the evolving bed and influence morphologic evolution.

Journal of Geophysical Research F: Earth Surface↗
USGS2015

The Landscape Evolution Observatory: a large-scale controllable infrastructure to study coupled Earth-surface processes

Zero-order drainage basins, and their constituent hillslopes, are the fundamental geomorphic unit comprising much of Earth's uplands. The convergent topography of these landscapes generates spatially variable substrate and moisture content, facilitating biological diversity and influencing how the landscape filters precipitation and sequesters atmospheric carbon dioxide. In light of these significant ecosystem services, refining our understanding of how these functions are affected by landscape evolution, weather variability, and long-term climate change is imperative. In this paper we introduce the Landscape Evolution Observatory (LEO): a large-scale controllable infrastructure consisting of three replicated artificial landscapes (each 330 m 2 surface area) within the climate-controlled Biosphere 2 facility in Arizona, USA. At LEO, experimental manipulation of rainfall, air temperature, relative humidity, and wind speed are possible at unprecedented scale. The Landscape Evolution Observatory was designed as a community resource to advance understanding of how topography, physical and chemical properties of soil, and biological communities coevolve, and how this coevolution affects water, carbon, and energy cycles at multiple spatial scales. With well-defined boundary conditions and an extensive network of sensors and samplers, LEO enables an iterative scientific approach that includes numerical model development and virtual experimentation, physical experimentation, data analysis, and model refinement. We plan to engage the broader scientific community through public dissemination of data from LEO, collaborative experimental design, and community-based model development.

Arizona↗
USGS2011

Spatially explicit shallow landslide susceptibility mapping over large areas

Recent advances in downscaling climate model precipitation predictions now yield spatially explicit patterns of rainfall that could be used to estimate shallow landslide susceptibility over large areas. In California, the United States Geological Survey is exploring community emergency response to the possible effects of a very large simulated storm event and to do so it has generated downscaled precipitation maps for the storm. To predict the corresponding pattern of shallow landslide susceptibility across the state, we have used the model Shalstab (a coupled steady state runoff and infinite slope stability model) which susceptibility spatially explicit estimates of relative potential instability. Such slope stability models that include the effects of subsurface runoff on potentially destabilizing pore pressure evolution require water routing and hence the definition of upslope drainage area to each potential cell. To calculate drainage area efficiently over a large area we developed a parallel framework to scale-up Shalstab and specifically introduce a new efficient parallel drainage area algorithm which produces seamless results. The single seamless shallow landslide susceptibility map for all of California was accomplished in a short run time, and indicates that much larger areas can be efficiently modelled. As landslide maps generally over predict the extent of instability for any given storm. Local empirical data on the fraction of predicted unstable cells that failed for observed rainfall intensity can be used to specify the likely extent of hazard for a given storm. This suggests that campaigns to collect local precipitation data and detailed shallow landslide location maps after major storms could be used to calibrate models and improve their use in hazard assessment for individual storms.

California↗
USGS2005

Field measurements of incision rates following bedrock exposure: Implications for process controls on the long profiles of valleys cut by rivers and debris flows

Until recently, published rates of incision of bedrock valleys came from indirect dating of incised surfaces. A small but growing literature based on direct measurement reports short-term bedrock lowering at geologically unsustainable rates. We report observations of bedrock lowering from erosion pins monitored over 1–7 yr in 10 valleys that cut indurated volcanic and sedimentary rocks in Washington, Oregon, California, and Taiwan. Most of these channels have historically been stripped of sediment. Their bedrock is exposed to bed-load abrasion, plucking, and seasonal wetting and drying that comminutes hard, intact rock into plates or equant fragments that are removed by higher flows. Consequent incision rates are proportional to the square of rock tensile strength, in agreement with experimental results of others. Measured rates up to centimeters per year far exceed regional long-term erosion-rate estimates, even for apparently minor sediment-transport rates. Cultural artifacts on adjoining strath terraces in Washington and Taiwan indicate at least several decades of lowering at these extreme rates. Lacking sediment cover, lithologies at these sites lower at rates that far exceed long-term rock-uplift rates. This rate disparity makes it unlikely that the long profiles of these rivers are directly adjusted to either bedrock hardness or rock-uplift rate in the manner predicted by the stream power law, despite the observation that their profiles are well fit by power-law plots of drainage area vs. slope. We hypothesize that the threshold of motion of a thin sediment mantle, rather than bedrock hardness or rock-uplift rate, controls channel slope in weak bedrock lithologies with tensile strengths below ∼3–5 MPa. To illustrate this hypothesis and to provide an alternative interpretation for power-law plots of area vs. slope, we combine Shields' threshold transport concept with measured hydraulic relationships and downstream fining rates. In contrast to fluvial reaches, none of the hundreds of erosion pins we installed in steep valleys recently scoured to bedrock by debris flows indicate any postevent fluvial lowering. These results are consistent with episodic debris flows as the primary agent of bedrock lowering in the steepest parts of the channel network above ∼0.03–0.10 slope.

California, Washington, Oregon↗