@article{HoffmannFeakinsBookhagenetal.2016, author = {Hoffmann, Bernd and Feakins, Sarah J. and Bookhagen, Bodo and Olen, Stephanie M. and Adhikari, Danda P. and Mainali, Janardan and Sachse, Dirk}, title = {Climatic and geomorphic drivers of plant organic matter transport in the Arun River, E Nepal}, series = {Earth \& planetary science letters}, volume = {452}, journal = {Earth \& planetary science letters}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0012-821X}, doi = {10.1016/j.epsl.2016.07.008}, pages = {104 -- 114}, year = {2016}, language = {en} } @misc{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{OlenBookhagen2018, author = {Olen, Stephanie M. and Bookhagen, Bodo}, title = {Mapping Damage-Affected Areas after Natural Hazard Events Using Sentinel-1 Coherence Time Series}, series = {remote sensing}, volume = {10}, journal = {remote sensing}, number = {8}, publisher = {Molecular Diversity Preservation International (MDPI)}, address = {Basel}, issn = {2072-4292}, doi = {10.3390/rs10081272}, pages = {1 -- 19}, year = {2018}, abstract = {The emergence of the Sentinel-1A and 1B satellites now offers freely available and widely accessible Synthetic Aperture Radar (SAR) data. Near-global coverage and rapid repeat time (6-12 days) gives Sentinel-1 data the potential to be widely used for monitoring the Earth's surface. Subtle land-cover and land surface changes can affect the phase and amplitude of the C-band SAR signal, and thus the coherence between two images collected before and after such changes. Analysis of SAR coherence therefore serves as a rapidly deployable and powerful tool to track both seasonal changes and rapid surface disturbances following natural disasters. An advantage of using Sentinel-1 C-band radar data is the ability to easily construct time series of coherence for a region of interest at low cost. In this paper, we propose a new method for Potentially Affected Area (PAA) detection following a natural hazard event. Based on the coherence time series, the proposed method (1) determines the natural variability of coherence within each pixel in the region of interest, accounting for factors such as seasonality and the inherent noise of variable surfaces; and (2) compares pixel-by-pixel syn-event coherence to temporal coherence distributions to determine where statistically significant coherence loss has occurred. The user can determine to what degree the syn-event coherence value (e.g., 1st, 5th percentile of pre-event distribution) constitutes a PAA, and integrate pertinent regional data, such as population density, to rank and prioritise PAAs. We apply the method to two case studies, Sarpol-e, Iran following the 2017 Iran-Iraq earthquake, and a landslide-prone region of NW Argentina, to demonstrate how rapid identification and interpretation of potentially affected areas can be performed shortly following a natural hazard event.}, language = {en} } @article{OlenBookhagen2020, author = {Olen, Stephanie M. and Bookhagen, Bodo}, title = {Applications of SAR interferometric coherence time series}, series = {Journal of geophysical research : Earth surface}, volume = {125}, journal = {Journal of geophysical research : Earth surface}, number = {3}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9003}, doi = {10.1029/2019JF005141}, pages = {22}, year = {2020}, abstract = {Sediment transport domains in mountain landscapes are characterized by fundamentally different processes and rates depending on several factors, including geology, climate, and biota. Accurately identifying where transitions between transport domains occur is an important step to quantify the past, present, and future contribution of varying erosion and sedimentation processes and enhance our predictive capabilities. We propose a new methodology based on time series of synthetic aperture radar (SAR) interferometric coherence images to map sediment transport regimes across arid and semiarid landscapes. Using 4 years of Sentinel-1 data, we analyze sediment transport regimes for the south-central Andes in northwestern Argentina characterized by steep topographic and climatic gradients. We observe seasonally low coherence during the regional wet season, particularly on hillslopes and in alluvial channels. The spatial distribution of coherence is compared to drainage areas extracted from digital topography to identify two distinct transitions within watersheds: (a) a hillslope-to-fluvial and (b) a fluvial-to-alluvial transition. While transitions within a given basin can be well-constrained, the relative role of each sediment transport domain varies widely over the climatic and topographic gradients. In semiarid regions, we observe larger relative contributions from hillslopes compared to arid regions. Across regional gradients, the range of coherence within basins positively correlates to previously published millennial catchment-wide erosion rates and to topographic metrics used to indicate long-term uplift. Our study suggests that a dense time series of interferometric coherence can be used as a proxy for surface sediment movement and landscape stability in vegetation-free settings at event to decadal timescales.}, language = {en} } @misc{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}, journal = {remote sensing}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-417766}, pages = {19}, year = {2018}, abstract = {The emergence of the Sentinel-1A and 1B satellites now offers freely available and widely accessible Synthetic Aperture Radar (SAR) data. Near-global coverage and rapid repeat time (6-12 days) gives Sentinel-1 data the potential to be widely used for monitoring the Earth's surface. Subtle land-cover and land surface changes can affect the phase and amplitude of the C-band SAR signal, and thus the coherence between two images collected before and after such changes. Analysis of SAR coherence therefore serves as a rapidly deployable and powerful tool to track both seasonal changes and rapid surface disturbances following natural disasters. An advantage of using Sentinel-1 C-band radar data is the ability to easily construct time series of coherence for a region of interest at low cost. In this paper, we propose a new method for Potentially Affected Area (PAA) detection following a natural hazard event. Based on the coherence time series, the proposed method (1) determines the natural variability of coherence within each pixel in the region of interest, accounting for factors such as seasonality and the inherent noise of variable surfaces; and (2) compares pixel-by-pixel syn-event coherence to temporal coherence distributions to determine where statistically significant coherence loss has occurred. The user can determine to what degree the syn-event coherence value (e.g., 1st, 5th percentile of pre-event distribution) constitutes a PAA, and integrate pertinent regional data, such as population density, to rank and prioritise PAAs. We apply the method to two case studies, Sarpol-e, Iran following the 2017 Iran-Iraq earthquake, and a landslide-prone region of NW Argentina, to demonstrate how rapid identification and interpretation of potentially affected areas can be performed shortly following a natural hazard event.}, language = {en} } @article{MeeseBookhagenOlenetal.2018, author = {Meese, Bernd and Bookhagen, Bodo and Olen, Stephanie M. and Barthold, Frauke Katrin and Sachse, Dirk}, title = {The effect of Indian Summer Monsoon rainfall on surface water delta D values in the central Himalaya}, series = {Hydrological processes}, volume = {32}, journal = {Hydrological processes}, number = {24}, publisher = {Wiley}, address = {Hoboken}, issn = {0885-6087}, doi = {10.1002/hyp.13281}, pages = {3662 -- 3674}, year = {2018}, abstract = {Stable isotope proxy records, such as speleothems, plant-wax biomarker records, and ice cores, are suitable archives for the reconstruction of regional palaeohydrologic conditions. But the interpretation of these records in the tropics, especially in the Indian Summer Monsoon (ISM) domain, is difficult due to differing moisture and water sources: precipitation from the ISM and Winter Westerlies, as well as snow- and glacial meltwater. In this study, we use interannual differences in ISM strength (2011-2012) to understand the stable isotopic composition of surface water in the Arun River catchment in eastern Nepal. We sampled main stem and tributary water (n = 204) for stable hydrogen and oxygen isotope analysis in the postmonsoon phase of two subsequent years with significantly distinct ISM intensities. In addition to the 2011/2012 sampling campaigns, we collected a 12-month time series of main stem waters (2012/2013, n = 105) in order to better quantify seasonal effects on the variability of surface water delta O-18/delta D. Furthermore, remotely sensed satellite data of rainfall, snow cover, glacial coverage, and evapotranspiration was evaluated. The comparison of datasets from both years revealed that surface waters of the main stem Arun and its tributaries were D-enriched by similar to 15 parts per thousand when ISM rainfall decreased by 20\%. This strong response emphasizes the importance of the ISM for surface water run-off in the central Himalaya. However, further spatio-temporal analysis of remote sensing data in combination with stream water d-excess revealed that most high-altitude tributaries and the Tibetan part of the Arun receive high portions of glacial melt water and likely Winter Westerly Disturbances precipitation. We make the following two implications: First, palaeohydrologic archives found in high-altitude tributaries and on the southern Tibetan Plateau record a mixture of past precipitation delta D values and variable amounts of additional water sources. Second, surface water isotope ratios of lower elevated tributaries strongly reflect the isotopic composition of ISM rainfall implying a suitable region for the analysis of potential delta D value proxy records.}, language = {en} } @phdthesis{Olen2016, author = {Olen, Stephanie M.}, title = {Understanding Himalayan denudation at the catchment and orogen scale}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-91423}, school = {Universit{\"a}t Potsdam}, pages = {xx, 174}, year = {2016}, abstract = {Understanding the rates and processes of denudation is key to unraveling the dynamic processes that shape active orogens. This includes decoding the roles of tectonic and climate-driven processes in the long-term evolution of high- mountain landscapes in regions with pronounced tectonic activity and steep climatic and surface-process gradients. Well-constrained denudation rates can be used to address a wide range of geologic problems. In steady-state landscapes, denudation rates are argued to be proportional to tectonic or isostatic uplift rates and provide valuable insight into the tectonic regimes underlying surface denudation. The use of denudation rates based on terrestrial cosmogenic nuclide (TCN) such as 10Beryllium has become a widely-used method to quantify catchment-mean denudation rates. Because such measurements are averaged over timescales of 102 to 105 years, they are not as susceptible to stochastic changes as shorter-term denudation rate estimates (e.g., from suspended sediment measurements) and are therefore considered more reliable for a comparison to long-term processes that operate on geologic timescales. However, the impact of various climatic, biotic, and surface processes on 10Be concentrations and the resultant denudation rates remains unclear and is subject to ongoing discussion. In this thesis, I explore the interaction of climate, the biosphere, topography, and geology in forcing and modulating denudation rates on catchment to orogen scales. There are many processes in highly dynamic active orogens that may effect 10Be concentrations in modern river sands and therefore impact 10Be-derived denudation rates. The calculation of denudation rates from 10Be concentrations, however, requires a suite of simplifying assumptions that may not be valid or applicable in many orogens. I investigate how these processes affect 10Be concentrations in the Arun Valley of Eastern Nepal using 34 new 10Be measurements from the main stem Arun River and its tributaries. The Arun Valley is characterized by steep gradients in climate and topography, with elevations ranging from <100 m asl in the foreland basin to >8,000 asl in the high sectors to the north. This is coupled with a five-fold increase in mean annual rainfall across strike of the orogen. Denudation rates from tributary samples increase toward the core of the orogen, from <0.2 to >5 mm/yr from the Lesser to Higher Himalaya. Very high denudation rates (>2 mm/yr), however, are likely the result of 10Be TCN dilution by surface and climatic processes, such as large landsliding and glaciation, and thus may not be representative of long-term denudation rates. Mainstem Arun denudation rates increase downstream from ~0.2 mm/yr at the border with Tibet to 0.91 mm/yr at its outlet into the Sapt Kosi. However, the downstream 10Be concentrations may not be representative of the entire upstream catchment. Instead, I document evidence for downstream fining of grains from the Tibetan Plateau, resulting in an order-of-magnitude apparent decrease in the measured 10Be concentration. In the Arun Valley and across the Himalaya, topography, climate, and vegetation are strongly interrelated. The observed increase in denudation rates at the transition from the Lesser to Higher Himalaya corresponds to abrupt increases in elevation, hillslope gradient, and mean annual rainfall. Thus, across strike (N-S), it is difficult to decipher the potential impacts of climate and vegetation cover on denudation rates. To further evaluate these relationships I instead took advantage of an along-strike west-to-east increase of mean annual rainfall and vegetation density in the Himalaya. An analysis of 136 published 10Be denudation rates from along strike of the revealed that median denudation rates do not vary considerably along strike of the Himalaya, ~1500 km E-W. However, the range of denudation rates generally decreases from west to east, with more variable denudation rates in the northwestern regions of the orogen than in the eastern regions. This denudation rate variability decreases as vegetation density increases (R=- 0.90), and increases proportionately to the annual seasonality of vegetation (R=0.99). Moreover, rainfall and vegetation modulate the relationship between topographic steepness and denudation rates such that in the wet, densely vegetated regions of the Himalaya, topography responds more linearly to changes in denudation rates than in dry, sparsely vegetated regions, where the response of topographic steepness to denudation rates is highly nonlinear. Understanding the relationships between denudation rates, topography, and climate is also critical for interpreting sedimentary archives. However, there is a lack of understanding of how terrestrial organic matter is transported out of orogens and into sedimentary archives. Plant wax lipid biomarkers derived from terrestrial and marine sedimentary records are commonly used as paleo- hydrologic proxy to help elucidate these problems. I address the issue of how to interpret the biomarker record by using the plant wax isotopic composition of modern suspended and riverbank organic matter to identify and quantify organic matter source regions in the Arun Valley. Topographic and geomorphic analysis, provided by the 10Be catchment-mean denudation rates, reveals that a combination of topographic steepness (as a proxy for denudation) and vegetation density is required to capture organic matter sourcing in the Arun River. My studies highlight the importance of a rigorous and careful interpretation of denudation rates in tectonically active orogens that are furthermore characterized by strong climatic and biotic gradients. Unambiguous information about these issues is critical for correctly decoding and interpreting the possible tectonic and climatic forces that drive erosion and denudation, and the manifestation of the erosion products in sedimentary archives.}, language = {en} } @article{OlenBookhagenHoffmannetal.2015, author = {Olen, Stephanie M. and Bookhagen, Bodo and Hoffmann, Bernd and Sachse, Dirk and Adhikari, Danda P. and Strecker, Manfred}, title = {Understanding erosion rates in the Himalayan orogen: A case study from the Arun Valley}, series = {Journal of geophysical research : Earth surface}, volume = {120}, journal = {Journal of geophysical research : Earth surface}, number = {10}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9003}, doi = {10.1002/2014JF003410}, pages = {2080 -- 2102}, year = {2015}, abstract = {Understanding the rates and pattern of erosion is a key aspect of deciphering the impacts of climate and tectonics on landscape evolution. Denudation rates derived from terrestrial cosmogenic nuclides (TCNs) are commonly used to quantify erosion and bridge tectonic (Myr) and climatic (up to several kiloyears) time scales. However, how the processes of erosion in active orogens are ultimately reflected in Be-10 TCN samples remains a topic of discussion. We investigate this problem in the Arun Valley of eastern Nepal with 34 new Be-10-derived catchment-mean denudation rates. The Arun Valley is characterized by steep north-south gradients in topography and climate. Locally, denudation rates increase northward, from <0.2mmyr(-1) to similar to 1.5mmyr(-1) in tributary samples, while main stem samples appear to increase downstream from similar to 0.2mmyr(-1) at the border with Tibet to 0.91mmyr(-1) in the foreland. Denudation rates most strongly correlate with normalized channel steepness (R-2=0.67), which has been commonly interpreted to indicate tectonic activity. Significant downstream decrease of Be-10 concentration in the main stem Arun suggests that upstream sediment grains are fining to the point that they are operationally excluded from the processed sample. This results in Be-10 concentrations and denudation rates that do not uniformly represent the upstream catchment area. We observe strong impacts on Be-10 concentrations from local, nonfluvial geomorphic processes, such as glaciation and landsliding coinciding with areas of peak rainfall rates, pointing toward climatic modulation of predominantly tectonically driven denudation rates.}, language = {en} } @article{OlenBookhagenStrecker2020, author = {Olen, Stephanie M. and Bookhagen, Bodo and Strecker, Manfred}, title = {Corrigendum to: Olen, Stephanie M.; Bookhagen, Bodo; Strecker, Manfred R. : Role of climate and vegetation density in modulating denudation rates in the Himalaya. - Earth and planetary science letters. - 445 (2016), S. 57 - 67. - doi: https://doi.org/10.1016/j.epsl.2016.03.047}, series = {Earth and planetary science letters}, volume = {540}, journal = {Earth and planetary science letters}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0012-821X}, doi = {10.1016/j.epsl.2020.116252}, pages = {5}, year = {2020}, abstract = {Vegetation has long been hypothesized to influence the nature and rates of surface processes. We test the possible impact of vegetation and climate on denudation rates at orogen scale by taking advantage of a pronounced along-strike gradient in rainfall and vegetation density in the Himalaya. We combine 12 new 10Be denudation rates from the Sutlej Valley and 123 published denudation rates from fluvially- dominated catchments in the Himalaya with remotely-sensed measures of vegetation density and rainfall metrics, and with tectonic and lithologic constraints. In addition, we perform topographic analyses to assess the contribution of vegetation and climate in modulating denudation rates along strike. We observe variations in denudation rates and the relationship between denudation and topography along strike that are most strongly controlled by local rainfall amount and vegetation density, and cannot be explained by along-strike differences in tectonics or lithology. A W-E along-strike decrease in denudation rate variability positively correlates with the seasonality of vegetation density (R = 0.95, p < 0.05), and negatively correlates with mean vegetation density (R = -0.84, p < 0.05). Vegetation density modulates the topographic response to changing denudation rates, such that the functional relationship between denudation rate and topographic steepness becomes increasingly linear as vegetation density increases. We suggest that while tectonic processes locally control the pattern of denudation rates across strike of the Himalaya (i.e., S-N), along strike of the orogen (i.e., E-W) climate exerts a measurable influence on how denudation rates scatter around long-term, tectonically-controlled erosion, and on the functional relationship between topography and denudation}, language = {en} } @article{OlenBookhagenStrecker2016, author = {Olen, Stephanie M. and Bookhagen, Bodo and Strecker, Manfred}, title = {Role of climate and vegetation density in modulating denudation rates in the Himalaya}, series = {Earth \& planetary science letters}, volume = {445}, journal = {Earth \& planetary science letters}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0012-821X}, doi = {10.1016/j.epsl.2016.03.047}, pages = {57 -- 67}, year = {2016}, abstract = {Vegetation has long been hypothesized to influence the nature and rates of surface processes. We test the possible impact of vegetation and climate on denudation rates at orogen scale by taking advantage of a pronounced along-strike gradient in rainfall and vegetation density in the Himalaya. We combine 12 new Be-10 denudation rates from the Sutlej Valley and 123 published denudation rates from fluvially-dominated catchments in the Himalaya with remotely-sensed measures of vegetation density and rainfall metrics, and with tectonic and lithologic constraints. In addition, we perform topographic analyses to assess the contribution of vegetation and climate in modulating denudation rates along strike. We observe variations in denudation rates and the relationship between denudation and topography along strike that are most strongly controlled by local rainfall amount and vegetation density, and cannot be explained by along-strike differences in tectonics or lithology. A W-E along-strike decrease in denudation rate variability positively correlates with the seasonality of vegetation density (R = 0.95, p < 0.05), and negatively correlates with mean vegetation density (R = -0.84, p < 0.05). Vegetation density modulates the topographic response to changing denudation rates, such that the functional relationship between denudation rate and topographic steepness becomes increasingly linear as vegetation density increases. We suggest that while tectonic processes locally control the pattern of denudation rates across strike of the Himalaya (i.e., S-N), along strike of the orogen (i.e., E-W) climate exerts a measurable influence on how denudation rates scatter around long-term, tectonically-controlled erosion, and on the functional relationship between topography and denudation. (C) 2016 Elsevier B.V. All rights reserved.}, language = {en} } @article{OlenEhlersDensmore2012, author = {Olen, Stephanie M. and Ehlers, Todd and Densmore, Mathew S.}, title = {Limits to reconstructing paleotopography from thermochronometer data}, series = {JOURNAL OF GEOPHYSICAL RESEARCH-EARTH SURFACE}, volume = {117}, journal = {JOURNAL OF GEOPHYSICAL RESEARCH-EARTH SURFACE}, publisher = {AMER GEOPHYSICAL UNION}, address = {WASHINGTON}, issn = {0148-0227}, doi = {10.1029/2011JF001985}, pages = {20}, year = {2012}, abstract = {Recent studies suggest that orogens can achieve a topographic steady state whereby equilibrium is reached between tectonics and erosion. However, steady state topography may not be the norm in many orogens experiencing large changes in climate or tectonics, which can produce topographic transients. The quantification of transient topography over geologic timescales requires reconstructing paleotopography, but this has proven difficult in many cases. This study investigates the utility of bedrock thermochronometer data to reconstruct orogen paleotopography over million year timescales. Apatite (U-Th)/He and fission track ages are integrated with a thermokinematic model for a single-parameter inversion of paleotopography. An iterative scheme is used that minimizes the misfit between predicted and observed cooling ages to identify the range of paleotopographies that could produce observed ages within sample uncertainty. Two approaches are considered. First, synthetic 2- D topographies are used to test the robustness of the approach. The following topographic evolution scenarios are considered: (1) lateral ridge migration, (2) topographic relief change, and (3) valley widening and deepening from glaciation. Second, the method is applied in three dimensions to existing data from the Coast Mountains of British Columbia, Canada. Results from both applications of the model suggest that (1) paleotopographic reconstruction will typically underpredict the magnitude of topographic change, especially relief change; (2) paleotopography is most successfully reconstructed after lateral ridge migration in long-wavelength topographies; and (3) reconstructed paleotopography from the Coast Mountains, British Columbia, suggests that glacial erosion may have the potential to remove drainage divides and laterally shift topographic ridges and peaks.}, language = {en} }