@article{BrunelloAndermannMarcetal.2020,
  author    = {Brunello, Camilla Francesca and Andermann, Christoff and Marc, Odin and Schneider, Katharina A. and Comiti, Francesco and Achleitner, Stefan and Hovius, Niels},
  title     = {Annually resolved monsoon onset and withdrawal dates across the Himalayas derived from local precipitation statistics},
  series = {Geophysical research letters},
  volume    = {47},
  journal   = {Geophysical research letters},
  number    = {23},
  publisher = {American Geophysical Union},
  address   = {Washington},
  issn      = {0094-8276},
  doi       = {10.1029/2020GL088420},
  pages     = {12},
  year      = {2020},
  abstract  = {A local and flexible definition of the monsoon season based on hydrological evidence is important for the understanding and management of Himalayan water resources. Here, we present an objective statistical method to retrieve seasonal hydrometeorological transitions. Applied to daily rainfall data (1951-2015), this method shows an average longitudinal delay of similar to 15 days, with later monsoon onset and earlier withdrawal in the western Himalaya, consistent with the continental progression of wet air masses. This delay leads to seasons of different length along the Himalaya and biased precipitation amounts when using uniform calendric monsoon boundaries. In the Central Himalaya annual precipitation has increased, due primarily to an increase of premonsoon precipitation. These findings highlight issues associated with a static definition of monsoon boundaries and call for a deeper understanding of nonmonsoonal precipitation over the Himalayan water tower. <br /> Plain Language Summary Precipitation in the Himalayas determines water availability for the Indian foreland with large socioeconomic implications. Despite its importance, spatial and temporal patterns of precipitation are poorly understood. Here, we estimate the long-term average and trends of seasonal precipitation at the scale of individual catchments draining the Himalayas. We apply a statistical method to detect the timing of hydrometeorological seasons from local precipitation measurements, focusing on monsoon onset and withdrawal. We identify longitudinal and latitudinal delays, resulting in seasons of different length along and across the Himalayas. These spatial patterns and the annual variability of the monsoon boundaries mean that oft-used, fixed calendric dates, for example, 1 June to 30 September, may be inadequate for retrieving monsoon rainfall totals. Moreover, we find that, despite its prominent contribution to annual rainfall totals, the Indian summer monsoon cannot explain the increase of the annual precipitation over the Central Himalayas. Instead, this appears to be mostly driven by changes in premonsoon and winter rainfall. So far, little attention has been paid to premonsoon precipitation, but governed by evaporative processes and surface water availability, it may be enhanced by irrigation and changed land use in the Gangetic foreland.},
  language  = {en}
}
@misc{EmbersonHoviusGalyetal.2016,
  author    = {Emberson, Robert and Hovius, Niels and Galy, Albert and Marc, Odin},
  title     = {Oxidation of sulfides and rapid weathering in recent landslides},
  series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe},
  journal   = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe},
  number    = {553},
  issn      = {1866-8372},
  doi       = {10.25932/publishup-41232},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-412326},
  pages     = {16},
  year      = {2016},
  abstract  = {Linking together the processes of rapid physical erosion and the resultant chemical dissolution of rock is a crucial step in building an overall deterministic understanding of weathering in mountain belts. Landslides, which are the most volumetrically important geomorphic process at these high rates of erosion, can generate extremely high rates of very localised weathering. To elucidate how this process works we have taken advantage of uniquely intense landsliding, resulting from Typhoon Morakot, in the T'aimali River and surrounds in southern Taiwan. Combining detailed analysis of landslide seepage chemistry with estimates of catchment-by-catchment landslide volumes, we demonstrate that in this setting the primary role of landslides is to introduce fresh, highly labile mineral phases into the surface weathering environment. There, rapid weathering is driven by the oxidation of pyrite and the resultant sulfuric-acid-driven dissolution of primarily carbonate rock. The total dissolved load correlates well with dissolved sulfate - the chief product of this style of weathering - in both landslides and streams draining the area (R-2 = 0.841 and 0.929 respectively; p < 0.001 in both cases), with solute chemistry in seepage from landslides and catchments affected by significant landsliding governed by the same weathering reactions. The predominance of coupled carbonate-sulfuric-acid-driven weathering is the key difference between these sites and previously studied landslides in New Zealand (Emberson et al., 2016), but in both settings increasing volumes of landslides drive greater overall solute concentrations in streams. Bedrock landslides, by excavating deep below saprolite-rock interfaces, create conditions for weathering in which all mineral phases in a lithology are initially unweathered within landslide deposits. As a result, the most labile phases dominate the weathering immediately after mobilisation and during a transient period of depletion. This mode of dissolution can strongly alter the overall output of solutes from catchments and their contribution to global chemical cycles if landslide-derived material is retained in catchments for extended periods after mass wasting.},
  language  = {en}
}
@article{EmbersonHoviusGalyetal.2016,
  author    = {Emberson, Robert and Hovius, Niels and Galy, Albert and Marc, Odin},
  title     = {Oxidation of sulfides and rapid weathering in recent landslides},
  series = {Earth surface dynamics},
  volume    = {4},
  journal   = {Earth surface dynamics},
  publisher = {Copernicus},
  address   = {G{\"o}ttingen},
  issn      = {2196-6311},
  doi       = {10.5194/esurf-4-727-2016},
  pages     = {727 -- 742},
  year      = {2016},
  abstract  = {Bedrock landslides, by excavating deep below saprolite-rock interfaces, create conditions for weathering in which all mineral phases in a lithology are initially unweathered within landslide deposits. As a result, the most labile phases dominate the weathering immediately after mobilisation and during a transient period of depletion. This mode of dissolution can strongly alter the overall output of solutes from catchments and their contribution to global chemical cycles if landslide-derived material is retained in catchments for extended periods after mass wasting.},
  language  = {en}
}
@article{EmbersonHoviusGalyetal.2016,
  author    = {Emberson, Robert and Hovius, Niels and Galy, Albert and Marc, Odin},
  title     = {Chemical weathering in active mountain belts controlled by stochastic bedrock landsliding},
  series = {Nature geoscience},
  volume    = {9},
  journal   = {Nature geoscience},
  publisher = {Nature Publ. Group},
  address   = {New York},
  issn      = {1752-0894},
  doi       = {10.1038/NGEO2600},
  pages     = {42 -- +},
  year      = {2016},
  abstract  = {A link between chemical weathering and physical erosion exists at the catchment scale over a wide range of erosion rates(1,2). However, in mountain environments, where erosion rates are highest, weathering may be kinetically limited(3-5) and therefore decoupled from erosion. In active mountain belts, erosion is driven by bedrock landsliding(6) at rates that depend strongly on the occurrence of extreme rainfall or seismicity(7). Although landslides affect only a small proportion of the landscape, bedrock landsliding can promote the collection and slow percolation of surface runoff in highly fragmented rock debris and create favourable conditions for weathering. Here we show from analysis of surface water chemistry in the Southern Alps of New Zealand that weathering in bedrock landslides controls the variability in solute load of these mountain rivers. We find that systematic patterns in surface water chemistry are strongly associated with landslide occurrence at scales from a single hillslope to an entire mountain belt, and that landslides boost weathering rates and river solute loads over decades. We conclude that landslides couple erosion and weathering in fast-eroding uplands and, thus, mountain weathering is a stochastic process that is sensitive to climatic and tectonic controls on mass wasting processes.},
  language  = {en}
}
@article{IllienSensSchoenfelderAndermannetal.2022,
  author    = {Illien, Luc and Sens-Sch{\"o}nfelder, Christoph and Andermann, Christoff and Marc, Odin and Cook, Kristen L. and Adhikari, Lok Bijaya and Hovius, Niels},
  title     = {Seismic velocity recovery in the subsurface},
  series = {Journal of geophysical research : Solid earth},
  volume    = {127},
  journal   = {Journal of geophysical research : Solid earth},
  number    = {2},
  publisher = {American Geophysical Union},
  address   = {Washington},
  issn      = {2169-9313},
  doi       = {10.1029/2021JB023402},
  pages     = {18},
  year      = {2022},
  abstract  = {Shallow earthquakes frequently disturb the hydrological and mechanical state of the subsurface, with consequences for hazard and water management. Transient post-seismic hydrological behavior has been widely reported, suggesting that the recovery of material properties (relaxation) following ground shaking may impact groundwater fluctuations. However, the monitoring of seismic velocity variations associated with earthquake damage and hydrological variations are often done assuming that both effects are independent. In a field site prone to highly variable hydrological conditions, we disentangle the different forcing of the relative seismic velocity variations delta v retrieved from a small dense seismic array in Nepal in the aftermath of the 2015 M-w 7.8 Gorkha earthquake. We successfully model transient damage effects by introducing a universal relaxation function that contains a unique maximum relaxation timescale for the main shock and the aftershocks, independent of the ground shaking levels. Next, we remove the modeled velocity from the raw data and test whether the corresponding residuals agree with a background hydrological behavior we inferred from a previously calibrated groundwater model. The fitting of the delta v data with this model is improved when we introduce transient hydrological properties in the phase immediately following the main shock. This transient behavior, interpreted as an enhanced permeability in the shallow subsurface, lasts for similar to 6 months and is shorter than the damage relaxation (similar to 1 yr). Thus, we demonstrate the capability of seismic interferometry to deconvolve transient hydrological properties after earthquakes from non-linear mechanical recovery.},
  language  = {en}
}
@article{KoberZeilingerHippeetal.2015,
  author    = {Kober, Florian and Zeilinger, Gerold and Hippe, Kristina and Marc, Odin and Lendzioch, Theodora and Grischott, Reto and Christl, Marcus and Kubik, Peter W. and Zola, Ramiro},
  title     = {Tectonic and lithological controls on denudation rates in the central Bolivian Andes},
  series = {Tectonophysics : international journal of geotectonics and the geology and physics of the interior of the earth},
  volume    = {657},
  journal   = {Tectonophysics : international journal of geotectonics and the geology and physics of the interior of the earth},
  publisher = {Elsevier},
  address   = {Amsterdam},
  issn      = {0040-1951},
  doi       = {10.1016/j.tecto.2015.06.037},
  pages     = {230 -- 244},
  year      = {2015},
  abstract  = {The topographic signature of a mountain belt depends on the interplay of tectonic, climatic and erosional processes, whose relative importance changes over times, while quantifying these processes and their rates at specific times remains a challenge. The eastern Andes of central Bolivia offer a natural laboratory in which such interplay has been debated. Here, we investigate the Rio Grande catchment which crosses orthogonally the eastern Andes orogen from the Eastern Cordillera into the Subandean Zone, exhibiting a catchment relief of up to 5000 m. Despite an enhanced tectonic activity in the Subandes, local relief, mean and modal slopes and channel steepness indices are largely similar compared to the Eastern Cordillera and the intervening Interandean Zone. Nevertheless, a dataset of 57 new cosmogenic 10Be and 26AI catchment wide denudation rates from the Rio Grande catchment reveals up to one order of magnitude higher denudation rates in the Subandean Zone (mean 0.8 mm/yr) compared to the upstream physiographic regions. We infer that tectonic activity in the thrusting dominated Subandean belt causes higher denudation rates based on cumulative rock uplift investigations and due to the absence of a pronounced climate gradient. Furthermore, the lower rock strength of the Subandean sedimentary units correlates with mean slopes similar to the ones of the Eastern Cordillera and Interandean Zone, highlighting the fact, that lithology and rock strength can control high denudation rates at low slopes. Low denudation rates measured at the outlet of the Rio Grande catchment (Abapo) are interpreted to be a result of a biased cosmogenic nuclide mixing that is dominated by headwater signals from the Eastern Cordillera and the Interandean zone and limited catchment sediment connectivity in the lower river reaches. Therefore, comparisons of short- (i.e., sediment yield) and millennial denudation rates require caution when postulating tectonic and/or climatic forcing without detailed studies. (C) 2015 The Authors. Published by Elsevier B.V.},
  language  = {en}
}
@phdthesis{Marc2016,
  author    = {Marc, Odin},
  title     = {Earthquake-induced landsliding},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-96808},
  school      = {Universit{\"a}t Potsdam},
  pages     = {xvi, 171},
  year      = {2016},
  abstract  = {Earthquakes deform Earth's surface, building long-lasting topographic features and contributing to landscape and mountain formation. However, seismic waves produced by earthquakes may also destabilize hillslopes, leading to large amounts of soil and bedrock moving downslope. Moreover, static deformation and shaking are suspected to damage the surface bedrock and therefore alter its future properties, affecting hydrological and erosional dynamics. Thus, earthquakes participate both in mountain building and stimulate directly or indirectly their erosion. Moreover, the impact of earthquakes on hillslopes has important implications for the amount of sediment and organic matter delivered to rivers, and ultimately to oceans, during episodic catastrophic seismic crises, the magnitude of life and property losses associated with landsliding, the perturbation and recovery of landscape properties after shaking, and the long term topographic evolution of mountain belts. Several of these aspects have been addressed recently through individual case studies but additional data compilation as well as theoretical or numerical modelling are required to tackle these issues in a more systematic and rigorous manner. This dissertation combines data compilation of earthquake characteristics, landslide mapping, and seismological data interpretation with physically-based modeling in order to address how earthquakes impact on erosional processes and landscape evolution. Over short time scales (10-100 s) and intermediate length scales (10 km), I have attempted to improve our understanding and ability to predict the amount of landslide debris triggered by seismic shaking in epicentral areas. Over long time scales (1-100 ky) and across a mountain belt (100 km) I have modeled the competition between erosional unloading and building of topography associated with earthquakes. Finally, over intermediate time scales (1-10 y) and at the hillslope scale (0.1-1 km) I have collected geomorphological and seismological data that highlight persistent effects of earthquakes on landscape properties and behaviour. First, I compiled a database on earthquakes that produced significant landsliding, including an estimate of the total landslide volume and area, and earthquake characteristics such as seismic moment and source depth. A key issue is the accurate conversion of landslide maps into volume estimates. Therefore I also estimated how amalgamation - when mapping errors lead to the bundling of multiple landslide into a single polygon - affects volume estimates from various earthquake-induced landslide inventories and developed an algorithm to automatically detect this artifact. The database was used to test a physically-based prediction of the total landslide area and volume caused by earthquakes, based on seismological scaling relationships and a statistical description of the landscape properties. The model outperforms empirical fits in accuracy, with 25 out of 40 cases well predicted, and allows interpretation of many outliers in physical terms. Apart from seismological complexities neglected by the model I found that exceptional rock strength properties or antecedent conditions may explain most outliers. Second, I assessed the geomorphic effects of large earthquakes on landscape dynamics by surveying the temporal evolution of precipitation-normalized landslide rate. I found strongly elevated landslide rates following earthquakes that progressively recover over 1 to 4 years, indicating that regolith strength drops and recovers. The relaxation is clearly non-linear for at least one case, and does not seem to correlate with coseismic landslide reactivation, water table level increase or tree root-system recovery. I suggested that shallow bedrock is damaged by the earthquake and then heals on annual timescales. Such variations in ground strength must be translated into shallow subsurface seismic velocities that are increasingly surveyed with ambient seismic noise correlations. With seismic noise autocorrelation I computed the seismic velocity in the epicentral areas of three earthquakes where I constrained a change in landslide rate. We found similar recovery dynamics and timescales, suggesting that seismic noise correlation techniques could be further developed to meaningfully assess ground strength variations for landscape dynamics. These two measurements are also in good agreement with the temporal dynamics of post-seismic surface displacement measured by GPS. This correlation suggests that the surface healing mechanism may be driven by tectonic deformation, and that the surface regolith and fractured bedrock may behave as a granular media that slowly compacts as it is sheared or vibrated. Last, I compared our model of earthquake-induced landsliding with a standard formulation of surface deformation caused by earthquakes to understand which parameters govern the competition between the building and destruction of topography caused by earthquakes. In contrast with previous studies I found that very large (Mw>8) earthquakes always increase the average topography, whereas only intermediate (Mw ~ 7) earthquakes in steep landscapes may reduce topography. Moreover, I illustrated how the net effect of earthquakes varies with depth or landscape steepness implying a complex and ambivalent role through the life of a mountain belt. Further I showed that faults producing a Gutenberg-Richter distribution of earthquake sizes, will limit topography over a larger range of fault sizes than faults producing repeated earthquakes with a characteristic size.},
  language  = {en}
}
@misc{MarcBehlingAndermannetal.2019,
  author    = {Marc, Odin and Behling, Robert and Andermann, Christoff and Turowski, Jens M. and Illien, Luc and Roessner, Sigrid and Hovius, Niels},
  title     = {Long-term erosion of the Nepal Himalayas by bedrock landsliding},
  series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe},
  journal   = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe},
  number    = {646},
  issn      = {1866-8372},
  doi       = {10.25932/publishup-42502},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-425022},
  pages     = {22},
  year      = {2019},
  abstract  = {In active mountain belts with steep terrain, bedrock landsliding is a major erosional agent. In the Himalayas, landsliding is driven by annual hydro-meteorological forcing due to the summer monsoon and by rarer, exceptional events, such as earthquakes. Independent methods yield erosion rate estimates that appear to increase with sampling time, suggesting that rare, high-magnitude erosion events dominate the erosional budget. Nevertheless, until now, neither the contribution of monsoon and earthquakes to landslide erosion nor the proportion of erosion due to rare, giant landslides have been quantified in the Himalayas. We address these challenges by combining and analysing earthquake- and monsoon-induced landslide inventories across different timescales. With time series of 5 m satellite images over four main valleys in central Nepal, we comprehensively mapped landslides caused by the monsoon from 2010 to 2018. We found no clear correlation between monsoon properties and landsliding and a similar mean landsliding rate for all valleys, except in 2015, where the valleys affected by the earthquake featured ∼ 5-8 times more landsliding than the pre-earthquake mean rate. The longterm size-frequency distribution of monsoon-induced landsliding (MIL) was derived from these inventories and from an inventory of landslides larger than ∼ 0.1 km 2 that occurred between 1972 and 2014. Using a published landslide inventory for the Gorkha 2015 earthquake, we derive the size-frequency distribution for earthquake-induced landsliding (EQIL). These two distributions are dominated by infrequent, large and giant landslides but under-predict an estimated Holocene frequency of giant landslides (> 1 km 3 ) which we derived from a literature compilation. This discrepancy can be resolved when modelling the effect of a full distribution of earthquakes of variable magnitude and when considering that a shallower earthquake may cause larger landslides. In this case, EQIL and MIL contribute about equally to a total long-term erosion of ∼ 2 ± 0.75 mm yr -1 in agreement with most thermo-chronological data. Independently of the specific total and relative erosion rates, the heavy-tailed size-frequency distribution from MIL and EQIL and the very large maximal landslide size in the Himalayas indicate that mean landslide erosion rates increase with sampling time, as has been observed for independent erosion estimates. Further, we find that the sampling timescale required to adequately capture the frequency of the largest landslides, which is necessary for deriving long-term mean erosion rates, is often much longer than the averaging time of cosmogenic 10 Be methods. This observation presents a strong caveat when interpreting spatial or temporal variability in erosion rates from this method. Thus, in areas where a very large, rare landslide contributes heavily to long-term erosion (as the Himalayas), we recommend 10 Be sample in catchments with source areas > 10 000 km 2 to reduce the method mean bias to below ∼ 20 \% of the long-term erosion.},
  language  = {en}
}
@article{MarcBehlingAndermannetal.2019,
  author    = {Marc, Odin and Behling, Robert and Andermann, Christoff and Turowski, Jens M. and Illien, Luc and Roessner, Sigrid and Hovius, Niels},
  title     = {Long-term erosion of the Nepal Himalayas by bedrock landsliding},
  series = {Earth surface dynamics},
  volume    = {7},
  journal   = {Earth surface dynamics},
  number    = {1},
  publisher = {Copernicus},
  address   = {G{\"o}ttingen},
  issn      = {2196-6311},
  doi       = {10.5194/esurf-7-107-2019},
  pages     = {107 -- 128},
  year      = {2019},
  abstract  = {In active mountain belts with steep terrain, bedrock landsliding is a major erosional agent. In the Himalayas, landsliding is driven by annual hydro-meteorological forcing due to the summer monsoon and by rarer, exceptional events, such as earthquakes. Independent methods yield erosion rate estimates that appear to increase with sampling time, suggesting that rare, high-magnitude erosion events dominate the erosional budget. Nevertheless, until now, neither the contribution of monsoon and earthquakes to landslide erosion nor the proportion of erosion due to rare, giant landslides have been quantified in the Himalayas. We address these challenges by combining and analysing earthquake- and monsoon-induced landslide inventories across different timescales. With time series of 5 m satellite images over four main valleys in central Nepal, we comprehensively mapped landslides caused by the monsoon from 2010 to 2018. We found no clear correlation between monsoon properties and landsliding and a similar mean landsliding rate for all valleys, except in 2015, where the valleys affected by the earthquake featured similar to 5-8 times more landsliding than the pre-earthquake mean rate. The longterm size-frequency distribution of monsoon-induced landsliding (MIL) was derived from these inventories and from an inventory of landslides larger than similar to 0.1 km(2) that occurred between 1972 and 2014. Using a published landslide inventory for the Gorkha 2015 earthquake, we derive the size-frequency distribution for earthquakeinduced landsliding (EQIL). These two distributions are dominated by infrequent, large and giant landslides but under-predict an estimated Holocene frequency of giant landslides (> 1 km(3)) which we derived from a literature compilation. This discrepancy can be resolved when modelling the effect of a full distribution of earthquakes of variable magnitude and when considering that a shallower earthquake may cause larger landslides. In this case, EQIL and MIL contribute about equally to a total long-term erosion of similar to 2 +/- 0.75 mm yr(-1) in agreement with most thermo-chronological data. Independently of the specific total and relative erosion rates, the heavy-tailed size-frequency distribution from MIL and EQIL and the very large maximal landslide size in the Himalayas indicate that mean landslide erosion rates increase with sampling time, as has been observed for independent erosion estimates. Further, we find that the sampling timescale required to adequately capture the frequency of the largest landslides, which is necessary for deriving long-term mean erosion rates, is often much longer than the averaging time of cosmogenic Be-10 methods. This observation presents a strong caveat when interpreting spatial or temporal variability in erosion rates from this method. Thus, in areas where a very large, rare landslide contributes heavily to long-term erosion (as the Himalayas), we recommend Be-10 sample in catchments with source areas > 10 000 km(2) to reduce the method mean bias to below similar to 20 \% of the long-term erosion.},
  language  = {en}
}
@article{MarcHovius2015,
  author    = {Marc, Odin and Hovius, Niels},
  title     = {Amalgamation in landslide maps},
  series = {Natural hazards and earth system sciences},
  volume    = {15},
  journal   = {Natural hazards and earth system sciences},
  number    = {4},
  publisher = {Copernicus},
  address   = {G{\"o}ttingen},
  issn      = {1561-8633},
  doi       = {10.5194/nhess-15-723-2015},
  pages     = {723 -- 733},
  year      = {2015},
  abstract  = {Inventories of individually delineated landslides are a key to understanding landslide physics and mitigating their impact. They permit assessment of area-frequency distributions and landslide volumes, and testing of statistical correlations between landslides and physical parameters such as topographic gradient or seismic strong motion. Amalgamation, i.e. the mapping of several adjacent landslides as a single polygon, can lead to potentially severe distortion of the statistics of these inventories. This problem can be especially severe in data sets produced by automated mapping. We present five inventories of earthquake-induced landslides mapped with different materials and techniques and affected by varying degrees of amalgamation. Errors on the total landslide volume and power-law exponent of the area-frequency distribution, resulting from amalgamation, may be up to 200 and 50 \%, respectively. We present an algorithm based on image and digital elevation model (DEM) analysis, for automatic identification of amalgamated polygons. On a set of about 2000 polygons larger than 1000 m(2), tracing landslides triggered by the 1994 Northridge earthquake, the algorithm performs well, with only 2.7-3.6\% incorrectly amalgamated landslides missed and 3.9-4.8\% correct polygons incorrectly identified as amalgams. This algorithm can be used broadly to check landslide inventories and allow faster correction by automating the identification of amalgamation.},
  language  = {en}
}