Refine
Year of publication
Document Type
- Article (37)
- Postprint (8)
- Habilitation Thesis (1)
- Review (1)
Language
- English (47)
Is part of the Bibliography
- yes (47)
Keywords
- hydropower (4)
- Graph theory (3)
- Himalayas (3)
- Nepal (3)
- glacial hazards (3)
- glacial lake outburst floods (3)
- water resources (3)
- Himalaya (2)
- Maritime Alps (2)
- WorldView-2 (2)
Increased rates of glacier retreat and thinning need accurate local estimates of glacier elevation change to predict future changes in glacier runoff and their contribution to sea level rise. Glacier elevation change is typically derived from digital elevation models (DEMs) tied to surface change analysis from satellite imagery. Yet, the rugged topography in mountain regions can cast shadows onto glacier surfaces, making it difficult to detect local glacier elevation changes in remote areas. A rather untapped resource comprises precise, time-stamped metadata on the solar position and angle in satellite images. These data are useful for simulating shadows from a given DEM. Accordingly, any differences in shadow length between simulated and mapped shadows in satellite images could indicate a change in glacier elevation relative to the acquisition date of the DEM. We tested this hypothesis at five selected glaciers with long-term monitoring programmes. For each glacier, we projected cast shadows onto the glacier surface from freely available DEMs and compared simulated shadows to cast shadows mapped from ∼40 years of Landsat images. W validated the relative differences with geodetic measurements of glacier elevation change where these shadows occurred. We find that shadow-derived glacier elevation changes are consistent with independent photogrammetric and geodetic surveys in shaded areas. Accordingly, a shadow cast on Baltoro Glacier (the Karakoram, Pakistan) suggests no changes in elevation between 1987 and 2020, while shadows on Great Aletsch Glacier (Switzerland) point to negative thinning rates of about 1 m yr−1 in our sample. Our estimates of glacier elevation change are tied to occurrence of mountain shadows and may help complement field campaigns in regions that are difficult to access. This information can be vital to quantify possibly varying elevation-dependent changes in the accumulation or ablation zone of a given glacier. Shadow-based retrieval of glacier elevation changes hinges on the precision of the DEM as the geometry of ridges and peaks constrains the shadow that we cast on the glacier surface. Future generations of DEMs with higher resolution and accuracy will improve our method, enriching the toolbox for tracking historical glacier mass balances from satellite and aerial images.
Channel steepness index, k(s), is a metric derived from the stream power model that, under certain conditions, scales with relative rock uplift rate. Channel steepness index is a property of rivers, which can be relatively easily extracted from digital elevation models (DEMs). As DEM data sets are widely available for Earth and are becoming more readily available for other planetary bodies, channel steepness index represents a powerful tool for interpreting tectonic processes. However, multiple approaches to calculate channel steepness index exist. From this several important questions arise; does choice of approach change the values of channel steepness index, can values be so different that choice of approach can influence the findings of a study, and are certain approaches better than others? With the aid of a synthetic river profile and a case study from the Sierra Nevada, California, we show that values of channel steepness index vary over orders of magnitude according to the methodology used in the calculation. We explore the limitations, advantages and disadvantages of the key approaches to calculating channel steepness index, and find that choosing an appropriate approach relies on the context of a study. Given these observations, it is important that authors acknowledge the methodology used to calculate channel steepness index, to ensure that results can be contextualised and reproduced.
A dataset of 2184 field measurements reported in the literature was used to evaluate the predictive capability of eight conventional flow resistance equations to predict the mean flow velocity in gravel-bed rivers. The results reveal considerable disagreement with the observed flow velocities for relative submergence less than 4 and for the non-uniformity of the bed material greater than 7.5 for all the equations. However, the predictions made using the Smart and Jaggi (1983), Ferguson (2007), and Rickenmann and Recking (2011) equations were closer to the observed values. Furthermore, bedload sediment transport also reduces the predictive capability of the equations considered in this study except for the Recking et al. (2008) equation, which was developed consid- ering active bedload transport. The performance of flow resistance equations improves when corrected by considering the geometric standard deviation of the bed material. Here we present an empirical approach using the whole dataset and its subsets for accounting for the additional energy losses occurring due to the wake vortices, spill losses, and free surface instabilities occurring due to the protrusions from the bed. The results obtained using the validation dataset shows the importance and usefulness of this approach to account for the additional energy losses, especially for the Strickler (1923) and Keulegan (1938) equations.
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.
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.
We propose a novel way to measure and analyze networks of drainage divides from digital elevation models. We developed an algorithm that extracts drainage divides based on the drainage basin boundaries defined by a stream network. In contrast to streams, there is no straightforward approach to order and classify divides, although it is intuitive that some divides are more important than others. A meaningful way of ordering divides is the average distance one would have to travel down on either side of a divide to reach a common stream location. However, because measuring these distances is computationally expensive and prone to edge effects, we instead sort divide segments based on their tree-like network structure, starting from endpoints at river confluences. The sorted nature of the network allows for assigning distances to points along the divides, which can be shown to scale with the average distance downslope to the common stream location. Furthermore, because divide segments tend to have characteristic lengths, an ordering scheme in which divide orders increase by 1 at junctions mimics these distances. We applied our new algorithm to the Big Tujunga catchment in the San Gabriel Mountains of southern California and studied the morphology of the drainage divide network. Our results show that topographic metrics, like the downstream flow distance to a stream and hillslope relief, attain characteristic values that depend on the drainage area threshold used to derive the stream network. Portions along the divide network that have lower than average relief or are closer than average to streams are often distinctly asymmetric in shape, suggesting that these divides are unstable. Our new and automated approach thus helps to objectively extract and analyze divide networks from digital elevation models.
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.
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.
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.
Mountain rivers respond to strong earthquakes by rapidly aggrading to accommodate excess sediment delivered by co-seismic landslides. Detailed sediment budgets indicate that rivers need several years to decades to recover from seismic disturbances, depending on how recovery is defined. We examine three principal proxies of river recovery after earthquake-induced sediment pulses around Pokhara, Nepal's second largest city. Freshly exhumed cohorts of floodplain trees in growth position indicate rapid and pulsed sedimentation that formed a fan covering 150 km2 in a Lesser Himalayan basin with tens of metres of debris between the 11th and 15th centuries AD. Radiocarbon dates of buried trees are consistent with those of nearby valley deposits linked to major medieval earthquakes, such that we can estimate average rates of re-incision since. We combine high-resolution digital elevation data, geodetic field surveys, aerial photos, and dated tree trunks to reconstruct geomorphic marker surfaces. The volumes of sediment relative to these surfaces require average net sediment yields of up to 4200 t km–2 yr–1 for the 650 years since the last inferred earthquake-triggered sediment pulse. The lithological composition of channel bedload differs from that of local bedrock, confirming that rivers are still mostly evacuating medieval valley fills, locally incising at rates of up to 0.2 m yr–1. Pronounced knickpoints and epigenetic gorges at tributary junctions further illustrate the protracted fluvial response; only the distal portions of the earthquake-derived sediment wedges have been cut to near their base. Our results challenge the notion that mountain rivers recover speedily from earthquakes within years to decades. The valley fills around Pokhara show that even highly erosive Himalayan rivers may need more than several centuries to adjust to catastrophic perturbations. Our results motivate some rethinking of post-seismic hazard appraisals and infrastructural planning in active mountain regions.
Plain Language Summary The 2015 Gorkha earthquake in Nepal caused severe losses in the hydropower sector. The country temporarily lost similar to 20% of its hydropower capacity, and >30 hydropower projects were damaged. The projects hit hardest were those that were affected by earthquake-triggered landslides. We show that these projects are located along very steep rivers with towering sidewalls that are prone to become unstable during strong seismic ground shaking. A statistical classification based on a topographic metric that expresses river steepness and earthquake ground acceleration is able to approximately predict hydropower damage during future earthquakes, based on successful testing of past cases. Thus, our model enables us to estimate earthquake damages to hydropower projects in other parts of the Himalayas. We find that >10% of the Himalayan drainage network may be unsuitable for hydropower infrastructure given high probabilities of high earthquake damages.
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.
Bumps in river profiles: uncertainty assessment and smoothing using quantile regression techniques
(2017)
The analysis of longitudinal river profiles is an important tool for studying landscape evolution. However, characterizing river profiles based on digital elevation models (DEMs) suffers from errors and artifacts that particularly prevail along valley bottoms. The aim of this study is to characterize uncertainties that arise from the analysis of river profiles derived from different, near-globally available DEMs. We devised new algorithms quantile carving and the CRS algorithm - that rely on quantile regression to enable hydrological correction and the uncertainty quantification of river profiles. We find that globally available DEMs commonly overestimate river elevations in steep topography. The distributions of elevation errors become increasingly wider and right skewed if adjacent hillslope gradients are steep. Our analysis indicates that the AW3D DEM has the highest precision and lowest bias for the analysis of river profiles in mountainous topography. The new 12m resolution TanDEM-X DEM has a very low precision, most likely due to the combined effect of steep valley walls and the presence of water surfaces in valley bottoms. Compared to the conventional approaches of carving and filling, we find that our new approach is able to reduce the elevation bias and errors in longitudinal river profiles.
Remote sensing technology serves as a powerful tool for analyzing geospatial characteristics of flood inundation events at various scales. However, the performance of remote sensing methods depends heavily on the flood characteristics and landscape settings. Difficulties might be encountered in mapping the extent of localized flooding with shallow water on riverine floodplain areas, where patches of herbaceous vegetation are interspersed with open water surfaces. To address the difficulties in mapping inundation on areas with complex water and vegetation compositions, a high spatial resolution dataset has to be used to reduce the problem of mixed pixels. The main objective of our study was to investigate the possibilities of using a single date WorldView-2 image of very high spatial resolution and supporting data to analyze spatial patterns of localized flooding on a riverine floodplain. We used a decision tree algorithm with various combinations of input variables including spectral bands of the WorldView-2 image, selected spectral indices dedicated to mapping water surfaces and vegetation, and topographic data. The overall accuracies of the twelve flood extent maps derived with the decision tree method and performed on both pixels and image objects ranged between 77% and 95%. The highest mapping overall accuracy was achieved with a method that utilized all available input data and the object-based image analysis. Our study demonstrates the possibility of using single date WorldView-2 data for analyzing flooding events at high spatial detail despite the absence of spectral bands from the short-waveform region that are frequently used in water related studies. Our study also highlights the importance of topographic data in inundation analyses. The greatest difficulties were met in mapping water surfaces under dense canopy herbaceous vegetation, due to limited water surface exposure and the dominance of vegetation reflectance.
Bumps in river profiles
(2017)
The analysis of longitudinal river profiles is an important tool for studying landscape evolution. However, characterizing river profiles based on digital elevation models (DEMs) suffers from errors and artifacts that particularly prevail along valley bottoms. The aim of this study is to characterize uncertainties that arise from the analysis of river profiles derived from different, near-globally available DEMs. We devised new algorithms quantile carving and the CRS algorithm - that rely on quantile regression to enable hydrological correction and the uncertainty quantification of river profiles. We find that globally available DEMs commonly overestimate river elevations in steep topography. The distributions of elevation errors become increasingly wider and right skewed if adjacent hillslope gradients are steep. Our analysis indicates that the AW3D DEM has the highest precision and lowest bias for the analysis of river profiles in mountainous topography. The new 12m resolution TanDEM-X DEM has a very low precision, most likely due to the combined effect of steep valley walls and the presence of water surfaces in valley bottoms. Compared to the conventional approaches of carving and filling, we find that our new approach is able to reduce the elevation bias and errors in longitudinal river profiles.