skip to main content

Title: Hydro‐geomorphic behavior of contrasting tropical landscapes and critical zone response to changing climate

The potentially important influence of climate change on landscape evolution and on critical zone processes is not sufficiently understood. The relative contribution of hydro‐climatic factors on hillslope erosion rates may significantly vary with topography at the watershed scale. The objective of this study is to quantify the hydro‐geomorphic behavior of two contrasting landscapes in response to different climate change scenarios in the Luquillo Critical Zone Observatory, a site of particular geomorphological interest, in terms of hillslope erosion and rainfall‐triggered landslides. We investigate the extent to which hillslope erosion and landslide occurrence remain relatively invariant with future hydro‐climatic perturbations. The adjacent Mameyes and Icacos watersheds are studied, which are underlain by contrasting lithologies. A high resolution coupled hydro‐geomorphic model based on tRIBS (Triangulated Irregular Network‐based Real‐time Integrated Basin Simulator) is used. Observations of landslide activity and hillslope erosion are used to evaluate the model performance. The process‐based model quantifies feedbacks among different hydrologic processes, landslide occurrence, and topsoil erosion and deposition. Simulations suggest that the propensity for landslide occurrence in the Luquillo Mountains is controlled by tropical storms, subsurface water flow, and by non‐climatic factors, and is expected to remain significant through 2099. The Icacos watershed, which is underlain by quartz diorite, is dominated by relatively large landslides. The relative frequency of smaller landslides is higher at the Mameyes watershed, which is underlain by volcaniclastic rock. While projections of precipitation decrease at the study site may lead to moderate decline in hillslope erosion rates, the simulated erosional potential of the two diverse landscapes likely remains significant. © 2018 John Wiley & Sons, Ltd.

