@article{vonSpechtOeztuerkVehetal.2019, author = {von Specht, Sebastian and {\"O}zt{\"u}rk, Ugur and Veh, Georg and Cotton, Fabrice and Korup, Oliver}, title = {Effects of finite source rupture on landslide triggering}, series = {Solid earth}, volume = {10}, journal = {Solid earth}, number = {2}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1869-9510}, doi = {10.5194/se-10-463-2019}, pages = {463 -- 486}, year = {2019}, abstract = {The propagation of a seismic rupture on a fault introduces spatial variations in the seismic wave field surrounding the fault. This directivity effect results in larger shaking amplitudes in the rupture propagation direction. Its seismic radiation pattern also causes amplitude variations between the strike-normal and strike-parallel components of horizontal ground motion. We investigated the landslide response to these effects during the 2016 Kumamoto earthquake (M-w 7.1) in central Kyushu (Japan). Although the distribution of some 1500 earthquake-triggered landslides as a function of rupture distance is consistent with the observed Arias intensity, the landslides were more concentrated to the northeast of the southwest-northeast striking rupture. We examined several landslide susceptibility factors: hillslope inclination, the median amplification factor (MAF) of ground shaking, lithology, land cover, and topographic wetness. None of these factors sufficiently explains the landslide distribution or orientation (aspect), although the landslide head scarps have an elevated hillslope inclination and MAF. We propose a new physics-based ground-motion model (GMM) that accounts for the seismic rupture effects, and we demonstrate that the low-frequency seismic radiation pattern is consistent with the overall landslide distribution. Its spatial pattern is influenced by the rupture directivity effect, whereas landslide aspect is influenced by amplitude variations between the fault-normal and fault-parallel motion at frequencies < 2 Hz. This azimuth dependence implies that comparable landslide concentrations can occur at different distances from the rupture. This quantitative link between the prevalent landslide aspect and the low-frequency seismic radiation pattern can improve coseismic landslide hazard assessment.}, language = {en} } @article{NakoudiGiannakakiDandouetal.2019, author = {Nakoudi, Konstantina and Giannakaki, Elina and Dandou, Aggeliki and Tombrou, Maria and Komppula, Mika}, title = {Planetary boundary layer height by means of lidar and numerical simulations over New Delhi, India}, series = {Atmospheric measurement techniques : an interactive open access journal of the European Geosciences Union}, volume = {12}, journal = {Atmospheric measurement techniques : an interactive open access journal of the European Geosciences Union}, number = {5}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1867-1381}, doi = {10.5194/amt-12-2595-2019}, pages = {2595 -- 2610}, year = {2019}, abstract = {In this work, the height of the planetary boundary layer (PBLH) is investigated over Gwal Pahari (Gual Pahari), New Delhi, for almost a year. To this end, ground-based measurements from a multiwavelength Raman lidar were used. The modified wavelet covariance transform (WCT) method was utilized for PBLH retrievals. Results were compared to data from Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) and the Weather Research and Forecasting (WRF) model. In order to examine the difficulties of PBLH detection from lidar, we analyzed three cases of PBLH diurnal evolution under different meteorological and aerosol load conditions. In the presence of multiple aerosol layers, the employed algorithm exhibited high efficiency (r = 0.9) in the attribution of PBLH, whereas weak aerosol gradients induced high variability in the PBLH. A sensitivity analysis corroborated the stability of the utilized methodology. The comparison with CALIPSO observations yielded satisfying results (r = 0.8), with CALIPSO slightly overestimating the PBLH. Due to the relatively warmer and drier winter and, correspondingly, colder and rainier pre-monsoon season, the seasonal PBLH cycle during the measurement period was slightly weaker than the cycle expected from long-term climate records.}, language = {en} } @article{PornsawadSapsakulBoeckmann2019, author = {Pornsawad, Pornsarp and Sapsakul, Nantawan and B{\"o}ckmann, Christine}, title = {A modified asymptotical regularization of nonlinear ill-posed problems}, series = {Mathematics}, volume = {7}, journal = {Mathematics}, edition = {5}, publisher = {MDPI}, address = {Basel, Schweiz}, issn = {2227-7390}, doi = {10.3390/math7050419}, pages = {19}, year = {2019}, abstract = {In this paper, we investigate the continuous version of modified iterative Runge-Kutta-type methods for nonlinear inverse ill-posed problems proposed in a previous work. The convergence analysis is proved under the tangential cone condition, a modified discrepancy principle, i.e., the stopping time T is a solution of ∥𝐹(𝑥𝛿(𝑇))-𝑦𝛿∥=𝜏𝛿+ for some 𝛿+>𝛿, and an appropriate source condition. We yield the optimal rate of convergence.}, language = {en} } @article{KoebschWinkelLiebneretal.2019, author = {Koebsch, Franziska and Winkel, Matthias and Liebner, Susanne and Liu, Bo and Westphal, Julia and Schmiedinger, Iris and Spitzy, Alejandro and Gehre, Matthias and Jurasinski, Gerald and K{\"o}hler, Stefan and Unger, Viktoria and Koch, Marian and Sachs, Torsten and B{\"o}ttcher, Michael E.}, title = {Sulfate deprivation triggers high methane production in a disturbed and rewetted coastal peatland}, series = {Biogeosciences}, volume = {16}, journal = {Biogeosciences}, number = {9}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1726-4170}, doi = {10.5194/bg-16-1937-2019}, pages = {1937 -- 1953}, year = {2019}, abstract = {In natural coastal wetlands, high supplies of marine sulfate suppress methanogenesis. Coastal wetlands are, however, often subject to disturbance by diking and drainage for agricultural use and can turn to potent methane sources when rewetted for remediation. This suggests that preceding land use measures can suspend the sulfate-related methane suppressing mechanisms. Here, we unravel the hydrological relocation and biogeochemical S and C transformation processes that induced high methane emissions in a disturbed and rewetted peatland despite former brackish impact. The underlying processes were investigated along a transect of increasing distance to the coastline using a combination of concentration patterns, stable isotope partitioning, and analysis of the microbial community structure. We found that diking and freshwater rewetting caused a distinct freshening and an efficient depletion of the brackish sulfate reservoir by dissimilatory sulfate reduction (DSR). Despite some legacy effects of brackish impact expressed as high amounts of sedimentary S and elevated electrical conductivities, contemporary metabolic processes operated mainly under sulfate-limited conditions. This opened up favorable conditions for the establishment of a prospering methanogenic community in the top 30-40 cm of peat, the structure and physiology of which resemble those of terrestrial organic-rich environments. Locally, high amounts of sulfate persisted in deeper peat layers through the inhibition of DSR, probably by competitive electron acceptors of terrestrial origin, for example Fe(III). However, as sulfate occurred only in peat layers below 30-40 cm, it did not interfere with high methane emissions on an ecosystem scale. Our results indicate that the climate effect of disturbed and remediated coastal wetlands cannot simply be derived by analogy with their natural counterparts. From a greenhouse gas perspective, the re-exposure of diked wetlands to natural coastal dynamics would literally open up the floodgates for a replenishment of the marine sulfate pool and therefore constitute an efficient measure to reduce methane emissions.}, language = {en} } @article{MitzscherlingHornWinterfeldetal.2019, author = {Mitzscherling, Julia and Horn, Fabian and Winterfeld, Maria and Mahler, Linda and Kallmeyer, Jens and Overduin, Pier Paul and Schirrmeister, Lutz and Winkel, Matthias and Grigoriev, Mikhail N. and Wagner, Dirk and Liebner, Susanne}, title = {Microbial community composition and abundance after millennia of submarine permafrost warming}, series = {Biogeosciences}, volume = {16}, journal = {Biogeosciences}, number = {19}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1726-4170}, doi = {10.5194/bg-16-3941-2019}, pages = {3941 -- 3958}, year = {2019}, abstract = {Warming of the Arctic led to an increase in permafrost temperatures by about 0.3 degrees C during the last decade. Permafrost warming is associated with increasing sediment water content, permeability, and diffusivity and could in the long term alter microbial community composition and abundance even before permafrost thaws. We studied the long-term effect (up to 2500 years) of submarine permafrost warming on microbial communities along an onshore-offshore transect on the Siberian Arctic Shelf displaying a natural temperature gradient of more than 10 degrees C. We analysed the in situ development of bacterial abundance and community composition through total cell counts (TCCs), quantitative PCR of bacterial gene abundance, and amplicon sequencing and correlated the microbial community data with temperature, pore water chemistry, and sediment physicochemical parameters. On timescales of centuries, permafrost warming coincided with an overall decreasing microbial abundance, whereas millennia after warming microbial abundance was similar to cold onshore permafrost. In addition, the dissolved organic carbon content of all cores was lowest in submarine permafrost after millennial-scale warming. Based on correlation analysis, TCC, unlike bacterial gene abundance, showed a significant rank-based negative correlation with increasing temperature, while bacterial gene copy numbers showed a strong negative correlation with salinity. Bacterial community composition correlated only weakly with temperature but strongly with the pore water stable isotopes delta O-18 and delta D, as well as with depth. The bacterial community showed substantial spatial variation and an overall dominance of Actinobacteria, Chloroflexi, Firmicutes, Gemmatimonadetes, and Proteobacteria, which are amongst the microbial taxa that were also found to be active in other frozen permafrost environments. We suggest that, millennia after permafrost warming by over 10 degrees C, microbial community composition and abundance show some indications for proliferation but mainly reflect the sedimentation history and paleoenvironment and not a direct effect through warming.}, language = {en} } @article{HaberPohlmeierToetzkeLehmannetal.2019, author = {Haber-Pohlmeier, Sabina and T{\"o}tzke, Christian and Lehmann, E. and Kardjilov, Nikolay and Pohlmeier, A. and Oswald, Sascha}, title = {Combination of magnetic resonance imaging and neutron computed tomography for three-dimensional rhizosphere imaging}, series = {Vadose zone journal}, volume = {18}, journal = {Vadose zone journal}, number = {1}, publisher = {Soil Science Society of America}, address = {Madison}, issn = {1539-1663}, doi = {10.2136/vzj2018.09.0166}, pages = {11}, year = {2019}, abstract = {Core Ideas 3D MRI relaxation time maps reflect water mobility in root, rhizosphere, and soil. 3D NCT water content maps of the same plant complement relaxation time maps. The relaxation time T1 decreases from soil to root, whereas water content increases. Parameters together indicate modification of rhizosphere pore space by gel phase. The zone of reduced T1 corresponds to the zone remaining dry after rewetting. In situ investigations of the rhizosphere require high-resolution imaging techniques, which allow a look into the optically opaque soil compartment. We present the novel combination of magnetic resonance imaging (MRI) and neutron computed tomography (NCT) to achieve synergistic information such as water mobility in terms of three-dimensional (3D) relaxation time maps and total water content maps. Besides a stationary MRI scanner for relaxation time mapping, we used a transportable MRI system on site in the NCT facility to capture rhizosphere properties before desiccation and after subsequent rewetting. First, we addressed two questions using water-filled test capillaries between 0.1 and 5 mm: which root diameters can still be detected by both methods, and to what extent are defined interfaces blurred by these imaging techniques? Going to real root system architecture, we demonstrated the sensitivity of the transportable MRI device by co-registration with NCT and additional validation using X-ray computed tomography. Under saturated conditions, we observed for the rhizosphere in situ a zone with shorter T1 relaxation time across a distance of about 1 mm that was not caused by reduced water content, as proven by successive NCT measurements. We conclude that the effective pore size in the pore network had changed, induced by a gel phase. After rewetting, NCT images showed a dry zone persisting while the MRI intensity inside the root increased considerably, indicating water uptake from the surrounding bulk soil through the still hydrophobic rhizosphere. Overall, combining NCT and MRI allows a more detailed analysis of the rhizosphere's functioning.}, language = {en} } @article{WeissQiuBarbotetal.2019, author = {Weiss, Jonathan R. and Qiu, Qiang and Barbot, Sylvain and Wright, Tim J. and Foster, James H. and Saunders, Alexander and Brooks, Benjamin A. and Bevis, Michael and Kendrick, Eric and Ericksen, Todd L. and Avery, Jonathan and Smalley, Robert and Cimbaro, Sergio R. and Lenzano, Luis Eduardo and Baron, Jorge and Carlos Baez, Juan and Echalar, Arturo}, title = {Illuminating subduction zone rheological properties in the wake of a giant earthquake}, series = {Science Advances}, volume = {5}, journal = {Science Advances}, number = {12}, publisher = {American Assoc. for the Advancement of Science}, address = {Washington}, issn = {2375-2548}, doi = {10.1126/sciadv.aax6720}, pages = {11}, year = {2019}, abstract = {Deformation associated with plate convergence at subduction zones is accommodated by a complex system involving fault slip and viscoelastic flow. These processes have proven difficult to disentangle. The 2010 M-w 8.8 Maule earthquake occurred close to the Chilean coast within a dense network of continuously recording Global Positioning System stations, which provide a comprehensive history of surface strain. We use these data to assemble a detailed picture of a structurally controlled megathrust fault frictional patchwork and the three-dimensional rheological and time-dependent viscosity structure of the lower crust and upper mantle, all of which control the relative importance of afterslip and viscoelastic relaxation during postseismic deformation. These results enhance our understanding of subduction dynamics including the interplay of localized and distributed deformation during the subduction zone earthquake cycle.}, language = {en} } @article{RaultRobertMarcetal.2019, author = {Rault, Claire and Robert, Alexandra and Marc, Odin and Hovius, Niels and Meunier, Patrick}, title = {Seismic and geologic controls on spatial clustering of landslides in three large earthquakes}, series = {Earth surface dynamics}, volume = {7}, journal = {Earth surface dynamics}, number = {3}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {2196-6311}, doi = {10.5194/esurf-7-829-2019}, pages = {829 -- 839}, year = {2019}, abstract = {The large, shallow earthquakes at Northridge, California (1994), Chi-Chi, Taiwan (1999), and Wenchuan, China (2008), each triggered thousands of landslides. We have determined the position of these landslides along hillslopes, normalizing for statistical bias. The landslide patterns have a co-seismic signature, with clustering at ridge crests and slope toes. A cross-check against rainfall-induced landslide inventories seems to confirm that crest clustering is specific to seismic triggering as observed in previous studies. In our three study areas, the seismic ground motion parameters and lithologic and topographic features used do not seem to exert a primary control on the observed patterns of landslide clustering. However, we show that at the scale of the epicentral area, crest and toe clustering occur in areas with specific geological features. Toe clustering of seismically induced landslides tends to occur along regional major faults. Crest clustering is concentrated at sites where the lithology along hillslopes is approximately uniform, or made of alternating soft and hard strata, and without strong overprint of geological structures. Although earthquake-induced landslides locate higher on hillslopes in a statistically significant way, geological features strongly modulate the landslide position along the hillslopes. As a result the observation of landslide clustering on topographic ridges cannot be used as a definite indicator of the topographic amplification of ground shaking.}, 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{JingHesseKumaretal.2019, author = {Jing, Miao and Hesse, Falk and Kumar, Rohini and Kolditz, Olaf and Kalbacher, Thomas and Attinger, Sabine}, title = {Influence of input and parameter uncertainty on the prediction of catchment-scale groundwater travel time distributions}, series = {Hydrology and earth system sciences : HESS}, volume = {23}, journal = {Hydrology and earth system sciences : HESS}, number = {1}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1027-5606}, doi = {10.5194/hess-23-171-2019}, pages = {171 -- 190}, year = {2019}, abstract = {Groundwater travel time distributions (TTDs) provide a robust description of the subsurface mixing behavior and hydrological response of a subsurface system. Lagrangian particle tracking is often used to derive the groundwater TTDs. The reliability of this approach is subjected to the uncertainty of external forcings, internal hydraulic properties, and the interplay between them. Here, we evaluate the uncertainty of catchment groundwater TTDs in an agricultural catchment using a 3-D groundwater model with an overall focus on revealing the relationship between external forcing, internal hydraulic properties, and TTD predictions. Eight recharge realizations are sampled from a high-resolution dataset of land surface fluxes and states. Calibration-constrained hydraulic conductivity fields (Ks fields) are stochastically generated using the null-space Monte Carlo (NSMC) method for each recharge realization. The random walk particle tracking (RWPT) method is used to track the pathways of particles and compute travel times. Moreover, an analytical model under the random sampling (RS) assumption is fit against the numerical solutions, serving as a reference for the mixing behavior of the model domain. The StorAge Selection (SAS) function is used to interpret the results in terms of quantifying the systematic preference for discharging young/old water. The simulation results reveal the primary effect of recharge on the predicted mean travel time (MTT). The different realizations of calibration-constrained Ks fields moderately magnify or attenuate the predicted MTTs. The analytical model does not properly replicate the numerical solution, and it underestimates the mean travel time. Simulated SAS functions indicate an overall preference for young water for all realizations. The spatial pattern of recharge controls the shape and breadth of simulated TTDs and SAS functions by changing the spatial distribution of particles' pathways. In conclusion, overlooking the spatial nonuniformity and uncertainty of input (forcing) will result in biased travel time predictions. We also highlight the worth of reliable observations in reducing predictive uncertainty and the good interpretability of SAS functions in terms of understanding catchment transport processes.}, language = {en} }