@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} }