more » « less
Author(s) / Creator(s):
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Earth Surface Processes and Landforms
Page Range / eLocation ID:
p. 641-654
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The dynamics of soil organic carbon (SOC) in tropical forests play an important role in the global carbon (C) cycle. Past attempts to quantify the net C exchange with the atmosphere in regional and global budgets do not systematically account for dynamic feedbacks among linked hydrological, geomorphological, and biogeochemical processes, which control the fate of SOC. Here we quantify effects of geomorphic perturbations on SOC oxidation and accumulation in two adjacent wet tropical forest watersheds underlain by contrasting lithology (volcaniclastic rock and quartz diorite) in the Luquillo Critical Zone Observatory. This study uses the spatially explicit and physically based model of SOC dynamics tRIBS‐ECO (Triangulated Irregular Network‐based Real‐time Integrated Basin Simulator‐Erosion and Carbon Oxidation) and measurements of SOC profiles and oxidation rates. Our results suggest that hillslope erosion at the two watersheds may drive C sequestration or CO2release to the atmosphere, depending on the forest type and land use. The net erosion‐induced C exchange with the atmosphere was controlled by the spatial distribution of forest types. The two watersheds were characterized by significant erosion and dynamic replacement of upland SOC stocks. Results suggest that the landscape underlain by volcaniclastic rock has reached a state close to geomorphic equilibrium, and the landscape underlain by quartz diorite is characterized by greater rates of denudation. These findings highlight the importance of the spatially explicit and physical representation of C erosion driven by local variation in lithological and geomorphological characteristics and in forest cover.

    more » « less
  2. Abstract

    The connection between topography and erosion rate is central to understanding landscape evolution and sediment hazards. However, investigation of this relationship in steep landscapes has been limited due to expectations of: (a) decoupling between erosion rate and “threshold” hillslope morphology; and (b) bias in detrital cosmogenic nuclide erosion rates due to deep‐seated landslides. Here we compile 120 new and published10Be erosion rates from catchments in the San Gabriel Mountains, California, and show that hillslope morphology and erosion rate are coupled for slopes approaching 50° due to progressive exposure of bare bedrock with increasing erosion rate. We find no evidence for drainage area dependence in10Be erosion rates in catchments as small as 0.09 km2, and we show that landslide deposits influence erosion rate estimates mainly by adding scatter. Our results highlight the potential and importance of sampling small catchments to better understand steep hillslope processes.

    more » « less
  3. Abstract

    Debris flows pose persistent hazards and shape high‐relief landscapes in diverse physiographic settings, but predicting the spatiotemporal occurrence of debris flows in postglacial topography remains challenging. To evaluate the debris flow process in high‐relief postglacial terrain, we conducted a geomorphic investigation to characterize geologic, glacial, volcanic, and land use contributions to landslide initiation across Southeast Alaska. To evaluate controls on landslide (esp. debris flow) occurrence in Sitka, we used field observation, geomorphic mapping, landslide characteristics as documented in the Tongass National Forest inventory, and a novel application of the shallow landslide model SHALSTAB to postglacial terrain. A complex geomorphic history of glaciation and volcanic activity provides a template for spatially heterogeneous landslide occurrence. Landslide density across the region is highly variable, but debris flow density is high on south‐ or southeast‐facing hillslopes where volcanic tephra soils are present and/or where timber harvest has occurred since 1900. High landslide density along the western coast of Baranof and Kruzof islands coincides with deposition of glacial sediment and thick tephra and exposure to extreme rainfall from atmospheric rivers on south‐facing aspects but the relative contributions of these controls are unclear. Timber harvest has also been identified as an important control on landslide occurrence in the region. Focusing on a subset of geo‐referenced landslides near Sitka, we used the SHALSTAB shallow landslide initiation model, which has been frequently applied in non‐glacial terrain, to identify areas of high landslide potential in steep, convergent terrain. In a validation against mapped landslide polygons, the model significantly outperformed random guessing, with area under the curve (AUC) = 0.709 on a performance classification curve of true positives vs. false positives. This successful application of SHALSTAB demonstrates practical utility for hazards analysis in postglacial landscapes to mitigate risk to people and infrastructure.

    more » « less
  4. Abstract

    Landslides are common natural disturbances in tropical montane forests. While the geomorphic drivers of landslides in the Andes have been studied, factors controlling post‐landslide forest recovery across the steep climatic and topographic gradients characteristic of tropical mountains are poorly understood.

    Here we use a LiDAR‐derived canopy height map coupled with a 25‐year landslide time‐series map to examine how landslide, topographic and biophysical factors, along with residual vegetation, affect canopy height and heterogeneity in regenerating landslides. We also calculate above‐ground biomass accumulation rates and estimate the time for landslides to recover to mature forest biomass levels.

    We find that age and elevation are the biggest determinants of forest recovery, and that the jump‐start in regeneration that residual vegetation provides lasts for at least 18 years. Our estimates of time to biomass recovery (31.6–37.1 years) are surprisingly rapid, and as a result we recommend that future research pair LiDAR with hyperspectral imagery to estimate forest above‐ground biomass in frequently disturbed landscapes.

    Synthesis. Using a high‐resolution LiDAR dataset and a time‐series inventory of 608 landslides distributed across a wide elevational gradient in Andean montane forest, we show that age and elevation are the most influential predictors of forest canopy height and canopy variability. Other features of landslides, in particular the presence of residual vegetation, shape post‐landslide regeneration trajectories. LiDAR allows for a detailed analysis of forest structural recovery across large landscapes and numbers of disturbances, and provides a reasonable upper bound on above‐ground biomass accumulation rates. However, because this method does not capture the effect of compositional change through succession on above‐ground biomass, wherein high‐wood density species gradually replace light‐wooded pioneer species, it overestimates above‐ground biomass. Given previously estimated stem turnover rates along this elevational gradient, we posit that above‐ground biomass recovery takes at least three times as long as our recovery time estimates based on LiDAR‐derived structure alone.

    more » « less
  5. Abstract

    Evidence suggests that following earthquakes, landslide activity, and susceptibility remain elevated for months to years after cessation of shaking. In landslides constrained within the soil mantle, enhanced post‐seismic activity may be attributed to persistence of progressive failure surface development that results in delayed triggering. With subsequent post‐seismic hydrologic disturbance, a commensurate decay in activity over time reflects evacuation of damaged hillslope materials and mechanical healing of colluvial media and vegetation recovery. In this study, we propose a hydromechanical model for representing the interplay of these processes including trends in timing and magnitude of post‐seismic shallow landslide activity. The hydromechanical model considers seismic or hydrologic disturbances as drivers for failure surface development in the soil mantle over a hillslope with commensurate weakening through accumulation of internal compression, shear softening, and loss of root strength. In the absence of well‐instrumented post‐seismic landslide measurement of evolving stress conditions and deformation, predictions were made against a unique field case study where hydrologically induced hillslope failure kinematics and internal compression were measured. The model accounts for seismic disturbances, demonstrating that certain peak ground accelerations directly reduce subsequent rainfall‐induced landslide‐triggering thresholds relative to undisturbed hillslopes. The model framework can be used to describe progressive, post‐seismic landslide activity, demonstrating that consideration of root or soil healing may suppress subsequent landslide triggering and contribute toward recovery to pre‐seismic disturbance levels of strength. Factors such as shear softening, vegetation, and healing are explored in the context of the temporal evolution of landslide rates after a seismic event. The results of this paper show mechanistically how (1) progressive development or cessation of a failure surface in active, post‐seismic landslides may describe elevated post‐seismic landslide rates and (2) that hillslope healing plays a modest role in the magnitude and timescale of attenuated landslide activity.

    more » « less