@article{KayaDupontNivetFrielingetal.2022, author = {Kaya, Mustafa Y{\"u}cel and Dupont-Nivet, Guillaume and Frieling, Joost and Fioroni, Chiara and Rohrmann, Alexander and Alt{\i}ner, Sevin{\c{c}} {\"O}zkan and Vardar, Ezgi and Tanyas, Hakan and Mamtimin, Mehmut and Zhaojie, Guo}, title = {The Eurasian epicontinental sea was an important carbon sink during the Palaeocene-Eocene thermal maximum}, series = {Communications earth and environment}, volume = {3}, journal = {Communications earth and environment}, number = {1}, publisher = {Springer Nature}, address = {London}, issn = {2662-4435}, doi = {10.1038/s43247-022-00451-4}, pages = {10}, year = {2022}, abstract = {The Palaeocene-Eocene Thermal Maximum (ca. 56 million years ago) offers a primary analogue for future global warming and carbon cycle recovery. Yet, where and how massive carbon emissions were mitigated during this climate warming event remains largely unknown. Here we show that organic carbon burial in the vast epicontinental seaways that extended over Eurasia provided a major carbon sink during the Palaeocene-Eocene Thermal Maximum. We coupled new and existing stratigraphic analyses to a detailed paleogeographic framework and using spatiotemporal interpolation calculated ca. 720-1300 Gt organic carbon excess burial, focused in the eastern parts of the Eurasian epicontinental seaways. A much larger amount (2160-3900 Gt C, and when accounting for the increase in inundated shelf area 7400-10300 Gt C) could have been sequestered in similar environments globally. With the disappearance of most epicontinental seas since the Oligocene-Miocene, an effective negative carbon cycle feedback also disappeared making the modern carbon cycle critically dependent on the slower silicate weathering feedback.}, language = {en} } @article{KloseGuillemoteauVignolietal.2023, author = {Klose, Tim and Guillemoteau, Julien and Vignoli, Giulio and Walter, Judith and Herrmann, Andreas and Tronicke, Jens}, title = {Structurally constrained inversion by means of a Minimum Gradient Support regularizer: examples of FD-EMI data inversion constrained by GPR reflection data}, series = {Geophysical journal international}, volume = {233}, journal = {Geophysical journal international}, number = {3}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggad041}, pages = {1938 -- 1949}, year = {2023}, abstract = {Many geophysical inverse problems are known to be ill-posed and, thus, requiring some kind of regularization in order to provide a unique and stable solution. A possible approach to overcome the inversion ill-posedness consists in constraining the position of the model interfaces. For a grid-based parameterization, such a structurally constrained inversion can be implemented by adopting the usual smooth regularization scheme in which the local weight of the regularization is reduced where an interface is expected. By doing so, sharp contrasts are promoted at interface locations while standard smoothness constraints keep affecting the other regions of the model. In this work, we present a structurally constrained approach and test it on the inversion of frequency-domain electromagnetic induction (FD-EMI) data using a regularization approach based on the Minimum Gradient Support stabilizer, which is capable to promote sharp transitions everywhere in the model, i.e., also in areas where no structural a prioriinformation is available. Using 1D and 2D synthetic data examples, we compare the proposed approach to a structurally constrained smooth inversion as well as to more standard (i.e., not structurally constrained) smooth and sharp inversions. Our results demonstrate that the proposed approach helps in finding a better and more reliable reconstruction of the subsurface electrical conductivity distribution, including its structural characteristics. Furthermore, we demonstrate that it allows to promote sharp parameter variations in areas where no structural information are available. Lastly, we apply our structurally constrained scheme to FD-EMI field data collected at a field site in Eastern Germany to image the thickness of peat deposits along two selected profiles. In this field example, we use collocated constant offset ground-penetrating radar (GPR) data to derive structural a priori information to constrain the inversion of the FD-EMI data. The results of this case study demonstrate the effectiveness and flexibility of the proposed approach.}, language = {en} } @article{CescaSuganRudzinskietal.2022, author = {Cesca, Simone and Sugan, Monica and Rudzinski, Lukasz and Vajedian, Sanaz and Niemz, Peter and Plank, Simon and Petersen, Gesa and Deng, Zhiguo and Rivalta, Eleonora and Vuan, Alessandro and Linares, Milton Percy Plasencia and Heimann, Sebastian and Dahm, Torsten}, title = {Massive earthquake swarm driven by magmatic intrusion at the Bransfield Strait, Antarctica}, series = {Communications earth and environment}, volume = {3}, journal = {Communications earth and environment}, number = {1}, publisher = {Springer Nature}, address = {London}, issn = {2662-4435}, doi = {10.1038/s43247-022-00418-5}, pages = {11}, year = {2022}, abstract = {An earthquake swarm affected the Bransfield Strait, Antarctica, a unique rift basin in transition from intra-arc rifting to ocean spreading. The swarm, counting similar to 85,000 volcano-tectonic earthquakes since August 2020, is located close to the Orca submarine volcano, previously considered inactive. Simultaneously, geodetic data reported up to similar to 11 cm north-westward displacement over King George Island. We use a broad variety of geophysical data and methods to reveal the complex migration of seismicity, accompanying the intrusion of 0.26-0.56 km(3) of magma. Strike-slip earthquakes mark the intrusion at depth, while shallower normal faulting the similar to 20 km long lateral growth of a dike. Seismicity abruptly decreased after a Mw 6.0 earthquake, suggesting the magmatic dike lost pressure with the slipping of a large fault. A seafloor eruption is likely, but not confirmed by sea surface temperature anomalies. The unrest documents episodic magmatic intrusion in the Bransfield Strait, providing unique insights into active continental rifting.}, language = {en} } @article{IrrgangBendixenFarquharsonetal.2022, author = {Irrgang, Anna M. and Bendixen, Mette and Farquharson, Louise M. and Baranskaya, Alisa and Erikson, Li H. and Gibbs, Ann E. and Ogorodov, Stanislav A. and Overduin, Pier Paul and Lantuit, Hugues and Grigoriev, Mikhail N. and Jones, Benjamin M.}, title = {Drivers, dynamics and impacts of changing Arctic coasts}, series = {Nature reviews earth and environment}, volume = {3}, journal = {Nature reviews earth and environment}, number = {1}, publisher = {Nature Research}, address = {London}, issn = {2662-138X}, doi = {10.1038/s43017-021-00232-1}, pages = {39 -- 54}, year = {2022}, abstract = {Arctic coasts are vulnerable to the effects of climate change, including rising sea levels and the loss of permafrost, sea ice and glaciers. Assessing the influence of anthropogenic warming on Arctic coastal dynamics, however, is challenged by the limited availability of observational, oceanographic and environmental data. Yet, with the majority of permafrost coasts being erosive, coupled with projected intensification of erosion and flooding, understanding these changes is critical. In this Review, we describe the morphological diversity of Arctic coasts, discuss important drivers of coastal change, explain the specific sensitivity of Arctic coasts to climate change and provide an overview of pan-Arctic shoreline change and its multifaceted impacts. Arctic coastal changes impact the human environment by threatening coastal settlements, infrastructure, cultural sites and archaeological remains. Changing sediment fluxes also impact the natural environment through carbon, nutrient and pollutant release on a magnitude that remains difficult to predict. Increasing transdisciplinary and interdisciplinary collaboration efforts will build the foundation for identifying sustainable solutions and adaptation strategies to reduce future risks for those living on, working at and visiting the rapidly changing Arctic coast.}, language = {en} } @article{BereswillGatzMillerSuetal.2023, author = {Bereswill, Sarah and Gatz-Miller, Hannah and Su, Danyang and T{\"o}tzke, Christian and Kardjilov, Nikolay and Oswald, Sascha and Mayer, Klaus Ulrich}, title = {Coupling non-invasive imaging and reactive transport modeling to investigate water and oxygen dynamics in the root zone}, series = {Vadose zone journal}, volume = {22}, journal = {Vadose zone journal}, number = {5}, publisher = {Wiley}, address = {Hoboken}, issn = {1539-1663}, doi = {10.1002/vzj2.20268}, pages = {19}, year = {2023}, abstract = {Oxygen (O-2) availability in soils is vital for plant growth and productivity. The transport and consumption of O-2 in the root zone is closely linked to soil moisture content, the spatial distribution of roots, as well as structure and heterogeneity of the surrounding soil. In this study, we measure three-dimensional root system architecture and the spatiotemporal dynamics of soil moisture (\& theta;) and O-2 concentrations in the root zone of maize (Zea mays) via non-invasive imaging, and then construct and parameterize a reactive transport model based on the experimental data. The combination of three non-invasive imaging methods allowed for a direct comparison of simulation results with observations at high spatial and temporal resolution. In three different modeling scenarios, we investigated how the results obtained for different levels of conceptual complexity in the model were able to match measured \& theta; and O-2 concentration patterns. We found that the modeling scenario that considers heterogeneous soil structure and spatial variability of hydraulic parameters (permeability, porosity, and van Genuchten \& alpha; and n), better reproduced the measured \& theta; and O-2 patterns relative to a simple model with a homogenous soil domain. The results from our combined imaging and modeling analysis reveal that experimental O-2 and water dynamics can be reproduced quantitatively in a reactive transport model, and that O-2 and water dynamics are best characterized when conditions unique to the specific system beyond the distribution of roots, such as soil structure and its effect on water saturation and macroscopic gas transport pathways, are considered.}, language = {en} } @article{StoltnowWeisKorges2023, author = {Stoltnow, Malte and Weis, Philipp and Korges, Maximilian}, title = {Hydrological controls on base metal precipitation and zoning at the porphyry-epithermal transition constrained by numerical modeling}, series = {Scientific reports}, volume = {13}, journal = {Scientific reports}, number = {1}, publisher = {Springer Nature}, address = {London}, issn = {2045-2322}, doi = {10.1038/s41598-023-30572-5}, pages = {15}, year = {2023}, abstract = {Ore precipitation in porphyry copper systems is generally characterized by metal zoning (Cu-Mo to Zn-Pb-Ag), which is suggested to be variably related to solubility decreases during fluid cooling, fluid-rock interactions, partitioning during fluid phase separation and mixing with external fluids. Here, we present new advances of a numerical process model by considering published constraints on the temperature- and salinity-dependent solubility of Cu, Pb and Zn in the ore fluid. We quantitatively investigate the roles of vapor-brine separation, halite saturation, initial metal contents, fluid mixing and remobilization as first-order controls of the physical hydrology on ore formation. The results show that the magmatic vapor and brine phases ascend with different residence times but as miscible fluid mixtures, with salinity increases generating metal-undersaturated bulk fluids. The release rates of magmatic fluids affect the location of the thermohaline fronts, leading to contrasting mechanisms for ore precipitation: higher rates result in halite saturation without significant metal zoning, lower rates produce zoned ore shells due to mixing with meteoric water. Varying metal contents can affect the order of the final metal precipitation sequence. Redissolution of precipitated metals results in zoned ore shell patterns in more peripheral locations and also decouples halite saturation from ore precipitation.}, language = {en} } @article{SchmidtFranckeGrosseetal.2023, author = {Schmidt, Lena Katharina and Francke, Till and Grosse, Peter Martin and Mayer, Christoph and Bronstert, Axel}, title = {Reconstructing five decades of sediment export from two glacierized high-alpine catchments in Tyrol, Austria, using nonparametric regression}, series = {Hydrology and earth system sciences : HESS}, volume = {27}, journal = {Hydrology and earth system sciences : HESS}, number = {9}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1027-5606}, doi = {10.5194/hess-27-1841-2023}, pages = {1841 -- 1863}, year = {2023}, abstract = {Knowledge on the response of sediment export to recent climate change in glacierized areas in the European Alps is limited, primarily because long-term records of suspended sediment concentrations (SSCs) are scarce. Here we tested the estimation of sediment export of the past five decades using quantile regression forest (QRF), a nonparametric, multivariate regression based on random forest. The regression builds on short-term records of SSCs and long records of the most important hydroclimatic drivers (discharge, precipitation and air temperature - QPT). We trained independent models for two nested and partially glacier-covered catchments, Vent (98 km(2)) and Vernagt (11.4 km(2)), in the upper otztal in Tyrol, Austria (1891 to 3772 m a.s.l.), where available QPT records start in 1967 and 1975. To assess temporal extrapolation ability, we used two 2-year SSC datasets at gauge Vernagt, which are almost 20 years apart, for a validation. For Vent, we performed a five-fold cross-validation on the 15 years of SSC measurements. Further, we quantified the number of days where predictors exceeded the range represented in the training dataset, as the inability to extrapolate beyond this range is a known limitation of QRF. Finally, we compared QRF performance to sediment rating curves (SRCs). We analyzed the modeled sediment export time series, the predictors and glacier mass balance data for trends (Mann-Kendall test and Sen's slope estimator) and step-like changes (using the widely applied Pettitt test and a complementary Bayesian approach).Our validation at gauge Vernagt demonstrated that QRF performs well in estimating past daily sediment export (Nash-Sutcliffe efficiency (NSE) of 0.73) and satisfactorily for SSCs (NSE of 0.51), despite the small training dataset. The temporal extrapolation ability of QRF was superior to SRCs, especially in periods with high-SSC events, which demonstrated the ability of QRF to model threshold effects. Days with high SSCs tended to be underestimated, but the effect on annual yields was small. Days with predictor exceedances were rare, indicating a good representativity of the training dataset. Finally, the QRF reconstruction models outperformed SRCs by about 20 percent points of the explained variance.Significant positive trends in the reconstructed annual suspended sediment yields were found at both gauges, with distinct step-like increases around 1981. This was linked to increased glacier melt, which became apparent through step-like increases in discharge at both gauges as well as change points in mass balances of the two largest glaciers in the Vent catchment. We identified exceptionally high July temperatures in 1982 and 1983 as a likely cause. In contrast, we did not find coinciding change points in precipitation. Opposing trends at the two gauges after 1981 suggest different timings of "peak sediment". We conclude that, given large-enough training datasets, the presented QRF approach is a promising tool with the ability to deepen our understanding of the response of high-alpine areas to decadal climate change.}, language = {en} } @article{CabiecesOlivar‐CastanoJunqueiraetal.2022, author = {Cabieces, Roberto and Olivar-Casta{\~n}o, Andr{\´e}s and Junqueira, Thiago C. and Relinque, Jes{\´u}s and Fernandez-Prieto, Luis M. and Vack{\´a}r, Jiř{\´i} and R{\"o}sler, Boris and Barco, Jaime and Pazos, Antonio and Garc{\´i}a-Mart{\´i}nez, Luz}, title = {Integrated Seismic Program (ISP): A new Python GUI-based software for earthquake seismology and seismic signal processing}, series = {Seismological research letters}, volume = {93}, journal = {Seismological research letters}, number = {3}, publisher = {Seismological Society of America}, address = {Albany}, issn = {0895-0695}, doi = {10.1785/0220210205}, pages = {1895 -- 1908}, year = {2022}, abstract = {Integrated Seismic Program (ISP) is a graphical user interface designed to facilitate and provide a user-friendly framework for performing diverse common and advanced tasks in seismological research. ISP is composed of five main modules for earthquake location, time-frequency analysis and advanced signal processing, implementation of array techniques to estimate the slowness vector, seismic moment tensor inversion, and receiver function computation and analysis. In addition, several support tools are available, allowing the user to create an event database, download data from International Federation of Digital Seismograph Networks services, inspect the background noise, and compute synthetic seismograms. ISP is written in Python3, supported by several open-source and/or publicly available tools. Its modular design allows for new features to be added in a collaborative development environment.}, language = {en} } @article{BlanchardPetitgirardLaurenzetal.2022, author = {Blanchard, Ingrid and Petitgirard, Sylvain and Laurenz, Vera and Miyajima, Nobuyoshi and Wilke, Max and Rubie, David C. and Lobanov, Sergey S. and Hennet, Louis and Morgenroth, Wolfgang and Tucoulou, R{\´e}mi and Bonino, Valentina and Zhao, Xuchao and Franchi, Ian}, title = {Chemical analysis of trace elements at the nanoscale in samples recovered from laser-heated diamond anvil cell experiments}, series = {Physics and chemistry of minerals}, volume = {49}, journal = {Physics and chemistry of minerals}, number = {6}, publisher = {Springer}, address = {New York}, issn = {0342-1791}, doi = {10.1007/s00269-022-01193-7}, pages = {16}, year = {2022}, abstract = {High pressure and high temperature experiments performed with laser-heated diamond anvil cells (LH-DAC) are being extensively used in geosciences to study matter at conditions prevailing in planetary interiors. Due to the size of the apparatus itself, the samples that are produced are extremely small, on the order of few tens of micrometers. There are several ways to analyze the samples and extract physical, chemical or structural information, using either in situ or ex situ methods. In this paper, we compare two nanoprobe techniques, namely nano-XRF and NanoSIMS, that can be used to analyze recovered samples synthetized in a LH-DAC. With these techniques, it is possible to extract the spatial distribution of chemical elements in the samples. We show the results for several standards and discuss the importance of proper calibration for the acquisition of quantifiable results. We used these two nanoprobe techniques to retrieve elemental ratios of dilute species (few tens of ppm) in quenched experimental molten samples relevant for the formation of the iron-rich core of the Earth. We finally discuss the applications of such probes to constrain the partitioning of trace elements between metal and silicate phases, with a focus on moderately siderophile elements, tungsten and molybdenum.}, 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{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{TimmermanKrmicekKrmičkovaetal.2023, author = {Timmerman, Martin Jan and Krmicek, Lukas and Krm{\´i}čkov{\´a}, Simona and Slama, Jiri and Sudo, Masafumi and Sobel, Edward}, title = {Tonian-Ediacaran evolution of the Brunovistulian microcontinent (Czech Republic) deciphered from LA-ICP-MS U-Pb zircon and 40Ar/39Ar muscovite ages}, series = {Precambrian research}, volume = {387}, journal = {Precambrian research}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0301-9268}, doi = {10.1016/j.precamres.2023.106981}, pages = {20}, year = {2023}, abstract = {Granitoids of the Slavkov Domain of the Brunovistulian microcontinent (BVM) in the Czech Republic have Ediacaran U-Pb zircon crystallization ages with the dominant magmatic activity occurring between ca. 597 and 595 Ma. The ages overlap published ages for the adjacent Thaya Domain, showing that both domains formed coevally in the same subduction setting. The data support published models in which the Slavkov Domain formed as arc crust. The main stage of magmatism stopped after ca. 595-590 Ma and was quickly followed by cooling accompanied by intrusion of small volumes of rhyolite dykes at ca. 594 Ma. Slavkov Domain metasedimentary rocks are dominated by Cryogenian-Ediacaran detrital zircon populations and their protoliths were locally derived erosional products of Cryogenian to Ediacaran arc rocks of the Thaya and Slavkov domains. Metasedi-mentary rocks from the NE part of the BVM contain younger, ca. 550 Ma zircons indicating that the BVM grew northeastward by accretion of progressively younger material derived from magmatic rocks with latest Ediacaran crystallization ages. In contrast to the Thaya and Slavkov domains, the Metavolcanic Zone that lies between them formed between ca. 740 and 725 Ma in the late Tonian to early Cryogenian. It predates the main stage magmatic activity in the BVM by 135 to 150 Ma and is probably a relic of older crust that formed during rifting of the Rodinia supercontinent. At ca. 552-551 Ma in the latest Ediacaran, parts of the BVM were exposed at the surface, during which time red, terrestrial siliciclastic sediments (Basal Clastics) were deposited. These largely had (very) proximal sources such as the main stage granitoids of the Thaya and Slavkov domains. Clasts of (meta)sandstones contain much older zircon populations and provide evidence that Neoarchaean and Palaeo-, meso- and early Neoproterozoic crustal rocks were exposed in erosional position nearby.}, language = {en} } @phdthesis{Pfalz2024, author = {Pfalz, Gregor}, title = {Advancing knowledge of Arctic lake system dynamics: A data-driven perspective on spatiotemporal patterns}, doi = {10.25932/publishup-63655}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-636554}, school = {Universit{\"a}t Potsdam}, pages = {xiii, 136}, year = {2024}, abstract = {Ecosystems play a pivotal role in addressing climate change but are also highly susceptible to drastic environmental changes. Investigating their historical dynamics can enhance our understanding of how they might respond to unprecedented future environmental shifts. With Arctic lakes currently under substantial pressure from climate change, lessons from the past can guide our understanding of potential disruptions to these lakes. However, individual lake systems are multifaceted and complex. Traditional isolated lake studies often fail to provide a global perspective because localized nuances—like individual lake parameters, catchment areas, and lake histories—can overshadow broader conclusions. In light of these complexities, a more nuanced approach is essential to analyze lake systems in a global context. A key to addressing this challenge lies in the data-driven analysis of sedimentological records from various northern lake systems. This dissertation emphasizes lake systems in the northern Eurasian region, particularly in Russia (n=59). For this doctoral thesis, we collected sedimentological data from various sources, which required a standardized framework for further analysis. Therefore, we designed a conceptual model for integrating and standardizing heterogeneous multi-proxy data into a relational database management system (PostgreSQL). Creating a database from the collected data enabled comparative numerical analyses between spatially separated lakes as well as between different proxies. When analyzing numerous lakes, establishing a common frame of reference was crucial. We achieved this by converting proxy values from depth dependency to age dependency. This required consistent age calculations across all lakes and proxies using one age-depth modeling software. Recognizing the broader implications and potential pitfalls of this, we developed the LANDO approach ("Linked Age and Depth Modelling"). LANDO is an innovative integration of multiple age-depth modeling software into a singular, cohesive platform (Jupyter Notebook). Beyond its ability to aggregate data from five renowned age-depth modeling software, LANDO uniquely empowers users to filter out implausible model outcomes using robust geoscientific data. Our method is not only novel but also significantly enhances the accuracy and reliability of lake analyses. Considering the preceding steps, this doctoral thesis further examines the relationship between carbon in sediments and temperature over the last 21,000 years. Initially, we hypothesized a positive correlation between carbon accumulation in lakes and modelled paleotemperature. Our homogenized dataset from heterogeneous lakes confirmed this association, even if the highest temperatures throughout our observation period do not correlate with the highest carbon values. We assume that rapid warming events contribute more to high accumulation, while sustained warming leads to carbon outgassing. Considering the current high concentration of carbon in the atmosphere and rising temperatures, ongoing climate change could cause northern lake systems to contribute to a further increase in atmospheric carbon (positive feedback loop). While our findings underscore the reliability of both our standardized data and the LANDO method, expanding our dataset might offer even greater assurance in our conclusions.}, language = {en} } @article{GluecklerGengGrimmetal.2022, author = {Gl{\"u}ckler, Ramesh and Geng, Rongwei and Grimm, Lennart and Baisheva, Izabella and Herzschuh, Ulrike and Stoof-Leichsenring, Kathleen R. and Kruse, Stefan and Andreev, Andrej Aleksandrovic and Pestryakova, Luidmila and Dietze, Elisabeth}, title = {Holocene wildfire and vegetation dynamics in Central Yakutia, Siberia, reconstructed from lake-sediment proxies}, 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.962906}, pages = {19}, year = {2022}, abstract = {Wildfires play an essential role in the ecology of boreal forests. In eastern Siberia, fire activity has been increasing in recent years, challenging the livelihoods of local communities. Intensifying fire regimes also increase disturbance pressure on the boreal forests, which currently protect the permafrost beneath from accelerated degradation. However, long-term relationships between changes in fire regime and forest structure remain largely unknown. We assess past fire-vegetation feedbacks using sedimentary proxy records from Lake Satagay, Central Yakutia, Siberia, covering the past c. 10,800 years. Results from macroscopic and microscopic charcoal analyses indicate high amounts of burnt biomass during the Early Holocene, and that the present-day, low-severity surface fire regime has been in place since c. 4,500 years before present. A pollen-based quantitative reconstruction of vegetation cover and a terrestrial plant record based on sedimentary ancient DNA metabarcoding suggest a pronounced shift in forest structure toward the Late Holocene. Whereas the Early Holocene was characterized by postglacial open larch-birch woodlands, forest structure changed toward the modern, mixed larch-dominated closed-canopy forest during the Mid-Holocene. We propose a potential relationship between open woodlands and high amounts of burnt biomass, as well as a mediating effect of dense larch forest on the climate-driven intensification of fire regimes. Considering the anticipated increase in forest disturbances (droughts, insect invasions, and wildfires), higher tree mortality may force the modern state of the forest to shift toward an open woodland state comparable to the Early Holocene. Such a shift in forest structure may result in a positive feedback on currently intensifying wildfires. These new long-term data improve our understanding of millennial-scale fire regime changes and their relationships to changes of vegetation in Central Yakutia, where the local population is already being confronted with intensifying wildfire seasons.}, language = {en} } @article{InceogluShpritsHeinemannetal.2022, author = {Inceoglu, Fadil and Shprits, Yuri and Heinemann, Stephan G. and Bianco, Stefano}, title = {Identification of coronal holes on AIA/SDO images using unsupervised machine learning}, series = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, volume = {930}, journal = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, number = {2}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1538-4357}, doi = {10.3847/1538-4357/ac5f43}, pages = {11}, year = {2022}, abstract = {Through its magnetic activity, the Sun governs the conditions in Earth's vicinity, creating space weather events, which have drastic effects on our space- and ground-based technology. One of the most important solar magnetic features creating the space weather is the solar wind that originates from the coronal holes (CHs). The identification of the CHs on the Sun as one of the source regions of the solar wind is therefore crucial to achieve predictive capabilities. In this study, we used an unsupervised machine-learning method, k-means, to pixel-wise cluster the passband images of the Sun taken by the Atmospheric Imaging Assembly on the Solar Dynamics Observatory in 171, 193, and 211 angstrom in different combinations. Our results show that the pixel-wise k-means clustering together with systematic pre- and postprocessing steps provides compatible results with those from complex methods, such as convolutional neural networks. More importantly, our study shows that there is a need for a CH database where a consensus about the CH boundaries is reached by observers independently. This database then can be used as the "ground truth," when using a supervised method or just to evaluate the goodness of the models.}, language = {en} } @article{GengAndreevKruseetal.2022, author = {Geng, Rongwei and Andreev, Andrei and Kruse, Stefan and Heim, Birgit and van Geffen, Femke and Pestryakova, Luidmila and Zakharov, Evgenii and Troeva, Elena I. and Shevtsova, Iuliia and Li, Furong and Zhao, Yan and Herzschuh, Ulrike}, title = {Modern pollen assemblages from lake sediments and soil in East Siberia and relative pollen productivity estimates for Major Taxa}, 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.837857}, pages = {17}, year = {2022}, abstract = {Modern pollen-vegetation-climate relationships underpin palaeovegetation and palaeoclimate reconstructions from fossil pollen records. East Siberia is an ideal area for investigating the relationships between modern pollen assemblages and near natural vegetation under cold continental climate conditions. Reliable pollen-based quantitative vegetation and climate reconstructions are still scarce due to the limited number of modern pollen datasets. Furthermore, differences in pollen representation of samples from lake sediments and soils are not well understood. Here, we present a new pollen dataset of 48 moss/soil and 24 lake surface-sediment samples collected in Chukotka and central Yakutia in East Siberia. The pollen-vegetation-climate relationships were investigated by ordination analyses. Generally, tundra and taiga vegetation types can be well distinguished in the surface pollen assemblages. Moss/soil and lake samples contain generally similar pollen assemblages as revealed by a Procrustes comparison with some exceptions. Overall, modern pollen assemblages reflect the temperature and precipitation gradients in the study areas as revealed by constrained ordination analysis. We estimate the relative pollen productivity (RPP) of major taxa and the relevant source area of pollen (RSAP) for moss/soil samples from Chukotka and central Yakutia using Extended R-Value (ERV) analysis. The RSAP of the tundra-forest transition area in Chukotka and taiga area in central Yakutia are ca. 1300 and 360 m, respectively. For Chukotka, RPPs relative to both Poaceae and Ericaceae were estimated while RPPs for central Yakutia were relative only to Ericaceae. Relative to Ericaceae (reference taxon, RPP = 1), Larix, Betula, Picea, and Pinus are overrepresented while Alnus, Cyperaceae, Poaceae, and Salix are underrepresented in the pollen spectra. Our estimates are in general agreement with previously published values and provide the basis for reliable quantitative reconstructions of East Siberian vegetation.}, 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{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{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{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} } @article{AstudilloSotomayorJaraMunozMelnicketal.2021, author = {Astudillo-Sotomayor, Luis and Jara Mu{\~n}oz, Julius and Melnick, Daniel and Cort{\´e}s-Aranda, Joaqu{\´i}n and Tassara, Andr{\´e}s and Strecker, Manfred}, title = {Fast Holocene slip and localized strain along the Liqui{\~n}e-Ofqui strike-slip fault system, Chile}, series = {Scientific reports}, volume = {11}, journal = {Scientific reports}, number = {1}, publisher = {Macmillan Publishers Limited, part of Springer Nature}, address = {London}, issn = {2045-2322}, doi = {10.1038/s41598-021-85036-5}, pages = {10}, year = {2021}, abstract = {In active tectonic settings dominated by strike-slip kinematics, slip partitioning across subparallel faults is a common feature; therefore, assessing the degree of partitioning and strain localization is paramount for seismic hazard assessments. Here, we estimate a slip rate of 18.8 +/- 2.0 mm/year over the past 9.0 +/- 0.1 ka for a single strand of the Liquirie-Ofqui Fault System, which straddles the Main Cordillera in Southern Chile. This Holocene rate accounts for similar to 82\% of the trench-parallel component of oblique plate convergence and is similar to million-year estimates integrated over the entire fault system. Our results imply that strain localizes on a single fault at millennial time scale but over longer time scales strain localization is not sustained. The fast millennial slip rate in the absence of historical Mw> 6.5 earthquakes along the Liquine-Ofqui Fault System implies either a component of aseismic slip or Mw similar to 7 earthquakes involving multi-trace ruptures and > 150-year repeat times. Our results have implications for the understanding of strike-slip fault system dynamics within volcanic arcs and seismic hazard assessments.}, language = {en} } @article{LopezCominoCescaNiemzetal.2021, author = {L{\´o}pez-Comino, Jos{\´e} {\´A}ngel and Cesca, Simone and Niemz, Peter and Dahm, Torsten and Zang, Arno}, title = {Rupture directivity in 3D inferred from acoustic emissions events in a mine-scale hydraulic fracturing experiment}, 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.670757}, pages = {9}, year = {2021}, abstract = {Rupture directivity, implying a predominant earthquake rupture propagation direction, is typically inferred upon the identification of 2D azimuthal patterns of seismic observations for weak to large earthquakes using surface-monitoring networks. However, the recent increase of 3D monitoring networks deployed in the shallow subsurface and underground laboratories toward the monitoring of microseismicity allows to extend the directivity analysis to 3D modeling, beyond the usual range of magnitudes. The high-quality full waveforms recorded for the largest, decimeter-scale acoustic emission (AE) events during a meter-scale hydraulic fracturing experiment in granites at similar to 410 m depth allow us to resolve the apparent durations observed at each AE sensor to analyze 3D-directivity effects. Unilateral and (asymmetric) bilateral ruptures are then characterized by the introduction of a parameter kappa, representing the angle between the directivity vector and the station vector. While the cloud of AE activity indicates the planes of the hydrofractures, the resolved directivity vectors show off-plane orientations, indicating that rupture planes of microfractures on a scale of centimeters have different geometries. Our results reveal a general alignment of the rupture directivity with the orientation of the minimum horizontal stress, implying that not only the slip direction but also the fracture growth produced by the fluid injections is controlled by the local stress conditions.}, language = {en} } @article{EsfahaniGholamiOhrnberger2020, author = {Esfahani, Reza Dokht Dolatabadi and Gholami, Ali and Ohrnberger, Matthias}, title = {An inexact augmented Lagrangian method for nonlinear dispersion-curve inversion using Dix-type global linear approximation}, series = {Geophysics : a journal of general and applied geophysics}, volume = {85}, journal = {Geophysics : a journal of general and applied geophysics}, number = {3}, publisher = {GeoScienceWorld}, address = {Tulsa, Okla.}, issn = {0016-8033}, doi = {10.1190/geo2019-0717.1}, pages = {EN77 -- EN85}, year = {2020}, abstract = {Dispersion-curve inversion of Rayleigh waves to infer subsurface shear-wave velocity is a long-standing problem in seismology. Due to nonlinearity and ill-posedness, sophisticated regularization techniques are required to solve the problem for a stable velocity model. We have formulated the problem as a minimization problem with nonlinear operator constraint and then solve it by using an inexact augmented Lagrangian method, taking advantage of the Haney-Tsai Dix-type relation (a global linear approximation of the nonlinear forward operator). This replaces the original regularized nonlinear problem with iterative minimization of a more tractable regularized linear problem followed by a nonlinear update of the phase velocity (data) in which the update can be performed accurately with any forward modeling engine, for example, the finite-element method. The algorithm allows discretizing the medium with thin layers (for the finite-element method) and thus omitting the layer thicknesses from the unknowns and also allows incorporating arbitrary regularizations to shape the desired velocity model. In this research, we use total variation regularization to retrieve the shear-wave velocity model. We use two synthetic and two real data examples to illustrate the performance of the inversion algorithm with total variation regularization. We find that the method is fast and stable, and it converges to the solution of the original nonlinear problem.}, 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} } @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} } @phdthesis{Koyan2024, author = {Koyan, Philipp}, title = {3D attribute analysis and classification to interpret ground-penetrating radar (GPR) data collected across sedimentary environments: Synthetic studies and field examples}, doi = {10.25932/publishup-63948}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-639488}, school = {Universit{\"a}t Potsdam}, pages = {xi, 115, A51}, year = {2024}, abstract = {Die Untersuchung des oberfl{\"a}chennahen Untergrundes erfolgt heutzutage bei Frage- stellungen aus den Bereichen des Bauwesens, der Arch{\"a}ologie oder der Geologie und Hydrologie oft mittels zerst{\"o}rungsfreier beziehungsweise zerst{\"o}rungsarmer Methoden der angewandten Geophysik. Ein Bereich, der eine immer zentralere Rolle in Forschung und Ingenieurwesen einnimmt, ist die Untersuchung von sediment{\"a}ren Umgebungen, zum Beispiel zur Charakterisierung oberfl{\"a}chennaher Grundwassersysteme. Ein in diesem Kontext h{\"a}ufig eingesetztes Verfahren ist das des Georadars (oftmals GPR - aus dem Englischen ground-penetrating radar). Dabei werden kurze elektromagnetische Impulse von einer Antenne in den Untergrund ausgesendet, welche dort wiederum an Kontrasten der elektromagnetischen Eigenschaften (wie zum Beispiel an der Grundwasseroberfl{\"a}che) reflektiert, gebrochen oder gestreut werden. Eine Empfangsantenne zeichnet diese Signale in Form derer Amplituden und Laufzeiten auf. Eine Analyse dieser aufgezeichneten Signale erm{\"o}glicht Aussagen {\"u}ber den Untergrund, beispielsweise {\"u}ber die Tiefenlage der Grundwasseroberfl{\"a}che oder die Lagerung und Charakteristika oberfl{\"a}chennaher Sedimentschichten. Dank des hohen Aufl{\"o}sungsverm{\"o}gens der GPR-Methode sowie stetiger technologischer Entwicklungen erfolgt heutzutage die Aufzeichnung von GPR- Daten immer h{\"a}ufiger in 3D. Trotz des hohen zeitlichen und technischen Aufwandes f{\"u}r die Datenakquisition und -bearbeitung werden die resultierenden 3D-Datens{\"a}tze, welche den Untergrund hochaufl{\"o}send abbilden, typischerweise von Hand interpretiert. Dies ist in der Regel ein {\"a}ußerst zeitaufwendiger Analyseschritt. Daher werden oft repr{\"a}sentative 2D-Schnitte aus dem 3D-Datensatz gew{\"a}hlt, in denen markante Reflektionsstrukuren markiert werden. Aus diesen Strukturen werden dann sich {\"a}hnelnde Bereiche im Untergrund als so genannte Radar-Fazies zusammengefasst. Die anhand von 2D-Schnitten erlangten Resultate werden dann als repr{\"a}sentativ f{\"u}r die gesamte untersuchte Fl{\"a}che angesehen. In dieser Form durchgef{\"u}hrte Interpretationen sind folglich oft unvollst{\"a}ndig sowie zudem in hohem Maße von der Expertise der Interpretierenden abh{\"a}ngig und daher in der Regel nicht reproduzierbar. Eine vielversprechende Alternative beziehungsweise Erg{\"a}nzung zur manuellen In- terpretation ist die Verwendung von so genannten GPR-Attributen. Dabei werden nicht die aufgezeichneten Daten selbst, sondern daraus abgeleitete Gr{\"o}ßen, welche die markanten Reflexionsstrukturen in 3D charakterisieren, zur Interpretation herangezogen. In dieser Arbeit wird anhand verschiedener Feld- und Modelldatens{\"a}tze untersucht, welche Attribute sich daf{\"u}r insbesondere eignen. Zudem zeigt diese Arbeit, wie ausgew{\"a}hlte Attribute mittels spezieller Bearbeitungs- und Klassifizierungsmethoden zur Erstellung von 3D-Faziesmodellen genutzt werden k{\"o}nnen. Dank der M{\"o}glichkeit der Erstellung so genannter attributbasierter 3D-GPR-Faziesmodelle k{\"o}nnen zuk{\"u}nftige Interpretationen zu gewissen Teilen automatisiert und somit effizienter durchgef{\"u}hrt werden. Weiterhin beschreiben die so erhaltenen Resultate den untersuchten Untergrund in reproduzierbarer Art und Weise sowie umf{\"a}nglicher als es bisher mittels manueller Interpretationsmethoden typischerweise m{\"o}glich war.}, 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{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{NooshiriBeanDahmetal.2021, author = {Nooshiri, Nima and Bean, Christopher J. and Dahm, Torsten and Grigoli, Francesco and Kristjansdottir, Sigriour and Obermann, Anne and Wiemer, Stefan}, title = {A multibranch, multitarget neural network for rapid point-source inversion in a microseismic environment}, series = {Geophysical journal international}, volume = {229}, journal = {Geophysical journal international}, number = {2}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggab511}, pages = {999 -- 1016}, year = {2021}, abstract = {Despite advanced seismological techniques, automatic source characterization for microseismic earthquakes remains difficult and challenging since current inversion and modelling of high-frequency signals are complex and time consuming. For real-time applications such as induced seismicity monitoring, the application of standard methods is often not fast enough for true complete real-time information on seismic sources. In this paper, we present an alternative approach based on recent advances in deep learning for rapid source-parameter estimation of microseismic earthquakes. The seismic inversion is represented in compact form by two convolutional neural networks, with individual feature extraction, and a fully connected neural network, for feature aggregation, to simultaneously obtain full moment tensor and spatial location of microseismic sources. Specifically, a multibranch neural network algorithm is trained to encapsulate the information about the relationship between seismic waveforms and underlying point-source mechanisms and locations. The learning-based model allows rapid inversion (within a fraction of second) once input data are available. A key advantage of the algorithm is that it can be trained using synthetic seismic data only, so it is directly applicable to scenarios where there are insufficient real data for training. Moreover, we find that the method is robust with respect to perturbations such as observational noise and data incompleteness (missing stations). We apply the new approach on synthesized and example recorded small magnitude (M <= 1.6) earthquakes at the Hellisheioi geothermal field in the Hengill area, Iceland. For the examined events, the model achieves excellent performance and shows very good agreement with the inverted solutions determined through standard methodology. In this study, we seek to demonstrate that this approach is viable for microseismicity real-time estimation of source parameters and can be integrated into advanced decision-support tools for controlling induced seismicity.}, language = {en} } @article{EzpeletaParraColloetal.2022, author = {Ezpeleta, Miguel and Parra, Mauricio and Collo, Gilda and Wunderlin, Cecilia and Borrego, Angeles G. and Sobel, Edward and Glodny, Johannes}, title = {Thermochronometry unveils ancient thermal regimes in the NW Pampean Ranges, Argentina}, series = {Basin research}, volume = {34}, journal = {Basin research}, number = {6}, publisher = {Wiley}, address = {Hoboken}, issn = {0950-091X}, doi = {10.1111/bre.12693}, pages = {1983 -- 2012}, year = {2022}, abstract = {Reconstructing thermal histories in thrust belts is commonly used to infer the age and rates of thrusting and hence the driving mechanisms of orogenesis. In areas where ancient basins have been incorporated into the orogenic wedge, a quantitative reconstruction of the thermal history helps distinguish among potential mechanisms responsible for heating events. We present such a reconstruction for the Ischigualasto-Villa Union basin in the western Pampean Ranges of Argentina, where Triassic rifting and late Cretaceous-Cenozoic retroarc foreland basin development has been widely documented, including Miocene flat-slab subduction. We report results of organic and inorganic thermal indicators acquired along three stratigraphic sections, including vitrinite reflectance and X-ray diffractometry in claystones and new thermochronological [(apatite fission-track and apatite and zircon [U-Th]/He)] analyses. Despite up to 5 km-thick Cenozoic overburden and unlike previously thought, the thermal peak in the basin is not due to Cenozoic burial but occurred in the Triassic, associated with a high heat flow of up to 90 mWm(-2) and <2 km of burial, which heated the base of the Triassic strata to similar to 160 degrees C. Following exhumation, attested by the development of an unconformity between the Triassic and Late-Cretaceous-Cenozoic sequences, Cenozoic re-burial increased the temperature to similar to 110 degrees C at the base of the Triassic section and only similar to 50 degrees C 7 km upsection, suggesting a dramatic decrease in the thermal gradient. The onset of Cenozoic cooling occurred at similar to 10(-8) Ma, concomitant with sediment accumulation and thus preceding the latest Miocene onset of thrusting that has been independently documented by stratigraphic-cross-cutting relationships. We argue that the onset of cooling is associated with lithospheric refrigeration following establishment of flat-slab subduction, leading to the eastward displacement of the asthenospheric wedge beneath the South American plate. Our study places time and temperature constraints on flat-slab cooling that calls for a careful interpretation of exhumation signals in thrustbelts inferred from thermochronology only.}, language = {en} } @article{FernandezPalominoHattermannKrysanovaetal.2022, author = {Fernandez-Palomino, Carlos Antonio and Hattermann, Fred and Krysanova, Valentina and Lobanova, Anastasia and Vega-Jacome, Fiorella and Lavado, Waldo and Santini, William and Aybar, Cesar and Bronstert, Axel}, title = {A novel high-resolution gridded precipitation dataset for peruvian and ecuadorian watersheds}, series = {Journal of hydrometeorology}, volume = {23}, journal = {Journal of hydrometeorology}, number = {3}, publisher = {American Meteorological Soc.}, address = {Boston}, issn = {1525-755X}, doi = {10.1175/JHM-D-20-0285.1}, pages = {309 -- 336}, year = {2022}, abstract = {A novel approach for estimating precipitation patterns is developed here and applied to generate a new hydrologically corrected daily precipitation dataset, called RAIN4PE (Rain for Peru and Ecuador), at 0.1 degrees spatial resolution for the period 1981-2015 covering Peru and Ecuador. It is based on the application of 1) the random forest method to merge multisource precipitation estimates (gauge, satellite, and reanalysis) with terrain elevation, and 2) observed and modeled streamflow data to first detect biases and second further adjust gridded precipitation by inversely applying the simulated results of the ecohydrological model SWAT (Soil and Water Assessment Tool). Hydrological results using RAIN4PE as input for the Peruvian and Ecuadorian catchments were compared against the ones when feeding other uncorrected (CHIRP and ERA5) and gauge-corrected (CHIRPS, MSWEP, and PISCO) precipitation datasets into the model. For that, SWAT was calibrated and validated at 72 river sections for each dataset using a range of performance metrics, including hydrograph goodness of fit and flow duration curve signatures. Results showed that gauge-corrected precipitation datasets outperformed uncorrected ones for streamflow simulation. However, CHIRPS, MSWEP, and PISCO showed limitations for streamflow simulation in several catchments draining into the Pacific Ocean and the Amazon River. RAIN4PE provided the best overall performance for streamflow simulation, including flow variability (low, high, and peak flows) and water budget closure. The overall good performance of RAIN4PE as input for hydrological modeling provides a valuable criterion of its applicability for robust countrywide hydrometeorological applications, including hydroclimatic extremes such as droughts and floods. Significance StatementWe developed a novel precipitation dataset RAIN4PE for Peru and Ecuador by merging multisource precipitation data (satellite, reanalysis, and ground-based precipitation) with terrain elevation using the random forest method. Furthermore, RAIN4PE was hydrologically corrected using streamflow data in watersheds with precipitation underestimation through reverse hydrology. The results of a comprehensive hydrological evaluation showed that RAIN4PE outperformed state-of-the-art precipitation datasets such as CHIRP, ERA5, CHIRPS, MSWEP, and PISCO in terms of daily and monthly streamflow simulations, including extremely low and high flows in almost all Peruvian and Ecuadorian catchments. This underlines the suitability of RAIN4PE for hydrometeorological applications in this region. Furthermore, our approach for the generation of RAIN4PE can be used in other data-scarce regions.}, language = {en} } @article{SchuelerCalabreseAttinger2021, author = {Sch{\"u}ler, Lennart and Calabrese, Justin M. and Attinger, Sabine}, title = {Data driven high resolution modeling and spatial analyses of the COVID-19 pandemic in Germany}, series = {PLoS one}, volume = {16}, journal = {PLoS one}, number = {8}, publisher = {PLoS}, address = {San Fransisco}, issn = {1932-6203}, doi = {10.1371/journal.pone.0254660}, pages = {14}, year = {2021}, abstract = {The SARS-CoV-2 virus has spread around the world with over 100 million infections to date, and currently many countries are fighting the second wave of infections. With neither sufficient vaccination capacity nor effective medication, non-pharmaceutical interventions (NPIs) remain the measure of choice. However, NPIs place a great burden on society, the mental health of individuals, and economics. Therefore the cost/benefit ratio must be carefully balanced and a target-oriented small-scale implementation of these NPIs could help achieve this balance. To this end, we introduce a modified SEIRD-class compartment model and parametrize it locally for all 412 districts of Germany. The NPIs are modeled at district level by time varying contact rates. This high spatial resolution makes it possible to apply geostatistical methods to analyse the spatial patterns of the pandemic in Germany and to compare the results of different spatial resolutions. We find that the modified SEIRD model can successfully be fitted to the COVID-19 cases in German districts, states, and also nationwide. We propose the correlation length as a further measure, besides the weekly incidence rates, to describe the current situation of the epidemic.}, language = {en} } @article{HuangStoofLeichsenringLiuetal.2021, author = {Huang, Sichao and Stoof-Leichsenring, Kathleen R. and Liu, Sisi and Courtin, Jeremy and Andreev, Andrej A. and Pestryakova, Luidmila. A. and Herzschuh, Ulrike}, title = {Plant sedimentary ancient DNA from Far East Russia covering the last 28,000 years reveals different assembly rules in cold and warm climates}, series = {Frontiers in Ecology and Evolution}, volume = {9}, journal = {Frontiers in Ecology and Evolution}, publisher = {Frontiers Media}, address = {Lausanne}, issn = {2296-701X}, doi = {10.3389/fevo.2021.763747}, pages = {15}, year = {2021}, abstract = {Woody plants are expanding into the Arctic in response to the warming climate. The impact on arctic plant communities is not well understood due to the limited knowledge about plant assembly rules. Records of past plant diversity over long time series are rare. Here, we applied sedimentary ancient DNA metabarcoding targeting the P6 loop of the chloroplast trnL gene to a sediment record from Lake Ilirney (central Chukotka, Far Eastern Russia) covering the last 28 thousand years. Our results show that forb-rich steppe-tundra and dwarf-shrub tundra dominated during the cold climate before 14 ka, while deciduous erect-shrub tundra was abundant during the warm period since 14 ka. Larix invasion during the late Holocene substantially lagged behind the likely warmest period between 10 and 6 ka, where the vegetation biomass could be highest. We reveal highest richness during 28-23 ka and a second richness peak during 13-9 ka, with both periods being accompanied by low relative abundance of shrubs. During the cold period before 14 ka, rich plant assemblages were phylogenetically clustered, suggesting low genetic divergence in the assemblages despite the great number of species. This probably originates from environmental filtering along with niche differentiation due to limited resources under harsh environmental conditions. In contrast, during the warmer period after 14 ka, rich plant assemblages were phylogenetically overdispersed. This results from a high number of species which were found to harbor high genetic divergence, likely originating from an erratic recruitment process in the course of warming. Some of our evidence may be of relevance for inferring future arctic plant assembly rules and diversity changes. By analogy to the past, we expect a lagged response of tree invasion. Plant richness might overshoot in the short term; in the long-term, however, the ongoing expansion of deciduous shrubs will eventually result in a phylogenetically more diverse community.}, 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{GhaffarJomaaMeonetal.2021, author = {Ghaffar, Salman and Jomaa, Seifeddine and Meon, G{\"u}nter and Rode, Michael}, title = {Spatial validation of a semi-distributed hydrological nutrient transport model}, series = {Journal of hydrology}, volume = {593}, journal = {Journal of hydrology}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0022-1694}, doi = {10.1016/j.jhydrol.2020.125818}, pages = {13}, year = {2021}, abstract = {Semi-distributed hydrological and water quality models are increasingly used as innovative and scientific-based management tools. However, their application is usually restricted to the gauging stations where they are originally calibrated, limiting their spatial capability. In this study, the semi-distributed hydrological water quality model HYPE (HYdrological Predictions for the Environment) was tested spatially to represent nitrate-N (NO3- N) and total phosphorus (TP) concentrations and loads of the nested and heterogeneous Selke catchment (463 km(2)) in central Germany. First, an automatic calibration procedure and uncertainty analysis were conducted using the DiffeRential Evolution Adaptive Metropolis (DREAM) tool to simulate discharge, NO3--N and TP concentrations. A multi-site and multi-objective calibration approach was applied using three main gauging stations, covering the most important hydro-meteorological and physiographical characteristics of the whole catchment. Second, the model's capability was tested to represent further internal stations, which were not initially considered for calibration. Results showed that discharge was well represented by the model at all three main stations during both calibration (1994-1998) and validation (1999-2014) periods with lowest Nash-Sutcliffe Efficiency (NSE) of 0.71 and maximum Percentage BIAS (PBIAS) of 18.0\%. The model was able to reproduce the seasonal dynamics of NO3--N and TP concentrations with low predictive uncertainty at the three main stations, reflected by PBIAS values in the ranges from 16.1\% to 6.4\% and from 20.0\% to 11.5\% for NO3--N and TP load simulations, respectively. At internal stations, the model could represent reasonably well the seasonal variation of nutrient concentrations with PBIAS values in the ranges from 9.0\% to 14.2\% for NO3--N and from 25.3\% to 34.3\% for TP concentration simulations. Overall, results suggested that the spatial validation of a nutrient transport model can be better ensured when a multi-site and multi-objective calibration approach using archetypical gauging stations is implemented. Further, results revealed that the delineation of sub-catchments should put more focus on hydro-meteorological conditions than on land-use features.}, language = {en} } @article{JanochaSmyrakSikoraSengeretal.2021, author = {Janocha, Julian and Smyrak-Sikora, Aleksandra and Senger, Kim and Birchall, Thomas}, title = {Seeing beyond the outcrop}, series = {Marine and petroleum geology}, volume = {125}, journal = {Marine and petroleum geology}, publisher = {Elsevier Science}, address = {Amsterdam [u.a.]}, issn = {0264-8172}, doi = {10.1016/j.marpetgeo.2020.104833}, pages = {16}, year = {2021}, abstract = {Paleokarst breccias are a common feature of sedimentary rift basins. The Billefjorden Trough in the High Arctic archipelago of Svalbard is an example of such a rift. Here the Carboniferous stratigraphy exhibits intervals of paleokarst breccias formed by gypsum dissolution. In this study we integrate digital outcrop models (DOMs) with a 2D ground penetrating radar (GPR) survey to extrapolate external irregular paleokarst geometries beyond the 2D outcrops. DOMs are obtained through combining a series of overlapping photographs with structure-frommotion photogrammetry, to create mmto dm-resolution georeferenced DOMs. GPR is typically used for surveying the shallow subsurface and relies on detecting the contrasts in electro-magnetic permittivity. We defined three geophysical facies based on their appearance in GPR. By integrating subsurface geophysical data with DOMs we were able to correlate reflection patterns in GPR with outcrop features. The chaotic nature of paleokarst breccias is seen both in outcrop and GPR. Key horizons in outcrop and the GPR profiles allow tying together observations between these methods. Furthermore, we show that this technique expands the twodimensional outcrop surface into a three-dimensional domain, thus complementing, strengthening and extending outcrop interpretations.}, language = {en} } @article{GomezZapataZafrirPittoreetal.2022, author = {Gomez Zapata, Juan Camilo and Zafrir, Raquel and Pittore, Massimiliano and Merino, Yvonne}, title = {Towards a sensitivity analysis in seismic risk with probabilistic building exposure models}, series = {ISPRS International Journal of Geo-Information}, volume = {11}, journal = {ISPRS International Journal of Geo-Information}, number = {2}, publisher = {MDPI}, address = {Basel}, issn = {2220-9964}, doi = {10.3390/ijgi11020113}, pages = {38}, year = {2022}, abstract = {Efforts have been made in the past to enhance building exposure models on a regional scale with increasing spatial resolutions by integrating different data sources. This work follows a similar path and focuses on the downscaling of the existing SARA exposure model that was proposed for the residential building stock of the communes of Valparaiso and Vina del Mar (Chile). Although this model allowed great progress in harmonising building classes and characterising their differential physical vulnerabilities, it is now outdated, and in any case, it is spatially aggregated over large administrative units. Hence, to more accurately consider the impact of future earthquakes on these cities, it is necessary to employ more reliable exposure models. For such a purpose, we propose updating this existing model through a Bayesian approach by integrating ancillary data that has been made increasingly available from Volunteering Geo-Information (VGI) activities. Its spatial representation is also optimised in higher resolution aggregation units that avoid the inconvenience of having incomplete building-by-building footprints. A worst-case earthquake scenario is presented to calculate direct economic losses and highlight the degree of uncertainty imposed by exposure models in comparison with other parameters used to generate the seismic ground motions within a sensitivity analysis. This example study shows the great potential of using increasingly available VGI to update worldwide building exposure models as well as its importance in scenario-based seismic risk assessment.}, language = {en} } @article{HannemannEulenfeldKruegeretal.2021, author = {Hannemann, Katrin and Eulenfeld, Tom and Kr{\"u}ger, Frank and Dahm, Torsten}, title = {Seismic scattering and absorption of oceanic lithospheric S waves in the Eastern North Atlantic}, series = {Geophysical journal international}, volume = {229}, journal = {Geophysical journal international}, number = {2}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggab493}, pages = {948 -- 961}, year = {2021}, abstract = {The scattering and absorption of high-frequency seismic waves in the oceanic lithosphere is to date only poorly constrained by observations. Such estimates would not only improve our understanding of the propagation of seismic waves, but also unravel the small-scale nature of the lithosphere and its variability. Our study benefits from two exceptional situations: (1) we deployed over 10 months a mid-aperture seismological array in the central part of the Eastern North Atlantic in 5 km water depth and (2) we could observe in total 340 high-frequency (up to 30 Hz) Po and So arrivals with tens to hundreds of seconds long seismic coda from local and regional earthquakes in a wide range of backazimuths and epicentral distances up to 850 km with a travel path in the oceanic lithosphere. Moreover, the array was located about 100 km north of the Gloria fault, defining the plate boundary between the Eurasian and African plates at this location which also allows an investigation of the influence of an abrupt change in lithospheric age (20 Ma in this case) on seismic waves. The waves travel with velocities indicating upper-mantle material. We use So waves and their coda of pre-selected earthquakes to estimate frequency-dependent seismic scattering and intrinsic attenuation parameters. The estimated scattering attenuation coefficients are between 10(-4) and 4 x 10(-5) m(-1) and are typical for the lithosphere or the upper mantle. Furthermore, the total quality factors for So waves below 5 Hz are between 20 and 500 and are well below estimates from previous modelling for observations in the Pacific Ocean. This implies that the Atlantic Ocean is more attenuative for So waves compared to the Pacific Ocean, which is inline with the expected behaviour for the lithospheric structures resulting from the slower spreading rates in the Atlantic Ocean. The results for the analysed events indicate that for frequencies above 3 Hz, intrinsic attenuation is equal to or slightly stronger than scattering attenuation and that the So-wave coda is weakly influenced by the oceanic crust. Both observations are in agreement with the proposed propagation mechanism of scattering in the oceanic mantle lithosphere. Furthermore, we observe an age dependence which shows that an increase in lithospheric age is associated with a decrease in attenuation. However, we also observe a trade-off of this age-dependent effect with either a change in lithospheric thickness or thermal variations, for example due to small-scale upwellings in the upper mantle in the southeast close to Madeira and the Canaries. Moreover, the influence of the nearby Gloria fault is visible in a reduction of the intrinsic attenuation below 3 Hz for estimates across the fault. This is the first study to estimate seismic scattering and absorption parameters of So waves for an area with several hundreds of kilometres radius centred in the Eastern North Atlantic and using them to characterize the nature of the oceanic lithosphere.}, language = {en} } @article{BenNsirJomaaYildirimetal.2022, author = {Ben Nsir, Siwar and Jomaa, Seifeddine and Yildirim, Umit and Zhou, Xiangqian and D'Oria, Marco and Rode, Michael and Khlifi, Slaheddine}, title = {Assessment of climate change impact on discharge of the lakhmass catchment (Northwest Tunisia)}, series = {Water}, volume = {14}, journal = {Water}, number = {14}, publisher = {MDPI}, address = {Basel}, issn = {2073-4441}, doi = {10.3390/w14142242}, pages = {17}, year = {2022}, abstract = {The Mediterranean region is increasingly recognized as a climate change hotspot but is highly underrepresented in hydrological climate change studies. This study aims to investigate the climate change effects on the hydrology of Lakhmass catchment in Tunisia. Lakhmass catchment is a part of the Medium Valley of Medjerda in northwestern Tunisia that drains an area of 126 km(2). First, the Hydrologiska Byrans Vattenbalansavdelning light (HBV-light) model was calibrated and validated successfully at a daily time step to simulate discharge during the 1981-1986 period. The Nash Sutcliffe Efficiency and Percent bias (NSE, PBIAS) were (0.80, +2.0\%) and (0.53, -9.5\%) for calibration (September 1982-August 1984) and validation (September 1984-August 1986) periods, respectively. Second, HBV-light model was considered as a predictive tool to simulate discharge in a baseline period (1981-2009) and future projections using data (precipitation and temperature) from thirteen combinations of General Circulation Models (GCMs) and Regional Climatic Models (RCMs). We used two trajectories of Representative Concentration Pathways, RCP4.5 and RCP8.5, suggested by the Intergovernmental Panel on Climate Change (IPCC). Each RCP is divided into three projection periods: near-term (2010-2039), mid-term (2040-2069) and long-term (2070-2099). For both scenarios, a decrease in precipitation and discharge will be expected with an increase in air temperature and a reduction in precipitation with almost 5\% for every +1 degrees C of global warming. By long-term (2070-2099) projection period, results suggested an increase in temperature with about 2.7 degrees C and 4 degrees C, and a decrease in precipitation of approximately 7.5\% and 15\% under RCP4.5 and RCP8.5, respectively. This will likely result in a reduction of discharge of 12.5\% and 36.6\% under RCP4.5 and RCP8.5, respectively. This situation calls for early climate change adaptation measures under a participatory approach, including multiple stakeholders and water users.}, 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{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{DrozdovAllisonShpritsetal.2022, author = {Drozdov, Alexander and Allison, Hayley J. and Shprits, Yuri and Usanova, Maria E. and Saikin, Anthony and Wang, Dedong}, title = {Depletions of Multi-MeV Electrons and their association to Minima in Phase Space Density}, series = {Geophysical research letters}, volume = {49}, journal = {Geophysical research letters}, number = {8}, publisher = {American Geophysical Union}, address = {Washington}, issn = {0094-8276}, doi = {10.1029/2021GL097620}, pages = {11}, year = {2022}, abstract = {Fast-localized electron loss, resulting from interactions with electromagnetic ion cyclotron (EMIC) waves, can produce deepening minima in phase space density (PSD) radial profiles. Here, we perform a statistical analysis of local PSD minima to quantify how readily these are associated with radiation belt depletions. The statistics of PSD minima observed over a year are compared to the Versatile Electron Radiation Belts (VERB) simulations, both including and excluding EMIC waves. The observed minima distribution can only be achieved in the simulation including EMIC waves, indicating their importance in the dynamics of the radiation belts. By analyzing electron flux depletions in conjunction with the observed PSD minima, we show that, in the heart of the outer radiation belt (L* < 5), on average, 53\% of multi-MeV electron depletions are associated with PSD minima, demonstrating that fast localized loss by interactions with EMIC waves are a common and crucial process for ultra-relativistic electron populations.}, language = {en} } @article{DahmHeimannMetzetal.2021, author = {Dahm, Torsten and Heimann, Sebastian and Metz, Malte and Isken, Marius Paul}, title = {A self-similar dynamic rupture model based on the simplified wave-rupture analogy}, series = {Geophysical journal international / the Royal Astronomical Society, the Deutsche Geophysikalische Gesellschaft and the European Geophysical Society}, volume = {225}, journal = {Geophysical journal international / the Royal Astronomical Society, the Deutsche Geophysikalische Gesellschaft and the European Geophysical Society}, number = {3}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggab045}, pages = {1586 -- 1604}, year = {2021}, abstract = {The investigation of stresses, faults, structure and seismic hazards requires a good understanding and mapping of earthquake rupture and slip. Constraining the finite source of earthquakes from seismic and geodetic waveforms is challenging because the directional effects of the rupture itself are small and dynamic numerical solutions often include a large number of free parameters. The computational effort is large and therefore difficult to use in an exploratory forward modelling or inversion approach. Here, we use a simplified self-similar fracture model with only a few parameters, where the propagation of the fracture front is decoupled from the calculation of the slip. The approximative method is flexible and computationally efficient. We discuss the strengths and limitations of the model with real-case examples of well-studied earthquakes. These include the M-w 8.3 2015 Illapel, Chile, megathrust earthquake at the plate interface of a subduction zone and examples of continental intraplate strike-slip earthquakes like the M-w 7.1 2016 Kumamoto, Japan, multisegment variable slip event or the M-w 7.5 2018 Palu, Indonesia, supershear earthquake. Despite the simplicity of the model, a large number of observational features ranging from different rupture-front isochrones and slip distributions to directional waveform effects or high slip patches are easy to model. The temporal evolution of slip rate and rise time are derived from the incremental growth of the rupture and the stress drop without imposing other constraints. The new model is fast and implemented in the open-source Python seismology toolbox Pyrocko, ready to study the physics of rupture and to be used in finite source inversions.}, language = {en} } @article{KrstulovićRosaBiedermannetal.2021, author = {Krstulović, Marija and Rosa, Angelika D. and Biedermann, Nicole and Irifune, Tetsuo and Wilke, Max}, title = {Structural changes in aluminosilicate glasses up to 164 GPa and the role of alkali, alkaline earth cations and alumina in the densification mechanism}, series = {Chemical geology : official journal of the European Association for Geochemistry}, volume = {560}, journal = {Chemical geology : official journal of the European Association for Geochemistry}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0009-2541}, doi = {10.1016/j.chemgeo.2020.119980}, pages = {14}, year = {2021}, abstract = {Pressure induced structural changes in silicate melts have a great impact on their physico-chemical properties and hence on their behaviour in the deep Earth's interior. In order to gain a deeper understanding we have studied the densification mechanism in multicomponent aluminosilicate glasses (albitic and albit-diopside composition) by means of extended X-ray absorption fine structure spectroscopy coupled to a diamond anvil cell up to 164 GPa. We have monitored the structural modifications from the network-former Ge as well as the network-modifier Sr. Notably, we tracked the evolution of Ge-O and Sr-O bond lengths (RGe-O, RSr-O) and their coordination number with pressure. We show that RGe-O increases strongly up to about 32 GPa, whereas RSr-O increases only slightly up to similar to 26 GPa. We assign these extensions to the increase of the coordination number from 4 to 6 (Ge) and from similar to 6 to at least 9 (Sr). Upon further compression RGe-O and RSr-O exhibit a continuous decrease to the highest probed pressure. These bond contractions, notably of RGe-O, that are continuous and exceed the one observed in pure SiO2 and GeO2, reflect a higher structural flexibility of multi-component glasses compared to those simple systems. Particularly, the high fraction of non-bridging oxygen atoms due to the presence of Na, Sr, Ca, Mg in the studied glasses, favours the simple compression of the highly-coordinated polyhedra of Si and Ge at pressure greater than 30 GPa. This is in strong contrast to pure oxides where cation polyhedral distortions govern the densification mechanism of the glass. The results of this study demonstrate that low field-strength alkali and alkaline earth cations, ubiquitous in deep Earth's melts, have a profound influence on the densification mechanism of glasses. Our results provide important constrains for interpreting the observed low velocity anomalies at the Earth's core-mantle boundary that have been, beyond others, referred to the presence of high-density melts. The hypothesis that non-buoyant melts at the Earth's core-mantle boundary can be formed by peculiar structural transformations in melts leading to higher coordination numbers compared to their crystalline equivalents is not supported from the present observations. The present results rather suggest that if velocity anomalies are to be explained by melts, these likely have considerable differences in chemical composition to the surrounding crystalline phase assemblage.}, 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{ZhuCottonKwaketal.2021, author = {Zhu, Chuanbin and Cotton, Fabrice and Kwak, Dong-Youp and Ji, Kun and Kawase, Hiroshi and Pilz, Marco}, title = {Within-site variability in earthquake site response}, series = {Geophysical journal international}, volume = {229}, journal = {Geophysical journal international}, number = {2}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggab481}, pages = {1268 -- 1281}, year = {2021}, abstract = {The within-site variability in site response is the randomness in site response at a given site from different earthquakes and is treated as aleatory variability in current seismic hazard/risk analyses. In this study, we investigate the single-station variability in linear site response at K-NET and KiK-net stations in Japan using a large number of earthquake recordings. We found that the standard deviation of the horizontal-to-vertical Fourier spectral ratio at individual sites, that is single-station horizontal-to-vertical spectral ratio (HVSR) sigma sigma(HV,s), approximates the within-site variability in site response quantified using surface-to-borehole spectral ratios (for oscillator frequencies higher than the site fundamental frequency) or empirical ground-motion models. Based on this finding, we then utilize the single-station HVSR sigma as a convenient tool to study the site-response variability at 697 KiK-net and 1169 K-NET sites. Our results show that at certain frequencies, stiff, rough and shallow sites, as well as small and local events tend to have a higher sigma(HV,s). However, when being averaged over different sites, the single-station HVSR sigma, that is sigma(HV), increases gradually with decreasing frequency. In the frequency range of 0.25-25 Hz, sigma(HV) is centred at 0.23-0.43 in ln scales (a linear scale factor of 1.26-1.54) with one standard deviation of less than 0.1. sigma(HV) is quite stable across different tectonic regions, and we present a constant, as well as earthquake magnitude- and distance-dependent sigma(HV) models.}, language = {en} } @inproceedings{SimonPapadopoulosGuillemoteauetal.2021, author = {Simon, Fran{\c{c}}ois-Xavier and Papadopoulos, Nikos and Guillemoteau, Julien and Oikonomou, Dimitris and Simirdanis, Kleanthis}, title = {Multi-frequency loop electromagnetic system measurement on shallow offshore archaeological site of Oulos}, series = {ArcheoSciences : revue d'arch{\´e}om{\´e}trie / Groupe des M{\´e}thodes Pluridisciplinaires Contribuant {\`a} l'Arch{\´e}ologie (GMPCA)}, volume = {45}, booktitle = {ArcheoSciences : revue d'arch{\´e}om{\´e}trie / Groupe des M{\´e}thodes Pluridisciplinaires Contribuant {\`a} l'Arch{\´e}ologie (GMPCA)}, number = {1}, publisher = {Presses Universitaires de Rennes}, address = {Rennes}, issn = {1960-1360}, doi = {10.4000/archeosciences.9690}, pages = {215 -- 218}, year = {2021}, language = {en} } @article{DaempflingMielkeKoellneretal.2022, author = {D{\"a}mpfling, Helge L. C. and Mielke, Christian and Koellner, Nicole and Lorenz, Melanie and Rogass, Christian and Altenberger, Uwe and Harlov, Daniel E. and Knoper, Michael}, title = {Automatic element and mineral detection in thin sections using hyperspectral transmittance imaging microscopy (HyperTIM)}, series = {European journal of mineralogy}, volume = {34}, journal = {European journal of mineralogy}, number = {3}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {0935-1221}, doi = {10.5194/ejm-34-275-2022}, pages = {275 -- 284}, year = {2022}, abstract = {In this study we present a novel method for the automatic detection of minerals and elements using hyperspectral transmittance imaging microscopy measurements of complete thin sections (HyperTIM). This is accomplished by using a hyperspectral camera system that operates in the visible and near-infrared (VNIR) range with a specifically designed sample holder, scanning setup, and a microscope lens. We utilize this method on a monazite ore thin section from Steenkampskraal (South Africa), which we analyzed for the rare earth element (REE)-bearing mineral monazite ((Ce,Nd,La)PO4), with high concentrations of Nd. The transmittance analyses with the hyperspectral VNIR camera can be used to identify REE minerals and Nd in thin sections. We propose a three-point band depth index, the Nd feature depth index (NdFD), and its related product the Nd band depth index (NdBDI), which enables automatic mineral detection and classification for the Nd-bearing monazites in thin sections. In combination with the average concentration of the relative Nd content, it permits a destruction-free, total concentration calculation for Nd across the entire thin section.}, language = {en} } @article{BarbosaCoelhoGusmaoetal.2022, author = {Barbosa, Luis Romero A. and Coelho, Victor Hugo R. and Gusmao, Ana Claudia V. L. F. and Fernandes, Lucila A. E. and da Silva, Bernardo B. and Galvao, Carlos de O. and Caicedo, Nelson O. L. and da Paz, Adriano R. and Xuan, Yunqing and Bertrand, Guillaume F. and Melo, Davi de C. D. and Montenegro, Suzana M. G. L. and Oswald, Sascha and Almeida, Cristiano das N.}, title = {A satellite-based approach to estimating spatially distributed groundwater recharge rates in a tropical wet sedimentary region despite cloudy conditions}, series = {Journal of hydrology}, volume = {607}, journal = {Journal of hydrology}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0022-1694}, doi = {10.1016/j.jhydrol.2022.127503}, pages = {15}, year = {2022}, abstract = {Groundwater recharge (GWR) is one of the most challenging water fluxes to estimate, as it relies on observed data that are often limited in many developing countries. This study developed an innovative water budget method using satellite products for estimating the spatially distributed GWR at monthly and annual scales in tropical wet sedimentary regions despite cloudy conditions. The distinctive features proposed in this study include the capacity to address 1) evapotranspiration estimations in tropical wet regions frequently overlaid by substantial cloud cover; and 2) seasonal root-zone water storage estimations in sedimentary regions prone to monthly variations. The method also utilises satellite-based information of the precipitation and surface runoff. The GWR was estimated and validated for the hydrologically contrasting years 2016 and 2017 over a tropical wet sedimentary region located in North-eastern Brazil, which has substantial potential for groundwater abstraction. This study showed that applying a cloud-cleaning procedure based on monthly compositions of biophysical data enables the production of a reasonable proxy for evapotranspiration able to estimate groundwater by the water budget method. The resulting GWR rates were 219 (2016) and 302 (2017) mm yr(-1), showing good correlations (CC = 0.68 to 0.83) and slight underestimations (PBIAS =-13 to-9\%) when compared with the referenced estimates obtained by the water table fluctuation method for 23 monitoring wells. Sensitivity analysis shows that water storage changes account for +19\% to-22\% of our monthly evaluation. The satellite-based approach consistently demonstrated that the consideration of cloud-cleaned evapotranspiration and root-zone soil water storage changes are essential for a proper estimation of spatially distributed GWR in tropical wet sedimentary regions because of their weather seasonality and cloudy conditions.}, language = {en} } @incollection{CardenasSchernthanner2022, author = {C{\´a}rdenas, Aura and Schernthanner, Harald}, title = {The role of livestock wastes in clean energy}, series = {Handbook of waste biorefinery}, booktitle = {Handbook of waste biorefinery}, editor = {Jacob-Lopes, Eduardo and Queiroz Zepka, Leila and Costa Depr{\´a}, Mariany}, publisher = {Springer}, address = {Cham}, isbn = {978-3-031-06561-3}, doi = {10.1007/978-3-031-06562-0_12}, pages = {337 -- 343}, year = {2022}, abstract = {Agricultural production worldwide has been increasing in the last decades at a very fast pace and with it the waste generation. Livestock activities are one of the largest producers of residues in the agricultural sector and contribute greatly to climate change. The present chapter gives an introduction and an in-depth analysis of the waste management of livestock for the conversion in a circular agriculture and economy based on research and experience in the sector conducted in the last decades. The conversion of animal waste into energy generation is an opportunity for farmers to obtain additional economic benefits, while contributing to the environment by preventing the release of GHGs into the atmosphere. The use of animal waste for energy generation through anaerobic digestion is a progressive technique and is being widely accepted in Europe, where Germany is the leading country in the use of biogas plants for energy production among others in the European Union. Economically speaking, the livestock industry faces the challenge of converting its production into a clean and more profitable production. The goal of this chapter is to analyze the economic benefit as well as the environmental contribution and future challenges of the use of livestock waste in the biorefineries sector from different perspectives, based on an intensive literature review. This review is accompanied by a geospatial analysis component, mapping biogas reactor hotspots and clusters in Germany, by means of methods of spatial statistics as analysis methods as kernel density estimations (KDE) and K-means clustering, based on volunteer geographic data. The applied methods easily can be transferred to other regions and allow a quick macroscopic overview over existing biogas reactors; furthermore, an identification of cluster and hotspots with a high biogas potential, that in a subsequent step can be analyzed in depth in larger scales.}, language = {en} }