E. Kirby

Aspect-dependent variations in regolith creep revealed by meteoric 10Be

Although variations in insolation and emergent feedbacks among soil moisture, vegetation, and soil cohesion are commonly invoked to explain topographic asymmetry that depends on aspect, few studies have directly quantified the efficiency of regolith transport along hillslopes of opposing aspect. We utilize meteoric 10 Be concentrations in regolith (n = 74) to determine mass flux along equatorial-facing and polar-facing hillslopes in three forested upland watersheds in and adjacent to the Susquehanna Shale Hills Critical Zone Observatory in central Pennsylvania (USA). In combination with regolith depth measurements and high-resolution topography, these fluxes allow us to evaluate transport rate laws and the efficiency of regolith creep. Concentrations of meteoric 10 Be in regolith along six separate transects imply that regolith flux is similar along all hillslopes, despite differences in topographic gradient and regolith thickness. Comparison of flux with regolith depth and topographic gradient reveals that transport depends on regolith depth, and that regolith creep is twice as efficient along low-gradient, south-facing slopes with thin regolith as compared to steep, north-facing slopes mantled with thicker regolith. We suggest that the observed topographic asymmetry in these watersheds has evolved over geologic time as a result of differences in the frequency of freeze-thaw events between hillslopes of opposing aspect.

Oxidative dissolution under the channel leads geomorphological evolution at the Shale Hills catchment

The hydrologic connectivity between hillslopes and streams impacts the geomorphological evolution of catchments. Here, we propose a conceptual model for hydrogeomorphological evolution of the Susquehanna Shale Hills Critical Zone Observatory (SSHCZO), a first-order catchment developed on shale in central Pennsylvania, U.S.A. At SSHCZO, the majority of available water (the difference between incoming meteoric water and outgoing evapotranspiration) flows laterally to the catchment outlet as interflow, while the rest is transported by regional groundwater flow. Interflow, shallow hillslope flow, is limited to the upper 5 to 8 m of highly fractured bedrock, thought to have formed during periglacial conditions in the late Pleistocene. In contrast, groundwater flowpaths are influenced by the primary permeability of the varying stratigraphic units. Both flowpaths respond to weathering-related secondary permeability. O2-rich interflow mixes with deep O2-poor groundwater under the catchment outlet at depths of 5–8 m. Penetration of this oxygenated interflow under the valley results in pyrite oxidation, release of sulfuric acid, dissolution of minerals, and weakening of bedrock. This is hypothesized to enhance channel incision and, in turn, to promote drainage of deep groundwater from the ridges. Drainage subsequently lowers the catchment water table, advancing the cascade of reactions that produce regolith. Weathering in the catchment is characterized by both sharp and diffuse reaction fronts. Relatively sharp fronts (pyrite, carbonate) mark where vertical, unsaturated flow changes to horizontal, saturated flow, while diffuse fronts (illite, chlorite, feldspar) mark where flow is largely vertical and unsaturated. According to this model, catchment morphology reflects subsurface pyrite reactions due to mixing of interflow and groundwater flow under the valley floor that ultimately results in clay weathering and regolith production nearer the land surface.

Pleistocene drainage reorganization driven by the isostatic response to deep incision into the northeastern Tibetan Plateau

Pleistocene drainage basin integration led to progressive excavation of Tertiary-Quaternary sedimentary basins along the Yellow River in the northeastern Tibetan Plateau. Cosmogenic burial dating of ancestral river deposits and basin fi ll from two key watershed divides confi rms a fl uvial connection between basins at 0.5-1.2 Ma, prior to excavation by the Yellow River. Pres- ervation of the relict depositional surface that represents the maximum height of basin fi ll allows reconstruction of the volume of eroded material across a broad region. We quantify the iso- static response to this erosional unloading using a two-dimensional fl exural model. Calculated maximum vertical displacements for different effective elastic thicknesses vary from ~160 m to ~260 m near the Pleistocene spillway from the Qinghai paleo-lake. We suggest that the isostatic response to fl excavation along the Yellow River defeated local tributaries, isolated Lake Qinghai, and led to the development of an internally drained basin in the past 0.5-1.2 m.y.

Rates and kinematics of active shortening along the eastern Qilian Shan, China, inferred from deformed fluvial terraces

