TY - GEN A1 - Atmani, Farid A1 - Bookhagen, Bodo A1 - Smith, Taylor T1 - Measuring Vegetation Heights and Their Seasonal Changes in the Western Namibian Savanna Using Spaceborne Lidars T2 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe N2 - 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. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 1275 KW - ICESat-2 KW - GEDI KW - canopy height KW - lidar KW - savanna Y1 - 2022 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-569915 SN - 1866-8372 IS - 1275 ER - TY - JOUR A1 - Atmani, Farid A1 - Bookhagen, Bodo A1 - Smith, Taylor T1 - Measuring vegetation heights and their seasonal changes in the Western Namibian Savanna using spaceborne lidars JF - Remote sensing / Molecular Diversity Preservation International (MDPI) N2 - 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. KW - ICESat-2 KW - GEDI KW - canopy height KW - lidar KW - savanna Y1 - 2022 U6 - https://doi.org/10.3390/rs14122928 SN - 2072-4292 VL - 14 IS - 12 SP - 1 EP - 20 PB - MDPI CY - Basel, Schweiz ET - 12 ER - TY - JOUR A1 - Regmi, Shakil A1 - Bookhagen, Bodo T1 - The spatial pattern of extreme precipitation from 40 years of gauge data in the central Himalaya JF - Weather and climate extremes N2 - The topography of the Himalaya exerts a substantial control on the spatial distribution of monsoonal rainfall, which is a vital water source for the regional economy and population. But the occurrence of short-lived and high-intensity precipitation results in socio-economic losses. This study relies on 40 years of daily data from 204 ground stations in Nepal to derive extreme precipitation thresholds, amounts, and days at the 95th percentile. We additionally determine the precipitation magnitude-frequency relation. We observe that extreme precipitation amounts follow an almost uniform band parallel to topographic contour lines in the southern Himalaya mountains in central and eastern Nepal but not in western Nepal. The relationship of extreme precipitation indices with topographic relief shows that extreme precipitation thresholds decrease with increasing elevation, but extreme precipitation days increase in higher elevation areas. Furthermore, stations above 1 km elevation exhibit a power-law relation in the rainfall magnitude-frequency framework. Stations at higher elevations generally have lower values of power-law exponents than low elevation areas. This suggests a fundamentally different behaviour of the rainfall distribution and an increased occurrence of extreme rainfall storms in the high elevation areas of Nepal. KW - Himalaya KW - Nepal KW - Indian summer monsoon KW - Precipitation KW - Extreme KW - precipitation Y1 - 2022 U6 - https://doi.org/10.1016/j.wace.2022.100470 SN - 2212-0947 VL - 37 PB - Elsevier CY - Amsterdam ER - TY - JOUR A1 - Erbello Doelesso, Asfaw A1 - Melnick, Daniel A1 - Zeilinger, Gerold A1 - Bookhagen, Bodo A1 - Pingel, Heiko A1 - Strecker, Manfred T1 - Geomorphic expression of a tectonically active rift-transfer zone in southern Ethiopia JF - Geomorphology : an international journal on pure and applied geomorphology N2 - 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. KW - rift transfer zone KW - Ethiopia rift KW - renya Rift KW - morphometric indices KW - knickpoints Y1 - 2022 U6 - https://doi.org/10.1016/j.geomorph.2022.108162 SN - 0169-555X SN - 1872-695X VL - 403 PB - Elsevier Science CY - Amsterdam [u.a.] ER - TY - JOUR A1 - Peña, Carlos A1 - Metzger, Sabrina A1 - Heidbach, Oliver A1 - Bedford, Jonathan A1 - Bookhagen, Bodo A1 - Moreno, Marcos A1 - Oncken, Onno A1 - Cotton, Fabrice T1 - Role of poroelasticity during the early postseismic deformation of the 2010 Maule megathrust earthquake JF - Geophysical research letters N2 - Megathrust earthquakes impose changes of differential stress and pore pressure in the lithosphere-asthenosphere system that are transiently relaxed during the postseismic period primarily due to afterslip, viscoelastic and poroelastic processes. Especially during the early postseismic phase, however, the relative contribution of these processes to the observed surface deformation is unclear. To investigate this, we use geodetic data collected in the first 48 days following the 2010 Maule earthquake and a poro-viscoelastic forward model combined with an afterslip inversion. This model approach fits the geodetic data 14% better than a pure elastic model. Particularly near the region of maximum coseismic slip, the predicted surface poroelastic uplift pattern explains well the observations. If poroelasticity is neglected, the spatial afterslip distribution is locally altered by up to +/- 40%. Moreover, we find that shallow crustal aftershocks mostly occur in regions of increased postseismic pore-pressure changes, indicating that both processes might be mechanically coupled. KW - Chilean subduction zone KW - poroelasticity KW - power-law rheology KW - afterslip inversion KW - InSAR KW - GNSS Y1 - 2022 U6 - https://doi.org/10.1029/2022GL098144 SN - 0094-8276 SN - 1944-8007 VL - 49 IS - 9 PB - Wiley CY - Hoboken, NJ ER - TY - JOUR A1 - Purinton, Benjamin A1 - Bookhagen, Bodo T1 - Tracking downstream variability in large grain-size distributions in the South-Central Andes JF - Journal of geophysical research : F, Earth surface N2 - 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. KW - digital grain sizing KW - downstream fining KW - fluvial geomorphology KW - grain-size distribution KW - pebblecounts Y1 - 2021 U6 - https://doi.org/10.1029/2021JF006260 SN - 2169-9003 SN - 2169-9011 VL - 126 IS - 8 SP - 1 EP - 29 PB - American Geophysical Union CY - Washington ER - TY - JOUR A1 - Smith, Taylor A1 - Rheinwalt, Aljoscha A1 - Bookhagen, Bodo T1 - Topography and climate in the upper Indus Basin BT - Mapping elevation-snow cover relationships JF - The science of the total environment : an international journal for scientific research into the environment and its relationship with man N2 - 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. KW - Snow-cover KW - Hierarchical clustering KW - Glaciers KW - Upper Indus Basin Y1 - 2021 U6 - https://doi.org/10.1016/j.scitotenv.2021.147363 SN - 0048-9697 SN - 1879-1026 VL - 786 PB - Elsevier CY - Amsterdam ER - TY - JOUR A1 - Smith, Taylor A1 - Bookhagen, Bodo T1 - Climatic and biotic controls on topographic asymmetry at the global scale JF - Journal of geophysical research : JGR, Earth surface N2 - 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. KW - erosion KW - freeze-thaw cycling KW - solar radiation KW - topographic asymmetry KW - topography KW - vegetation cover Y1 - 2021 U6 - https://doi.org/10.1029/2020JF005692 SN - 2169-9003 SN - 2169-9011 VL - 126 IS - 1 PB - American Geophysical Union CY - Washington ER - TY - JOUR A1 - Purinton, Benjamin A1 - Bookhagen, Bodo T1 - Beyond Vertical Point Accuracy BT - Assessing Inter-pixel Consistency in 30 m Global DEMs for the Arid Central Andes JF - Frontiers in Earth Science N2 - 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. KW - DEM noise KW - Fourier analysis KW - TanDEM-X KW - ASTER GDEM KW - Copernicus DEM KW - WorldDEM KW - SRTM KW - ALOS World 3D Y1 - 2021 U6 - https://doi.org/10.3389/feart.2021.758606 SN - 2296-6463 SP - 1 EP - 24 PB - Frontiers Media CY - Lausanne, Schweiz ER - TY - JOUR A1 - Ramezani Ziarani, Maryam A1 - Bookhagen, Bodo A1 - Schmidt, Torsten A1 - Wickert, Jens A1 - de la Torre, Alejandro A1 - Deng, Zhiguo A1 - Calori, Andrea T1 - A model for the relationship between rainfall, GNSS-derived integrated water vapour, and CAPE in the eastern central Andes JF - Remote Sensing N2 - 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. KW - Global Navigation Satellite System (GNSS) KW - GNSS-integrated water vapour KW - convective available potential energy (CAPE) KW - extreme rainfall KW - TRMM Y1 - 2021 U6 - https://doi.org/10.3390/rs13183788 SN - 2072-4292 VL - 13 IS - 18 PB - MDPI CY - Basel ER - TY - GEN A1 - Ramezani Ziarani, Maryam A1 - Bookhagen, Bodo A1 - Schmidt, Torsten A1 - Wickert, Jens A1 - de la Torre, Alejandro A1 - Deng, Zhiguo A1 - Calori, Andrea T1 - A model for the relationship between rainfall, GNSS-derived integrated water vapour, and CAPE in the eastern central Andes T2 - Postprints der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe N2 - 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. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 1172 KW - Global Navigation Satellite System (GNSS) KW - GNSS-integrated water vapour KW - convective available potential energy (CAPE) KW - extreme rainfall KW - TRMM Y1 - 2021 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-523256 SN - 1866-8372 IS - 1172 ER - TY - GEN A1 - Purinton, Benjamin A1 - Bookhagen, Bodo T1 - Beyond Vertical Point Accuracy BT - Assessing Inter-pixel Consistency in 30 m Global DEMs for the Arid Central Andes T2 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe N2 - 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. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 1242 KW - DEM noise KW - Fourier analysis KW - TanDEM-X KW - ASTER GDEM KW - Copernicus DEM KW - WorldDEM KW - SRTM KW - ALOS World 3D Y1 - 2022 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-549805 SN - 1866-8372 SP - 1 EP - 24 PB - Universitätsverlag Potsdam CY - Potsdam ER - TY - JOUR A1 - Müting, Friederike Ariane A1 - Bookhagen, Bodo A1 - Strecker, Manfred T1 - Identification of debris-flow channels using high-resolution topographic data BT - a case study in the Quebrada del Toro, NW Argentina JF - Journal of geophysical research : JGR, Earth surface N2 - Resolving Earth's surface at the meter scale is essential for an improved understanding of the dynamics of mass-movement processes. In this study, we explore the applicability and potential of digital elevation models (DEMs) derived from stereophotogrammetry to detect debris-flow channels in the Quebrada del Toro in the northwestern Argentine Andes. Our analysis relies on a high-resolution (3 m) DEM created from SPOT-7 tri-stereo satellite data. We carefully validated DEM quality with ∼6,000 differential GPS points and identified optimal parameters for DEM generation in high-relief terrain. After multiple processing steps, we achieved an accuracy of 0.051 ± 1.915 m (1σ) using n = 3,139 control points with cm precision. Previous studies have used the drainage area and slope framework to identify topographic signatures of debris flows within a catchment. We built upon this and investigated individual river-channel segments using connected-component (CC) analysis on meter-scale topographic data. We define CC as segments of similar slope along the channel profile. Based on seven manually identified debris-flow catchments, we developed a debris-flow similarity index using component length and mean channel-segment slope and identified channel segments that have likely been shaped by debris flows. The presented approach has the potential to resolve intra-catchment variability of transport processes, allows to constrain the extent of debris-flow channels more precisely than slope-area analysis, and highlights the versatility of combined space- and field-based observations for natural-hazard assessments. Y1 - 2021 U6 - https://doi.org/10.1029/2021JF006330 SN - 2169-9003 SN - 2169-9011 VL - 126 IS - 12 PB - American Geophysical Union CY - Washington ER - TY - JOUR A1 - Olen, Stephanie M. A1 - Bookhagen, Bodo T1 - Applications of SAR interferometric coherence time series BT - satiotemporal dynamics of geomorphic transitions in the South-Central Andes JF - Journal of geophysical research : Earth surface N2 - 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. KW - Copernicus KW - SAR KW - critical infrastructure resilience KW - early warning KW - landslides Y1 - 2020 U6 - https://doi.org/10.1029/2019JF005141 SN - 2169-9003 SN - 2169-9011 VL - 125 IS - 3 PB - American Geophysical Union CY - Washington ER - TY - GEN A1 - Smith, Taylor A1 - Bookhagen, Bodo T1 - Assessing Multi-Temporal Snow-Volume Trends in High Mountain Asia From 1987 to 2016 Using High-Resolution Passive Microwave Data T2 - Postprints der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe N2 - 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. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 1020 KW - snow KW - glacier KW - climate change KW - passive microwave KW - special sensor microwave imager KW - special sensor microwave imager/sounder Y1 - 2020 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-484176 SN - 1866-8372 IS - 1020 ER - TY - JOUR A1 - Smith, Taylor A1 - Bookhagen, Bodo T1 - Assessing Multi-Temporal Snow-Volume Trends in High Mountain Asia From 1987 to 2016 Using High-Resolution Passive Microwave Data JF - Frontiers in Earth Science N2 - 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. KW - snow KW - glacier KW - climate change KW - passive microwave KW - special sensor microwave imager KW - special sensor microwave imager/sounder Y1 - 2020 U6 - https://doi.org/10.3389/feart.2020.559175 SN - 2296-6463 VL - 8 PB - Frontiers Media CY - Lausanne ER - TY - GEN A1 - Milewski, Robert A1 - Chabrillat, Sabine A1 - Bookhagen, Bodo T1 - Analyses of Namibian Seasonal Salt Pan Crust Dynamics and Climatic Drivers Using Landsat 8 Time-Series and Ground Data T2 - Postprints der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe N2 - 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. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 988 KW - salt pan KW - playa KW - spectral analysis KW - crust KW - saline pan cycle KW - evaporites KW - time-series mapping Y1 - 2020 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-475685 SN - 1866-8372 IS - 988 ER - TY - JOUR A1 - Voss, Katalyn A. A1 - Bookhagen, Bodo A1 - Sachse, Dirk A1 - Chadwick, Oliver A. T1 - Variation of deuterium excess in surface waters across a 5000-m elevation gradient in eastern Nepal JF - Journal of hydrology N2 - 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. KW - stable isotopes KW - Himalaya KW - glacier KW - snow KW - precipitation KW - seasonality Y1 - 2020 U6 - https://doi.org/10.1016/j.jhydrol.2020.124802 SN - 0022-1694 SN - 1879-2707 VL - 586 PB - Elsevier CY - Amsterdam ER - TY - JOUR A1 - Castino, Fabiana A1 - Bookhagen, Bodo A1 - De la Torre, Alejandro T1 - Atmospheric dynamics of extreme discharge events from 1979 to 2016 in the southern Central Andes JF - Climate dynamics : observational, theoretical and computational research on the climate system N2 - 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. KW - South American monsoon system KW - cold surges KW - orographic barrier KW - mesoscale convective systems KW - extreme discharge KW - daily-discharge time KW - series ERA-interim Y1 - 2019 U6 - https://doi.org/10.1007/s00382-020-05458-1 SN - 0930-7575 SN - 1432-0894 VL - 55 IS - 11-12 SP - 3485 EP - 3505 PB - Springer CY - Berlin ; Heidelberg [u.a.] ER - TY - JOUR A1 - Milewski, Robert A1 - Chabrillat, Sabine A1 - Bookhagen, Bodo T1 - Analyses of Namibian Seasonal Salt Pan Crust Dynamics and Climatic Drivers Using Landsat 8 Time-Series and Ground Data JF - Remote Sensing N2 - 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. KW - salt pan KW - playa KW - spectral analysis KW - crust KW - saline pan cycle KW - evaporites KW - time-series mapping Y1 - 2019 U6 - https://doi.org/10.3390/rs12030474 SN - 2072-4292 IS - 3 PB - MDPI CY - Basel ER -