An official website of the United States government
Abstract The advance of a chemical weathering front into the bedrock of a hillslope is often limited by the rate weathering products that can be carried away, maintaining chemical disequilibrium. If the weathering front is within the saturated zone, groundwater flow downslope may affect the rate of transport and weathering—however, weathering also modifies the rock permeability and the subsurface potential gradient that drives lateral groundwater flow. This feedback may help explain why there tends to be neither “runaway weathering” to great depth nor exposed bedrock covering much of the earth and may provide a mechanism for weathering front advance to keep pace with incision of adjacent streams into bedrock. This is the second of a two‐part paper exploring the coevolution of bedrock weathering and lateral flow in hillslopes using a simple low‐dimensional model based on hydraulic groundwater theory. Here, we show how a simplified kinetic model of 1‐D rock weathering can be extended to consider lateral flow in a 2‐D hillslope. Exact and approximate analytical solutions for the location and thickness of weathering within the hillslope are obtained for a number of cases. A location for the weathering front can be found such that lateral flow is able to export weathering products at the rate required to keep pace with stream incision at steady state. Three pathways of solute export are identified: “diffusing up,” where solutes diffuse up and away from the weathering front into the laterally flowing aquifer; “draining down,” where solutes are advected primarily downward into the unweathered bedrock; and “draining along,” where solutes travel laterally within the weathering zone. For each pathway, a different subsurface topography and overall relief of unweathered bedrock within the hillslope is needed to remove solutes at steady state. The relief each pathway requires depends on the rate of stream incision raised to a different power, such that at a given incision rate, one pathway requires minimal relief and, therefore, likely determines the steady‐state hillslope profile.
more »
« less
Modeling the Shape and Evolution of Normal‐Fault Facets
Facets formed along the footwalls of active normal‐fault blocks display a variety of longitudinal profile forms, with variations in gradient, shape, degree of soil cover, and presence or absence of a slope break at the fault trace. We show that a two‐dimensional, process‐oriented cellular automaton model of facet profile evolution can account for the observed morphologic diversity. The model uses two dimensionless parameters to represent fault slip, progressive rock weathering, and downslope colluvial‐soil transport driven by gravity and stochastic disturbance events. The parameters represent rock weathering and soil disturbance rates, respectively, scaled by fault slip rate; both can be derived from field‐estimated rate coefficients. In the model's transport‐limited regime, slope gradient depends on the ratio of disturbance to slip rate, with a maximum that represents the angle of repose for colluvium. In this regime, facet evolution is consistent with nonlinear diffusion models of soil‐mantled hillslope evolution. Under the weathering‐limited regime, bedrock becomes partly exposed but microtopography helps trap some colluvium even when facet gradient exceeds the threshold angle. Whereas the model predicts a continuous gradient from footwall to colluvial wedge under transport‐limited behavior, fully weathering‐limited facets tend to develop a slope break between footwall and basal colluvium as a result of reduced transport efficiency on the rocky footwall slope. To the extent that the model provides a reasonable analogy for natural facets, its behavior suggests that facet profile morphology can provide useful constraints on relative potential rates of rock weathering, soil disturbance, and fault slip.
more »
« less
- Award ID(s):
- 1831623
- PAR ID:
- 10458110
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Earth Surface
- Volume:
- 125
- Issue:
- 3
- ISSN:
- 2169-9003
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract It is widely recognized that fluid injection can trigger aseismic fault slip. However, the processes by which the fluid‐rock interactions facilitate or inhibit slip are poorly understood and some are oversimplified in most models of injection‐induced slip. In this study, we perform a 2D anti‐plane shear investigation of aseismic slip that occurs in response to fluid injection into a permeable fault governed by rate‐and‐state friction. We account for porosity and permeability changes that accompany slip, including dilatancy, and quantify how these processes affect pore pressure diffusion, which couples to aseismic slip. Fault response to injection has two phases. In the first phase, slip is negligible and pore pressure closely follows the standard linear diffusion model. Pressurization eventually triggers aseismic slip close to the injection site. In the second phase, aseismic slip front expands outward and dilatancy causes pore pressure to depart from the linear diffusion model. We quantify how prestress, injection rate, permeability and other fluid transport properties affect the slip front migration rate, finding rates ranging from 10 to 1,000 m/day for typical parameters. The migration rate is strongly influenced by the fault's closeness to failure and injection rate. The total slip on the fault, on the other hand, is primarily determined by the injected volume, with minimal sensitivity to injection rate. Additionally, we show that when dilatancy is neglected, slip front migration rate and total slip can be several times higher. Our modeling demonstrates that porosity and permeability evolution, especially dilatancy, fundamentally alters how faults respond to fluid injection.more » « less
-
In the Earth sciences, weathering encompasses all the physical, chemical, and biological processes that break down rocks in place. Rock weathering takes decades to millions of years and impacts climate and soil formation. In our two-part lesson, students develop an understanding of weathering and how it can influence climate and human society through hands-on experiments. Lesson 1 focuses on how rock weathering impacts climate; students investigate how changing the temperature and acidity of weathering agents affects the rate of rock weathering. Lesson 2 focuses on how weathering impacts human society; students perform experiments simulating weathering of mudstone and granite via shaking rocks in containers; students observe that these rocks weather at different rates and produce different-sized particles because of physical weathering. Students relate their experimental observations to the process of soil formation and then apply this knowledge to societally relevant topics. These lessons bring rock weathering into the classroom with crosscutting concepts and connect the Earth, climate, and human society together in an interactive way.more » « less
-
The geometry and magnitude of finite strain in the ductile footwalls of metamorphic core complexes are important parameters for testing the predictions of models of extension, yet are often difficult to quantify due to the rare preservation of deformed markers. The footwall of the Northern Snake Range core complex in eastern Nevada preserves a coherent stratigraphy of ductilely thinned Neoproterozoic-Cambrian metasedimentary rocks that are exposed over a 30 km transport-parallel distance, and thus provides an important opportunity to quantify footwall strain. We measured strain ellipsoids from stretched detrital quartz grains and ribbons in 45 samples that span the full exposed distance of the ductilely sheared footwall of the master detachment fault (the Northern Snake Range décollement), and we combined our data with 11 published strain ellipsoids. On the eastern side of the range, where recrystallization limits the preservation of detrital quartz grains, we estimated finite strain by comparing the attenuated thicknesses of Neoproterozoic-Cambrian rock units to their regional stratigraphic thicknesses. Our data demonstrate a dramatic gradient in ductile strain in the transport direction, from 39% subhorizontal extension and 32% subvertical thinning at the western flank of the Northern Snake Range to 450-1440% extension and 81-94% thinning at the eastern flank. The footwall underwent 18-20 km of cumulative ductile extension, which is equivalent to 38-43% of the 47 km of total extension accommodated on brittle structures. Kinematic vorticity estimates from published quartz petrofabrics define an eastward-increasing component of top-to-the-ESE simple shear. Our data are compatible with a rolling hinge model of extension, where displacement on a low-angle, upper-crustal, brittle detachment fault system was fed downward to a zone of distributed, simple shear-dominant, top-down-to-ESE ductile shearing beneath the quartz crystal-plastic transition. The progressive eastward translation and brittle thinning of the hanging wall resulted in the eastward migration of exhumation of footwall rocks. Migrating exhumation may in part be responsible for the eastward-increasing finite strain gradient, as footwall rocks in the eastern part of the range experienced a longer strain history.more » « less
-
Landscape disturbance events (e.g., earthquakes, slope failures) play key roles in landscape evolution in tectonically active areas. Along the Teton fault, fault scarps vary in height by up to tens of meters. LiDAR-based mapping indicates that scarp height is affected by glacial geomorphology, slope failure, and alluvial processes. LiDAR data, digital and field mapping were used to characterize fault scarps and slope failure deposits along the Teton fault zone. Based on vertical separation (VS; the vertical offset between faulted surfaces) across fault scarps and the expected behavior of normal faults, we propose a four-section model of the Teton fault. At a broad scale, VS is greatest along the southern fault zone. At a finer scale, VS is least at the ends of the fault and at three areas within the central fault zone. Transitions between these four sections may represent segment boundaries with potentially important implications for geohazards assessment.more » « less