In the eastern Qilian Shan, a flight of fluvial terraces developed along the Jinta River valley are deformed across the Nanying anticline. Four individual fluvial terraces are preserved at different elevations above the river, and higher terrace treads are draped by systematically thicker aeolian loess. Optically stimulated luminescence dating of deposits at the base of the loess provides constraints on the timing of surface abandonment; terraces were abandoned at 69 +/- 4 ka B.P. (T4), 57 +/- 4 ka B.P. (T3), and between 34 +/- 3 ka B.P. (T2), respectively. Differential GPS measurement of the terrace profile across the anticline allows reconstruction of subsurface fault geometry; we model terrace deformation above a listric thrust fault with a tip line at 2.2 +/- 0.1 km depth and whose dip shallows systematically to 23 +/- 3 degrees at depth of 5.8 +/- 1.1 km. Combining terrace ages with this model of fault geometry, we estimate a shortening rate of 0.8 +/- 0.2 mm/a across the Nanying fold and a shortening rate of similar to 0.1 mm/a across the mountain front fault since similar to 70 ka B.P. This rate suggests that the frontal fault system along the eastern Qilian Shan accomplishes crustal shortening at rates of approximately 0.9 +/- 0.3 mm/a during late Pleistocene time.

Holocene slip rate along the Gyaring Co Fault, central Tibet

Although geodetic measurements of interseismic deformation in interior Tibet suggest slow strain accumulation, active slip along the right-lateral Gyaring Co Fault is suggested to be between 8 and 21 mm/yr. Reliable geologic constraints on the slip rate along this fault are sparse. Here we document 12 ± 2 m of right-lateral displacement of lacustrine shorelines across the Gyaring Co Fault. Optically stimulated luminescence ages of the shorelines are tightly clustered between 4.1 and 4.4 ka. These data require an average slip rate of 2.2–3.2 mm/yr along the central Gyaring Co Fault during the latter half of the Holocene. Consideration of seismic cycle effects allows the possibility of slightly higher average slip rates, up to 2.2–4.5 mm/yr. Overall, our results suggest that the slip rate along the Gyaring Co Fault is similar to other strike-slip faults in interior Tibet, supporting the notion that active deformation in this region is distributed among numerous, slowly moving faults.

Late Miocene erosion and evolution of topography along the western slope of the Colorado Rockies

In the Colorado Rocky Mountains, the association of high topography and low seismic velocity in the underlying mantle suggests that recent changes in lithospheric buoyancy may have been associated with surface uplift of the range. This paper examines the relationships among late Cenozoic fluvial incision, channel steepness, and mantle velocity domains along the western slope of the northern Colorado Rockies. New 40 Ar/ 39 Ar ages on basalts capping the Tertiary Browns Park Formation range from ca. 11 to 6 Ma and provide markers from which we reconstruct incision along the White, Yampa, and Little Snake rivers. The magnitude of post–10 Ma incision varies systematically from north to south, increasing from ∼500 m along the Little Snake River to ∼1500 m along the Colorado River. Spatial variations in the amount of late Cenozoic incision are matched by metrics of channel steepness; the upper Colorado River and its tributaries (e.g., Gunnison and Dolores rivers) are two to three times steeper than the Yampa and White rivers, and these variations are independent of both discharge and lithologic substrate. The coincidence of steep river profiles with deep incision suggests that the fluvial systems are dynamically adjusting to an external forcing but is not readily explained by a putative increase in erosivity associated with late Cenozoic climate change. Rather, channel steepness correlates with the position of the channels relative to low-velocity mantle. We suggest that the history of late Miocene–present incision and channel adjustment reflects long-wavelength tilting across the western slope of the Rocky Mountains.

Kinematics of late Quaternary slip along the Yabrai fault: Implications for Cenozoic tectonics across the Gobi Alashan block, China

The Yabrai range-front fault accommodates deformation within the middle Gobi Alashan block between the Tibetan Plateau and the Ordos block. As such, it provides the opportunity to examine the transition between contractional deformation associated with the growth of the Tibetan Plateau and extensional deformation across North China. Geomorphic mapping of the active fault trace and trench investigations reveal that the Yabrai range-front fault is composed of three segments of varying fault strike, but for which the sense of motion, scarp height, and slip history appear to be kinematically compatible along the fault. Displaced Holocene and late Pleistocene alluvial deposits indicate that the southwestern segment is characterized by oblique-normal displacement with a minor sinistral component, whereas the middle segment appears to exhibit nearly dip-slip normal displacement. In contrast, slip along the northeastern segment appears to be primarily sinistral strike-slip with a minor reverse component. Geomorphically fresh fault scarps are developed within late Pleistocene–Holocene alluvial fans and terraces along the southwestern and northeastern segments, whereas the middle segment of the fault defines the bedrock-alluvial contact along the range front. The 10Be exposure ages of displaced alluvial fans along the southwestern segment yield a throw rate of ∼0.1 mm/yr over late Pleistocene time. Lateral slip rates along the northeastern fault segment range between 0.23 ± 0.02 and 0.78 ± 0.12 mm/yr. Regionally, the orientation and sense of motion along the Yabrai range-front fault are consistent with NE-SW shortening, and we suggest that recent activity along this fault system reflects incipient deformation of the foreland at the northeastern margin of the Tibetan Plateau.

