@article{BeckHollowaySchwanghart2013, author = {Beck, Jan and Holloway, Jeremy D. and Schwanghart, Wolfgang}, title = {Undersampling and the measurement of beta diversity}, series = {Methods in ecology and evolution : an official journal of the British Ecological Society}, volume = {4}, journal = {Methods in ecology and evolution : an official journal of the British Ecological Society}, number = {4}, publisher = {Wiley-Blackwell}, address = {Hoboken}, issn = {2041-210X}, doi = {10.1111/2041-210x.12023}, pages = {370 -- 382}, year = {2013}, abstract = {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.}, language = {en} } @article{BernhardtSchwanghart2021, author = {Bernhardt, Anne and Schwanghart, Wolfgang}, title = {Where and why do submarine canyons remain connected to the shore during sea-level rise?}, series = {Geophysical research letters : GRL / American Geophysical Union}, volume = {48}, journal = {Geophysical research letters : GRL / American Geophysical Union}, number = {10}, publisher = {American Geophysical Union}, address = {Washington}, issn = {0094-8276}, doi = {10.1029/2020GL092234}, pages = {15}, year = {2021}, abstract = {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.}, language = {en} } @article{BernhardtSchwanghartHebbelnetal.2017, author = {Bernhardt, Anne and Schwanghart, Wolfgang and Hebbeln, Dierk and Stuut, Jan-Berend W. and Strecker, Manfred}, title = {Immediate propagation of deglacial environmental change to deep-marine turbidite systems along the Chile convergent margin}, series = {Earth \& planetary science letters}, volume = {473}, journal = {Earth \& planetary science letters}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0012-821X}, doi = {10.1016/j.epsl.2017.05.017}, pages = {190 -- 204}, year = {2017}, abstract = {Understanding how Earth-surface processes respond to past climatic perturbations is crucial for making informed predictions about future impacts of climate change on sediment "uxes. Sedimentary records provide the archives for inferring these processes, but their interpretation is compromised by our incomplete understanding of how sediment-routing systems respond to millennial-scale climate cycles. We analyzed seven sediment cores recovered from marine turbidite depositional sites along the Chile continental margin. The sites span a pronounced arid-to-humid gradient with variable relief and related sediment connectivity of terrestrial and marine environments. These sites allowed us to study event related depositional processes in different climatic and geomorphic settings from the Last Glacial Maximum to the present day. The three sites reveal a steep decline of turbidite deposition during deglaciation. High rates of sea-level rise postdate the decline in turbidite deposition. Comparison with paleoclimate proxies documents that the spatio-temporal sedimentary pattern rather mirrors the deglacial humidity decrease and concomitant warming with no resolvable lag times. Our results let us infer that declining deglacial humidity decreased "uvial sediment supply. This signal propagated rapidly through the highly connected systems into the marine sink in north-central Chile. In contrast, in south-central Chile, connectivity between the Andean erosional zone and the "uvial transfer zone probably decreased abruptly by sediment trapping in piedmont lakes related to deglaciation, resulting in a sudden decrease of sediment supply to the ocean. Additionally, reduced moisture supply may have contributed to the rapid decline of turbidite deposition. These different causes result in similar depositional patterns in the marine sinks. We conclude that turbiditic strata may constitute reliable recorders of climate change across a wide range of climatic zones and geomorphic conditions. However, the underlying causes for similar signal manifestations in the sinks may differ, ranging from maintained high system connectivity to abrupt connectivity loss. (C) 2017 Elsevier B.V. All rights reserved.}, language = {en} } @article{BloetheKorupSchwanghart2015, author = {Bl{\"o}the, Jan Henrik and Korup, Oliver and Schwanghart, Wolfgang}, title = {Large landslides lie low: Excess topography in the Himalaya-Karakoram ranges}, series = {Geology}, volume = {43}, journal = {Geology}, number = {6}, publisher = {American Institute of Physics}, address = {Boulder}, issn = {0091-7613}, doi = {10.1130/G36527.1}, pages = {523 -- 526}, year = {2015}, abstract = {Mass wasting is an important process for denuding hillslopes and lowering ridge crests in active mountain belts such as the Himalaya-Karakoram ranges (HKR). Such a high-relief landscape is likely to be at its mechanical threshold, maintained by competing rapid rock uplift, river incision, and pervasive slope failure. We introduce excess topography, Z(E), for quantifying potentially unstable rock-mass volumes inclined at angles greater than a specified threshold angle. We find that Z(E) peaks along major fluvial and glacial inner gorges, which is also where the majority of 492 large (>0.1 km(2)) rock-slope failures occur in the Himalaya's largest cluster of documented Pleistocene to Holocene bedrock landslides. Our data reveal that bedrock landslides in the HKR chiefly detached from near or below the median elevation, whereas glaciers and rock glaciers occupy higher-elevation bands almost exclusively. Less than 10\% of the area of the HKR is upslope of glaciers, such that possible censoring of evidence of large bedrock landslides above the permanent snow line barely affects this finding. Bedrock landslides appear to preferentially undermine topographic relief in response to fluvial and glacial incision along inner gorges, unless more frequent and smaller undetected failures, or rigorous (peri-)glacial erosion, compensate for this role at higher elevation. Either way, the distinct patterns of excess topography and large bedrock landsliding in the HKR juxtapose two stacked domains of landslide and (peri-)glacial erosion that may respond to different time scales of perturbation. Our findings call for more detailed analysis of vertical erosional domains and their geomorphic coupling in active mountain belts.}, language = {en} } @article{BrickerSchwanghartAdhikarietal.2017, author = {Bricker, Jeremy D. and Schwanghart, Wolfgang and Adhikari, Basanta Raj and Moriguchi, Shuji and Roeber, Volker and Giri, Sanjay}, title = {Performance of Models for Flash Flood Warning and Hazard Assessment}, series = {Mountain research and development}, volume = {37}, journal = {Mountain research and development}, number = {1}, publisher = {International Mountain Society}, address = {Lawrence}, issn = {0276-4741}, doi = {10.1659/MRD-JOURNAL-D-16-00043.1}, pages = {5 -- 15}, year = {2017}, abstract = {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.}, language = {en} } @article{CampfortsSchwanghartGovers2017, author = {Campforts, Benjamin and Schwanghart, Wolfgang and Govers, Gerard}, title = {Accurate simulation of transient landscape evolution by eliminating numerical diffusion}, series = {Earth surface dynamics}, volume = {5}, journal = {Earth surface dynamics}, number = {1}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {2196-6311}, doi = {10.5194/esurf-5-47-2017}, pages = {47 -- 66}, year = {2017}, abstract = {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.}, language = {en} } @misc{CampfortsSchwanghartGovers2017, author = {Campforts, Benjamin and Schwanghart, Wolfgang and Govers, Gerard}, title = {Accurate simulation of transient landscape evolution by eliminating numerical diffusion}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {664}, issn = {1866-8372}, doi = {10.25932/publishup-41878}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-418784}, pages = {20}, year = {2017}, abstract = {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.}, language = {en} } @misc{ClasonMairNienowetal.2015, author = {Clason, Caroline C. and Mair, D. W. F. and Nienow, P. W. and Bartholomew, I. D. and Sole, Andrew and Palmer, Steven and Schwanghart, Wolfgang}, title = {Modelling the transfer of supraglacial meltwater to the bed of Leverett Glacier, Southwest Greenland}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {513}, issn = {1866-8372}, doi = {10.25932/publishup-40905}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-409053}, pages = {16}, year = {2015}, abstract = {Meltwater delivered to the bed of the Greenland Ice Sheet is a driver of variable ice-motion through changes in effective pressure and enhanced basal lubrication. Ice surface velocities have been shown to respond rapidly both to meltwater production at the surface and to drainage of supraglacial lakes, suggesting efficient transfer of meltwater from the supraglacial to subglacial hydrological systems. Although considerable effort is currently being directed towards improved modelling of the controlling surface and basal processes, modelling the temporal and spatial evolution of the transfer of melt to the bed has received less attention. Here we present the results of spatially distributed modelling for prediction of moulins and lake drainages on the Leverett Glacier in Southwest Greenland. The model is run for the 2009 and 2010 ablation seasons, and for future increased melt scenarios. The temporal pattern of modelled lake drainages are qualitatively comparable with those documented from analyses of repeat satellite imagery. The modelled timings and locations of delivery of meltwater to the bed also match well with observed temporal and spatial patterns of ice surface speed-ups. This is particularly true for the lower catchment (< 1000 m a.s.l.) where both the model and observations indicate that the development of moulins is the main mechanism for the transfer of surface meltwater to the bed. At higher elevations (e.g. 1250-1500 m a.s.l.) the development and drainage of supraglacial lakes becomes increasingly important. At these higher elevations, the delay between modelled melt generation and subsequent delivery of melt to the bed matches the observed delay between the peak air temperatures and subsequent velocity speed-ups, while the instantaneous transfer of melt to the bed in a control simulation does not. Although both moulins and lake drainages are predicted to increase in number for future warmer climate scenarios, the lake drainages play an increasingly important role in both expanding the area over which melt accesses the bed and in enabling a greater proportion of surface melt to reach the bed.}, language = {en} } @article{ClasonMairNienowetal.2015, author = {Clason, Caroline C. and Mair, D. W. F. and Nienow, P. W. and Bartholomew, I. D. and Sole, Andrew and Palmer, Steven and Schwanghart, Wolfgang}, title = {Modelling the transfer of supraglacial meltwater to the bed of Leverett Glacier, Southwest Greenland}, series = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, volume = {9}, journal = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, number = {1}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1994-0416}, doi = {10.5194/tc-9-123-2015}, pages = {123 -- 138}, year = {2015}, abstract = {Meltwater delivered to the bed of the Greenland Ice Sheet is a driver of variable ice-motion through changes in effective pressure and enhanced basal lubrication. Ice surface velocities have been shown to respond rapidly both to meltwater production at the surface and to drainage of supraglacial lakes, suggesting efficient transfer of meltwater from the supraglacial to subglacial hydrological systems. Although considerable effort is currently being directed towards improved modelling of the controlling surface and basal processes, modelling the temporal and spatial evolution of the transfer of melt to the bed has received less attention. Here we present the results of spatially distributed modelling for prediction of moulins and lake drainages on the Leverett Glacier in Southwest Greenland. The model is run for the 2009 and 2010 ablation seasons, and for future increased melt scenarios. The temporal pattern of modelled lake drainages are qualitatively comparable with those documented from analyses of repeat satellite imagery. The modelled timings and locations of delivery of meltwater to the bed also match well with observed temporal and spatial patterns of ice surface speed-ups. This is particularly true for the lower catchment (< 1000 m a.s.l.) where both the model and observations indicate that the development of moulins is the main mechanism for the transfer of surface meltwater to the bed. At higher elevations (e.g. 1250-1500 m a.s.l.) the development and drainage of supraglacial lakes becomes increasingly important. At these higher elevations, the delay between modelled melt generation and subsequent delivery of melt to the bed matches the observed delay between the peak air temperatures and subsequent velocity speed-ups, while the instantaneous transfer of melt to the bed in a control simulation does not. Although both moulins and lake drainages are predicted to increase in number for future warmer climate scenarios, the lake drainages play an increasingly important role in both expanding the area over which melt accesses the bed and in enabling a greater proportion of surface melt to reach the bed.}, language = {en} } @article{GreeneThirumalaiKearneyetal.2019, author = {Greene, Chad A. and Thirumalai, Kaustubh and Kearney, Kelly A. and Delgado, Jose Miguel Martins and Schwanghart, Wolfgang and Wolfenbarger, Natalie S. and Thyng, Kristen M. and Gwyther, David E. and Gardner, Alex S. and Blankenship, Donald D.}, title = {The Climate Data Toolbox for MATLAB}, series = {Geochemistry, geophysics, geosystems}, volume = {20}, journal = {Geochemistry, geophysics, geosystems}, number = {7}, publisher = {American Geophysical Union}, address = {Washington}, issn = {1525-2027}, doi = {10.1029/2019GC008392}, pages = {3774 -- 3781}, year = {2019}, abstract = {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.}, language = {en} }