@article{FlovenzWangHersiretal.2022, author = {Fl{\´o}venz, {\´O}lafur G. and Wang, Rongjiang and Hersir, Gylfi P{\´a}ll and Dahm, Torsten and Hainzl, Sebastian and Vassileva, Magdalena and Drouin, Vincent and Heimann, Sebastian and Isken, Marius Paul and Gudnason, Egill {\´A}. and {\´A}g{\´u}stsson, Kristj{\´a}n and {\´A}g{\´u}stsd{\´o}ttir, Thorbj{\"o}rg and Hor{\´a}lek, Josef and Motagh, Mahdi and Walter, Thomas R. and Rivalta, Eleonora and Jousset, Philippe and Krawczyk, Charlotte M. and Milkereit, Claus}, title = {Cyclical geothermal unrest as a precursor to Iceland's 2021 Fagradalsfjall eruption}, series = {Nature geoscience}, volume = {15}, journal = {Nature geoscience}, number = {5}, publisher = {Nature Research}, address = {Berlin}, issn = {1752-0894}, doi = {10.1038/s41561-022-00930-5}, pages = {397 -- 404}, year = {2022}, abstract = {Understanding and constraining the source of geodetic deformation in volcanic areas is an important component of hazard assessment. Here, we analyse deformation and seismicity for one year before the March 2021 Fagradalsfjall eruption in Iceland. We generate a high-resolution catalogue of 39,500 earthquakes using optical cable recordings and develop a poroelastic model to describe three pre-eruptional uplift and subsidence cycles at the Svartsengi geothermal field, 8 km west of the eruption site. We find the observed deformation is best explained by cyclic intrusions into a permeable aquifer by a fluid injected at 4 km depth below the geothermal field, with a total volume of 0.11 ± 0.05 km3 and a density of 850 ± 350 kg m-3. We therefore suggest that ingression of magmatic CO2 can explain the geodetic, gravity and seismic data, although some contribution of magma cannot be excluded.}, language = {en} } @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}, 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{SharmaHainzlZoeller2023, author = {Sharma, Shubham and Hainzl, Sebastian and Z{\"o}ller, Gert}, title = {Seismicity parameters dependence on main shock-induced co-seismic stress}, series = {Geophysical journal international}, volume = {235}, journal = {Geophysical journal international}, number = {1}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggad201}, pages = {509 -- 517}, year = {2023}, abstract = {The Gutenberg-Richter (GR) and the Omori-Utsu (OU) law describe the earthquakes' energy release and temporal clustering and are thus of great importance for seismic hazard assessment. Motivated by experimental results, which indicate stress-dependent parameters, we consider a combined global data set of 127 main shock-aftershock sequences and perform a systematic study of the relationship between main shock-induced stress changes and associated seismicity patterns. For this purpose, we calculate space-dependent Coulomb Stress (\& UDelta;CFS) and alternative receiver-independent stress metrics in the surrounding of the main shocks. Our results indicate a clear positive correlation between the GR b-value and the induced stress, contrasting expectations from laboratory experiments and suggesting a crucial role of structural heterogeneity and strength variations. Furthermore, we demonstrate that the aftershock productivity increases nonlinearly with stress, while the OU parameters c and p systematically decrease for increasing stress changes. Our partly unexpected findings can have an important impact on future estimations of the aftershock hazard.}, 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} }