Characterizing the transient geomorphic response to base‐level fall in the northeastern Tibetan Plateau

Analysis of hillslope gradient, landscape relief, and channel steepness in the Daxia River basin provides evidence of a transient geomorphic response to base-level fall on the northeastern Tibetan Plateau. Low-gradient channels and gentle hillslopes of the upper watershed are separated from a steeper, high-relief landscape by a series of convex knickzones along channel longitudinal profiles. Downstream projection of the “relict” portions of the profiles implies ~800–850 m of incision, consistent with geologic and geomorphic records of post ~1.7 Ma incision in the lower watershed. We combine optically stimulated luminescence dating of fluvial terrace deposits to constrain incision rates downstream of knickpoints with catchment-averaged 10Be concentrations in modern sediment to estimate erosion rates in tributary basins both above and below knickpoints. Both sources of data imply landscape lowering rates of ~300 m Ma−1 below the knickpoint and ~50–100 m Ma−1 above. Field measurements of channel width (n = 48) and calculations of bankfull discharge (n = 9) allow determination of scaling relations among channel hydraulic geometry, discharge, and contributing area that we employ to estimate the patterns of basal shear stress, unit stream power, and bed load transport rate throughout the channel network. Our results imply a clear downstream increase of incision potential; this result would appear to be consistent with a detachment-limited response to imposed base-level fall, in which steepening of channels drives an increase in erosion rates. In contrast, however, we do not observe apparent narrowing of channels across the transition from slowly eroding to rapidly eroding portions of the watershed. Rather, the present-day channel morphology as well as its scaling of hydraulic geometry imply that the river is primarily adjusted to transport its sediment load and suggest that channel morphology may not always reflect the presence of knickpoints and differences in landscape relief.

Evaluating the size and extent of paleolakes in central Tibet during the late Pleistocene

Subhorizontal lake shorelines allow a geodynamic test of the size and extent of a hypothesized paleo-lake in central Tibet, the East Qiangtang Lake (EQL), during the last interglacial period (Marine Isotope Stage [MIS] 5e). Reconstructions based on relict lake deposits suggest that the EQL would have been ~400 m deep and over ~66,000 km2. Models of flexural rebound driven by lake recession predict that shorelines near the EQL center, at the present-day location of Siling Co, would have rebounded 60-90 m above their initial elevation. New 36Cl chronology of the highest relict shorelines around Siling Co indicates that they reflect lake levels between 110-190 ka. These shorelines, however, are presently >300 m below their predicted elevations, implying a substantially smaller water load. Our results reveal that the expansion of Tibetan lakes during MIS 5e was relatively limited. Instead, individual lakes were supplied by river networks, much as they are today.

Slip inversion along inner fore-arc faults, eastern Tohoku, Japan

Funding was provided by a grant from the National Science Foundation Tectonics Program grant EAR-0809939 to D.M.F. and E.K., Geologic Society of America Graduate Research Grants, and the P.D. Krynine Memorial Fund. The authors thank Gaku Kimura, Kyoko Tonegawa, Hiroko Watanabe, Jun Kameda, and Asuka Yamaguchi for scientific and logistical support, and Kristin Morell for comments on early versions of the manuscript. We also thank Yuzuru Yamamoto and Kohtaro Ujiie for their detailed reviews and suggestions for improvement to the manuscript. The authors acknowledge the use of the Move Software Suite granted by Midland Valleys Academic Software Initiative. Geologic, structural, stratigraphic, and chronologic data used herein are accessible in manuscript figures, and in the citations therein. Input geologic data for trishear kinematic modeling can be accessed in Table 1 and in the supporting information. (EAR-0809939 - National Science Foundation Tectonics Program grant; Geologic Society of America Graduate Research Grants; P.D. Krynine Memorial Fund)