@article{RacanoJaraMunozCosentinoetal.2020, author = {Racano, Simone and Jara Mu{\~n}oz, Julius and Cosentino, Domenico and Melnick, Daniel}, title = {Variable quaternary uplift along the Southern Margin of the Central Anatolian Plateau inferred from modeling Marine Terrace sequences}, series = {Tectonics}, volume = {39}, journal = {Tectonics}, number = {12}, publisher = {American Geophysical Union}, address = {Washington}, issn = {0278-7407}, doi = {10.1029/2019TC005921}, pages = {22}, year = {2020}, abstract = {The southern margin of the Central Anatolian Plateau (CAP) records a strong uplift phase after the early Middle Pleistocene, which has been related to the slab break-off of the subducting Arabian plate beneath the Anatolian microplate. During the last 450 kyr the area underwent an uplift phase at a mean rate of similar to 3.2 m/kyr, as suggested by Middle Pleistocene marine sediments exposed at similar to 1,500 m above sea level. These values are significantly higher than the 1.0-1.5 m/kyr estimated since the Late Pleistocene, suggesting temporal variations in uplift rate. To estimate changes in uplift rate during the Pleistocene we studied the marine terraces along the CAP southern margin, mapping the remnants of the platforms and their associated deposits in the field, and used the TerraceM software to identify the position and elevation of associated shoreline angles. We used shoreline angles and the timing of Quaternary marine sedimentation as constrains for a Landscape Evolution Model that simulates wave erosion of an uplifting coast. We applied random optimization algorithms and minimization statistics to find the input parameters that better reproduce the morphology of CAP marine terraces. The best-fitting uplift rate history suggests a significative increase from 1.9 to 3.5 m/kyr between 500 and 200 kyr, followed by an abrupt decrease to 1.4 m/kyr until the present. Our results agree with slab break-off models, which suggest a strong uplift pulse during slab rupture followed by a smoother decrease.}, language = {en} } @article{ZiebarthvonSpechtHeidbachetal.2020, author = {Ziebarth, Malte J. and von Specht, Sebastian and Heidbach, Oliver and Cotton, Fabrice and Anderson, John G.}, title = {Applying conservation of energy to estimate earthquake frequencies from strain rates and stresses}, series = {Journal of geophysical research : Solid earth}, volume = {125}, journal = {Journal of geophysical research : Solid earth}, number = {8}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9313}, doi = {10.1029/2020JB020186}, pages = {25}, year = {2020}, abstract = {Estimating earthquake occurrence rates from the accumulation rate of seismic moment is an established tool of seismic hazard analysis. We propose an alternative, fault-agnostic approach based on the conservation of energy: the Energy-Conserving Seismicity Framework (ENCOS). Working in energy space has the advantage that the radiated energy is a better predictor of the damage potential of earthquake waves than the seismic moment release. In a region, ENCOS balances the stationary power available to cause earthquakes with the long-term seismic energy release represented by the energy-frequency distribution's first moment. Accumulation and release are connected through the average seismic efficiency, by which we mean the fraction of released energy that is converted into seismic waves. Besides measuring earthquakes in energy, ENCOS differs from moment balance essentially in that the energy accumulation rate depends on the total stress in addition to the strain rate tensor. To validate ENCOS, we exemplarily model the energy-frequency distribution around Southern California. We estimate the energy accumulation rate due to tectonic loading assuming poroelasticity and hydrostasis. Using data from the World Stress Map and assuming the frictional limit to estimate the stress tensor, we obtain a power of 0.8 GW. The uncertainty range, 0.3-2.0GW, originates mainly from the thickness of the seismogenic crust, the friction coefficient on preexisting faults, and models of Global Positioning System (GPS) derived strain rates. Based on a Gutenberg-Richter magnitude-frequency distribution, this power can be distributed over a range of energies consistent with historical earthquake rates and reasonable bounds on the seismic efficiency.}, language = {en} } @article{ChengMilsch2020, author = {Cheng, Chaojie and Milsch, Harald}, title = {Permeability variations in illite-bearing sandstone}, series = {Journal of geophysical research : Solid earth}, volume = {125}, journal = {Journal of geophysical research : Solid earth}, number = {9}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9313}, doi = {10.1029/2020JB020122}, pages = {21}, year = {2020}, abstract = {Temperature changes and variations in pore fluid salinity may negatively affect the permeability of clay-bearing sandstones with implications for natural fluid flow and geotechnical applications alike. In this study these factors are investigated for a sandstone dominated by illite as the clay phase. Systematic long-term flow-through experiments were conducted and complemented with comprehensive microstructural investigations and the application of Derjaguin-Landau-Verwey-Overbeek (DLVO) theory to explain mechanistically the observed permeability changes. Initially, sample permeability was not affected by low pore fluid salinity indicating strong attraction of the illite particles to the pore walls as supported by electron microprobe analysis (EMPA). Increasing temperature up to 145 degrees C resulted in an irreversible permeability decrease by 1.5 orders of magnitude regardless of the pore fluid composition (i.e., deionized water and 2 M NaCl solution). Subsequently diluting the high salinity pore fluid to below 0.5 M yielded an additional permeability decline by 1.5 orders of magnitude, both at 145 degrees C and after cooling to room temperature. By applying scanning electron microscopy (SEM) and mercury intrusion porosimetry (MIP) thermo-mechanical pore throat closure and illite particle migration were identified as independently operating mechanisms responsible for observed permeability changes during heating and dilution, respectively. These observations indicate that permeability of illite-bearing sandstones will be impaired by heating and exposure to low salinity pore fluids. In addition, chemically induced permeability variations proved to be path dependent with respect to the applied succession of fluid salinity changes.}, language = {en} } @article{PilzFranckeBaronietal.2020, author = {Pilz, Tobias and Francke, Till and Baroni, Gabriele and Bronstert, Axel}, title = {How to Tailor my process-based hydrological model?}, series = {Water resources research}, volume = {56}, journal = {Water resources research}, number = {8}, publisher = {American Geophysical Union}, address = {Washington}, issn = {0043-1397}, doi = {10.1029/2020WR028042}, pages = {24}, year = {2020}, abstract = {In the field of hydrological modeling, many alternative representations of natural processes exist. Choosing specific process formulations when building a hydrological model is therefore associated with a high degree of ambiguity and subjectivity. In addition, the numerical integration of the underlying differential equations and parametrization of model structures influence model performance. Identifiability analysis may provide guidance by constraining the a priori range of alternatives based on observations. In this work, a flexible simulation environment is used to build an ensemble of semidistributed, process-based hydrological model configurations with alternative process representations, numerical integration schemes, and model parametrizations in an integrated manner. The flexible simulation environment is coupled with an approach for dynamic identifiability analysis. The objective is to investigate the applicability of the framework to identify the most adequate model. While an optimal model configuration could not be clearly distinguished, interesting results were obtained when relating model identifiability with hydro-meteorological boundary conditions. For instance, we tested the Penman-Monteith and Shuttleworth \& Wallace evapotranspiration models and found that the former performs better under wet and the latter under dry conditions. Parametrization of model structures plays a dominant role as it can compensate for inadequate process representations and poor numerical solvers. Therefore, it was found that numerical solvers of high order of accuracy do often, though not necessarily, lead to better model performance. The proposed coupled framework proved to be a straightforward diagnostic tool for model building and hypotheses testing and shows potential for more in-depth analysis of process implementations and catchment functioning.}, language = {en} } @article{DahmStillerMechieetal.2020, author = {Dahm, Torsten and Stiller, Manfred and Mechie, James and Heimann, Sebastian and Hensch, Martin and Woith, Heiko and Schmidt, Bernd and Gabriel, Gerald and Weber, Michael}, title = {Seismological and geophysical signatures of the deep crustal magma systems of the cenozoic volcanic fields Beneath the Eifel, Germany}, series = {Geochemistry, geophysics, geosystems}, volume = {21}, journal = {Geochemistry, geophysics, geosystems}, number = {9}, publisher = {American Geophysical Union}, address = {Washington}, issn = {1525-2027}, doi = {10.1029/2020GC009062}, pages = {21}, year = {2020}, abstract = {The Quaternary volcanic fields of the Eifel (Rhineland-Palatinate, Germany) had their last eruptions less than 13,000 years ago. Recently, deep low-frequency (DLF) earthquakes were detected beneath one of the volcanic fields showing evidence of ongoing magmatic activity in the lower crust and upper mantle. In this work, seismic wide- and steep-angle experiments from 1978/1979 and 1987/1988 are compiled, partially reprocessed and interpreted, together with other data to better determine the location, size, shape, and state of magmatic reservoirs in the Eifel region near the crust-mantle boundary. We discuss seismic evidence for a low-velocity gradient layer from 30-36 km depth, which has developed over a large region under all Quaternary volcanic fields of the Rhenish Massif and can be explained by the presence of partial melts. We show that the DLF earthquakes connect the postulated upper mantle reservoir with the upper crust at a depth of about 8 km, directly below one of the youngest phonolitic volcanic centers in the Eifel, where CO(2)originating from the mantle is massively outgassing. A bright spot in the West Eifel between 6 and 10 km depth represents a Tertiary magma reservoir and is seen as a model for a differentiated reservoir beneath the young phonolitic center today. We find that the distribution of volcanic fields is controlled by the Variscan lithospheric structures and terrane boundaries as a whole, which is reflected by an offset of the Moho depth, a wedge-shaped transparent zone in the lower crust and the system of thrusts over about 120 km length.}, 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{SchoppaSiegVogeletal.2020, author = {Schoppa, Lukas and Sieg, Tobias and Vogel, Kristin and Z{\"o}ller, Gert and Kreibich, Heidi}, title = {Probabilistic flood loss models for companies}, series = {Water resources research}, volume = {56}, journal = {Water resources research}, number = {9}, publisher = {American Geophysical Union}, address = {Washington}, issn = {0043-1397}, doi = {10.1029/2020WR027649}, pages = {19}, year = {2020}, abstract = {Flood loss modeling is a central component of flood risk analysis. Conventionally, this involves univariable and deterministic stage-damage functions. Recent advancements in the field promote the use of multivariable and probabilistic loss models, which consider variables beyond inundation depth and account for prediction uncertainty. Although companies contribute significantly to total loss figures, novel modeling approaches for companies are lacking. Scarce data and the heterogeneity among companies impede the development of company flood loss models. We present three multivariable flood loss models for companies from the manufacturing, commercial, financial, and service sector that intrinsically quantify prediction uncertainty. Based on object-level loss data (n = 1,306), we comparatively evaluate the predictive capacity of Bayesian networks, Bayesian regression, and random forest in relation to deterministic and probabilistic stage-damage functions, serving as benchmarks. The company loss data stem from four postevent surveys in Germany between 2002 and 2013 and include information on flood intensity, company characteristics, emergency response, private precaution, and resulting loss to building, equipment, and goods and stock. We find that the multivariable probabilistic models successfully identify and reproduce essential relationships of flood damage processes in the data. The assessment of model skill focuses on the precision of the probabilistic predictions and reveals that the candidate models outperform the stage-damage functions, while differences among the proposed models are negligible. Although the combination of multivariable and probabilistic loss estimation improves predictive accuracy over the entire data set, wide predictive distributions stress the necessity for the quantification of uncertainty.}, language = {en} } @article{MuldashevSobolev2020, author = {Muldashev, Iskander A. and Sobolev, Stephan}, title = {What controls maximum magnitudes of giant subduction earthquakes?}, series = {Geochemistry, geophysics, geosystems}, volume = {21}, journal = {Geochemistry, geophysics, geosystems}, number = {9}, publisher = {American Geophysical Union}, address = {Washington}, issn = {1525-2027}, doi = {10.1029/2020GC009145}, pages = {15}, year = {2020}, abstract = {Giant earthquakes with magnitudes above 8.5 occur only in subduction zones. Despite the developments made in observing large subduction zone earthquakes with geophysical instruments, the factors controlling the maximum size of these earthquakes are still poorly understood. Previous studies have suggested the importance of slab shape, roughness of the plate interface contact, state of the strain in the upper plate, thickness of sediments filling the trenches, and subduction rate. Here, we present 2-D cross-scale numerical models of seismic cycles for subduction zones with various geometries, subduction channel friction configurations, and subduction rates. We found that low-angle subduction and thick sediments in the subduction channel are the necessary conditions for generating giant earthquakes, while the subduction rate has a negligible effect. We suggest that these key parameters determine the maximum magnitude of a subduction earthquake by controlling the seismogenic zone width and smoothness of the subduction interface. This interpretation supports previous studies that are based upon observations and scaling laws. Our modeling results also suggest that low static friction in the sediment-filled subduction channel results in neutral or moderate compressive deformation in the overriding plate for low-angle subduction zones hosting giant earthquakes. These modeling results agree well with observations for the largest earthquakes. Based on our models we predict maximum magnitudes of subduction earthquakes worldwide, demonstrating the fit to magnitudes of all giant earthquakes of the 20th and 21st centuries and good agreement with the predictions based on statistical analyses of observations.}, language = {en} } @article{SharmaHainzlZoelleretal.2020, author = {Sharma, Shubham and Hainzl, Sebastian and Z{\"o}ller, Gert and Holschneider, Matthias}, title = {Is Coulomb stress the best choice for aftershock forecasting?}, series = {Journal of geophysical research : Solid earth}, volume = {125}, journal = {Journal of geophysical research : Solid earth}, number = {9}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9313}, doi = {10.1029/2020JB019553}, pages = {12}, year = {2020}, abstract = {The Coulomb failure stress (CFS) criterion is the most commonly used method for predicting spatial distributions of aftershocks following large earthquakes. However, large uncertainties are always associated with the calculation of Coulomb stress change. The uncertainties mainly arise due to nonunique slip inversions and unknown receiver faults; especially for the latter, results are highly dependent on the choice of the assumed receiver mechanism. Based on binary tests (aftershocks yes/no), recent studies suggest that alternative stress quantities, a distance-slip probabilistic model as well as deep neural network (DNN) approaches, all are superior to CFS with predefined receiver mechanism. To challenge this conclusion, which might have large implications, we use 289 slip inversions from SRCMOD database to calculate more realistic CFS values for a layered half-space and variable receiver mechanisms. We also analyze the effect of the magnitude cutoff, grid size variation, and aftershock duration to verify the use of receiver operating characteristic (ROC) analysis for the ranking of stress metrics. The observations suggest that introducing a layered half-space does not improve the stress maps and ROC curves. However, results significantly improve for larger aftershocks and shorter time periods but without changing the ranking. We also go beyond binary testing and apply alternative statistics to test the ability to estimate aftershock numbers, which confirm that simple stress metrics perform better than the classic Coulomb failure stress calculations and are also better than the distance-slip probabilistic model.}, language = {en} } @article{KayaDupontNivetProustetal.2020, author = {Kaya, Mustafa Yuecel and Dupont-Nivet, Guillaume and Proust, Jean-No{\"e}l and Roperch, Pierrick and Meijer, Niels and Frieling, Joost and Fioroni, Chiara and Altiner, Sevin{\c{c}} {\"O}zkan and Stoica, Marius and Aminov, Jovid and Mamtimin, Mehmut and Guo, Zhaojie}, title = {Cretaceous evolution of the Central Asian Proto-Paratethys Sea}, series = {Tectonics}, volume = {39}, journal = {Tectonics}, number = {9}, publisher = {American Geophysical Union}, address = {Washington}, issn = {0278-7407}, doi = {10.1029/2019TC005983}, pages = {27}, year = {2020}, abstract = {The timing and mechanisms of the Cretaceous sea incursions into Central Asia are still poorly constrained. We provide a new chronostratigraphic framework based on biostratigraphy and magnetostratigraphy together with detailed paleoenvironmental analyses of Cretaceous records of the proto-Paratethys Sea fluctuations in the Tajik and Tarim basins. The Early Cretaceous marine incursion in the western Tajik Basin was followed by major marine incursions during the Cenomanian (ca. 100 Ma) and Santonian (ca. 86 Ma) that reached far into the eastern Tajik and Tarim basins. These marine incursions were separated by a Turonian-Coniacian (ca. 92-86 Ma) regression. Basin-wide tectonic subsidence analyses imply that the Early Cretaceous sea incursion into the Tajik Basin was related to increased Pamir tectonism. We find that thrusting along the northern edge of the Pamir at ca. 130-90 Ma resulted in increased subsidence in a retro-arc basin setting. This tectonic event and coeval eustatic highstand resulted in the maximum observed geographic extent of the sea during the Cenomanian (ca. 100 Ma). The following Turonian-Coniacian (ca. 92-86 Ma) major regression, driven by eustasy, coincides with a sharp slowdown in tectonic subsidence during the late orogenic unloading period with limited thrusting. The Santonian (ca. 86 Ma) major sea incursion was likely controlled by eustasy as evidenced by the coeval fluctuations in the west Siberian Basin. An early Maastrichtian cooling (ca. 71-70 Ma), potentially connected to global Late Cretaceous trends, is inferred from the replacement of mollusk-rich limestones by bryozoan- and echinoderm-rich limestones.}, language = {en} } @article{LadeiraMarwanDestroFilhoetal.2020, author = {Ladeira, Guenia and Marwan, Norbert and Destro-Filho, Joao-Batista and Ramos, Camila Davi and Lima, Gabriela}, title = {Frequency spectrum recurrence analysis}, series = {Scientific reports}, volume = {10}, journal = {Scientific reports}, number = {1}, publisher = {Nature portfolio}, address = {Berlin}, issn = {2045-2322}, doi = {10.1038/s41598-020-77903-4}, pages = {9}, year = {2020}, abstract = {In this paper, we present the new frequency spectrum recurrence analysis technique by means of electro-encephalon signals (EES) analyses. The technique is suitable for time series analysis with noise and disturbances. EES were collected, and alpha waves of the occipital region were analysed by comparing the signals from participants in two states, eyes open and eyes closed. Firstly, EES were characterized and analysed by means of techniques already known to compare with the results of the innovative technique that we present here. We verified that, standard recurrence quantification analysis by means of EES time series cannot statistically distinguish the two states. However, the new frequency spectrum recurrence quantification exhibit quantitatively whether the participants have their eyes open or closed. In sequence, new quantifiers are created for analysing the recurrence concentration on frequency bands. These analyses show that EES with similar frequency spectrum have different recurrence levels revealing different behaviours of the nervous system. The technique can be used to deepen the study on depression, stress, concentration level and other neurological issues and also can be used in any complex system.}, language = {en} } @article{SoaresYamazakiCnossenetal.2020, author = {Soares, Gabriel Brando and Yamazaki, Yosuke and Cnossen, Ingrid and Matzka, J{\"u}rgen and Pinheiro, Katia J. and Morschhauser, Achim and Alken, Patrick and Stolle, Claudia}, title = {Evolution of the geomagnetic daily variation at Tatuoca, Brazil, From 1957 to 2019}, series = {Journal of geophysical research : Space physics}, volume = {125}, journal = {Journal of geophysical research : Space physics}, number = {9}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9380}, doi = {10.1029/2020JA028109}, pages = {20}, year = {2020}, abstract = {The magnetic equator in the Brazilian region has moved over 1,100 km northward since 1957, passing the geomagnetic observatory Tatuoca (TTB), in northern Brazil, around 2013. We recovered and processed TTB hourly mean values of the geomagnetic field horizontal (H) component from 1957 until 2019, allowing the investigation of long-term changes in the daily variation due to the influence of secular variation, solar activity, season, and lunar phase. The H day-to-day variability and the occurrence of the counter electrojet at TTB were also investigated. Until the 1990s, ionospheric solar quiet currents dominated the quiet-time daily variation at TTB. After 2000, the magnitude of the daily variation became appreciably greater due to the equatorial electrojet (EEJ) contribution. The H seasonal and day-to-day variability increased as the magnetic equator approached, but their amplitudes normalized to the average daily variation remained at similar levels. Meanwhile, the amplitude of the lunar variation, normalized in the same way, increased from 5\% to 12\%. Within the EEJ region, the occurrence rate of the morning counter electrojet (MCEJ) increased with proximity to the magnetic equator, while the afternoon counter electrojet (ACEJ) did not. EEJ currents derived from CHAMP and Swarm satellite data revealed that the MCEJ rate varies with magnetic latitude within the EEJ region while the ACEJ rate is largely constant. Simulations with the Thermosphere-Ionosphere-Electrodynamics General Circulation Model based on different geomagnetic main field configurations suggest that long-term changes in the geomagnetic daily variation at TTB can be attributed to the main field secular variation.}, language = {en} } @article{RodriguezPicedaScheckWenderothGomezDacaletal.2020, author = {Rodriguez Piceda, Constanza and Scheck Wenderoth, Magdalena and Gomez Dacal, Maria Laura and Bott, Judith and Prezzi, Claudia Beatriz and Strecker, Manfred}, title = {Lithospheric density structure of the southern Central Andes constrained by 3D data-integrative gravity modelling}, series = {International journal of earth sciences}, volume = {110}, journal = {International journal of earth sciences}, number = {7}, publisher = {Springer}, address = {New York}, issn = {1437-3254}, doi = {10.1007/s00531-020-01962-1}, pages = {2333 -- 2359}, year = {2020}, abstract = {The southern Central Andes (SCA) (between 27 degrees S and 40 degrees S) is bordered to the west by the convergent margin between the continental South American Plate and the oceanic Nazca Plate. The subduction angle along this margin is variable, as is the deformation of the upper plate. Between 33 degrees S and 35 degrees S, the subduction angle of the Nazca plate increases from sub-horizontal (< 5 degrees) in the north to relatively steep (similar to 30 degrees) in the south. The SCA contain inherited lithological and structural heterogeneities within the crust that have been reactivated and overprinted since the onset of subduction and associated Cenozoic deformation within the Andean orogen. The distribution of the deformation within the SCA has often been attributed to the variations in the subduction angle and the reactivation of these inherited heterogeneities. However, the possible influence that the thickness and composition of the continental crust have had on both short-term and long-term deformation of the SCA is yet to be thoroughly investigated. For our investigations, we have derived density distributions and thicknesses for various layers that make up the lithosphere and evaluated their relationships with tectonic events that occurred over the history of the Andean orogeny and, in particular, investigated the short- and long-term nature of the present-day deformation processes. We established a 3D model of lithosphere beneath the orogen and its foreland (29 degrees S-39 degrees S) that is consistent with currently available geological and geophysical data, including the gravity data. The modelled crustal configuration and density distribution reveal spatial relationships with different tectonic domains: the crystalline crust in the orogen (the magmatic arc and the main orogenic wedge) is thicker (similar to 55 km) and less dense (similar to 2900 kg/m(3)) than in the forearc (similar to 35 km, similar to 2975 kg/m(3)) and foreland (similar to 30 km, similar to 3000 kg/m(3)). Crustal thickening in the orogen probably occurred as a result of stacking of low-density domains, while density and thickness variations beneath the forearc and foreland most likely reflect differences in the tectonic evolution of each area following crustal accretion. No clear spatial relationship exists between the density distribution within the lithosphere and previously proposed boundaries of crustal terranes accreted during the early Paleozoic. Areas with ongoing deformation show a spatial correlation with those areas that have the highest topographic gradients and where there are abrupt changes in the average crustal-density contrast. This suggests that the short-term deformation within the interior of the Andean orogen and its foreland is fundamentally influenced by the crustal composition and the relative thickness of different crustal layers. A thicker, denser, and potentially stronger lithosphere beneath the northern part of the SCA foreland is interpreted to have favoured a strong coupling between the Nazca and South American plates, facilitating the development of a sub-horizontal slab.}, language = {en} } @article{GhaniSobelZeilingeretal.2020, author = {Ghani, Humaad and Sobel, Edward and Zeilinger, Gerold and Glodny, Johannes and Zapata, Sebastian and Irum, Irum}, title = {Palaeozoic and Pliocene tectonic evolution of the Salt Range constrained by low-temperature thermochronology}, series = {Terra nova}, volume = {33}, journal = {Terra nova}, number = {3}, publisher = {Wiley}, address = {Hoboken}, issn = {0954-4879}, doi = {10.1111/ter.12515}, pages = {293 -- 305}, year = {2020}, abstract = {The Salt Range in Pakistan exposes Precambrian to Pleistocene strata outcropping along the Salt Range Thrust (SRT). To better understand the in-situ Cambrian and Pliocene tectonic evolution of the Pakistan Subhimalaya, we have conducted low-temperature thermochronological analysis using apatite (U-Th-Sm)/He and fission track dating. We combine cooling ages from different samples located along the thrust front of the SRT into a thermal model that shows two major cooling events associated with rifting and regional erosion in the Late Palaeozoic and SRT activity since the Pliocene. Our results suggest that the SRT maintained a long-term average shortening rate of similar to 5-6 mm/yr and a high exhumation rate above the SRT ramp since similar to 4 Ma.}, language = {en} } @article{StichMartinMoralesetal.2020, author = {Stich, Daniel and Martin, Rosa and Morales, Jose and Lopez-Comino, Jose Angel and Mancilla, Flor de Lis}, title = {Slip partitioning in the 2016 Alboran Sea earthquake sequence (western Mediterranean)}, series = {Frontiers in Earth Science}, volume = {8}, journal = {Frontiers in Earth Science}, publisher = {Frontiers Media}, address = {Lausanne}, issn = {2296-6463}, doi = {10.3389/feart.2020.587356}, pages = {19}, year = {2020}, abstract = {AM(W)= 5.1 earthquake on January 21st, 2016 marked the beginning of a significant seismic sequence in the southern Alboran Sea, culminating in aM(W)= 6.3 earthquake on January 25th, and continuing with further moderate magnitude earthquakes until March. We use data from 35 seismic broadband stations in Spain, Morocco and Portugal to relocate the seismicity, estimate seismic moment tensors, and isolate regional apparent source time functions for the main earthquake. Relocation and regional moment tensor inversion consistently yield very shallow depths for the majority of events. We obtain 50 moment tensors for the sequence, showing a mixture of strike-slip faulting for the foreshock and the main event and reverse faulting for the major aftershocks. The leading role of reverse focal mechanisms among the aftershocks may be explained by the geometry of the fault network. The mainshock nucleates at a bend along the left-lateral Al-Idrisi fault, introducing local transpression within the transtensional Alboran Basin. The shallow depths of the 2016 Alboran Sea earthquakes may favor slip-partitioning on the involved faults. Apparent source durations for the main event suggest a similar to 21 km long, asymmetric rupture that propagates primarily toward NE into the restraining fault segment, with fast rupture speed of similar to 3.0 km/s. Consistently, the inversion for laterally variable fault displacement situates the main slip in the restraining segment. The partitioning into strike-slip rupture and dip-slip aftershocks confirms a non-optimal orientation of this segment, and suggests that the 2016 event settled a slip deficit from previous ruptures that could not propagate into the stronger restraining segment.}, 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{OuKulikovaYuetal.2020, author = {Ou, Qi and Kulikova, Galina and Yu, Jingxing and Elliott, Austin and Parsons, Bethany and Walker, Richard}, title = {Magnitude of the 1920 Haiyuan earthquake reestimated using seismological and geomorphological methods}, series = {Journal of geophysical research : Solid earth}, volume = {125}, journal = {Journal of geophysical research : Solid earth}, number = {8}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9313}, doi = {10.1029/2019JB019244}, pages = {28}, year = {2020}, abstract = {Reported magnitudes ranging between 7.8 and 8.7 highlight a confusion about the true size of the 1920 Haiyuan earthquake, the largest earthquake recorded in the northeast Tibetan Plateau. We compiled a global data set of previously unlooked-at historical seismograms and performed modern computational analyses on the digitized seismic records. We found the epicenter to be near Haiyuan town and obtained a moment magnitude of M-W=7.90.2. Following traditional approaches, we obtained m(B)=7.90.3 with data from 21 stations and M-S(20)=8.10.2 with data from three stations. Geomorphologically, we mapped the surface rupture and horizontal offsets on high-resolution Pleiades satellite and drone imagery that covered the entire western and middle sections of the 1920 Haiyuan earthquake rupture and compiled offsets reported on the eastern section from field measurements in the 1980s. Careful discrimination between single-event and cumulative offsets suggests average horizontal slips of 3.01.0m on the western section, 4.51.5m on the middle section, and 3.5 +/- 0.5m on the eastern section, indicating a total moment magnitude of M-W=7.8 +/- 0.1. Thus, the seismological and geomorphological results agree within the uncertainties, a weighted average giving a moment magnitude of M-W=7.9 +/- 0.2 for the 1920 Haiyuan earthquake. It is likely that earthquake magnitudes based on the historical M were systematically overestimated.
Plain Language Summary Earthquakes are the main mechanism by which elastic energy accumulating due to tectonic motion is released. As the earthquake magnitude scale is logarithmic, major earthquakes control the bulk of this energy budget and are often the most destructive, like the 1920 Haiyuan earthquake with similar to 230,000 casualties. However, major earthquakes tend to have recurrence periods of several hundred years, longer than our instrumental records. To obtain knowledge of historic major earthquakes, paleoseismologists measure geomorphic offsets and map surface ruptures left by past events and estimate the shaking intensity from historical writings. However, in the case of the Haiyuan earthquake, which happened in the late historic, early instrumental period, the magnitudes reported from these two communities differed significantly. In order to constrain the magnitude of this earthquake for seismic hazard assessment and to reconcile the differences between published magnitudes, we reestimated its magnitude from both newly compiled and digitized seismological records and modern satellite and drone imagery. The results show that the early seismological magnitudes were overestimated, which may affect historical magnitudes systematically. The 1920 Haiyuan earthquake was of a similar magnitude to the 2001 Kokoxili and 2008 Wenchuan earthquakes that also occurred in and around the Tibetan Plateau, instead of more than half a magnitude larger.}, language = {en} } @article{GiarollaVeigaNobreetal.2020, author = {Giarolla, Emanuel and Veiga, Sandro F. and Nobre, Paulo and Silva, Manoel B. and Capistrano, Vinicius B. and Callegare, Andyara O.}, title = {Sea surface height trends in the southern hemisphere oceans simulated by the Brazilian Earth System Model under RCP4.5 and RCP8.5 scenarios}, series = {Journal of southern hemisphere earth systems science}, volume = {70}, journal = {Journal of southern hemisphere earth systems science}, number = {1}, publisher = {CSIRO}, address = {Clayton}, issn = {2206-5865}, doi = {10.1071/ES19042}, pages = {280 -- 289}, year = {2020}, abstract = {The Brazilian Earth System Model (BESM-OA2.5), while simulating the historical period proposed by the fifth phase of the Coupled Model Intercomparison Project (CMIP5), detects an increasing trend in the sea surface height (SSH) on the southern hemisphere oceans relative to that of the pre-industrial era. The increasing trend is accentuated in the CMIP5 RCP4.5 and RCP8.5 future scenarios with higher concentrations of greenhouse gases in the atmosphere. This study sheds light on the sources of such trends in these regions. The results suggest an association with the thermal expansion of the oceans in the upper 700 m due to a gradual warming inflicted by those future scenarios. BESM-OA2.5 presents a surface height increase of 0.11 m in the historical period of 1850-2005. Concerning future projections, BESM-OA2.5 projects SSH increases of 0.14 and 0.23 m (relative to the historical 2005 value) for RCP4.5 and RCP8.5, respectively, by the end of 2100. These increases are predominantly in a band of latitude within 35-60 degrees S in the Atlantic and Indian oceans. The reproducibility of the trend signal detected in the BESM-OA2.5 simulations is confirmed by the results of three other CMIP5 models.}, language = {en} } @article{MerzKuhlickeKunzetal.2020, author = {Merz, Bruno and Kuhlicke, Christian and Kunz, Michael and Pittore, Massimiliano and Babeyko, Andrey and Bresch, David N. and Domeisen, Daniela I. and Feser, Frauke and Koszalka, Inga and Kreibich, Heidi and Pantillon, Florian and Parolai, Stefano and Pinto, Joaquim G. and Punge, Heinz J{\"u}rgen and Rivalta, Eleonora and Schr{\"o}ter, Kai and Strehlow, Karen and Weisse, Ralf and Wurpts, Andreas}, title = {Impact forecasting to support emergency management of natural hazards}, series = {Reviews of geophysics}, volume = {58}, journal = {Reviews of geophysics}, number = {4}, publisher = {American Geophysical Union}, address = {Washington}, issn = {8755-1209}, doi = {10.1029/2020RG000704}, pages = {52}, year = {2020}, abstract = {Forecasting and early warning systems are important investments to protect lives, properties, and livelihood. While early warning systems are frequently used to predict the magnitude, location, and timing of potentially damaging events, these systems rarely provide impact estimates, such as the expected amount and distribution of physical damage, human consequences, disruption of services, or financial loss. Complementing early warning systems with impact forecasts has a twofold advantage: It would provide decision makers with richer information to take informed decisions about emergency measures and focus the attention of different disciplines on a common target. This would allow capitalizing on synergies between different disciplines and boosting the development of multihazard early warning systems. This review discusses the state of the art in impact forecasting for a wide range of natural hazards. We outline the added value of impact-based warnings compared to hazard forecasting for the emergency phase, indicate challenges and pitfalls, and synthesize the review results across hazard types most relevant for Europe.}, language = {en} } @article{MarciszJasseyKosakyanetal.2020, author = {Marcisz, Katarzyna and Jassey, Vincent E. J. and Kosakyan, Anush and Krashevska, Valentyna and Lahr, Daniel J. G. and Lara, Enrique and Lamentowicz, Lukasz and Lamentowicz, Mariusz and Macumber, Andrew and Mazei, Yuri and Mitchell, Edward A. D. and Nasser, Nawaf A. and Patterson, R. Timothy and Roe, Helen M. and Singer, David and Tsyganov, Andrey N. and Fournier, Bertrand}, title = {Testate amoeba functional traits and their use in paleoecology}, series = {Frontiers in Ecology and Evolution}, volume = {8}, journal = {Frontiers in Ecology and Evolution}, publisher = {Frontiers Media}, address = {Lausanne}, issn = {2296-701X}, doi = {10.3389/fevo.2020.575966}, pages = {28}, year = {2020}, abstract = {This review provides a synthesis of current knowledge on the morphological and functional traits of testate amoebae, a polyphyletic group of protists commonly used as proxies of past hydrological changes in paleoecological investigations from peatland, lake sediment and soil archives. A trait-based approach to understanding testate amoebae ecology and paleoecology has gained in popularity in recent years, with research showing that morphological characteristics provide complementary information to the commonly used environmental inferences based on testate amoeba (morpho-)species data. We provide a broad overview of testate amoeba morphological and functional traits and trait-environment relationships in the context of ecology, evolution, genetics, biogeography, and paleoecology. As examples we report upon previous ecological and paleoecological studies that used trait-based approaches, and describe key testate amoebae traits that can be used to improve the interpretation of environmental studies. We also highlight knowledge gaps and speculate on potential future directions for the application of trait-based approaches in testate amoeba research.}, language = {en} } @article{RajaonarisonStampsFishwicketal.2020, author = {Rajaonarison, Tahiry A. and Stamps, D. Sarah and Fishwick, Stewart and Brune, Sascha and Glerum, Anne and Hu, Jiashun}, title = {Numerical modeling of mantle flow beneath Madagascar to constrain upper mantle rheology beneath continental regions}, series = {Journal of geophysical research : Solid earth}, volume = {125}, journal = {Journal of geophysical research : Solid earth}, number = {2}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9313}, doi = {10.1029/2019JB018560}, pages = {23}, year = {2020}, abstract = {Over the past few decades, azimuthal seismic anisotropy measurements have been widely used proxy to study past and present-day deformation of the lithosphere and to characterize convection in the mantle. Beneath continental regions, distinguishing between shallow and deep sources of anisotropy remains difficult due to poor depth constraints of measurements and a lack of regional-scale geodynamic modeling. Here, we constrain the sources of seismic anisotropy beneath Madagascar where a complex pattern cannot be explained by a single process such as absolute plate motion, global mantle flow, or geology. We test the hypotheses that either Edge-Driven Convection (EDC) or mantle flow derived from mantle wind interactions with lithospheric topography is the dominant source of anisotropy beneath Madagascar. We, therefore, simulate two sets of mantle convection models using regional-scale 3-D computational modeling. We then calculate Lattice Preferred Orientation that develops along pathlines of the mantle flow models and use them to calculate synthetic splitting parameters. Comparison of predicted with observed seismic anisotropy shows a good fit in northern and southern Madagascar for the EDC model, but the mantle wind case only fits well in northern Madagascar. This result suggests the dominant control of the measured anisotropy may be from EDC, but the role of localized fossil anisotropy in narrow shear zones cannot be ruled out in southern Madagascar. Our results suggest that the asthenosphere beneath northern and southern Madagascar is dominated by dislocation creep. Dislocation creep rheology may be dominant in the upper asthenosphere beneath other regions of continental lithosphere.}, language = {en} } @article{TofeldeBernhardtGueritetal.2020, author = {Tofelde, Stefanie and Bernhardt, Anne and Guerit, Laure and Romans, Brian W.}, title = {Times Associated With Source-to-Sink Propagation of Environmental Signals During Landscape Transience}, series = {Frontiers in Earth Science}, volume = {9}, journal = {Frontiers in Earth Science}, publisher = {Frontiers Media}, address = {Lausanne, Schweiz}, issn = {2296-6463}, doi = {10.3389/feart.2021.628315}, pages = {1 -- 26}, year = {2020}, abstract = {Sediment archives in the terrestrial and marine realm are regularly analyzed to infer changes in climate, tectonic, or anthropogenic boundary conditions of the past. However, contradictory observations have been made regarding whether short period events are faithfully preserved in stratigraphic archives; for instance, in marine sediments offshore large river systems. On the one hand, short period events are hypothesized to be non-detectable in the signature of terrestrially derived sediments due to buffering during sediment transport along large river systems. On the other hand, several studies have detected signals of short period events in marine records offshore large river systems. We propose that this apparent discrepancy is related to the lack of a differentiation between different types of signals and the lack of distinction between river response times and signal propagation times. In this review, we (1) expand the definition of the term 'signal' and group signals in sub-categories related to hydraulic grain size characteristics, (2) clarify the different types of 'times' and suggest a precise and consistent terminology for future use, and (3) compile and discuss factors influencing the times of signal transfer along sediment routing systems and how those times vary with hydraulic grain size characteristics. Unraveling different types of signals and distinctive time periods related to signal propagation addresses the discrepancies mentioned above and allows a more comprehensive exploration of event preservation in stratigraphy - a prerequisite for reliable environmental reconstructions from terrestrially derived sedimentary records.}, language = {en} } @article{ZapataSobelDelPapaetal.2020, author = {Zapata, Sebastian and Sobel, Edward and Del Papa, Cecilia and Glodny, Johannes}, title = {Upper Plate Controls on the Formation of Broken Foreland Basins in the Andean Retroarc Between 26°S and 28°S}, series = {Geochemistry, Geophysics, Geosystems}, volume = {21}, journal = {Geochemistry, Geophysics, Geosystems}, number = {7}, publisher = {John Wiley \& Sons, Inc.}, address = {New Jersey}, pages = {22}, year = {2020}, abstract = {Marked along-strike changes in stratigraphy, mountain belt morphology, basement exhumation, and deformation styles characterize the Andean retroarc; these changes have previously been related to spatiotemporal variations in the subduction angle. We modeled new apatite fission track and apatite (U-Th-Sm)/He data from nine ranges located between 26 degrees S and 28 degrees S. Using new and previously published data, we constructed a Cretaceous to Pliocene paleogeographic model that delineates a four-stage tectonic evolution: extensional tectonics during the Cretaceous (120-75 Ma), the formation of a broken foreland basin between 55 and 30 Ma, reheating due to burial beneath sedimentary rocks (18-13 Ma), and deformation, exhumation, and surface uplift during the Late Miocene and the Pliocene (13-3 Ma). Our model highlights how preexisting upper plate structures control the deformation patterns of broken foreland basins. Because retroarc deformation predates flat-slab subduction, we propose that slab anchoring may have been the precursor of Eocene-Oligocene compression in the Andean retroarc. Our model challenges models which consider broken foreland basins and retroarc deformation in the NW Argentinian Andes to be directly related to Miocene flat subduction.}, language = {en} } @article{LuDewaldKoutsodendrisetal.2020, author = {Lu, Yin and Dewald, Nico and Koutsodendris, Andreas and Kaboth-Bahr, Stefanie and R{\"o}sler, Wolfgang and Fang, Xiaomin and Pross, J{\"o}rg and Appel, Erwin and Friedrich, Oliver}, title = {Sedimentological evidence for pronounced glacial-interglacial climate fluctuations in NE Tibet in the latest Pliocene to early Pleistocene}, series = {Paleoceanography and Paleoclimatology}, volume = {35}, journal = {Paleoceanography and Paleoclimatology}, number = {5}, publisher = {John Wiley \& Sons, Inc.}, address = {New Jersey}, pages = {16}, year = {2020}, abstract = {The intensification of Northern Hemisphere glaciation (iNHG) and uplift of the Tibetan Plateau have been argued to be among the main drivers of climate change in midlatitude Central Asia during the Pliocene/Pleistocene. While most proxy records that support this hypothesis are from regions outside the Tibetan Plateau (such as from the Chinese Loess Plateau), detailed paleoclimatic information for the plateau itself during that time has yet remained elusive. Here we present a temporally highly resolved (similar to 500 years) sedimentological record from the Qaidam Basin situated on the northeastern Tibetan Plateau that shows pronounced glacial-interglacial climate variability during the interval from 2.7 to 2.1 Ma. Glacial (interglacial) intervals are generally characterized by coarser (finer) grain size, minima (maxima) in organic matter content, and maxima (minima) in carbonate content. Comparison of our results with Earth's orbital parameters and proxy records from the Chinese Loess Plateau suggests that the observed climate fluctuations were mainly driven by changes in the Siberian High/East Asian winter monsoon system as a response to the iNHG. They are further proposed to be enhanced by the topography of the Tibetan Plateau and its impact on the position and intensity of the westerlies.}, language = {en} } @article{SmirnovKronbergDalyetal.2020, author = {Smirnov, Artem G. and Kronberg, Elena A. and Daly, Patrick W. and Aseev, Nikita and Shprits, Yuri and Kellerman, Adam C.}, title = {Adiabatic Invariants Calculations for Cluster Mission: A Long-Term Product for Radiation Belts Studies}, series = {Journal of Geophysical Research: Space Physics}, volume = {125}, journal = {Journal of Geophysical Research: Space Physics}, number = {2}, publisher = {John Wiley \& Sons, Inc.}, address = {New Jersey}, pages = {12}, year = {2020}, abstract = {The Cluster mission has produced a large data set of electron flux measurements in the Earth's magnetosphere since its launch in late 2000. Electron fluxes are measured using Research with Adaptive Particle Imaging Detector (RAPID)/Imaging Electron Spectrometer (IES) detector as a function of energy, pitch angle, spacecraft position, and time. However, no adiabatic invariants have been calculated for Cluster so far. In this paper we present a step-by-step guide to calculations of adiabatic invariants and conversion of the electron flux to phase space density (PSD) in these coordinates. The electron flux is measured in two RAPID/IES energy channels providing pitch angle distribution at energies 39.2-50.5 and 68.1-94.5 keV in nominal mode since 2004. A fitting method allows to expand the conversion of the differential fluxes to the range from 40 to 150 keV. Best data coverage for phase space density in adiabatic invariant coordinates can be obtained for values of second adiabatic invariant, K, similar to 10(2), and values of the first adiabatic invariant mu in the range approximate to 5-20 MeV/G. Furthermore, we describe the production of a new data product "LSTAR," equivalent to the third adiabatic invariant, available through the Cluster Science Archive for years 2001-2018 with 1-min resolution. The produced data set adds to the availability of observations in Earth's radiation belts region and can be used for long-term statistical purposes.}, language = {en} } @article{ZeedenObrehtVeresetal.2020, author = {Zeeden, Christian and Obreht, Igor and Veres, Daniel and Kaboth-Bahr, Stefanie and Hošek, Jan and Marković, Slobodan B. and B{\"o}sken, Janina and Lehmkuhl, Frank and Rolf, Christian and Hambach, Ulrich}, title = {Smoothed millennial-scale palaeoclimatic reference data as unconventional comparison targets: Application to European loess records}, series = {Scientific Reports}, volume = {10}, journal = {Scientific Reports}, number = {1}, publisher = {Springer Nature}, address = {Berlin}, pages = {13}, year = {2020}, abstract = {Millennial-scale palaeoclimate variability has been documented in various terrestrial and marine palaeoclimate proxy records throughout the Northern Hemisphere for the last glacial cycle. Its clear expression and rapid shifts between different states of climate (Greenland Interstadials and Stadials) represents a correlation tool beyond the resolution of e.g. luminescence dating, especially relevant for terrestrial deposits. Usually, comparison of terrestrial proxy datasets and the Greenland ice cores indicates a complex expression of millennial-scale climate variability as recorded in terrestrial geoarchives including loess. Loess is the most widespread terrestrial geoarchive of the Quaternary and especially widespread over Eurasia. However, loess often records a smoothed representation of millennial-scale variability without all fidelity when compared to the Greenland data, this being a relevant limiting feature in integrating loess with other palaeoclimate records. To better understand the loess proxy-response to millennial-scale climate variability, we simulate a proxy signal smoothing by natural processes through application of low-pass filters of delta O-18 data from Greenland, a high-resolution palaeoclimate reference record, alongside speleothem isotope records from the Black Sea-Mediterranean region. We show that low-pass filters represent rather simple models for better constraining the expression of millennial-scale climate variability in low sedimentation environments, and in sediments where proxy-response signals are most likely affected by natural smoothing (by e.g. bioturbation). Interestingly, smoothed datasets from Greenland and the Black Sea-Mediterranean region are most similar in the last similar to 15 ka and between similar to 50-30 ka. Between similar to 30-15 ka, roughly corresponding to the Last Glacial Maximum and the deglaciation, the records show dissimilarities, challenging the construction of robust correlative time-scales in this age range. From our analysis it becomes apparent that patterns of palaeoclimate signals in loess-palaeosol sequences often might be better explained by smoothed Greenland reference data than the original high-resolution Greenland dataset, or other reference data. This opens the possibility to better assess the temporal resolution and palaeoclimate potential of loess-palaeosol sequences in recording supra-regional climate patterns, as well as to securely integrate loess with other chronologically better-resolved palaeoclimate records.}, language = {en} } @article{NaliboffGlerumBruneetal.2020, author = {Naliboff, John B. and Glerum, Anne and Brune, Sascha and P{\´e}ron-Pinvidic, G. and Wrona, Thilo}, title = {Development of 3-D rift heterogeneity through fault network evolution}, series = {Geophysical Research Letters}, volume = {47}, journal = {Geophysical Research Letters}, number = {13}, publisher = {John Wiley \& Sons, Inc.}, address = {New Jersey}, pages = {11}, year = {2020}, abstract = {Observations of rift and rifted margin architecture suggest that significant spatial and temporal structural heterogeneity develops during the multiphase evolution of continental rifting. Inheritance is often invoked to explain this heterogeneity, such as preexisting anisotropies in rock composition, rheology, and deformation. Here, we use high-resolution 3-D thermal-mechanical numerical models of continental extension to demonstrate that rift-parallel heterogeneity may develop solely through fault network evolution during the transition from distributed to localized deformation. In our models, the initial phase of distributed normal faulting is seeded through randomized initial strength perturbations in an otherwise laterally homogeneous lithosphere extending at a constant rate. Continued extension localizes deformation onto lithosphere-scale faults, which are laterally offset by tens of km and discontinuous along-strike. These results demonstrate that rift- and margin-parallel heterogeneity of large-scale fault patterns may in-part be a natural byproduct of fault network coalescence.}, language = {en} } @article{CoesfeldKuesterKuechlyetal.2020, author = {Coesfeld, Jacqueline and Kuester, Theres and Kuechly, Helga U. and Kyba, Christopher C. M.}, title = {Reducing variability and removing natural light from nighttime satellite imagery: A case study using the VIIRS DNB}, series = {Sensors}, volume = {20}, journal = {Sensors}, publisher = {MDPI}, address = {Basel}, pages = {13}, year = {2020}, abstract = {Temporal variation of natural light sources such as airglow limits the ability of night light sensors to detect changes in small sources of artificial light (such as villages). This study presents a method for correcting for this effect globally, using the satellite radiance detected from regions without artificial light emissions. We developed a routine to define an approximate grid of locations worldwide that do not have regular light emission. We apply this method with a 5 degree equally spaced global grid (total of 2016 individual locations), using data from the Visible Infrared Imaging Radiometer Suite (VIIRS) Day-Night Band (DNB). This code could easily be adapted for other future global sensors. The correction reduces the standard deviation of data in the Earth Observation Group monthly DNB composites by almost a factor of two. The code and datasets presented here are available under an open license by GFZ Data Services, and are implemented in the Radiance Light Trends web application.}, language = {en} } @article{RolfesWilhelm2020, author = {Rolfes, Manfred and Wilhelm, Jan Lorenz}, title = {Urban complexity and urban change}, series = {Urban Change and Citizenship in Times of Crisis}, journal = {Urban Change and Citizenship in Times of Crisis}, publisher = {Routledge}, address = {London}, isbn = {978-0-429-26226-5}, pages = {135 -- 155}, year = {2020}, abstract = {The chapter illustrates how cities can be understood from a system-theory perspective as complex social systems. It argues that the classical and linear intervention methods are often no (longer) suitable for the complex structure, temporal dynamics and multifaceted processuality of urban development. It offers a systemic and systems theory-inspired method as an alternative approach, which allows for extended possibilities that are more appropriate for dealing with urban development processes. The method was developed on the basis of practical experience and theoretical insights. The approach should demonstrate for local decision-makers potential areas of activity for organising urban changes through co-design.}, language = {en} } @article{EngelsMedeirosAxfordetal.2020, author = {Engels, Stefan and Medeiros, Andrew S. and Axford, Yarrow and Brooks, Steve and Heiri, Oliver and Luoto, Tomi P. and Nazarova, Larisa B. and Porinchu, David F. and Quinlan, Roberto and Self, Angela E.}, title = {Temperature change as a driver of spatial patterns and long-term trends in chironomid (Insecta: Diptera) diversity}, series = {Global change biology}, volume = {26}, journal = {Global change biology}, number = {3}, publisher = {Wiley}, address = {Hoboken}, issn = {1354-1013}, doi = {10.1111/gcb.14862}, pages = {1155 -- 1169}, year = {2020}, abstract = {Anthropogenic activities have led to a global decline in biodiversity, and monitoring studies indicate that both insect communities and wetland ecosystems are particularly affected. However, there is a need for long-term data (over centennial or millennial timescales) to better understand natural community dynamics and the processes that govern the observed trends. Chironomids (Insecta: Diptera: Chironomidae) are often the most abundant insects in lake ecosystems, sensitive to environmental change, and, because their larval exoskeleton head capsules preserve well in lake sediments, they provide a unique record of insect community dynamics through time. Here, we provide the results of a metadata analysis of chironomid diversity across a range of spatial and temporal scales. First, we analyse spatial trends in chironomid diversity using Northern Hemispheric data sets overall consisting of 837 lakes. Our results indicate that in most of our data sets, summer temperature (T-jul) is strongly associated with spatial trends in modern-day chironomid diversity. We observe a strong increase in chironomid alpha diversity with increasing T-jul in regions with present-day T-jul between 2.5 and 14 degrees C. In some areas with T-jul > 14 degrees C, chironomid diversity stabilizes or declines. Second, we demonstrate that the direction and amplitude of change in alpha diversity in a compilation of subfossil chironomid records spanning the last glacial-interglacial transition (similar to 15,000-11,000 years ago) are similar to those observed in our modern data. A compilation of Holocene records shows that during phases when the amplitude of temperature change was small, site-specific factors had a greater influence on the chironomid fauna obscuring the chironomid diversity-temperature relationship. Our results imply expected overall chironomid diversity increases in colder regions such as the Arctic under sustained global warming, but with complex and not necessarily predictable responses for individual sites.}, language = {en} } @article{RepaschWittmannScheingrossetal.2020, author = {Repasch, Marisa and Wittmann, Hella and Scheingross, Joel S. and Sachse, Dirk and Szupiany, Ricardo and Orfeo, Oscar and Fuchs, Margret and Hovius, Niels}, title = {Sediment Transit Time and Floodplain Storage Dynamics in Alluvial Rivers Revealed by Meteoric 10Be}, series = {Journal of Geophysical Research: Earth Surface}, volume = {125}, journal = {Journal of Geophysical Research: Earth Surface}, publisher = {Wiley}, address = {Hoboken, NJ}, issn = {2169-9011}, doi = {10.1029/2019JF005419}, pages = {19}, year = {2020}, abstract = {Quantifying the time scales of sediment transport and storage through river systems is fundamental for understanding weathering processes, biogeochemical cycling, and improving watershed management, but measuring sediment transit time is challenging. Here we provide the first systematic test of measuring cosmogenic meteoric Beryllium-10 (10Bem) in the sediment load of a large alluvial river to quantify sediment transit times. We take advantage of a natural experiment in the Rio Bermejo, a lowland alluvial river traversing the east Andean foreland basin in northern Argentina. This river has no tributaries along its trunk channel for nearly 1,300 km downstream from the mountain front. We sampled suspended sediment depth profiles along the channel and measured the concentrations of 10Bem in the chemically extracted grain coatings. We calculated depth-integrated 10Bem concentrations using sediment flux data and found that 10Bem concentrations increase 230\% from upstream to downstream, indicating a mean total sediment transit time of 8.4 ± 2.2 kyr. Bulk sediment budget-based estimates of channel belt and fan storage times suggest that the 10Bem tracer records mixing of old and young sediment reservoirs. On a reach scale, 10Bem transit times are shorter where the channel is braided and superelevated above the floodplain, and longer where the channel is incised and meandering, suggesting that transit time is controlled by channel morphodynamics. This is the first systematic application of 10Bem as a sediment transit time tracer and highlights the method's potential for inferring sediment routing and storage dynamics in large river systems.}, language = {en} } @article{SmithBookhagen2020, author = {Smith, Taylor and Bookhagen, Bodo}, title = {Assessing Multi-Temporal Snow-Volume Trends in High Mountain Asia From 1987 to 2016 Using High-Resolution Passive Microwave Data}, series = {Frontiers in Earth Science}, volume = {8}, journal = {Frontiers in Earth Science}, publisher = {Frontiers Media}, address = {Lausanne}, issn = {2296-6463}, doi = {10.3389/feart.2020.559175}, pages = {13}, year = {2020}, abstract = {High Mountain Asia (HMA) is dependent upon both the amount and timing of snow and glacier meltwater. Previous model studies and coarse resolution (0.25° × 0.25°, ∼25 km × 25 km) passive microwave assessments of trends in the volume and timing of snowfall, snowmelt, and glacier melt in HMA have identified key spatial and seasonal heterogeneities in the response of snow to changes in regional climate. Here we use recently developed, continuous, internally consistent, and high-resolution passive microwave data (3.125 km × 3.125 km, 1987-2016) from the special sensor microwave imager instrument family to refine and extend previous estimates of changes in the snow regime of HMA. We find an overall decline in snow volume across HMA; however, there exist spatially contiguous regions of increasing snow volume—particularly during the winter season in the Pamir, Karakoram, Hindu Kush, and Kunlun Shan. Detailed analysis of changes in snow-volume trends through time reveal a large step change from negative trends during the period 1987-1997, to much more positive trends across large regions of HMA during the periods 1997-2007 and 2007-2016. We also find that changes in high percentile monthly snow-water volume exhibit steeper trends than changes in low percentile snow-water volume, which suggests a reduction in the frequency of high snow-water volumes in much of HMA. Regions with positive snow-water storage trends generally correspond to regions of positive glacier mass balances.}, language = {en} } @article{RungeGrosse2020, author = {Runge, Alexandra and Grosse, Guido}, title = {Mosaicking Landsat and Sentinel-2 Data to Enhance LandTrendr Time Series Analysis in Northern High Latitude Permafrost Regions}, series = {Remote Sensing}, volume = {12}, journal = {Remote Sensing}, number = {15}, publisher = {MDPI}, address = {Basel}, issn = {2072-4292}, doi = {10.3390/rs12152471}, pages = {23}, year = {2020}, abstract = {Permafrost is warming in the northern high latitudes, inducing highly dynamic thaw-related permafrost disturbances across the terrestrial Arctic. Monitoring and tracking of permafrost disturbances is important as they impact surrounding landscapes, ecosystems and infrastructure. Remote sensing provides the means to detect, map, and quantify these changes homogeneously across large regions and time scales. Existing Landsat-based algorithms assess different types of disturbances with similar spatiotemporal requirements. However, Landsat-based analyses are restricted in northern high latitudes due to the long repeat interval and frequent clouds, in particular at Arctic coastal sites. We therefore propose to combine Landsat and Sentinel-2 data for enhanced data coverage and present a combined annual mosaic workflow, expanding currently available algorithms, such as LandTrendr, to achieve more reliable time series analysis. We exemplary test the workflow for twelve sites across the northern high latitudes in Siberia. We assessed the number of images and cloud-free pixels, the spatial mosaic coverage and the mosaic quality with spectral comparisons. The number of available images increased steadily from 1999 to 2019 but especially from 2016 onward with the addition of Sentinel-2 images. Consequently, we have an increased number of cloud-free pixels even under challenging environmental conditions, which then serve as the input to the mosaicking process. In a comparison of annual mosaics, the Landsat+Sentinel-2 mosaics always fully covered the study areas (99.9-100 \%), while Landsat-only mosaics contained data-gaps in the same years, only reaching coverage percentages of 27.2 \%, 58.1 \%, and 69.7 \% for Sobo Sise, East Taymyr, and Kurungnakh in 2017, respectively. The spectral comparison of Landsat image, Sentinel-2 image, and Landsat+Sentinel-2 mosaic showed high correlation between the input images and mosaic bands (e.g., for Kurungnakh 0.91-0.97 between Landsat and Landsat+Sentinel-2 mosaic and 0.92-0.98 between Sentinel-2 and Landsat+Sentinel-2 mosaic) across all twelve study sites, testifying good quality mosaic results. Our results show that especially the results for northern, coastal areas was substantially improved with the Landsat+Sentinel-2 mosaics. By combining Landsat and Sentinel-2 data we accomplished to create reliably high spatial resolution input mosaics for time series analyses. Our approach allows to apply a high temporal continuous time series analysis to northern high latitude permafrost regions for the first time, overcoming substantial data gaps, and assess permafrost disturbance dynamics on an annual scale across large regions with algorithms such as LandTrendr by deriving the location, timing and progression of permafrost thaw disturbances}, 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{MilewskiChabrillatBookhagen2020, author = {Milewski, Robert and Chabrillat, Sabine and Bookhagen, Bodo}, title = {Analyses of Namibian Seasonal Salt Pan Crust Dynamics and Climatic Drivers Using Landsat 8 Time-Series and Ground Data}, series = {Remote Sensing}, journal = {Remote Sensing}, number = {3}, publisher = {MDPI}, address = {Basel}, issn = {2072-4292}, doi = {10.3390/rs12030474}, pages = {24}, year = {2020}, abstract = {Salt pans are highly dynamic environments that are difficult to study by in situ methods because of their harsh climatic conditions and large spatial areas. Remote sensing can help to elucidate their environmental dynamics and provide important constraints regarding their sedimentological, mineralogical, and hydrological evolution. This study utilizes spaceborne multitemporal multispectral optical data combined with spectral endmembers to document spatial distribution of surface crust types over time on the Omongwa pan located in the Namibian Kalahari. For this purpose, 49 surface samples were collected for spectral and mineralogical characterization during three field campaigns (2014-2016) reflecting different seasons and surface conditions of the salt pan. An approach was developed to allow the spatiotemporal analysis of the salt pan crust dynamics in a dense time-series consisting of 77 Landsat 8 cloud-free scenes between 2014 and 2017, covering at least three major wet-dry cycles. The established spectral analysis technique Sequential Maximum Angle Convex Cone (SMACC) extraction method was used to derive image endmembers from the Landsat time-series stack. Evaluation of the extracted endmember set revealed that the multispectral data allowed the differentiation of four endmembers associated with mineralogical mixtures of the crust's composition in dry conditions and three endmembers associated with flooded or muddy pan conditions. The dry crust endmember spectra have been identified in relation to visible, near infrared, and short-wave infrared (VNIR-SWIR) spectroscopy and X-ray diffraction (XRD) analyses of the collected surface samples. According these results, the spectral endmembers are interpreted as efflorescent halite crust, mixed halite-gypsum crust, mixed calcite quartz sepiolite crust, and gypsum crust. For each Landsat scene the spatial distribution of these crust types was mapped with the Spectral Angle Mapper (SAM) method and significant spatiotemporal dynamics of the major surface crust types were observed. Further, the surface crust dynamics were analyzed in comparison with the pan's moisture regime and other climatic parameters. The results show that the crust dynamics are mainly driven by flooding events in the wet season, but are also influenced by temperature and aeolian activity in the dry season. The approach utilized in this study combines the advantages of multitemporal satellite data for temporal event characterization with advantages from hyperspectral methods for the image and ground data analyses that allow improved mineralogical differentiation and characterization.}, language = {en} } @article{NardiniRybackiKrauseetal.2020, author = {Nardini, Livia and Rybacki, Erik and Krause, Michael and Morales, Luiz F.G. and Dresen, Georg}, title = {Control of the geometric arrangement of material heterogeneities on strain localization at the brittle-to-ductile transition in experimentally deformed carbonate rocks}, series = {Journal of Structural Geology}, volume = {135}, journal = {Journal of Structural Geology}, publisher = {Pergamon Press}, address = {Oxford ; Frankfurt, M.}, issn = {0191-8141}, doi = {10.1016/j.jsg.2020.104038}, year = {2020}, abstract = {Triaxial high temperature (900 °C) deformation experiments were conducted at constant strain rate in a Paterson-type deformation apparatus on cylinders of Carrara marble with two right or left stepping, non-overlapping weak inclusions of Solnhofen limestone, oriented at 45° to the cylinders' longitudinal axes. Applying different values of confinement (30, 50, 100 and 300 MPa) we induced various amounts of brittle deformation in the marble matrix and investigated the effect of brittle precursors on the initiation and development of heterogeneity-induced high temperature shear zones. Viscosity contrast between the matrix and the inclusions induces local stress concentration at the tips of these latter. The initial arrangement of the inclusions results in either an overpressured (contractional) or underpressured (extensional) domain in the step-over region of the sample. At low confinement (30 and 50 MPa) abundant brittle deformation is observed, but the spatial distribution of microfractures is dependent on the kinematics of the step-over region: microcracks occur either along the shearing plane between inclusions (in extensional bridge samples), or broadly distributed outside the step-over region (contractional bridge samples). Accordingly, ductile deformation localizes along the inclusions plane in the extensional bridge samples as opposed to distributing over large areas of the matrix in the contractional bridge samples. If microcracking is suppressed (high confinement), strain is accommodated by viscous creep and strain progressively de-localizes in extensional bridge samples. Our experiments demonstrate that brittle precursors enhance the degree of localization in the ductile deformation regime, but only if the interaction of pre-existing heterogeneities induces an extensional mean stress regime in between.}, language = {en} } @article{DommainAndamaMcDonoughetal.2020, author = {Dommain, Ren{\´e} and Andama, Morgan and McDonough, Molly M. and Prado, Natalia A. and Goldhammer, Tobias and Potts, Richard and Maldonado, Jes{\´u}s E. and Nkurunungi, John Bosco and Campana, Michael G.}, title = {The Challenges of Reconstructing Tropical Biodiversity With Sedimentary Ancient DNA}, series = {Frontiers in Ecology and Evolution}, volume = {8}, journal = {Frontiers in Ecology and Evolution}, publisher = {Frontiers Media}, address = {Lausanne}, issn = {2296-701X}, doi = {10.3389/fevo.2020.00218}, pages = {26}, year = {2020}, abstract = {Sedimentary ancient DNA has been proposed as a key methodology for reconstructing biodiversity over time. Yet, despite the concentration of Earth's biodiversity in the tropics, this method has rarely been applied in this region. Moreover, the taphonomy of sedimentary DNA, especially in tropical environments, is poorly understood. This study elucidates challenges and opportunities of sedimentary ancient DNA approaches for reconstructing tropical biodiversity. We present shotgun-sequenced metagenomic profiles and DNA degradation patterns from multiple sediment cores from Mubwindi Swamp, located in Bwindi Impenetrable Forest (Uganda), one of the most diverse forests in Africa. We describe the taxonomic composition of the sediments covering the past 2200 years and compare the sedimentary DNA data with a comprehensive set of environmental and sedimentological parameters to unravel the conditions of DNA degradation. Consistent with the preservation of authentic ancient DNA in tropical swamp sediments, DNA concentration and mean fragment length declined exponentially with age and depth, while terminal deamination increased with age. DNA preservation patterns cannot be explained by any environmental parameter alone, but age seems to be the primary driver of DNA degradation in the swamp. Besides degradation, the presence of living microbial communities in the sediment also affects DNA quantity. Critically, 92.3\% of our metagenomic data of a total 81.8 million unique, merged reads cannot be taxonomically identified due to the absence of genomic references in public databases. Of the remaining 7.7\%, most of the data (93.0\%) derive from Bacteria and Archaea, whereas only 0-5.8\% are from Metazoa and 0-6.9\% from Viridiplantae, in part due to unbalanced taxa representation in the reference data. The plant DNA record at ordinal level agrees well with local pollen data but resolves less diversity. Our animal DNA record reveals the presence of 41 native taxa (16 orders) including Afrotheria, Carnivora, and Ruminantia at Bwindi during the past 2200 years. Overall, we observe no decline in taxonomic richness with increasing age suggesting that several-thousand-year-old information on past biodiversity can be retrieved from tropical sediments. However, comprehensive genomic surveys of tropical biota need prioritization for sedimentary DNA to be a viable methodology for future tropical biodiversity studies.}, 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{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{CrisologoHeistermann2020, author = {Crisologo, Irene and Heistermann, Maik}, title = {Using ground radar overlaps to verify the retrieval of calibration bias estimates from spaceborne platforms}, series = {Atmospheric Measurement Techniques}, volume = {13}, journal = {Atmospheric Measurement Techniques}, number = {2}, publisher = {Copernicus Publications}, address = {G{\"o}ttingen}, issn = {1867-1381}, doi = {10.5194/amt-13-645-2020}, pages = {645 -- 659}, year = {2020}, abstract = {Many institutions struggle to tap into the potential of their large archives of radar reflectivity: these data are often affected by miscalibration, yet the bias is typically unknown and temporally volatile. Still, relative calibration techniques can be used to correct the measurements a posteriori. For that purpose, the usage of spaceborne reflectivity observations from the Tropical Rainfall Measuring Mission (TRMM) and Global Precipitation Measurement (GPM) platforms has become increasingly popular: the calibration bias of a ground radar (GR) is estimated from its average reflectivity difference to the spaceborne radar (SR). Recently, Crisologo et al. (2018) introduced a formal procedure to enhance the reliability of such estimates: each match between SR and GR observations is assigned a quality index, and the calibration bias is inferred as a quality-weighted average of the differences between SR and GR. The relevance of quality was exemplified for the Subic S-band radar in the Philippines, which is greatly affected by partial beam blockage. The present study extends the concept of quality-weighted averaging by accounting for path-integrated attenuation (PIA) in addition to beam blockage. This extension becomes vital for radars that operate at the C or X band. Correspondingly, the study setup includes a C-band radar that substantially overlaps with the S-band radar. Based on the extended quality-weighting approach, we retrieve, for each of the two ground radars, a time series of calibration bias estimates from suitable SR overpasses. As a result of applying these estimates to correct the ground radar observations, the consistency between the ground radars in the region of overlap increased substantially. Furthermore, we investigated if the bias estimates can be interpolated in time, so that ground radar observations can be corrected even in the absence of prompt SR overpasses. We found that a moving average approach was most suitable for that purpose, although limited by the absence of explicit records of radar maintenance operations.}, language = {en} }