Refine
Has Fulltext
- no (37)
Year of publication
Document Type
- Article (37) (remove)
Language
- English (37)
Is part of the Bibliography
- yes (37)
Keywords
- Nepal (3)
- hydropower (3)
- Graph theory (2)
- Himalaya (2)
- Himalayas (2)
- glacial hazards (2)
- glacial lake outburst floods (2)
- water resources (2)
- 3D CAVE (1)
- ALS (1)
Institute
Landscape evolution models (LEMs) allow the study of earth surface responses to changing climatic and tectonic forcings. While much effort has been devoted to the development of LEMs that simulate a wide range of processes, the numerical accuracy of these models has received less attention. Most LEMs use first-order accurate numerical methods that suffer from substantial numerical diffusion. Numerical diffusion particularly affects the solution of the advection equation and thus the simulation of retreating landforms such as cliffs and river knickpoints. This has potential consequences for the integrated response of the simulated landscape. Here we test a higher-order flux-limiting finite volume method that is total variation diminishing (TVD-FVM) to solve the partial differential equations of river incision and tectonic displacement. We show that using the TVD-FVM to simulate river incision significantly influences the evolution of simulated landscapes and the spatial and temporal variability of catchment-wide erosion rates. Furthermore, a two-dimensional TVD-FVM accurately simulates the evolution of landscapes affected by lateral tectonic displacement, a process whose simulation was hitherto largely limited to LEMs with flexible spatial discretization. We implement the scheme in TTLEM (TopoToolbox Landscape Evolution Model), a spatially explicit, raster-based LEM for the study of fluvially eroding landscapes in TopoToolbox 2.
The 2015 magnitude 7.8 Gorkha earthquake and its aftershocks weakened mountain slopes in Nepal. Co- and postseismic landsliding and the formation of landslide-dammed lakes along steeply dissected valleys were widespread, among them a landslide that dammed the Kali Gandaki River. Overtopping of the landslide dam resulted in a flash flood downstream, though casualties were prevented because of timely evacuation of low-lying areas. We hindcast the flood using the BREACH physically based dam-break model for upstream hydrograph generation, and compared the resulting maximum flow rate with those resulting from various empirical formulas and a simplified hydrograph based on published observations. Subsequent modeling of downstream flood propagation was compromised by a coarse-resolution digital elevation model with several artifacts. Thus, we used a digital-elevation-model preprocessing technique that combined carving and smoothing to derive topographic data. We then applied the 1-dimensional HEC-RAS model for downstream flood routing, and compared it to the 2-dimensional Delft-FLOW model. Simulations were validated using rectified frames of a video recorded by a resident during the flood in the village of Beni, allowing estimation of maximum flow depth and speed. Results show that hydrological smoothing is necessary when using coarse topographic data (such as SRTM or ASTER), as using raw topography underestimates flow depth and speed and overestimates flood wave arrival lag time. Results also show that the 2-dimensional model produces more accurate results than the 1-dimensional model but the 1-dimensional model generates a more conservative result and can be run in a much shorter time. Therefore, a 2-dimensional model is recommended for hazard assessment and planning, whereas a 1-dimensional model would facilitate real-time warning declaration.
The evaluation and verification of landscape evolution models (LEMs) has long been limited by a lack of suitable observational data and statistical measures which can fully capture the complexity of landscape changes. This lack of data limits the use of objective function based evaluation prolific in other modelling fields, and restricts the application of sensitivity analyses in the models and the consequent assessment of model uncertainties. To overcome this deficiency, a novel model function approach has been developed, with each model function representing an aspect of model behaviour, which allows for the application of sensitivity analyses. The model function approach is used to assess the relative sensitivity of the CAESAR-Lisflood LEM to a set of model parameters by applying the Morris method sensitivity analysis for two contrasting catchments. The test revealed that the model was most sensitive to the choice of the sediment transport formula for both catchments, and that each parameter influenced model behaviours differently, with model functions relating to internal geomorphic changes responding in a different way to those relating to the sediment yields from the catchment outlet. The model functions proved useful for providing a way of evaluating the sensitivity of LEMs in the absence of data and methods for an objective function approach.
Natural catchments are likely to show the existence of knickpoints in their river networks. The origin and genesis of the knickpoints can be manifold, considering that the present morphology is the result of the interactions of different factors such as tectonic movements, quaternary glaciations, river captures, variable lithology, and base-level changes. We analyzed the longitudinal profiles of the river channels in the Stura di Demonte Valley (Maritime Alps) to identify the knickpoints of such an alpine setting and to characterize their origins. The distribution and the geometry of stream profiles were used to identify the possible causes of the changes in stream gradients and to define zones with genetically linked knickpoints. Knickpoints are key geomorphological features for reconstructing the evolution of fluvial dissected basins, when the different perturbing factors affecting the ideally graded fluvial system have been detected. This study shows that even in a regionally small area, perturbations of river profiles are caused by multiple factors. Thus, attributing (automatically)-extracted knickpoints solely to one factor, can potentially lead to incomplete interpretations of catchment evolution.
Climate science is highly interdisciplinary by nature, so understanding interactions between Earth processes inherently warrants the use of analytical software that can operate across the disciplines of Earth science. Toward this end, we present the Climate Data Toolbox for MATLAB, which contains more than 100 functions that span the major climate-related disciplines of Earth science. The toolbox enables streamlined, entirely scriptable workflows that are intuitive to write and easy to share. Included are functions to evaluate uncertainty, perform matrix operations, calculate climate indices, and generate common data displays. Documentation is presented pedagogically, with thorough explanations of how each function works and tutorials showing how the toolbox can be used to replicate results of published studies. As a well-tested, well-documented platform for interdisciplinary collaborations, the Climate Data Toolbox for MATLAB aims to reduce time spent writing low-level code, let researchers focus on physics rather than coding and encourage more efficacious code sharing. Plain Language Summary This article describes a collection of computer code that has recently been released to help scientists analyze many types of Earth science data. The code in this toolbox makes it easy to investigate things like global warming, El Nino, or other major climate-related processes such as how winds affect ocean circulation. Although the toolbox was designed to be used by expert climate scientists, its instruction manual is well written, and beginners may be able to learn a great deal about coding and Earth science, simply by following along with the provided examples. The toolbox is intended to help scientists save time, help them ensure their analysis is accurate, and make it easy for other scientists to repeat the results of previous studies.
Knickpoints in longitudinal river profiles are proxies for the climatic and tectonic history of active mountains. The analysis of river profiles commonly relies on the assumption that drainage network configurations are stable. Here, we show that this assumption must be made cautiously if changes in contributing area are fast relative to knickpoint migration rates. We studied the Parachute Creek basin in the Roan Plateau, Colorado, United States, where knickpoint retreat occurs in horizontally uniform lithology so that drainage area is the sole governing variable. In this basin, we identified an anomalous catchment in the degree to which a stream power-based model predicted knickpoint locations. The catchment is experiencing area loss as the plateau edge is eroded by cliff migration in proximity to the Colorado River. Model predictions improve if the plateau edge is assumed to have migrated over the time scale of knickpoint retreat. Finally, a Lagrangian model of knickpoint migration enabled us to study the kinematic links between drainage area loss and knickpoint migration and offered constraints on the temporal aspects of area loss. Modeled onset and amount of area loss are consistent with cliff retreat rates along the margin of the Roan Plateau inferred from the incisional history of the upper Colorado River.
Drainage divide networks
(2020)
Drainage divides are organized into tree-like networks that may record information about drainage divide mobility. However, views diverge about how to best assess divide mobility. Here, we apply a new approach of automatically extracting and ordering drainage divide networks from digital elevation models to results from landscape evolution model experiments. We compared landscapes perturbed by strike-slip faulting and spatiotemporal variations in erodibility to a reference model to assess which topographic metrics (hillslope relief, flow distance, and chi) are diagnostic of divide mobility. Results show that divide segments that are a minimum distance of similar to 5 km from river confluences strive to attain constant values of hillslope relief and flow distance to the nearest stream. Disruptions of such patterns can be related to mobile divides that are lower than stable divides, closer to streams, and often asymmetric in shape. In general, we observe that drainage divides high up in the network, i.e., at great distances from river confluences, are more susceptible to disruptions than divides closer to these confluences and are thus more likely to record disturbance for a longer time period. We found that across-divide differences in hillslope relief proved more useful for assessing divide migration than other tested metrics. However, even stable drainage divide networks exhibit across-divide differences in any of the studied topographic metrics. Finally, we propose a new metric to quantify the connectivity of divide junctions.
The efficiency of sediment routing from land to the ocean depends on the position of submarine canyon heads with regard to terrestrial sediment sources. We aim to identify the main controls on whether a submarine canyon head remains connected to terrestrial sediment input during Holocene sea-level rise. Globally, we identified 798 canyon heads that are currently located at the 120m-depth contour (the Last Glacial Maximum shoreline) and 183 canyon heads that are connected to the shore (within a distance of 6 km) during the present-day highstand. Regional hotspots of shore-connected canyons are the Mediterranean active margin and the Pacific coast of Central and South America. We used 34 terrestrial and marine predictor variables to predict shore-connected canyon occurrence using Bayesian regression. Our analysis shows that steep and narrow shelves facilitate canyon-head connectivity to the shore. Moreover, shore-connected canyons occur preferentially along active margins characterized by resistant bedrock and high river-water discharge.
Purpose: Dry land vegetation is expected to respond sensitively to climate change and the projected variability of rainfall events. Rainfall as a water source is an obvious factor for the water supply of vegetation. However, the interaction of water and surface on rocky desert slopes with a patchy soil cover is also vital for vegetation in drylands. In particular, runoff on rocky surfaces and infiltration capacity of soil patches determine plant available water. Process-based studies into rock-soil interaction benefit from rainfall simulation, but require an approach accounting for the micro-scale heterogeneity of the slope surfaces. This study therefore aims at developing a suitable procedure for examining rock-soil interaction and the relevance of soil volume for storing plant available water in the northern Negev, Israel.
Materials and methods: To determine the amount of rainfall required to fill the available soil water storage capacity rainfall simulation experiments were conducted. The design of the rainfall-simulator and the selection of the plots aimed specifically at observing infiltration into small soil patches on a micro-scale relevant for the prevalent vegetation cover.
Results and discussion: The preliminary results of the study in the Negev Desert indicate that the ratio between soil volume and frequency of rainfall events determine the effect of climate change on plant available water and thus ultimately vegetation cover.
Conclusions: Based on the experiments examining runoff and soil moisture the qualitative understanding of hillslope ecohydrology in a rocky desert environment can be expanded into a quantitative assessment of the potential impact of varying rainfall conditions. The study also illustrates the contribution of rainfall simulation experiments for studies on the impact of climate change.
Beta diversity is a conceptual link between diversity at local and regional scales. Various additional methodologies of quantifying this and related phenomena have been applied. Among them, measures of pairwise (dis)similarity of sites are particularly popular. Undersampling, i.e. not recording all taxa present at a site, is a common situation in ecological data. Bias in many metrics related to beta diversity must be expected, but only few studies have explicitly investigated the properties of various measures under undersampling conditions. On the basis of an empirical data set, representing near-complete local inventories of the Lepidoptera from an isolated Pacific island, as well as simulated communities with varying properties, we mimicked different levels of undersampling. We used 14 different approaches to quantify beta diversity, among them dataset-wide multiplicative partitioning (i.e. true beta diversity') and pairwise site x site dissimilarities. We compared their values from incomplete samples to true results from the full data. We used these comparisons to quantify undersampling bias and we calculated correlations of the dissimilarity measures of undersampled data with complete data of sites. Almost all tested metrics showed bias and low correlations under moderate to severe undersampling conditions (as well as deteriorating precision, i.e. large chance effects on results). Measures that used only species incidence were very sensitive to undersampling, while abundance-based metrics with high dependency on the distribution of the most common taxa were particularly robust. Simulated data showed sensitivity of results to the abundance distribution, confirming that data sets of high evenness and/or the application of metrics that are strongly affected by rare species are particularly sensitive to undersampling. The class of beta measure to be used should depend on the research question being asked as different metrics can lead to quite different conclusions even without undersampling effects. For each class of metric, there is a trade-off between robustness to undersampling and sensitivity to rare species. In consequence, using incidence-based metrics carries a particular risk of false conclusions when undersampled data are involved. Developing bias corrections for such metrics would be desirable.