skip to main content


Title: The Importance of Hillslope Scale in Responses of Chemical Erosion Rate to Changes in Tectonics and Climate
Abstract

Chemical erosion is of wide interest due to its influence on topography, nutrient supply to streams and soils, sediment composition, and Earth's climate. While controls on chemical erosion rate have been studied extensively in steady‐state models, few studies have explored the controls on chemical erosion rate during transient responses to external perturbations. Here we develop a numerical model for the coevolution of soil‐mantled topography, soil thickness, and soil mineralogy, and we use it to simulate responses to step changes in rates of rock uplift, soil production, soil transport, and mineral dissolution. These simulations suggest that tectonic and climatic perturbations can generate responses in soil chemical erosion rate that differ in speed, magnitude, and spatial pattern and that climatic and tectonic perturbations may impart distinct signatures on hillslope mass fluxes, soil chemistry, and sediment composition. The response time of chemical erosion rate is dominantly controlled by hillslope length and is secondarily modulated by rates of rock uplift, soil production, transport, and mineral dissolution. This strong dependence on drainage density implies that a landscape's chemical erosion response should depend on the relative efficiencies of river incision and soil transport and thus may be mediated by climatic and biological factors. The simulations further suggest that the timescale of the hillslope response may be long relative to that of river channel profiles, implying that chemical erosion response times may be limited more by the sluggishness of the hillslopes than by the rate of signal propagation through river channel profiles.

 
more » « less
Award ID(s):
1946762
NSF-PAR ID:
10449122
Author(s) / Creator(s):
 ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Earth Surface
