@article{OeztuerkBozzolanHolcombeetal.2022, author = {{\"O}zt{\"u}rk, Ugur and Bozzolan, Elisa and Holcombe, Elizabeth A. and Shukla, Roopam and Pianosi, Francesca and Wagener, Thorsten}, title = {How climate change and unplanned urban sprawl bring more landslides}, series = {Nature : the international weekly journal of science}, volume = {608}, journal = {Nature : the international weekly journal of science}, number = {7922}, publisher = {Nature portfolio}, address = {Berlin}, issn = {0028-0836}, doi = {10.1038/d41586-022-02141-9}, pages = {262 -- 265}, year = {2022}, abstract = {More settlements will suffer as heavy rains and unregulated construction destabilize slopes in the tropics, models show.}, language = {en} } @article{ZuhrDolmanHoetal.2022, author = {Zuhr, Alexandra M. and Dolman, Andrew M. and Ho, Sze Ling and Groeneveld, Jeroen and Loewemark, Ludvig and Grotheer, Hendrik and Su, Chih-Chieh and Laepple, Thomas}, title = {Age-heterogeneity in marine sediments revealed by three-dimensional high-resolution radiocarbon measurements}, series = {Frontiers in Earth Science}, volume = {10}, journal = {Frontiers in Earth Science}, publisher = {Frontiers Media}, address = {Lausanne}, issn = {2296-6463}, doi = {10.3389/feart.2022.871902}, pages = {15}, year = {2022}, abstract = {Marine sedimentary archives are routinely used to reconstruct past environmental changes. In many cases, bioturbation and sedimentary mixing affect the proxy time-series and the age-depth relationship. While idealized models of bioturbation exist, they usually assume homogeneous mixing, thus that a single sample is representative for the sediment layer it is sampled from. However, it is largely unknown to which extent this assumption holds for sediments used for paleoclimate reconstructions. To shed light on 1) the age-depth relationship and its full uncertainty, 2) the magnitude of mixing processes affecting the downcore proxy variations, and 3) the representativity of the discrete sample for the sediment layer, we designed and performed a case study on South China Sea sediment material which was collected using a box corer and which covers the last glacial cycle. Using the radiocarbon content of foraminiferal tests as a tracer of time, we characterize the spatial age-heterogeneity of sediments in a three-dimensional setup. In total, 118 radiocarbon measurements were performed on defined small- and large-volume bulk samples ( similar to 200 specimens each) to investigate the horizontal heterogeneity of the sediment. Additionally, replicated measurements on small numbers of specimens (10 x 5 specimens) were performed to assess the heterogeneity within a sample volume. Visual assessment of X-ray images and a quantitative assessment of the mixing strength show typical mixing from bioturbation corresponding to around 10 cm mixing depth. Notably, our 3D radiocarbon distribution reveals that the horizontal heterogeneity (up to 1,250 years), contributing to the age uncertainty, is several times larger than the typically assumed radiocarbon based age-model error (single errors up to 250 years). Furthermore, the assumption of a perfectly bioturbated layer with no mixing underneath is not met. Our analysis further demonstrates that the age-heterogeneity might be a function of sample size; smaller samples might contain single features from the incomplete mixing and are thus less representative than larger samples. We provide suggestions for future studies, optimal sampling strategies for quantitative paleoclimate reconstructions and realistic uncertainty in age models, as well as discuss possible implications for the interpretation of paleoclimate records.}, language = {en} } @article{ZhangCaoYuanetal.2022, author = {Zhang, Di and Cao, Kai and Yuan, Xiaoping and Wang, Guocan and van der Beek, Pieter A.}, title = {Late Oligocene-early Miocene origin of the First Bend of the Yangtze River explained by thrusting-induced river reorganization}, series = {Geomorphology}, volume = {411}, journal = {Geomorphology}, publisher = {Elsevier Science}, address = {Amsterdam [u.a.]}, issn = {0169-555X}, doi = {10.1016/j.geomorph.2022.108303}, pages = {13}, year = {2022}, abstract = {The origin of the First Bend of the Yangtze River is key to understanding the birth of the modern Yangtze River. Despite considerable efforts, the timing and mechanism of formation of the First Bend remain highly debated. Inverse river-profile modeling of three tributaries (Chongjiang, Lima, and Gudu) of the Jinsha River, integrated with regional tectonic and geomorphic interpretations, allows the onset of incision at the First Bend to be constrained to 28-20 Ma. The spatio-temporal coincidence of initial river incision and activity of Yulong Thrust Belt in southeastern Tibet highlights thrusting to be fundamental in reshaping the pre-existing stream network at the First Bend. These results enable us to reinterpret a change in sedimentary environment from a braided river to a swamp-like lake in the Jianchuan Basin south of the First Bend, recording the destruction of the hypothesized southwards-flowing paleo-Jinsha and Shuiluo Rivers at ~36-35 Ma by magmatism. During the late Oligoceneearly Miocene, the paleo-Shuiluo River was diverted to the north by focused rock uplift due to thrusting along the Yulong Thrust Belt, which also led to exhumation of the Jianchuan Basin. Diversion of the paleo-Shuiluo River can be explained by capture from a downstream river in the footwall of the Yulong Thrust Belt. Subsequent rapid headward erosion, that was caused by thrusting-induced drop of local base level, is recorded by upstream younging ages for the onset of incision and led to the formation of the First Bend. The combination of new ages for the onset of incision at 28-20 Ma at the First Bend and younger ages upstream indicates northwards expansion of the Jinsha River at a rate of 62 +/- 18 mm/yr. Our results suggest that the origin of the First Bend was likely triggered by thrusting at 28-20 Ma, after which the Yangtze River formed.}, language = {en} } @article{YuanJiaoDupontNivetetal.2022, author = {Yuan, Xiaoping P. and Jiao, Ruohong and Dupont-Nivet, Guillaume and Shen, Xiaoming}, title = {Southeastern Tibetan Plateau growth revealed by inverse analysis of landscape evolution model}, series = {Geophysical research letters}, volume = {49}, journal = {Geophysical research letters}, number = {10}, publisher = {American Geophysical Union}, address = {Washington}, issn = {0094-8276}, doi = {10.1029/2021GL097623}, pages = {10}, year = {2022}, abstract = {The Cenozoic history of the Tibetan Plateau topography is critical for understanding the evolution of the Indian-Eurasian collision, climate, and biodiversity. However, the long-term growth and landscape evolution of the Tibetan Plateau remain ambiguous, it remains unclear if plateau uplift occurred soon after the India-Asia collision in the Paleogene (similar to 50-25 Ma) or later in the Neogene (similar to 20-5 Ma). Here, we reproduce the uplift history of the southeastern Tibetan Plateau using a 2D landscape evolution model, which simultaneously solves fluvial erosion and sediment transport processes in the drainage basins of the Three Rivers region (Yangtze, Mekong, and Salween Rivers). Our model was optimized through a formal inverse analysis with 20,000 forward simulations, which aims to reconcile the transient states of the present-day river profiles. The results, compared to existing paleoelevation and thermochronologic data, suggest initially low elevations (similar to 300-500 m) during the Paleogene, followed by a gradual southeastward propagation of topographic uplift of the plateau margin.}, language = {en} } @article{YenvonSpechtLinetal.2022, author = {Yen, Ming-Hsuan and von Specht, Sebastian and Lin, Yen-Yu and Cotton, Fabrice and Ma, Kuo-Fong}, title = {Within- and between-event variabilities of strong-velocity pulses of moderate earthquakes within dense seismic arrays}, series = {Bulletin of the Seismological Society of America}, volume = {112}, journal = {Bulletin of the Seismological Society of America}, number = {1}, publisher = {Seismological Society of America}, address = {El Cerito, Calif.}, issn = {0037-1106}, doi = {10.1785/0120200376}, pages = {361 -- 380}, year = {2022}, abstract = {Ground motion with strong-velocity pulses can cause significant damage to buildings and structures at certain periods; hence, knowing the period and velocity amplitude of such pulses is critical for earthquake structural engineering. However, the physical factors relating the scaling of pulse periods with magnitude are poorly understood. In this study, we investigate moderate but damaging earthquakes (M-w 6-7) and characterize ground- motion pulses using the method of Shahi and Baker (2014) while considering the potential static-offset effects. We confirm that the within-event variability of the pulses is large. The identified pulses in this study are mostly from strike-slip-like earthquakes. We further perform simulations using the freq uency-wavenumber algorithm to investigate the causes of the variability of the pulse periods within and between events for moderate strike-slip earthquakes. We test the effect of fault dips, and the impact of the asperity locations and sizes. The simulations reveal that the asperity properties have a high impact on the pulse periods and amplitudes at nearby stations. Our results emphasize the importance of asperity characteristics, in addition to earthquake magnitudes for the occurrence and properties of pulses produced by the forward directivity effect. We finally quantify and discuss within- and between-event variabilities of pulse properties at short distances.}, language = {en} } @article{WutzlerHudsonThieken2022, author = {Wutzler, Bianca and Hudson, Paul and Thieken, Annegret}, title = {Adaptation strategies of flood-damaged businesses in Germany}, series = {Frontiers in water}, volume = {4}, journal = {Frontiers in water}, publisher = {Frontiers Media}, address = {Lausanne}, issn = {2624-9375}, doi = {10.3389/frwa.2022.932061}, pages = {13}, year = {2022}, abstract = {Flood risk management in Germany follows an integrative approach in which both private households and businesses can make an important contribution to reducing flood damage by implementing property-level adaptation measures. While the flood adaptation behavior of private households has already been widely researched, comparatively less attention has been paid to the adaptation strategies of businesses. However, their ability to cope with flood risk plays an important role in the social and economic development of a flood-prone region. Therefore, using quantitative survey data, this study aims to identify different strategies and adaptation drivers of 557 businesses damaged by a riverine flood in 2013 and 104 businesses damaged by pluvial or flash floods between 2014 and 2017. Our results indicate that a low perceived self-efficacy may be an important factor that can reduce the motivation of businesses to adapt to flood risk. Furthermore, property-owners tended to act more proactively than tenants. In addition, high experience with previous flood events and low perceived response costs could strengthen proactive adaptation behavior. These findings should be considered in business-tailored risk communication.}, language = {en} } @article{WolfHuismansBraunetal.2022, author = {Wolf, Sebastian G. and Huismans, Ritske S. and Braun, Jean and Yuan, Xiaoping}, title = {Topography of mountain belts controlled by rheology and surface processes}, series = {Nature : the international weekly journal of science}, volume = {606}, journal = {Nature : the international weekly journal of science}, number = {7914}, publisher = {Nature portfolio}, address = {Berlin}, issn = {0028-0836}, doi = {10.1038/s41586-022-04700-6}, pages = {516 -- 521}, year = {2022}, abstract = {It is widely recognized that collisional mountain belt topography is generated by crustal thickening and lowered by river bedrock erosion, linking climate and tectonics(1-4). However, whether surface processes or lithospheric strength control mountain belt height, shape and longevity remains uncertain. Additionally, how to reconcile high erosion rates in some active orogens with long-term survival of mountain belts for hundreds of millions of years remains enigmatic. Here we investigate mountain belt growth and decay using a new coupled surface process(5,6) and mantle-scale tectonic model(7). End-member models and the new non-dimensional Beaumont number, Bm, quantify how surface processes and tectonics control the topographic evolution of mountain belts, and enable the definition of three end-member types of growing orogens: type 1, non-steady state, strength controlled (Bm > 0.5); type 2, flux steady state(8), strength controlled (Bm approximate to 0.4-0.5); and type 3, flux steady state, erosion controlled (Bm < 0.4). Our results indicate that tectonics dominate in Himalaya-Tibet and the Central Andes (both type 1), efficient surface processes balance high convergence rates in Taiwan (probably type 2) and surface processes dominate in the Southern Alps of New Zealand (type 3). Orogenic decay is determined by erosional efficiency and can be subdivided into two phases with variable isostatic rebound characteristics and associated timescales. The results presented here provide a unified framework explaining how surface processes and lithospheric strength control the height, shape, and longevity of mountain belts.}, language = {en} } @article{WilhelmsBoersigYangetal.2022, author = {Wilhelms, Andre and B{\"o}rsig, Nicolas and Yang, Jingwei and Holbach, Andreas and Norra, Stefan}, title = {Insights into phytoplankton dynamics and water quality monitoring with the BIOFISH at the Elbe River, Germany}, series = {Water}, volume = {14}, journal = {Water}, number = {13}, publisher = {MDPI}, address = {Basel}, issn = {2073-4441}, doi = {10.3390/w14132078}, pages = {20}, year = {2022}, abstract = {Understanding the key factors influencing the water quality of large river systems forms an important basis for the assessment and protection of cross-regional ecosystems and the implementation of adapted water management concepts. However, identifying these factors requires in-depth comprehension of the unique environmental systems, which can only be achieved by detailed water quality monitoring. Within the scope of the joint science and sports event "Elbschwimmstaffel" (swimming relay on the river Elbe) in June/July 2017 organized by the German Ministry of Education and Research, water quality data were acquired along a 550 km long stretch of the Elbe River in Germany. During the survey, eight physiochemical water quality parameters were recorded in high spatial and temporal resolution with the BIOFISH multisensor system. Multivariate statistical methods were applied to identify and delineate processes influencing the water quality. The BIOFISH dataset revealed that phytoplankton activity has a major impact on the water quality of the Elbe River in the summer months. The results suggest that phytoplankton biomass constitutes a substantial proportion of the suspended particles and that photosynthetic activity of phytoplankton is closely related to significant temporal changes in pH and oxygen saturation. An evaluation of the BIOFISH data based on the combination of statistical analysis with weather and discharge data shows that the hydrological and meteorological history of the sampled water body was the main driver of phytoplankton dynamics. This study demonstrates the capacity of longitudinal river surveys with the BIOFISH or similar systems for water quality assessment, the identification of pollution sources and their utilization for online in situ monitoring of rivers.}, language = {en} } @article{WestRosolemMacDonaldetal.2022, author = {West, Charles and Rosolem, Rafael and MacDonald, Alan M. and Cuthbert, Mark O. and Wagener, Thorsten}, title = {Understanding process controls on groundwater recharge variability across Africa through recharge landscapes}, series = {Journal of hydrology}, volume = {612}, journal = {Journal of hydrology}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0022-1694}, doi = {10.1016/j.jhydrol.2022.127967}, pages = {19}, year = {2022}, abstract = {Groundwater is critical in supporting current and future reliable water supply throughout Africa. Although continental maps of groundwater storage and recharge have been developed, we currently lack a clear understanding on how the controls on groundwater recharge vary across the entire continent. Reviewing the existing literature, we synthesize information on reported groundwater recharge controls in Africa. We find that 15 out of 22 of these controls can be characterised using global datasets. We develop 11 descriptors of climatic, topographic, vegetation, soil and geologic properties using global datasets, to characterise groundwater recharge controls in Africa. These descriptors cluster Africa into 15 Recharge Landscape Units for which we expect recharge controls to be similar. Over 80\% of the continents land area is organized by just nine of these units. We also find that aggregating the Units by similarity into four broader Recharge Landscapes (Desert, Dryland, Wet tropical and Wet tropical forest) provides a suitable level of landscape organisation to explain differences in ground-based long-term mean annual recharge and recharge ratio (annual recharge / annual precipitation) estimates. Furthermore, wetter Recharge Landscapes are more efficient in converting rainfall to recharge than drier Recharge Landscapes as well as having higher annual recharge rates. In Dryland Recharge Landscapes, we found that annual recharge rates largely varied according to mean annual precipitation, whereas recharge ratio estimates increase with increasing monthly variability in P-PET. However, we were unable to explain why ground based estimates of recharge signatures vary across other Recharge Landscapes, in which there are fewer ground based recharge estimates, using global datasets alone. Even in dryland regions, there is still considerable unexplained variability in the estimates of annual recharge and recharge ratio, stressing the limitations of global datasets for investigating ground-based information.}, language = {en} } @article{WeidleWiesenbergElSharkawyetal.2022, author = {Weidle, Christian and Wiesenberg, Lars and El-Sharkawy, Amr and Kr{\"u}ger, Frank and Scharf, Andreas and Agard, Philippe and Meier, Thomas}, title = {A 3-D crustal shear wave velocity model and Moho map below the Semail Ophiolite, eastern Arabia}, series = {Geophysical journal international}, volume = {231}, journal = {Geophysical journal international}, number = {2}, publisher = {Oxford University Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggac223}, pages = {817 -- 834}, year = {2022}, abstract = {The Semail Ophiolite in eastern Arabia is the largest and best-exposed slice of oceanic lithosphere on land. Detailed knowledge of the tectonic evolution of the shallow crust, in particular during and after ophiolite obduction in Late Cretaceous times is contrasted by few constraints on physical and compositional properties of the middle and lower continental crust below the obducted units. The role of inherited, pre-obduction crustal architecture remains therefore unaccounted for in our understanding of crustal evolution and the present-day geology. Based on seismological data acquired during a 27-month campaign in northern Oman, Ambient Seismic Noise Tomography and Receiver Function analysis provide for the first time a 3-D radially anisotropic shear wave velocity (V-S) model and a consistent Moho map below the iconic Semail Ophiolite. The model highlights deep crustal boundaries that segment the eastern Arabian basement in two distinct units. The previously undescribed Western Jabal Akhdar Zone separates Arabian crust with typical continental properties and a thickness of similar to 40-45 km in the northwest from a compositionally different terrane in the southeast that is interpreted as a terrane accreted during the Pan-African orogeny in Neoproterozoic times. East of the Ibra Zone, another deep crustal boundary, crustal thickness decreases to 30-35 km and very high lower crustal V-S suggest large-scale mafic intrusions into, and possible underplating of the Arabian continental crust that occurred most likely during Permian breakup of Pangea. Mafic reworking is sharply bounded by the (upper crustal) Semail Gap Fault Zone, northwest of which no such high velocities are found in the crust. Topography of the Oman Mountains is supported by a mild crustal root and Moho depth below the highest topography, the Jabal Akhdar Dome, is similar to 42 km. Radial anisotropy is robustly resolved in the upper crust and aids in discriminating dipping allochthonous units from autochthonous sedimentary rocks that are indistinguishable by isotropic V-S alone. Lateral thickness variations of the ophiolite highlight the Haylayn Ophiolite Massif on the northern flank of Jabal Akhdar Dome and the Hawasina Window as the deepest reaching unit. Ophiolite thickness is similar to 10 km in the southern and northern massifs, and <= 5 km elsewhere.}, language = {en} } @article{ViltresNobileVasyuraBathkeetal.2022, author = {Viltres, Renier and Nobile, Adriano and Vasyura-Bathke, Hannes and Trippanera, Daniele and Xu, Wenbin and J{\´o}nsson, Sigurj{\´o}n}, title = {Transtensional rupture within a diffuse plate boundary zone during the 2020 M-w 6.4 Puerto Rico earthquake}, series = {Seismological research letters}, volume = {93}, journal = {Seismological research letters}, number = {2A}, publisher = {Seismological Society of America}, address = {Boulder, Colo.}, issn = {0895-0695}, doi = {10.1785/0220210261}, pages = {567 -- 583}, year = {2022}, abstract = {On 7 January 2020, an M-w 6.4 earthquake occurred in the northeastern Caribbean, a few kilometers offshore of the island of Puerto Rico. It was the mainshock of a complex seismic sequence, characterized by a large number of energetic earthquakes illuminating an east-west elongated area along the southwestern coast of Puerto Rico. Deformation fields constrained by Interferometric Synthetic Aperture Radar and Global Navigation Satellite System data indicate that the coseismic movements affected only the western part of the island. To assess the mainshock's source fault parameters, we combined the geodetically derived coseismic deformation with teleseismic waveforms using Bayesian inference. The results indicate a roughly east-west oriented fault, dipping northward and accommodating similar to 1.4 m of transtensional motion. Besides, the determined location and orientation parameters suggest an offshore continuation of the recently mapped North Boqueron Bay-Punta Montalva fault in southwest Puerto Rico. This highlights the existence of unmapped faults with moderate-to-large earthquake potential within the Puerto Rico region.}, language = {en} } @article{VehLuetzowKharlamovaetal.2022, author = {Veh, Georg and L{\"u}tzow, Natalie and Kharlamova, Varvara and Petrakov, Dmitry and Hugonnet, Romain and Korup, Oliver}, title = {Trends, breaks, and biases in the frequency of reported glacier lake outburst floods}, series = {Earth's future}, volume = {10}, journal = {Earth's future}, number = {3}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2328-4277}, doi = {10.1029/2021EF002426}, pages = {14}, year = {2022}, abstract = {Thousands of glacier lakes have been forming behind natural dams in high mountains following glacier retreat since the early 20th century. Some of these lakes abruptly released pulses of water and sediment with disastrous downstream consequences. Yet it remains unclear whether the reported rise of these glacier lake outburst floods (GLOFs) has been fueled by a warming atmosphere and enhanced meltwater production, or simply a growing research effort. Here we estimate trends and biases in GLOF reporting based on the largest global catalog of 1,997 dated glacier-related floods in six major mountain ranges from 1901 to 2017. We find that the positive trend in the number of reported GLOFs has decayed distinctly after a break in the 1970s, coinciding with independently detected trend changes in annual air temperatures and in the annual number of field-based glacier surveys (a proxy of scientific reporting). We observe that GLOF reports and glacier surveys decelerated, while temperature rise accelerated in the past five decades. Enhanced warming alone can thus hardly explain the annual number of reported GLOFs, suggesting that temperature-driven glacier lake formation, growth, and failure are weakly coupled, or that outbursts have been overlooked. Indeed, our analysis emphasizes a distinct geographic and temporal bias in GLOF reporting, and we project that between two to four out of five GLOFs on average might have gone unnoticed in the early to mid-20th century. We recommend that such biases should be considered, or better corrected for, when attributing the frequency of reported GLOFs to atmospheric warming.}, language = {en} } @article{vanGeffenHeimBriegeretal.2022, author = {van Geffen, Femke and Heim, Birgit and Brieger, Frederic and Geng, Rongwei and Shevtsova, Iuliia and Schulte, Luise and Stuenzi, Simone M. and Bernhardt, Nadine and Troeva, Elena I. and Pestryakova, Luidmila Agafyevna and Zakharov, Evgenii S. and Pflug, Bringfried and Herzschuh, Ulrike and Kruse, Stefan}, title = {SiDroForest: a comprehensive forest inventory of Siberian boreal forest investigations including drone-based point clouds, individually labeled trees, synthetically generated tree crowns, and Sentinel-2 labeled image patches}, series = {Earth system science data}, volume = {14}, journal = {Earth system science data}, number = {11}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1866-3508}, doi = {10.5194/essd-14-4967-2022}, pages = {4967 -- 4994}, year = {2022}, abstract = {The SiDroForest (Siberian drone-mapped forest inventory) data collection is an attempt to remedy the scarcity of forest structure data in the circumboreal region by providing adjusted and labeled tree-level and vegetation plot-level data for machine learning and upscaling purposes. We present datasets of vegetation composition and tree and plot level forest structure for two important vegetation transition zones in Siberia, Russia; the summergreen-evergreen transition zone in Central Yakutia and the tundra-taiga transition zone in Chukotka (NE Siberia). The SiDroForest data collection consists of four datasets that contain different complementary data types that together support in-depth analyses from different perspectives of Siberian Forest plot data for multi-purpose applications. i. Dataset 1 provides unmanned aerial vehicle (UAV)-borne data products covering the vegetation plots surveyed during fieldwork (Kruse et al., 2021, ). The dataset includes structure-from-motion (SfM) point clouds and red-green-blue (RGB) and red-green-near-infrared (RGN) orthomosaics. From the orthomosaics, point-cloud products were created such as the digital elevation model (DEM), canopy height model (CHM), digital surface model (DSM) and the digital terrain model (DTM). The point-cloud products provide information on the three-dimensional (3D) structure of the forest at each plot. Dataset 2 contains spatial data in the form of point and polygon shapefiles of 872 individually labeled trees and shrubs that were recorded during fieldwork at the same vegetation plots (van Geffen et al., 2021c, ). The dataset contains information on tree height, crown diameter, and species type. These tree and shrub individually labeled point and polygon shapefiles were generated on top of the RGB UVA orthoimages. The individual tree information collected during the expedition such as tree height, crown diameter, and vitality are provided in table format. This dataset can be used to link individual information on trees to the location of the specific tree in the SfM point clouds, providing for example, opportunity to validate the extracted tree height from the first dataset. The dataset provides unique insights into the current state of individual trees and shrubs and allows for monitoring the effects of climate change on these individuals in the future. Dataset 3 contains a synthesis of 10 000 generated images and masks that have the tree crowns of two species of larch ( and ) automatically extracted from the RGB UAV images in the common objects in context (COCO) format (van Geffen et al., 2021a, ). As machine-learning algorithms need a large dataset to train on, the synthetic dataset was specifically created to be used for machine-learning algorithms to detect Siberian larch species. Larix gmeliniiLarix cajanderiDataset 4 contains Sentinel-2 (S-2) Level-2 bottom-of-atmosphere processed labeled image patches with seasonal information and annotated vegetation categories covering the vegetation plots (van Geffen et al., 2021b, ). The dataset is created with the aim of providing a small ready-to-use validation and training dataset to be used in various vegetation-related machine-learning tasks. It enhances the data collection as it allows classification of a larger area with the provided vegetation classes. The SiDroForest data collection serves a variety of user communities.
The detailed vegetation cover and structure information in the first two datasets are of use for ecological applications, on one hand for summergreen and evergreen needle-leaf forests and also for tundra-taiga ecotones. Datasets 1 and 2 further support the generation and validation of land cover remote-sensing products in radar and optical remote sensing. In addition to providing information on forest structure and vegetation composition of the vegetation plots, the third and fourth datasets are prepared as training and validation data for machine-learning purposes. For example, the synthetic tree-crown dataset is generated from the raw UAV images and optimized to be used in neural networks. Furthermore, the fourth SiDroForest dataset contains S-2 labeled image patches processed to a high standard that provide training data on vegetation class categories for machine-learning classification with JavaScript Object Notation (JSON) labels provided. The SiDroForest data collection adds unique insights into remote hard-to-reach circumboreal forest regions.}, language = {en} } @article{ValenzuelaMalebranCescaLopezCominoetal.2022, author = {Valenzuela-Malebran, Carla and Cesca, Simone and Lopez-Comino, Jos{\´e} {\´A}ngel and Zeckra, Martin and Kr{\"u}ger, F. and Dahm, Torsten}, title = {Source mechanisms and rupture processes of the Jujuy seismic nest, Chile-Argentina border}, series = {Journal of South American earth sciences}, volume = {117}, journal = {Journal of South American earth sciences}, publisher = {Elsevier}, address = {Oxford}, issn = {0895-9811}, doi = {10.1016/j.jsames.2022.103887}, pages = {13}, year = {2022}, abstract = {The Altiplano-Puna plateau, in Central Andes, is the second-largest continental plateau on Earth, extending between 22 degrees and 27 degrees S at an average altitude of 4400 m. The Puna plateau has been formed in consequence of the subduction of the oceanic Nazca Plate beneath the continental South American plate, which has an average crustal thickness of 50 km at this location. A large seismicity cluster, the Jujuy cluster, is observed at depth of 150-250 km beneath the central region of the Puna plateau. The cluster is seismically very active, with hundreds of earthquakes reported and a peak magnitude MW 6.6 on 25th August 2006. The cluster is situated in one of three band of intermediate-depth focus seismicity, which extend parallel to the trench roughly North to South. It has been hypothesized that the Jujuy cluster could be a seismic nest, a compact seismogenic region characterized by a high stationary activity relative to its surroundings. In this study, we collected more than 40 years of data from different catalogs and proof that the cluster meets the three conditions of a seismic nest. Compared to other known intermediate depth nests at Hindu Kush (Afganisthan) or Bucaramanga (Colombia), the Jujuy nest presents an outstanding seismicity rate, with more than 100 M4+ earthquakes per year. We additionally performed a detailed analysis of the rupture process of some of the largest earthquakes in the nest, by means of moment tensor inversion and directivity analysis. We focused on the time period 2017-2018, where the seismic monitoring was the most extended. Our results show that earthquakes in the nest take place within the eastward subducting oceanic plate, but rupture along sub-horizontal planes dipping westward. We suggest that seismicity at Jujuy nest is controlled by dehydration processes, which are also responsible for the generation of fluids ascending to the crust beneath the Puna volcanic region. We use the rupture plane and nest geometry to provide a constraint to maximal expected magnitude, which we estimate as MW -6.7.}, language = {en} } @article{TuerkerCottonPilzetal.2022, author = {T{\"u}rker, Elif and Cotton, Fabrice and Pilz, Marco and Weatherill, Graeme}, title = {Analysis of the 2019 Mw 5.8 Silivri earthquake ground motions}, series = {Seismological research letters}, volume = {93}, journal = {Seismological research letters}, number = {2A}, publisher = {Seismological Society of America}, address = {Boulder, Colo.}, issn = {0895-0695}, doi = {10.1785/0220210168}, pages = {693 -- 705}, year = {2022}, abstract = {The main Marmara fault (MMF) extends for 150 km through the Sea of Marmara and forms the only portion of the North Anatolian fault zone that has not ruptured in a large event (Mw >7) for the last 250 yr. Accordingly, this portion is potentially a major source contributing to the seismic hazard of the Istanbul region. On 26 September 2019, a sequence of moderate-sized events started along the MMF only 20 km south of Istanbul and were widely felt by the population. The largest three events, 26 September Mw 5.8 (10:59 UTC), 26 September 2019 Mw 4.1 (11:26 UTC), and 20 January 2020 Mw 4.7 were recorded by numerous strong-motion seismic stations and the resulting ground motions were compared to the predicted means resulting from a set of the most recent ground-motion prediction equations (GMPEs). The estimated residuals were used to investigate the spatial variation of ground motion across the Marmara region. Our results show a strong azimuthal trend in ground-motion residuals, which might indicate systematically repeating directivity effects toward the eastern Marmara region.}, language = {en} } @article{TrauthMarwan2022, author = {Trauth, Martin H. and Marwan, Norbert}, title = {Introduction-time series analysis for Earth, climate and life interactions}, series = {Quaternary science reviews : the international multidisciplinary research and review journal}, volume = {284}, journal = {Quaternary science reviews : the international multidisciplinary research and review journal}, publisher = {Elsevier}, address = {Amsterdam [u.a.]}, issn = {0277-3791}, doi = {10.1016/j.quascirev.2022.107475}, pages = {3}, year = {2022}, language = {en} } @book{Trauth2022, author = {Trauth, Martin H.}, title = {Python Recipes for Earth Sciences}, series = {Springer Textbooks in Earth Sciences, Geography and Environment}, journal = {Springer Textbooks in Earth Sciences, Geography and Environment}, publisher = {Springer}, address = {Cham}, isbn = {978-3-031-07719-7}, issn = {2510-1307}, doi = {10.1007/978-3-031-07719-7}, pages = {453}, year = {2022}, abstract = {Python is used in a wide range of geoscientific applications, such as in processing images for remote sensing, in generating and processing digital elevation models, and in analyzing time series. This book introduces methods of data analysis in the geosciences using Python that include basic statistics for univariate, bivariate, and multivariate data sets, time series analysis, and signal processing; the analysis of spatial and directional data; and image analysis. The text includes numerous examples that demonstrate how Python can be used on data sets from the earth sciences. The supplementary electronic material (available online through Springer Link) contains the example data as well as recipes that include all the Python commands featured in the book.}, language = {en} } @phdthesis{Tranter2022, author = {Tranter, Morgan Alan}, title = {Numerical quantification of barite reservoir scaling and the resulting injectivity loss in geothermal systems}, doi = {10.25932/publishup-56113}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-561139}, school = {Universit{\"a}t Potsdam}, pages = {131}, year = {2022}, abstract = {Due to the major role of greenhouse gas emissions in global climate change, the development of non-fossil energy technologies is essential. Deep geothermal energy represents such an alternative, which offers promising properties such as a high base load capability and a large untapped potential. The present work addresses barite precipitation within geothermal systems and the associated reduction in rock permeability, which is a major obstacle to maintaining high efficiency. In this context, hydro-geochemical models are essential to quantify and predict the effects of precipitation on the efficiency of a system. The objective of the present work is to quantify the induced injectivity loss using numerical and analytical reactive transport simulations. For the calculations, the fractured-porous reservoirs of the German geothermal regions North German Basin (NGB) and Upper Rhine Graben (URG) are considered. Similar depth-dependent precipitation potentials could be determined for both investigated regions (2.8-20.2 g/m3 fluid). However, the reservoir simulations indicate that the injectivity loss due to barite deposition in the NGB is significant (1.8\%-6.4\% per year) and the longevity of the system is affected as a result; this is especially true for deeper reservoirs (3000 m). In contrast, simulations of URG sites indicate a minor role of barite (< 0.1\%-1.2\% injectivity loss per year). The key differences between the investigated regions are reservoir thicknesses and the presence of fractures in the rock, as well as the ionic strength of the fluids. The URG generally has fractured-porous reservoirs with much higher thicknesses, resulting in a greater distribution of precipitates in the subsurface. Furthermore, ionic strengths are higher in the NGB, which accelerates barite precipitation, causing it to occur more concentrated around the wellbore. The more concentrated the precipitates occur around the wellbore, the higher the injectivity loss. In this work, a workflow was developed within which numerical and analytical models can be used to estimate and quantify the risk of barite precipitation within the reservoir of geothermal systems. A key element is a newly developed analytical scaling score that provides a reliable estimate of induced injectivity loss. The key advantage of the presented approach compared to fully coupled reservoir simulations is its simplicity, which makes it more accessible to plant operators and decision makers. Thus, in particular, the scaling score can find wide application within geothermal energy, e.g., in the search for potential plant sites and the estimation of long-term efficiency.}, language = {en} } @article{TofeldeBufeTurowski2022, author = {Tofelde, Stefanie and Bufe, Aaron and Turowski, Jens M.}, title = {Hillslope Sediment Supply Limits Alluvial Valley Width}, series = {AGU Advances}, journal = {AGU Advances}, publisher = {American Geophysical Union (AGU); Wiley}, address = {Hoboken, New Jersey, USA}, issn = {2576-604X}, doi = {10.1029/2021AV000641}, pages = {20}, year = {2022}, abstract = {River-valley morphology preserves information on tectonic and climatic conditions that shape landscapes. Observations suggest that river discharge and valley-wall lithology are the main controls on valley width. Yet, current models based on these observations fail to explain the full range of cross-sectional valley shapes in nature, suggesting hitherto unquantified controls on valley width. In particular, current models cannot explain the existence of paired terrace sequences that form under cyclic climate forcing. Paired river terraces are staircases of abandoned floodplains on both valley sides, and hence preserve past valley widths. Their formation requires alternating phases of predominantly river incision and predominantly lateral planation, plus progressive valley narrowing. While cyclic Quaternary climate changes can explain shifts between incision and lateral erosion, the driving mechanism of valley narrowing is unknown. Here, we extract valley geometries from climatically formed, alluvial river-terrace sequences and show that across our dataset, the total cumulative terrace height (here: total valley height) explains 90\%-99\% of the variance in valley width at the terrace sites. This finding suggests that valley height, or a parameter that scales linearly with valley height, controls valley width in addition to river discharge and lithology. To explain this valley-width-height relationship, we reformulate existing valley-width models and suggest that, when adjusting to new boundary conditions, alluvial valleys evolve to a width at which sediment removal from valley walls matches lateral sediment supply from hillslope erosion. Such a hillslope-channel coupling is not captured in current valley-evolution models. Our model can explain the existence of paired terrace sequences under cyclic climate forcing and relates valley width to measurable field parameters. Therefore, it facilitates the reconstruction of past climatic and tectonic conditions from valley topography.}, language = {en} } @misc{TofeldeBufeTurowski2022, author = {Tofelde, Stefanie and Bufe, Aaron and Turowski, Jens M.}, title = {Hillslope Sediment Supply Limits Alluvial Valley Width}, series = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {1289}, issn = {1866-8372}, doi = {10.25932/publishup-57287}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-572879}, pages = {20}, year = {2022}, abstract = {River-valley morphology preserves information on tectonic and climatic conditions that shape landscapes. Observations suggest that river discharge and valley-wall lithology are the main controls on valley width. Yet, current models based on these observations fail to explain the full range of cross-sectional valley shapes in nature, suggesting hitherto unquantified controls on valley width. In particular, current models cannot explain the existence of paired terrace sequences that form under cyclic climate forcing. Paired river terraces are staircases of abandoned floodplains on both valley sides, and hence preserve past valley widths. Their formation requires alternating phases of predominantly river incision and predominantly lateral planation, plus progressive valley narrowing. While cyclic Quaternary climate changes can explain shifts between incision and lateral erosion, the driving mechanism of valley narrowing is unknown. Here, we extract valley geometries from climatically formed, alluvial river-terrace sequences and show that across our dataset, the total cumulative terrace height (here: total valley height) explains 90\%-99\% of the variance in valley width at the terrace sites. This finding suggests that valley height, or a parameter that scales linearly with valley height, controls valley width in addition to river discharge and lithology. To explain this valley-width-height relationship, we reformulate existing valley-width models and suggest that, when adjusting to new boundary conditions, alluvial valleys evolve to a width at which sediment removal from valley walls matches lateral sediment supply from hillslope erosion. Such a hillslope-channel coupling is not captured in current valley-evolution models. Our model can explain the existence of paired terrace sequences under cyclic climate forcing and relates valley width to measurable field parameters. Therefore, it facilitates the reconstruction of past climatic and tectonic conditions from valley topography.}, language = {en} } @article{TellaWinterleitnerMutti2022, author = {Tella, Timothy Oluwatobi and Winterleitner, Gerd and Mutti, Maria}, title = {Investigating the role of differential biotic production on carbonate geometries through stratigraphic forward modelling and sensitivity analysis}, series = {Petroleum geoscience}, volume = {28}, journal = {Petroleum geoscience}, number = {2}, publisher = {Geological Soc. Publ. House}, address = {Bath}, issn = {1354-0793}, doi = {10.1144/petgeo2021-053}, pages = {20}, year = {2022}, abstract = {The geometry of carbonate platforms reflects the interaction of several factors. However, the impact of carbonate-producing organisms has been poorly investigated so far. This study applies stratigraphic forward modelling (SFM) and sensitivity analysis to examine, referenced to the Miocene Llucmajor Platform, the effect of changes of dominant biotic production in the oligophotic and euphotic zones on platform geometry. Our results show that the complex interplay of carbonate production rates, bathymetry and variations in accommodation space control the platform geometry. The main driver of progradation is the oligophotic production of rhodalgal sediments during the lowstands. This study demonstrates that platform geometry and internal architecture varies significantly according to the interaction of the predominant carbonate-producing biotas. The input parameters for this study are based on well-understood Miocene carbonate biotas with characteristic euphotic, oligophotic and photo-independent carbonate production in which it is crucial that each carbonate-producing class is modelled explicitly within the simulation run and not averaged with a single carbonate production-depth profile. This is important in subsurface exploration studies based on stratigraphic forward models where the overall platform geometry may be approximated through calibration runs, and constrained by seismic surveys and wellbores. However, the internal architecture is likely to be oversimplified without an in-depth understanding of the target carbonate system and a transfer to forward modelling parameters.}, language = {en} } @article{TawfikOndrakWinterleitneretal.2022, author = {Tawfik, Ahmed Y. and Ondrak, Robert and Winterleitner, Gerd and Mutti, Maria}, title = {Source rock evaluation and petroleum system modeling of the East Beni Suef Basin, north Eastern Desert, Egypt}, series = {Journal of African earth sciences}, volume = {193}, journal = {Journal of African earth sciences}, publisher = {Elsevier}, address = {Oxford}, issn = {1464-343X}, doi = {10.1016/j.jafrearsci.2022.104575}, pages = {21}, year = {2022}, abstract = {This study deals with the East Beni Suef Basin (Eastern Desert, Egypt) and aims to evaluate the source-generative potential, reconstruct the burial and thermal history, examine the most influential parameters on thermal maturity modeling, and improve on the models already published for the West Beni Suef to ultimately formulate a complete picture of the whole basin evolution. Source rock evaluation was carried out based on TOC, Rock-Eval pyrolysis, and visual kerogen petrography analyses. Three kerogen types (II, II/III, and III) are distinguished in the East Beni Suef Basin, where the Abu Roash "F" Member acts as the main source rock with good to excellent source potential, oil-prone mainly type II kerogen, and immature to marginal maturity levels. The burial history shows four depositional and erosional phases linked with the tectonic evolution of the basin. A hiatus (due to erosion or non-deposition) has occurred during the Late Eocene-Oligocene in the East Beni Suef Basin, while the West Beni Suef Basin has continued subsiding. Sedimentation began later (Middle to Late Albian) with lower rates in the East Beni Suef Basin compared with the West Beni Suef Basin (Early Albian). The Abu Roash "F" source rock exists in the early oil window with a present-day transformation ratio of about 19\% and 21\% in the East and West Beni Suef Basin, respectively, while the Lower Kharita source rock, which is only recorded in the West Beni Suef Basin, has reached the late oil window with a present-day transformation ratio of about 70\%. The magnitude of erosion and heat flow have proportional and mutual effects on thermal maturity. We present three possible scenarios of basin modeling in the East Beni Suef Basin concerning the erosion from the Apollonia and Dabaa formations. Results of this work can serve as a basis for subsequent 2D and/or 3D basin modeling, which are highly recommended to further investigate the petroleum system evolution of the Beni Suef Basin.}, language = {en} } @article{SvennevigHermannsKeidingetal.2022, author = {Svennevig, Kristian and Hermanns, Reginald L. and Keiding, Marie and Binder, Daniel and Citterio, Michele and Dahl-Jensen, Trine and Mertl, Stefan and S{\o}rensen, Erik Vest and Voss, Peter Henrik}, title = {A large frozen debris avalanche entraining warming permafrost ground-the June 2021 Assapaat landslibe, West Greenland}, series = {Landslides}, volume = {19}, journal = {Landslides}, publisher = {Springer}, address = {Heidelberg}, issn = {1612-510X}, doi = {10.1007/s10346-022-01922-7}, pages = {2549 -- 2567}, year = {2022}, abstract = {A large landslide (frozen debris avalanche) occurred at Assapaat on the south coast of the Nuussuaq Peninsula in Central West Greenland on June 13, 2021, at 04:04 local time. We present a compilation of available data from field observations, photos, remote sensing, and seismic monitoring to describe the event. Analysis of these data in combination with an analysis of pre- and post-failure digital elevation models results in the first description of this type of landslide. The frozen debris avalanche initiated as a 6.9 * 10(6) m(3) failure of permafrozen talus slope and underlying colluvium and till at 600-880 m elevation. It entrained a large volume of permafrozen colluvium along its 2.4 km path in two subsequent entrainment phases accumulating a total volume between 18.3 * 10(6) and 25.9 * 10(6) m(3). About 3.9 * 10(6) m(3) is estimated to have entered the Vaigat strait; however, no tsunami was reported, or is evident in the field. This is probably because the second stage of entrainment along with a flattening of slope angle reduced the mobility of the frozen debris avalanche. We hypothesise that the initial talus slope failure is dynamically conditioned by warming of the ice matrix that binds the permafrozen talus slope. When the slope ice temperature rises to a critical level, its shear resistance is reduced, resulting in an unstable talus slope prone to failure. Likewise, we attribute the large-scale entrainment to increasing slope temperature and take the frozen debris avalanche as a strong sign that the permafrost in this region is increasingly at a critical state. Global warming is enhanced in the Arctic and frequent landslide events in the past decade in Western Greenland let us hypothesise that continued warming will lead to an increase in the frequency and magnitude of these types of landslides. Essential data for critical arctic slopes such as precipitation, snowmelt, and ground and surface temperature are still missing to further test this hypothesis. It is thus strongly required that research funds are made available to better predict the change of landslide threat in the Arctic.}, language = {en} } @article{SudibyoEiblHainzletal.2022, author = {Sudibyo, Maria R. P. and Eibl, Eva P. S. and Hainzl, Sebastian and Hersir, Gylfi P{\´a}ll}, title = {Eruption Forecasting of Strokkur Geyser, Iceland, Using Permutation Entropy}, series = {Journal of geophysical research : Solid earth}, volume = {127}, journal = {Journal of geophysical research : Solid earth}, number = {10}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9313}, doi = {10.1029/2022JB024840}, pages = {15}, year = {2022}, abstract = {A volcanic eruption is usually preceded by seismic precursors, but their interpretation and use for forecasting the eruption onset time remain a challenge. A part of the eruptive processes in open conduits of volcanoes may be similar to those encountered in geysers. Since geysers erupt more often, they are useful sites for testing new forecasting methods. We tested the application of Permutation Entropy (PE) as a robust method to assess the complexity in seismic recordings of the Strokkur geyser, Iceland. Strokkur features several minute-long eruptive cycles, enabling us to verify in 63 recorded cycles whether PE behaves consistently from one eruption to the next one. We performed synthetic tests to understand the effect of different parameter settings in the PE calculation. Our application to Strokkur shows a distinct, repeating PE pattern consistent with previously identified phases in the eruptive cycle. We find a systematic increase in PE within the last 15 s before the eruption, indicating that an eruption will occur. We quantified the predictive power of PE, showing that PE performs better than seismic signal strength or quiescence when it comes to forecasting eruptions.}, language = {en} } @article{StuenziKruseBoikeetal.2022, author = {Stuenzi, Simone Maria and Kruse, Stefan and Boike, Julia and Herzschuh, Ulrike and Oehme, Alexander and Pestryakova, Luidmila A. and Westermann, Sebastian and Langer, Moritz}, title = {Thermohydrological impact of forest disturbances on ecosystem-protected permafrost}, series = {Journal of geophysical research : Biogeosciences}, volume = {127}, journal = {Journal of geophysical research : Biogeosciences}, number = {5}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-8953}, doi = {10.1029/2021JG006630}, pages = {24}, year = {2022}, abstract = {Boreal forests cover over half of the global permafrost area and protect underlying permafrost. Boreal forest development, therefore, has an impact on permafrost evolution, especially under a warming climate. Forest disturbances and changing climate conditions cause vegetation shifts and potentially destabilize the carbon stored within the vegetation and permafrost. Disturbed permafrost-forest ecosystems can develop into a dry or swampy bush- or grasslands, shift toward broadleaf- or evergreen needleleaf-dominated forests, or recover to the pre-disturbance state. An increase in the number and intensity of fires, as well as intensified logging activities, could lead to a partial or complete ecosystem and permafrost degradation. We study the impact of forest disturbances (logging, surface, and canopy fires) on the thermal and hydrological permafrost conditions and ecosystem resilience. We use a dynamic multilayer canopy-permafrost model to simulate different scenarios at a study site in eastern Siberia. We implement expected mortality, defoliation, and ground surface changes and analyze the interplay between forest recovery and permafrost. We find that forest loss induces soil drying of up to 44\%, leading to lower active layer thicknesses and abrupt or steady decline of a larch forest, depending on disturbance intensity. Only after surface fires, the most common disturbances, inducing low mortality rates, forests can recover and overpass pre-disturbance leaf area index values. We find that the trajectory of larch forests after surface fires is dependent on the precipitation conditions in the years after the disturbance. Dryer years can drastically change the direction of the larch forest development within the studied period.}, language = {en} } @article{StoltnowLuedersGraafetal.2022, author = {Stoltnow, Malte and L{\"u}ders, Volker and Graaf, Stefan de and Niedermann, Samuel}, title = {A geochemical study of the Sweet Home mine, Colorado Mineral Belt, USA}, series = {Mineralium deposita : international journal for geology, mineralogy and geochemistry of mineral deposits}, volume = {57}, journal = {Mineralium deposita : international journal for geology, mineralogy and geochemistry of mineral deposits}, number = {5}, publisher = {Springer}, address = {Berlin ; Heidelberg}, issn = {0026-4598}, doi = {10.1007/s00126-022-01102-6}, pages = {801 -- 825}, year = {2022}, abstract = {Deep hydrothermal Mo, W, and base metal mineralization at the Sweet Home mine (Detroit City portal) formed in response to magmatic activity during the Oligocene. Microthermometric data of fluid inclusions trapped in greisen quartz and fluorite suggest that the early-stage mineralization at the Sweet Home mine precipitated from low- to medium-salinity (1.5-11.5 wt\% equiv. NaCl), CO2-bearing fluids at temperatures between 360 and 415 degrees C and at depths of at least 3.5 km. Stable isotope and noble gas isotope data indicate that greisen formation and base metal mineralization at the Sweet Home mine was related to fluids of different origins. Early magmatic fluids were the principal source for mantle-derived volatiles (CO2, H2S/SO2, noble gases), which subsequently mixed with significant amounts of heated meteoric water. Mixing of magmatic fluids with meteoric water is constrained by delta H-2(w)-delta O-18(w) relationships of fluid inclusions. The deep hydrothermal mineralization at the Sweet Home mine shows features similar to deep hydrothermal vein mineralization at Climax-type Mo deposits or on their periphery. This suggests that fluid migration and the deposition of ore and gangue minerals in the Sweet Home mine was triggered by a deep-seated magmatic intrusion. The findings of this study are in good agreement with the results of previous fluid inclusion studies of the mineralization of the Sweet Home mine and from Climax-type Mo porphyry deposits in the Colorado Mineral Belt.}, language = {en} } @article{StolpmannMollenhauerMorgensternetal.2022, author = {Stolpmann, Lydia and Mollenhauer, Gesine and Morgenstern, Anne and Hammes, Jens S. and Boike, Julia and Overduin, Pier Paul and Grosse, Guido}, title = {Origin and pathways of dissolved organic carbon in a small catchment in the Lena River Delta}, series = {Frontiers in Earth Science}, volume = {9}, journal = {Frontiers in Earth Science}, publisher = {Frontiers Media}, address = {Lausanne}, issn = {2296-6463}, doi = {10.3389/feart.2021.759085}, pages = {15}, year = {2022}, abstract = {The Arctic is rich in aquatic systems and experiences rapid warming due to climate change. The accelerated warming causes permafrost thaw and the mobilization of organic carbon. When dissolved organic carbon is mobilized, this DOC can be transported to aquatic systems and degraded in the water bodies and further downstream. Here, we analyze the influence of different landscape components on DOC concentrations and export in a small (6.45 km(2)) stream catchment in the Lena River Delta. The catchment includes lakes and ponds, with the flow path from Pleistocene yedoma deposits across Holocene non-yedoma deposits to the river outlet. In addition to DOC concentrations, we use radiocarbon dating of DOC as well as stable oxygen and hydrogen isotopes (delta O-18 and delta D) to assess the origin of DOC. We find significantly higher DOC concentrations in the Pleistocene yedoma area of the catchment compared to the Holocene non-yedoma area with medians of 5 and 4.5 mg L-1 (p < 0.05), respectively. When yedoma thaw streams with high DOC concentration reach a large yedoma thermokarst lake, we observe an abrupt decrease in DOC concentration, which we attribute to dilution and lake processes such as mineralization. The DOC ages in the large thermokarst lake (between 3,428 and 3,637 C-14 y BP) can be attributed to a mixing of mobilized old yedoma and Holocene carbon. Further downstream after the large thermokarst lake, we find progressively younger DOC ages in the stream water to its mouth, paired with decreasing DOC concentrations. This process could result from dilution with leaching water from Holocene deposits and/or emission of ancient yedoma carbon to the atmosphere. Our study shows that thermokarst lakes and ponds may act as DOC filters, predominantly by diluting incoming waters of higher DOC concentrations or by re-mineralizing DOC to CO2 and CH4. Nevertheless, our results also confirm that the small catchment still contributes DOC on the order of 1.2 kg km(-2) per day from a permafrost landscape with ice-rich yedoma deposits to the Lena River.}, language = {en} } @article{StevanatoBaroniOswaldetal.2022, author = {Stevanato, Luca and Baroni, Gabriele and Oswald, Sascha and Lunardon, Marcello and Mareš, Vratislav and Marinello, Francesco and Moretto, Sandra and Polo, Matteo and Sartori, Paolo and Schattan, Paul and R{\"u}hm, Werner}, title = {An alternative incoming correction for cosmic-ray neutron sensing observations using local muon measurement}, series = {Geophysical research letters}, volume = {49}, journal = {Geophysical research letters}, number = {6}, publisher = {American Geophysical Union}, address = {Washington}, issn = {0094-8276}, doi = {10.1029/2021GL095383}, pages = {9}, year = {2022}, abstract = {Measuring the variability of incoming neutrons locally would be usefull for the cosmic-ray neutron sensing (CRNS) method. As the measurement of high energy neutrons is not so easy, alternative particles can be considered for such purpose. Among them, muons are particles created from the same cascade of primary cosmic-ray fluxes that generate neutrons at the ground. In addition, they can be easily detected by small and relatively inexpensive detectors. For these reasons they could provide a suitable local alternative to incoming corrections based on remote neutron monitor data. The reported measurements demonstrated that muon detection system can detect incoming cosmic-ray variations locally. Furthermore the precision of this measurement technique is considered adequate for many CRNS applications.}, language = {en} } @article{SteirouGerlitzSunetal.2022, author = {Steirou, Eva and Gerlitz, Lars and Sun, Xun and Apel, Heiko and Agarwal, Ankit and Totz, Sonja Juliana and Merz, Bruno}, title = {Towards seasonal forecasting of flood probabilities in Europe using climate and catchment information}, series = {Scientific reports}, volume = {12}, journal = {Scientific reports}, number = {1}, publisher = {Nature portfolio}, address = {Berlin}, issn = {2045-2322}, doi = {10.1038/s41598-022-16633-1}, pages = {10}, year = {2022}, abstract = {We investigate whether the distribution of maximum seasonal streamflow is significantly affected by catchment or climate state of the season/month ahead. We fit the Generalized Extreme Value (GEV) distribution to extreme seasonal streamflow for around 600 stations across Europe by conditioning the GEV location and scale parameters on 14 indices, which represent the season-ahead climate or catchment state. The comparison of these climate-informed models with the classical GEV distribution, with time-constant parameters, suggests that there is a substantial potential for seasonal forecasting of flood probabilities. The potential varies between seasons and regions. Overall, the season-ahead catchment wetness shows the highest potential, although climate indices based on large-scale atmospheric circulation, sea surface temperature or sea ice concentration also show some skill for certain regions and seasons. Spatially coherent patterns and a substantial fraction of climate-informed models are promising signs towards early alerts to increase flood preparedness already a season ahead.}, language = {en} } @phdthesis{Steding2022, author = {Steding, Svenja}, title = {Geochemical and Hydraulic Modeling of Cavernous Structures in Potash Seams}, doi = {10.25932/publishup-54818}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-548182}, school = {Universit{\"a}t Potsdam}, pages = {IX, 104}, year = {2022}, abstract = {Salt deposits offer a variety of usage types. These include the mining of rock salt and potash salt as important raw materials, the storage of energy in man-made underground caverns, and the disposal of hazardous substances in former mines. The most serious risk with any of these usage types comes from the contact with groundwater or surface water. It causes an uncontrolled dissolution of salt rock, which in the worst case can result in the flooding or collapse of underground facilities. Especially along potash seams, cavernous structures can spread quickly, because potash salts show a much higher solubility than rock salt. However, as their chemical behavior is quite complex, previous models do not account for these highly soluble interlayers. Therefore, the objective of the present thesis is to describe the evolution of cavernous structures along potash seams in space and time in order to improve hazard mitigation during the utilization of salt deposits. The formation of cavernous structures represents an interplay of chemical and hydraulic processes. Hence, the first step is to systematically investigate the dissolution and precipitation reactions that occur when water and potash salt come into contact. For this purpose, a geochemical reaction model is used. The results show that the minerals are only partially dissolved, resulting in a porous sponge like structure. With the saturation of the solution increasing, various secondary minerals are formed, whose number and type depend on the original rock composition. Field data confirm a correlation between the degree of saturation and the distance from the center of the cavern, where solution is entering. Subsequently, the reaction model is coupled with a flow and transport code and supplemented by a novel approach called 'interchange'. The latter enables the exchange of solution and rock between areas of different porosity and mineralogy, and thus ultimately the growth of the cavernous structure. By means of several scenario analyses, cavern shape, growth rate and mineralogy are systematically investigated, taking also heterogeneous potash seams into account. The results show that basically four different cases can be distinguished, with mixed forms being a frequent occurrence in nature. The classification scheme is based on the dimensionless numbers P{\´e}clet and Damk{\"o}hler, and allows for a first assessment of the hazard potential. In future, the model can be applied to any field case, using measurement data for calibration. The presented research work provides a reactive transport model that is able to spatially and temporally characterize the propagation of cavernous structures along potash seams for the first time. Furthermore, it allows to determine thickness and composition of transition zones between cavern center and unaffected salt rock. The latter is particularly important in potash mining, so that natural cavernous structures can be located at an early stage and the risk of mine flooding can thus be reduced. The models may also contribute to an improved hazard prevention in the construction of storage caverns and the disposal of hazardous waste in salt deposits. Predictions regarding the characteristics and evolution of cavernous structures enable a better assessment of potential hazards, such as integrity or stability loss, as well as of suitable mitigation measures.}, language = {en} } @article{SpallanzaniKogaCichyetal.2022, author = {Spallanzani, Roberta and Koga, Kenneth T. and Cichy, Sarah B. and Wiedenbeck, Michael and Schmidt, Burkhard C. and Oelze, Marcus and Wilke, Max}, title = {Lithium and boron diffusivity and isotopic fractionation in hydrated rhyolitic melts}, series = {Contributions to mineralogy and petrology}, volume = {177}, journal = {Contributions to mineralogy and petrology}, number = {8}, publisher = {Springer}, address = {New York}, issn = {0010-7999}, doi = {10.1007/s00410-022-01937-2}, pages = {17}, year = {2022}, abstract = {Lithium and boron are trace components of magmas, released during exsolution of a gas phase during volcanic activity. In this study, we determine the diffusivity and isotopic fractionation of Li and B in hydrous silicate melts. Two glasses were synthesized with the same rhyolitic composition (4.2 wt\% water), having different Li and B contents; these were studied in diffusion-couple experiments that were performed using an internally heated pressure vessel, operated at 300 MPa in the temperature range 700-1250 degrees C for durations from 0 s to 24 h. From this we determined activation energies for Li and B diffusion of 57 +/- 4 kJ/mol and 152 +/- 15 kJ/mol with pre-exponential factors of 1.53 x 10(-7) m(2)/s and 3.80 x 10(-8) m(2)/s, respectively. Lithium isotopic fractionation during diffusion gave beta values between 0.15 and 0.20, whereas B showed no clear isotopic fractionation. Our Li diffusivities and isotopic fractionation results differ somewhat from earlier published values, but overall confirm that Li diffusivity increases with water content. Our results on B diffusion show that similarly to Li, B mobility increases in the presence of water. By applying the Eyring relation, we confirm that B diffusivity is limited by viscous flow in silicate melts. Our results on Li and B diffusion present a new tool for understanding degassing-related processes, offering a potential geospeedometer to measure volcanic ascent rates.}, language = {en} } @phdthesis{Spallanzani2022, author = {Spallanzani, Roberta}, title = {Li and B in ascending magmas: an experimental study on their mobility and isotopic fractionation}, doi = {10.25932/publishup-56061}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-560619}, school = {Universit{\"a}t Potsdam}, pages = {xiv, 131}, year = {2022}, abstract = {This research study focuses on the behaviour of Li and B during magmatic ascent, and decompression-driven degassing related to volcanic systems. The main objective of this dissertation is to determine whether it is possible to use the diffusion properties of the two trace elements as a tool to trace magmatic ascent rate. With this objective, diffusion-couple and decompression experiments have been performed in order to study Li and B mobility in intra-melt conditions first, and then in an evolving system during decompression-driven degassing. Synthetic glasses were prepared with rhyolitic composition and an initial water content of 4.2 wt\%, and all the experiments were performed using an internally heated pressure vessel, in order to ensure a precise control on the experimental parameters such as temperature and pressure. Diffusion-couple experiments were performed with a fix pressure 300 MPa. The temperature was varied in the range of 700-1250 °C with durations between 0 seconds and 24 hours. The diffusion-couple results show that Li diffusivity is very fast and starts already at very low temperature. Significant isotopic fractionation occurs due to the faster mobility of 6Li compared to 7Li. Boron diffusion is also accelerated by the presence of water, but the results of the isotopic ratios are unclear, and further investigation would be necessary to well constrain the isotopic fractionation process of boron in hydrous silicate melts. The isotopic ratios results show that boron isotopic fractionation might be affected by the speciation of boron in the silicate melt structure, as 10B and 11B tend to have tetrahedral and trigonal coordination, respectively. Several decompression experiments were performed at 900 °C and 1000 °C, with pressures going from 300 MPa to 71-77 MPa and durations of 30 minutes, two, five and ten hours, in order to trigger water exsolution and the formation of vesicles in the sample. Textural observations and the calculation of the bubble number density confirmed that the bubble size and distribution after decompression is directly proportional to the decompression rate. The overall SIMS results of Li and B show that the two trace elements tend to progressively decrease their concentration with decreasing decompression rates. This is explained because for longer decompression times, the diffusion of Li and B into the bubbles has more time to progress and the melt continuously loses volatiles as the bubbles expand their volumes. For fast decompression, Li and B results show a concentration increase with a δ7Li and δ11B decrease close to the bubble interface, related to the sudden formation of the gas bubble, and the occurrence of a diffusion process in the opposite direction, from the bubble meniscus to the unaltered melt. When the bubble growth becomes dominant and Li and B start to exsolve into the gas phase, the silicate melt close to the bubble gets depleted in Li and B, because of a stronger diffusion of the trace elements into the bubble. Our data are being applied to different models, aiming to combine the dynamics of bubble nucleation and growth with the evolution of trace elements concentration and isotopic ratios. Here, first considerations on these models will be presented, giving concluding remarks on this research study. All in all, the final remarks constitute a good starting point for further investigations. These results are a promising base to continue to study this process, and Li and B can indeed show clear dependences on decompression-related magma ascent rates in volcanic systems.}, language = {en} } @article{SongJieGaoetal.2022, author = {Song, Lina and Jie, Dongmei and Gao, Guizai and Liu, Lidan and Liu, Hongyan and Li, Dehui and Liu, Ying}, title = {Application of a topsoil phytolith dataset to quantitative paleoclimate reconstruction in Northeast China}, series = {Palaeogeography, palaeoclimatology, palaeoecology : an international journal for the geo-sciences}, volume = {601}, journal = {Palaeogeography, palaeoclimatology, palaeoecology : an international journal for the geo-sciences}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0031-0182}, doi = {10.1016/j.palaeo.2022.111108}, pages = {12}, year = {2022}, abstract = {Although phytoliths are recognized as an important proxy for paleoenvironmental reconstruction, the quantitative relationship between phytoliths and climate is still debated. In order to provide an improved basis for phytolith-based paleoclimate reconstructions, a representative modern phytolith dataset is essential. Here, we synthesize a modern topsoil phytolith dataset for Northeast China, analyze its climatic significance, and apply it to a fossil phytolith series from the Hani peat core in Northeast China. The dataset comprises 660 topsoil phytolith assemblages from 289 sample sites. We compiled modern meteorological data to assess the quantitative relationship between the phytolith assemblages and climatic variables. Detrended correspondence analysis (DCA) and Redundancy analysis (RDA) were used to determine the dominant climatic variable influencing the phytolith distributions. The results showed that mean annual temperature (MAT) is the dominant variable controlling the spatial distribution of phytoliths, accounting for 8.91\% of the total variance. Transfer function based on inverse deshrinking locally-weighted weighted averaging (LWWA_Inv) was developed for MAT (R-_boot(2) = 0.86, RMSEP = 1.02 degrees C). Applying the LWWA_Inv transfer function to fossil phytolith records from the Hani peat core enables quantitative inferences to be made about Holocene climate changes in Northeast China. Overall, combined with the LWWA_Inv method, the topsoil phytolith dataset of Northeast China can be used for reliable quantitative MAT reconstruction.}, language = {en} } @article{SmithTraxlBoers2022, author = {Smith, Taylor and Traxl, Dominik and Boers, Niklas}, title = {Empirical evidence for recent global shifts in vegetation resilience}, series = {Nature climate change}, volume = {12}, journal = {Nature climate change}, number = {5}, publisher = {Nature Publ. Group}, address = {London}, issn = {1758-678X}, doi = {10.1038/s41558-022-01352-2}, pages = {477 -- 484}, year = {2022}, abstract = {The authors demonstrate that a vegetation system's ability to recover from disturbances-its resilience-can be estimated from its natural variability. Global patterns of resilience loss and gains since the early 1990s reveal shifts towards widespread resilience loss since the early 2000s. The character and health of ecosystems worldwide is tightly coupled to changes in Earth's climate. Theory suggests that ecosystem resilience-the ability of ecosystems to resist and recover from external shocks such as droughts and fires-can be inferred from their natural variability. Here, we quantify vegetation resilience globally with complementary metrics based on two independent long-term satellite records. We first empirically confirm that the recovery rates from large perturbations can be closely approximated from internal vegetation variability across vegetation types and climate zones. On the basis of this empirical relationship, we quantify vegetation resilience continuously and globally from 1992 to 2017. Long-term vegetation resilience trends are spatially heterogeneous, with overall increasing resilience in the tropics and decreasing resilience at higher latitudes. Shorter-term trends, however, reveal a marked shift towards a global decline in vegetation resilience since the early 2000s, particularly in the equatorial rainforest belt.}, language = {en} } @article{SinghSinhaMishraetal.2022, author = {Singh, Manudeo and Sinha, Rajiv and Mishra, Arjit and Babu, Suresh}, title = {Wetlandscape (dis)connectivity and fragmentation in a large wetland (Haiderpur) in west Ganga plains, India}, series = {Earth surface processes and landforms : the journal of the British Geomorphological Research Group}, volume = {47}, journal = {Earth surface processes and landforms : the journal of the British Geomorphological Research Group}, number = {7}, publisher = {Wiley}, address = {New York, NY [u.a.]}, issn = {0197-9337}, doi = {10.1002/esp.5352}, pages = {1872 -- 1887}, year = {2022}, abstract = {Wetlands are dynamic ecosystems that require continuous monitoring and assessment of degradation status to design strategies for their sustainable management. While hydrology provides the primary functional control for the wetland ecosystem, the loss of landscape connectivity influences wetland degradation in a major way as it leads to fragmentation. This article aims to integrate hydrogeomorphic and ecological concepts for the assessment of degradation status and its causal factors for a large wetland in the western Ganga plains, India, the Haiderpur, using a wetlandscape approach. We have used a remote-sensing-based approach, which offers a powerful tool for assessing and linking cross-scale structures, functions, and controls in a wetlandscape. The Haiderpur, a Ramsar site since December 2021, is an artificial wetland located on the right bank of the Ganga River wherein the inflows are controlled by a barrage constructed on the Ganga River apart from smaller tributaries flowing in from the north. A novel aspect of this work is the integration of river dynamics and its connectivity to the wetlandscape to understand the spatiotemporal variability in the waterspread area in the wetland. In this work, we have developed an integrated wetlandscape assessment approach by evaluating wetland's geomorphic and hydrological connectivity status for the period 1993-2019 (25 years) across three different spatial scales - regional, catchment, and wetland. We have highlighted the ecological implications of connectivity and patch dynamics for developing sustainable wetland management plans.}, language = {en} } @article{SiegmundFunkSommeretal.2022, author = {Siegmund, Nicole and Funk, Roger and Sommer, Michael and Avecilla, Fernando and Esteban Panebianco, Juan and Iturri, Laura Antonela and Buschiazzo, Daniel}, title = {Horizontal and vertical fluxes of particulate matter during wind erosion on arable land in the province La Pampa, Argentina}, series = {International journal of sediment research}, volume = {37}, journal = {International journal of sediment research}, number = {5}, publisher = {IRTCES}, address = {Beijing}, issn = {1001-6279}, doi = {10.1016/j.ijsrc.2022.01.004}, pages = {539 -- 552}, year = {2022}, abstract = {A detailed analysis of horizontal and vertical particulate matter (PM) fluxes during wind erosion has been done, based on measurements of PM smaller than 10, 2.5, and 1.0 mu mm, at windward and leeward positions on a measuring field. The three fractions of PM measurement are differently influenced by the increasing wind and shear velocities of the wind. The measured concentrations of the coarser fractions of the fine dust, PM10, and PM2.5, increase with wind and shear velocity, whereas the PM1.0 concentrations show no clear correlation to the shear velocity. The share of PM2.5 on PM10 depends on the measurement height and wind speed and varies between 4 and 12 m/s at the 1 m height ranging from 25\% to 7\% (average 10\%), and at the 4 m height from 39\% to 23\% (average 30\%). Although general relationships between wind speed, PM concentration, and horizontal and vertical fluxes could be found, the contribution of the measuring field was very low, as balances of incoming and outgoing fluxes show. Consequently, the measured PM concentrations are determined from a variety of sources, such as traffic on unpaved roads, cattle drives, tillage operations, and wind erosion, and thus, represent all components of land use and landscape structure in the near and far surroundings of the measuring field. The current results may reflect factors from the landscape scale rather than the influence of field-related variables. The measuring devices used to monitor PM concentrations showed differences of up to 20\%, which led to considerable deviations when determining total balances. Differences up to 67\% between the calculated fluxes prove the necessity of a previous calibration of the devices used. (c) 2022 International Research and Training Centre on Erosion and Sedimentation/the World Association for Sedimentation and Erosion Research.}, language = {en} } @phdthesis{Siegmund2022, author = {Siegmund, Nicole}, title = {Wind driven soil particle uptake Quantifying drivers of wind erosion across the particle size spectrum}, doi = {10.25932/publishup-57489}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-574897}, school = {Universit{\"a}t Potsdam}, pages = {ix, 56}, year = {2022}, abstract = {Among the multitude of geomorphological processes, aeolian shaping processes are of special character, Pedogenic dust is one of the most important sources of atmospheric aerosols and therefore regarded as a key player for atmospheric processes. Soil dust emissions, being complex in composition and properties, influence atmospheric processes and air quality and has impacts on other ecosystems. In this because even though their immediate impact can be considered low (exceptions exist), their constant and large-scale force makes them a powerful player in the earth system. dissertation, we unravel a novel scientific understanding of this complex system based on a holistic dataset acquired during a series of field experiments on arable land in La Pampa, Argentina. The field experiments as well as the generated data provide information about topography, various soil parameters, the atmospheric dynamics in the very lower atmosphere (4m height) as well as measurements regarding aeolian particle movement across a wide range of particle size classes between 0.2μm up to the coarse sand. The investigations focus on three topics: (a) the effects of low-scale landscape structures on aeolian transport processes of the coarse particle fraction, (b) the horizontal and vertical fluxes of the very fine particles and (c) the impact of wind gusts on particle emissions. Among other considerations presented in this thesis, it could in particular be shown, that even though the small-scale topology does have a clear impact on erosion and deposition patterns, also physical soil parameters need to be taken into account for a robust statistical modelling of the latter. Furthermore, specifically the vertical fluxes of particulate matter have different characteristics for the particle size classes. Finally, a novel statistical measure was introduced to quantify the impact of wind gusts on the particle uptake and its application on the provided data set. The aforementioned measure shows significantly increased particle concentrations during points in time defined as gust event. With its holistic approach, this thesis further contributes to the fundamental understanding of how atmosphere and pedosphere are intertwined and affect each other.}, language = {en} } @article{SieberYaxleyHermann2022, author = {Sieber, Melanie Jutta and Yaxley, Greg and Hermann, J{\"o}rg}, title = {COH-fluid induced metasomatism of peridotites in the forearc mantle}, series = {Contributions to Mineralogy and Petrology}, volume = {177}, journal = {Contributions to Mineralogy and Petrology}, number = {4}, publisher = {Springer}, address = {New York}, issn = {0010-7999}, doi = {10.1007/s00410-022-01905-w}, pages = {22}, year = {2022}, abstract = {Devolatilization of subducting lithologies liberates COH-fluids. These may become partially sequestered in peridotites in the slab and the overlying forearc mantle, affecting the cycling of volatiles and fluid mobile elements in subduction zones. Here we assess the magnitudes, timescales and mechanism of channelized injection of COH-fluids doped with Ca-aq(2+), Sr-aq(2+) and Ba-aq(2+) into the dry forearc mantle by performing piston cylinder experiments between 1-2.5 GPa and 600-700 degrees C. Cylindrical cores of natural spinel-bearing harzburgites were used as starting materials. Based on mineral assemblage and composition three reaction zones are distinguishable from the rim towards the core of primary olivine and orthopyroxene grains. Zone 1 contains carbonates + quartz +/- kyanite and zone 2 contains carbonates + talc +/- chlorite. Olivine is further replaced in zone 3 by either antigorite+ magnesite or magnesite +talc within or above antigorite stability, respectively. Orthopyroxene is replaced in zone 3 by talc + chlorite. Mineral assemblages and the compositions of secondary minerals depend on fluid composition and the replaced primary silicate. The extent of alteration depends on fluid CO2 content and fluid/rock-ratio, and is further promoted by fluid permeable reaction zones and reaction driven cracking. Our results show that COH-fluid induced metasomatism of the forearc mantle is self-perpetuating and efficient at sequestering Ca-aq(2+), Sr-aq(2+), Ba-aq(2+) and CO2aq into newly formed carbonates. This process is fast with 90\% of the available C sequestered and nearly 50\% of the initial minerals altered at 650 degrees C, 2 GPa within 55 h. The dissolution of primary silicates under high COH-fluid/rock-ratios, as in channelized fluid flow, enriches SiO2aq in the fluid, while CO2aq is sequestered into carbonates. In an open system, the remaining CO2-depleted, Si-enriched aqueous fluid may cause Si-metasomatism in the forearc further away from the injection of the COH-fluid into peridotite.}, language = {en} } @article{ShpritsAllisonWangetal.2022, author = {Shprits, Yuri Y. and Allison, Hayley J. and Wang, Dedong and Drozdov, Alexander and Szabo-Roberts, Matyas and Zhelavskaya, Irina and Vasile, Ruggero}, title = {A new population of ultra-relativistic electrons in the outer radiation zone}, series = {Journal of geophysical research : Space physics}, volume = {127}, journal = {Journal of geophysical research : Space physics}, number = {5}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9380}, doi = {10.1029/2021JA030214}, pages = {34}, year = {2022}, abstract = {Van Allen Probes measurements revealed the presence of the most unusual structures in the ultra-relativistic radiation belts. Detailed modeling, analysis of pitch angle distributions, analysis of the difference between relativistic and ultra-realistic electron evolution, along with theoretical studies of the scattering and wave growth, all indicate that electromagnetic ion cyclotron (EMIC) waves can produce a very efficient loss of the ultra-relativistic electrons in the heart of the radiation belts. Moreover, a detailed analysis of the profiles of phase space densities provides direct evidence for localized loss by EMIC waves. The evolution of multi-MeV fluxes shows dramatic and very sudden enhancements of electrons for selected storms. Analysis of phase space density profiles reveals that growing peaks at different values of the first invariant are formed at approximately the same radial distance from the Earth and show the sequential formation of the peaks from lower to higher energies, indicating that local energy diffusion is the dominant source of the acceleration from MeV to multi-MeV energies. Further simultaneous analysis of the background density and ultra-relativistic electron fluxes shows that the acceleration to multi-MeV energies only occurs when plasma density is significantly depleted outside of the plasmasphere, which is consistent with the modeling of acceleration due to chorus waves.}, language = {en} } @phdthesis{Schoenfeldt2022, author = {Sch{\"o}nfeldt, Elisabeth}, title = {Giant landslides in Patagonia, Argentina}, pages = {XXII, 156}, year = {2022}, language = {en} } @phdthesis{Schuster2022, author = {Schuster, Valerian}, title = {Mechanical and hydraulic properties of Opalinus Clay}, doi = {10.25932/publishup-56678}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-566786}, school = {Universit{\"a}t Potsdam}, year = {2022}, abstract = {Deep geological repositories represent a promising solution for the final disposal of nuclear waste. Due to its low permeability, high sorption capacity and self-sealing potential, Opalinus Clay (OPA) is considered a suitable host rock formation for the long-term storage of nuclear waste in Switzerland and Germany. However, the clay formation is characterized by compositional and structural variabilities including the occurrence of carbonate- and quartz-rich layers, pronounced bedding planes as well as tectonic elements such as pre-existing fault zones and fractures, suggesting heterogeneous rock mass properties. Characterizing the heterogeneity of host rock properties is therefore essential for safety predictions of future repositories. This includes a detailed understanding of the mechanical and hydraulic properties, deformation behavior and the underlying deformation processes for an improved assessment of the sealing integrity and long-term safety of a deep repository in OPA. Against this background, this thesis presents the results of deformation experiments performed on intact and artificially fractured specimens of the quartz-rich, sandy and clay-rich, shaly facies of OPA. The experiments focus on the influence of mineralogical composition on the deformation behavior as well as the reactivation and sealing properties of pre-existing faults and fractures at different boundary conditions (e.g., pressure, temperature, strain rate). The anisotropic mechanical properties of the sandy facies of OPA are presented in the first section, which were determined from triaxial deformation experiments using dried and resaturated samples loaded at 0°, 45° and 90° to the bedding plane orientation. A Paterson-type deformation apparatus was used that allowed to investigate how the deformation behavior is influenced by the variation of confining pressure (50 - 100 MPa), temperature (25 - 200 °C), and strain rate (1 × 10-3 - 5 × 10-6 s-1). Constant strain rate experiments revealed brittle to semi-brittle deformation behavior of the sandy facies at the applied conditions. Deformation behavior showed a strong dependence on confining pressure, degree of water saturation as well as bedding orientation, whereas the variation of temperature and strain rate had no significant effect on deformation. Furthermore, the sandy facies displays higher strength and stiffness compared to the clay-rich shaly facies deformed at similar conditions by N{\"u}esch (1991). From the obtained results it can be concluded that cataclastic mechanisms dominate the short-term deformation behavior of dried samples from both facies up to elevated pressure (<200 MPa) and temperature (<200 °C) conditions. The second part presents triaxial deformation tests that were performed to investigate how structural discontinuities affect the deformation behavior of OPA and how the reactivation of preexisting faults is influenced by mineral composition and confining pressure. To this end, dried cylindrical samples of the sandy and shaly facies of OPA were used, which contained a saw-cut fracture oriented at 30° to the long axis. After hydrostatic pre-compaction at 50 MPa, constant strain rate deformation tests were performed at confining pressures of 5, 20 or 35 MPa. With increasing confinement, a gradual transition from brittle, highly localized fault slip including a stress drop at fault reactivation to semi-brittle deformation behavior, characterized by increasing delocalization and non-linear strain hardening without dynamic fault reactivation, can be observed. Brittle localization was limited by the confining pressure at which the fault strength exceeded the matrix yield strength, above which strain partitioning between localized fault slip and distributed matrix deformation occurred. The sandy facies displayed a slightly higher friction coefficient (≈0.48) compared to the shaly facies (≈0.4). In addition, slide-hold-slide tests were conducted, revealing negative or negligible frictional strengthening, which suggests stable creep and long-term weakness of faults in both facies of OPA. The conducted experiments demonstrate that dilatant brittle fault reactivation in OPA may be favored at high overconsolidation ratios and shallow depths, increasing the risk of seismic hazard and the creation of fluid pathways. The final section illustrates how the sealing capacity of fractures in OPA is affected by mineral composition. Triaxial flow-through experiments using Argon-gas were performed with dried samples from the sandy and shaly facies of OPA containing a roughened, artificial fracture. Slate, graywacke, quartzite, natural fault gouge, and granite samples were also tested to highlight the influence of normal stress, mineralogy and diagenesis on the sustainability of fracture transmissivity. With increasing normal stress, a non-linear decrease of fracture transmissivity can be observed that resulted in a permanent reduction of transmissivity after stress release. The transmissivity of rocks with a high portion of strong minerals (e.g., quartz) and high unconfined compressive strength was less sensitive to stress changes. In accordance with this, the sandy facies of OPA displayed a higher initial transmissivity that was less sensitive to stress changes compared to the shaly facies. However, transmissivity of rigid slate was less sensitive to stress changes than the sandy facies of OPA, although the slate is characterized by a higher phyllosilicate content. This demonstrates that in addition to mineral composition, other factors such as the degree of metamorphism, cementation and consolidation have to be considered when evaluating the sealing capacity of phyllosilicate-rich rocks. The results of this thesis highlighted the role of confining pressure on the failure behavior of intact and artificially fractured OPA. Although the quartz-rich sandy facies may be considered as being more favorable for underground constructions due to its higher shear strength and stiffness than the shaly facies, the results indicate that when fractures develop in the sandy facies, they are more conductive and remain more permeable compared to fractures in the clay-dominated shaly facies at a given stress. The results may provide the basis for constitutive models to predict the integrity and evolution of a future repository. Clearly, the influence of composition and consolidation, e.g., by geological burial and uplift, on the mechanical sealing behavior of OPA highlights the need for a detailed site-specific material characterization for a future repository.}, language = {en} } @article{ScholzVoigt2022, author = {Scholz, Carolin and Voigt, Christian C.}, title = {Diet analysis of bats killed at wind turbines suggests large-scale losses of trophic interactions}, series = {Conservation science and practice}, volume = {4}, journal = {Conservation science and practice}, number = {7}, publisher = {Wiley}, address = {Hoboken}, issn = {2578-4854}, doi = {10.1111/csp2.12744}, pages = {12}, year = {2022}, abstract = {Agricultural practice has led to landscape simplification and biodiversity decline, yet recently, energy-producing infrastructures, such as wind turbines, have been added to these simplified agroecosystems, turning them into multi-functional energy-agroecosystems. Here, we studied the trophic interactions of bats killed at wind turbines using a DNA metabarcoding approach to shed light on how turbine-related bat fatalities may possibly affect local habitats. Specifically, we identified insect DNA in the stomachs of common noctule bats (Nyctalus noctula) killed by wind turbines in Germany to infer in which habitats these bats hunted. Common noctule bats consumed a wide variety of insects from different habitats, ranging from aquatic to terrestrial ecosystems (e.g., wetlands, farmland, forests, and grasslands). Agricultural and silvicultural pest insects made up about 20\% of insect species consumed by the studied bats. Our study suggests that the potential damage of wind energy production goes beyond the loss of bats and the decline of bat populations. Bat fatalities at wind turbines may lead to the loss of trophic interactions and ecosystem services provided by bats, which may add to the functional simplification and impaired crop production, respectively, in multi-functional ecosystems.}, language = {en} } @misc{ScholzVoigt2022, author = {Scholz, Carolin and Voigt, Christian C.}, title = {Diet analysis of bats killed at wind turbines suggests large-scale losses of trophic interactions}, series = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {7}, issn = {1866-8372}, doi = {10.25932/publishup-59156}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-591568}, pages = {14}, year = {2022}, abstract = {Agricultural practice has led to landscape simplification and biodiversity decline, yet recently, energy-producing infrastructures, such as wind turbines, have been added to these simplified agroecosystems, turning them into multi-functional energy-agroecosystems. Here, we studied the trophic interactions of bats killed at wind turbines using a DNA metabarcoding approach to shed light on how turbine-related bat fatalities may possibly affect local habitats. Specifically, we identified insect DNA in the stomachs of common noctule bats (Nyctalus noctula) killed by wind turbines in Germany to infer in which habitats these bats hunted. Common noctule bats consumed a wide variety of insects from different habitats, ranging from aquatic to terrestrial ecosystems (e.g., wetlands, farmland, forests, and grasslands). Agricultural and silvicultural pest insects made up about 20\% of insect species consumed by the studied bats. Our study suggests that the potential damage of wind energy production goes beyond the loss of bats and the decline of bat populations. Bat fatalities at wind turbines may lead to the loss of trophic interactions and ecosystem services provided by bats, which may add to the functional simplification and impaired crop production, respectively, in multi-functional ecosystems.}, language = {en} } @misc{SchmidtFranckeRottleretal.2022, author = {Schmidt, Lena Katharina and Francke, Till and Rottler, Erwin and Blume, Theresa and Sch{\"o}ber, Johannes}, title = {Suspended sediment and discharge dynamics in a glaciated alpine environment: identifying crucial areas and time periods on several spatial and temporal scales in the {\"O}tztal, Austria}, series = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {1296}, issn = {1866-8372}, doi = {10.25932/publishup-57656}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-576564}, pages = {653 -- 669}, year = {2022}, abstract = {Glaciated high-alpine areas are fundamentally altered by climate change, with well-known implications for hydrology, e.g., due to glacier retreat, longer snow-free periods, and more frequent and intense summer rainstorms. While knowledge on how these hydrological changes will propagate to suspended sediment dynamics is still scarce, it is needed to inform mitigation and adaptation strategies. To understand the processes and source areas most relevant to sediment dynamics, we analyzed discharge and sediment dynamics in high temporal resolution as well as their patterns on several spatial scales, which to date few studies have done. We used a nested catchment setup in the Upper {\"O}tztal in Tyrol, Austria, where high-resolution (15 min) time series of discharge and suspended sediment concentrations are available for up to 15 years (2006-2020). The catchments of the gauges in Vent, S{\"o}lden and Tumpen range from 100 to almost 800 km2 with 10 \% to 30 \% glacier cover and span an elevation range of 930 to 3772 m a.s.l. We analyzed discharge and suspended sediment yields (SSY), their distribution in space, their seasonality and spatial differences therein, and the relative importance of short-term events. We complemented our analysis by linking the observations to satellite-based snow cover maps, glacier inventories, mass balances and precipitation data. Our results indicate that the areas above 2500 m a.s.l., characterized by glacier tongues and the most recently deglaciated areas, are crucial for sediment generation in all sub-catchments. This notion is supported by the synchronous spring onset of sediment export at the three gauges, which coincides with snowmelt above 2500 m but lags behind spring discharge onsets. This points at a limitation of suspended sediment supply as long as the areas above 2500 m are snow-covered. The positive correlation of annual SSY with glacier cover (among catchments) and glacier mass balances (within a catchment) further supports the importance of the glacier-dominated areas. The analysis of short-term events showed that summer precipitation events were associated with peak sediment concentrations and yields but on average accounted for only 21 \% of the annual SSY in the headwaters. These results indicate that under current conditions, thermally induced sediment export (through snow and glacier melt) is dominant in the study area. Our results extend the scientific knowledge on current hydro-sedimentological conditions in glaciated high-alpine areas and provide a baseline for studies on projected future changes in hydro-sedimentological system dynamics.}, language = {en} } @article{SchmidPetersenHooftetal.2022, author = {Schmid, Florian and Petersen, Gesa M. and Hooft, Emilie E. E. and Paulatto, Michele and Chrapkiewicz, Kajetan and Hensch, Martin and Dahm, Torsten}, title = {Heralds of future volcanism: Swarms of microseismicity beneath the submarine Kolumbo volcano indicate opening of near-vertical fractures exploited by ascending melts}, series = {Geochemistry, geophysics, geosystems}, volume = {23}, journal = {Geochemistry, geophysics, geosystems}, number = {7}, publisher = {American Geophysical Union}, address = {Washington}, issn = {1525-2027}, doi = {10.1029/2022GC010420}, pages = {21}, year = {2022}, abstract = {The Kolumbo submarine volcano in the southern Aegean (Greece) is associated with repeated seismic unrest since at least two decades and the causes of this unrest are poorly understood. We present a ten-month long microseismicity data set for the period 2006-2007. The majority of earthquakes cluster in a cone-shaped portion of the crust below Kolumbo. The tip of this cone coincides with a low Vp-anomaly at 2-4 km depth, which is interpreted as a crustal melt reservoir. Our data set includes several earthquake swarms, of which we analyze the four with the highest events numbers in detail. Together the swarms form a zone of fracturing elongated in the SW-NE direction, parallel to major regional faults. All four swarms show a general upward migration of hypocenters and the cracking front propagates unusually fast, compared to swarms in other volcanic areas. We conclude that the swarm seismicity is most likely triggered by a combination of pore-pressure perturbations and the re-distribution of elastic stresses. Fluid pressure perturbations are induced likely by obstructions in the melt conduits in a rheologically strong layer between 6 and 9 km depth. We conclude that the zone of fractures below Kolumbo is exploited by melts ascending from the mantle and filling the crustal melt reservoir. Together with the recurring seismic unrest, our study suggests that a future eruption is probable and monitoring of the Kolumbo volcanic system is highly advisable.}, language = {en} } @article{SchildgenvanderBeekD'Arcyetal.2022, author = {Schildgen, Taylor F. and van der Beek, Pieter A. and D'Arcy, Mitch and Roda-Boluda, Duna N. and Orr, Elizabeth N. and Wittmann, Hella}, title = {Quantifying drainage-divide migration from orographic rainfall over geologic timescales}, series = {Earth \& planetary science letters}, volume = {579}, journal = {Earth \& planetary science letters}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0012-821X}, doi = {10.1016/j.epsl.2021.117345}, pages = {13}, year = {2022}, abstract = {Drainage-divide migration, controlled by rock-uplift and rainfall patterns, may play a major role in the geomorphic evolution of mountain ranges. However, divide-migration rates over geologic timescales have only been estimated by theoretical studies and remain empirically poorly constrained. Geomorphological evidence suggests that the Sierra de Aconquija, on the eastern side of the southern Central Andes, northwest Argentina, is undergoing active westward drainage-divide migration. The mountain range has been subjected to steep rock trajectories and pronounced orographic rainfall for the last several million years, presenting an ideal setting for using low-temperature thermochronometric data to explore its topographic evolution. We perform three-dimensional thermal-kinematic modeling of previously published thermochronometric data spanning the windward and leeward sides of the range to explore the most likely structural and topographic evolution of the range. We find that the data can be explained by scenarios involving drainage-divide migration alone, or by scenarios that also involve changes in the structures that have accommodated deformation through time. By combining new Be-10-derived catchment-average denudation rates with geomorphic constraints on probable fault activity, we conclude that the evolution of the range was likely dominated by west-vergent faulting on a high-angle reverse fault underlying the range, together with westward drainage-divide migration at a rate of several km per million years. Our findings place new constraints on the magnitudes and rates of drainage-divide migration in real landscapes, quantify the effects of orographic rainfall and erosion on the topographic evolution of a mountain range, and highlight the importance of considering drainage-divide migration when interpreting thermochronometer age patterns.}, language = {en} } @article{SchifferleLobanov2022, author = {Schifferle, Lukas and Lobanov, Sergey S.}, title = {Evolution of chemical bonding and spin-pairing energy in ferropericlase across Its spin transition}, series = {ACS Earth and Space Chemistry}, volume = {6}, journal = {ACS Earth and Space Chemistry}, number = {3}, publisher = {American Chemical Society}, address = {Washington}, issn = {2472-3452}, doi = {10.1021/acsearthspacechem.2c00014}, pages = {788 -- 799}, year = {2022}, abstract = {The evolution of chemical bonding in ferropericlase, (Mg,Fe)O, with pressure may affect the physical and chemical properties of the Earth's lower mantle. Here, we report high-pressure optical absorption spectra of single-crystalline ferropericlase ((Mg0.87Fe0.13)O) up to 135 GPa. Combined with a re-evaluation of published partial fluorescence yield X-ray absorption spectroscopy data, we show that the covalency of the Fe-O bond increases with pressure, but the iron spin transition at 57-76.5 GPa reverses this trend. The qualitative crossover in chemical bonding suggests that the spin-pairing transition weakens the Fe-O bond in ferropericlase. We find, that the spin transition in ferropericlase is caused by both the increase of the ligand field-splitting energy and the decrease in the spin-pairing energy of high-spin Fe2+.}, language = {en} } @article{SarrazinKumarBasuetal.2022, author = {Sarrazin, Fanny J. and Kumar, Rohini and Basu, Nandita B. and Musolff, Andreas and Weber, Michael and Van Meter, Kimberly J. and Attinger, Sabine}, title = {Characterizing catchment-scale nitrogen legacies and constraining their uncertainties}, series = {Water resources research}, volume = {58}, journal = {Water resources research}, number = {4}, publisher = {American Geophysical Union}, address = {Washington}, issn = {0043-1397}, doi = {10.1029/2021WR031587}, pages = {32}, year = {2022}, abstract = {Improving nitrogen (N) status in European water bodies is a pressing issue. N levels depend not only on current but also past N inputs to the landscape, that have accumulated through time in legacy stores (e.g., soil, groundwater). Catchment-scale N models, that are commonly used to investigate in-stream N levels, rarely examine the magnitude and dynamics of legacy components. This study aims to gain a better understanding of the long-term fate of the N inputs and its uncertainties, using a legacy-driven N model (ELEMeNT) in Germany's largest national river basin (Weser; 38,450 km(2)) over the period 1960-2015. We estimate the nine model parameters based on a progressive constraining strategy, to assess the value of different observational data sets. We demonstrate that beyond in-stream N loading, soil N content and in-stream N concentration allow to reduce the equifinality in model parameterizations. We find that more than 50\% of the N surplus denitrifies (1480-2210 kg ha(-1)) and the stream export amounts to around 18\% (410-640 kg ha(-1)), leaving behind as much as around 230-780 kg ha(-1) of N in the (soil) source zone and 10-105 kg ha(-1) in the subsurface. A sensitivity analysis reveals the importance of different factors affecting the residual uncertainties in simulated N legacies, namely hydrologic travel time, denitrification rates, a coefficient characterizing the protection of organic N in source zone and N surplus input. Our study calls for proper consideration of uncertainties in N legacy characterization, and discusses possible avenues to further reduce the equifinality in water quality modeling.}, language = {en} } @article{SarrDonnadieuBoltonetal.2022, author = {Sarr, Anta-Clarisse and Donnadieu, Yannick and Bolton, Clara T. and Ladant, Jean-Baptiste and Licht, Alexis and Fluteau, Fr{\´e}d{\´e}ric and Laugi{\´e}, Marie and Tardif, Delphine and Dupont-Nivet, Guillaume}, title = {Neogene South Asian monsoon rainfall and wind histories diverged due to topographic effects}, series = {Nature geoscience}, volume = {15}, journal = {Nature geoscience}, number = {4}, publisher = {Nature Research}, address = {Berlin}, issn = {1752-0894}, doi = {10.1038/s41561-022-00919-0}, pages = {314 -- 319}, year = {2022}, abstract = {The drivers of the evolution of the South Asian Monsoon remain widely debated. An intensification of monsoonal rainfall recorded in terrestrial and marine sediment archives from the earliest Miocene (23-20 million years ago (Ma)) is generally attributed to Himalayan uplift. However, Indian Ocean palaeorecords place the onset of a strong monsoon around 13 Ma, linked to strengthening of the southwesterly winds of the Somali Jet that also force Arabian Sea upwelling. Here we reconcile these divergent records using Earth system model simulations to evaluate the interactions between palaeogeography and ocean-atmosphere dynamics. We show that factors forcing the South Asian Monsoon circulation versus rainfall are decoupled and diachronous. Himalayan and Tibetan Plateau topography predominantly controlled early Miocene rainfall patterns, with limited impact on ocean-atmosphere circulation. The uplift of the East African and Middle Eastern topography played a pivotal role in the establishment of the modern Somali Jet structure above the western Indian Ocean, while strong upwelling initiated as a direct consequence of the emergence of the Arabian Peninsula and the onset of modern-like atmospheric circulation. Our results emphasize that although elevated rainfall seasonality was probably a persistent feature since the India-Asia collision in the Paleogene, modern-like monsoonal atmospheric circulation only emerged in the late Neogene.}, language = {en} } @article{RolphOverduinRavensetal.2022, author = {Rolph, Rebecca and Overduin, Pier Paul and Ravens, Thomas and Lantuit, Hugues and Langer, Moritz}, title = {ArcticBeach v1.0}, series = {Frontiers in Earth Science}, volume = {10}, journal = {Frontiers in Earth Science}, publisher = {Frontiers Media}, address = {Lausanne}, issn = {2296-6463}, doi = {10.3389/feart.2022.962208}, pages = {19}, year = {2022}, abstract = {In the Arctic, air temperatures are increasing and sea ice is declining, resulting in larger waves and a longer open water season, all of which intensify the thaw and erosion of ice-rich coasts. Climate change has been shown to increase the rate of Arctic coastal erosion, causing problems for Arctic cultural heritage, existing industrial, military, and civil infrastructure, as well as changes in nearshore biogeochemistry. Numerical models that reproduce historical and project future Arctic erosion rates are necessary to understand how further climate change will affect these problems, and no such model yet exists to simulate the physics of erosion on a pan-Arctic scale. We have coupled a bathystrophic storm surge model to a simplified physical erosion model of a permafrost coastline. This Arctic erosion model, called ArcticBeach v1.0, is a first step toward a physical parameterization of Arctic shoreline erosion for larger-scale models. It is forced by wind speed and direction, wave period and height, sea surface temperature, all of which are masked during times of sea ice cover near the coastline. Model tuning requires observed historical retreat rates (at least one value), as well as rough nearshore bathymetry. These parameters are already available on a pan-Arctic scale. The model is validated at three study sites at 1) Drew Point (DP), Alaska, 2) Mamontovy Khayata (MK), Siberia, and 3) Veslebogen Cliffs, Svalbard. Simulated cumulative retreat rates for DP and MK respectively (169 and 170 m) over the time periods studied at each site (2007-2016, and 1995-2018) are found to the same order of magnitude as observed cumulative retreat (172 and 120 m). The rocky Veslebogen cliffs have small observed cumulative retreat rates (0.05 m over 2014-2016), and our model was also able to reproduce this same order of magnitude of retreat (0.08 m). Given the large differences in geomorphology between the study sites, this study provides a proof-of-concept that ArcticBeach v1.0 can be applied on very different permafrost coastlines. ArcticBeach v1.0 provides a promising starting point to project retreat of Arctic shorelines, or to evaluate historical retreat in places that have had few observations.}, language = {en} } @article{RodriguezMoskwaOsesetal.2022, author = {Rodriguez, Victoria and Moskwa, Lisa-Marie and Oses, Romulo and K{\"u}hn, Peter and Riveras-Mu{\~n}oz, Nicol{\´a}s and Seguel, Oscar and Scholten, Thomas and Wagner, Dirk}, title = {Impact of climate and slope aspects on the composition of soil bacterial communities involved in pedogenetic processes along the chilean coastal cordillera}, series = {Microorganisms}, volume = {10}, journal = {Microorganisms}, number = {5}, publisher = {MDPI}, address = {Basel}, issn = {2076-2607}, doi = {10.3390/microorganisms10050847}, pages = {20}, year = {2022}, abstract = {Soil bacteria play a fundamental role in pedogenesis. However, knowledge about both the impact of climate and slope aspects on microbial communities and the consequences of these items in pedogenesis is lacking. Therefore, soil-bacterial communities from four sites and two different aspects along the climate gradient of the Chilean Coastal Cordillera were investigated. Using a combination of microbiological and physicochemical methods, soils that developed in arid, semi-arid, mediterranean, and humid climates were analyzed. Proteobacteria, Acidobacteria, Chloroflexi, Verrucomicrobia, and Planctomycetes were found to increase in abundance from arid to humid climates, while Actinobacteria and Gemmatimonadetes decreased along the transect. Bacterial-community structure varied with climate and aspect and was influenced by pH, bulk density, plant-available phosphorus, clay, and total organic-matter content. Higher bacterial specialization was found in arid and humid climates and on the south-facing slope and was likely promoted by stable microclimatic conditions. The presence of specialists was associated with ecosystem-functional traits, which shifted from pioneers that accumulated organic matter in arid climates to organic decomposers in humid climates. These findings provide new perspectives on how climate and slope aspects influence the composition and functional capabilities of bacteria, with most of these capabilities being involved in pedogenetic processes.}, language = {en} } @article{RodriguezPicedaScheckWenderothCacaceetal.2022, author = {Rodriguez Piceda, Constanza and Scheck-Wenderoth, Magdalena and Cacace, Mauro and Bott, Judith and Strecker, Manfred}, title = {Long-Term Lithospheric Strength and Upper-Plate Seismicity in the Southern Central Andes, 29 degrees-39 degrees S}, series = {Geochemistry, geophysics, geosystems}, volume = {23}, journal = {Geochemistry, geophysics, geosystems}, number = {3}, publisher = {American Geophysical Union}, address = {Washington}, issn = {1525-2027}, doi = {10.1029/2021GC010171}, pages = {22}, year = {2022}, abstract = {We examined the relationship between the mechanical strength of the lithosphere and the distribution of seismicity within the overriding continental plate of the southern Central Andes (SCA, 29 degrees-39 degrees S), where the oceanic Nazca Plate changes its subduction angle between 33 degrees S and 35 degrees S, from subhorizontal in the north (<5 degrees) to steep in the south (similar to 30 degrees). We computed the long-term lithospheric strength based on an existing 3D model describing variations in thickness, density, and temperature of the main geological units forming the lithosphere of the SCA and adjacent forearc and foreland regions. The comparison between our results and seismicity within the overriding plate (upper-plate seismicity) shows that most of the events occur within the modeled brittle domain of the lithosphere. The depth where the deformation mode switches from brittle frictional to thermally activated ductile creep provides a conservative lower bound to the seismogenic zone in the overriding plate of the study area. We also found that the majority of upper-plate earthquakes occurs within the realm of first-order contrasts in integrated strength (12.7-13.3 log Pam in the Andean orogen vs. 13.5-13.9 log Pam in the forearc and the foreland). Specific conditions characterize the mechanically strong northern foreland of the Andes, where seismicity is likely explained by the effects of slab steepening.}, language = {en} } @article{RodriguezPicedaScheckWenderothBottetal.2022, author = {Rodriguez Piceda, Constanza and Scheck-Wenderoth, Magdalena and Bott, Judith and Gomez Dacal, Maria Laura and Cacace, Mauro and Pons, Michael and Prezzi, Claudia and Strecker, Manfred}, title = {Controls of the Lithospheric Thermal Field of an Ocean-Continent Subduction Zone}, series = {Lithosphere / Geological Society of America}, volume = {2022}, journal = {Lithosphere / Geological Society of America}, number = {1}, publisher = {GeoScienceWorld}, address = {McLean}, issn = {1941-8264}, doi = {10.2113/2022/2237272}, pages = {26}, year = {2022}, abstract = {In an ocean-continent subduction zone, the assessment of the lithospheric thermal state is essential to determine the controls of the deformation within the upper plate and the dip angle of the subducting lithosphere. In this study, we evaluate the degree of influence of both the configuration of the upper plate (i.e., thickness and composition of the rock units) and variations of the subduction angle on the lithospheric thermal field of the southern Central Andes (29 degrees-39 degrees S). Here, the subduction angle increases from subhorizontal (5 degrees) north of 33 degrees S to steep (similar to 30 degrees) in the south. We derived the 3D temperature and heat flow distribution of the lithosphere in the southern Central Andes considering conversion of S wave tomography to temperatures together with steady-state conductive thermal modeling. We found that the orogen is overall warmer than the forearc and the foreland and that the lithosphere of the northern part of the foreland appears colder than its southern counterpart. Sedimentary blanketing and the thickness of the radiogenic crust exert the main control on the shallow thermal field (<50km depth). Specific conditions are present where the oceanic slab is relatively shallow (<85 km depth) and the radiogenic crust is thin. This configuration results in relatively colder temperatures compared to regions where the radiogenic crust is thick and the slab is steep. At depths >50km, the temperatures of the overriding plate are mainly controlled by the mantle heat input and the subduction angle. The thermal field of the upper plate likely preserves the flat subduction angle and influences the spatial distribution of shortening.}, language = {en} } @phdthesis{RodriguezPiceda2022, author = {Rodriguez Piceda, Constanza}, title = {Thermomechanical state of the southern Central Andes}, doi = {10.25932/publishup-54927}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-549275}, school = {Universit{\"a}t Potsdam}, pages = {xx, 228}, year = {2022}, abstract = {The Andes are a ~7000 km long N-S trending mountain range developed along the South American western continental margin. Driven by the subduction of the oceanic Nazca plate beneath the continental South American plate, the formation of the northern and central parts of the orogen is a type case for a non-collisional orogeny. In the southern Central Andes (SCA, 29°S-39°S), the oceanic plate changes the subduction angle between 33°S and 35°S from almost horizontal (< 5° dip) in the north to a steeper angle (~30° dip) in the south. This sector of the Andes also displays remarkable along- and across- strike variations of the tectonic deformation patterns. These include a systematic decrease of topographic elevation, of crustal shortening and foreland and orogenic width, as well as an alternation of the foreland deformation style between thick-skinned and thin-skinned recorded along- and across the strike of the subduction zone. Moreover, the SCA are a very seismically active region. The continental plate is characterized by a relatively shallow seismicity (< 30 km depth) which is mainly focussed at the transition from the orogen to the lowland areas of the foreland and the forearc; in contrast, deeper seismicity occurs below the interiors of the northern foreland. Additionally, frequent seismicity is also recorded in the shallow parts of the oceanic plate and in a sector of the flat slab segment between 31°S and 33°S. The observed spatial heterogeneity in tectonic and seismic deformation in the SCA has been attributed to multiple causes, including variations in sediment thickness, the presence of inherited structures and changes in the subduction angle of the oceanic slab. However, there is no study that inquired the relationship between the long-term rheological configuration of the SCA and the spatial deformation patterns. Moreover, the effects of the density and thickness configuration of the continental plate and of variations in the slab dip angle in the rheological state of the lithosphere have been not thoroughly investigated yet. Since rheology depends on composition, pressure and temperature, a detailed characterization of the compositional, structural and thermal fields of the lithosphere is needed. Therefore, by using multiple geophysical approaches and data sources, I constructed the following 3D models of the SCA lithosphere: (i) a seismically-constrained structural and density model that was tested against the gravity field; (ii) a thermal model integrating the conversion of mantle shear-wave velocities to temperature with steady-state conductive calculations in the uppermost lithosphere (< 50 km depth), validated by temperature and heat-flow measurements; and (iii) a rheological model of the long-term lithospheric strength using as input the previously-generated models. The results of this dissertation indicate that the present-day thermal and rheological fields of the SCA are controlled by different mechanisms at different depths. At shallow depths (< 50 km), the thermomechanical field is modulated by the heterogeneous composition of the continental lithosphere. The overprint of the oceanic slab is detectable where the oceanic plate is shallow (< 85 km depth) and the radiogenic crust is thin, resulting in overall lower temperatures and higher strength compared to regions where the slab is steep and the radiogenic crust is thick. At depths > 50 km, largest temperatures variations occur where the descending slab is detected, which implies that the deep thermal field is mainly affected by the slab dip geometry. The outcomes of this thesis suggests that long-term thermomechanical state of the lithosphere influences the spatial distribution of seismic deformation. Most of the seismicity within the continental plate occurs above the modelled transition from brittle to ductile conditions. Additionally, there is a spatial correlation between the location of these events and the transition from the mechanically strong domains of the forearc and foreland to the weak domain of the orogen. In contrast, seismicity within the oceanic plate is also detected where long-term ductile conditions are expected. I therefore analysed the possible influence of additional mechanisms triggering these earthquakes, including the compaction of sediments in the subduction interface and dehydration reactions in the slab. To that aim, I carried out a qualitative analysis of the state of hydration in the mantle using the ratio between compressional- and shear-wave velocity (vp/vs ratio) from a previous seismic tomography. The results from this analysis indicate that the majority of the seismicity spatially correlates with hydrated areas of the slab and overlying continental mantle, with the exception of the cluster within the flat slab segment. In this region, earthquakes are likely triggered by flexural processes where the slab changes from a flat to a steep subduction angle. First-order variations in the observed tectonic patterns also seem to be influenced by the thermomechanical configuration of the lithosphere. The mechanically strong domains of the forearc and foreland, due to their resistance to deformation, display smaller amounts of shortening than the relatively weak orogenic domain. In addition, the structural and thermomechanical characteristics modelled in this dissertation confirm previous analyses from geodynamic models pointing to the control of the observed heterogeneities in the orogen and foreland deformation style. These characteristics include the lithospheric and crustal thickness, the presence of weak sediments and the variations in gravitational potential energy. Specific conditions occur in the cold and strong northern foreland, which is characterized by active seismicity and thick-skinned structures, although the modelled crustal strength exceeds the typical values of externally-applied tectonic stresses. The additional mechanisms that could explain the strain localization in a region that should resist deformation are: (i) increased tectonic forces coming from the steepening of the slab and (ii) enhanced weakening along inherited structures from pre-Andean deformation events. Finally, the thermomechanical conditions of this sector of the foreland could be a key factor influencing the preservation of the flat subduction angle at these latitudes of the SCA.}, language = {en} } @article{RiedlMelnickNjueetal.2022, author = {Riedl, Simon and Melnick, Daniel and Njue, Lucy and Sudo, Masafumi and Strecker, Manfred}, title = {Mid-Pleistocene to recent crustal extension in the inner graben of the Northern Kenya Rift}, series = {Geochemistry, geophysics, geosystems}, volume = {23}, journal = {Geochemistry, geophysics, geosystems}, number = {3}, publisher = {American Geophysical Union}, address = {Washington}, issn = {1525-2027}, doi = {10.1029/2021GC010123}, pages = {25}, year = {2022}, abstract = {Magmatic continental rifts often constitute nascent plate boundaries, yet long-term extension rates and transient rate changes associated with these early stages of continental breakup remain difficult to determine. Here, we derive a time-averaged minimum extension rate for the inner graben of the Northern Kenya Rift (NKR) of the East African Rift System for the last 0.5 m.y. We use the TanDEM-X science digital elevation model to evaluate fault-scarp geometries and determine fault throws across the volcano-tectonic axis of the inner graben of the NKR. Along rift-perpendicular profiles, amounts of cumulative extension are determined, and by integrating four new Ar-40/Ar-39 radiometric dates for the Silali volcano into the existing geochronology of the faulted volcanic units, time-averaged extension rates are calculated. This study reveals that in the inner graben of the NKR, the long-term extension rate based on mid-Pleistocene to recent brittle deformation has minimum values of 1.0-1.6 mm yr(-1), locally with values up to 2.0 mm yr(-1). A comparison with the decadal, geodetically determined extension rate reveals that at least 65\% of the extension must be accommodated within a narrow, 20-km-wide zone of the inner rift. In light of virtually inactive border faults of the NKR, we show that extension is focused in the region of the active volcano-tectonic axis in the inner graben, thus highlighting the maturing of continental rifting in the NKR.}, language = {en} } @phdthesis{Richter2022, author = {Richter, Maximilian Jacob Enzo Amandus}, title = {Continental rift dynamics across the scales}, doi = {10.25932/publishup-55060}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-550606}, school = {Universit{\"a}t Potsdam}, pages = {129}, year = {2022}, abstract = {Localisation of deformation is a ubiquitous feature in continental rift dynamics and observed across drastically different time and length scales. This thesis comprises one experimental and two numerical modelling studies investigating strain localisation in (1) a ductile shear zone induced by a material heterogeneity and (2) in an active continental rift setting. The studies are related by the fact that the weakening mechanisms on the crystallographic and grain size scale enable bulk rock weakening, which fundamentally enables the formation of shear zones, continental rifts and hence plate tectonics. Aiming to investigate the controlling mechanisms on initiation and evolution of a shear zone, the torsion experiments of the experimental study were conducted in a Patterson type apparatus with strong Carrara marble cylinders with a weak, planar Solnhofen limestone inclusion. Using state-of-the-art numerical modelling software, the torsion experiments were simulated to answer questions regarding localisation procedure like stress distribution or the impact of rheological weakening. 2D numerical models were also employed to integrate geophysical and geological data to explain characteristic tectonic evolution of the Southern and Central Kenya Rift. Key elements of the numerical tools are a randomized initial strain distribution and the usage of strain softening. During the torsion experiments, deformation begins to localise at the limestone inclusion tips in a process zone, which propagates into the marble matrix with increasing deformation until a ductile shear zone is established. Minor indicators for coexisting brittle deformation are found close to the inclusion tip and presumed to slightly facilitate strain localisation besides the dominant ductile deformation processes. The 2D numerical model of the torsion experiment successfully predicts local stress concentration and strain rate amplification ahead of the inclusion in first order agreement with the experimental results. A simple linear parametrization of strain weaking enables high accuracy reproduction of phenomenological aspects of the observed weakening. The torsion experiments suggest that loading conditions do not affect strain localisation during high temperature deformation of multiphase material with high viscosity contrasts. A numerical simulation can provide a way of analysing the process zone evolution virtually and extend the examinable frame. Furthermore, the nested structure and anastomosing shape of an ultramylonite band was mimicked with an additional second softening step. Rheological weakening is necessary to establish a shear zone in a strong matrix around a weak inclusion and for ultramylonite formation. Such strain weakening laws are also incorporated into the numerical models of the Southern and Central Kenya Rift that capture the characteristic tectonic evolution. A three-stage early rift evolution is suggested that starts with (1) the accommodation of strain by a single border fault and flexure of the hanging-wall crust, after which (2) faulting in the hanging-wall and the basin centre increases before (3) the early-stage asymmetry is lost and basinward localisation of deformation occurs. Along-strike variability of rifts can be produced by modifying the initial random noise distribution. In summary, the three studies address selected aspects of the broad range of mechanisms and processes that fundamentally enable the deformation of rock and govern the localisation patterns across the scales. In addition to the aforementioned results, the first and second manuscripts combined, demonstrate a procedure to find new or improve on existing numerical formulations for specific rheologies and their dynamic weakening. These formulations are essential in addressing rock deformation from the grain to the global scale. As within the third study of this thesis, where geodynamic controls on the evolution of a rift were examined and acquired by the integration of geological and geophysical data into a numerical model.}, language = {en} } @article{RibackiTrumbullLopezDeLuchietal.2022, author = {Ribacki, Enrico and Trumbull, Robert B. and Lopez De Luchi, Monica Graciela and Altenberger, Uwe}, title = {The chemical and B-Isotope composition of Tourmaline from intra-granitic Pegmatites in the Las Chacras-Potrerillos Batholith, Argentina}, series = {The Canadian mineralogist : journal of the Mineralogical Association of Canada}, volume = {60}, journal = {The Canadian mineralogist : journal of the Mineralogical Association of Canada}, number = {1}, publisher = {Association of Canada}, address = {Ottawa}, issn = {0008-4476}, doi = {10.3749/canmin.2100036}, pages = {49 -- 66}, year = {2022}, abstract = {The Devonian Las Chacras-Potrerillos batholith comprises six nested monzonitic to granitic intrusions with metaluminous to weakly peraluminous composition and a Sr-Nd isotopic signature indicating a dominantly juvenile mantle-derived source. The chemically most evolved units in the southern batholith contain a large number of intra-granitic, pod-shaped tourmaline-bearing pegmatites. This study uses in situ chemical and boron isotopic analyses of tourmaline from nine of these pegmatites to discuss their relationship to the respective host intrusions and the implications of their B-isotope composition for the source and evolution of the magmas. The tourmalines reveal a diversity in element composition (e.g., FeO, MgO, TiO2, CaO, MnO, F) which distinguishes individual pegmatites from one another. However, all have a narrow 5 11 B range of -13.7 to -10.5\%0 (n = 100) which indicates a relatively uniform magmatic system and similar temperature conditions during tourmaline crystallization. The average delta(11) B value of -11.7\%0 is typical for S-type granites and is within the range reported for peraluminous granites. pegmatites, and metamorphic units of the Ordovician basement into which the Las Chacras-Potrerillos batholith intruded. The B-isotope evidence argues for a crustal boron source like that of the Ordovician basement, in contrast to the metaluminous to weakly peraluminous composition and juvenile initial Sr and Nd isotope ratios of the Las Chacras-Potrerillos batholith magmas. We propose that the boron was not derived from the magma source region but was incorporated from dehydration melting of elastic metasedimentary rocks higher up in the crustal column.}, language = {en} } @misc{RembeSobelKleyetal.2022, author = {Rembe, Johannes and Sobel, Edward and Kley, Jonas and Terbishalieva, Baiansulu and Musiol, Antje and Chen, Jie and Zhou, Renjie}, title = {Geochronology, Geochemistry, and Geodynamic Implications of Permo-Triassic Back-Arc Basin Successions in the North Pamir, Central Asia}, series = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {1309}, issn = {1866-8372}, doi = {10.25932/publishup-58331}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-583318}, pages = {21}, year = {2022}, abstract = {The Permo-Triassic period marks the time interval between Hercynian (Variscan) orogenic events in the Tien Shan and the North Pamir, and the Cimmerian accretion of the Gondwana-derived Central and South Pamir to the southern margin of the Paleo-Asian continent. A well-preserved Permo-Triassic volcano-sedimentary sequence from the Chinese North Pamir yields important information on the geodynamic evolution of Asia's pre-Cimmerian southern margin. The oldest volcanic rocks from that section are dated to the late Guadalupian epoch by a rhyolite and a dacitic dike that gave zircon U-Pb ages of ~260 Ma. Permian volcanism was largely pyroclastic and mafic to intermediate. Upsection, a massive ignimbritic crystal tuff in the Chinese Qimgan valley was dated to 244.1 +/- 1.1 Ma, a similar unit in the nearby Gez valley to 245 +/- 11 Ma, and an associated rhyolite to 233.4 +/- 1.1 Ma. Deposition of the locally ~200 m thick crystal tuff unit follows an unconformity and marks the onset of intense, mainly mafic to intermediate, calc-alkaline magmatic activity. Triassic volcanic activity in the North Pamir was coeval with the major phase of Cimmerian intrusive activity in the Karakul-Mazar arc-accretionary complex to the south, caused by northward subduction of the Paleo-Tethys. It also coincided with the emplacement of basanitic and carbonatitic dikes and a thermal event in the South Tien Shan, to the north of our study area. Evidence for arc-related magmatic activity in a back-arc position provides strong arguments for back-arc extension or transtension and basin formation. This puts the Qimgan succession in line with a more than 1000 km long realm of extensional Triassic back-arc basins known from the North Pamir in the Kyrgyz Altyn Darya valley (Myntekin formation), the North Pamir of Tajikistan and Afghanistan, and the Afghan Hindukush (Doab formation) and further west from the Paropamisus and Kopet Dag (Aghdarband, NE Iran).}, language = {en} } @article{RembeSobelKleyetal.2022, author = {Rembe, Johannes and Sobel, Edward and Kley, Jonas and Terbishalieva, Baiansulu and Musiol, Antje and Chen, Jie and Zhou, Renjie}, title = {Geochronology, Geochemistry, and Geodynamic Implications of Permo-Triassic Back-Arc Basin Successions in the North Pamir, Central Asia}, series = {Lithosphere}, volume = {2022}, journal = {Lithosphere}, number = {1}, publisher = {GeoScienceWorld, Geological Society of America}, address = {Boulder, Colorado, USA}, issn = {1947-4253}, doi = {10.2113/2022/7514691}, pages = {21}, year = {2022}, abstract = {The Permo-Triassic period marks the time interval between Hercynian (Variscan) orogenic events in the Tien Shan and the North Pamir, and the Cimmerian accretion of the Gondwana-derived Central and South Pamir to the southern margin of the Paleo-Asian continent. A well-preserved Permo-Triassic volcano-sedimentary sequence from the Chinese North Pamir yields important information on the geodynamic evolution of Asia's pre-Cimmerian southern margin. The oldest volcanic rocks from that section are dated to the late Guadalupian epoch by a rhyolite and a dacitic dike that gave zircon U-Pb ages of ~260 Ma. Permian volcanism was largely pyroclastic and mafic to intermediate. Upsection, a massive ignimbritic crystal tuff in the Chinese Qimgan valley was dated to 244.1 +/- 1.1 Ma, a similar unit in the nearby Gez valley to 245 +/- 11 Ma, and an associated rhyolite to 233.4 +/- 1.1 Ma. Deposition of the locally ~200 m thick crystal tuff unit follows an unconformity and marks the onset of intense, mainly mafic to intermediate, calc-alkaline magmatic activity. Triassic volcanic activity in the North Pamir was coeval with the major phase of Cimmerian intrusive activity in the Karakul-Mazar arc-accretionary complex to the south, caused by northward subduction of the Paleo-Tethys. It also coincided with the emplacement of basanitic and carbonatitic dikes and a thermal event in the South Tien Shan, to the north of our study area. Evidence for arc-related magmatic activity in a back-arc position provides strong arguments for back-arc extension or transtension and basin formation. This puts the Qimgan succession in line with a more than 1000 km long realm of extensional Triassic back-arc basins known from the North Pamir in the Kyrgyz Altyn Darya valley (Myntekin formation), the North Pamir of Tajikistan and Afghanistan, and the Afghan Hindukush (Doab formation) and further west from the Paropamisus and Kopet Dag (Aghdarband, NE Iran).}, language = {en} } @article{RegmiBookhagen2022, author = {Regmi, Shakil and Bookhagen, Bodo}, title = {The spatial pattern of extreme precipitation from 40 years of gauge data in the central Himalaya}, series = {Weather and climate extremes}, volume = {37}, journal = {Weather and climate extremes}, publisher = {Elsevier}, address = {Amsterdam}, issn = {2212-0947}, doi = {10.1016/j.wace.2022.100470}, pages = {14}, year = {2022}, abstract = {The topography of the Himalaya exerts a substantial control on the spatial distribution of monsoonal rainfall, which is a vital water source for the regional economy and population. But the occurrence of short-lived and high-intensity precipitation results in socio-economic losses. This study relies on 40 years of daily data from 204 ground stations in Nepal to derive extreme precipitation thresholds, amounts, and days at the 95th percentile. We additionally determine the precipitation magnitude-frequency relation. We observe that extreme precipitation amounts follow an almost uniform band parallel to topographic contour lines in the southern Himalaya mountains in central and eastern Nepal but not in western Nepal. The relationship of extreme precipitation indices with topographic relief shows that extreme precipitation thresholds decrease with increasing elevation, but extreme precipitation days increase in higher elevation areas. Furthermore, stations above 1 km elevation exhibit a power-law relation in the rainfall magnitude-frequency framework. Stations at higher elevations generally have lower values of power-law exponents than low elevation areas. This suggests a fundamentally different behaviour of the rainfall distribution and an increased occurrence of extreme rainfall storms in the high elevation areas of Nepal.}, language = {en} } @article{RasigrafWagner2022, author = {Rasigraf, Olivia and Wagner, Dirk}, title = {Landslides}, series = {Earth science reviews : the international geological journal bridging the gap between research articles and textbooks}, volume = {231}, journal = {Earth science reviews : the international geological journal bridging the gap between research articles and textbooks}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0012-8252}, doi = {10.1016/j.earscirev.2022.104064}, pages = {16}, year = {2022}, abstract = {Erosion by landslides is a common phenomenon in mountain regions around the globe, affecting all climatic zones. Landslides facilitate bedrock weathering, pedogenesis and ecological succession, being key drivers of biodiversity. Landslide chronosequences have long been used for studies of vegetation succession in initial ecosystems, but they further offer ideal model systems for studies of soil development and microbial community succession. In this review we synthesize the state of knowledge on the role of landslides in ecosystems, their influence on element cycles and interactions with biota. Further, we discuss feedback mechanisms between global warming, landslide activity and greenhouse gas emissions. In the view of increasing anthropogenic influence and climate change, soils are becoming a critical resource. Due to their ubiquity, landslide chronosequences have the potential to provide critical insights into soil development under different climates and thereby contribute to future soil restoration efforts.}, language = {en} } @article{RamachandranRupakhetiCherianetal.2022, author = {Ramachandran, Srikanthan and Rupakheti, Maheswar and Cherian, R. and Lawrence, Mark}, title = {Climate Benefits of Cleaner Energy Transitions in East and South Asia Through Black Carbon Reduction}, series = {Frontiers in environmental science}, volume = {10}, journal = {Frontiers in environmental science}, publisher = {Frontiers Media}, address = {Lausanne}, issn = {2296-665X}, doi = {10.3389/fenvs.2022.842319}, pages = {16}, year = {2022}, abstract = {The state of air pollution has historically been tightly linked to how we produce and use energy. Air pollutant emissions over Asia are now changing rapidly due to cleaner energy transitions; however, magnitudes of benefits for climate and air quality remain poorly quantified. The associated risks involve adverse health impacts, reduced agricultural yields, reduced freshwater availability, contributions to climate change, and economic costs. We focus particularly on climate benefits of energy transitions by making first-time use of two decades of high quality observations of atmospheric loading of light-absorbing black carbon (BC) over Kanpur (South Asia) and Beijing (East Asia) and relating these observations to changing energy, emissions, and economic trends in India and China. Our analysis reveals that absorption aerosol optical depth (AAOD) due to BC has decreased substantially, by 40\% over Kanpur and 60\% over Beijing between 2001 and 2017, and thus became decoupled from regional economic growth. Furthermore, the resultant decrease in BC emissions and BC AAOD over Asia is regionally coherent and occurs primarily due to transitions into cleaner energies (both renewables and fossil fuels) and not due to the decrease in primary energy supply or decrease in use of fossil use and biofuels and waste. Model simulations show that BC aerosols alone contribute about half of the surface temperature change (warming) of the total forcing due to greenhouse gases, natural and internal variability, and aerosols, thus clearly revealing the climate benefits due to a reduction in BC emissions, which would significantly reduce global warming. However, this modeling study excludes responses from natural variability, circulation, and sea ice responses, which cause relatively strong temperature fluctuations that may mask signals from BC aerosols. Our findings show additional benefits for climate (beyond benefits of CO2 reduction) and for several other issues of sustainability over South and East Asia, provide motivation for ongoing cleaner energy production, and consumption transitions, especially when they are associated with reduced emissions of air pollutants. Such an analysis connecting the trends in energy transitions and aerosol absorption loading, unavailable so far, is crucial for simulating the aerosol climate impacts over Asia which is quite uncertain.}, language = {en} } @article{RadosavljevicLantuitKnoblauchetal.2022, author = {Radosavljevic, Boris and Lantuit, Hugues and Knoblauch, Christian and Couture, Nicole and Herzschuh, Ulrike and Fritz, Michael}, title = {Arctic nearshore sediment dynamics - an example from Herschel Island - Qikiqtaruk, Canada}, series = {Journal of marine science and engineering}, volume = {10}, journal = {Journal of marine science and engineering}, number = {11}, publisher = {MDPI}, address = {Basel}, issn = {2077-1312}, doi = {10.3390/jmse10111589}, pages = {18}, year = {2022}, abstract = {Increasing arctic coastal erosion rates imply a greater release of sediments and organic matter into the coastal zone. With 213 sediment samples taken around Herschel Island-Qikiqtaruk, Canadian Beaufort Sea, we aimed to gain new insights on sediment dynamics and geochemical properties of a shallow arctic nearshore zone. Spatial characteristics of nearshore sediment texture (moderately to poorly sorted silt) are dictated by hydrodynamic processes, but ice-related processes also play a role. We determined organic matter (OM) distribution and inferred the origin and quality of organic carbon by C/N ratios and stable carbon isotopes delta C-13. The carbon content was higher offshore and in sheltered areas (mean: 1.0 wt.\%., S.D.: 0.9) and the C/N ratios also showed a similar spatial pattern (mean: 11.1, S.D.: 3.1), while the delta C-13 (mean: -26.4 parts per thousand VPDB, S.D.: 0.4) distribution was more complex. We compared the geochemical parameters of our study with terrestrial and marine samples from other studies using a bootstrap approach. Sediments of the current study contained 6.5 times and 1.8 times less total organic carbon than undisturbed and disturbed terrestrial sediments, respectively. Therefore, degradation of OM and separation of carbon pools take place on land and continue in the nearshore zone, where OM is leached, mineralized, or transported beyond the study area.}, language = {en} } @article{ProcykLovejoyHebert2022, author = {Procyk, Roman and Lovejoy, Shaun and H{\´e}bert, Raphaёl}, title = {The fractional energy balance equation for climate projections through 2100}, series = {Earth system dynamics / European Geosciences Union}, volume = {13}, journal = {Earth system dynamics / European Geosciences Union}, number = {1}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {2190-4979}, doi = {10.5194/esd-13-81-2022}, pages = {81 -- 107}, year = {2022}, abstract = {We produce climate projections through the 21st century using the fractional energy balance equation (FEBE): a generalization of the standard energy balance equation (EBE). The FEBE can be derived from Budyko-Sellers models or phenomenologically through the application of the scaling symmetry to energy storage processes, easily implemented by changing the integer order of the storage (derivative) term in the EBE to a fractional value. The FEBE is defined by three parameters: a fundamental shape parameter, a timescale and an amplitude, corresponding to, respectively, the scaling exponent h, the relaxation time tau and the equilibrium climate sensitivity (ECS). Two additional parameters were needed for the forcing: an aerosol recalibration factor alpha to account for the large aerosol uncertainty and a volcanic intermittency correction exponent upsilon. A Bayesian framework based on historical temperatures and natural and anthropogenic forcing series was used for parameter estimation. Significantly, the error model was not ad hoc but rather predicted by the model itself: the internal variability response to white noise internal forcing. The 90 \% credible interval (CI) of the exponent and relaxation time were h = [0.33, 0.44] (median = 0.38) and tau = [2.4, 7.0] (median = 4.7) years compared to the usual EBE h = 1, and literature values of tau typically in the range 2-8 years. Aerosol forcings were too strong, requiring a decrease by an average factor alpha = [0.2, 1.0] (median = 0.6); the volcanic intermittency correction exponent was upsilon = [0.15, 0.41] (median = 0.28) compared to standard values alpha = upsilon = 1. The overpowered aerosols support a revision of the global modern (2005) aerosol forcing 90 \% CI to a narrower range [ -1.0, -0.2] W m(-2). The key parameter ECS in comparison to IPCC AR5 (and to the CMIP6 MME), the 90 \% CI range is reduced from [1.5, 4.5] K ([2.0, 5.5] K) to [1.6, 2.4] K ([1.5, 2.2] K), with median value lowered from 3.0 K (3.7 K) to 2.0 K (1.8 K) Similarly we found for the transient climate response (TCR), the 90 \% CI range shrinks from [1.0, 2.5] K ([1.2, 2.8] K) to [1.2, 1.8] K ([1.1, 1.6] K) and the median estimate decreases from 1.8 K (2.0 K) to 1.5 K (1.4 K). As often seen in other observational-based studies, the FEBE values for climate sensitivities are therefore somewhat lower but still consistent with those in IPCC AR5 and the CMIP6 MME.
Using these parameters, we made projections to 2100 using both the Representative Concentration Pathway (RCP) and Shared Socioeconomic Pathway (SSP) scenarios, and compared them to the corresponding CMIP5 and CMIP6 multi-model ensembles (MMEs). The FEBE historical reconstructions (1880-2020) closely follow observations, notably during the 1998-2014 slowdown ("hiatus"). We also reproduce the internal variability with the FEBE and statistically validate this against centennial-scale temperature observations. Overall, the FEBE projections were 10 \%-15 \% lower but due to their smaller uncertainties, their 90 \% CIs lie completely within the GCM 90 \% CIs. This agreement means that the FEBE validates the MME, and vice versa.}, language = {en} } @article{PorȩbaRacioppiGarbarinoetal.2022, author = {Porȩba, Tomasz and Racioppi, Stefano and Garbarino, Gaston and Morgenroth, Wolfgang and Mezouar, Mohamed}, title = {Investigating the structural symmetrization of CsI3 at high pressures through combined X-ray diffraction experiments and theoretical analysis}, series = {Inorganic chemistry}, volume = {61}, journal = {Inorganic chemistry}, number = {28}, publisher = {American Chemical Society}, address = {Washington}, issn = {0020-1669}, doi = {10.1021/acs.inorgchem.2c01690}, pages = {10977 -- 10985}, year = {2022}, abstract = {ABSTRACT: Structural evolution of cesium triiodide at high pressures has been revealed by synchrotron single-crystal X-ray diffraction. Cesium triiodide undergoes a first-order phase transition above 1.24(3) GPa from an orthorhombic to a trigonal system. This transition is coupled with severe reorganization of the polyiodide network from a layered to three-dimensional architecture. Quantum chemical calculations show that even though the two polymorphic phases are nearly isoenergetic under ambient conditions, the PV term is decisive in stabilizing the trigonal polymorph above the transition point. Phonon calculations using a non-local correlation functional that accounts for dispersion interactions confirm that this polymorph is dynamically unstable under ambient conditions. The high-pressure behavior of crystalline CsI3 can be correlated with other alkali metal trihalides, which undergo a similar sequence of structural changes upon load.}, language = {en} } @article{PlatzWeckmannPeketal.2022, author = {Platz, Anna and Weckmann, Ute and Pek, Josef and Kovacikova, Svetlana and Klanica, Radek and Mair, Johannes and Aleid, Basel}, title = {3D imaging of the subsurface electrical resistivity structure in West Bohemia/Upper Palatinate covering mofettes and quaternary volcanic structures by using magnetotellurics}, series = {Tectonophysics : international journal of geotectonics and the geology and physics of the interior of the earth}, volume = {833}, journal = {Tectonophysics : international journal of geotectonics and the geology and physics of the interior of the earth}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0040-1951}, doi = {10.1016/j.tecto.2022.229353}, pages = {20}, year = {2022}, abstract = {The region of West Bohemia and Upper Palatinate belongs to the West Bohemian Massif. The study area is situated at the junction of three different Variscan tectonic units and hosts the ENE-WSW trending Ohre Rift as well as many different fault systems. The entire region is characterized by ongoing magmatic processes in the intra-continental lithospheric mantle expressed by a series of phenomena, including e.g. the occurrence of repeated earthquake swarms and massive degassing of mantle derived CO2 in form of mineral springs and mofettes. Ongoing active tectonics is mainly manifested by Cenozoic volcanism represented by different Quaternary volcanic structures. All these phenomena make the Ohre Rift a unique target area for European intra-continental geo-scientific research. With magnetotelluric (MT) measurements we image the subsurface distribution of the electrical resistivity and map possible fluid pathways. Two-dimensional (2D) inversion results by Munoz et al. (2018) reveal a conductive channel in the vicinity of the earthquake swarm region that extends from the lower crust to the surface forming a pathway for fluids into the region of the mofettes. A second conductive channel is present in the south of their model; however, their 2D inversions allow ambiguous interpretations of this feature. Therefore, we conducted a large 3D MT field experiment extending the study area towards the south. The 3D inversion result matches well with the known geology imaging different fluid/magma reservoirs at crust-mantle depth and mapping possible fluid pathways from the reservoirs to the surface feeding known mofettes and spas. A comparison of 3D and 2D inversion results suggests that the 2D inversion results are considerably characterized by 3D and off-profile structures. In this context, the new results advocate for the swarm earthquakes being located in the resistive host rock surrounding the conductive channels; a finding in line with observations e.g. at the San Andreas Fault, California.}, language = {en} } @article{PenaMetzgerHeidbachetal.2022, author = {Pe{\~n}a, Carlos and Metzger, Sabrina and Heidbach, Oliver and Bedford, Jonathan and Bookhagen, Bodo and Moreno, Marcos and Oncken, Onno and Cotton, Fabrice}, title = {Role of poroelasticity during the early postseismic deformation of the 2010 Maule megathrust earthquake}, series = {Geophysical research letters}, volume = {49}, journal = {Geophysical research letters}, number = {9}, publisher = {Wiley}, address = {Hoboken, NJ}, issn = {0094-8276}, doi = {10.1029/2022GL098144}, pages = {11}, year = {2022}, abstract = {Megathrust earthquakes impose changes of differential stress and pore pressure in the lithosphere-asthenosphere system that are transiently relaxed during the postseismic period primarily due to afterslip, viscoelastic and poroelastic processes. Especially during the early postseismic phase, however, the relative contribution of these processes to the observed surface deformation is unclear. To investigate this, we use geodetic data collected in the first 48 days following the 2010 Maule earthquake and a poro-viscoelastic forward model combined with an afterslip inversion. This model approach fits the geodetic data 14\% better than a pure elastic model. Particularly near the region of maximum coseismic slip, the predicted surface poroelastic uplift pattern explains well the observations. If poroelasticity is neglected, the spatial afterslip distribution is locally altered by up to +/- 40\%. Moreover, we find that shallow crustal aftershocks mostly occur in regions of increased postseismic pore-pressure changes, indicating that both processes might be mechanically coupled.}, language = {en} } @phdthesis{Patyniak2022, author = {Patyniak, Magda}, title = {Seismotectonic segmentation, paleoearthquakes and style of deformation along the Northern Pamir thrust system, South Kyrgyzstan}, doi = {10.25932/publishup-57714}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-577141}, school = {Universit{\"a}t Potsdam}, pages = {xxii, 165}, year = {2022}, abstract = {The Pamir Frontal Thrust (PFT) located in the Trans Alai range in Central Asia is the principal active fault of the intracontinental India-Eurasia convergence zone and constitutes the northernmost boundary of the Pamir orogen at the NW edge of this collision zone. Frequent seismic activity and ongoing crustal shortening reflect the northward propagation of the Pamir into the intermontane Alai Valley. Quaternary deposits are being deformed and uplifted by the advancing thrust front of the Trans Alai range. The Alai Valley separates the Pamir range front from the Tien Shan mountains in the north; the Alai Valley is the vestige of a formerly contiguous basin that linked the Tadjik Depression in the west with the Tarim Basin in the east. GNSS measurements across the Central Pamir document a shortening rate of ~25 mm/yr, with a dramatic decrease of ~10-15 mm over a short distance across the northernmost Trans Alai range. This suggests that almost half of the shortening in the greater Pamir - Tien Shan collision zone is absorbed along the PFT. The short-term (geodetic) and long-term (geologic) shortening rates across the northern Pamir appear to be at odds with an apparent slip-rate discrepancy along the frontal fault system of the Pamir. Moreover, the present-day seismicity and historical records have not revealed great Mw > 7 earthquakes that might be expected with such a significant slip accommodation. In contrast, recent and historic earthquakes exhibit complex rupture patterns within and across seismotectonic segments bounding the Pamir mountain front, challenging our understanding of fault interaction and the seismogenic potential of this area, and leaving the relationships between seismicity and the geometry of the thrust front not well understood. In this dissertation I employ different approaches to assess the seismogenic behavior along the PFT. Firstly, I provide paleoseismic data from five trenches across the central PFT segment (cPFT) and compute a segment-wide earthquake chronology over the past 16 kyr. This novel dataset provides important insights into the recurrence, magnitude, and rupture extent of past earthquakes along the cPFT. I interpret five, possibly six paleoearthquakes that have ruptured the Pamir mountain front since ∼7 ka and 16 ka, respectively. My results indicate that at least three major earthquakes ruptured the full-segment length and possibly crossed segment boundaries with a recurrence interval of ∼1.9 kyr and potential magnitudes of up to Mw 7.4. Importantly, I did not find evidence for great (i.e., Mw ≥8) earthquakes. Secondly, I combine my paleoseimic results with morphometric analyses to establish a segment-wide distribution of the cumulative vertical separation along offset fluvial terraces and I model a long-term slip rate for the cPFT. My investigations reveal discrepancies between the extents of slip and rupture during apparent partial segment ruptures in the western half of the cPFT. Combined with significantly higher fault scarp offsets in this sector of the cPFT, the observations indicate a more mature fault section with a potential for future fault linkage. I estimate an average rate of horizontal motion for the cPFT of 4.1 ± 1.5 mm/yr during the past ∼5 kyr, which does not fully match the GNSS-derived present-day shortening rate of ∼10 mm/yr. This suggests a complex distribution of strain accumulation and potential slip partitioning between the cPFT and additional faults and folds within the Pamir that may be associated with a partially locked regional d{\´e}collement. The third part of the thesis provides new insights regarding the surface rupture of the 2008 Mw 6.6 Nura earthquake that ruptured along the eastern PFT sector. I explore this rupture in the context of its structural complexity by combining extensive field observations with high-resolution digital surface models. I provide a map of the rupture extent, net slip measurements, and updated regional geological observations. Based on this data I propose a tectonic model in this area associated with secondary flexural-slip faulting along steeply dipping bedding of folded Paleogene sedimentary strata that is related to deformation along a deeper blind thrust. Here, the strain release seems to be transferred from the PFT towards older inherited basement structures within the area of advanced Pamir-Tien Shan collision zone. The extensive research of my dissertation results in a paleoseismic database of the past 16 ~kyr, which contributes to the understanding of the seismogenic behavior of the PFT, but also to that of segmented thrust-fault systems in active collisional settings. My observations underscore the importance of combining different methodological approaches in the geosciences, especially in structurally complex tectonic settings like the northern Pamir. Discrepancy between GNSS-derived present-day deformation rates and those from different geological archives in the central part, as well as the widespread distribution of the deformation due to earthquake triggered strain transfer in the eastern part reveals the complexity of this collision zone and calls for future studies involving multi-temporal and interdisciplinary approaches.}, language = {en} } @phdthesis{Pan2022, author = {Pan, Mengdi}, title = {Systematic studies on the thermodynamic properties of gas hydrates and their formation/dissociation/transformation behaviors}, doi = {10.25932/publishup-55476}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-554760}, school = {Universit{\"a}t Potsdam}, pages = {XX, 192}, year = {2022}, abstract = {Gas hydrates are ice-like crystalline compounds made of water cavities that retain various types of guest molecules. Natural gas hydrates are CH4-rich but also contain higher hydrocarbons as well as CO2, H2S, etc. They are highly dependent of local pressure and temperature conditions. Considering the high energy content, natural gas hydrates are artificially dissociated for the production of methane gas. Besides, they may also dissociate in response to global warming. It is therefore crucial to investigate the hydrate nucleation and growth process at a molecular level. The understanding of how guest molecules in the hydrate cavities respond to warming climate or gas injection is also of great importance. This thesis is concerned with a systematic investigation of simple and mixed gas hydrates at conditions relevant to the natural hydrate reservoir in Qilian Mountain permafrost, China. A high-pressure cell that integrated into the confocal Raman spectroscopy ensured a precise and continuous characterization of the hydrate phase during formation/dissociation/transformation processes with a high special and spectral resolution. By applying laboratory experiments, the formation of mixed gas hydrates containing other hydrocarbons besides methane was simulated in consideration of the effects from gas supply conditions and sediments. The results revealed a preferential enclathration of different guest molecules in hydrate cavities and further refute the common hypothesis of the coexistence of hydrate phases due to a changing feed gas phase. However, the presence of specific minerals and organic compounds in sediments may have significant impacts on the coexisting solid phases. With regard to the dissociation, the formation damage caused by fines mobilization and migration during hydrate decomposition was reported for the first time, illustrating the complex interactions between fine grains and hydrate particles. Gas hydrates, starting from simple CH4 hydrates to binary CH4—C3H8 hydrates and multi-component mixed hydrates were decomposed by thermal stimulation mimicking global warming. The mechanisms of guest substitution in hydrate structures were studied through the experimental data obtained from CH4—CO2, CH4—mixed gas hydrates and mixed gas hydrates—CO2 systems. For the first time, a second transformation behavior was documented during the transformation process from CH4 hydrates to CO2-rich mixed hydrates. Most of the crystals grew or maintained when exposed to CO2 gas while some others decreased in sizes and even disappeared over time. The highlight of the two last experimental simulations was to visualize and characterize the hydrate crystals which were at different structural transition stages. These experimental simulations enhanced our knowledge about the mixed gas hydrates in natural reservoirs and improved our capability to assess the response to global warming.}, language = {en} } @article{OuDaoutWeissetal.2022, author = {Ou, Qi and Daout, Simon and Weiss, Jonathan R. and Shen, Lin and Lazecky, Milan and Wright, Tim J. and Parsons, Barry E.}, title = {Large-Scale interseismic strain mapping of the NE Tibetan Plateau from Sentinel-1 Interferometry}, series = {Journal of geophysical research : Solid earth}, volume = {127}, journal = {Journal of geophysical research : Solid earth}, number = {6}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9313}, doi = {10.1029/2022JB024176}, pages = {29}, year = {2022}, abstract = {The launches of the Sentinel-1 synthetic aperture radar satellites in 2014 and 2016 started a new era of high-resolution velocity and strain rate mapping for the continents. However, multiple challenges exist in tying independently processed velocity data sets to a common reference frame and producing high-resolution strain rate fields. We analyze Sentinel-1 data acquired between 2014 and 2019 over the northeast Tibetan Plateau, and develop new methods to derive east and vertical velocities with similar to 100 m resolution and similar to 1 mm/yr accuracy across an area of 440,000 km(2). By implementing a new method of combining horizontal gradients of filtered east and interpolated north velocities, we derive the first similar to 1 km resolution strain rate field for this tectonically active region. The strain rate fields show concentrated shear strain along the Haiyuan and East Kunlun Faults, and local contractional strain on fault junctions, within the Qilianshan thrusts, and around the Longyangxia Reservoir. The Laohushan-Jingtai creeping section of the Haiyuan Fault is highlighted in our data set by extremely rapid strain rates. Strain across unknown portions of the Haiyuan Fault system, including shear on the eastern extension of the Dabanshan Fault and contraction at the western flank of the Quwushan, highlight unmapped tectonic structures. In addition to the uplift across most of the lowlands, the vertical velocities also contain climatic, hydrological or anthropogenic-related deformation signals. We demonstrate the enhanced view of large-scale active tectonic processes provided by high-resolution velocities and strain rates derived from Sentinel-1 data and highlight associated wide-ranging research applications.}, language = {en} } @article{OrtizSaezAlvaradoetal.2022, author = {Ortiz, Gustavo and Saez, Mauro and Alvarado, Patricia and Rivas, Carolina and Garc{\´i}a, V{\´i}ctor Hugo and Alonso, Ricardo and Zullo, Fernando Morales}, title = {Seismotectonic characterization of the 1948 (M-W 6.9) Anta earthquake Santa Barbara System, central Andes broken foreland of northwestern Argentina}, series = {Journal of South American earth sciences}, volume = {116}, journal = {Journal of South American earth sciences}, publisher = {Elsevier}, address = {Oxford}, issn = {0895-9811}, doi = {10.1016/j.jsames.2022.103822}, pages = {15}, year = {2022}, abstract = {The region of the Andean back-arc of northwestern Argentina has been struck by several magnitude >= 6 crustal earthquakes since the first historically recorded event in 1692. One of these events corresponds to the Anta earthquake on 25 August 1948, with epicenter in the Santa Barbara System causing three deaths and severe damage in Salta and Jujuy provinces with maximum Modified Mercalli seismic intensities (MMI) of IX. We collected and digitized analog seismograms of this earthquake from worldwide seismic observatories in order to perform first-motion analysis and modeling of long-period teleseismic P-waveforms. Our results indicate a simple seismic source of M0 = 2.85 x 1019 N m consistent with a moment magnitude Mw = 6.9. We have also tested for the focal depth determining a shallow source at 8 km with a reverse focal mechanism solution with a minor dextral strike-slip component (strike 20 degrees, dip 30 degrees, rake 120 degrees) from the best fit of waveforms. Using magnitude size empirical relationships, the comparison of the obtained Mw 6.9 magnitude value and the ca. 10,000 km2 area of MMI >= IX from our seismic intensity map, which was obtained from newspaper and many historical reports, indicates a rupture length of 42 +/- 8 km for the Anta earthquake. We show our results in a 3D geological model around the epicentral area, which integrates modern seismicity, geological data, and information of a previously studied east-west cross section located a few kilometers south of the 1948 epicenter. The integration of all available information provides evidence of the re-activation of the Pie de la Sierra del Gallo fault during the 1948 Mw 6.9 shallow earthquake; this thrust fault bounds the Santa Barbara System along its western foothill.}, language = {en} } @article{NomosatryoTjallingiiHennyetal.2022, author = {Nomosatryo, Sulung and Tjallingii, Rik and Henny, Cynthia and Ridwansyah, Iwan and Wagner, Dirk and Tom{\´a}s, Sara and Kallmeyer, Jens}, title = {Surface sediment composition and depositional environments in tropical Lake Sentani, Papua Province, Indonesia}, series = {Journal of Paleolimnology}, journal = {Journal of Paleolimnology}, publisher = {Springer}, address = {Dordrecht}, issn = {0921-2728}, doi = {10.1007/s10933-022-00259-4}, pages = {20}, year = {2022}, abstract = {Tropical Lake Sentani in the Indonesian Province Papua consists of four separate basins and is surrounded by a catchment with a very diverse geology. We characterized the surface sediment (upper 5 cm) of the lake's four sub-basins based on multivariate statistical analyses (principal component analysis, hierarchical clustering) of major element compositions obtained by X-ray fluorescence scanning. Three types of sediment are identified based on distinct compositional differences between rivers, shallow/proximal and deep/distal lake sediments. The different sediment types are mainly characterized by the correlation of elements associated with redox processes (S, Mn, Fe), carbonates (Ca), and detrital input (Ti, Al, Si, K) derived by river discharge. The relatively coarse-grained river sediments mainly derive form the mafic catchment geology and contribution of the limestone catchment geology is only limited. Correlation of redox sensitive and detrital elements are used to reveal oxidation conditions, and indicate oxic conditions in river samples and reducing conditions for lake sediments. Organic carbon (TOC) generally correlates with redox sensitive elements, although a correlation between TOC and individual elements change strongly between the three sediment types. Pyrite is the quantitatively dominant reduced sulfur mineral, monosulfides only reach appreciable concentrations in samples from rivers draining mafic and ultramafic catchments. Our study shows large spatial heterogeneity within the lake's sub-basins that is mainly caused by catchment geology and topography, river runoff as well as the bathymetry and the depth of the oxycline. We show that knowledge about lateral heterogeneity is crucial for understanding the geochemical and sedimentological variations recorded by these sediments. The highly variable conditions make Lake Sentani a natural laboratory, with its different sub-basins representing different depositional environments under identical tropical climate conditions.}, language = {en} } @article{NievasPilzPrehnetal.2022, author = {Nievas, Cecilia and Pilz, Marco and Prehn, Karsten and Schorlemmer, Danijel and Weatherill, Graeme and Cotton, Fabrice}, title = {Calculating earthquake damage building by building}, series = {Bulletin of earthquake engineering : official publication of the European Association for Earthquake Engineering}, volume = {20}, journal = {Bulletin of earthquake engineering : official publication of the European Association for Earthquake Engineering}, number = {3}, publisher = {Springer}, address = {Dordrecht}, issn = {1570-761X}, doi = {10.1007/s10518-021-01303-w}, pages = {1519 -- 1565}, year = {2022}, abstract = {The creation of building exposure models for seismic risk assessment is frequently challenging due to the lack of availability of detailed information on building structures. Different strategies have been developed in recent years to overcome this, including the use of census data, remote sensing imagery and volunteered graphic information (VGI). This paper presents the development of a building-by-building exposure model based exclusively on openly available datasets, including both VGI and census statistics, which are defined at different levels of spatial resolution and for different moments in time. The initial model stemming purely from building-level data is enriched with statistics aggregated at the neighbourhood and city level by means of a Monte Carlo simulation that enables the generation of full realisations of damage estimates when using the exposure model in the context of an earthquake scenario calculation. Though applicable to any other region of interest where analogous datasets are available, the workflow and approach followed are explained by focusing on the case of the German city of Cologne, for which a scenario earthquake is defined and the potential damage is calculated. The resulting exposure model and damage estimates are presented, and it is shown that the latter are broadly consistent with damage data from the 1978 Albstadt earthquake, notwithstanding the differences in the scenario. Through this real-world application we demonstrate the potential of VGI and open data to be used for exposure modelling for natural risk assessment, when combined with suitable knowledge on building fragility and accounting for the inherent uncertainties.}, language = {en} } @phdthesis{Niemz2022, author = {Niemz, Peter}, title = {Imaging and modeling of hydraulic fractures in crystalline rock via induced seismic activity}, doi = {10.25932/publishup-55659}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-556593}, school = {Universit{\"a}t Potsdam}, pages = {135}, year = {2022}, abstract = {Enhanced geothermal systems (EGS) are considered a cornerstone of future sustainable energy production. In such systems, high-pressure fluid injections break the rock to provide pathways for water to circulate in and heat up. This approach inherently induces small seismic events that, in rare cases, are felt or can even cause damage. Controlling and reducing the seismic impact of EGS is crucial for a broader public acceptance. To evaluate the applicability of hydraulic fracturing (HF) in EGS and to improve the understanding of fracturing processes and the hydromechanical relation to induced seismicity, six in-situ, meter-scale HF experiments with different injection schemes were performed under controlled conditions in crystalline rock in a depth of 410 m at the {\"A}sp{\"o} Hard Rock Laboratory (Sweden). I developed a semi-automated, full-waveform-based detection, classification, and location workflow to extract and characterize the acoustic emission (AE) activity from the continuous recordings of 11 piezoelectric AE sensors. Based on the resulting catalog of 20,000 AEs, with rupture sizes of cm to dm, I mapped and characterized the fracture growth in great detail. The injection using a novel cyclic injection scheme (HF3) had a lower seismic impact than the conventional injections. HF3 induced fewer AEs with a reduced maximum magnitude and significantly larger b-values, implying a decreased number of large events relative to the number of small ones. Furthermore, HF3 showed an increased fracture complexity with multiple fractures or a fracture network. In contrast, the conventional injections developed single, planar fracture zones (Publication 1). An independent, complementary approach based on a comparison of modeled and observed tilt exploits transient long-period signals recorded at the horizontal components of two broad-band seismometers a few tens of meters apart from the injections. It validated the efficient creation of hydraulic fractures and verified the AE-based fracture geometries. The innovative joint analysis of AEs and tilt signals revealed different phases of the fracturing process, including the (re-)opening, growth, and aftergrowth of fractures, and provided evidence for the reactivation of a preexisting fault in one of the experiments (Publication 2). A newly developed network-based waveform-similarity analysis applied to the massive AE activity supports the latter finding. To validate whether the reduction of the seismic impact as observed for the cyclic injection schemes during the {\"A}sp{\"o} mine-scale experiments is transferable to other scales, I additionally calculated energy budgets for injection experiments from previously conducted laboratory tests and from a field application. Across all three scales, the cyclic injections reduce the seismic impact, as depicted by smaller maximum magnitudes, larger b-values, and decreased injection efficiencies (Publication 3).}, language = {en} } @article{NguyenKumarMusolffetal.2022, author = {Nguyen, Tam and Kumar, Rohini and Musolff, Andreas and Lutz, Stefanie R. and Sarrazin, Fanny and Attinger, Sabine and Fleckenstein, Jan H.}, title = {Disparate Seasonal Nitrate Export From Nested Heterogeneous Subcatchments Revealed With StorAge Selection Functions}, series = {Water resources research}, volume = {58}, journal = {Water resources research}, number = {3}, publisher = {American Geophysical Union}, address = {Washington}, issn = {0043-1397}, doi = {10.1029/2021WR030797}, pages = {20}, year = {2022}, abstract = {Understanding catchment controls on catchment solute export is a prerequisite for water quality management. StorAge Selection (SAS) functions encapsulate essential information about catchment functioning in terms of discharge selection preference and solute export dynamics. However, they lack information on the spatial origin of solutes when applied at the catchment scale, thereby limiting our understanding of the internal (subcatchment) functioning. Here, we parameterized SAS functions in a spatially explicit way to understand the internal catchment responses and transport dynamics of reactive dissolved nitrate (N-NO3). The model was applied in a nested mesoscale catchment (457 km(2)), consisting of a mountainous partly forested, partly agricultural subcatchment, a middle-reach forested subcatchment, and a lowland agricultural subcatchment. The model captured flow and nitrate concentration dynamics not only at the catchment outlet but also at internal gauging stations. Results reveal disparate subsurface mixing dynamics and nitrate export among headwater and lowland subcatchments. The headwater subcatchment has high seasonal variation in subsurface mixing schemes and younger water in discharge, while the lowland subcatchment has less pronounced seasonality in subsurface mixing and much older water in discharge. Consequently, nitrate concentration in discharge from the headwater subcatchment shows strong seasonality, whereas that from the lowland subcatchment is stable in time. The temporally varying responses of headwater and lowland subcatchments alternate the dominant contribution to nitrate export in high and low-flow periods between subcatchments. Overall, our results demonstrate that the spatially explicit SAS modeling provides useful information about internal catchment functioning, helping to develop or evaluate spatial management practices.}, language = {en} } @article{NeuharthBruneWronaetal.2022, author = {Neuharth, Derek and Brune, Sascha and Wrona, Thilo and Glerum, Anne and Braun, Jean and Yuan, Xiaoping}, title = {Evolution of rift systems and their fault networks in response to surface processes}, series = {Tectonics}, volume = {41}, journal = {Tectonics}, number = {3}, publisher = {American Geophysical Union}, address = {Washington}, issn = {0278-7407}, doi = {10.1029/2021TC007166}, pages = {22}, year = {2022}, abstract = {Continental rifting is responsible for the generation of major sedimentary basins, both during rift inception and during the formation of rifted continental margins. Geophysical and field studies revealed that rifts feature complex networks of normal faults but the factors controlling fault network properties and their evolution are still matter of debate. Here, we employ high-resolution 2D geodynamic models (ASPECT) including two-way coupling to a surface processes (SP) code (FastScape) to conduct 12 models of major rift types that are exposed to various degrees of erosion and sedimentation. We further present a novel quantitative fault analysis toolbox (Fatbox), which allows us to isolate fault growth patterns, the number of faults, and their length and displacement throughout rift history. Our analysis reveals that rift fault networks may evolve through five major phases: (a) distributed deformation and coalescence, (b) fault system growth, (c) fault system decline and basinward localization, (d) rift migration, and (e) breakup. These phases can be correlated to distinct rifted margin domains. Models of asymmetric rifting suggest rift migration is facilitated through both ductile and brittle deformation within a weak exhumation channel that rotates subhorizontally and remains active at low angles. In sedimentation-starved settings, this channel satisfies the conditions for serpentinization. We find that SP are not only able to enhance strain localization and to increase fault longevity but that they also reduce the total length of the fault system, prolong rift phases and delay continental breakup.}, language = {en} } @phdthesis{Neuharth2022, author = {Neuharth, Derek}, title = {Evolution of divergent and strike-slip boundaries in response to surface processes}, doi = {10.25932/publishup-54940}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-549403}, school = {Universit{\"a}t Potsdam}, pages = {xiii, 108}, year = {2022}, abstract = {Plate tectonics describes the movement of rigid plates at the surface of the Earth as well as their complex deformation at three types of plate boundaries: 1) divergent boundaries such as rift zones and mid-ocean ridges, 2) strike-slip boundaries where plates grind past each other, such as the San Andreas Fault, and 3) convergent boundaries that form large mountain ranges like the Andes. The generally narrow deformation zones that bound the plates exhibit complex strain patterns that evolve through time. During this evolution, plate boundary deformation is driven by tectonic forces arising from Earth's deep interior and from within the lithosphere, but also by surface processes, which erode topographic highs and deposit the resulting sediment into regions of low elevation. Through the combination of these factors, the surface of the Earth evolves in a highly dynamic way with several feedback mechanisms. At divergent boundaries, for example, tensional stresses thin the lithosphere, forcing uplift and subsequent erosion of rift flanks, which creates a sediment source. Meanwhile, the rift center subsides and becomes a topographic low where sediments accumulate. This mass transfer from foot- to hanging wall plays an important role during rifting, as it prolongs the activity of individual normal faults. When rifting continues, continents are eventually split apart, exhuming Earth's mantle and creating new oceanic crust. Because of the complex interplay between deep tectonic forces that shape plate boundaries and mass redistribution at the Earth's surface, it is vital to understand feedbacks between the two domains and how they shape our planet. In this study I aim to provide insight on two primary questions: 1) How do divergent and strike-slip plate boundaries evolve? 2) How is this evolution, on a large temporal scale and a smaller structural scale, affected by the alteration of the surface through erosion and deposition? This is done in three chapters that examine the evolution of divergent and strike-slip plate boundaries using numerical models. Chapter 2 takes a detailed look at the evolution of rift systems using two-dimensional models. Specifically, I extract faults from a range of rift models and correlate them through time to examine how fault networks evolve in space and time. By implementing a two-way coupling between the geodynamic code ASPECT and landscape evolution code FastScape, I investigate how the fault network and rift evolution are influenced by the system's erosional efficiency, which represents many factors like lithology or climate. In Chapter 3, I examine rift evolution from a three-dimensional perspective. In this chapter I study linkage modes for offset rifts to determine when fast-rotating plate-boundary structures known as continental microplates form. Chapter 4 uses the two-way numerical coupling between tectonics and landscape evolution to investigate how a strike-slip boundary responds to large sediment loads, and whether this is sufficient to form an entirely new type of flexural strike-slip basin.}, language = {en} } @article{NagakuraSchubertWagneretal.2022, author = {Nagakura, Toshiki and Schubert, Florian and Wagner, Dirk and Kallmeyer, Jens}, title = {Biological sulfate reduction in deep subseafloor sediment of Guaymas Basin}, series = {Frontiers in microbiology}, volume = {13}, journal = {Frontiers in microbiology}, publisher = {Frontiers Media}, address = {Lausanne}, organization = {IODP Exp 385 Shipboard Sci Party}, issn = {1664-302X}, doi = {10.3389/fmicb.2022.845250}, pages = {12}, year = {2022}, abstract = {Sulfate reduction is the quantitatively most important process to degrade organic matter in anoxic marine sediment and has been studied intensively in a variety of settings. Guaymas Basin, a young marginal ocean basin, offers the unique opportunity to study sulfate reduction in an environment characterized by organic-rich sediment, high sedimentation rates, and high geothermal gradients (100-958 degrees C km(-1)). We measured sulfate reduction rates (SRR) in samples taken during the International Ocean Discovery Program (IODP) Expedition 385 using incubation experiments with radiolabeled (SO42-)-S-35 carried out at in situ pressure and temperature. The highest SRR (387 nmol cm(-3) d(-1)) was recorded in near-surface sediments from Site U1548C, which had the steepest geothermal gradient (958 degrees C km(-1)). At this site, SRR were generally over an order of magnitude higher than at similar depths at other sites (e.g., 387-157 nmol cm(-3) d(-1) at 1.9 mbsf from Site U1548C vs. 46-1.0 nmol cm(-3) d(-1) at 2.1 mbsf from Site U1552B). Site U1546D is characterized by a sill intrusion, but it had already reached thermal equilibrium and SRR were in the same range as nearby Site U1545C, which is minimally affected by sills. The wide temperature range observed at each drill site suggests major shifts in microbial community composition with very different temperature optima but awaits confirmation by molecular biological analyses. At the transition between the mesophilic and thermophilic range around 40 degrees C-60 degrees C, sulfate-reducing activity appears to be decreased, particularly in more oligotrophic settings, but shows a slight recovery at higher temperatures.}, language = {en} } @article{MuellerNeugebauerBenDoretal.2022, author = {M{\"u}ller, Daniela and Neugebauer, Ina and Ben Dor, Yoav and Enzel, Yehouda and Schwab, Markus Julius and Tjallingii, Rik and Brauer, Achim}, title = {Phases of stability during major hydroclimate change ending the Last Glacial in the Levant}, series = {Scientific reports}, volume = {12}, journal = {Scientific reports}, number = {1}, publisher = {Macmillan Publishers Limited, part of Springer Nature}, address = {London}, issn = {2045-2322}, doi = {10.1038/s41598-022-10217-9}, pages = {12}, year = {2022}, abstract = {In-depth understanding of the reorganization of the hydrological cycle in response to global climate change is crucial in highly sensitive regions like the eastern Mediterranean, where water availability is a major factor for socioeconomic and political development. The sediments of Lake Lisan provide a unique record of hydroclimatic change during the last glacial to Holocene transition (ca. 24-11 ka) with its tremendous water level drop of similar to 240 m that finally led to its transition into the present hypersaline water body-the Dead Sea. Here we utilize high-resolution sedimentological analyses from the marginal terraces and deep lake to reconstruct an unprecedented seasonal record of the last millennia of Lake Lisan. Aragonite varve formation in intercalated intervals of our record demonstrates that a stepwise long-term lake level decline was interrupted by almost one millennium of rising or stable water level. Even periods of pronounced water level drops indicated by gypsum deposition were interrupted by decades of positive water budgets. Our results thus highlight that even during major climate change at the end of the last glacial, decadal to millennial periods of relatively stable or positive moisture supply occurred which could have been an important premise for human sedentism.}, language = {en} } @phdthesis{Mueller2022, author = {M{\"u}ller, Daniela}, title = {Abrupt climate changes and extreme events in two different varved lake sediment records}, doi = {10.25932/publishup-55833}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-558331}, school = {Universit{\"a}t Potsdam}, pages = {XVIII, 209}, year = {2022}, abstract = {Different lake systems might reflect different climate elements of climate changes, while the responses of lake systems are also divers, and are not completely understood so far. Therefore, a comparison of lakes in different climate zones, during the high-amplitude and abrupt climate fluctuations of the Last Glacial to Holocene transition provides an exceptional opportunity to investigate distinct natural lake system responses to different abrupt climate changes. The aim of this doctoral thesis was to reconstruct climatic and environmental fluctuations down to (sub-) annual resolution from two different lake systems during the Last Glacial-Interglacial transition (~17 and 11 ka). Lake Gościąż, situated in the temperate central Poland, developed in the Aller{\o}d after recession of the Last Glacial ice sheets. The Dead Sea is located in the Levant (eastern Mediterranean) within a steep gradient from sub-humid to hyper-arid climate, and formed in the mid-Miocene. Despite their differences in sedimentation processes, both lakes form annual laminations (varves), which are crucial for studies of abrupt climate fluctuations. This doctoral thesis was carried out within the DFG project PALEX-II (Paleohydrology and Extreme Floods from the Dead Sea ICDP Core) that investigates extreme hydro-meteorological events in the ICDP core in relation to climate changes, and ICLEA (Virtual Institute of Integrated Climate and Landscape Evolution Analyses) that intends to better the understanding of climate dynamics and landscape evolutions in north-central Europe since the Last Glacial. Further, it contributes to the Helmholtz Climate Initiative REKLIM (Regional Climate Change and Humans) Research Theme 3 "Extreme events across temporal and spatial scales" that investigates extreme events using climate data, paleo-records and model-based simulations. The three main aims were to (1) establish robust chronologies of the lakes, (2) investigate how major and abrupt climate changes affect the lake systems, and (3) to compare the responses of the two varved lakes to these hemispheric-scale climate changes. Robust chronologies are a prerequisite for high-resolved climate and environmental reconstructions, as well as for archive comparisons. Thus, addressing the first aim, the novel chronology of Lake Gościąż was established by microscopic varve counting and Bayesian age-depth modelling in Bacon for a non-varved section, and was corroborated by independent age constrains from 137Cs activity concentration measurements, AMS radiocarbon dating and pollen analysis. The varve chronology reaches from the late Aller{\o}d until AD 2015, revealing more Holocene varves than a previous study of Lake Gościąż suggested. Varve formation throughout the complete Younger Dryas (YD) even allowed the identification of annually- to decadal-resolved leads and lags in proxy responses at the YD transitions. The lateglacial chronology of the Dead Sea (DS) was thus far mainly based on radiocarbon and U/Th-dating. In the unique ICDP core from the deep lake centre, continuous search for cryptotephra has been carried out in lateglacial sediments between two prominent gypsum deposits - the Upper and Additional Gypsum Units (UGU and AGU, respectively). Two cryptotephras were identified with glass analyses that correlate with tephra deposits from the S{\"u}phan and Nemrut volcanoes indicating that the AGU is ~1000 years younger than previously assumed, shifting it into the YD, and the underlying varved interval into the B{\o}lling/Aller{\o}d, contradicting previous assumptions. Using microfacies analyses, stable isotopes and temperature reconstructions, the second aim was achieved at Lake Gościąż. The YD lake system was dynamic, characterized by higher aquatic bioproductivity, more re-suspended material and less anoxia than during the Aller{\o}d and Early Holocene, mainly influenced by stronger water circulation and catchment erosion due to stronger westerly winds and less lake sheltering. Cooling at the YD onset was ~100 years longer than the final warming, while environmental proxies lagged the onset of cooling by ~90 years, but occurred contemporaneously during the termination of the YD. Chironomid-based temperature reconstructions support recent studies indicating mild YD summer temperatures. Such a comparison of annually-resolved proxy responses to both abrupt YD transitions is rare, because most European lake archives do not preserve varves during the YD. To accomplish the second aim at the DS, microfacies analyses were performed between the UGU (~17 ka) and Holocene onset (~11 ka) in shallow- (Masada) and deep-water (ICDP core) environments. This time interval is marked by a huge but fluctuating lake level drop and therefore the complete transition into the Holocene is only recorded in the deep-basin ICDP core. In this thesis, this transition was investigated for the first time continuously and in detail. The final two pronounced lake level drops recorded by deposition of the UGU and AGU, were interrupted by one millennium of relative depositional stability and a positive water budget as recorded by aragonite varve deposition interrupted by only a few event layers. Further, intercalation of aragonite varves between the gypsum beds of the UGU and AGU shows that these generally dry intervals were also marked by decadal- to centennial-long rises in lake level. While continuous aragonite varves indicate decadal-long stable phases, the occurrence of thicker and more frequent event layers suggests general more instability during the gypsum units. These results suggest a pattern of complex and variable hydroclimate at different time scales during the Lateglacial at the DS. The third aim was accomplished based on the individual studies above that jointly provide an integrated picture of different lake responses to different climate elements of hemispheric-scale abrupt climate changes during the Last Glacial-Interglacial transition. In general, climatically-driven facies changes are more dramatic in the DS than at Lake Gościąż. Further, Lake Gościąż is characterized by continuous varve formation nearly throughout the complete profile, whereas the DS record is widely characterized by extreme event layers, hampering the establishment of a continuous varve chronology. The lateglacial sedimentation in Lake Gościąż is mainly influenced by westerly winds and minor by changes in catchment vegetation, whereas the DS is primarily influenced by changes in winter precipitation, which are caused by temperature variations in the Mediterranean. Interestingly, sedimentation in both archives is more stable during the B{\o}lling/Aller{\o}d and more dynamic during the YD, even when sedimentation processes are different. In summary, this doctoral thesis presents seasonally-resolved records from two lake archives during the Lateglacial (ca 17-11 ka) to investigate the impact of abrupt climate changes in different lake systems. New age constrains from the identification of volcanic glass shards in the lateglacial sediments of the DS allowed the first lithology-based interpretation of the YD in the DS record and its comparison to Lake Gościąż. This highlights the importance of the construction of a robust chronology, and provides a first step for synchronization of the DS with other eastern Mediterranean archives. Further, climate reconstructions from the lake sediments showed variability on different time scales in the different archives, i.e. decadal- to millennial fluctuations in the lateglacial DS, and even annual variations and sub-decadal leads and lags in proxy responses during the rapid YD transitions in Lake Gościąż. This showed the importance of a comparison of different lake archives to better understand the regional and local impacts of hemispheric-scale climate variability. An unprecedented example is demonstrated here of how different lake systems show different lake responses and also react to different climate elements of abrupt climate changes. This further highlights the importance of the understanding of the respective lake system for climate reconstructions.}, language = {en} } @article{MukherjeeAdhikariNicolietal.2022, author = {Mukherjee, Shreya and Adhikari, Avishek and Nicoli, Gautier and Vadlamani, Ravikant}, title = {Neoarchean (similar to 2.73-2.70 Ga) accretionary history of the eastern Dharwar Craton, India}, series = {Precambrian research}, volume = {375}, journal = {Precambrian research}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0301-9268}, doi = {10.1016/j.precamres.2022.106657}, pages = {23}, year = {2022}, abstract = {Cratonic mid-crustal plutons may contain supracrustal enclaves that preserve evidence of an earlier growth history. The Eastern Dharwar craton records Neoarchean two-stage accretionary sequential growth (2.70 and 2.55 Ga) and a chronology of their enclaves could refine orogenic models. To test whether the metamorphic history of their enclaves was related to any of these stages, phase equilibria modelling and combined Lu-Hf and Sm-Nd geochronology on garnet were conducted on metapsammite, now preserved as garnet-orthopyroxene-cordierite gneiss. Phase equilibria modelling indicates peak metamorphic conditions, similar to 850 degrees C and similar to 8.5 kbar (M1a), were followed by near isothermal decompression to 5-6 kbar (M1b) and isobaric cooling to similar to 800 degrees C (M1c). The thermobaric gradient related to peak metamorphic conditions, similar to 30 degrees C kbar(-1), is typical of collisional orogens. Regression of the whole-rock and garnet, for sample S17b, yield Lu-Hf isochron ages of 2733 +/- 29 Ma, and for sample S18, 2724 +/- 13 Ma. A Lu-Hf weighted mean age for the porphyroblastic garnet suggests growth at 2725.5 +/- 11.9 Ma during the M1a-M1b stages. In contrast, the whole-rock sample S17b and the garnet fractions yield a Sm-Nd isochron age of 2696 +/- 10 Ma. From sample S18 the whole rock, garnet fractions, and orthopyroxene yield an isochron age of 2683 +/- 15 Ma. The garnet Sm-Nd weighted mean age at 2692.0 +/- 8.3 Ma constrains the M1b-M1c stages. We suggest that the protoliths to these supracrustal enclaves were deposited in an arc tectonic setting and underwent thickening followed by heating during peeled-back lithospheric convergence. Therefore, the earliest of the craton-forming accretionary stages is preserved as the similar to 2.73 Ga granulite-facies enclaves, marginally older than the 2.70-2.65 Ga cratonic greenstone volcanism. Tectonic exhumation of these mid-crustal granulite enclaves was in response to the late-Proterozoic (similar to 1.7 Ga) Bhopalpatnam orogeny.}, language = {en} } @article{MonhonvalStraussThomasetal.2022, author = {Monhonval, Arthur and Strauss, Jens and Thomas, Maxime and Hirst, Catherine and Titeux, Hugues and Louis, Justin and Gilliot, Alexia and D'Aische, Eleonore du Bois and Pereira, Benoit and Vandeuren, Aubry and Grosse, Guido and Schirrmeister, Lutz and Jongejans, Loeka Laura and Ulrich, Mathias and Opfergelt, Sophie}, title = {Thermokarst processes increase the supply of stabilizing surfaces and elements (Fe, Mn, Al, and Ca) for mineral-organic carbon interactions}, series = {Permafrost and periglacial processes}, volume = {33}, journal = {Permafrost and periglacial processes}, number = {4}, publisher = {Wiley}, address = {Hoboken}, issn = {1045-6740}, doi = {10.1002/ppp.2162}, pages = {452 -- 469}, year = {2022}, abstract = {The stabilizing properties of mineral-organic carbon (OC) interactions have been studied in many soil environments (temperate soils, podzol lateritic soils, and paddy soils). Recently, interest in their role in permafrost regions is increasing as permafrost was identified as a hotspot of change. In thawing ice-rich permafrost regions, such as the Yedoma domain, 327-466 Gt of frozen OC is buried in deep sediments. Interactions between minerals and OC are important because OC is located very near the mineral matrix. Mineral surfaces and elements could mitigate recent and future greenhouse gas emissions through physical and/or physicochemical protection of OC. The dynamic changes in redox and pH conditions associated with thermokarst lake formation and drainage trigger metal-oxide dissolution and precipitation, likely influencing OC stabilization and microbial mineralization. However, the influence of thermokarst processes on mineral-OC interactions remains poorly constrained. In this study, we aim to characterize Fe, Mn, Al, and Ca minerals and their potential protective role for OC. Total and selective extractions were used to assess the crystalline and amorphous oxides or complexed metal pools as well as the organic acids found within these pools. We analyzed four sediment cores from an ice-rich permafrost area in Central Yakutia, which were drilled (i) in undisturbed Yedoma uplands, (ii) beneath a recent lake formed within Yedoma deposits, (iii) in a drained thermokarst lake basin, and (iv) beneath a mature thermokarst lake from the early Holocene period. We find a decrease in the amount of reactive Fe, Mn, Al, and Ca in the deposits on lake formation (promoting reduction reactions), and this was largely balanced by an increase in the amount of reactive metals in the deposits on lake drainage (promoting oxidation reactions). We demonstrate an increase in the metal to C molar ratio on thermokarst process, which may indicate an increase in metal-C bindings and could provide a higher protective role against microbial mineralization of organic matter. Finally, we find that an increase in mineral-OC interactions corresponded to a decrease in CO2 and CH4 gas emissions on thermokarst process. Mineral-OC interactions could mitigate greenhouse gas production from permafrost thaw as soon as lake drainage occurs.}, language = {en} } @article{MitzscherlingMacLeanLipusetal.2022, author = {Mitzscherling, Julia and MacLean, Joana and Lipus, Daniel and Bartholom{\"a}us, Alexander and Mangelsdorf, Kai and Lipski, Andr{\´e} and Roddatis, Vladimir and Liebner, Susanne and Wagner, Dirk}, title = {Nocardioides alcanivorans sp. nov., a novel hexadecane-degrading species isolated from plastic waste}, series = {International journal of systematic and evolutionary microbiology}, volume = {72}, journal = {International journal of systematic and evolutionary microbiology}, number = {4}, publisher = {Microbiology Society}, address = {London}, issn = {1466-5026}, doi = {10.1099/ijsem.0.005319}, pages = {11}, year = {2022}, abstract = {Strain NGK65(T), a novel hexadecane degrading, non-motile, Gram-positive, rod-to-coccus shaped, aerobic bacterium, was isolated from plastic polluted soil sampled at a landfill. Strain NGK65(T) hydrolysed casein, gelatin, urea and was catalase-positive. It optimally grew at 28 degrees C. in 0-1\% NaCl and at pH 7.5-8.0. Glycerol, D-glucose, arbutin, aesculin, salicin, potassium 5-ketogluconate. sucrose, acetate, pyruvate and hexadecane were used as sole carbon sources. The predominant membrane fatty acids were iso-C-16:0 followed by iso-C(17:)0 and C-18:1 omega 9c. The major polar lipids were phosphatidylglycerol, phosphatidylethanolamine, phosphatidylinositol and hydroxyphosphatidylinositol. The cell-wall peptidoglycan type was A3 gamma, with LL-diaminopimelic acid and glycine as the diagnostic amino acids. MK 8 (H-4) was the predominant menaquinone. Phylogenetic analysis based on 16S rRNA gene sequences indicated that strain NGK65(T) belongs to the genus Nocardioides (phylum Actinobacteria). appearing most closely related to Nocardioides daejeonensis MJ31(T) (98.6\%) and Nocardioides dubius KSL-104(T) (98.3\%). The genomic DNA G+C content of strain NGK65(T) was 68.2\%. Strain NGK65(T) and the type strains of species involved in the analysis had average nucleotide identity values of 78.3-71.9\% as well as digital DNA-DNA hybridization values between 22.5 and 19.7\%, which clearly indicated that the isolate represents a novel species within the genus Nocardioides. Based on phenotypic and molecular characterization, strain NGK65(T) can clearly be differentiated from its phylogenetic neighbours to establish a novel species, for which the name Nocardioides alcanivorans sp. nov. is proposed. The type strain is NGK65(T) (=DSM 113112(T)=NCCB 100846(T)).}, language = {en} } @article{MerzBassoFischeretal.2022, author = {Merz, Bruno and Basso, Stefano and Fischer, Svenja and Lun, David and Bloeschl, Guenter and Merz, Ralf and Guse, Bjorn and Viglione, Alberto and Vorogushyn, Sergiy and Macdonald, Elena and Wietzke, Luzie and Schumann, Andreas}, title = {Understanding heavy tails of flood peak distributions}, series = {Water resources research}, volume = {58}, journal = {Water resources research}, number = {6}, publisher = {American Geophysical Union}, address = {Washington}, issn = {0043-1397}, doi = {10.1029/2021WR030506}, pages = {37}, year = {2022}, abstract = {Statistical distributions of flood peak discharge often show heavy tail behavior, that is, extreme floods are more likely to occur than would be predicted by commonly used distributions that have exponential asymptotic behavior. This heavy tail behavior may surprise flood managers and citizens, as human intuition tends to expect light tail behavior, and the heaviness of the tails is very difficult to predict, which may lead to unnecessarily high flood damage. Despite its high importance, the literature on the heavy tail behavior of flood distributions is rather fragmented. In this review, we provide a coherent overview of the processes causing heavy flood tails and the implications for science and practice. Specifically, we propose nine hypotheses on the mechanisms causing heavy tails in flood peak distributions related to processes in the atmosphere, the catchment, and the river system. We then discuss to which extent the current knowledge supports or contradicts these hypotheses. We also discuss the statistical conditions for the emergence of heavy tail behavior based on derived distribution theory and relate them to the hypotheses and flood generation mechanisms. We review the degree to which the heaviness of the tails can be predicted from process knowledge and data. Finally, we recommend further research toward testing the hypotheses and improving the prediction of heavy tails.}, language = {en} } @article{MelchertWischhoeferKnoblauchetal.2022, author = {Melchert, Jan Olaf and Wischh{\"o}fer, Philipp and Knoblauch, Christian and Eckhardt, Tim and Liebner, Susanne and Rethemeyer, Janet}, title = {Sources of CO2 Produced in Freshly Thawed Pleistocene-Age Yedoma Permafrost}, series = {Frontiers in Earth Science}, volume = {9}, journal = {Frontiers in Earth Science}, publisher = {Frontiers Media}, address = {Lausanne}, issn = {2296-6463}, doi = {10.3389/feart.2021.737237}, pages = {13}, year = {2022}, abstract = {The release of greenhouse gases from the large organic carbon stock in permafrost deposits in the circumarctic regions may accelerate global warming upon thaw. The extent of this positive climate feedback is thought to be largely controlled by the microbial degradability of the organic matter preserved in these sediments. In addition, weathering and oxidation processes may release inorganic carbon preserved in permafrost sediments as CO2, which is generally not accounted for. We used C-13 and C-14 analysis and isotopic mass balances to differentiate and quantify organic and inorganic carbon released as CO2 in the field from an active retrogressive thaw slump of Pleistocene-age Yedoma and during a 1.5-years incubation experiment. The results reveal that the dominant source of the CO2 released from freshly thawed Yedoma exposed as thaw mound is Pleistocene-age organic matter (48-80\%) and to a lesser extent modern organic substrate (3-34\%). A significant portion of the CO2 originated from inorganic carbon in the Yedoma (17-26\%). The mixing of young, active layer material with Yedoma at a site on the slump floor led to the preferential mineralization of this young organic carbon source. Admixtures of younger organic substrates in the Yedoma thaw mound were small and thus rapidly consumed as shown by lower contributions to the CO2 produced during few weeks of aerobic incubation at 4 degrees C corresponding to approximately one thaw season. Future CO2 fluxes from the freshly thawed Yedoma will contain higher proportions of ancient inorganic (22\%) and organic carbon (61-78\%) as suggested by the results at the end, after 1.5 years of incubation. The increasing contribution of inorganic carbon during the incubation is favored by the accumulation of organic acids from microbial organic matter degradation resulting in lower pH values and, in consequence, in inorganic carbon dissolution. Because part of the inorganic carbon pool is assumed to be of pedogenic origin, these emissions would ultimately not alter carbon budgets. The results of this study highlight the preferential degradation of younger organic substrates in freshly thawed Yedoma, if available, and a substantial release of CO2 from inorganic sources.}, language = {en} } @article{McCoolCoddingVernonetal.2022, author = {McCool, Weston C. and Codding, Brian F. and Vernon, Kenneth B. and Wilson, Kurt M. and Yaworsky, Peter M. and Marwan, Norbert and Kennett, Douglas J.}, title = {Climate change-induced population pressure drives high rates of lethal violence in the Prehispanic central Andes}, series = {Proceedings of the National Academy of Sciences of the United States of America : PNAS}, volume = {119}, journal = {Proceedings of the National Academy of Sciences of the United States of America : PNAS}, number = {17}, publisher = {National Acad. of Sciences}, address = {Washington}, issn = {0027-8424}, doi = {10.1073/pnas.2117556119}, pages = {8}, year = {2022}, abstract = {Understanding the influence of climate change and population pressure on human conflict remains a critically important topic in the social sciences. Long-term records that evaluate these dynamics across multiple centuries and outside the range of modern climatic variation are especially capable of elucidating the relative effect of-and the interaction between-climate and demography. This is crucial given that climate change may structure population growth and carrying capacity, while both climate and population influence per capita resource availability. This study couples paleoclimatic and demographic data with osteological evaluations of lethal trauma from 149 directly accelerator mass spectrometry C-14-dated individuals from the Nasca highland region of Peru. Multiple local and supraregional precipitation proxies are combined with a summed probability distribution of 149 C-14 dates to estimate population dynamics during a 700-y study window. Counter to previous findings, our analysis reveals a precipitous increase in violent deaths associated with a period of productive and stable climate, but volatile population dynamics. We conclude that favorable local climate conditions fostered population growth that put pressure on the marginal and highly circumscribed resource base, resulting in violent resource competition that manifested in over 450 y of internecine warfare. These findings help support a general theory of intergroup violence, indicating that relative resource scarcity-whether driven by reduced resource abundance or increased competition-can lead to violence in subsistence societies when the outcome is lower per capita resource availability.}, language = {en} } @article{MarUngerWalderdorffetal.2022, author = {Mar, Kathleen A. and Unger, Charlotte and Walderdorff, Ludmila and Butler, Tim}, title = {Beyond CO2 equivalence}, series = {Environmental science \& policy}, volume = {134}, journal = {Environmental science \& policy}, publisher = {Elsevier}, address = {Oxford}, issn = {1462-9011}, doi = {10.1016/j.envsci.2022.03.027}, pages = {127 -- 136}, year = {2022}, abstract = {In this article we review the physical and chemical properties of methane (CH4) relevant to impacts on climate, ecosystems, and air pollution, and examine the extent to which this is reflected in climate and air pollution governance. Although CH4 is governed under the UNFCCC climate regime, its treatment there is limited to the ways in which it acts as a "CO2 equivalent" climate forcer on a 100-year time frame. The UNFCCC framework neglects the impacts that CH4 has on near-term climate, as well its impacts on human health and ecosystems, which are primarily mediated by methane's role as a precursor to tropospheric ozone. Frameworks for air quality governance generally address tropospheric ozone as a pollutant, but do not regulate CH4 itself. Methane's climate and air quality impacts, together with its alarming rise in atmospheric concentrations in recent years, make it clear that mitigation of CH4 emissions needs to be accelerated globally. We examine challenges and opportunities for further progress on CH4 mitigation within the international governance landscapes for climate change and air pollution.}, language = {en} } @article{MacdonaldMerzGuseetal.2022, author = {Macdonald, Elena and Merz, Bruno and Guse, Bj{\"o}rn and Wietzke, Luzie and Ullrich, Sophie and Kemter, Matthias and Ahrens, Bodo and Vorogushyn, Sergiy}, title = {Event and catchment controls of heavy tail behavior of floods}, series = {Water resources research}, volume = {58}, journal = {Water resources research}, number = {6}, publisher = {American Geophysical Union}, address = {Washington}, issn = {0043-1397}, doi = {10.1029/2021WR031260}, pages = {25}, year = {2022}, abstract = {In some catchments, the distribution of annual maximum streamflow shows heavy tail behavior, meaning the occurrence probability of extreme events is higher than if the upper tail decayed exponentially. Neglecting heavy tail behavior can lead to an underestimation of the likelihood of extreme floods and the associated risk. Partly contradictory results regarding the controls of heavy tail behavior exist in the literature and the knowledge is still very dispersed and limited. To better understand the drivers, we analyze the upper tail behavior and its controls for 480 catchments in Germany and Austria over a period of more than 50 years. The catchments span from quickly reacting mountain catchments to large lowland catchments, allowing for general conclusions. We compile a wide range of event and catchment characteristics and investigate their association with an indicator of the tail heaviness of flood distributions, namely the shape parameter of the GEV distribution. Following univariate analyses of these characteristics, along with an evaluation of different aggregations of event characteristics, multiple linear regression models, as well as random forests, are constructed. A novel slope indicator, which represents the relation between the return period of flood peaks and event characteristics, captures the controls of heavy tails best. Variables describing the catchment response are found to dominate the heavy tail behavior, followed by event precipitation, flood seasonality, and catchment size. The pre-event moisture state in a catchment has no relevant impact on the tail heaviness even though it does influence flood magnitudes.}, language = {en} } @article{LilienkampvonSpechtWeatherilletal.2022, author = {Lilienkamp, Henning and von Specht, Sebastian and Weatherill, Graeme and Caire, Giuseppe and Cotton, Fabrice}, title = {Ground-Motion modeling as an image processing task}, series = {Bulletin of the Seismological Society of America}, volume = {112}, journal = {Bulletin of the Seismological Society of America}, number = {3}, publisher = {Seismological Society of America}, address = {Albany}, issn = {0037-1106}, doi = {10.1785/0120220008}, pages = {1565 -- 1582}, year = {2022}, abstract = {We construct and examine the prototype of a deep learning-based ground-motion model (GMM) that is both fully data driven and nonergodic. We formulate ground-motion modeling as an image processing task, in which a specific type of neural network, the U-Net, relates continuous, horizontal maps of earthquake predictive parameters to sparse observations of a ground-motion intensity measure (IM). The processing of map-shaped data allows the natural incorporation of absolute earthquake source and observation site coordinates, and is, therefore, well suited to include site-, source-, and path-specific amplification effects in a nonergodic GMM. Data-driven interpolation of the IM between observation points is an inherent feature of the U-Net and requires no a priori assumptions. We evaluate our model using both a synthetic dataset and a subset of observations from the KiK-net strong motion network in the Kanto basin in Japan. We find that the U-Net model is capable of learning the magnitude???distance scaling, as well as site-, source-, and path-specific amplification effects from a strong motion dataset. The interpolation scheme is evaluated using a fivefold cross validation and is found to provide on average unbiased predictions. The magnitude???distance scaling as well as the site amplification of response spectral acceleration at a period of 1 s obtained for the Kanto basin are comparable to previous regional studies.}, language = {en} } @article{LiSpangenbergSchicksetal.2022, author = {Li, Zhen and Spangenberg, Erik and Schicks, Judith Maria and Kempka, Thomas}, title = {Numerical simulation of hydrate formation in the LArge-Scale Reservoir Simulator (LARS)}, series = {Energies : open-access journal of related scientific research, technology development and studies in policy and management}, volume = {15}, journal = {Energies : open-access journal of related scientific research, technology development and studies in policy and management}, number = {6}, publisher = {MDPI}, address = {Basel}, issn = {1996-1073}, doi = {10.3390/en15061974}, pages = {27}, year = {2022}, abstract = {The LArge-scale Reservoir Simulator (LARS) has been previously developed to study hydrate dissociation in hydrate-bearing systems under in-situ conditions. In the present study, a numerical framework of equations of state describing hydrate formation at equilibrium conditions has been elaborated and integrated with a numerical flow and transport simulator to investigate a multi-stage hydrate formation experiment undertaken in LARS. A verification of the implemented modeling framework has been carried out by benchmarking it against another established numerical code. Three-dimensional (3D) model calibration has been performed based on laboratory data available from temperature sensors, fluid sampling, and electrical resistivity tomography. The simulation results demonstrate that temperature profiles, spatial hydrate distribution, and bulk hydrate saturation are consistent with the observations. Furthermore, our numerical framework can be applied to calibrate geophysical measurements, optimize post-processing workflows for monitoring data, improve the design of hydrate formation experiments, and investigate the temporal evolution of sub-permafrost methane hydrate reservoirs.}, language = {en} } @article{LiSpangenbergSchicksetal.2022, author = {Li, Zhen and Spangenberg, Erik and Schicks, Judith Maria and Kempka, Thomas}, title = {Numerical Simulation of Coastal Sub-Permafrost Gas Hydrate Formation in the Mackenzie Delta, Canadian Arctic}, series = {Energies}, volume = {15}, journal = {Energies}, number = {14}, publisher = {MDPI}, address = {Basel}, issn = {1996-1073}, doi = {10.3390/en15144986}, pages = {25}, year = {2022}, abstract = {The Mackenzie Delta (MD) is a permafrost-bearing region along the coasts of the Canadian Arctic which exhibits high sub-permafrost gas hydrate (GH) reserves. The GH occurring at the Mallik site in the MD is dominated by thermogenic methane (CH4), which migrated from deep conventional hydrocarbon reservoirs, very likely through the present fault systems. Therefore, it is assumed that fluid flow transports dissolved CH4 upward and out of the deeper overpressurized reservoirs via the existing polygonal fault system and then forms the GH accumulations in the Kugmallit-Mackenzie Bay Sequences. We investigate the feasibility of this mechanism with a thermo-hydraulic-chemical numerical model, representing a cross section of the Mallik site. We present the first simulations that consider permafrost formation and thawing, as well as the formation of GH accumulations sourced from the upward migrating CH4-rich formation fluid. The simulation results show that temperature distribution, as well as the thickness and base of the ice-bearing permafrost are consistent with corresponding field observations. The primary driver for the spatial GH distribution is the permeability of the host sediments. Thus, the hypothesis on GH formation by dissolved CH4 originating from deeper geological reservoirs is successfully validated. Furthermore, our results demonstrate that the permafrost has been substantially heated to 0.8-1.3 degrees C, triggered by the global temperature increase of about 0.44 degrees C and further enhanced by the Arctic Amplification effect at the Mallik site from the early 1970s to the mid-2000s.}, language = {en} } @article{LiTianRudayaetal.2022, author = {Li, Wenjia and Tian, Fang and Rudaya, Natalya A. and Herzschuh, Ulrike and Cao, Xianyong}, title = {Pollen-based holocene thawing-history of permafrost in Northern Asia and its potential impacts on climate change}, series = {Frontiers in Ecology and Evolution}, volume = {10}, journal = {Frontiers in Ecology and Evolution}, publisher = {Frontiers Media}, address = {Lausanne}, issn = {2296-701X}, doi = {10.3389/fevo.2022.894471}, pages = {13}, year = {2022}, abstract = {As the recent permafrost thawing of northern Asia proceeds due to anthropogenic climate change, precise and detailed palaeoecological records from past warm periods are essential to anticipate the extent of future permafrost variations. Here, based on the modern relationship between permafrost and vegetation (represented by pollen assemblages), we trained a Random Forest model using pollen and permafrost data and verified its reliability to reconstruct the history of permafrost in northern Asia during the Holocene. An early Holocene (12-8 cal ka BP) strong thawing trend, a middle-to-late Holocene (8-2 cal ka BP) relatively slow thawing trend, and a late Holocene freezing trend of permafrost in northern Asia are consistent with climatic proxies such as summer solar radiation and Northern Hemisphere temperature. The extensive distribution of permafrost in northern Asia inhibited the spread of evergreen coniferous trees during the early Holocene warming and might have decelerated the enhancement of the East Asian summer monsoon (EASM) by altering hydrological processes and albedo. Based on these findings, we suggest that studies of the EASM should consider more the state of permafrost and vegetation in northern Asia, which are often overlooked and may have a profound impact on climate change in this region.}, language = {en} } @article{LescesenSrajBasarinetal.2022, author = {Lescesen, Igor and Sraj, Mojca and Basarin, Biljana and Pavic, Dragoslav and Mesaros, Minucer and Mudelsee, Manfred}, title = {Regional flood frequency analysis of the sava river in south-eastern Europe}, series = {Sustainability}, volume = {14}, journal = {Sustainability}, number = {15}, publisher = {MDPI}, address = {Basel}, issn = {2071-1050}, doi = {10.3390/su14159282}, pages = {19}, year = {2022}, abstract = {Regional flood frequency analysis (RFFA) is a powerful method for interrogating hydrological series since it combines observational time series from several sites within a region to estimate risk-relevant statistical parameters with higher accuracy than from single-site series. Since RFFA extreme value estimates depend on the shape of the selected distribution of the data-generating stochastic process, there is need for a suitable goodness-of-distributional-fit measure in order to optimally utilize given data. Here we present a novel, least-squares-based measure to select the optimal fit from a set of five distributions, namely Generalized Extreme Value (GEV), Generalized Logistic, Gumbel, Log-Normal Type III and Log-Pearson Type III. The fit metric is applied to annual maximum discharge series from six hydrological stations along the Sava River in South-eastern Europe, spanning the years 1961 to 2020. Results reveal that (1) the Sava River basin can be assessed as hydrologically homogeneous and (2) the GEV distribution provides typically the best fit. We offer hydrological-meteorological insights into the differences among the six stations. For the period studied, almost all stations exhibit statistically insignificant trends, which renders the conclusions about flood risk as relevant for hydrological sciences and the design of regional flood protection infrastructure.}, language = {en} } @article{LandisSaikinZhelavskayaetal.2022, author = {Landis, Daji August and Saikin, Anthony and Zhelavskaya, Irina and Drozdov, Alexander and Aseev, Nikita and Shprits, Yuri Y. and Pfitzer, Maximilian F. and Smirnov, Artem G.}, title = {NARX Neural Network Derivations of the Outer Boundary Radiation Belt Electron Flux}, series = {Space Weather: the international journal of research and applications}, volume = {20}, journal = {Space Weather: the international journal of research and applications}, number = {5}, publisher = {American Geophysical Union}, address = {Washington}, issn = {1542-7390}, doi = {10.1029/2021SW002774}, pages = {18}, year = {2022}, abstract = {We present two new empirical models of radiation belt electron flux at geostationary orbit. GOES-15 measurements of 0.8 MeV electrons were used to train a Nonlinear Autoregressive with Exogenous input (NARX) neural network for both modeling GOES-15 flux values and an upper boundary condition scaling factor (BF). The GOES-15 flux model utilizes an input and feedback delay of 2 and 2 time steps (i.e., 5 min time steps) with the most efficient number of hidden layers set to 10. Magnetic local time, Dst, Kp, solar wind dynamic pressure, AE, and solar wind velocity were found to perform as predicative indicators of GOES-15 flux and therefore were used as the exogenous inputs. The NARX-derived upper boundary condition scaling factor was used in conjunction with the Versatile Electron Radiation Belt (VERB) code to produce reconstructions of the radiation belts during the period of July-November 1990, independent of in-situ observations. Here, Kp was chosen as the sole exogenous input to be more compatible with the VERB code. This Combined Release and Radiation Effects Satellite-era reconstruction showcases the potential to use these neural network-derived boundary conditions as a method of hindcasting the historical radiation belts. This study serves as a companion paper to another recently published study on reconstructing the radiation belts during Solar Cycles 17-24 (Saikin et al., 2021, ), for which the results featured in this paper were used.}, language = {en} } @article{KuehnHainzlDahmetal.2022, author = {K{\"u}hn, Daniela and Hainzl, Sebastian and Dahm, Torsten and Richter, Gudrun and Vera Rodriguez, Ismael}, title = {A review of source models to further the understanding of the seismicity of the Groningen field}, series = {Netherlands journal of geosciences : NJG}, volume = {101}, journal = {Netherlands journal of geosciences : NJG}, publisher = {Cambridge Univ. Press}, address = {Cambridge}, issn = {0016-7746}, doi = {10.1017/njg.2022.7}, pages = {12}, year = {2022}, abstract = {The occurrence of felt earthquakes due to gas production in Groningen has initiated numerous studies and model attempts to understand and quantify induced seismicity in this region. The whole bandwidth of available models spans the range from fully deterministic models to purely empirical and stochastic models. In this article, we summarise the most important model approaches, describing their main achievements and limitations. In addition, we discuss remaining open questions and potential future directions of development.}, language = {en} } @article{KumarGuntuAgarwaletal.2022, author = {Kumar, Satish and Guntu, Ravi Kumar and Agarwal, Ankit and Villuri, Vasant Govind Kumar and Pasupuleti, Srinivas and Kaushal, Deo Raj and Gosian, Ashwin Kumar and Bronstert, Axel}, title = {Multi-objective optimization for stormwater management by green-roofs and infiltration trenches to reduce urban flooding in central Delhi}, series = {Journal of hydrology}, volume = {606}, journal = {Journal of hydrology}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0022-1694}, doi = {10.1016/j.jhydrol.2022.127455}, pages = {16}, year = {2022}, abstract = {Urban surface runoff management via best management practices (BMP) and low impact development (LID) has earned significant recognition owing to positive environmental and ecological impacts. However, due to the complexity of the parameters involved, the estimation of LID efficiency in attenuating the urban surface runoff at the watershed scale is challenging. A planning analysis of employing Green Roofs and Infiltration Trenches as BMPs/LIDs practices for urban surface runoff control is presented in this study. A multi-objective optimization decision-making framework is established by coupling SWMM (Storm Water Management Model) with NSGA-II models to check the performance of BMPs/LIDs concerning the cost-benefit analysis of LID at the watershed scale. Two urbanized areas belonging to Central Delhi in India were used as case studies. The results showed that the SWMM model is useful in simulating optimization problems for managing urban surface runoff. The optimum scenarios efficiently minimized the urban runoff volume while maintaining the BMPs/LIDs implementation costs and size. With BMPs/LIDs implementation, the reduction in runoff volume increases as expenses increase initially; however, there is no noticeable reduction in flood volume after a certain threshold. Contrasted with the haphazard arrangement of BMPs/LIDs, the proposed approach demonstrates 22\%-24\% runoff reductions for the same expenditures in watershed 1 and 23\%-26\% in watershed 2. The result of the study provides insights into planning and management of the urban surface runoff control with LID practices. The proposed framework assists the hydrologists in optimum selection and placements of BMPs/LIDs practices to acquire the most extreme ecological advantages with the least expenses.}, language = {en} } @article{KraemerGelbrechtPavithranetal.2022, author = {Kr{\"a}mer, Hauke Kai and Gelbrecht, Maximilian and Pavithran, Induja and Sujith, Ravindran and Marwan, Norbert}, title = {Optimal state space reconstruction via Monte Carlo decision tree search}, series = {Nonlinear Dynamics}, volume = {108}, journal = {Nonlinear Dynamics}, number = {2}, publisher = {Springer}, address = {Dordrecht}, issn = {0924-090X}, doi = {10.1007/s11071-022-07280-2}, pages = {1525 -- 1545}, year = {2022}, abstract = {A novel idea for an optimal time delay state space reconstruction from uni- and multivariate time series is presented. The entire embedding process is considered as a game, in which each move corresponds to an embedding cycle and is subject to an evaluation through an objective function. This way the embedding procedure can be modeled as a tree, in which each leaf holds a specific value of the objective function. By using a Monte Carlo ansatz, the proposed algorithm populates the tree with many leafs by computing different possible embedding paths and the final embedding is chosen as that particular path, which ends at the leaf with the lowest achieved value of the objective function. The method aims to prevent getting stuck in a local minimum of the objective function and can be used in a modular way, enabling practitioners to choose a statistic for possible delays in each embedding cycle as well as a suitable objective function themselves. The proposed method guarantees the optimization of the chosen objective function over the parameter space of the delay embedding as long as the tree is sampled sufficiently. As a proof of concept, we demonstrate the superiority of the proposed method over the classical time delay embedding methods using a variety of application examples. We compare recurrence plot-based statistics inferred from reconstructions of a Lorenz-96 system and highlight an improved forecast accuracy for map-like model data as well as for palaeoclimate isotope time series. Finally, we utilize state space reconstruction for the detection of causality and its strength between observables of a gas turbine type thermoacoustic combustor.}, language = {en} } @article{KrstulovićRosaFerreiraSanchezetal.2022, author = {Krstulović, Marija and Rosa, Angelika D. and Ferreira Sanchez, Dario and Libon, L{\´e}lia and Albers, Christian and Merkulova, Margarita and Grolimund, Daniel and Irifune, Tetsuo and Wilke, Max}, title = {Effect of temperature on the densification of silicate melts to lower earth's mantle conditions}, series = {Physics of the earth and planetary interiors}, volume = {323}, journal = {Physics of the earth and planetary interiors}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0031-9201}, doi = {10.1016/j.pepi.2021.106823}, pages = {13}, year = {2022}, abstract = {Physical properties of silicate melts play a key role for global planetary dynamics, controlling for example volcanic eruption styles, mantle convection and elemental cycling in the deep Earth. They are significantly modified by structural changes at the atomic scale due to external parameters such as pressure and temperature or due to chemistry. Structural rearrangements such as 4- to 6-fold coordination change of Si with increasing depth may profoundly influence melt properties, but have so far mostly been studied at ambient temperature due to experimental difficulties. In order to investigate the structural properties of silicate melts and their densification mechanisms at conditions relevant to the deep Earth's interior, we studied haplo basaltic glasses and melts (albite-diopside composition) at high pressure and temperature conditions in resistively and laser-heated diamond anvil cells using X-ray absorption near edge structure spectroscopy. Samples were doped with 10 wt\% of Ge, which is accessible with this experimental technique and which commonly serves as a structural analogue for the network forming cation Si. We acquired spectra on the Ge K edge up to 48 GPa and 5000 K and derived the average Ge-O coordination number NGe-O, and bond distance RGe-O as functions of pressure. Our results demonstrate a continuous transformation from tetrahedral to octahedral coordination between ca. 5 and 30 GPa at ambient temperature. Above 1600 K the data reveal a reduction of the pressure needed to complete conversion to octahedral coordination by ca. 30 \%. The results allow us to determine the influence of temperature on the Si coordination number changes in natural melts in the Earth's interior. We propose that the complete transition to octahedral coordination in basaltic melts is reached at about 40 GPa, corresponding to a depth of ca. 1200 km in the uppermost lower mantle. At the core-mantle boundary (2900 km, 130 GPa, 3000 K) the existence of non-buoyant melts has been proposed to explain observed low seismic wave velocity features. Our results highlight that the melt composition can affect the melt density at such extreme conditions and may strongly influence the structural response.}, language = {en} } @article{KothaWeatherillBindietal.2022, author = {Kotha, Sreeram Reddy and Weatherill, Graeme and Bindi, Dino and Cotton, Fabrice}, title = {Near-source magnitude scaling of spectral accelerations}, series = {Bulletin of earthquake engineering : official publication of the European Association for Earthquake Engineering}, volume = {20}, journal = {Bulletin of earthquake engineering : official publication of the European Association for Earthquake Engineering}, number = {3}, publisher = {Springer}, address = {Dordrecht}, issn = {1570-761X}, doi = {10.1007/s10518-021-01308-5}, pages = {1343 -- 1370}, year = {2022}, abstract = {Ground-motion models (GMMs) are often used to predict the random distribution of Spectral accelerations (SAs) at a site due to a nearby earthquake. In probabilistic seismic hazard and risk assessment, large earthquakes occurring close to a site are considered as critical scenarios. GMMs are expected to predict realistic SAs with low within-model uncertainty (sigma(mu)) for such rare scenarios. However, the datasets used to regress GMMs are usually deficient of data from critical scenarios. The (Kotha et al., A Regionally Adaptable Ground-Motion Model for Shallow Crustal Earthquakes in Europe Bulletin of Earthquake Engineering 18:4091-4125, 2020) GMM developed from the Engineering strong motion (ESM) dataset was found to predict decreasing short-period SAs with increasing M-W >= M-h = 6.2, and with large sigma(mu) at near-source distances <= 30km. In this study, we updated the parametrisation of the GMM based on analyses of ESM and the Near source strong motion (NESS) datasets. With M-h = 5.7, we could rectify the M-W scaling issue, while also reducing sigma(mu). at M-W >= M-h. We then evaluated the GMM against NESS data, and found that the SAs from a few large, thrust-faulting events in California, New Zealand, Japan, and Mexico are significantly higher than GMM median predictions. However, recordings from these events were mostly made on soft-soil geology, and contain anisotropic pulse-like effects. A more thorough non-ergodic treatment of NESS was not possible because most sites sampled unique events in very diverse tectonic environments. We provide an updated set of GMM coefficients,sigma(mu), and heteroscedastic variance models; while also cautioning against its application for M-W <= 4 in low-moderate seismicity regions without evaluating the homogeneity of M-W estimates between pan-European ESM and regional datasets.}, language = {en} } @article{KloseGuillemoteauVignolietal.2022, author = {Klose, Tim and Guillemoteau, Julien and Vignoli, Giulio and Tronicke, Jens}, title = {Laterally constrained inversion (LCI) of multi-configuration EMI data with tunable sharpness}, series = {Journal of applied geophysics}, volume = {196}, journal = {Journal of applied geophysics}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0926-9851}, doi = {10.1016/j.jappgeo.2021.104519}, pages = {9}, year = {2022}, abstract = {Frequency-domain electromagnetic (FDEM) data are commonly inverted to characterize subsurface geoelectrical properties using smoothness constraints in 1D inversion schemes assuming a layered medium. Smoothness constraints are suitable for imaging gradual transitions of subsurface geoelectrical properties caused, for example, by varying sand, clay, or fluid content. However, such inversion approaches are limited in characterizing sharp interfaces. Alternative regularizations based on the minimum gradient support (MGS) stabilizers can, instead, be used to promote results with different levels of smoothness/sharpness selected by simply acting on the so-called focusing parameter. The MGS regularization has been implemented for different kinds of geophysical data inversion strategies. However, concerning FDEM data, the MGS regularization has only been implemented for vertically constrained inversion (VCI) approaches but not for laterally constrained inversion (LCI) approaches. We present a novel LCI approach for FDEM data using the MGS regularization for the vertical and lateral direction. Using synthetic and field data examples, we demonstrate that our approach can efficiently and automatically provide a set of model solutions characterized by different levels of sharpness and variable lateral consistencies. In terms of data misfit, the obtained set of solutions contains equivalent models allowing us also to investigate the non-uniqueness of FDEM data inversion.}, language = {en} }