@article{AgarwalMarwanMaheswaranetal.2020, author = {Agarwal, Ankit and Marwan, Norbert and Maheswaran, Rathinasamy and {\"O}zt{\"u}rk, Ugur and Kurths, J{\"u}rgen and Merz, Bruno}, title = {Optimal design of hydrometric station networks based on complex network analysis}, series = {Hydrology and Earth System Sciences}, volume = {24}, journal = {Hydrology and Earth System Sciences}, number = {5}, publisher = {Copernicus Publ.}, address = {G{\"o}ttingen}, issn = {1027-5606}, doi = {10.5194/hess-24-2235-2020}, pages = {2235 -- 2251}, year = {2020}, abstract = {Hydrometric networks play a vital role in providing information for decision-making in water resource management. They should be set up optimally to provide as much information as possible that is as accurate as possible and, at the same time, be cost-effective. Although the design of hydrometric networks is a well-identified problem in hydrometeorology and has received considerable attention, there is still scope for further advancement. In this study, we use complex network analysis, defined as a collection of nodes interconnected by links, to propose a new measure that identifies critical nodes of station networks. The approach can support the design and redesign of hydrometric station networks. The science of complex networks is a relatively young field and has gained significant momentum over the last few years in different areas such as brain networks, social networks, technological networks, or climate networks. The identification of influential nodes in complex networks is an important field of research. We propose a new node-ranking measure - the weighted degree-betweenness (WDB) measure - to evaluate the importance of nodes in a network. It is compared to previously proposed measures used on synthetic sample networks and then applied to a real-world rain gauge network comprising 1229 stations across Germany to demonstrate its applicability. The proposed measure is evaluated using the decline rate of the network efficiency and the kriging error. The results suggest that WDB effectively quantifies the importance of rain gauges, although the benefits of the method need to be investigated in more detail.}, language = {en} } @article{AlbrechtWinkelmannLevermann2020, author = {Albrecht, Torsten and Winkelmann, Ricarda and Levermann, Anders}, title = {Glacial-cycle simulations of the Antarctic Ice Sheet with the Parallel Ice Sheet Model (PISM)}, series = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, volume = {14}, journal = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, number = {2}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1994-0416}, doi = {10.5194/tc-14-599-2020}, pages = {599 -- 632}, year = {2020}, abstract = {Simulations of the glacial-interglacial history of the Antarctic Ice Sheet provide insights into dynamic threshold behavior and estimates of the ice sheet's contributions to global sea-level changes for the past, present and future. However, boundary conditions are weakly constrained, in particular at the interface of the ice sheet and the bedrock. Also climatic forcing covering the last glacial cycles is uncertain, as it is based on sparse proxy data.
We use the Parallel Ice Sheet Model (PISM) to investigate the dynamic effects of different choices of input data, e.g., for modern basal heat flux or reconstructions of past changes of sea level and surface temperature. As computational resources are limited, glacial-cycle simulations are performed using a comparably coarse model grid of 16 km and various parameterizations, e.g., for basal sliding, iceberg calving, or for past variations in precipitation and ocean temperatures. In this study we evaluate the model's transient sensitivity to corresponding parameter choices and to different boundary conditions over the last two glacial cycles and provide estimates of involved uncertainties. We also discuss isolated and combined effects of climate and sea-level forcing. Hence, this study serves as a "cookbook" for the growing community of PISM users and paleo-ice sheet modelers in general.
For each of the different model uncertainties with regard to climatic forcing, ice and Earth dynamics, and basal processes, we select one representative model parameter that captures relevant uncertainties and motivates corresponding parameter ranges that bound the observed ice volume at present. The four selected parameters are systematically varied in a parameter ensemble analysis, which is described in a companion paper.}, language = {en} } @article{AllroggenBeiterTronicke2020, author = {Allroggen, Niklas and Beiter, Daniel and Tronicke, Jens}, title = {Ground-penetrating radar monitoring of fast subsurface processes}, series = {Geophysics}, volume = {85}, journal = {Geophysics}, number = {3}, publisher = {Society of Exploration Geophysicists}, address = {Tulsa}, issn = {0016-8033}, doi = {10.1190/GEO2019-0737.1}, pages = {A19 -- A23}, year = {2020}, abstract = {Earth and environmental sciences rely on detailed information about subsurface processes. Whereas geophysical techniques typically provide highly resolved spatial images, monitoring subsurface processes is often associated with enormous effort and, therefore, is usually limited to point information in time or space. Thus, the development of spatial and temporal continuous field monitoring methods is a major challenge for the understanding of subsurface processes. We have developed a novel method for ground-penetrating-radar (GPR) reflection monitoring of subsurface flow processes under unsaturated conditions and applied it to a hydrological infiltration experiment performed across a periglacial slope deposit in northwest Luxembourg. Our approach relies on a spatial and temporal quasicontinuous data recording and processing, followed by an attribute analysis based on analyzing differences between individual time steps. The results demonstrate the ability of time-lapse GPR monitoring to visualize the spatial and temporal dynamics of preferential flow processes with a spatial resolution in the order of a few decimeters and temporal resolution in the order of a few minutes. We observe excellent agreement with water table information originating from different boreholes. This demonstrates the potential of surface-based GPR reflection monitoring to observe the spatiotemporal dynamics of water movements in the subsurface. It provides valuable, and so far not accessible, information for example in the field of hydrology and pedology that allows studying the actual subsurface processes rather than deducing them from point information.}, language = {en} } @article{AngelopoulosOverduinWestermannetal.2020, author = {Angelopoulos, Michael and Overduin, Pier Paul and Westermann, Sebastian and Tronicke, Jens and Strauss, Jens and Schirrmeister, Lutz and Biskaborn, Boris and Liebner, Susanne and Maksimov, Georgii and Grigoriev, Mikhail N. and Grosse, Guido}, title = {Thermokarst lake to lagoon transitions in Eastern Siberia}, series = {Journal of geophysical research : Earth surface}, volume = {125}, journal = {Journal of geophysical research : Earth surface}, number = {10}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9003}, doi = {10.1029/2019JF005424}, pages = {21}, year = {2020}, abstract = {As the Arctic coast erodes, it drains thermokarst lakes, transforming them into lagoons, and, eventually, integrates them into subsea permafrost. Lagoons represent the first stage of a thermokarst lake transition to a marine setting and possibly more saline and colder upper boundary conditions. In this research, borehole data, electrical resistivity surveying, and modeling of heat and salt diffusion were carried out at Polar Fox Lagoon on the Bykovsky Peninsula, Siberia. Polar Fox Lagoon is a seasonally isolated water body connected to Tiksi Bay through a channel, leading to hypersaline waters under the ice cover. The boreholes in the center of the lagoon revealed floating ice and a saline cryotic bed underlain by a saline cryotic talik, a thin ice-bearing permafrost layer, and unfrozen ground. The bathymetry showed that most of the lagoon had bedfast ice in spring. In bedfast ice areas, the electrical resistivity profiles suggested that an unfrozen saline layer was underlain by a thick layer of refrozen talik. The modeling showed that thermokarst lake taliks can refreeze when submerged in saltwater with mean annual bottom water temperatures below or slightly above 0 degrees C. This occurs, because the top-down chemical degradation of newly formed ice-bearing permafrost is slower than the refreezing of the talik. Hence, lagoons may precondition taliks with a layer of ice-bearing permafrost before encroachment by the sea, and this frozen layer may act as a cap on gas migration out of the underlying talik.}, language = {en} } @article{ArnousZeckraVenerdinietal.2020, author = {Arnous, Ahmad and Zeckra, Martin and Venerdini, Agostina and Alvarado, Patricia and Arrowsmith, Ram{\´o}n and Guillemoteau, Julien and Landgraf, Angela and Guti{\´e}rrez, Adolfo Antonio and Strecker, Manfred}, title = {Neotectonic Activity in the Low-Strain Broken Foreland (Santa B{\´a}rbara System) of the North-Western Argentinean Andes (26°S)}, series = {Lithosphere}, volume = {2020}, journal = {Lithosphere}, number = {1}, publisher = {GSA}, address = {Boulder, Colo.}, issn = {1947-4253}, doi = {10.2113/2020/8888588}, pages = {1 -- 25}, year = {2020}, abstract = {Uplift in the broken Andean foreland of the Argentine Santa B{\´a}rbara System (SBS) is associated with the contractional reactivation of basement anisotropies, similar to those reported from the thick-skinned Cretaceous-Eocene Laramide province of North America. Fault scarps, deformed Quaternary deposits and landforms, disrupted drainage patterns, and medium-sized earthquakes within the SBS suggest that movement along these structures may be a recurring phenomenon, with yet to be defined repeat intervals and rupture lengths. In contrast to the Subandes thrust belt farther north, where eastward-migrating deformation has generated a well-defined thrust front, the SBS records spatiotemporally disparate deformation along structures that are only known to the first order. We present herein the results of geomorphic desktop analyses, structural field observations, and 2D electrical resistivity tomography and seismic-refraction tomography surveys and an interpretation of seismic reflection profiles across suspected fault scarps in the sedimentary basins adjacent to the Candelaria Range (CR) basement uplift, in the south-central part of the SBS. Our analysis in the CR piedmont areas reveals consistency between the results of near-surface electrical resistivity and seismic-refraction tomography surveys, the locations of prominent fault scarps, and structural geometries at greater depth imaged by seismic reflection data. We suggest that this deformation is driven by deep-seated blind thrusting beneath the CR and associated regional warping, while shortening involving Mesozoic and Cenozoic sedimentary strata in the adjacent basins was accommodated by layer-parallel folding and flexural-slip faults that cut through Quaternary landforms and deposits at the surface.}, language = {en} } @article{AyzelSchefferHeistermann2020, author = {Ayzel, Georgy and Scheffer, Tobias and Heistermann, Maik}, title = {RainNet v1.0}, series = {Geoscientific Model Development}, volume = {13}, journal = {Geoscientific Model Development}, number = {6}, publisher = {Copernicus Publ.}, address = {G{\"o}ttingen}, issn = {1991-959X}, doi = {10.5194/gmd-13-2631-2020}, pages = {2631 -- 2644}, year = {2020}, abstract = {In this study, we present RainNet, a deep convolutional neural network for radar-based precipitation nowcasting. Its design was inspired by the U-Net and SegNet families of deep learning models, which were originally designed for binary segmentation tasks. RainNet was trained to predict continuous precipitation intensities at a lead time of 5min, using several years of quality-controlled weather radar composites provided by the German Weather Service (DWD). That data set covers Germany with a spatial domain of 900km × 900km and has a resolution of 1km in space and 5min in time. Independent verification experiments were carried out on 11 summer precipitation events from 2016 to 2017. In order to achieve a lead time of 1h, a recursive approach was implemented by using RainNet predictions at 5min lead times as model inputs for longer lead times. In the verification experiments, trivial Eulerian persistence and a conventional model based on optical flow served as benchmarks. The latter is available in the rainymotion library and had previously been shown to outperform DWD's operational nowcasting model for the same set of verification events. RainNet significantly outperforms the benchmark models at all lead times up to 60min for the routine verification metrics mean absolute error (MAE) and the critical success index (CSI) at intensity thresholds of 0.125, 1, and 5mm h⁻¹. However, rainymotion turned out to be superior in predicting the exceedance of higher intensity thresholds (here 10 and 15mm h⁻¹). The limited ability of RainNet to predict heavy rainfall intensities is an undesirable property which we attribute to a high level of spatial smoothing introduced by the model. At a lead time of 5min, an analysis of power spectral density confirmed a significant loss of spectral power at length scales of 16km and below. Obviously, RainNet had learned an optimal level of smoothing to produce a nowcast at 5min lead time. In that sense, the loss of spectral power at small scales is informative, too, as it reflects the limits of predictability as a function of spatial scale. Beyond the lead time of 5min, however, the increasing level of smoothing is a mere artifact - an analogue to numerical diffusion - that is not a property of RainNet itself but of its recursive application. In the context of early warning, the smoothing is particularly unfavorable since pronounced features of intense precipitation tend to get lost over longer lead times. Hence, we propose several options to address this issue in prospective research, including an adjustment of the loss function for model training, model training for longer lead times, and the prediction of threshold exceedance in terms of a binary segmentation task. Furthermore, we suggest additional input data that could help to better identify situations with imminent precipitation dynamics. The model code, pretrained weights, and training data are provided in open repositories as an input for such future studies.}, language = {en} } @article{BaesSobolevGeryaetal.2020, author = {Baes, Marzieh and Sobolev, Stephan and Gerya, Taras V. and Brune, Sascha}, title = {Plume-induced subduction initiation}, series = {Geochemistry, geophysics, geosystems}, volume = {21}, journal = {Geochemistry, geophysics, geosystems}, number = {2}, publisher = {American Geophysical Union}, address = {Washington}, issn = {1525-2027}, doi = {10.1029/2019GC008663}, pages = {19}, year = {2020}, abstract = {Initiation of subduction following the impingement of a hot buoyant mantle plume is one of the few scenarios that allow breaking the lithosphere and recycling a stagnant lid without requiring any preexisting weak zones. Here, we investigate factors controlling the number and shape of retreating subducting slabs formed by plume-lithosphere interaction. Using 3-D thermomechanical models we show that the deformation regime, which defines formation of single-slab or multi-slab subduction, depends on several parameters such as age of oceanic lithosphere, thickness of the crust and large-scale lithospheric extension rate. Our model results indicate that on present-day Earth multi-slab plume-induced subduction is initiated only if the oceanic lithosphere is relatively young (<30-40 Myr, but >10 Myr), and the crust has a typical thickness of 8 km. In turn, development of single-slab subduction is facilitated by older lithosphere and pre-imposed extensional stresses. In early Earth, plume-lithosphere interaction could have led to formation of either episodic short-lived circular subduction when the oceanic lithosphere was young or to multi-slab subduction when the lithosphere was old.}, language = {en} } @article{BaesSobolevGeryaetal.2020, author = {Baes, Marzieh and Sobolev, Stephan V. and Gerya, Taras V. and Brune, Sascha}, title = {Subduction initiation by Plume-Plateau interaction}, series = {Geochemistry, geophysics, geosystems}, volume = {21}, journal = {Geochemistry, geophysics, geosystems}, number = {8}, publisher = {American Geophysical Union}, address = {Washington}, issn = {1525-2027}, doi = {10.1029/2020GC009119}, pages = {17}, year = {2020}, abstract = {It has recently been demonstrated that the interaction of a mantle plume with sufficiently old oceanic lithosphere can initiate subduction. However, the existence of large lithospheric heterogeneities, such as a buoyant plateau, in proximity to a rising plume head may potentially hinder the formation of a new subduction zone. Here, we investigate this scenario by means of 3-D numerical thermomechanical modeling. We explore how plume-lithosphere interaction is affected by lithospheric age, relative location of plume head and plateau border, and the strength of the oceanic crust. Our numerical experiments suggest four different geodynamic regimes: (a) oceanic trench formation, (b) circular oceanic-plateau trench formation, (c) plateau trench formation, and (d) no trench formation. We show that regardless of the age and crustal strength of the oceanic lithosphere, subduction can initiate when the plume head is either below the plateau border or at a distance less than the plume radius from the plateau edge. Crustal heterogeneity facilitates subduction initiation of old oceanic lithosphere. High crustal strength hampers the formation of a new subduction zone when the plume head is located below a young lithosphere containing a thick and strong plateau. We suggest that plume-plateau interaction in the western margin of the Caribbean could have resulted in subduction initiation when the plume head impinged onto the oceanic lithosphere close to the border between plateau and oceanic crust.}, language = {en} } @article{BahrKolberKabothBahretal.2020, author = {Bahr, Andr{\´e} and Kolber, Gilles and Kaboth-Bahr, Stefanie and Reinhardt, Lutz and Friedrich, Oliver and Pross, J{\"o}rg}, title = {Mega-monsoon variability during the late Triassic}, series = {Sedimentology : the journal of the International Association of Sedimentologists}, volume = {67}, journal = {Sedimentology : the journal of the International Association of Sedimentologists}, number = {2}, publisher = {Wiley-Blackwell}, address = {Oxford}, issn = {0037-0746}, doi = {10.1111/sed.12668}, pages = {951 -- 970}, year = {2020}, abstract = {The formation of the supercontinent Pangaea during the Permo-Triassic gave rise to an extreme monsoonal climate (often termed 'mega-monsoon') that has been documented by numerous palaeo-records. However, considerable debate exists about the role of orbital forcing in causing humid intervals in an otherwise arid climate. To shed new light on the forcing of monsoonal variability in subtropical Pangaea, this study focuses on sediment facies and colour variability of playa and alluvial fan deposits in an outcrop from the late Carnian (ca 225 Ma) in the southern Germanic Basin, south-western Germany. The sediments were deposited against a background of increasingly arid conditions following the humid Carnian Pluvial Event (ca 234 to 232 Ma). The ca 2 center dot 4 Myr long sedimentary succession studied shows a tripartite long-term evolution, starting with a distal mud-flat facies deposited under arid conditions. This phase was followed by a highly variable playa-lake environment that documents more humid conditions and finally a regression of the playa-lake due to a return of arid conditions. The red-green (a*) and lightness (L*) records show that this long-term variability was overprinted by alternating wet/dry cycles driven by orbital precession and ca 405 kyr eccentricity, without significant influence of obliquity. The absence of obliquity in this record indicates that high-latitude forcing played only a minor role in the southern Germanic Basin during the late Carnian. This is different from the subsequent Norian when high-latitude signals became more pronounced, potentially related to the northward drift of the Germanic Basin. The recurring pattern of pluvial events during the late Triassic demonstrates that orbital forcing, in particular eccentricity, stimulated the occurrence and intensity of wet phases. It also highlights the possibility that the Carnian Pluvial Event, although most likely triggered by enhanced volcanic activity, may also have been modified by an orbital stimulus.}, language = {en} } @article{BalischewskiBehrensZehbeetal.2020, author = {Balischewski, Christian and Behrens, Karsten and Zehbe, Kerstin and G{\"u}nter, Christina and Mies, Stefan and Sperlich, Eric and Kelling, Alexandra and Taubert, Andreas}, title = {Ionic liquids with more than one metal}, series = {Chemistry - a European journal}, volume = {26}, journal = {Chemistry - a European journal}, number = {72}, publisher = {Wiley-VCH}, address = {Weinheim}, issn = {0947-6539}, doi = {10.1002/chem.202003097}, pages = {17504 -- 17513}, year = {2020}, abstract = {Thirteen N-butylpyridinium salts, including three monometallic [C4Py](2)[MCl4], nine bimetallic [C4Py](2)[(M1-xMxCl4)-M-a-Cl-b] and one trimetallic compound [C4Py](2)[(M1-y-zMyMz (c) Cl4)-M-a-M-b] (M=Co, Cu, Mn; x=0.25, 0.50 or 0.75 and y=z=0.33), were synthesized and their structure and thermal and electrochemical properties were studied. All compounds are ionic liquids (ILs) with melting points between 69 and 93 degrees C. X-ray diffraction proves that all ILs are isostructural. The conductivity at room temperature is between 10(-4) and 10(-8) S cm(-1). Some Cu-based ILs reach conductivities of 10(-2) S cm(-1), which is, however, probably due to IL dec. This correlates with the optical bandgap measurements indicating the formation of large bandgap semiconductors. At elevated temperatures approaching the melting points, the conductivities reach up to 1.47x10(-1) S cm(-1) at 70 degrees C. The electrochemical stability windows of the ILs are between 2.5 and 3.0 V.}, language = {en} }