TY - JOUR A1 - Tofelde, Stefanie A1 - Schildgen, Taylor F. A1 - Savi, Sara A1 - Pingel, Heiko A1 - Wickert, Andrew D. A1 - Bookhagen, Bodo A1 - Wittmann, Hella A1 - Alonso, Ricardo N. A1 - Cottle, John A1 - Strecker, Manfred T1 - 100 kyr fluvial cut-and-fill terrace cycles since the Middle Pleistocene in the southern Central Andes, NW Argentina JF - Earth & planetary science letters N2 - Fluvial fill terraces in intermontane basins are valuable geomorphic archives that can record tectonically and/or climatically driven changes of the Earth-surface process system. However, often the preservation of fill terrace sequences is incomplete and/or they may form far away from their source areas, complicating the identification of causal links between forcing mechanisms and landscape response, especially over multi-millennial timescales. The intermontane Toro Basin in the southern Central Andes exhibits at least five generations of fluvial terraces that have been sculpted into several-hundred-meter-thick Quaternary valley-fill conglomerates. New surface-exposure dating using nine cosmogenic Be-10 depth profiles reveals the successive abandonment of these terraces with a 100 kyr cyclicity between 75 +/- 7 and 487 +/- 34 ka. Depositional ages of the conglomerates, determined by four Al-26/Be-10 burial samples and U-Pb zircon ages of three intercalated volcanic ash beds, range from 18 +/- 141 to 936 +/- 170 ka, indicating that there were multiple cut-and-fill episodes. Although the initial onset of aggradation at similar to 1 Ma and the overall net incision since ca. 500 ka can be linked to tectonic processes at the narrow basin outlet, the superimposed 100 kyr cycles of aggradation and incision are best explained by eccentricity-driven climate change. Within these cycles, the onset of river incision can be correlated with global cold periods and enhanced humid phases recorded in paleoclimate archives on the adjacent Bolivian Altiplano, whereas deposition occurred mainly during more arid phases on the Altiplano and global interglacial periods. We suggest that enhanced runoff during global cold phases - due to increased regional precipitation rates, reduced evapotranspiration, or both - resulted in an increased sediment-transport capacity in the Toro Basin, which outweighed any possible increases in upstream sediment supply and thus triggered incision. Compared with two nearby basins that record precessional (21-kyr) and long-eccentricity (400-kyr) forcing within sedimentary and geomorphic archives, the recorded cyclicity scales with the square of the drainage basin length. (C) 2017 Elsevier B.V. All rights reserved. KW - Be-10 depth-profiles KW - surface inflation KW - aggradation-incision cycles KW - glacial-interglacial cycles KW - landscape response to climate change KW - Eastern Cordillera Y1 - 2017 U6 - https://doi.org/10.1016/j.epsl.2017.06.001 SN - 0012-821X SN - 1385-013X VL - 473 SP - 141 EP - 153 PB - Elsevier CY - Amsterdam ER - TY - JOUR A1 - Brell, Maximilian A1 - Segl, Karl A1 - Guanter, Luis A1 - Bookhagen, Bodo T1 - 3D hyperspectral point cloud generation BT - Fusing airborne laser scanning and hyperspectral imaging sensors for improved object-based information extraction JF - ISPRS journal of photogrammetry and remote sensing : official publication of the International Society for Photogrammetry and Remote Sensing N2 - 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. KW - Lidar KW - Multispectral point cloud KW - Laser return intensity KW - Unmixing KW - Sharpening KW - Imaging spectroscopy KW - In-flight KW - Pixel level KW - Sensor fusion KW - Data fusion KW - Preprocessing KW - Point cloud segmentation KW - Semantic labeling Y1 - 2019 U6 - https://doi.org/10.1016/j.isprsjprs.2019.01.022 SN - 0924-2716 SN - 1872-8235 VL - 149 SP - 200 EP - 214 PB - Elsevier CY - Amsterdam 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 - JOUR A1 - Rheinwahlt, Aljoscha A1 - Goswami, Bedartha A1 - Bookhagen, Bodo T1 - A network-based flow accumulation algorithm for point clouds BT - Facet-Flow Networks (FFNs) JF - Journal of geophysical research : Earth surface N2 - Flow accumulation algorithms estimate the steady state of flow on real or modeled topographic surfaces and are crucial for hydrological and geomorphological assessments, including delineation of river networks, drainage basins, and sediment transport processes. Existing flow accumulation algorithms are typically designed to compute flows on regular grids and are not directly applicable to arbitrarily sampled topographic data such as lidar point clouds. In this study we present a random sampling scheme that generates homogeneous point densities, in combination with a novel flow path tracing approach-the Facet-Flow Network (FFN)-that estimates flow accumulation in terms of specific catchment area (SCA) on triangulated surfaces. The random sampling minimizes biases due to spatial sampling and the FFN allows for direct flow estimation from point clouds. We validate our approach on a Gaussian hill surface and study the convergence of its SCA compared to the analytical solution. Here, our algorithm outperforms the multiple flow direction algorithm, which is optimized for divergent surfaces. We also compute the SCA of a 6-km(2)-steep, vegetated catchment on Santa Cruz Island, California, based on airborne lidar point-cloud data. Point-cloud-based SCA values estimated by our method compare well with those estimated by the D-infinity or multiple flow direction algorithm on gridded data. The advantage of computing SCA from point clouds becomes relevant especially for divergent topography and for small drainage areas: These are depicted with much more detail due to the higher sampling density of point clouds. KW - point clouds KW - drainage networks KW - lidar KW - tin KW - surface runoff KW - spatial sampling Y1 - 2019 U6 - https://doi.org/10.1029/2018JF004827 SN - 2169-9003 SN - 2169-9011 VL - 124 IS - 7 SP - 2013 EP - 2033 PB - American Geophysical Union CY - Washington ER - TY - JOUR A1 - Bookhagen, Bodo A1 - Thiede, Rasmus Christoph A1 - Strecker, Manfred T1 - Abnormal monsoon years and their control on erosion and sediment flux in the high, and northwest Himalaya N2 - The interplay between topography and Indian summer monsoon circulation profoundly controls precipitation distribution, sediment transport, and river discharge along the Southern Himalayan Mountain Front (SHF). The Higher Himalayas form a major orographic barrier that separates humid sectors to the south and and regions to the north. During the Indian summer monsoon, vortices transport moisture from the Bay of Bengal, swirl along the SHF to the northwest, and cause heavy rainfall when colliding with the mountain front. In the eastern and central parts of the Himalaya, precipitation measurements derived from passive microwave analysis (SSM/I) show a strong gradient, with high values at medium elevations and extensive penetration of moisture along major river valleys into the orogen. The end of the monsoonal conveyer belt is near the Sutlej Valley in the NW Himalaya, where precipitation is lower and rainfall maxima move to lower elevations. This region thus comprises a climatic transition zone that is very sensitive to changes in Indian summer monsoon strength. To constrain magnitude, temporal, and spatial distribution of precipitation, we analyzed high-resolution passive microwave data from the last decade and identified an abnormal monsoon year (AMY) in 2002. During the 2002 AMY, violent rainstorms conquered orographic barriers and penetrated far into otherwise and regions in the northwest Himalaya at elevations in excess of 3 km asl. While precipitation in these regions was significantly increased and triggered extensive erosional processes (i.e., debris flows) on sparsely vegetated, steep hillslopes, mean rainfall along the low to medium elevations was not significantly greater in magnitude. This shift may thus play an important role in the overall sediment flux toward the Himalayan foreland. Using extended precipitation and sediment flux records for the last century, we show that these events have a decadal recurrence interval during the present-day monsoon circulation. Hence, episodically occurring AMYs control geomorphic processes primarily in the high-elevation and sectors of the orogen, while annual recurring monsoonal rainfall distribution dominates erosion in the low- to medium- elevation parts along the SHF. (C) 2004 Elsevier B.V. All rights reserved Y1 - 2005 ER - TY - GEN A1 - Brieger, Frederic A1 - Herzschuh, Ulrike A1 - Pestryakova, Luidmila Agafyevna A1 - Bookhagen, Bodo A1 - Zakharov, Evgenii S. A1 - Kruse, Stefan T1 - Advances in the derivation of Northeast Siberian forest metrics using high-resolution UAV-based photogrammetric point clouds T2 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe N2 - 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. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 1337 KW - UAV KW - photogrammetry KW - remote sensing KW - structure from motion KW - tundra–taiga ecotone KW - point cloud KW - forest structure Y1 - 2019 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-473318 SN - 1866-8372 IS - 1337 ER - TY - JOUR A1 - Brieger, Frederic A1 - Herzschuh, Ulrike A1 - Pestryakova, Luidmila Agafyevna A1 - Bookhagen, Bodo A1 - Zakharov, Evgenii S. A1 - Kruse, Stefan T1 - Advances in the Derivation of Northeast Siberian Forest Metrics Using High-Resolution UAV-Based Photogrammetric Point Clouds JF - Remote sensing N2 - 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. KW - UAV KW - photogrammetry KW - remote sensing KW - structure from motion KW - tundra-taiga ecotone KW - point cloud KW - forest structure Y1 - 2019 U6 - https://doi.org/10.3390/rs11121447 SN - 2072-4292 VL - 11 IS - 12 PB - MDPI CY - Basel ER - TY - JOUR A1 - Harvey, Jonathan E. A1 - Burbank, Douglas W. A1 - Bookhagen, Bodo T1 - Along-strike changes in Himalayan thrust geometry: Topographic and tectonic discontinuities in western Nepal JF - Lithosphere N2 - 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. Y1 - 2015 U6 - https://doi.org/10.1130/L444.1 SN - 1941-8264 SN - 1947-4253 VL - 7 IS - 5 SP - 511 EP - 518 PB - American Institute of Physics CY - Boulder ER - TY - JOUR A1 - Neely, Alexander B. A1 - Bookhagen, Bodo A1 - Burbank, Douglas W. T1 - An automated knickzone selection algorithm (KZ-Picker) to analyze transient landscapes: Calibration and validation JF - Journal of geophysical research : Earth surface N2 - Streams commonly respond to base-level fall by localizing erosion within steepened, convex knickzone reaches. Localized incision causes knickzone reaches to migrate upstream. Such migrating knickzones dictate the pace of landscape response to changes in tectonics or erosional efficiency and can help quantify the timing and source of base-level fall. Identification of knickzones typically requires individual selection of steepened reaches: a process that is tedious and subjective and has no efficient means to measure knickzone size. We construct an algorithm to automate this procedure by selecting the bounds of knickzone reaches in a -space (drainage-area normalized) framework. An automated feature calibrates algorithm parameters to a subset of knickzones handpicked by the user. The algorithm uses these parameters as consistent criteria to identify knickzones objectively, and then the algorithm measures the height, length, and slope of each knickzone reach. We test the algorithm on 1, 10, and 30m resolution digital elevation models (DEMs) of six catchments (trunk-stream lengths: 2.1-5.4km) on Santa Cruz Island, southern California. On the 1m DEM, algorithm-selected knickzones confirm 93% of handpicked knickzone positions (n=178) to a spatial accuracy of 100m, 88% to an accuracy within 50m, and 46% to an accuracy within 10m. Using 10 and 30m DEMs, accuracy is similar: 88-86% to 100m and 82% to 50m (n=38 and 36, respectively). The algorithm enables efficient regional comparison of the size and location of knickzones with geologic structures, mapped landforms, and hillslope morphology, thereby facilitating approaches to characterize the dynamics of transient landscapes. Plain Language Summary The shape of rivers reflects the environments that they flow through and the environments that they link together: mountains and oceans. Anywhere along the length of a river, changes in environmental conditions are propagated upstream and downstream as the river changes its morphology to match the new environmental conditions. Commonly, rivers steepen as land uplifts faster in regions of high tectonic convergence. The steepening of river gradients is propagated upstream and can be mapped to trace zones of high tectonic activity across landscapes and estimate the source and timing of environmental change. Such insights may indicate regions where earthquakes have become more frequent in the recent past and how rivers respond to these changes. In this submission, we detail an algorithm that can use digital topographic data (similar to google earth), to automatically map and measure anomalously steep river reaches across continental scales. This technology can highlight areas that have experienced recent sustained changes in environmental conditions, evident by changes in the morphology of rivers. Such environmental conditions could be changes in tectonic uplift and earthquake activity, changes in sea level, changes in land-use, or changes in climate, all factors that can produce measurable differences in river morphology over time. KW - knickpoint KW - transient KW - knickzone KW - incision KW - relict landscape KW - Santa Cruz Island Y1 - 2017 U6 - https://doi.org/10.1002/2017JF004250 SN - 2169-9003 SN - 2169-9011 VL - 122 SP - 1236 EP - 1261 PB - American Geophysical Union CY - Washington ER -