@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{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{JobeLiBookhagenetal.2018, author = {Jobe, Jessica Ann Thompson and Li, Tao and Bookhagen, Bodo and Chen, Jie and Burbank, Douglas W.}, title = {Dating growth strata and basin fill by combining 26Al/10Be burial dating and magnetostratigraphy}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {1044}, issn = {1866-8372}, doi = {10.25932/publishup-46806}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-468067}, pages = {806 -- 828}, year = {2018}, abstract = {Cosmogenic burial dating enables dating of coarse-grained, Pliocene-Pleistocene sedimentary units that are typically difficult to date with traditional methods, such as magnetostratigraphy. In the actively deforming western Tarim Basin in NW China, Pliocene-Pleistocene conglomerates were dated at eight sites, integrating Al-26/Be-10 burial dating with previously published magnetostratigraphic sections. These samples were collected from growth strata on the flanks of growing folds and from sedimentary units beneath active faults to place timing constraints on the initiation of deformation of structures within the basin and on shortening rates on active faults. These new basin-fill and growthstrata ages document the late Neogene and Quaternary growth of the Pamir and Tian Shan orogens between >5 and 1 Ma and delineate the eastward propagation of deformation at rates up to 115 km/m.y. and basinward growth of both mountain belts at rates up to 12 km/m.y.}, language = {en} } @misc{TostEhmelHeidmannetal.2018, author = {Tost, Jordi and Ehmel, Fabian and Heidmann, Frank and Olen, Stephanie M. and Bookhagen, Bodo}, title = {Hazards and accessibility}, series = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, number = {710}, issn = {1866-8372}, doi = {10.25932/publishup-42785}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-427853}, pages = {8}, year = {2018}, abstract = {The assessment of natural hazards and risk has traditionally been built upon the estimation of threat maps, which are used to depict potential danger posed by a particular hazard throughout a given area. But when a hazard event strikes, infrastructure is a significant factor that can determine if the situation becomes a disaster. The vulnerability of the population in a region does not only depend on the area's local threat, but also on the geographical accessibility of the area. This makes threat maps by themselves insufficient for supporting real-time decision-making, especially for those tasks that involve the use of the road network, such as management of relief operations, aid distribution, or planning of evacuation routes, among others. To overcome this problem, this paper proposes a multidisciplinary approach divided in two parts. First, data fusion of satellite-based threat data and open infrastructure data from OpenStreetMap, introducing a threat-based routing service. Second, the visualization of this data through cartographic generalization and schematization. This emphasizes critical areas along roads in a simple way and allows users to visually evaluate the impact natural hazards may have on infrastructure. We develop and illustrate this methodology with a case study of landslide threat for an area in Colombia.}, language = {en} } @article{PurintonBookhagen2018, author = {Purinton, Benjamin and Bookhagen, Bodo}, title = {Measuring decadal vertical land-level changes from SRTM-C (2000) and TanDEM-X (∼ 2015) in the south-central Andes}, series = {Earth Surface Dynamics}, volume = {6}, journal = {Earth Surface Dynamics}, publisher = {Copernicus Publ.}, address = {G{\"o}ttingen}, issn = {2196-6311}, doi = {10.5194/esurf-6-971-2018}, pages = {971 -- 987}, year = {2018}, abstract = {In the arctic and high mountains it is common to measure vertical changes of ice sheets and glaciers via digital elevation model (DEM) differencing. This requires the signal of change to outweigh the noise associated with the datasets. Excluding large landslides, on the ice-free earth the land-level change is smaller in vertical magnitude and thus requires more accurate DEMs for differencing and identification of change. Previously, this has required meter to submeter data at small spatial scales. Following careful corrections, we are able to measure land-level changes in gravel-bed channels and steep hillslopes in the south-central Andes using the SRTM-C (collected in 2000) and the TanDEM-X (collected from 2010 to 2015) near-global 12-30m DEMs. Long-standing errors in the SRTM-C are corrected using the TanDEM-X as a control surface and applying cosine-fit co-registration to remove ∼ 1∕10 pixel (∼ 3m) shifts, fast Fourier transform (FFT) and filtering to remove SRTM-C short- and long-wavelength stripes, and blocked shifting to remove remaining complex biases. The datasets are then differenced and outlier pixels are identified as a potential signal for the case of gravel-bed channels and hillslopes. We are able to identify signals of incision and aggradation (with magnitudes down to ∼ 3m in the best case) in two  > 100km river reaches, with increased geomorphic activity downstream of knickpoints. Anthropogenic gravel excavation and piling is prominently measured, with magnitudes exceeding ±5m (up to  > 10m for large piles). These values correspond to conservative average rates of 0.2 to > 0.5myr-1 for vertical changes in gravel-bed rivers. For hillslopes, since we require stricter cutoffs for noise, we are only able to identify one major landslide in the study area with a deposit volume of 16±0.15×106m3. Additional signals of change can be garnered from TanDEM-X auxiliary layers; however, these are more difficult to quantify. The methods presented can be extended to any region of the world with SRTM-C and TanDEM-X coverage where vertical land-level changes are of interest, with the caveat that remaining vertical uncertainties in primarily the SRTM-C limit detection in steep and complex topography.}, language = {en} } @article{OlenBookhagen2018, author = {Olen, Stephanie M. and Bookhagen, Bodo}, title = {Mapping Damage-Affected Areas after Natural Hazard Events Using Sentinel-1 Coherence Time Series}, series = {remote sensing}, volume = {10}, journal = {remote sensing}, number = {8}, publisher = {Molecular Diversity Preservation International (MDPI)}, address = {Basel}, issn = {2072-4292}, doi = {10.3390/rs10081272}, pages = {1 -- 19}, year = {2018}, abstract = {The emergence of the Sentinel-1A and 1B satellites now offers freely available and widely accessible Synthetic Aperture Radar (SAR) data. Near-global coverage and rapid repeat time (6-12 days) gives Sentinel-1 data the potential to be widely used for monitoring the Earth's surface. Subtle land-cover and land surface changes can affect the phase and amplitude of the C-band SAR signal, and thus the coherence between two images collected before and after such changes. Analysis of SAR coherence therefore serves as a rapidly deployable and powerful tool to track both seasonal changes and rapid surface disturbances following natural disasters. An advantage of using Sentinel-1 C-band radar data is the ability to easily construct time series of coherence for a region of interest at low cost. In this paper, we propose a new method for Potentially Affected Area (PAA) detection following a natural hazard event. Based on the coherence time series, the proposed method (1) determines the natural variability of coherence within each pixel in the region of interest, accounting for factors such as seasonality and the inherent noise of variable surfaces; and (2) compares pixel-by-pixel syn-event coherence to temporal coherence distributions to determine where statistically significant coherence loss has occurred. The user can determine to what degree the syn-event coherence value (e.g., 1st, 5th percentile of pre-event distribution) constitutes a PAA, and integrate pertinent regional data, such as population density, to rank and prioritise PAAs. We apply the method to two case studies, Sarpol-e, Iran following the 2017 Iran-Iraq earthquake, and a landslide-prone region of NW Argentina, to demonstrate how rapid identification and interpretation of potentially affected areas can be performed shortly following a natural hazard event.}, language = {en} } @misc{PurintonBookhagen2018, author = {Purinton, Benjamin and Bookhagen, Bodo}, title = {Measuring decadal vertical land-level changes from SRTM-C (2000) and TanDEM-X (∼ 2015) in the south-central Andes}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {480}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-420487}, pages = {16}, year = {2018}, abstract = {In the arctic and high mountains it is common to measure vertical changes of ice sheets and glaciers via digital elevation model (DEM) differencing. This requires the signal of change to outweigh the noise associated with the datasets. Excluding large landslides, on the ice-free earth the land-level change is smaller in vertical magnitude and thus requires more accurate DEMs for differencing and identification of change. Previously, this has required meter to submeter data at small spatial scales. Following careful corrections, we are able to measure land-level changes in gravel-bed channels and steep hillslopes in the south-central Andes using the SRTM-C (collected in 2000) and the TanDEM-X (collected from 2010 to 2015) near-global 12-30m DEMs. Long-standing errors in the SRTM-C are corrected using the TanDEM-X as a control surface and applying cosine-fit co-registration to remove ∼ 1∕10 pixel (∼ 3m) shifts, fast Fourier transform (FFT) and filtering to remove SRTM-C short- and long-wavelength stripes, and blocked shifting to remove remaining complex biases. The datasets are then differenced and outlier pixels are identified as a potential signal for the case of gravel-bed channels and hillslopes. We are able to identify signals of incision and aggradation (with magnitudes down to ∼ 3m in the best case) in two  > 100km river reaches, with increased geomorphic activity downstream of knickpoints. Anthropogenic gravel excavation and piling is prominently measured, with magnitudes exceeding ±5m (up to  > 10m for large piles). These values correspond to conservative average rates of 0.2 to > 0.5myr-1 for vertical changes in gravel-bed rivers. For hillslopes, since we require stricter cutoffs for noise, we are only able to identify one major landslide in the study area with a deposit volume of 16±0.15×106m3. Additional signals of change can be garnered from TanDEM-X auxiliary layers; however, these are more difficult to quantify. The methods presented can be extended to any region of the world with SRTM-C and TanDEM-X coverage where vertical land-level changes are of interest, with the caveat that remaining vertical uncertainties in primarily the SRTM-C limit detection in steep and complex topography.}, language = {en} }