@article{GuentherSchueleZurelletal.2023, author = {G{\"u}nther, Oliver and Sch{\"u}le, Manja and Zurell, Damaris and Jeltsch, Florian and Roeleke, Manuel and Kampe, Heike and Zimmermann, Matthias and Scholz, Jana and Engbert, Ralf and Elsner, Birgit and Schlangen, David and Agrofylax, Luisa and Georgi, Doreen and Weymar, Mathias and Wagener, Thorsten and Bookhagen, Bodo and Eibl, Eva P. S. and Korup, Oliver and Oswald, Sascha and Thieken, Annegret and van der Beek, Peter}, title = {Portal Wissen = Exzellenz}, series = {Portal Wissen: Das Forschungsmagazin der Universit{\"a}t Potsdam}, journal = {Portal Wissen: Das Forschungsmagazin der Universit{\"a}t Potsdam}, number = {02/2023}, issn = {2194-4245}, doi = {10.25932/publishup-61144}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-611440}, pages = {98}, year = {2023}, abstract = {Was nicht nur gut oder sehr gut ist, nennen wir gern exzellent. Aber was meint das eigentlich? Vom lateinischen „excellere" kommend, beschreibt es Dinge, Personen oder Handlungen, die „hervor-" oder „herausragen" aus der Menge, sich „auszeichnen" gegen{\"u}ber anderen. Mehr geht nicht. Exzellenz ist das Mittel der Wahl, wenn es darum geht, der Erste oder Beste zu sein. Und das macht auch vor der Forschung nicht halt. Wer auf die Universit{\"a}t Potsdam schaut, findet zahlreiche ausgezeichnete Forschende, hervorragende Projekte und immer wieder auch aufsehenerregende Erkenntnisse, Ver{\"o}ffentlichungen und Ergebnisse. Aber ist die UP auch exzellent? Eine Frage, die 2023 ganz sicher andere Wellen schl{\"a}gt als vielleicht vor 20 Jahren. Denn seit dem Start der Exzellenzinitiative 2005 gelten als - w{\"o}rtlich - exzellent jene Hochschulen, denen es gelingt, in dem umfangreichsten F{\"o}rderprogramm f{\"u}r Wissenschaft in Deutschland einen Zuschlag zu erhalten. Egal ob in Form von Graduiertenschulen, Forschungsclustern oder - seit Fortsetzung des Programms ab 2019 unter dem Titel „Exzellenzstrategie" - ganzen Exzellenzuniversit{\"a}ten: Wer im Kreis der Forschungsuniversit{\"a}ten zu den Besten geh{\"o}ren will, braucht das Siegel der Exzellenz. In der gerade eingel{\"a}uteten neuen Wettbewerbsrunde der „Exzellenzstrategie des Bundes und der L{\"a}nder" bewirbt sich die Universit{\"a}t Potsdam mit drei Clusterskizzen um F{\"o}rderung. Ein Antrag kommt aus der {\"O}kologie- und Biodiversit{\"a}tsforschung. Ziel ist es, ein komplexes Bild {\"o}kologischer Prozesse zu zeichnen - und dabei die Rolle von einzelnen Individuen ebenso zu betrachten wie das Zusammenwirken vieler Arten in einem {\"O}kosystem, um die Funktion der Artenvielfalt genauer zu bestimmen. Eine zweite Skizze haben die Kognitionswissenschaften eingereicht. Hier soll das komplexe Nebeneinander von Sprache und Kognition, Entwicklung und Lernen sowie Motivation und Verhalten als dynamisches Miteinander erforscht werden - wobei auch mit den Erziehungswissenschaften kooperiert wird, um verkn{\"u}pfte Lernund Bildungsprozesse stets mitzudenken. Der dritte Antrag aus den Geo- und Umweltwissenschaften nimmt extreme und besonders folgenschwere Naturgefahren und -prozesse wie {\"U}berschwemmungen und D{\"u}rren in den Blick. Die Forschenden untersuchen die Extremereignisse mit besonderem Fokus auf deren Wechselwirkung mit der Gesellschaft, um mit ihnen einhergehende Risiken und Sch{\"a}den besser einsch{\"a}tzen sowie k{\"u}nftig rechtzeitig Maßnahmen einleiten zu k{\"o}nnen. „Alle drei Antr{\"a}ge zeichnen ein hervorragendes Bild unserer Leistungsf{\"a}higkeit", betont der Pr{\"a}sident der Universit{\"a}t, Prof. Oliver G{\"u}nther, Ph.D. „Die Skizzen dokumentieren eindrucksvoll unser Engagement, vorhandene Forschungsexzellenz sowie die Potenziale der Universit{\"a}t Potsdam insgesamt. Allein die Tatsache, dass sich drei schlagkr{\"a}ftige Konsortien in ganz unterschiedlichen Themenbereichen zusammengefunden haben, zeigt, dass wir auf unserem Weg in die Spitzengruppe der deutschen Universit{\"a}ten einen guten Schritt vorangekommen sind." In diesem Heft schauen wir, was sich in und hinter diesen Antr{\"a}gen verbirgt: Wir haben mit den Wissenschaftlerinnen und Wissenschaftlern gesprochen, die sie geschrieben haben, und sie gefragt, was sie sich vornehmen, sollten sie den Zuschlag erhalten und ein Cluster an die Universit{\"a}t holen. Wir haben aber auch auf die Forschung geschaut, die zu den Antr{\"a}gen gef{\"u}hrt hat und die schon l{\"a}nger das Profil der Universit{\"a}t pr{\"a}gt und ihr national wie international Anerkennung eingebracht hat. Wir stellen eine kleine Auswahl an Projekten, Methoden und Forschenden vor, um zu zeigen, warum in diesen Antr{\"a}gen tats{\"a}chlich exzellente Forschung steckt! {\"U}brigens: Auch „Exzellenz" ist nicht das Ende der Fahnenstange. Immerhin l{\"a}sst sich das Adjektiv exzellent sogar steigern. In diesem Sinne w{\"u}nschen wir exzellentestes Vergn{\"u}gen beim Lesen!}, language = {de} } @article{GuentherSchueleZurelletal.2023, author = {G{\"u}nther, Oliver and Sch{\"u}le, Manja and Zurell, Damaris and Jeltsch, Florian and Roeleke, Manuel and Kampe, Heike and Zimmermann, Matthias and Scholz, Jana and Mikulla, Stefanie and Engbert, Ralf and Elsner, Birgit and Schlangen, David and Agrofylax, Luisa and Georgi, Doreen and Weymar, Mathias and Wagener, Thorsten and Bookhagen, Bodo and Eibl, Eva P. S. and Korup, Oliver and Oswald, Sascha and Thieken, Annegret and van der Beek, Peter}, title = {Portal Wissen = Excellence}, series = {Portal Wissen: The research magazine of the University of Potsdam}, journal = {Portal Wissen: The research magazine of the University of Potsdam}, number = {02/2023}, issn = {2198-9974}, doi = {10.25932/publishup-61145}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-611456}, pages = {58}, year = {2023}, abstract = {When something is not just good or very good, we often call it excellent. But what does that really mean? Coming from the Latin word "excellere," it describes things, persons, or actions that are outstanding or superior and distinguish themselves from others. It cannot get any better. Excellence is the top choice for being the first or the best. Research is no exception. At the university, you will find numerous exceptional researchers, outstanding projects, and, time and again, sensational findings, publications, and results. But is the University of Potsdam also excellent? A question that will certainly create a different stir in 2023 than it did perhaps 20 years ago. Since the launch of the Excellence Initiative in 2005, universities that succeed in winning the most comprehensive funding program for research in Germany have been considered - literally - excellent. Whether in the form of graduate schools, research clusters, or - since the program was continued in 2019 under the title "Excellence Strategy" - entire universities of excellence: Anyone who wants to be among the best research universities needs the seal of excellence. The University of Potsdam is applying for funding with three cluster proposals in the recently launched new round of the "Excellence Strategy of the German Federal and State Governments." One proposal comes from ecology and biodiversity research. The aim is to paint a comprehensive picture of ecological processes by examining the role of single individuals as well as the interactions among many species in an ecosystem to precisely determine the function of biodiversity. A second proposal has been submitted by the cognitive sciences. Here, the complex coexistence of language and cognition, development and learning, as well as motivation and behavior will be researched as a dynamic interrelation. The projects will include cooperation with the educational sciences to constantly consider linked learning and educational processes. The third proposal from the geo and environmental sciences concentrates on extreme and particularly devastating natural hazards and processes such as floods and droughts. The researchers examine these extreme events, focusing on their interaction with society, to be able to better assess the risks and damages they might involve and to initiate timely measures in the future. "All three proposals highlight the excellence of our performance," emphasizes University President Prof. Oliver G{\"u}nther, Ph.D. "The outlines impressively document our commitment, existing research excellence, and the potential of the University of Potsdam as a whole. The fact that three powerful consortia have come together in different subject areas shows that we have taken a good step forward on our way to becoming one of the top German universities." In this issue, we are looking at what is in and behind these proposals: We talked to the researchers who wrote them. We asked them about their plans in case their proposals are successful and they bring a cluster of excellence to the university. But we also looked at the research that has led to the proposals, has long shaped the university's profile, and earned it national and international recognition. We present a small selection of projects, methods, and researchers to illustrate why there really is excellent research in these proposals! By the way, "excellence" is also not the end of the flagpole. After all, the adjective "excellent" even has a comparative and a superlative. With this in mind, I wish you the most excellent pleasure reading this issue!}, language = {en} } @misc{AtmaniBookhagenSmith2022, author = {Atmani, Farid and Bookhagen, Bodo and Smith, Taylor}, title = {Measuring Vegetation Heights and Their Seasonal Changes in the Western Namibian Savanna Using Spaceborne Lidars}, series = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {1275}, issn = {1866-8372}, doi = {10.25932/publishup-56991}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-569915}, pages = {20}, year = {2022}, abstract = {The Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) with its land and vegetation height data product (ATL08), and Global Ecosystem Dynamics Investigation (GEDI) with its terrain elevation and height metrics data product (GEDI Level 2A) missions have great potential to globally map ground and canopy heights. Canopy height is a key factor in estimating above-ground biomass and its seasonal changes; these satellite missions can also improve estimated above-ground carbon stocks. This study presents a novel Sparse Vegetation Detection Algorithm (SVDA) which uses ICESat-2 (ATL03, geolocated photons) data to map tree and vegetation heights in a sparsely vegetated savanna ecosystem. The SVDA consists of three main steps: First, noise photons are filtered using the signal confidence flag from ATL03 data and local point statistics. Second, we classify ground photons based on photon height percentiles. Third, tree and grass photons are classified based on the number of neighbors. We validated tree heights with field measurements (n = 55), finding a root-mean-square error (RMSE) of 1.82 m using SVDA, GEDI Level 2A (Geolocated Elevation and Height Metrics product): 1.33 m, and ATL08: 5.59 m. Our results indicate that the SVDA is effective in identifying canopy photons in savanna ecosystems, where ATL08 performs poorly. We further identify seasonal vegetation height changes with an emphasis on vegetation below 3 m; widespread height changes in this class from two wet-dry cycles show maximum seasonal changes of 1 m, possibly related to seasonal grass-height differences. Our study shows the difficulties of vegetation measurements in savanna ecosystems but provides the first estimates of seasonal biomass changes.}, language = {en} } @article{AtmaniBookhagenSmith2022, author = {Atmani, Farid and Bookhagen, Bodo and Smith, Taylor}, title = {Measuring vegetation heights and their seasonal changes in the Western Namibian Savanna using spaceborne lidars}, series = {Remote sensing / Molecular Diversity Preservation International (MDPI)}, volume = {14}, journal = {Remote sensing / Molecular Diversity Preservation International (MDPI)}, number = {12}, edition = {12}, publisher = {MDPI}, address = {Basel, Schweiz}, issn = {2072-4292}, doi = {10.3390/rs14122928}, pages = {1 -- 20}, year = {2022}, abstract = {The Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) with its land and vegetation height data product (ATL08), and Global Ecosystem Dynamics Investigation (GEDI) with its terrain elevation and height metrics data product (GEDI Level 2A) missions have great potential to globally map ground and canopy heights. Canopy height is a key factor in estimating above-ground biomass and its seasonal changes; these satellite missions can also improve estimated above-ground carbon stocks. This study presents a novel Sparse Vegetation Detection Algorithm (SVDA) which uses ICESat-2 (ATL03, geolocated photons) data to map tree and vegetation heights in a sparsely vegetated savanna ecosystem. The SVDA consists of three main steps: First, noise photons are filtered using the signal confidence flag from ATL03 data and local point statistics. Second, we classify ground photons based on photon height percentiles. Third, tree and grass photons are classified based on the number of neighbors. We validated tree heights with field measurements (n = 55), finding a root-mean-square error (RMSE) of 1.82 m using SVDA, GEDI Level 2A (Geolocated Elevation and Height Metrics product): 1.33 m, and ATL08: 5.59 m. Our results indicate that the SVDA is effective in identifying canopy photons in savanna ecosystems, where ATL08 performs poorly. We further identify seasonal vegetation height changes with an emphasis on vegetation below 3 m; widespread height changes in this class from two wet-dry cycles show maximum seasonal changes of 1 m, possibly related to seasonal grass-height differences. Our study shows the difficulties of vegetation measurements in savanna ecosystems but provides the first estimates of seasonal biomass changes.}, language = {en} } @article{RegmiBookhagen2022, author = {Regmi, Shakil and Bookhagen, Bodo}, title = {The spatial pattern of extreme precipitation from 40 years of gauge data in the central Himalaya}, series = {Weather and climate extremes}, volume = {37}, journal = {Weather and climate extremes}, publisher = {Elsevier}, address = {Amsterdam}, issn = {2212-0947}, doi = {10.1016/j.wace.2022.100470}, pages = {14}, year = {2022}, abstract = {The topography of the Himalaya exerts a substantial control on the spatial distribution of monsoonal rainfall, which is a vital water source for the regional economy and population. But the occurrence of short-lived and high-intensity precipitation results in socio-economic losses. This study relies on 40 years of daily data from 204 ground stations in Nepal to derive extreme precipitation thresholds, amounts, and days at the 95th percentile. We additionally determine the precipitation magnitude-frequency relation. We observe that extreme precipitation amounts follow an almost uniform band parallel to topographic contour lines in the southern Himalaya mountains in central and eastern Nepal but not in western Nepal. The relationship of extreme precipitation indices with topographic relief shows that extreme precipitation thresholds decrease with increasing elevation, but extreme precipitation days increase in higher elevation areas. Furthermore, stations above 1 km elevation exhibit a power-law relation in the rainfall magnitude-frequency framework. Stations at higher elevations generally have lower values of power-law exponents than low elevation areas. This suggests a fundamentally different behaviour of the rainfall distribution and an increased occurrence of extreme rainfall storms in the high elevation areas of Nepal.}, language = {en} } @article{ErbelloDoelessoMelnickZeilingeretal.2022, author = {Erbello Doelesso, Asfaw and Melnick, Daniel and Zeilinger, Gerold and Bookhagen, Bodo and Pingel, Heiko and Strecker, Manfred}, title = {Geomorphic expression of a tectonically active rift-transfer zone in southern Ethiopia}, series = {Geomorphology : an international journal on pure and applied geomorphology}, volume = {403}, journal = {Geomorphology : an international journal on pure and applied geomorphology}, publisher = {Elsevier Science}, address = {Amsterdam [u.a.]}, issn = {0169-555X}, doi = {10.1016/j.geomorph.2022.108162}, pages = {20}, year = {2022}, abstract = {The Gofa Province and the Chew Bahir Basin of southern Ethiopia constitute tectonically active regions, where the Southern Main Ethiopian Rift converges with the Northern Kenya Rift through a wide zone of extensional deformation with several north to northeast-trending, left-stepping en-e \& PRIME;chelon basins. This sector of the Southern Main Ethiopian Rift is characterized by a semi-arid climate and a largely uniform lithology, and thus provides ideal conditions for studying the different parameters that define the tectonic and geomorphic features of this complex kinematic transfer zone. In this study, the degree of tectonic activity, spatiotemporal variations in extension, and the nature of kinematic linkage between different fault systems of the transfer zone are constrained by detailed quantitative geomorphic analysis of river catchments and focused field work. We analyzed fluvial and landscape morphometric characteristics in combination with structural, seismicity, and climatic data to better evaluate the tectono-geomorphic history of this transfer zone. Our data reveal significant north-south variations in the degree of extension from the Sawula Basin in the north (mature) to the Chew Bahir Basin in the south (juvenile). First, normalized channel-steepness indices and the spatial arrangement of knickpoints in footwall-draining streams suggest a gradual, southward shift in extensional deformation and recent tectonic activity. Second, based on 1-k(m) radius local relief and mean-hillslope maximum values that are consistent with ksn anomalies, we confirm strain localization within zones of fault interaction. Third, morphometric indices such as hypsometry, basin asymmetry factor, and valley floor width to valley height ratio also indicate a north to south gradient in tectonic activity, highlighting the importance of such a wide transfer zone with diffuse extension linking different rift segments during the break-up of continental crust.}, language = {en} } @article{PenaMetzgerHeidbachetal.2022, author = {Pe{\~n}a, Carlos and Metzger, Sabrina and Heidbach, Oliver and Bedford, Jonathan and Bookhagen, Bodo and Moreno, Marcos and Oncken, Onno and Cotton, Fabrice}, title = {Role of poroelasticity during the early postseismic deformation of the 2010 Maule megathrust earthquake}, series = {Geophysical research letters}, volume = {49}, journal = {Geophysical research letters}, number = {9}, publisher = {Wiley}, address = {Hoboken, NJ}, issn = {0094-8276}, doi = {10.1029/2022GL098144}, pages = {11}, year = {2022}, abstract = {Megathrust earthquakes impose changes of differential stress and pore pressure in the lithosphere-asthenosphere system that are transiently relaxed during the postseismic period primarily due to afterslip, viscoelastic and poroelastic processes. Especially during the early postseismic phase, however, the relative contribution of these processes to the observed surface deformation is unclear. To investigate this, we use geodetic data collected in the first 48 days following the 2010 Maule earthquake and a poro-viscoelastic forward model combined with an afterslip inversion. This model approach fits the geodetic data 14\% better than a pure elastic model. Particularly near the region of maximum coseismic slip, the predicted surface poroelastic uplift pattern explains well the observations. If poroelasticity is neglected, the spatial afterslip distribution is locally altered by up to +/- 40\%. Moreover, we find that shallow crustal aftershocks mostly occur in regions of increased postseismic pore-pressure changes, indicating that both processes might be mechanically coupled.}, language = {en} } @article{PurintonBookhagen2021, author = {Purinton, Benjamin and Bookhagen, Bodo}, title = {Tracking downstream variability in large grain-size distributions in the South-Central Andes}, series = {Journal of geophysical research : F, Earth surface}, volume = {126}, journal = {Journal of geophysical research : F, Earth surface}, number = {8}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9003}, doi = {10.1029/2021JF006260}, pages = {1 -- 29}, year = {2021}, abstract = {Mixed sand- and gravel-bed rivers record erosion, transport, and fining signals in their bedload size distributions. Thus, grain-size data are imperative for studying these processes. However, collecting hundreds to thousands of pebble measurements in steep and dynamic high-mountain river settings remains challenging. Using the recently published digital grain-sizing algorithm PebbleCounts, we were able to survey seven large (>= 1,000 m2) channel cross-sections and measure thousands to tens-of-thousands of grains per survey along a 100-km stretch of the trunk stream of the Toro Basin in Northwest Argentina. The study region traverses a steep topographic and environmental gradient on the eastern margin of the Central Andean Plateau. Careful counting and validation allows us to identify measurement errors and constrain percentile uncertainties using large sample sizes. In the coarse >= 2.5 cm fraction of bedload, only the uppermost size percentiles (>= 95th) vary significantly downstream, whereas the 50th and 84th percentiles show less variability. We note a relation between increases in these upper percentiles and along-channel junctions with large, steep tributaries. This signal is strongly influenced by lithology and geologic structures, and mixed with local hillslope input. In steep catchments like the Toro Basin, we suggest nonlinear relationships between geomorphic metrics and grain size, whereby the steepest parts of the landscape exert primary control on the upper grain-size percentiles. Thus, average or median metrics that do not apply weights or thresholds to steeper topography may be less predictive of grain-size distributions in such settings.}, language = {en} } @article{SmithRheinwaltBookhagen2021, author = {Smith, Taylor and Rheinwalt, Aljoscha and Bookhagen, Bodo}, title = {Topography and climate in the upper Indus Basin}, series = {The science of the total environment : an international journal for scientific research into the environment and its relationship with man}, volume = {786}, journal = {The science of the total environment : an international journal for scientific research into the environment and its relationship with man}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0048-9697}, doi = {10.1016/j.scitotenv.2021.147363}, pages = {11}, year = {2021}, abstract = {The Upper Indus Basin (UIB), which covers a wide range of climatic and topographic settings, provides an ideal venue to explore the relationship between climate and topography. While the distribution of snow and glaciers is spatially and temporally heterogeneous, there exist regions with similar elevation-snow relationships. In this work, we construct elevation-binned snow-cover statistics to analyze 3415 watersheds and 7357 glaciers in the UIB region. We group both glaciers and watersheds using a hierarchical clustering approach and find that (1) watershed clusters mirror large-scale moisture transport patterns and (2) are highly dependent on median watershed elevation. (3) Glacier clusters are spatially heterogeneous and are less strongly controlled by elevation, but rather by local topographic parameters that modify solar insolation. Our clustering approach allows us to clearly define self-similar snow-topographic regions. Eastern watersheds in the UIB show a steep snow cover-elevation relationship whereas watersheds in the central and western UIB have moderately sloped relationships, but cluster in distinct groups. We highlight this snow-cover-topographic transition zone and argue that these watersheds have different hydrologic responses than other regions. Our hierarchical clustering approach provides a potential new framework to use in defining climatic zones in the cyrosphere based on empirical data.}, language = {en} } @article{SmithBookhagen2021, author = {Smith, Taylor and Bookhagen, Bodo}, title = {Climatic and biotic controls on topographic asymmetry at the global scale}, series = {Journal of geophysical research : JGR, Earth surface}, volume = {126}, journal = {Journal of geophysical research : JGR, Earth surface}, number = {1}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9003}, doi = {10.1029/2020JF005692}, pages = {24}, year = {2021}, abstract = {Insolation differences play a primary role in controlling microclimate and vegetation cover, which together influence the development of topography. Topographic asymmetry (TA), or slope differences between terrain aspects, has been well documented in small-scale, field-based, and modeling studies. Here we combine a suite of environmental (e.g., vegetation, temperature, solar insolation) and topographic (e.g., elevation, drainage network) data to explore the driving mechanisms and markers of TA on a global scale. Using a novel empirical TA analysis method, we find that (1) steeper terrain has higher TA magnitudes, (2) globally, pole-facing terrain is on average steeper than equator-facing terrain, especially in mid-latitude, tectonically quiescent, and vegetated landscapes, and (3) high-elevation and low-temperature regions tend to have terrain steepened toward the equator. We further show that there are distinct differences in climate and vegetation cover across terrain aspects, and that TA is reflected in the size and form of fluvial drainage networks. Our work supports the argument that insolation asymmetries engender differences in local microclimates and vegetation on opposing terrain aspects, which broadly encourage the development of asymmetric topography across a range of lithologic, tectonic, geomorphic, and climatic settings.}, language = {en} } @article{PurintonBookhagen2021, author = {Purinton, Benjamin and Bookhagen, Bodo}, title = {Beyond Vertical Point Accuracy}, series = {Frontiers in Earth Science}, journal = {Frontiers in Earth Science}, publisher = {Frontiers Media}, address = {Lausanne, Schweiz}, issn = {2296-6463}, doi = {10.3389/feart.2021.758606}, pages = {1 -- 24}, year = {2021}, abstract = {Quantitative geomorphic research depends on accurate topographic data often collected via remote sensing. Lidar, and photogrammetric methods like structure-from-motion, provide the highest quality data for generating digital elevation models (DEMs). Unfortunately, these data are restricted to relatively small areas, and may be expensive or time-consuming to collect. Global and near-global DEMs with 1 arcsec (∼30 m) ground sampling from spaceborne radar and optical sensors offer an alternative gridded, continuous surface at the cost of resolution and accuracy. Accuracy is typically defined with respect to external datasets, often, but not always, in the form of point or profile measurements from sources like differential Global Navigation Satellite System (GNSS), spaceborne lidar (e.g., ICESat), and other geodetic measurements. Vertical point or profile accuracy metrics can miss the pixel-to-pixel variability (sometimes called DEM noise) that is unrelated to true topographic signal, but rather sensor-, orbital-, and/or processing-related artifacts. This is most concerning in selecting a DEM for geomorphic analysis, as this variability can affect derivatives of elevation (e.g., slope and curvature) and impact flow routing. We use (near) global DEMs at 1 arcsec resolution (SRTM, ASTER, ALOS, TanDEM-X, and the recently released Copernicus) and develop new internal accuracy metrics to assess inter-pixel variability without reference data. Our study area is in the arid, steep Central Andes, and is nearly vegetation-free, creating ideal conditions for remote sensing of the bare-earth surface. We use a novel hillshade-filtering approach to detrend long-wavelength topographic signals and accentuate short-wavelength variability. Fourier transformations of the spatial signal to the frequency domain allows us to quantify: 1) artifacts in the un-projected 1 arcsec DEMs at wavelengths greater than the Nyquist (twice the nominal resolution, so > 2 arcsec); and 2) the relative variance of adjacent pixels in DEMs resampled to 30-m resolution (UTM projected). We translate results into their impact on hillslope and channel slope calculations, and we highlight the quality of the five DEMs. We find that the Copernicus DEM, which is based on a carefully edited commercial version of the TanDEM-X, provides the highest quality landscape representation, and should become the preferred DEM for topographic analysis in areas without sufficient coverage of higher-quality local DEMs.}, language = {en} } @article{RamezaniZiaraniBookhagenSchmidtetal.2021, author = {Ramezani Ziarani, Maryam and Bookhagen, Bodo and Schmidt, Torsten and Wickert, Jens and de la Torre, Alejandro and Deng, Zhiguo and Calori, Andrea}, title = {A model for the relationship between rainfall, GNSS-derived integrated water vapour, and CAPE in the eastern central Andes}, series = {Remote Sensing}, volume = {13}, journal = {Remote Sensing}, number = {18}, publisher = {MDPI}, address = {Basel}, issn = {2072-4292}, doi = {10.3390/rs13183788}, pages = {19}, year = {2021}, abstract = {Atmospheric water vapour content is a key variable that controls the development of deep convective storms and rainfall extremes over the central Andes. Direct measurements of water vapour are challenging; however, recent developments in microwave processing allow the use of phase delays from L-band radar to measure the water vapour content throughout the atmosphere: Global Navigation Satellite System (GNSS)-based integrated water vapour (IWV) monitoring shows promising results to measure vertically integrated water vapour at high temporal resolutions. Previous works also identified convective available potential energy (CAPE) as a key climatic variable for the formation of deep convective storms and rainfall in the central Andes. Our analysis relies on GNSS data from the Argentine Continuous Satellite Monitoring Network, Red Argentina de Monitoreo Satelital Continuo (RAMSAC) network from 1999 to 2013. CAPE is derived from version 2.0 of the ECMWF's (European Centre for Medium-Range Weather Forecasts) Re-Analysis (ERA-interim) and rainfall from the TRMM (Tropical Rainfall Measuring Mission) product. In this study, we first analyse the rainfall characteristics of two GNSS-IWV stations by comparing their complementary cumulative distribution function (CCDF). Second, we separately derive the relation between rainfall vs. CAPE and GNSS-IWV. Based on our distribution fitting analysis, we observe an exponential relation of rainfall to GNSS-IWV. In contrast, we report a power-law relationship between the daily mean value of rainfall and CAPE at the GNSS-IWV station locations in the eastern central Andes that is close to the theoretical relationship based on parcel theory. Third, we generate a joint regression model through a multivariable regression analysis using CAPE and GNSS-IWV to explain the contribution of both variables in the presence of each other to extreme rainfall during the austral summer season. We found that rainfall can be characterised with a higher statistical significance for higher rainfall quantiles, e.g., the 0.9 quantile based on goodness-of-fit criterion for quantile regression. We observed different contributions of CAPE and GNSS-IWV to rainfall for each station for the 0.9 quantile. Fourth, we identify the temporal relation between extreme rainfall (the 90th, 95th, and 99th percentiles) and both GNSS-IWV and CAPE at 6 h time steps. We observed an increase before the rainfall event and at the time of peak rainfall—both for GNSS-integrated water vapour and CAPE. We show higher values of CAPE and GNSS-IWV for higher rainfall percentiles (99th and 95th percentiles) compared to the 90th percentile at a 6-h temporal scale. Based on our correlation analyses and the dynamics of the time series, we show that both GNSS-IWV and CAPE had comparable magnitudes, and we argue to consider both climatic variables when investigating their effect on rainfall extremes.}, language = {en} } @misc{RamezaniZiaraniBookhagenSchmidtetal.2021, author = {Ramezani Ziarani, Maryam and Bookhagen, Bodo and Schmidt, Torsten and Wickert, Jens and de la Torre, Alejandro and Deng, Zhiguo and Calori, Andrea}, title = {A model for the relationship between rainfall, GNSS-derived integrated water vapour, and CAPE in the eastern central Andes}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {1172}, issn = {1866-8372}, doi = {10.25932/publishup-52325}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-523256}, pages = {21}, year = {2021}, abstract = {Atmospheric water vapour content is a key variable that controls the development of deep convective storms and rainfall extremes over the central Andes. Direct measurements of water vapour are challenging; however, recent developments in microwave processing allow the use of phase delays from L-band radar to measure the water vapour content throughout the atmosphere: Global Navigation Satellite System (GNSS)-based integrated water vapour (IWV) monitoring shows promising results to measure vertically integrated water vapour at high temporal resolutions. Previous works also identified convective available potential energy (CAPE) as a key climatic variable for the formation of deep convective storms and rainfall in the central Andes. Our analysis relies on GNSS data from the Argentine Continuous Satellite Monitoring Network, Red Argentina de Monitoreo Satelital Continuo (RAMSAC) network from 1999 to 2013. CAPE is derived from version 2.0 of the ECMWF's (European Centre for Medium-Range Weather Forecasts) Re-Analysis (ERA-interim) and rainfall from the TRMM (Tropical Rainfall Measuring Mission) product. In this study, we first analyse the rainfall characteristics of two GNSS-IWV stations by comparing their complementary cumulative distribution function (CCDF). Second, we separately derive the relation between rainfall vs. CAPE and GNSS-IWV. Based on our distribution fitting analysis, we observe an exponential relation of rainfall to GNSS-IWV. In contrast, we report a power-law relationship between the daily mean value of rainfall and CAPE at the GNSS-IWV station locations in the eastern central Andes that is close to the theoretical relationship based on parcel theory. Third, we generate a joint regression model through a multivariable regression analysis using CAPE and GNSS-IWV to explain the contribution of both variables in the presence of each other to extreme rainfall during the austral summer season. We found that rainfall can be characterised with a higher statistical significance for higher rainfall quantiles, e.g., the 0.9 quantile based on goodness-of-fit criterion for quantile regression. We observed different contributions of CAPE and GNSS-IWV to rainfall for each station for the 0.9 quantile. Fourth, we identify the temporal relation between extreme rainfall (the 90th, 95th, and 99th percentiles) and both GNSS-IWV and CAPE at 6 h time steps. We observed an increase before the rainfall event and at the time of peak rainfall—both for GNSS-integrated water vapour and CAPE. We show higher values of CAPE and GNSS-IWV for higher rainfall percentiles (99th and 95th percentiles) compared to the 90th percentile at a 6-h temporal scale. Based on our correlation analyses and the dynamics of the time series, we show that both GNSS-IWV and CAPE had comparable magnitudes, and we argue to consider both climatic variables when investigating their effect on rainfall extremes.}, language = {en} } @misc{PurintonBookhagen2021, author = {Purinton, Benjamin and Bookhagen, Bodo}, title = {Beyond Vertical Point Accuracy}, series = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, publisher = {Universit{\"a}tsverlag Potsdam}, address = {Potsdam}, issn = {1866-8372}, doi = {10.25932/publishup-54980}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-549805}, pages = {1 -- 24}, year = {2021}, abstract = {Quantitative geomorphic research depends on accurate topographic data often collected via remote sensing. Lidar, and photogrammetric methods like structure-from-motion, provide the highest quality data for generating digital elevation models (DEMs). Unfortunately, these data are restricted to relatively small areas, and may be expensive or time-consuming to collect. Global and near-global DEMs with 1 arcsec (∼30 m) ground sampling from spaceborne radar and optical sensors offer an alternative gridded, continuous surface at the cost of resolution and accuracy. Accuracy is typically defined with respect to external datasets, often, but not always, in the form of point or profile measurements from sources like differential Global Navigation Satellite System (GNSS), spaceborne lidar (e.g., ICESat), and other geodetic measurements. Vertical point or profile accuracy metrics can miss the pixel-to-pixel variability (sometimes called DEM noise) that is unrelated to true topographic signal, but rather sensor-, orbital-, and/or processing-related artifacts. This is most concerning in selecting a DEM for geomorphic analysis, as this variability can affect derivatives of elevation (e.g., slope and curvature) and impact flow routing. We use (near) global DEMs at 1 arcsec resolution (SRTM, ASTER, ALOS, TanDEM-X, and the recently released Copernicus) and develop new internal accuracy metrics to assess inter-pixel variability without reference data. Our study area is in the arid, steep Central Andes, and is nearly vegetation-free, creating ideal conditions for remote sensing of the bare-earth surface. We use a novel hillshade-filtering approach to detrend long-wavelength topographic signals and accentuate short-wavelength variability. Fourier transformations of the spatial signal to the frequency domain allows us to quantify: 1) artifacts in the un-projected 1 arcsec DEMs at wavelengths greater than the Nyquist (twice the nominal resolution, so > 2 arcsec); and 2) the relative variance of adjacent pixels in DEMs resampled to 30-m resolution (UTM projected). We translate results into their impact on hillslope and channel slope calculations, and we highlight the quality of the five DEMs. We find that the Copernicus DEM, which is based on a carefully edited commercial version of the TanDEM-X, provides the highest quality landscape representation, and should become the preferred DEM for topographic analysis in areas without sufficient coverage of higher-quality local DEMs.}, language = {en} } @article{MuetingBookhagenStrecker2021, author = {M{\"u}ting, Friederike Ariane and Bookhagen, Bodo and Strecker, Manfred}, title = {Identification of debris-flow channels using high-resolution topographic data}, series = {Journal of geophysical research : JGR, Earth surface}, volume = {126}, journal = {Journal of geophysical research : JGR, Earth surface}, number = {12}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9003}, doi = {10.1029/2021JF006330}, pages = {25}, year = {2021}, abstract = {Resolving Earth's surface at the meter scale is essential for an improved understanding of the dynamics of mass-movement processes. In this study, we explore the applicability and potential of digital elevation models (DEMs) derived from stereophotogrammetry to detect debris-flow channels in the Quebrada del Toro in the northwestern Argentine Andes. Our analysis relies on a high-resolution (3 m) DEM created from SPOT-7 tri-stereo satellite data. We carefully validated DEM quality with ∼6,000 differential GPS points and identified optimal parameters for DEM generation in high-relief terrain. After multiple processing steps, we achieved an accuracy of 0.051 ± 1.915 m (1σ) using n = 3,139 control points with cm precision. Previous studies have used the drainage area and slope framework to identify topographic signatures of debris flows within a catchment. We built upon this and investigated individual river-channel segments using connected-component (CC) analysis on meter-scale topographic data. We define CC as segments of similar slope along the channel profile. Based on seven manually identified debris-flow catchments, we developed a debris-flow similarity index using component length and mean channel-segment slope and identified channel segments that have likely been shaped by debris flows. The presented approach has the potential to resolve intra-catchment variability of transport processes, allows to constrain the extent of debris-flow channels more precisely than slope-area analysis, and highlights the versatility of combined space- and field-based observations for natural-hazard assessments.}, language = {en} } @article{PurintonBookhagen2020, author = {Purinton, Benjamin and Bookhagen, Bodo}, title = {Multiband (X, C, L) radar amplitude analysis for a mixed sand- and gravel-bed river in the eastern Central Andes}, series = {Remote sensing of environment : an interdisciplinary journal}, volume = {246}, journal = {Remote sensing of environment : an interdisciplinary journal}, publisher = {Elsevier}, address = {New York}, issn = {0034-4257}, doi = {10.1016/j.rse.2020.111799}, pages = {16}, year = {2020}, abstract = {Synthetic Aperture Radar (SAR) amplitude measurements from spaceborne sensors are sensitive to surface roughness conditions near their radar wavelength. These backscatter signals are often exploited to assess the roughness of plowed agricultural fields and water surfaces, and less so to complex, heterogeneous geological surfaces. The bedload of mixed sand- and gravel-bed rivers can be considered a mixture of smooth (compacted sand) and rough (gravel) surfaces. Here, we assess backscatter gradients over a large high-mountain alluvial river in the eastern Central Andes with aerially exposed sand and gravel bedload using X-band TerraSAR-X/TanDEM-X, C-band Sentinel-1, and L-band ALOS-2 PALSAR-2 radar scenes. In a first step, we present theory and hypotheses regarding radar response to an alluvial channel bed. We test our hypotheses by comparing backscatter responses over vegetation-free endmember surfaces from inside and outside of the active channel-bed area. We then develop methods to extract smoothed backscatter gradients downstream along the channel using kernel density estimates. In a final step, the local variability of sand-dominated patches is analyzed using Fourier frequency analysis, by fitting stretched-exponential and power-law regression models to the 2-D power spectrum of backscatter amplitude. We find a large range in backscatter depending on the heterogeneity of contiguous smooth- and rough-patches of bedload material. The SAR amplitude signal responds primarily to the fraction of smooth-sand bedload, but is further modified by gravel elements. The sensitivity to gravel is more apparent in longer wavelength L-band radar, whereas C- and X-band is sensitive only to sand variability. Because the spatial extent of smooth sand patches in our study area is typically< 50 m, only higher resolution sensors (e.g., TerraSAR-X/TanDEM-X) are useful for power spectrum analysis. Our results show the potential for mapping sand-gravel transitions and local geomorphic complexity in alluvial rivers with aerially exposed bedload using SAR amplitude.}, language = {en} } @article{OlenBookhagen2020, author = {Olen, Stephanie M. and Bookhagen, Bodo}, title = {Applications of SAR interferometric coherence time series}, series = {Journal of geophysical research : Earth surface}, volume = {125}, journal = {Journal of geophysical research : Earth surface}, number = {3}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9003}, doi = {10.1029/2019JF005141}, pages = {22}, year = {2020}, abstract = {Sediment transport domains in mountain landscapes are characterized by fundamentally different processes and rates depending on several factors, including geology, climate, and biota. Accurately identifying where transitions between transport domains occur is an important step to quantify the past, present, and future contribution of varying erosion and sedimentation processes and enhance our predictive capabilities. We propose a new methodology based on time series of synthetic aperture radar (SAR) interferometric coherence images to map sediment transport regimes across arid and semiarid landscapes. Using 4 years of Sentinel-1 data, we analyze sediment transport regimes for the south-central Andes in northwestern Argentina characterized by steep topographic and climatic gradients. We observe seasonally low coherence during the regional wet season, particularly on hillslopes and in alluvial channels. The spatial distribution of coherence is compared to drainage areas extracted from digital topography to identify two distinct transitions within watersheds: (a) a hillslope-to-fluvial and (b) a fluvial-to-alluvial transition. While transitions within a given basin can be well-constrained, the relative role of each sediment transport domain varies widely over the climatic and topographic gradients. In semiarid regions, we observe larger relative contributions from hillslopes compared to arid regions. Across regional gradients, the range of coherence within basins positively correlates to previously published millennial catchment-wide erosion rates and to topographic metrics used to indicate long-term uplift. Our study suggests that a dense time series of interferometric coherence can be used as a proxy for surface sediment movement and landscape stability in vegetation-free settings at event to decadal timescales.}, language = {en} } @misc{SmithBookhagen2020, author = {Smith, Taylor and Bookhagen, Bodo}, title = {Assessing Multi-Temporal Snow-Volume Trends in High Mountain Asia From 1987 to 2016 Using High-Resolution Passive Microwave Data}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {1020}, issn = {1866-8372}, doi = {10.25932/publishup-48417}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-484176}, pages = {15}, year = {2020}, abstract = {High Mountain Asia (HMA) is dependent upon both the amount and timing of snow and glacier meltwater. Previous model studies and coarse resolution (0.25° × 0.25°, ∼25 km × 25 km) passive microwave assessments of trends in the volume and timing of snowfall, snowmelt, and glacier melt in HMA have identified key spatial and seasonal heterogeneities in the response of snow to changes in regional climate. Here we use recently developed, continuous, internally consistent, and high-resolution passive microwave data (3.125 km × 3.125 km, 1987-2016) from the special sensor microwave imager instrument family to refine and extend previous estimates of changes in the snow regime of HMA. We find an overall decline in snow volume across HMA; however, there exist spatially contiguous regions of increasing snow volume—particularly during the winter season in the Pamir, Karakoram, Hindu Kush, and Kunlun Shan. Detailed analysis of changes in snow-volume trends through time reveal a large step change from negative trends during the period 1987-1997, to much more positive trends across large regions of HMA during the periods 1997-2007 and 2007-2016. We also find that changes in high percentile monthly snow-water volume exhibit steeper trends than changes in low percentile snow-water volume, which suggests a reduction in the frequency of high snow-water volumes in much of HMA. Regions with positive snow-water storage trends generally correspond to regions of positive glacier mass balances.}, language = {en} } @article{SmithBookhagen2020, author = {Smith, Taylor and Bookhagen, Bodo}, title = {Assessing Multi-Temporal Snow-Volume Trends in High Mountain Asia From 1987 to 2016 Using High-Resolution Passive Microwave Data}, series = {Frontiers in Earth Science}, volume = {8}, journal = {Frontiers in Earth Science}, publisher = {Frontiers Media}, address = {Lausanne}, issn = {2296-6463}, doi = {10.3389/feart.2020.559175}, pages = {13}, year = {2020}, abstract = {High Mountain Asia (HMA) is dependent upon both the amount and timing of snow and glacier meltwater. Previous model studies and coarse resolution (0.25° × 0.25°, ∼25 km × 25 km) passive microwave assessments of trends in the volume and timing of snowfall, snowmelt, and glacier melt in HMA have identified key spatial and seasonal heterogeneities in the response of snow to changes in regional climate. Here we use recently developed, continuous, internally consistent, and high-resolution passive microwave data (3.125 km × 3.125 km, 1987-2016) from the special sensor microwave imager instrument family to refine and extend previous estimates of changes in the snow regime of HMA. We find an overall decline in snow volume across HMA; however, there exist spatially contiguous regions of increasing snow volume—particularly during the winter season in the Pamir, Karakoram, Hindu Kush, and Kunlun Shan. Detailed analysis of changes in snow-volume trends through time reveal a large step change from negative trends during the period 1987-1997, to much more positive trends across large regions of HMA during the periods 1997-2007 and 2007-2016. We also find that changes in high percentile monthly snow-water volume exhibit steeper trends than changes in low percentile snow-water volume, which suggests a reduction in the frequency of high snow-water volumes in much of HMA. Regions with positive snow-water storage trends generally correspond to regions of positive glacier mass balances.}, language = {en} } @misc{MilewskiChabrillatBookhagen2020, author = {Milewski, Robert and Chabrillat, Sabine and Bookhagen, Bodo}, title = {Analyses of Namibian Seasonal Salt Pan Crust Dynamics and Climatic Drivers Using Landsat 8 Time-Series and Ground Data}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {988}, issn = {1866-8372}, doi = {10.25932/publishup-47568}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-475685}, pages = {26}, year = {2020}, abstract = {Salt pans are highly dynamic environments that are difficult to study by in situ methods because of their harsh climatic conditions and large spatial areas. Remote sensing can help to elucidate their environmental dynamics and provide important constraints regarding their sedimentological, mineralogical, and hydrological evolution. This study utilizes spaceborne multitemporal multispectral optical data combined with spectral endmembers to document spatial distribution of surface crust types over time on the Omongwa pan located in the Namibian Kalahari. For this purpose, 49 surface samples were collected for spectral and mineralogical characterization during three field campaigns (2014-2016) reflecting different seasons and surface conditions of the salt pan. An approach was developed to allow the spatiotemporal analysis of the salt pan crust dynamics in a dense time-series consisting of 77 Landsat 8 cloud-free scenes between 2014 and 2017, covering at least three major wet-dry cycles. The established spectral analysis technique Sequential Maximum Angle Convex Cone (SMACC) extraction method was used to derive image endmembers from the Landsat time-series stack. Evaluation of the extracted endmember set revealed that the multispectral data allowed the differentiation of four endmembers associated with mineralogical mixtures of the crust's composition in dry conditions and three endmembers associated with flooded or muddy pan conditions. The dry crust endmember spectra have been identified in relation to visible, near infrared, and short-wave infrared (VNIR-SWIR) spectroscopy and X-ray diffraction (XRD) analyses of the collected surface samples. According these results, the spectral endmembers are interpreted as efflorescent halite crust, mixed halite-gypsum crust, mixed calcite quartz sepiolite crust, and gypsum crust. For each Landsat scene the spatial distribution of these crust types was mapped with the Spectral Angle Mapper (SAM) method and significant spatiotemporal dynamics of the major surface crust types were observed. Further, the surface crust dynamics were analyzed in comparison with the pan's moisture regime and other climatic parameters. The results show that the crust dynamics are mainly driven by flooding events in the wet season, but are also influenced by temperature and aeolian activity in the dry season. The approach utilized in this study combines the advantages of multitemporal satellite data for temporal event characterization with advantages from hyperspectral methods for the image and ground data analyses that allow improved mineralogical differentiation and characterization.}, language = {en} } @article{VossBookhagenSachseetal.2020, author = {Voss, Katalyn A. and Bookhagen, Bodo and Sachse, Dirk and Chadwick, Oliver A.}, title = {Variation of deuterium excess in surface waters across a 5000-m elevation gradient in eastern Nepal}, series = {Journal of hydrology}, volume = {586}, journal = {Journal of hydrology}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0022-1694}, doi = {10.1016/j.jhydrol.2020.124802}, pages = {17}, year = {2020}, abstract = {The strong elevation gradient of the Himalaya allows for investigation of altitude and orographic impacts on surface water delta O-18 and delta D stable isotope values. This study differentiates the time- and altitude-variable contributions of source waters to the Arun River in eastern Nepal. It provides isotope data along a 5000-m gradient collected from tributaries as well as groundwater, snow, and glacial-sourced surface waters and time-series data from April to October 2016. We find nonlinear trends in delta O-18 and delta D lapse rates with high-elevation lapse rates (4000-6000 masl) 5-7 times more negative than low-elevation lapse rates (1000-3000 masl). A distinct seasonal signal in delta O-18 and delta D lapse rates indicates time-variable source-water contributions from glacial and snow meltwater as well as precipitation transitions between the Indian Summer Monsoon and Winter Westerly Disturbances. Deuterium excess correlates with the extent of snowpack and tracks melt events during the Indian Summer Monsoon season. Our analysis identifies the influence of snow and glacial melt waters on river composition during low-flow conditions before the monsoon (April/May 2016) followed by a 5-week transition to the Indian Summer Monsoon-sourced rainfall around mid-June 2016. In the post-monsoon season, we find continued influence from glacial melt waters as well as ISM-sourced groundwater.}, language = {en} } @article{CastinoBookhagenDelaTorre2020, author = {Castino, Fabiana and Bookhagen, Bodo and De la Torre, Alejandro}, title = {Atmospheric dynamics of extreme discharge events from 1979 to 2016 in the southern Central Andes}, series = {Climate dynamics : observational, theoretical and computational research on the climate system}, volume = {55}, journal = {Climate dynamics : observational, theoretical and computational research on the climate system}, number = {11-12}, publisher = {Springer}, address = {Berlin ; Heidelberg [u.a.]}, issn = {0930-7575}, doi = {10.1007/s00382-020-05458-1}, pages = {3485 -- 3505}, year = {2020}, abstract = {During the South-American Monsoon season, deep convective systems occur at the eastern flank of the Central Andes leading to heavy rainfall and flooding. We investigate the large- and meso-scale atmospheric dynamics associated with extreme discharge events (> 99.9th percentile) observed in two major river catchments meridionally stretching from humid to semi-arid conditions in the southern Central Andes. Based on daily gauge time series and ERA-Interim reanalysis, we made the following three key observations: (1) for the period 1940-2016 daily discharge exhibits more pronounced variability in the southern, semi-arid than in the northern, humid catchments. This is due to a smaller ratio of discharge magnitudes between intermediate (0.2 year return period) and rare events (20 year return period) in the semi-arid compared to the humid areas; (2) The climatological composites of the 40 largest discharge events showed characteristic atmospheric features of cold surges based on 5-day time-lagged sequences of geopotential height at different levels in the troposphere; (3) A subjective classification revealed that 80\% of the 40 largest discharge events are mainly associated with the north-northeastward migration of frontal systems and 2/3 of these are cold fronts, i.e. cold surges. This work highlights the importance of cold surges and their related atmospheric processes for the generation of heavy rainfall events and floods in the southern Central Andes.}, language = {en} } @article{MilewskiChabrillatBookhagen2020, author = {Milewski, Robert and Chabrillat, Sabine and Bookhagen, Bodo}, title = {Analyses of Namibian Seasonal Salt Pan Crust Dynamics and Climatic Drivers Using Landsat 8 Time-Series and Ground Data}, series = {Remote Sensing}, journal = {Remote Sensing}, number = {3}, publisher = {MDPI}, address = {Basel}, issn = {2072-4292}, doi = {10.3390/rs12030474}, pages = {24}, year = {2020}, abstract = {Salt pans are highly dynamic environments that are difficult to study by in situ methods because of their harsh climatic conditions and large spatial areas. Remote sensing can help to elucidate their environmental dynamics and provide important constraints regarding their sedimentological, mineralogical, and hydrological evolution. This study utilizes spaceborne multitemporal multispectral optical data combined with spectral endmembers to document spatial distribution of surface crust types over time on the Omongwa pan located in the Namibian Kalahari. For this purpose, 49 surface samples were collected for spectral and mineralogical characterization during three field campaigns (2014-2016) reflecting different seasons and surface conditions of the salt pan. An approach was developed to allow the spatiotemporal analysis of the salt pan crust dynamics in a dense time-series consisting of 77 Landsat 8 cloud-free scenes between 2014 and 2017, covering at least three major wet-dry cycles. The established spectral analysis technique Sequential Maximum Angle Convex Cone (SMACC) extraction method was used to derive image endmembers from the Landsat time-series stack. Evaluation of the extracted endmember set revealed that the multispectral data allowed the differentiation of four endmembers associated with mineralogical mixtures of the crust's composition in dry conditions and three endmembers associated with flooded or muddy pan conditions. The dry crust endmember spectra have been identified in relation to visible, near infrared, and short-wave infrared (VNIR-SWIR) spectroscopy and X-ray diffraction (XRD) analyses of the collected surface samples. According these results, the spectral endmembers are interpreted as efflorescent halite crust, mixed halite-gypsum crust, mixed calcite quartz sepiolite crust, and gypsum crust. For each Landsat scene the spatial distribution of these crust types was mapped with the Spectral Angle Mapper (SAM) method and significant spatiotemporal dynamics of the major surface crust types were observed. Further, the surface crust dynamics were analyzed in comparison with the pan's moisture regime and other climatic parameters. The results show that the crust dynamics are mainly driven by flooding events in the wet season, but are also influenced by temperature and aeolian activity in the dry season. The approach utilized in this study combines the advantages of multitemporal satellite data for temporal event characterization with advantages from hyperspectral methods for the image and ground data analyses that allow improved mineralogical differentiation and characterization.}, language = {en} } @article{OlenBookhagenStrecker2020, author = {Olen, Stephanie M. and Bookhagen, Bodo and Strecker, Manfred}, title = {Corrigendum to: Olen, Stephanie M.; Bookhagen, Bodo; Strecker, Manfred R. : Role of climate and vegetation density in modulating denudation rates in the Himalaya. - Earth and planetary science letters. - 445 (2016), S. 57 - 67. - doi: https://doi.org/10.1016/j.epsl.2016.03.047}, series = {Earth and planetary science letters}, volume = {540}, journal = {Earth and planetary science letters}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0012-821X}, doi = {10.1016/j.epsl.2020.116252}, pages = {5}, year = {2020}, abstract = {Vegetation has long been hypothesized to influence the nature and rates of surface processes. We test the possible impact of vegetation and climate on denudation rates at orogen scale by taking advantage of a pronounced along-strike gradient in rainfall and vegetation density in the Himalaya. We combine 12 new 10Be denudation rates from the Sutlej Valley and 123 published denudation rates from fluvially- dominated catchments in the Himalaya with remotely-sensed measures of vegetation density and rainfall metrics, and with tectonic and lithologic constraints. In addition, we perform topographic analyses to assess the contribution of vegetation and climate in modulating denudation rates along strike. We observe variations in denudation rates and the relationship between denudation and topography along strike that are most strongly controlled by local rainfall amount and vegetation density, and cannot be explained by along-strike differences in tectonics or lithology. A W-E along-strike decrease in denudation rate variability positively correlates with the seasonality of vegetation density (R = 0.95, p < 0.05), and negatively correlates with mean vegetation density (R = -0.84, p < 0.05). Vegetation density modulates the topographic response to changing denudation rates, such that the functional relationship between denudation rate and topographic steepness becomes increasingly linear as vegetation density increases. We suggest that while tectonic processes locally control the pattern of denudation rates across strike of the Himalaya (i.e., S-N), along strike of the orogen (i.e., E-W) climate exerts a measurable influence on how denudation rates scatter around long-term, tectonically-controlled erosion, and on the functional relationship between topography and denudation}, language = {en} } @article{LoiblBookhagenValadeetal.2019, author = {Loibl, David and Bookhagen, Bodo and Valade, Sebastien and Schneider, Christoph}, title = {OSARIS, the "Open Source SAR Investigation System" for Automatized Parallel InSAR Processing of Sentinel-1 Time Series Data With Special Emphasis on Cryosphere Applications}, series = {Frontiers in Earth Science}, volume = {7}, journal = {Frontiers in Earth Science}, publisher = {Frontiers Media}, address = {Lausanne}, issn = {2296-6463}, doi = {10.3389/feart.2019.00172}, pages = {20}, year = {2019}, abstract = {With the advent of the two Sentinel-1 (S1) satellites, Synthetic Aperture Radar (SAR) data with high temporal and spatial resolution are freely available. This provides a promising framework to facilitate detailed investigations of surface instabilities and movements on large scales with high temporal resolution, but also poses substantial processing challenges because of storage and computation requirements. Methods are needed to efficiently detect short term changes in dynamic environments. Approaches considering pair-wise processing of a series of consecutive scenes to retain maximum temporal resolution in conjunction with time series analyses are required. Here we present OSARIS, the "Open Source SAR Investigation System," as a framework to process large stacks of S1 data on high-performance computing clusters. Based on Generic Mapping Tools SAR, shell scripts, and the workload manager Slurm, OSARIS provides an open and modular framework combining parallelization of high-performance C programs, flexible processing schemes, convenient configuration, and generation of geocoded stacks of analysis-ready base data, including amplitude, phase, coherence, and unwrapped interferograms. Time series analyses can be conducted by applying automated modules to the data stacks. The capabilities of OSARIS are demonstrated in a case study from the northwestern Tien Shan, Central Asia. After merging of slices, a total of 80 scene pairs were processed from 174 total input scenes. The coherence time series exhibits pronounced seasonal variability, with relatively high coherence values prevailing during the summer months in the nival zone. As an example of a time series analysis module, we present OSARIS' "Unstable Coherence Metric" which identifies pixels affected by significant drops from high to low coherence values. Measurements of motion provided by LOSD measurements require careful evaluation because interferometric phase unwrapping is prone to errors. Here, OSARIS provides a series of modules to detect and mask unwrapping errors, correct for atmospheric disturbances, and remove large-scale trends. Wall clock processing time for the case study (area ~9,000 km2) was ~12 h 4 min on a machine with 400 cores and 2 TB RAM. In total, ~12 d 10 h 44 min (~96\%) were saved through parallelization. A comparison of selected OSARIS datasets to results from two state-of-the-art SAR processing suites, ISCE and SNAP, shows that OSARIS provides products of competitive quality despite its high level of automatization. OSARIS thus facilitates efficient S1-based region-wide investigations of surface movement events over multiple years.}, language = {en} } @article{WeldeabRuehlemannBookhagenetal.2019, author = {Weldeab, Syee and R{\"u}hlemann, Carsten and Bookhagen, Bodo and Pausata, Francesco S. R. and Perez-Lua, Fabiola M.}, title = {Enhanced Himalayan glacial melting during YD and H1 recorded in the Northern Bay of Bengal}, series = {Geochemistry, geophysics, geosystems}, volume = {20}, journal = {Geochemistry, geophysics, geosystems}, number = {5}, publisher = {American Geophysical Union}, address = {Washington}, issn = {1525-2027}, doi = {10.1029/2018GC008065}, pages = {2449 -- 2461}, year = {2019}, abstract = {Ocean-land thermal feedback mechanisms in the Indian Summer Monsoon (ISM) domain are an important but not well understood component of regional climate dynamics. Here we present a O-18 record analyzed in the mixed-layer dwelling planktonic foraminifer Globigerinoides ruber (sensu stricto) from the northernmost Bay of Bengal (BoB). The O-18 time series provides a spatially integrated measure of monsoonal precipitation and Himalayan meltwater runoff into the northern BoB and reveals two brief episodes of anomalously low O-18 values between 16.30.4 and 160.5 and 12.60.4 and 12.30.4 thousand years before present. The timing of these events is centered at Heinrich event 1 and the Younger Dryas, well-known phases of weak northern hemisphere monsoon systems. Numerical climate model experiments, simulating Heinrich event-like conditions, suggest a surface warming over the monsoon-dominated Himalaya and foreland in response to ISM weakening. Corroborating the simulation results, our analysis of published moraine exposure ages in the monsoon-dominated Himalaya indicates enhanced glacier retreats that, considering age model uncertainties, coincide and overlap with the episodes of anomalously low O-18 values in the northernmost BoB. Our climate proxy and simulation results provide insights into past regional climate dynamics, suggesting reduced cloud cover, increased solar radiation, and air warming of the Himalaya and foreland areas and, as a result, glacier mass losses in response to weakened ISM. Plain Language Summary Indian Summer Monsoon rainfall and Himalayan glacier/snow melts constitute the main water source for the densely populated Indian subcontinent. Better understanding of how future climate changes will affect the monsoon rainfall and Himalayan glaciers requires a long climate record. In this study, we create a 13,000-year-long climate record that allows us to better understand the response of Indian Summer Monsoon rainfall and Himalayan glaciers to past climate changes. The focus of our study is the time window between 9,000 and 22,000 years ago, an episode where the global climate experienced large and rapid changes. Our sediment record from the northern Bay of Bengal and climate change simulation indicate that during episodes of weak monsoon, the melting of the Himalayan glaciers increases substantially significantly. This is because the weakening of the monsoon results in less cloud cover and, as a result, the surface receives more sunlight and causes glacier melting.}, language = {en} } @article{MerchelGaertnerBeutneretal.2019, author = {Merchel, Silke and G{\"a}rtner, Andreas and Beutner, Sabrina and Bookhagen, Bodo and Chabilan, Amelie}, title = {Attempts to understand potential deficiencies in chemical procedures for AMS: Cleaning and dissolving quartz for Be-10 and Al-26 analysis}, series = {Nuclear instruments \& methods in physics research : a journal on accelerators, instrumentation and techniques applied to research in nuclear and atomic physics, materials science and related fields in physics, Section B, Beam interactions with materials and atoms}, volume = {455}, journal = {Nuclear instruments \& methods in physics research : a journal on accelerators, instrumentation and techniques applied to research in nuclear and atomic physics, materials science and related fields in physics, Section B, Beam interactions with materials and atoms}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0168-583X}, doi = {10.1016/j.nimb.2019.02.007}, pages = {293 -- 299}, year = {2019}, abstract = {The purity of the analysed samples (e.g. quartz) with respect to chemical composition and radionuclide contamination is essential for geomorphologic applications using so-called terrestrial cosmogenic nuclides (TCNs). To guarantee this, numerous cleaning and dissolution procedures have been developed. At the DREsden Accelerator Mass Spectrometry (DREAMS) facility, we also work on enhancing the chemical quartz-enrichment methodology from bulk rock and dissolution of quartz. Repeated exposure of the bulk material to acid mixtures (HCl/H2SiF6) at room temperature for cleaning and its monitoring by optical microscopy works for most quartz-rich samples. The quartz dissolution in HF under rather mild conditions (at room temperature on a shaker-table) has the advantage to leave difficult-to-dissolve minerals (e.g., tourmaline, zircon, rutile, sillimanite, kyanite, chromite, corundum), not separated by other physical methods before, as residue. Our comparison with a high-temperature dissolution method (in a microwave) indicates an additional amount of interfering elements, such as in average about 3 mg of Ti, more than 7 mg of Al, and about 22 mu g of Be (for 50 g SiO2), is added to the sample, hence showing the superiority of our mild method. This way, we reduce problems for chemistry and AMS, but also ensure better comparability to production rates of cleaner stoichiometric quartz from calibration sites.}, language = {en} } @article{ClubbBookhagenRheinwalt2019, author = {Clubb, Fiona J. and Bookhagen, Bodo and Rheinwalt, Aljoscha}, title = {Clustering river profiles to classify geomorphic domains}, series = {Journal of geophysical research : Earth surface}, volume = {124}, journal = {Journal of geophysical research : Earth surface}, number = {6}, publisher = {American Geophysical Union}, address = {Hoboken}, issn = {2169-9003}, doi = {10.1029/2019JF005025}, pages = {1417 -- 1439}, year = {2019}, abstract = {The structure and organization of river networks has been used for decades to investigate the influence of climate and tectonics on landscapes. The majority of these studies either analyze rivers in profile view by extracting channel steepness or calculate planform metrics such as drainage density. However, these techniques rely on the assumption of homogeneity: that intrinsic and external factors are spatially or temporally invariant over the measured profile. This assumption is violated for the majority of Earth's landscapes, where variations in uplift rate, rock strength, climate, and geomorphic process are almost ubiquitous. We propose a method for classifying river profiles to identify landscape regions with similar characteristics by adapting hierarchical clustering algorithms developed for time series data. We first test our clustering on two landscape evolution scenarios and find that we can successfully cluster regions with different erodibility and detect the transient response to sudden base level fall. We then test our method in two real landscapes: first in Bitterroot National Forest, Idaho, where we demonstrate that our method can detect transient incision waves and the topographic signature of fluvial and debris flow process regimes; and second, on Santa Cruz Island, California, where our technique identifies spatial patterns in lithology not detectable through normalized channel steepness analysis. By calculating channel steepness separately for each cluster, our method allows the extraction of more reliable steepness metrics than if calculated for the landscape as a whole. These examples demonstrate the method's ability to disentangle fluvial morphology in complex lithological and tectonic settings.}, language = {en} } @article{RheinwahltGoswamiBookhagen2019, author = {Rheinwahlt, Aljoscha and Goswami, Bedartha and Bookhagen, Bodo}, title = {A network-based flow accumulation algorithm for point clouds}, series = {Journal of geophysical research : Earth surface}, volume = {124}, journal = {Journal of geophysical research : Earth surface}, number = {7}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9003}, doi = {10.1029/2018JF004827}, pages = {2013 -- 2033}, year = {2019}, abstract = {Flow accumulation algorithms estimate the steady state of flow on real or modeled topographic surfaces and are crucial for hydrological and geomorphological assessments, including delineation of river networks, drainage basins, and sediment transport processes. Existing flow accumulation algorithms are typically designed to compute flows on regular grids and are not directly applicable to arbitrarily sampled topographic data such as lidar point clouds. In this study we present a random sampling scheme that generates homogeneous point densities, in combination with a novel flow path tracing approach-the Facet-Flow Network (FFN)-that estimates flow accumulation in terms of specific catchment area (SCA) on triangulated surfaces. The random sampling minimizes biases due to spatial sampling and the FFN allows for direct flow estimation from point clouds. We validate our approach on a Gaussian hill surface and study the convergence of its SCA compared to the analytical solution. Here, our algorithm outperforms the multiple flow direction algorithm, which is optimized for divergent surfaces. We also compute the SCA of a 6-km(2)-steep, vegetated catchment on Santa Cruz Island, California, based on airborne lidar point-cloud data. Point-cloud-based SCA values estimated by our method compare well with those estimated by the D-infinity or multiple flow direction algorithm on gridded data. The advantage of computing SCA from point clouds becomes relevant especially for divergent topography and for small drainage areas: These are depicted with much more detail due to the higher sampling density of point clouds.}, language = {en} } @misc{RamezaniZiaraniBookhagenSchmidtetal.2019, author = {Ramezani Ziarani, Maryam and Bookhagen, Bodo and Schmidt, Torsten and Wickert, Jens and de la Torre, Alejandro and Hierro, Rodrigo}, title = {Using Convective Available Potential Energy (CAPE) and Dew-Point Temperature to Characterize Rainfall-Extreme Events in the South-Central Andes}, series = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, number = {771}, issn = {1866-8372}, doi = {10.25932/publishup-43886}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-438865}, pages = {22}, year = {2019}, abstract = {The interactions between atmosphere and steep topography in the eastern south-central Andes result in complex relations with inhomogenous rainfall distributions. The atmospheric conditions leading to deep convection and extreme rainfall and their spatial patterns—both at the valley and mountain-belt scales—are not well understood. In this study, we aim to identify the dominant atmospheric conditions and their spatial variability by analyzing the convective available potential energy (CAPE) and dew-point temperature (Td). We explain the crucial effect of temperature on extreme rainfall generation along the steep climatic and topographic gradients in the NW Argentine Andes stretching from the low-elevation eastern foreland to the high-elevation central Andean Plateau in the west. Our analysis relies on version 2.0 of the ECMWF's (European Centre for Medium-RangeWeather Forecasts) Re-Analysis (ERA-interim) data and TRMM (Tropical Rainfall Measuring Mission) data. We make the following key observations: First, we observe distinctive gradients along and across strike of the Andes in dew-point temperature and CAPE that both control rainfall distributions. Second, we identify a nonlinear correlation between rainfall and a combination of dew-point temperature and CAPE through a multivariable regression analysis. The correlation changes in space along the climatic and topographic gradients and helps to explain controlling factors for extreme-rainfall generation. Third, we observe more contribution (or higher importance) of Td in the tropical low-elevation foreland and intermediate-elevation areas as compared to the high-elevation central Andean Plateau for 90th percentile rainfall. In contrast, we observe a higher contribution of CAPE in the intermediate-elevation area between low and high elevation, especially in the transition zone between the tropical and subtropical areas for the 90th percentile rainfall. Fourth, we find that the parameters of the multivariable regression using CAPE and Td can explain rainfall with higher statistical significance for the 90th percentile compared to lower rainfall percentiles. Based on our results, the spatial pattern of rainfall-extreme events during the past ∼16 years can be described by a combination of dew-point temperature and CAPE in the south-central Andes.}, language = {en} } @article{RamezaniZiaraniBookhagenSchmidtetal.2019, author = {Ramezani Ziarani, Maryam and Bookhagen, Bodo and Schmidt, Torsten and Wickert, Jens and de la Torre, Alejandro and Hierro, Rodrigo}, title = {Using Convective Available Potential Energy (CAPE) and Dew-Point Temperature to Characterize Rainfall-Extreme Events in the South-Central Andes}, series = {Atmosphere}, volume = {10}, journal = {Atmosphere}, number = {7}, publisher = {MDPI}, address = {Basel}, issn = {2073-4433}, doi = {10.3390/atmos10070379}, pages = {22}, year = {2019}, abstract = {The interactions between atmosphere and steep topography in the eastern south-central Andes result in complex relations with inhomogenous rainfall distributions. The atmospheric conditions leading to deep convection and extreme rainfall and their spatial patterns—both at the valley and mountain-belt scales—are not well understood. In this study, we aim to identify the dominant atmospheric conditions and their spatial variability by analyzing the convective available potential energy (CAPE) and dew-point temperature (Td). We explain the crucial effect of temperature on extreme rainfall generation along the steep climatic and topographic gradients in the NW Argentine Andes stretching from the low-elevation eastern foreland to the high-elevation central Andean Plateau in the west. Our analysis relies on version 2.0 of the ECMWF's (European Centre for Medium-RangeWeather Forecasts) Re-Analysis (ERA-interim) data and TRMM (Tropical Rainfall Measuring Mission) data. We make the following key observations: First, we observe distinctive gradients along and across strike of the Andes in dew-point temperature and CAPE that both control rainfall distributions. Second, we identify a nonlinear correlation between rainfall and a combination of dew-point temperature and CAPE through a multivariable regression analysis. The correlation changes in space along the climatic and topographic gradients and helps to explain controlling factors for extreme-rainfall generation. Third, we observe more contribution (or higher importance) of Td in the tropical low-elevation foreland and intermediate-elevation areas as compared to the high-elevation central Andean Plateau for 90th percentile rainfall. In contrast, we observe a higher contribution of CAPE in the intermediate-elevation area between low and high elevation, especially in the transition zone between the tropical and subtropical areas for the 90th percentile rainfall. Fourth, we find that the parameters of the multivariable regression using CAPE and Td can explain rainfall with higher statistical significance for the 90th percentile compared to lower rainfall percentiles. Based on our results, the spatial pattern of rainfall-extreme events during the past ∼16 years can be described by a combination of dew-point temperature and CAPE in the south-central Andes.}, language = {en} } @article{BoersGoswamiRheinwaltetal.2019, author = {Boers, Niklas and Goswami, Bedartha and Rheinwalt, Aljoscha and Bookhagen, Bodo and Hoskins, Brian and Kurths, J{\"u}rgen}, title = {Complex networks reveal global pattern of extreme-rainfall teleconnections}, series = {Nature : the international weekly journal of science}, volume = {566}, journal = {Nature : the international weekly journal of science}, number = {7744}, publisher = {Nature Publ. Group}, address = {London}, issn = {0028-0836}, doi = {10.1038/s41586-018-0872-x}, pages = {373 -- 377}, year = {2019}, abstract = {Climatic observables are often correlated across long spatial distances, and extreme events, such as heatwaves or floods, are typically assumed to be related to such teleconnections(1,2). Revealing atmospheric teleconnection patterns and understanding their underlying mechanisms is of great importance for weather forecasting in general and extreme-event prediction in particular(3,4), especially considering that the characteristics of extreme events have been suggested to change under ongoing anthropogenic climate change(5-8). Here we reveal the global coupling pattern of extreme-rainfall events by applying complex-network methodology to high-resolution satellite data and introducing a technique that corrects for multiple-comparison bias in functional networks. We find that the distance distribution of significant connections (P < 0.005) around the globe decays according to a power law up to distances of about 2,500 kilometres. For longer distances, the probability of significant connections is much higher than expected from the scaling of the power law. We attribute the shorter, power-law-distributed connections to regional weather systems. The longer, super-power-law-distributed connections form a global rainfall teleconnection pattern that is probably controlled by upper-level Rossby waves. We show that extreme-rainfall events in the monsoon systems of south-central Asia, east Asia and Africa are significantly synchronized. Moreover, we uncover concise links between south-central Asia and the European and North American extratropics, as well as the Southern Hemisphere extratropics. Analysis of the atmospheric conditions that lead to these teleconnections confirms Rossby waves as the physical mechanism underlying these global teleconnection patterns and emphasizes their crucial role in dynamical tropical-extratropical couplings. Our results provide insights into the function of Rossby waves in creating stable, global-scale dependencies of extreme-rainfall events, and into the potential predictability of associated natural hazards.}, language = {en} } @article{BrellSeglGuanteretal.2019, author = {Brell, Maximilian and Segl, Karl and Guanter, Luis and Bookhagen, Bodo}, title = {3D hyperspectral point cloud generation}, series = {ISPRS journal of photogrammetry and remote sensing : official publication of the International Society for Photogrammetry and Remote Sensing}, volume = {149}, journal = {ISPRS journal of photogrammetry and remote sensing : official publication of the International Society for Photogrammetry and Remote Sensing}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0924-2716}, doi = {10.1016/j.isprsjprs.2019.01.022}, pages = {200 -- 214}, year = {2019}, abstract = {Remote Sensing technologies allow to map biophysical, biochemical, and earth surface parameters of the land surface. Of especial interest for various applications in environmental and urban sciences is the combination of spectral and 3D elevation information. However, those two data streams are provided separately by different instruments, namely airborne laser scanner (ALS) for elevation and a hyperspectral imager (HSI) for high spectral resolution data. The fusion of ALS and HSI data can thus lead to a single data entity consistently featuring rich structural and spectral information. In this study, we present the application of fusing the first pulse return information from ALS data at a sub-decimeter spatial resolution with the lower-spatial resolution hyperspectral information available from the HSI into a hyperspectral point cloud (HSPC). During the processing, a plausible hyperspectral spectrum is assigned to every first-return ALS point. We show that the complementary implementation of spectral and 3D information at the point-cloud scale improves object-based classification and information extraction schemes. This improvements have great potential for numerous land cover mapping and environmental applications.}, language = {en} } @misc{SmithRheinwaltBookhagen2019, author = {Smith, Taylor and Rheinwalt, Aljoscha and Bookhagen, Bodo}, title = {Determining the optimal grid resolution for topographic analysis on an airborne lidar dataset}, series = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, number = {725}, issn = {1866-8372}, doi = {10.25932/publishup-43016}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-430165}, pages = {475 -- 489}, year = {2019}, abstract = {Digital elevation models (DEMs) are a gridded representation of the surface of the Earth and typically contain uncertainties due to data collection and processing. Slope and aspect estimates on a DEM contain errors and uncertainties inherited from the representation of a continuous surface as a grid (referred to as truncation error; TE) and from any DEM uncertainty. We analyze in detail the impacts of TE and propagated elevation uncertainty (PEU) on slope and aspect. Using synthetic data as a control, we define functions to quantify both TE and PEU for arbitrary grids. We then develop a quality metric which captures the combined impact of both TE and PEU on the calculation of topographic metrics. Our quality metric allows us to examine the spatial patterns of error and uncertainty in topographic metrics and to compare calculations on DEMs of different sizes and accuracies. Using lidar data with point density of ∼10 pts m-2 covering Santa Cruz Island in southern California, we are able to generate DEMs and uncertainty estimates at several grid resolutions. Slope (aspect) errors on the 1 m dataset are on average 0.3∘ (0.9∘) from TE and 5.5∘ (14.5∘) from PEU. We calculate an optimal DEM resolution for our SCI lidar dataset of 4 m that minimizes the error bounds on topographic metric calculations due to the combined influence of TE and PEU for both slope and aspect calculations over the entire SCI. Average slope (aspect) errors from the 4 m DEM are 0.25∘ (0.75∘) from TE and 5∘ (12.5∘) from PEU. While the smallest grid resolution possible from the high-density SCI lidar is not necessarily optimal for calculating topographic metrics, high point-density data are essential for measuring DEM uncertainty across a range of resolutions.}, language = {en} } @article{PurintonBookhagen2019, author = {Purinton, Benjamin and Bookhagen, Bodo}, title = {Introducing PebbleCounts}, series = {Earth Surface Dynamics}, volume = {2019}, journal = {Earth Surface Dynamics}, number = {7}, publisher = {Copernicus Publ}, address = {G{\"o}ttingen}, issn = {2196-6311}, doi = {10.5194/esurf-7-859-2019}, pages = {859 -- 877}, year = {2019}, abstract = {Grain-size distributions are a key geomorphic metric of gravel-bed rivers. Traditional measurement methods include manual counting or photo sieving, but these are achievable only at the 1-10 ㎡ scale. With the advent of drones and increasingly high-resolution cameras, we can now generate orthoimagery over hectares at millimeter to centimeter resolution. These scales, along with the complexity of high-mountain rivers, necessitate different approaches for photo sieving. As opposed to other image segmentation methods that use a watershed approach, our open-source algorithm, PebbleCounts, relies on k-means clustering in the spatial and spectral domain and rapid manual selection of well-delineated grains. This improves grain-size estimates for complex riverbed imagery, without post-processing. We also develop a fully automated method, PebbleCountsAuto, that relies on edge detection and filtering suspect grains, without the k-means clustering or manual selection steps. The algorithms are tested in controlled indoor conditions on three arrays of pebbles and then applied to 12 × 1 ㎡ orthomosaic clips of high-energy mountain rivers collected with a camera-on-mast setup (akin to a low-flying drone). A 20-pixel b-axis length lower truncation is necessary for attaining accurate grain-size distributions. For the k-means PebbleCounts approach, average percentile bias and precision are 0.03 and 0.09 ψ, respectively, for ∼1.16 mm pixel⁻¹ images, and 0.07 and 0.05 ψ for one 0.32 mm pixel⁻¹ image. The automatic approach has higher bias and precision of 0.13 and 0.15 ψ, respectively, for ∼1.16 mm pixel⁻¹ images, but similar values of -0.06 and 0.05 ψ for one 0.32 mm pixel⁻¹ image. For the automatic approach, only at best 70 \% of the grains are correct identifications, and typically around 50 \%. PebbleCounts operates most effectively at the 1 ㎡ patch scale, where it can be applied in ∼5-10 min on many patches to acquire accurate grain-size data over 10-100 ㎡ areas. These data can be used to validate PebbleCountsAuto, which may be applied at the scale of entire survey sites (102-104 ㎡ ). We synthesize results and recommend best practices for image collection, orthomosaic generation, and grain-size measurement using both algorithms.}, language = {en} } @article{SmithRheinwaltBookhagen2019, author = {Smith, Taylor and Rheinwalt, Aljoscha and Bookhagen, Bodo}, title = {Determining the optimal grid resolution for topographic analysis on an airborne lidar dataset}, series = {Earth Surface Dynamics}, volume = {7}, journal = {Earth Surface Dynamics}, publisher = {Copernicus Publ.}, address = {G{\"o}ttingen}, issn = {2196-6311}, doi = {10.5194/esurf-7-475-2019}, pages = {475 -- 489}, year = {2019}, abstract = {Digital elevation models (DEMs) are a gridded representation of the surface of the Earth and typically contain uncertainties due to data collection and processing. Slope and aspect estimates on a DEM contain errors and uncertainties inherited from the representation of a continuous surface as a grid (referred to as truncation error; TE) and from any DEM uncertainty. We analyze in detail the impacts of TE and propagated elevation uncertainty (PEU) on slope and aspect. Using synthetic data as a control, we define functions to quantify both TE and PEU for arbitrary grids. We then develop a quality metric which captures the combined impact of both TE and PEU on the calculation of topographic metrics. Our quality metric allows us to examine the spatial patterns of error and uncertainty in topographic metrics and to compare calculations on DEMs of different sizes and accuracies. Using lidar data with point density of ∼10 pts m-2 covering Santa Cruz Island in southern California, we are able to generate DEMs and uncertainty estimates at several grid resolutions. Slope (aspect) errors on the 1 m dataset are on average 0.3∘ (0.9∘) from TE and 5.5∘ (14.5∘) from PEU. We calculate an optimal DEM resolution for our SCI lidar dataset of 4 m that minimizes the error bounds on topographic metric calculations due to the combined influence of TE and PEU for both slope and aspect calculations over the entire SCI. Average slope (aspect) errors from the 4 m DEM are 0.25∘ (0.75∘) from TE and 5∘ (12.5∘) from PEU. While the smallest grid resolution possible from the high-density SCI lidar is not necessarily optimal for calculating topographic metrics, high point-density data are essential for measuring DEM uncertainty across a range of resolutions.}, language = {en} } @misc{PurintonBookhagen2019, author = {Purinton, Benjamin and Bookhagen, Bodo}, title = {Introducing PebbleCounts}, series = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, number = {783}, issn = {1866-8372}, doi = {10.25932/publishup-43946}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-439468}, pages = {21}, year = {2019}, abstract = {Grain-size distributions are a key geomorphic metric of gravel-bed rivers. Traditional measurement methods include manual counting or photo sieving, but these are achievable only at the 1-10 ㎡ scale. With the advent of drones and increasingly high-resolution cameras, we can now generate orthoimagery over hectares at millimeter to centimeter resolution. These scales, along with the complexity of high-mountain rivers, necessitate different approaches for photo sieving. As opposed to other image segmentation methods that use a watershed approach, our open-source algorithm, PebbleCounts, relies on k-means clustering in the spatial and spectral domain and rapid manual selection of well-delineated grains. This improves grain-size estimates for complex riverbed imagery, without post-processing. We also develop a fully automated method, PebbleCountsAuto, that relies on edge detection and filtering suspect grains, without the k-means clustering or manual selection steps. The algorithms are tested in controlled indoor conditions on three arrays of pebbles and then applied to 12 × 1 ㎡ orthomosaic clips of high-energy mountain rivers collected with a camera-on-mast setup (akin to a low-flying drone). A 20-pixel b-axis length lower truncation is necessary for attaining accurate grain-size distributions. For the k-means PebbleCounts approach, average percentile bias and precision are 0.03 and 0.09 ψ, respectively, for ∼1.16 mm pixel⁻¹ images, and 0.07 and 0.05 ψ for one 0.32 mm pixel⁻¹ image. The automatic approach has higher bias and precision of 0.13 and 0.15 ψ, respectively, for ∼1.16 mm pixel⁻¹ images, but similar values of -0.06 and 0.05 ψ for one 0.32 mm pixel⁻¹ image. For the automatic approach, only at best 70 \% of the grains are correct identifications, and typically around 50 \%. PebbleCounts operates most effectively at the 1 ㎡ patch scale, where it can be applied in ∼5-10 min on many patches to acquire accurate grain-size data over 10-100 ㎡ areas. These data can be used to validate PebbleCountsAuto, which may be applied at the scale of entire survey sites (102-104 ㎡ ). We synthesize results and recommend best practices for image collection, orthomosaic generation, and grain-size measurement using both algorithms.}, language = {en} } @article{TeshebaevaEchtlerBookhagenetal.2019, author = {Teshebaeva, Kanayim and Echtler, Helmut and Bookhagen, Bodo and Strecker, Manfred}, title = {Deep-seated gravitational slope deformation (DSGSD) and slow-moving landslides in the southern Tien Shan Mountains: new insights from InSAR, tectonic and geomorphic analysis}, series = {Earth surface processes and landforms : the journal of the British Geomorphological Research Group}, volume = {44}, journal = {Earth surface processes and landforms : the journal of the British Geomorphological Research Group}, number = {12}, publisher = {Wiley}, address = {Hoboken}, issn = {0197-9337}, doi = {10.1002/esp.4648}, pages = {2333 -- 2348}, year = {2019}, abstract = {We investigated deep-seated gravitational slope deformation (DSGSD) and slow mass movements in the southern Tien Shan Mountains front using synthetic aperture radar (SAR) time-series data obtained by the ALOS/PALSAR satellite. DSGSD evolves with a variety of geomorphological changes (e.g. valley erosion, incision of slope drainage networks) over time that affect earth surfaces and, therefore, often remain unexplored. We analysed 118 interferograms generated from 20 SAR images that covered about 900 km(2). To understand the spatial pattern of the slope movements and to identify triggering parameters, we correlated surface dynamics with the tectono-geomorphic processes and lithologic conditions of the active front of the Alai Range. We observed spatially continuous, constant hillslope movements with a downslope speed of approximately 71 mm year(-1) velocity. Our findings suggest that the lithological and structural framework defined by protracted deformation was the main controlling factor for sustained relief and, consequently, downslope mass movements. The analysed structures revealed integration of a geological/structural setting with the superposition of Cretaceous-Paleogene alternating carbonatic and clastic sedimentary structures as the substratum for younger, less consolidated sediments. This type of structural setting causes the development of large-scale, gravity-driven DSGSD and slow mass movement. Surface deformations with clear scarps and multiple crest lines triggered planes for large-scale deep mass creeps, and these were related directly to active faults and folds in the geologic structures. Our study offers a new combination of InSAR techniques and structural field observations, along with morphometric and seismologic correlations, to identify and quantify slope instability phenomena along a tectonically active mountain front. These results contribute to an improved natural risk assessment in these structures.}, language = {en} } @misc{BriegerHerzschuhPestryakovaetal.2019, author = {Brieger, Frederic and Herzschuh, Ulrike and Pestryakova, Luidmila Agafyevna and Bookhagen, Bodo and Zakharov, Evgenii S. and Kruse, Stefan}, title = {Advances in the derivation of Northeast Siberian forest metrics using high-resolution UAV-based photogrammetric point clouds}, series = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {1337}, issn = {1866-8372}, doi = {10.25932/publishup-47331}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-473318}, pages = {24}, year = {2019}, abstract = {Forest structure is a crucial component in the assessment of whether a forest is likely to act as a carbon sink under changing climate. Detailed 3D structural information about the tundra-taiga ecotone of Siberia is mostly missing and still underrepresented in current research due to the remoteness and restricted accessibility. Field based, high-resolution remote sensing can provide important knowledge for the understanding of vegetation properties and dynamics. In this study, we test the applicability of consumer-grade Unmanned Aerial Vehicles (UAVs) for rapid calculation of stand metrics in treeline forests. We reconstructed high-resolution photogrammetric point clouds and derived canopy height models for 10 study sites from NE Chukotka and SW Yakutia. Subsequently, we detected individual tree tops using a variable-window size local maximum filter and applied a marker-controlled watershed segmentation for the delineation of tree crowns. With this, we successfully detected 67.1\% of the validation individuals. Simple linear regressions of observed and detected metrics show a better correlation (R2) and lower relative root mean square percentage error (RMSE\%) for tree heights (mean R2 = 0.77, mean RMSE\% = 18.46\%) than for crown diameters (mean R2 = 0.46, mean RMSE\% = 24.9\%). The comparison between detected and observed tree height distributions revealed that our tree detection method was unable to representatively identify trees <2 m. Our results show that plot sizes for vegetation surveys in the tundra-taiga ecotone should be adapted to the forest structure and have a radius of >15-20 m to capture homogeneous and representative forest stands. Additionally, we identify sources of omission and commission errors and give recommendations for their mitigation. In summary, the efficiency of the used method depends on the complexity of the forest's stand structure.}, language = {en} } @article{BriegerHerzschuhPestryakovaetal.2019, author = {Brieger, Frederic and Herzschuh, Ulrike and Pestryakova, Luidmila Agafyevna and Bookhagen, Bodo and Zakharov, Evgenii S. and Kruse, Stefan}, title = {Advances in the Derivation of Northeast Siberian Forest Metrics Using High-Resolution UAV-Based Photogrammetric Point Clouds}, series = {Remote sensing}, volume = {11}, journal = {Remote sensing}, number = {12}, publisher = {MDPI}, address = {Basel}, issn = {2072-4292}, doi = {10.3390/rs11121447}, pages = {24}, year = {2019}, abstract = {Forest structure is a crucial component in the assessment of whether a forest is likely to act as a carbon sink under changing climate. Detailed 3D structural information about the tundra-taiga ecotone of Siberia is mostly missing and still underrepresented in current research due to the remoteness and restricted accessibility. Field based, high-resolution remote sensing can provide important knowledge for the understanding of vegetation properties and dynamics. In this study, we test the applicability of consumer-grade Unmanned Aerial Vehicles (UAVs) for rapid calculation of stand metrics in treeline forests. We reconstructed high-resolution photogrammetric point clouds and derived canopy height models for 10 study sites from NE Chukotka and SW Yakutia. Subsequently, we detected individual tree tops using a variable-window size local maximum filter and applied a marker-controlled watershed segmentation for the delineation of tree crowns. With this, we successfully detected 67.1\% of the validation individuals. Simple linear regressions of observed and detected metrics show a better correlation (R2) and lower relative root mean square percentage error (RMSE\%) for tree heights (mean R2 = 0.77, mean RMSE\% = 18.46\%) than for crown diameters (mean R2 = 0.46, mean RMSE\% = 24.9\%). The comparison between detected and observed tree height distributions revealed that our tree detection method was unable to representatively identify trees <2 m. Our results show that plot sizes for vegetation surveys in the tundra-taiga ecotone should be adapted to the forest structure and have a radius of >15-20 m to capture homogeneous and representative forest stands. Additionally, we identify sources of omission and commission errors and give recommendations for their mitigation. In summary, the efficiency of the used method depends on the complexity of the forest's stand structure.}, language = {en} } @misc{OlenBookhagen2018, author = {Olen, Stephanie M. and Bookhagen, Bodo}, title = {Mapping Damage-Affected Areas after Natural Hazard Events Using Sentinel-1 Coherence Time Series}, series = {remote sensing}, journal = {remote sensing}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-417766}, pages = {19}, year = {2018}, abstract = {The emergence of the Sentinel-1A and 1B satellites now offers freely available and widely accessible Synthetic Aperture Radar (SAR) data. Near-global coverage and rapid repeat time (6-12 days) gives Sentinel-1 data the potential to be widely used for monitoring the Earth's surface. Subtle land-cover and land surface changes can affect the phase and amplitude of the C-band SAR signal, and thus the coherence between two images collected before and after such changes. Analysis of SAR coherence therefore serves as a rapidly deployable and powerful tool to track both seasonal changes and rapid surface disturbances following natural disasters. An advantage of using Sentinel-1 C-band radar data is the ability to easily construct time series of coherence for a region of interest at low cost. In this paper, we propose a new method for Potentially Affected Area (PAA) detection following a natural hazard event. Based on the coherence time series, the proposed method (1) determines the natural variability of coherence within each pixel in the region of interest, accounting for factors such as seasonality and the inherent noise of variable surfaces; and (2) compares pixel-by-pixel syn-event coherence to temporal coherence distributions to determine where statistically significant coherence loss has occurred. The user can determine to what degree the syn-event coherence value (e.g., 1st, 5th percentile of pre-event distribution) constitutes a PAA, and integrate pertinent regional data, such as population density, to rank and prioritise PAAs. We apply the method to two case studies, Sarpol-e, Iran following the 2017 Iran-Iraq earthquake, and a landslide-prone region of NW Argentina, to demonstrate how rapid identification and interpretation of potentially affected areas can be performed shortly following a natural hazard event.}, language = {en} } @misc{SmithBookhagen2018, author = {Smith, Taylor and Bookhagen, Bodo}, title = {Using passive microwave data to understand spatio-temporal trends and dynamics in snow-water storage in High Mountain Asia}, series = {active and passive microwave remote sensing for environmental monitoring II}, volume = {10788}, journal = {active and passive microwave remote sensing for environmental monitoring II}, publisher = {SPIE-INT Soc Optical Engineering}, address = {Bellingham}, isbn = {978-1-5106-2160-2}, issn = {0277-786X}, doi = {10.1117/12.2323827}, pages = {8}, year = {2018}, abstract = {High Mountain Asia provides water for more than a billion downstream users. Many catchments receive the majority of their yearly water budget in the form of snow - the vast majority of which is not monitored by sparse weather networks. We leverage passive microwave data from the SSMI series of satellites (SSMI, SSMI/S, 1987-2016), reprocessed to 3.125 km resolution, to examine trends in the volume and spatial distribution of snow-water equivalent (SWE) in the Indus Basin. We find that the majority of the Indus has seen an increase in snow-water storage. There exists a strong elevation-trend relationship, where high-elevation zones have more positive SWE trends. Negative trends are confined to the Himalayan foreland and deeply-incised valleys which run into the Upper Indus. This implies a temperature-dependent cutoff below which precipitation increases are not translated into increased SWE. Earlier snowmelt or a higher percentage of liquid precipitation could both explain this cutoff.(1) Earlier work 2 found a negative snow-water storage trend for the entire Indus catchment over the time period 1987-2009 (-4 x 10(-3) mm/yr). In this study based on an additional seven years of data, the average trend reverses to 1.4 x 10(-3). This implies that the decade since the mid-2000s was likely wetter, and positively impacted long-term SWE trends. This conclusion is supported by an analysis of snowmelt onset and end dates which found that while long-term trends are negative, more recent (since 2005) trends are positive (moving later in the year).(3)}, language = {en} } @article{OezcanBookhagenMusaoglu2018, author = {{\"O}zcan, Orkan and Bookhagen, Bodo and Musaoglu, Nebiye}, title = {Impact of the Ataturk Dam Lake on Agro-Meteorological Aspects of the Southeastern Anatolia Region, Turkey}, series = {Journal of the Indian Society of Remote Sensing}, volume = {46}, journal = {Journal of the Indian Society of Remote Sensing}, number = {3}, publisher = {Springer}, address = {New York}, issn = {0255-660X}, doi = {10.1007/s12524-017-0703-9}, pages = {471 -- 481}, year = {2018}, abstract = {In this study, the spatial and temporal impacts of the Ataturk Dam on agro-meteorological aspects of the Southeastern Anatolia region have been investigated. Change detection and environmental impacts due to water-reserve changes in Ataturk Dam Lake have been determined and evaluated using multi-temporal Landsat satellite imageries and meteorological datasets within a period of 1984-2011. These time series have been evaluated for three time periods. Dam construction period constitutes the first part of the study. Land cover/use changes especially on agricultural fields under the Ataturk Dam Lake and its vicinity have been identified between the periods of 1984-1992. The second period comprises the 10-year period after the completion of filling up the reservoir in 1992. At this period, Landsat and meteorological time-series analyses are examined to assess the impact of the Ataturk Dam Lake on selected irrigated agricultural areas. For the last 9-year period from 2002 to 2011, the relationships between seasonal water-reserve changes and irrigated plains under changing climatic factors primarily driving vegetation activity (monthly, seasonal, and annual fluctuations of rainfall rate, air temperature, humidity) on the watershed have been investigated using a 30-year meteorological time series. The results showed that approximately 368 km(2) of agricultural fields have been affected because of inundation due to the Ataturk Dam Lake. However, irrigated agricultural fields have been increased by 56.3\% of the total area (1552 of 2756 km(2)) on Harran Plain within the period of 1984-2011.}, language = {en} } @misc{RheinwaltBookhagen2018, author = {Rheinwalt, Aljoscha and Bookhagen, Bodo}, title = {Network-based flow accumulation for point clouds}, series = {Remote Sensing for Agriculture, Ecosystems, and Hydrology XX}, volume = {10783}, journal = {Remote Sensing for Agriculture, Ecosystems, and Hydrology XX}, publisher = {SPIE-INT Society of Photo-Optical Instrumentation Engineers}, address = {Bellingham}, isbn = {978-1-5106-2150-3}, issn = {0277-786X}, doi = {10.1117/12.2318424}, pages = {12}, year = {2018}, abstract = {Point clouds provide high-resolution topographic data which is often classified into bare-earth, vegetation, and building points and then filtered and aggregated to gridded Digital Elevation Models (DEMs) or Digital Terrain Models (DTMs). Based on these equally-spaced grids flow-accumulation algorithms are applied to describe the hydrologic and geomorphologic mass transport on the surface. In this contribution, we propose a stochastic point-cloud filtering that, together with a spatial bootstrap sampling, allows for a flow accumulation directly on point clouds using Facet-Flow Networks (FFN). Additionally, this provides a framework for the quantification of uncertainties in point-cloud derived metrics such as Specific Catchment Area (SCA) even though the flow accumulation itself is deterministic.}, language = {en} } @article{SmithBookhagen2018, author = {Smith, Taylor and Bookhagen, Bodo}, title = {Changes in seasonal snow water equivalent distribution in High Mountain Asia (1987 to 2009)}, series = {Science Advances}, volume = {4}, journal = {Science Advances}, number = {1}, publisher = {American Assoc. for the Advancement of Science}, address = {Washington}, issn = {2375-2548}, doi = {10.1126/sciadv.1701550}, pages = {8}, year = {2018}, abstract = {Snow meltwaters account for most of the yearly water budgets of many catchments in High Mountain Asia (HMA). We examine trends in snow water equivalent (SWE) using passive microwave data (1987 to 2009). We find an overall decrease in SWE in HMA, despite regions of increased SWE in the Pamir, Kunlun Shan, Eastern Himalaya, and Eastern Tien Shan. Although the average decline in annual SWE across HMA (contributing area, 2641 x 10(3) km(2)) is low (average, -0.3\%), annual SWE losses conceal distinct seasonal and spatial heterogeneities across the study region. For example, the Tien Shan has seen both strong increases in winter SWE and sharp declines in spring and summer SWE. In the majority of catchments, the most negative SWE trends are found in mid-elevation zones, which often correspond to the regions of highest snow-water storage and are somewhat distinct from glaciated areas. Negative changes in SWE storage in these mid-elevation zones have strong implications for downstream water availability.}, language = {en} } @article{MeeseBookhagenOlenetal.2018, author = {Meese, Bernd and Bookhagen, Bodo and Olen, Stephanie M. and Barthold, Frauke Katrin and Sachse, Dirk}, title = {The effect of Indian Summer Monsoon rainfall on surface water delta D values in the central Himalaya}, series = {Hydrological processes}, volume = {32}, journal = {Hydrological processes}, number = {24}, publisher = {Wiley}, address = {Hoboken}, issn = {0885-6087}, doi = {10.1002/hyp.13281}, pages = {3662 -- 3674}, year = {2018}, abstract = {Stable isotope proxy records, such as speleothems, plant-wax biomarker records, and ice cores, are suitable archives for the reconstruction of regional palaeohydrologic conditions. But the interpretation of these records in the tropics, especially in the Indian Summer Monsoon (ISM) domain, is difficult due to differing moisture and water sources: precipitation from the ISM and Winter Westerlies, as well as snow- and glacial meltwater. In this study, we use interannual differences in ISM strength (2011-2012) to understand the stable isotopic composition of surface water in the Arun River catchment in eastern Nepal. We sampled main stem and tributary water (n = 204) for stable hydrogen and oxygen isotope analysis in the postmonsoon phase of two subsequent years with significantly distinct ISM intensities. In addition to the 2011/2012 sampling campaigns, we collected a 12-month time series of main stem waters (2012/2013, n = 105) in order to better quantify seasonal effects on the variability of surface water delta O-18/delta D. Furthermore, remotely sensed satellite data of rainfall, snow cover, glacial coverage, and evapotranspiration was evaluated. The comparison of datasets from both years revealed that surface waters of the main stem Arun and its tributaries were D-enriched by similar to 15 parts per thousand when ISM rainfall decreased by 20\%. This strong response emphasizes the importance of the ISM for surface water run-off in the central Himalaya. However, further spatio-temporal analysis of remote sensing data in combination with stream water d-excess revealed that most high-altitude tributaries and the Tibetan part of the Arun receive high portions of glacial melt water and likely Winter Westerly Disturbances precipitation. We make the following two implications: First, palaeohydrologic archives found in high-altitude tributaries and on the southern Tibetan Plateau record a mixture of past precipitation delta D values and variable amounts of additional water sources. Second, surface water isotope ratios of lower elevated tributaries strongly reflect the isotopic composition of ISM rainfall implying a suitable region for the analysis of potential delta D value proxy records.}, language = {en} } @article{LunaBookhagenNiedermannetal.2018, author = {Luna, Lisa Victoria and Bookhagen, Bodo and Niedermann, Samuel and Rugel, Georg and Scharf, Andreas and Merchel, Silke}, title = {Glacial chronology and production rate cross-calibration of five cosmogenic nuclide and mineral systems from the southern Central Andean Plateau}, series = {Earth \& planetary science letters}, volume = {500}, journal = {Earth \& planetary science letters}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0012-821X}, doi = {10.1016/j.epsl.2018.07.034}, pages = {242 -- 253}, year = {2018}, abstract = {Glacial deposits on the high-altitude, arid southern Central Andean Plateau (CAP), the Puna in northwestern Argentina, document past changes in climate, but the associated geomorphic features have rarely been directly dated. This study provides direct age control of glacial moraine deposits from the central Puna (24 degrees S) at elevations of 3900-5000 m through surface exposure dating with cosmogenic nuclides. Our results show that the most extensive glaciations occurred before 95 ka and an additional major advance occurred between 46 and 39 ka. The latter period is synchronous with the highest lake levels in the nearby Pozuelos basin and the Minchin (Inca Huasi) wet phase on the Altiplano in the northern CAP. None of the dated moraines produced boulder ages corresponding to the Tauca wet phase (24-15 ka). Additionally, the volcanic lithologies of the deposits allow us to establish production ratios at low latitude and high elevation for five different nuclide and mineral systems: Be-10, Ne-21, and Al-26 from quartz (11 or 12 samples) and He-3 and Ne-21 from pyroxene (10 samples). We present production ratios for all combinations of the measured nuclides and cross-calibrated production rates for 21Ne in pyroxene and quartz for the high, (sub-)tropical Andes. The production rates are based on our Be-10-normalized production ratios and a weighted mean of reference 10Be production rates calibrated in the high, tropical Andes (4.02 +/- 0.12 at g(-1) yr(-1)). These are, Ne-21(qtz): 18.1 +/- 1.2 at g(-1) yr(-1) and Ne-21(px): 36.6 +/- 1.8 at g(-1) yr(-1) (En(88-94)) scaled to sea level and high latitude using the Lal/Stone scheme, with 1 sigma uncertainties. As He-3 and Al-26 have been directly calibrated in the tropical Andes, we recommend using those rates. Finally, we compare exposure ages calculated using all measured cosmogenic nuclides from each sample, including 11 feldspar samples measured for Cl-36, and a suite of previously published production rates. (C) 2018 Published by Elsevier B.V.}, language = {en} } @article{ThompsonChenYangetal.2018, author = {Thompson, Jessica A. and Chen, Jie and Yang, Huili and Li, Tao and Bookhagen, Bodo and Burbank, Douglas}, title = {Coarse- versus fine-grain quartz OSL and cosmogenic Be-10 dating of deformed fluvial terraces on the northeast Pamir margin, northwest China}, series = {Quaternary geochronology : the international research and review journal on advances in quaternary dating techniques}, volume = {46}, journal = {Quaternary geochronology : the international research and review journal on advances in quaternary dating techniques}, publisher = {Elsevier}, address = {Oxford}, issn = {1871-1014}, doi = {10.1016/j.quageo.2018.01.002}, pages = {1 -- 15}, year = {2018}, abstract = {Along the NE Pamir margin, flights of late Quaternary fluvial terraces span actively deforming fault-related folds. We present detailed results on two terraces dated using optically stimulated luminescence (OSL) and cosmogenic radionuclide Be-10 (CRN) techniques. Quartz OSL dating of two different grain sizes (4-11 mu m and 90-180 mu m) revealed the fine-grain quartz fraction may overestimate the terrace ages by up to a factor of ten. Two-mm, small-aliquot, coarse-grain quartz OSL ages, calculated using the minimum age model, yielded stratigraphically consistent ages within error and dated times of terrace deposition to similar to 9 and similar to 16 ka. We speculate that, in this arid environment, fine-grain samples can be transported and deposited in single, turbid, and (sometimes) night-time floods that prevent thorough bleaching and, thereby, can lead to relatively large residual OSL signals. In contrast, sand in the fluvial system is likely to have a much longer residence time during transport, thereby providing greater opportunities for thorough bleaching. CRN Be-10 depth profiles date the timing of terrace abandonment to similar to 8 and similar to 14 ka: ages that generally agree with the coarse-grain quartz OSL ages. Our new terrace age of similar to 13-14 ka is broadly consistent with other terraces in the region that indicate terrace deposition and subsequent abandonment occurred primarily during glacial-interglacial transitions, thereby suggesting a climatic control on the formation of these terraces on the margins of the Tarim Basin. Furthermore, tectonic shortening rates calculated from these deformed terraces range from similar to 1.2 to similar to 4.6 mm/a and, when combined with shortening rates from other structures in the region, illuminate the late Quaternary basinward migration of deformation to faults and folds along the Pamir-Tian Shan collisional interface.}, language = {en} } @article{NennewitzThiedeBookhagen2018, author = {Nennewitz, Markus and Thiede, Rasmus C. and Bookhagen, Bodo}, title = {Fault activity, tectonic segmentation, and deformation pattern of the western Himalaya on Ma timescales inferred from landscape morphology}, series = {Lithosphere}, volume = {10}, journal = {Lithosphere}, number = {5}, publisher = {American Institute of Physics}, address = {Boulder}, issn = {1941-8264}, doi = {10.1130/L681.1}, pages = {632 -- 640}, year = {2018}, abstract = {The location and magnitude of Himalayan tectonic activity has been debated for decades, and several aspects remain unknown. For instance, the spatial distribution of crustal shortening that ultimately sustains Himalayan topography and the activity of major fault zones remain unknown at Ma timescales. In this study, we address the spatial deformation pattern in the data-scarce western Himalaya. We calculated catchment averaged, normalized river-steepness indices of non-glaciated drainage basins with tributary catchment areas between 5 and 200 km(2) (n = 2138). We analyzed the spatial distribution of the relative change of river steepness both along and across strike to gain information about the regional distribution of differential uplift pattern and relate this to the activity of distinctive fault segments. For our study area, we observe a positive correlation of averaged k(sn) values with long-term exhumation rates derived from previously published thermochronologic datasets combined with thermal modeling as well as with millennial timescale denudation rates based on cosmogenic nuclide dating. Our results indicate three tectono-geomorphic segments with distinctive landscape morphology, structural architecture, and fault geometry along the western Himalaya: Garhwal-Sutlej, Chamba, and Kashmir Himalaya (from east to west). Moreover, our data recognize distinctive fault segments showing varying thrust activity along strike of the Main Frontal Thrust, the Main Boundary Thrust, and in the vicinity of the steep topographic transition between the Lesser and Greater Himalaya. In this region, we relate out-of-sequence deformation along major basement thrust ramps, such as the Munsiari Thrust with deformation along a mid-crustal ramp along the basal decollement. We suggest that during the Quaternary, all major fault zones in the Western Himalaya experienced out-of-sequence faulting and have accommodated some portion of crustal shortening.}, language = {en} } @article{GrujicGovinBarrieretal.2018, author = {Grujic, Djordje and Govin, Gwladys and Barrier, Laurie and Bookhagen, Bodo and Coutand, Isabelle and Cowan, Beth and Hren, Michael T. and Najman, Yani}, title = {Formation of a Rain Shadow}, series = {Geochemistry, geophysics, geosystems}, volume = {19}, journal = {Geochemistry, geophysics, geosystems}, number = {9}, publisher = {American Geophysical Union}, address = {Washington}, issn = {1525-2027}, doi = {10.1029/2017GC007254}, pages = {3430 -- 3447}, year = {2018}, abstract = {We measure the oxygen and hydrogen stable isotope composition of authigenic clays from Himalayan foreland sediments (Siwalik Group), and from present day small stream waters in eastern Bhutan to explore the impact of uplift of the Shillong Plateau on rain shadow formation over the Himalayan foothills. Stable isotope data from authigenic clay minerals (<2 μm) suggest the presence of three paleoclimatic periods during deposition of the Siwalik Group, between ∼7 and ∼1 Ma. The mean δ18O value in paleometeoric waters, which were in equilibrium with clay minerals, is ∼2.5 per mille lower than in modern meteoric and stream waters at the elevation of the foreland basin. We discuss the factors that could have changed the isotopic composition of water over time and we conclude that (a) the most likely and significant cause for the increase in meteoric water δ18O values over time is the "amount effect," specifically, a decrease in mean annual precipitation. (b) The change in mean annual precipitation over the foreland basin and foothills of the Himalaya is the result of orographic effect caused by the Shillong Plateau's uplift. The critical elevation of the Shillong Plateau required to induce significant orographic precipitation was attained after ∼1.2 Ma. (c) By applying scale analysis, we estimate that the mean annual precipitation over the foreland basin of the eastern Bhutan Himalayas has decreased by a factor of 1.7-2.5 over the last 1-3 million years.}, language = {en} }