@article{GuenthelDonisKirillinetal.2019, author = {G{\"u}nthel, Marco and Donis, Daphne and Kirillin, Georgiy and Ionescu, Danny and Bizic, Mina and McGinnis, Daniel F. and Grossart, Hans-Peter and Tang, Kam W.}, title = {Contribution of oxic methane production to surface methane emission in lakes and its global importance}, series = {Nature Communications}, volume = {10}, journal = {Nature Communications}, publisher = {Nature Publishing Group UK}, address = {London}, issn = {2041-1723}, doi = {10.1038/s41467-019-13320-0}, pages = {10}, year = {2019}, abstract = {Recent discovery of oxic methane production in sea and lake waters, as well as wetlands, demands re-thinking of the global methane cycle and re-assessment of the contribution of oxic waters to atmospheric methane emission. Here we analysed system-wide sources and sinks of surface-water methane in a temperate lake. Using a mass balance analysis, we show that internal methane production in well-oxygenated surface water is an important source for surface-water methane during the stratified period. Combining our results and literature reports, oxic methane contribution to emission follows a predictive function of littoral sediment area and surface mixed layer volume. The contribution of oxic methane source(s) is predicted to increase with lake size, accounting for the majority (>50\%) of surface methane emission for lakes with surface areas >1 km(2).}, 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} } @article{HeimannVasyuraBathkeSudhausetal.2019, author = {Heimann, Sebastian and Vasyura-Bathke, Hannes and Sudhaus, Henriette and Isken, Marius Paul and Kriegerowski, Marius and Steinberg, Andreas and Dahm, Torsten}, title = {A Python framework for efficient use of pre-computed Green's functions in seismological and other physical forward and inverse source problems}, series = {Solid earth}, volume = {10}, journal = {Solid earth}, number = {6}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1869-9510}, doi = {10.5194/se-10-1921-2019}, pages = {1921 -- 1935}, year = {2019}, abstract = {The computation of such synthetic GFs is computationally and operationally demanding. As a consequence, the onthe-fly recalculation of synthetic GFs in each iteration of an optimisation is time-consuming and impractical. Therefore, the pre-calculation and efficient storage of synthetic GFs on a dense grid of source to receiver combinations enables the efficient lookup and utilisation of GFs in time-critical scenarios. We present a Python-based framework and toolkit - Pyrocko-GF - that enables the pre-calculation of synthetic GF stores, which are independent of their numerical calculation method and GF transfer function. The framework aids in the creation of such GF stores by interfacing a suite of established numerical forward modelling codes in seismology (computational back ends). So far, interfaces to back ends for layered Earth model cases have been provided; however, the architecture of Pyrocko-GF is designed to cover back ends for other geometries (e.g. full 3-D heterogeneous media) and other physical quantities (e.g. gravity, pressure, tilt). Therefore, Pyrocko-GF defines an extensible GF storage format suitable for a wide range of GF types, especially handling elasticity and wave propagation problems. The framework assists with visualisations, quality control, and the exchange of GF stores, which is supported through an online platform that provides many pre-calculated GF stores for local, regional, and global studies. The Pyrocko-GF toolkit comes with a well-documented application programming interface (API) for the Python programming language to efficiently facilitate forward modelling of geophysical processes, e.g. synthetic waveforms or static displacements for a wide range of source models.}, language = {en} } @article{AlHalbouniHolohanTaherietal.2019, author = {Al-Halbouni, Djamil and Holohan, Eoghan P. and Taheri, Abbas and Watson, Robert A. and Polom, Ulrich and Schoepfer, Martin P. J. and Emam, Sacha and Dahm, Torsten}, title = {Distinct element geomechanical modelling of the formation of sinkhole clusters within large-scale karstic depressions}, series = {Solid earth}, volume = {10}, journal = {Solid earth}, number = {4}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1869-9510}, doi = {10.5194/se-10-1219-2019}, pages = {1219 -- 1241}, year = {2019}, abstract = {The 2-D distinct element method (DEM) code (PFC2D_V5) is used here to simulate the evolution of subsidence-related karst landforms, such as single and clustered sinkholes, and associated larger-scale depressions. Subsurface material in the DEM model is removed progressively to produce an array of cavities; this simulates a network of subsurface groundwater conduits growing by chemical/mechanical erosion. The growth of the cavity array is coupled mechanically to the gravitationally loaded surroundings, such that cavities can grow also in part by material failure at their margins, which in the limit can produce individual collapse sinkholes. Two end-member growth scenarios of the cavity array and their impact on surface subsidence were examined in the models: (1) cavity growth at the same depth level and growth rate; (2) cavity growth at progressively deepening levels with varying growth rates. These growth scenarios are characterised by differing stress patterns across the cavity array and its overburden, which are in turn an important factor for the formation of sinkholes and uvalalike depressions. For growth scenario (1), a stable compression arch is established around the entire cavity array, hindering sinkhole collapse into individual cavities and favouring block-wise, relatively even subsidence across the whole cavity array. In contrast, for growth scenario (2), the stress system is more heterogeneous, such that local stress concentrations exist around individual cavities, leading to stress interactions and local wall/overburden fractures. Consequently, sinkhole collapses occur in individual cavities, which results in uneven, differential subsidence within a larger-scale depression. Depending on material properties of the cavity-hosting material and the overburden, the larger-scale depression forms either by sinkhole coalescence or by widespread subsidence linked geometrically to the entire cavity array. The results from models with growth scenario (2) are in close agreement with surface morphological and subsurface geophysical observations from an evaporite karst area on the eastern shore of the Dead Sea.}, language = {en} } @article{DeschlerNeherSchmidtMende2019, author = {Deschler, Felix and Neher, Dieter and Schmidt-Mende, Lukas}, title = {Perovskite semiconductors for next generation optoelectronic applications}, series = {APL Materials}, volume = {7}, journal = {APL Materials}, number = {8}, publisher = {American Institute of Physics}, address = {Melville}, issn = {2166-532X}, doi = {10.1063/1.5119744}, pages = {3}, year = {2019}, language = {en} } @article{NikolisMischokSiegmundetal.2019, author = {Nikolis, Vasileios C. and Mischok, Andreas and Siegmund, Bernhard and Kublitski, Jonas and Jia, Xiangkun and Benduhn, Johannes and H{\"o}rmann, Ulrich and Neher, Dieter and Gather, Malte C. and Spoltore, Donato and Vandewal, Koen}, title = {Strong light-matter coupling for reduced photon energy losses in organic photovoltaics}, series = {Nature Communications}, volume = {10}, journal = {Nature Communications}, publisher = {Nature Publ. Group}, address = {London}, issn = {2041-1723}, doi = {10.1038/s41467-019-11717-5}, pages = {8}, year = {2019}, abstract = {Strong light-matter coupling can re-arrange the exciton energies in organic semiconductors. Here, we exploit strong coupling by embedding a fullerene-free organic solar cell (OSC) photo-active layer into an optical microcavity, leading to the formation of polariton peaks and a red-shift of the optical gap. At the same time, the open-circuit voltage of the device remains unaffected. This leads to reduced photon energy losses for the low-energy polaritons and a steepening of the absorption edge. While strong coupling reduces the optical gap, the energy of the charge-transfer state is not affected for large driving force donor-acceptor systems. Interestingly, this implies that strong coupling can be exploited in OSCs to reduce the driving force for electron transfer, without chemical or microstructural modifications of the photoactive layer. Our work demonstrates that the processes determining voltage losses in OSCs can now be tuned, and reduced to unprecedented values, simply by manipulating the device architecture.}, language = {en} } @article{BiskabornSmithNoetzlietal.2019, author = {Biskaborn, Boris and Smith, Sharon L. and Noetzli, Jeannette and Matthes, Heidrun and Vieira, Goncalo and Streletskiy, Dmitry A. and Schoeneich, Philippe and Romanovsky, Vladimir E. and Lewkowicz, Antoni G. and Abramov, Andrey and Allard, Michel and Boike, Julia and Cable, William L. and Christiansen, Hanne H. and Delaloye, Reynald and Diekmann, Bernhard and Drozdov, Dmitry and Etzelmueller, Bernd and Grosse, Guido and Guglielmin, Mauro and Ingeman-Nielsen, Thomas and Isaksen, Ketil and Ishikawa, Mamoru and Johansson, Margareta and Johannsson, Halldor and Joo, Anseok and Kaverin, Dmitry and Kholodov, Alexander and Konstantinov, Pavel and Kroeger, Tim and Lambiel, Christophe and Lanckman, Jean-Pierre and Luo, Dongliang and Malkova, Galina and Meiklejohn, Ian and Moskalenko, Natalia and Oliva, Marc and Phillips, Marcia and Ramos, Miguel and Sannel, A. Britta K. and Sergeev, Dmitrii and Seybold, Cathy and Skryabin, Pavel and Vasiliev, Alexander and Wu, Qingbai and Yoshikawa, Kenji and Zheleznyak, Mikhail and Lantuit, Hugues}, title = {Permafrost is warming at a global scale}, series = {Nature Communications}, volume = {10}, journal = {Nature Communications}, publisher = {Nature Publ. Group}, address = {London}, issn = {2041-1723}, doi = {10.1038/s41467-018-08240-4}, pages = {11}, year = {2019}, abstract = {Permafrost warming has the potential to amplify global climate change, because when frozen sediments thaw it unlocks soil organic carbon. Yet to date, no globally consistent assessment of permafrost temperature change has been compiled. Here we use a global data set of permafrost temperature time series from the Global Terrestrial Network for Permafrost to evaluate temperature change across permafrost regions for the period since the International Polar Year (2007-2009). During the reference decade between 2007 and 2016, ground temperature near the depth of zero annual amplitude in the continuous permafrost zone increased by 0.39 +/- 0.15 degrees C. Over the same period, discontinuous permafrost warmed by 0.20 +/- 0.10 degrees C. Permafrost in mountains warmed by 0.19 +/- 0.05 degrees C and in Antarctica by 0.37 +/- 0.10 degrees C. Globally, permafrost temperature increased by 0.29 +/- 0.12 degrees C. The observed trend follows the Arctic amplification of air temperature increase in the Northern Hemisphere. In the discontinuous zone, however, ground warming occurred due to increased snow thickness while air temperature remained statistically unchanged.}, language = {en} } @article{BiskabornNazarovaPestryakovaetal.2019, author = {Biskaborn, Boris and Nazarova, Larisa B. and Pestryakova, Luidmila Agafyevna and Syrykh, Liudmila and Funck, Kim and Meyer, Hanno and Chapligin, Bernhard and Vyse, Stuart Andrew and Gorodnichev, Ruslan and Zakharov, Evgenii and Wang, Rong and Schwamborn, Georg and Bailey, Hannah L. and Diekmann, Bernhard}, title = {Spatial distribution of environmental indicators in surface sediments of Lake Bolshoe Toko, Yakutia, Russia}, series = {Biogeosciences}, volume = {16}, journal = {Biogeosciences}, number = {20}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1726-4170}, doi = {10.5194/bg-16-4023-2019}, pages = {4023 -- 4049}, year = {2019}, abstract = {Rapidly changing climate in the Northern Hemisphere and associated socio-economic impacts require reliable understanding of lake systems as important freshwater resources and sensitive sentinels of environmental change. To better understand time-series data in lake sediment cores, it is necessary to gain information on within-lake spatial variabilities of environmental indicator data. Therefore, we retrieved a set of 38 samples from the sediment surface along spatial habitat gradients in the boreal, deep, and yet pristine Lake Bolshoe Toko in southern Yakutia, Russia. Our methods comprise laboratory analyses of the sediments for multiple proxy parameters, including diatom and chironomid taxonomy, oxygen isotopes from diatom silica, grain-size distributions, elemental compositions (XRF), organic carbon content, and mineralogy (XRD). We analysed the lake water for cations, anions, and isotopes. Our results show that the diatom assemblages are strongly influenced by water depth and dominated by planktonic species, i.e. Pliocaenicus bolshetokoensis. Species richness and diversity are higher in the northern part of the lake basin, associated with the availability of benthic, i.e. periphytic, niches in shallower waters. delta O-18(diatom) values are higher in the deeper south-western part of the lake, probably related to water temperature differences. The highest amount of the chironomid taxa underrepresented in the training set used for palaeoclimate inference was found close to the Utuk River and at southern littoral and profundal sites. Abiotic sediment components are not symmetrically distributed in the lake basin, but vary along restricted areas of differential environmental forcing. Grain size and organic matter are mainly controlled by both river input and water depth. Mineral (XRD) data distributions are influenced by the methamorphic lithology of the Stanovoy mountain range, while elements (XRF) are intermingled due to catchment and diagenetic differences. We conclude that the lake represents a valuable archive for multiproxy environmental reconstruction based on diatoms (including oxygen isotopes), chironomids, and sediment-geochemical parameters. Our analyses suggest multiple coring locations preferably at intermediate depth in the northern basin and the deep part in the central basin, to account for representative bioindicator distributions and higher temporal resolution, respectively.}, language = {en} }