@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} } @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} } @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{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} } @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} } @article{JobeLiBookhagenetal.2018, author = {Jobe, Jessica Ann Thompson and Li, Tao and Bookhagen, Bodo and Chen, Jie and Burbank, Douglas W.}, title = {Dating growth strata and basin fill by combining Al-26/Be-10 burial dating and magnetostratigraphy}, series = {Lithosphere}, volume = {10}, journal = {Lithosphere}, number = {6}, publisher = {American Institute of Physics}, address = {Boulder}, issn = {1941-8264}, doi = {10.1130/L727.1}, pages = {806 -- 828}, year = {2018}, abstract = {Cosmogenic burial dating enables dating of coarse-grained, Pliocene-Pleistocene sedimentary units that are typically difficult to date with traditional methods, such as magnetostratigraphy. In the actively deforming western Tarim Basin in NW China, Pliocene-Pleistocene conglomerates were dated at eight sites, integrating Al-26/Be-10 burial dating with previously published magnetostratigraphic sections. These samples were collected from growth strata on the flanks of growing folds and from sedimentary units beneath active faults to place timing constraints on the initiation of deformation of structures within the basin and on shortening rates on active faults. These new basin-fill and growthstrata ages document the late Neogene and Quaternary growth of the Pamir and Tian Shan orogens between >5 and 1 Ma and delineate the eastward propagation of deformation at rates up to 115 km/m.y. and basinward growth of both mountain belts at rates up to 12 km/m.y.}, 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{SmithBookhagenCannon2015, author = {Smith, Taylor and Bookhagen, Bodo and Cannon, Forest}, title = {Improving semi-automated glacier mapping with a multi-method approach: applications in central Asia}, series = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, volume = {9}, journal = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, number = {5}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1994-0416}, doi = {10.5194/tc-9-1747-2015}, pages = {1747 -- 1759}, year = {2015}, abstract = {Studies of glaciers generally require precise glacier outlines. Where these are not available, extensive manual digitization in a geographic information system (GIS) must be performed, as current algorithms struggle to delineate glacier areas with debris cover or other irregular spectral profiles. Although several approaches have improved upon spectral band ratio delineation of glacier areas, none have entered wide use due to complexity or computational intensity. In this study, we present and apply a glacier mapping algorithm in Central Asia which delineates both clean glacier ice and debris-covered glacier tongues. The algorithm is built around the unique velocity and topographic characteristics of glaciers and further leverages spectral and spatial relationship data. We found that the algorithm misclassifies between 2 and 10\% of glacier areas, as compared to a similar to 750 glacier control data set, and can reliably classify a given Landsat scene in 3-5 min. The algorithm does not completely solve the difficulties inherent in classifying glacier areas from remotely sensed imagery but does represent a significant improvement over purely spectral-based classification schemes, such as the band ratio of Landsat 7 bands three and five or the normalized difference snow index. The main caveats of the algorithm are (1) classification errors at an individual glacier level, (2) reliance on manual intervention to separate connected glacier areas, and (3) dependence on fidelity of the input Landsat data.}, language = {en} } @article{ThompsonBurbankLietal.2015, author = {Thompson, Jessica A. and Burbank, Douglas W. and Li, Tao and Chen, Jie and Bookhagen, Bodo}, title = {Late Miocene northward propagation of the northeast Pamir thrust system, northwest China}, series = {Tectonics}, volume = {34}, journal = {Tectonics}, number = {3}, publisher = {American Geophysical Union}, address = {Washington}, issn = {0278-7407}, doi = {10.1002/2014TC003690}, pages = {510 -- 534}, year = {2015}, abstract = {Piggyback basins on the margins of growing orogens commonly serve as sensitive recorders of the onset of thrust deformation and changes in source areas. The Bieertuokuoyi piggyback basin, located in the hanging wall of the Pamir Frontal Thrust, provides an unambiguous record of the outward growth of the northeast Pamir margin in northwest China from the Miocene through the Quaternary. To reconstruct the deformation along the margin, we synthesized structural mapping, stratigraphy, magnetostratigraphy, and cosmogenic burial dating of basin fill and growth strata. The Bieertuokuoyi basin records the initiation of the Pamir Frontal Thrust and the Takegai Thrust similar to 5-6Ma, as well as clast provenance and paleocurrent changes resulting from the Pliocene-to-Recent uplift and exhumation of the Pamir to the south. Our results show that coeval deformation was accommodated on the major structures on the northeast Pamir margin throughout the Miocene to Recent. Furthermore, our data support a change in the regional kinematics around the Miocene-Pliocene boundary (similar to 5-6Ma). Rapid exhumation of NE Pamir extensional domes, coupled with cessation of the Kashgar-Yecheng Transfer System on the eastern margin of the Pamir, accelerated the outward propagation of the northeastern Pamir margin and the southward propagation of the Kashi-Atushi fold-and-thrust belt in the southern Tian Shan. This coeval deformation signifies the coupling of the Pamir and Tarim blocks and the transfer of shortening north to the Pamir frontal faults and across the quasi-rigid Tarim Basin to the southern Tian Shan Kashi-Atushi fold-and-thrust system.}, language = {en} } @article{BoersBarbosaBookhagenetal.2015, author = {Boers, Niklas and Barbosa, Henrique M. J. and Bookhagen, Bodo and Marengo, Jose A. and Marwan, Norbert and Kurths, J{\"u}rgen}, title = {Propagation of Strong Rainfall Events from Southeastern South America to the Central Andes}, series = {Journal of climate}, volume = {28}, journal = {Journal of climate}, number = {19}, publisher = {American Meteorological Soc.}, address = {Boston}, issn = {0894-8755}, doi = {10.1175/JCLI-D-15-0137.1}, pages = {7641 -- 7658}, year = {2015}, abstract = {Based on high-spatiotemporal-resolution data, the authors perform a climatological study of strong rainfall events propagating from southeastern South America to the eastern slopes of the central Andes during the monsoon season. These events account for up to 70\% of total seasonal rainfall in these areas. They are of societal relevance because of associated natural hazards in the form of floods and landslides, and they form an intriguing climatic phenomenon, because they propagate against the direction of the low-level moisture flow from the tropics. The responsible synoptic mechanism is analyzed using suitable composites of the relevant atmospheric variables with high temporal resolution. The results suggest that the low-level inflow from the tropics, while important for maintaining sufficient moisture in the area of rainfall, does not initiate the formation of rainfall clusters. Instead, alternating low and high pressure anomalies in midlatitudes, which are associated with an eastward-moving Rossby wave train, in combination with the northwestern Argentinean low, create favorable pressure and wind conditions for frontogenesis and subsequent precipitation events propagating from southeastern South America toward the Bolivian Andes.}, language = {en} } @article{StolbovaMartinBookhagenetal.2014, author = {Stolbova, Veronika and Martin, P. and Bookhagen, Bodo and Marwan, Norbert and Kurths, J{\"u}rgen}, title = {Topology and seasonal evolution of the network of extreme precipitation over the Indian subcontinent and Sri Lanka}, series = {Nonlinear processes in geophysics}, volume = {21}, journal = {Nonlinear processes in geophysics}, number = {4}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1023-5809}, doi = {10.5194/npg-21-901-2014}, pages = {901 -- 917}, year = {2014}, abstract = {This paper employs a complex network approach to determine the topology and evolution of the network of extreme precipitation that governs the organization of extreme rainfall before, during, and after the Indian Summer Monsoon (ISM) season. We construct networks of extreme rainfall events during the ISM (June-September), post-monsoon (October-December), and pre-monsoon (March-May) periods from satellite-derived (Tropical Rainfall Measurement Mission, TRMM) and rain-gauge interpolated (Asian Precipitation Highly Resolved Observational Data Integration Towards the Evaluation of Water Resources, APHRODITE) data sets. The structure of the networks is determined by the level of synchronization of extreme rainfall events between different grid cells throughout the Indian subcontinent. Through the analysis of various complex-network metrics, we describe typical repetitive patterns in North Pakistan (NP), the Eastern Ghats (EG), and the Tibetan Plateau (TP). These patterns appear during the pre-monsoon season, evolve during the ISM, and disappear during the post-monsoon season. These are important meteorological features that need further attention and that may be useful in ISM timing and strength prediction.}, language = {en} } @article{HarveyBurbankBookhagen2015, author = {Harvey, Jonathan E. and Burbank, Douglas W. and Bookhagen, Bodo}, title = {Along-strike changes in Himalayan thrust geometry: Topographic and tectonic discontinuities in western Nepal}, series = {Lithosphere}, volume = {7}, journal = {Lithosphere}, number = {5}, publisher = {American Institute of Physics}, address = {Boulder}, issn = {1941-8264}, doi = {10.1130/L444.1}, pages = {511 -- 518}, year = {2015}, abstract = {Geodetic and seismologic studies support a tectonic model for the central Himalaya wherein similar to 2 cm/yr of Indo-Asian convergence is accommodated along the primary decollement under the range, the Main Himalayan thrust. A steeper midcrustal ramp in the Main Himalayan thrust is commonly invoked as driving rapid rock uplift along a range-parallel band in the Greater Himalaya. This tectonic model, developed primarily from studies in central Nepal, is commonly assumed to project along strike with little lateral variation in Main Himalayan thrust geometry or associated rock uplift patterns. Here, we synthesize multiple lines of evidence for a major discontinuity in the Main Himalayan thrust in western Nepal. Analysis of topography and seismicity indicates that west of similar to 82.5 degrees E, the single band of steep topography and seismicity along the Main Himalayan thrust ramp in central Nepal bifurcates around a high-elevation, low-relief landscape, resulting in a two-step topographic front along an similar to 150 km segment of the central Himalaya. Although multiple models could explain this bifurcation, the full suite of data appears to be most consistent with a northward bend to the Main Himalayan thrust ramp and activation of a young duplex horse to the south. This poorly documented segmentation of the Main Himalayan thrust has important implications for the seismogenic potential of the western Nepal seismic gap and for models of the ongoing evolution of the orogen.}, 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{HartmanBookhagenChadwick2016, author = {Hartman, Brett D. and Bookhagen, Bodo and Chadwick, Oliver A.}, title = {The effects of check dams and other erosion control structures on the restoration of Andean bofedal ecosystems}, series = {Restoration Ecology}, volume = {24}, journal = {Restoration Ecology}, publisher = {Wiley-Blackwell}, address = {Hoboken}, issn = {1061-2971}, doi = {10.1111/rec.12402}, pages = {761 -- 772}, year = {2016}, abstract = {Restoring degraded lands in rural environments that are heavily managed to meet subsistence needs is a challenge due to high rates of disturbance and resource extraction. This study investigates the efficacy of erosion control structures (ECSs) as restoration tools in the context of a watershed rehabilitation and wet meadow (bofedal) restoration program in the Bolivian Andes. In an effort to enhance water security and increase grazing stability, Aymara indigenous communities built over 15,000 check dams, 9,100 terraces, 5,300 infiltration ditches, and 35 pasture improvement trials. Communities built ECSs at different rates, and we compared vegetation change in the highest restoration management intensity, lowest restoration management intensity, and nonproject control communities. We used line transects to measure changes in vegetation cover and standing water in gullies with check dams and without check dams, and related these ground measurements to a time series (1986-2009) of normalized difference vegetation index derived from Landsat TM5 images. Evidence suggests that check dams increase bofedal vegetation and standing water at a local scale, and lead to increased greenness at a basin scale when combined with other ECSs. Watershed rehabilitation enhances ecosystem services significant to local communities (grazing stability, water security), which creates important synergies when conducting land restoration in rural development settings.}, 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} } @article{HoffmannFeakinsBookhagenetal.2016, author = {Hoffmann, Bernd and Feakins, Sarah J. and Bookhagen, Bodo and Olen, Stephanie M. and Adhikari, Danda P. and Mainali, Janardan and Sachse, Dirk}, title = {Climatic and geomorphic drivers of plant organic matter transport in the Arun River, E Nepal}, series = {Earth \& planetary science letters}, volume = {452}, journal = {Earth \& planetary science letters}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0012-821X}, doi = {10.1016/j.epsl.2016.07.008}, pages = {104 -- 114}, year = {2016}, language = {en} } @misc{JobeLiBookhagenetal.2018, author = {Jobe, Jessica Ann Thompson and Li, Tao and Bookhagen, Bodo and Chen, Jie and Burbank, Douglas W.}, title = {Dating growth strata and basin fill by combining 26Al/10Be burial dating and magnetostratigraphy}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {1044}, issn = {1866-8372}, doi = {10.25932/publishup-46806}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-468067}, pages = {806 -- 828}, year = {2018}, abstract = {Cosmogenic burial dating enables dating of coarse-grained, Pliocene-Pleistocene sedimentary units that are typically difficult to date with traditional methods, such as magnetostratigraphy. In the actively deforming western Tarim Basin in NW China, Pliocene-Pleistocene conglomerates were dated at eight sites, integrating Al-26/Be-10 burial dating with previously published magnetostratigraphic sections. These samples were collected from growth strata on the flanks of growing folds and from sedimentary units beneath active faults to place timing constraints on the initiation of deformation of structures within the basin and on shortening rates on active faults. These new basin-fill and growthstrata ages document the late Neogene and Quaternary growth of the Pamir and Tian Shan orogens between >5 and 1 Ma and delineate the eastward propagation of deformation at rates up to 115 km/m.y. and basinward growth of both mountain belts at rates up to 12 km/m.y.}, 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{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 Eric and Thieken, Annegret Henriette 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} } @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{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} } @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} } @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} } @misc{TostEhmelHeidmannetal.2018, author = {Tost, Jordi and Ehmel, Fabian and Heidmann, Frank and Olen, Stephanie M. and Bookhagen, Bodo}, title = {Hazards and accessibility}, series = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, number = {710}, issn = {1866-8372}, doi = {10.25932/publishup-42785}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-427853}, pages = {8}, year = {2018}, abstract = {The assessment of natural hazards and risk has traditionally been built upon the estimation of threat maps, which are used to depict potential danger posed by a particular hazard throughout a given area. But when a hazard event strikes, infrastructure is a significant factor that can determine if the situation becomes a disaster. The vulnerability of the population in a region does not only depend on the area's local threat, but also on the geographical accessibility of the area. This makes threat maps by themselves insufficient for supporting real-time decision-making, especially for those tasks that involve the use of the road network, such as management of relief operations, aid distribution, or planning of evacuation routes, among others. To overcome this problem, this paper proposes a multidisciplinary approach divided in two parts. First, data fusion of satellite-based threat data and open infrastructure data from OpenStreetMap, introducing a threat-based routing service. Second, the visualization of this data through cartographic generalization and schematization. This emphasizes critical areas along roads in a simple way and allows users to visually evaluate the impact natural hazards may have on infrastructure. We develop and illustrate this methodology with a case study of landslide threat for an area in Colombia.}, language = {en} } @article{SmithBookhagenRheinwalt2017, author = {Smith, Taylor and Bookhagen, Bodo and Rheinwalt, Aljoscha}, title = {Spatiotemporal patterns of High Mountain Asia's snowmelt season identified with an automated snowmelt detection algorithm, 1987-2016}, series = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, volume = {11}, journal = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, issn = {1994-0416}, doi = {10.5194/tc-11-2329-2017}, pages = {2329 -- 2343}, year = {2017}, abstract = {High Mountain Asia (HMA) - encompassing the Tibetan Plateau and surrounding mountain ranges - is the primary water source for much of Asia, serving more than a billion downstream users. Many catchments receive the majority of their yearly water budget in the form of snow, which is poorly monitored by sparse in situ weather networks. Both the timing and volume of snowmelt play critical roles in downstream water provision, as many applications - such as agriculture, drinking-water generation, and hydropower - rely on consistent and predictable snowmelt runoff. Here, we examine passive microwave data across HMA with five sensors (SSMI, SSMIS, AMSR-E, AMSR2, and GPM) from 1987 to 2016 to track the timing of the snowmelt season - defined here as the time between maximum passive microwave signal separation and snow clearance. We validated our method against climate model surface temperatures, optical remote-sensing snow-cover data, and a manual control dataset (n = 2100, 3 variables at 25 locations over 28 years); our algorithm is generally accurate within 3-5 days. Using the algorithm-generated snowmelt dates, we examine the spatiotemporal patterns of the snowmelt season across HMA. The climatically short (29-year) time series, along with complex interannual snowfall variations, makes determining trends in snowmelt dates at a single point difficult. We instead identify trends in snowmelt timing by using hierarchical clustering of the passive microwave data to determine trends in self-similar regions. We make the following four key observations. (1) The end of the snowmelt season is trending almost universally earlier in HMA (negative trends). Changes in the end of the snowmelt season are generally between 2 and 8 days decade 1 over the 29-year study period (5-25 days total). The length of the snowmelt season is thus shrinking in many, though not all, regions of HMA. Some areas exhibit later peak signal separation (positive trends), but with generally smaller magnitudes than trends in snowmelt end. (2) Areas with long snowmelt periods, such as the Tibetan Plateau, show the strongest compression of the snowmelt season (negative trends). These trends are apparent regardless of the time period over which the regression is performed. (3) While trends averaged over 3 decades indicate generally earlier snowmelt seasons, data from the last 14 years (2002-2016) exhibit positive trends in many regions, such as parts of the Pamir and Kunlun Shan. Due to the short nature of the time series, it is not clear whether this change is a reversal of a long-term trend or simply interannual variability. (4) Some regions with stable or growing glaciers - such as the Karakoram and Kunlun Shan - see slightly later snowmelt seasons and longer snowmelt periods. It is likely that changes in the snowmelt regime of HMA account for some of the observed heterogeneity in glacier response to climate change. While the decadal increases in regional temperature have in general led to earlier and shortened melt seasons, changes in HMA's cryosphere have been spatially and temporally heterogeneous.}, language = {en} } @article{PurintonBookhagen2017, author = {Purinton, Benjamin and Bookhagen, Bodo}, title = {Validation of digital elevation models (DEMs) and comparison of geomorphic metrics on the southern Central Andean Plateau}, series = {Earth surface dynamics}, volume = {5}, journal = {Earth surface dynamics}, number = {2}, publisher = {Copernicus Publ.}, address = {G{\"o}ttingen}, issn = {2196-632X}, doi = {10.5194/esurf-5-211-2017}, pages = {211 -- 237}, year = {2017}, abstract = {In this study, we validate and compare elevation accuracy and geomorphic metrics of satellite-derived digital elevation models (DEMs) on the southern Central Andean Plateau. The plateau has an average elevation of 3.7 km and is characterized by diverse topography and relief, lack of vegetation, and clear skies that create ideal conditions for remote sensing. At 30m resolution, SRTM-C, ASTER GDEM2, stacked ASTER L1A stereopair DEM, ALOS World 3D, and TanDEM-X have been analyzed. The higher-resolution datasets include 12m TanDEM-X, 10m single-CoSSC TerraSAR-X/TanDEM-X DEMs, and 5m ALOS World 3D. These DEMs are state of the art for optical (ASTER and ALOS) and radar (SRTM-C and TanDEM-X) spaceborne sensors. We assessed vertical accuracy by comparing standard deviations of the DEM elevation versus 307 509 differential GPS measurements across 4000m of elevation. For the 30m DEMs, the ASTER datasets had the highest vertical standard deviation at > 6.5 m, whereas the SRTM-C, ALOS World 3D, and TanDEM-X were all < 3.5 m. Higher-resolution DEMs generally had lower uncertainty, with both the 12m TanDEM-X and 5m ALOSWorld 3D having < 2m vertical standard deviation. Analysis of vertical uncertainty with respect to terrain elevation, slope, and aspect revealed the low uncertainty across these attributes for SRTM-C (30 m), TanDEM-X (12-30 m), and ALOS World 3D (5-30 m). Single-CoSSC TerraSAR-X/TanDEM-X 10m DEMs and the 30m ASTER GDEM2 displayed slight aspect biases, which were removed in their stacked counterparts (TanDEM-X and ASTER Stack). Based on low vertical standard deviations and visual inspection alongside optical satellite data, we selected the 30m SRTM-C, 12-30m TanDEM-X, 10m single-CoSSC TerraSAR-X/TanDEM-X, and 5m ALOS World 3D for geomorphic metric comparison in a 66 km2 catchment with a distinct river knickpoint. Consistent m=n values were found using chi plot channel profile analysis, regardless of DEM type and spatial resolution. Slope, curvature, and drainage area were calculated and plotting schemes were used to assess basin-wide differences in the hillslope-to-valley transition related to the knickpoint. While slope and hillslope length measurements vary little between datasets, curvature displays higher magnitude measurements with fining resolution. This is especially true for the optical 5m ALOS World 3D DEM, which demonstrated high-frequency noise in 2-8 pixel steps through a Fourier frequency analysis. The improvements in accurate space-radar DEMs (e.g., TanDEM-X) for geomorphometry are promising, but airborne or terrestrial data are still necessary for meter-scale analysis.}, language = {en} } @article{PurintonBookhagen2018, author = {Purinton, Benjamin and Bookhagen, Bodo}, title = {Measuring decadal vertical land-level changes from SRTM-C (2000) and TanDEM-X (∼ 2015) in the south-central Andes}, series = {Earth Surface Dynamics}, volume = {6}, journal = {Earth Surface Dynamics}, publisher = {Copernicus Publ.}, address = {G{\"o}ttingen}, issn = {2196-6311}, doi = {10.5194/esurf-6-971-2018}, pages = {971 -- 987}, year = {2018}, abstract = {In the arctic and high mountains it is common to measure vertical changes of ice sheets and glaciers via digital elevation model (DEM) differencing. This requires the signal of change to outweigh the noise associated with the datasets. Excluding large landslides, on the ice-free earth the land-level change is smaller in vertical magnitude and thus requires more accurate DEMs for differencing and identification of change. Previously, this has required meter to submeter data at small spatial scales. Following careful corrections, we are able to measure land-level changes in gravel-bed channels and steep hillslopes in the south-central Andes using the SRTM-C (collected in 2000) and the TanDEM-X (collected from 2010 to 2015) near-global 12-30m DEMs. Long-standing errors in the SRTM-C are corrected using the TanDEM-X as a control surface and applying cosine-fit co-registration to remove ∼ 1∕10 pixel (∼ 3m) shifts, fast Fourier transform (FFT) and filtering to remove SRTM-C short- and long-wavelength stripes, and blocked shifting to remove remaining complex biases. The datasets are then differenced and outlier pixels are identified as a potential signal for the case of gravel-bed channels and hillslopes. We are able to identify signals of incision and aggradation (with magnitudes down to ∼ 3m in the best case) in two  > 100km river reaches, with increased geomorphic activity downstream of knickpoints. Anthropogenic gravel excavation and piling is prominently measured, with magnitudes exceeding ±5m (up to  > 10m for large piles). These values correspond to conservative average rates of 0.2 to > 0.5myr-1 for vertical changes in gravel-bed rivers. For hillslopes, since we require stricter cutoffs for noise, we are only able to identify one major landslide in the study area with a deposit volume of 16±0.15×106m3. Additional signals of change can be garnered from TanDEM-X auxiliary layers; however, these are more difficult to quantify. The methods presented can be extended to any region of the world with SRTM-C and TanDEM-X coverage where vertical land-level changes are of interest, with the caveat that remaining vertical uncertainties in primarily the SRTM-C limit detection in steep and complex topography.}, language = {en} } @article{AlonzoBookhagenMcFaddenetal.2015, author = {Alonzo, Michael and Bookhagen, Bodo and McFadden, Joseph P. and Sun, Alex and Roberts, Dar A.}, title = {Mapping urban forest leaf area index with airborne lidar using penetration metrics and allometry}, series = {Remote sensing of environment : an interdisciplinary journal}, volume = {162}, journal = {Remote sensing of environment : an interdisciplinary journal}, publisher = {Elsevier}, address = {New York}, issn = {0034-4257}, doi = {10.1016/j.rse.2015.02.025}, pages = {141 -- 153}, year = {2015}, abstract = {In urban areas, leaf area index (LAI) is a key ecosystem structural attribute with implications for energy and water balance, gas exchange, and anthropogenic energy use. In this study, we estimated LAI spatially using airborne lidar in downtown Santa Barbara, California, USA. We implemented two different modeling approaches. First, we directly estimated effective LAI (LAIe) using scan angle- and clump-corrected lidar laser penetration metrics (LPM). Second, we adapted existing allometric equations to estimate crown structural metrics including tree height and crown base height using lidar. The latter approach allowed for LAI estimates at the individual tree-crown scale. The LPM method, at both high and decimated point densities, resulted in good linear agreement with estimates from ground-based hemispherical photography (r(2) = 0.82, y = 0.99x) using a model that assumed a spherical leaf angle distribution. Within individual tree crown segments, the lidar estimates of crown structure closely paralleled field measurements (e.g., r(2) = 0.87 for crown length). LAI estimates based on the lidar crown measurements corresponded well with estimates from field measurements (r(2) = 0.84, y = 0.97x + 0.10). Consistency of the LPM and allometric lidar methods was also strong at 71 validation plots (r(2) = 0.88) and at 450 additional sample locations across the entire study area (r(2) = 0.72). This level of correspondence exceeded that of the canopy hemispherical photography and allometric, ground-based estimates (r(2) = 0.53). The first-order alignment of these two disparate methods may indicate that the error bounds for mapping LAI in cities are small enough to pursue large scale, spatially explicit estimation. (C) 2015 Elsevier Inc All rights reserved.}, language = {en} } @article{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}, volume = {10}, journal = {remote sensing}, number = {8}, publisher = {Molecular Diversity Preservation International (MDPI)}, address = {Basel}, issn = {2072-4292}, doi = {10.3390/rs10081272}, pages = {1 -- 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} } @article{BoersBookhagenMarengoetal.2015, author = {Boers, Niklas and Bookhagen, Bodo and Marengo, Jose and Marwan, Norbert and von Storch, Jin-Song and Kurths, J{\"u}rgen}, title = {Extreme Rainfall of the South American Monsoon System: A Dataset Comparison Using Complex Networks}, series = {Journal of climate}, volume = {28}, journal = {Journal of climate}, number = {3}, publisher = {American Meteorological Soc.}, address = {Boston}, issn = {0894-8755}, doi = {10.1175/JCLI-D-14-00340.1}, pages = {1031 -- 1056}, year = {2015}, abstract = {In this study, the authors compare six different rainfall datasets for South America with a focus on their representation of extreme rainfall during the monsoon season (December February): the gauge-calibrated TRMM 3B42 V7 satellite product; the near-real-time TRMM 3B42 V7 RT, the GPCP 1 degrees daily (1DD) V1.2 satellite gauge combination product, the Interim ECMWF Re-Analysis (ERA-Interim) product; output of a high-spatial-resolution run of the ECHAM6 global circulation model; and output of the regional climate model Eta. For the latter three, this study can be understood as a model evaluation. In addition to statistical values of local rainfall distributions, the authors focus on the spatial characteristics of extreme rainfall covariability. Since traditional approaches based on principal component analysis are not applicable in the context of extreme events, they apply and further develop methods based on complex network theory. This way, the authors uncover substantial differences in extreme rainfall patterns between the different datasets: (i) The three model-derived datasets yield very different results than the satellite gauge combinations regarding the main climatological propagation pathways of extreme events as well as the main convergence zones of the monsoon system. (ii) Large discrepancies are found for the development of mesoscale convective systems in southeastern South America. (iii) Both TRMM datasets and ECHAM6 indicate a linkage of extreme rainfall events between the central Amazon basin and the eastern slopes of the central Andes, but this pattern is not reproduced by the remaining datasets. The authors' study suggests that none of the three model-derived datasets adequately captures extreme rainfall patterns in South America.}, language = {en} } @misc{PurintonBookhagen2018, author = {Purinton, Benjamin and Bookhagen, Bodo}, title = {Measuring decadal vertical land-level changes from SRTM-C (2000) and TanDEM-X (∼ 2015) 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 = {480}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-420487}, pages = {16}, year = {2018}, abstract = {In the arctic and high mountains it is common to measure vertical changes of ice sheets and glaciers via digital elevation model (DEM) differencing. This requires the signal of change to outweigh the noise associated with the datasets. Excluding large landslides, on the ice-free earth the land-level change is smaller in vertical magnitude and thus requires more accurate DEMs for differencing and identification of change. Previously, this has required meter to submeter data at small spatial scales. Following careful corrections, we are able to measure land-level changes in gravel-bed channels and steep hillslopes in the south-central Andes using the SRTM-C (collected in 2000) and the TanDEM-X (collected from 2010 to 2015) near-global 12-30m DEMs. Long-standing errors in the SRTM-C are corrected using the TanDEM-X as a control surface and applying cosine-fit co-registration to remove ∼ 1∕10 pixel (∼ 3m) shifts, fast Fourier transform (FFT) and filtering to remove SRTM-C short- and long-wavelength stripes, and blocked shifting to remove remaining complex biases. The datasets are then differenced and outlier pixels are identified as a potential signal for the case of gravel-bed channels and hillslopes. We are able to identify signals of incision and aggradation (with magnitudes down to ∼ 3m in the best case) in two  > 100km river reaches, with increased geomorphic activity downstream of knickpoints. Anthropogenic gravel excavation and piling is prominently measured, with magnitudes exceeding ±5m (up to  > 10m for large piles). These values correspond to conservative average rates of 0.2 to > 0.5myr-1 for vertical changes in gravel-bed rivers. For hillslopes, since we require stricter cutoffs for noise, we are only able to identify one major landslide in the study area with a deposit volume of 16±0.15×106m3. Additional signals of change can be garnered from TanDEM-X auxiliary layers; however, these are more difficult to quantify. The methods presented can be extended to any region of the world with SRTM-C and TanDEM-X coverage where vertical land-level changes are of interest, with the caveat that remaining vertical uncertainties in primarily the SRTM-C limit detection in steep and complex topography.}, language = {en} } @misc{SmithBookhagenRheinwalt2017, author = {Smith, Taylor and Bookhagen, Bodo and Rheinwalt, Aljoscha}, title = {Spatiotemporal patterns of High Mountain Asia's snowmelt season identified with an automated snowmelt detection algorithm, 1987-2016}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-403911}, pages = {15}, year = {2017}, abstract = {High Mountain Asia (HMA) - encompassing the Tibetan Plateau and surrounding mountain ranges - is the primary water source for much of Asia, serving more than a billion downstream users. Many catchments receive the majority of their yearly water budget in the form of snow, which is poorly monitored by sparse in situ weather networks. Both the timing and volume of snowmelt play critical roles in downstream water provision, as many applications - such as agriculture, drinking-water generation, and hydropower - rely on consistent and predictable snowmelt runoff. Here, we examine passive microwave data across HMA with five sensors (SSMI, SSMIS, AMSR-E, AMSR2, and GPM) from 1987 to 2016 to track the timing of the snowmelt season - defined here as the time between maximum passive microwave signal separation and snow clearance. We validated our method against climate model surface temperatures, optical remote-sensing snow-cover data, and a manual control dataset (n = 2100, 3 variables at 25 locations over 28 years); our algorithm is generally accurate within 3-5 days. Using the algorithm-generated snowmelt dates, we examine the spatiotemporal patterns of the snowmelt season across HMA. The climatically short (29-year) time series, along with complex interannual snowfall variations, makes determining trends in snowmelt dates at a single point difficult. We instead identify trends in snowmelt timing by using hierarchical clustering of the passive microwave data to determine trends in self-similar regions. We make the following four key observations. (1) The end of the snowmelt season is trending almost universally earlier in HMA (negative trends). Changes in the end of the snowmelt season are generally between 2 and 8 days decade 1 over the 29-year study period (5-25 days total). The length of the snowmelt season is thus shrinking in many, though not all, regions of HMA. Some areas exhibit later peak signal separation (positive trends), but with generally smaller magnitudes than trends in snowmelt end. (2) Areas with long snowmelt periods, such as the Tibetan Plateau, show the strongest compression of the snowmelt season (negative trends). These trends are apparent regardless of the time period over which the regression is performed. (3) While trends averaged over 3 decades indicate generally earlier snowmelt seasons, data from the last 14 years (2002-2016) exhibit positive trends in many regions, such as parts of the Pamir and Kunlun Shan. Due to the short nature of the time series, it is not clear whether this change is a reversal of a long-term trend or simply interannual variability. (4) Some regions with stable or growing glaciers - such as the Karakoram and Kunlun Shan - see slightly later snowmelt seasons and longer snowmelt periods. It is likely that changes in the snowmelt regime of HMA account for some of the observed heterogeneity in glacier response to climate change. While the decadal increases in regional temperature have in general led to earlier and shortened melt seasons, changes in HMA's cryosphere have been spatially and temporally heterogeneous.}, language = {en} } @article{PurintonBookhagen2017, author = {Purinton, Benjamin and Bookhagen, Bodo}, title = {Validation of digital elevation models (DEMs) and comparison of geomorphic metrics on the southern Central Andean Plateau}, series = {Earth surface dynamics}, volume = {5}, journal = {Earth surface dynamics}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {2196-6311}, doi = {10.5194/esurf-5-211-2017}, pages = {211 -- 237}, year = {2017}, abstract = {In this study, we validate and compare elevation accuracy and geomorphic metrics of satellite-derived digital elevation models (DEMs) on the southern Central Andean Plateau. The plateau has an average elevation of 3.7 km and is characterized by diverse topography and relief, lack of vegetation, and clear skies that create ideal conditions for remote sensing. At 30m resolution, SRTM-C, ASTER GDEM2, stacked ASTER L1A stereopair DEM, ALOS World 3D, and TanDEM-X have been analyzed. The higher-resolution datasets include 12m TanDEM-X, 10m single-CoSSC TerraSAR-X/TanDEM-X DEMs, and 5m ALOS World 3D. These DEMs are state of the art for optical (ASTER and ALOS) and radar (SRTM-C and TanDEM-X) spaceborne sensors. We assessed vertical accuracy by comparing standard deviations of the DEM elevation versus 307 509 differential GPS measurements across 4000m of elevation. For the 30m DEMs, the ASTER datasets had the highest vertical standard deviation at > 6.5 m, whereas the SRTM-C, ALOS World 3D, and TanDEM-X were all < 3.5 m. Higher-resolution DEMs generally had lower uncertainty, with both the 12m TanDEM-X and 5m ALOSWorld 3D having < 2m vertical standard deviation. Analysis of vertical uncertainty with respect to terrain elevation, slope, and aspect revealed the low uncertainty across these attributes for SRTM-C (30 m), TanDEM-X (12-30 m), and ALOS World 3D (5-30 m). Single-CoSSC TerraSAR-X/TanDEM-X 10m DEMs and the 30m ASTER GDEM2 displayed slight aspect biases, which were removed in their stacked counterparts (TanDEM-X and ASTER Stack). Based on low vertical standard deviations and visual inspection alongside optical satellite data, we selected the 30m SRTM-C, 12-30m TanDEM-X, 10m single-CoSSC TerraSAR-X/TanDEM-X, and 5m ALOS World 3D for geomorphic metric comparison in a 66 km2 catchment with a distinct river knickpoint. Consistent m = n values were found using chi plot channel profile analysis, regardless of DEM type and spatial resolution. Slope, curvature, and drainage area were calculated and plotting schemes were used to assess basin-wide differences in the hillslope-to-valley transition related to the knickpoint. While slope and hillslope length measurements vary little between datasets, curvature displays higher magnitude measurements with fining resolution. This is especially true for the optical 5m ALOS World 3D DEM, which demonstrated high-frequency noise in 2-8 pixel steps through a Fourier frequency analysis. The improvements in accurate space-radar DEMs (e. g., TanDEM-X) for geomorphometry are promising, but airborne or terrestrial data are still necessary for meter-scale analysis.}, 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{BrellSeglGuanteretal.2017, author = {Brell, Maximilian and Segl, Karl and Guanter, Luis and Bookhagen, Bodo}, title = {Hyperspectral and Lidar Intensity Data Fusion: A Framework for the Rigorous Correction of Illumination, Anisotropic Effects, and Cross Calibration}, series = {IEEE transactions on geoscience and remote sensing}, volume = {55}, journal = {IEEE transactions on geoscience and remote sensing}, publisher = {Inst. of Electr. and Electronics Engineers}, address = {Piscataway}, issn = {0196-2892}, doi = {10.1109/TGRS.2017.2654516}, pages = {2799 -- 2810}, year = {2017}, abstract = {The fusion of hyperspectral imaging (HSI) sensor and airborne lidar scanner (ALS) data provides promising potential for applications in environmental sciences. Standard fusion approaches use reflectance information from the HSI and distance measurements from the ALS to increase data dimen-sionality and geometric accuracy. However, the potential for data fusion based on the respective intensity information of the complementary active and passive sensor systems is high and not yet fully exploited. Here, an approach for the rigorous illumination correction of HSI data, based on the radiometric cross-calibrated return intensity information of ALS data, is presented. The cross calibration utilizes a ray tracing-based fusion of both sensor measurements by intersecting their particular beam shapes. The developed method is capable of compensating for the drawbacks of passive HSI systems, such as cast and cloud shadowing effects, illumination changes over time, across track illumination, and partly anisotropy effects. During processing, spatial and temporal differences in illumination patterns are detected and corrected over the entire HSI wavelength domain. The improvement in the classification accuracy of urban and vegetation surfaces demonstrates the benefit and potential of the proposed HSI illumination correction. The presented approach is the first step toward the rigorous in-flight fusion of passive and active system characteristics, enabling new capabilities for a variety of applications.}, language = {en} } @misc{PurintonBookhagen2017, author = {Purinton, Benjamin and Bookhagen, Bodo}, title = {Validation of digital elevation models (DEMs) and comparison of geomorphic metrics on the southern Central Andean Plateau}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-396277}, pages = {27}, year = {2017}, abstract = {In this study, we validate and compare elevation accuracy and geomorphic metrics of satellite-derived digital elevation models (DEMs) on the southern Central Andean Plateau. The plateau has an average elevation of 3.7 km and is characterized by diverse topography and relief, lack of vegetation, and clear skies that create ideal conditions for remote sensing. At 30m resolution, SRTM-C, ASTER GDEM2, stacked ASTER L1A stereopair DEM, ALOS World 3D, and TanDEM-X have been analyzed. The higher-resolution datasets include 12m TanDEM-X, 10m single-CoSSC TerraSAR-X/TanDEM-X DEMs, and 5m ALOS World 3D. These DEMs are state of the art for optical (ASTER and ALOS) and radar (SRTM-C and TanDEM-X) spaceborne sensors. We assessed vertical accuracy by comparing standard deviations of the DEM elevation versus 307 509 differential GPS measurements across 4000m of elevation. For the 30m DEMs, the ASTER datasets had the highest vertical standard deviation at > 6.5 m, whereas the SRTM-C, ALOS World 3D, and TanDEM-X were all < 3.5 m. Higher-resolution DEMs generally had lower uncertainty, with both the 12m TanDEM-X and 5m ALOSWorld 3D having < 2m vertical standard deviation. Analysis of vertical uncertainty with respect to terrain elevation, slope, and aspect revealed the low uncertainty across these attributes for SRTM-C (30 m), TanDEM-X (12-30 m), and ALOS World 3D (5-30 m). Single-CoSSC TerraSAR-X/TanDEM-X 10m DEMs and the 30m ASTER GDEM2 displayed slight aspect biases, which were removed in their stacked counterparts (TanDEM-X and ASTER Stack). Based on low vertical standard deviations and visual inspection alongside optical satellite data, we selected the 30m SRTM-C, 12-30m TanDEM-X, 10m single-CoSSC TerraSAR-X/TanDEM-X, and 5m ALOS World 3D for geomorphic metric comparison in a 66 km2 catchment with a distinct river knickpoint. Consistent m=n values were found using chi plot channel profile analysis, regardless of DEM type and spatial resolution. Slope, curvature, and drainage area were calculated and plotting schemes were used to assess basin-wide differences in the hillslope-to-valley transition related to the knickpoint. While slope and hillslope length measurements vary little between datasets, curvature displays higher magnitude measurements with fining resolution. This is especially true for the optical 5m ALOS World 3D DEM, which demonstrated high-frequency noise in 2-8 pixel steps through a Fourier frequency analysis. The improvements in accurate space-radar DEMs (e.g., TanDEM-X) for geomorphometry are promising, but airborne or terrestrial data are still necessary for meter-scale analysis.}, 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{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{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{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} } @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{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} } @incollection{Bookhagen2017, author = {Bookhagen, Bodo}, title = {The influence of hydrology and glaciology on wetlands in the Himalayas}, series = {Bird migration across the Himalayas : wetland functioning amidst mountains and glaciers}, booktitle = {Bird migration across the Himalayas : wetland functioning amidst mountains and glaciers}, editor = {Prins, Herbert H.T. and Namgail, Tsewang}, publisher = {Cambridge University Press}, address = {Cambridge}, isbn = {978-1-107-11471-5}, doi = {10.1017/9781316335420}, pages = {175 -- 188}, year = {2017}, abstract = {Birds migrating across the Himalayan region fly over the highest peaks in the world, facing immense physiological and climatic challenges. The authors show the different strategies used by birds to cope with these challenges. Many wetland avian species are seen in the high-altitude lakes of the Himalayas and the adjoining Tibetan Plateau, such as Bar-Headed Geese. Ringing programmes have generated information about origins and destinations, and this book is the first to present information on the bird's exact migratory paths. Capitalising on knowledge generated through satellite telemetry, the authors describe the migratory routes of a multitude of birds flying over or skirting the Himalayas. The myriad of threats to migratory birds and the wetland system in the Central Asian Flyway are discussed, with ways to mitigate them. This volume will inform and persuade policy-makers and conservation practitioners to take appropriate measures for the long-term survival of this unique migration}, 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{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} }