Volume:
125
Issue:
9
ISSN:
2169-9003
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Temporal and spatial variations of tectonic rock uplift are generally thought to be the main controls on long-term erosion rates in various landscapes. However, rivers continuously lengthen and capture drainages in strike-slip fault systems due to ongoing motion across the fault, which can induce changes in landscape forms, drainage networks, and local erosion rates. Located along the restraining bend of the San Andreas Fault, the San Bernardino Mountains provide a suitable location for assessing the influence of topographic disequilibrium from perturbations by tectonic forcing and channel reorganization on measured erosion rates. In this study, we measured 17 new basin-averaged erosion rates using cosmogenic 10Be in river sands (hereafter, 10Be-derived erosion rates) and compiled 31 10Be-derived erosion rates from previous work. We quantify the degree of topographic disequilibrium using topographic analysis by examining hillslope and channel decoupling, the areal extent of pre-uplift surface, and drainage divide asymmetry across various landscapes. Similar to previous work, we find that erosion rates generally increase from north to south across the San Bernardino Mountains, reflecting a southward increase in tectonic activity. However, a comparison between 10Be-derived erosion rates and various topographic metrics in the southern San Bernardino Mountains suggests that the presence of transient landscape features such as relict topography and drainage-divide migration may explain local variations in 10Be-derived erosion rates. Our work shows that coupled analysis of erosion rates and topographic metrics provides tools for assessing the influence of tectonic uplift and channel reorganization on landscape evolution and 10Be-derived erosion rates in an evolving strike-slip restraining bend.

     
    more » « less
  2. Abstract

    River profiles are shaped by climatic and tectonic history, lithology, and internal feedbacks between flow hydraulics, sediment transport and erosion. In steep channels, waterfalls may self‐form without changes in external forcing (i.e., autogenic formation) and erode at rates faster or slower than an equivalent channel without waterfalls. We use a 1‐D numerical model to investigate how self‐formed waterfalls alter the morphology of bedrock river longitudinal profiles. We modify the standard stream power model to include a slope threshold above which waterfalls spontaneously form and a rate constant allowing waterfalls to erode faster or slower than other fluvial processes. Using this model, we explore how waterfall formation alters both steady state and transient longitudinal profile forms. Our model predicts that fast waterfalls create km‐scale reaches in a dynamic equilibrium with channel slope held approximately constant at the threshold slope for waterfall formation, while slow waterfalls can create local channel slope maxima at the location of slow waterfall development. Furthermore, slow waterfall profiles integrate past base level histories, leading to multiple possible profile forms, even at steady‐state. Consistency between our model predictions and field observations of waterfall‐rich rivers in the Kings and Kaweah drainages in the southern Sierra Nevada, California, supports the hypothesis that waterfall formation can modulate river profiles in nature. Our findings may help identify how bedrock channels are influenced by waterfall erosion and aid in distinguishing between signatures of external and internal perturbations, thereby strengthening our ability to interpret past climate and tectonic changes from river longitudinal profiles.

     
    more » « less
  3. This study examines dissolved rhenium (Re) concentrations as a function of water runoff using river samples from two contrasting mountainous watersheds, the Eel and Umpqua Rivers in the Pacific Northwest, USA. These watersheds share many key characteristics in terms of size, discharge, climate, and vegetation, but they have a 15-fold difference in sediment yield due to differences in their tectonic setting and uplift and erosion rates. We evaluate concentration-runoff (C-R) relationships and ratios of coefficients of variation (CVC/CVR) for major cations, anions, dissolved inorganic carbon, selected trace elements including Re, and 87Sr/86Sr ratios. Recent research outlines the potential of Re to serve as a tracer for the oxidation of ancient/fossil organic matter because of its close association with petrogenic carbon (OCpetro) in rocks. In both the Eel and Umpqua Rivers, our measurements show that Re behaves similarly to major weathering derived-solutes corrected for atmospheric input, such as Ca2+*, Mg2+*, and Na+* with modest dilution across all tributaries with increasing runoff. Rhenium behaves dissimilarly from other trace elements, such as Mo and U, and is also dissimilar to biologically-cycled nutrients, such as NO3 – , PO4 3 , and K+*, suggesting differences in sources, solute generation mechanisms, and flowpaths. Rhenium behavior is also distinct from that of colloids, which have increasing concentrations with increasing runoff. We find that Re and sulfate corrected for atmospheric input (SO4 2 *) have distinct CR relationships, in which SO4 2 * undergoes greater dilution with increasing runoff. This implies that Re is not dominantly sourced from sulfide weathering, which leaves primary bedrock minerals and OCpetro hosted in bedrock of these watersheds as the likely dominant sources of dissolved Re release. At mean discharge, Re concentration in the Eel river (3.5 pmol L-1) is more than two times greater than Re concentrations in the Umpqua River (1.5 pmol L-1). Furthermore, comparison of two tributary watersheds with similar bedrock but marked differences in erosion rates show higher Re concentrations in Bull Creek (erosion rate of 0.5 mm yr 1) relative to Elder Creek (erosion rate of 0.2 mm yr 1). The results of this study suggest that dissolved Re in the Eel and Umpqua River basins is likely derived from primary mineral dissolution or OCpetro oxidation, and Re fluxes are higher in areas with higher erosion rates, suggesting that tectonic setting is one factor that controls Re release and therefore OCpetro oxidation. 
    more » « less
  4. null (Ed.)
    The Indian (southwest) summer monsoon is one of the most intense climatic phenomena on Earth. Its long-term development has been linked to the growth of high topography in South and Central Asia. The Indian continental margin, adjoining the Arabian Sea, offers a unique opportunity to investigate tectonic–climatic interactions and the net impact of these processes on weathering and erosion of the western Himalaya. During International Ocean Discovery Program Expedition 355, two sites (U1456 and U1457) were drilled in Laxmi Basin in the eastern Arabian Sea to document the coevolution of mountain building, weathering, erosion, and climate over a range of timescales. In addition, recovering basement from the eastern Arabian Sea provides constraints on the early rifting history of the western continental margin of India with special emphasis on continental breakup between India and the Seychelles and its relationship to the plume-related volcanism of the Deccan Plateau. Drilling and coring operations during Expedition 355 recovered sediment from Sites U1456 and U1457 in the Laxmi Basin, penetrating 1109.4 and 1108.6 m below seafloor (mbsf), respectively. Drilling reached sediment dated to 13.5–17.7 Ma (late early to early middle Miocene) at Site U1456, although with a large hiatus between the lowermost sediment and overlying deposits dated to <10.9 Ma. At Site U1457, a much longer hiatus occurs near the base of the cored section, spanning from 10.9 to ~62 Ma. At both sites, hiatuses span ~8.2–9.2 and ~3.6–5.6 Ma, with a possible condensed section spanning ~2.0–2.6 Ma, although the total duration for each hiatus is slightly different between the two sites. A major submarine fan draining the western Himalaya and Karakoram must have been supplying sediment to the eastern Arabian Sea since at least ~17 Ma. Sand mineral assemblages indicate that the Greater Himalayan Crystalline Sequence was fully exposed to the surface by this time. Most of the recovered sediment appears to be derived from the Indus River and includes minerals that are unique to the Indus Suture Zone, in particular glaucophane and hypersthene, most likely originating from the structural base of the Kohistan arc. Pliocene sandy intervals at Site U1456 were deposited in lower fan “sheet lobe” settings, with intervals of basin plain turbidites separated by hemipelagic muddy sections deposited during the Miocene. Site U1457 is more distal in facies, reflecting its more marginal setting. No major active lobe appears to have affected the Laxmi Basin since the Middle Pleistocene (~1.2 Ma). We succeeded in recovering sections spanning the 8 Ma climatic transition, when monsoon intensity is believed to have changed strongly, although the nature of this change awaits postcruise analysis. We also recovered sediment from a large mass transport deposit measuring ~330 and ~190 m thick at Sites U1456 and U1457, respectively. This section includes an upper sequence of slump-folded muddy and silty rocks, as well as underlying calcarenites and limestone breccias, together with smaller amounts of volcanic clasts, all of which are likely derived from the western Indian continental shelf. Identification of similar facies on the regional seismic lines in Laxmi Basin suggests that these deposits form parts of one of the world’s largest mass transport deposits. Coring of igneous basement was successful at Site U1457. Recovery of massive basalt and associated volcaniclastic sediment at this site should address the key questions related to rifting and volcanism associated with formation of Laxmi Basin. Geochemical analysis is required to understand the petrogenesis and thus the tectonic setting of volcanism that will reveal whether it is oceanic basalt or volcanic rock contaminated by underlying continental crust or continental flood basalt. However, the fact that the lavas are massive and have few vesicles implies water depths of eruption likely deeper than 2000 m. This precludes opening of the basin in the presence of a major mantle thermal anomaly, such as that associated with the Deccan Large Igneous Province. Other observations made at the two sites during Expedition 355 provide vital constraints on the rift history of this margin. Heat flow measurements at the two drill sites were calculated to be ~57 and ~60 mW/m2. Such heat flow values are compatible with those observed in average oceanic crust of 63–84 Ma age, as well as with the presence of highly extended continental crust. Postcruise analyses of the more than ~1722 m of core will provide further information about the nature of tectonic–climatic interactions in this global type area for such studies. 
    more » « less
  5. Abstract

    How tectonic forcing, expressed as base level change, is encoded in the stratigraphic and geomorphic records of coupled source‐to‐sink systems remains uncertain. Using sedimentological, geochronological and geomorphic approaches, we describe the relationship between transient topographic change and sediment deposition for a low‐storage system forced by rapid rock uplift. We present five new luminescence ages and two terrestrial cosmogenic nuclide paleo‐erosion rates for the late Pleistocene Pagliara fan‐delta complex and we model corresponding base level fall history and erosion of the source catchment located on the Ionian flank of the Peloritani Mountains (NE‐Sicily, Italy). The Pagliara delta complex is part of the broader Messina Gravel‐and‐Sands lithostratigraphic unit that outcrops along the Peloritani coastal belt as extensional basins have been recently inverted by both normal faults and regional uplift at the Messina Straits. The deltas exposed at the mouth of the Pagliara River have constructional tops at ca. 300 m a.s.l. and onlap steeply east‐dipping bedrock at the coast to thickness between ca. 100 and 200 m. Five infrared‐stimulated luminescence (IRSL) ages collected from the delta range in age from ca. 327 to 208 ka and indicate a vertical long‐term sediment accumulation rate as rapid as ca. 2.2 cm/yr during MIS 7. Two cosmogenic10Be concentrations measured in samples of delta sediment indicate paleo‐erosion rates during MIS 8–7 near or slightly higher than the modern rates of ca. 1 mm/yr. Linear inversion of Pagliara fluvial topography indicates an unsteady base level fall history in phase with eustasy that is superimposed on a longer, tectonically driven trend that doubled in rate from ca. 0.95 to 1.8 mm/yr in the past 150 ky. The combination of footwall uplift rate and eustasy determines the accommodation space history to trap the fan‐deltas at the Peloritani coast in hanging wall basins, which are now inverted, uplifted and exposed hundreds of metres above the sea level.

     
    more » « less