ALBERT

Description: We performed seismic ambient noise tomography to investigate the shallow crustal structure around the Ivrea geophysical body (IGB) in the Ivrea-Verbano Zone (IVZ). We achieved higher resolution with respect to previous tomographic works covering the Western Alps, by processing seismic data collected by both permanent and temporary seismic networks (61 broad-band seismic stations in total). This included IvreaArray, a temporary, passive seismic experiment designed to investigate the IVZ crustal structure. Starting from continuous seismic ambient noise recordings, we measured and inverted the dispersion of the group velocity of surface Rayleigh waves (fundamental mode) in the period range 4–25 s. We obtained a new, 3-D vS model of the IVZ crust via the stochastic neighbourhood algorithm (NA), with the highest resolution between 3 to 40 km depth. The fast and shallow shear wave velocity anomaly associated with the IGB presents velocities of 3.6 km s−1 directly at the surface, in remarkable agreement with the location of the exposed lower-to-middle crustal and mantle outcrops. This suggests a continuity between the surface geological observations and the subsurface geophysical anomalies. The fast IGB structure reaches vS of 4 km s−1 at 20–25 km depth, at the boundary between the European and Adriatic tectonic plates, and in correspondence with the earlier identified Moho jump in the same area. The interpretation of a very shallow reaching IGB is further supported by the comparison of our new results with recent geophysical investigations, based on receiver functions and gravity anomaly data. By combining the new geophysical constraints and the geological observations at the surface, we provide a new structural interpretation of the IGB, which features lower crustal and mantle rocks at upper crustal depths. The comparison of the obtained vS values with the physical properties from laboratory analysis of local rock samples suggests that the bulk of the IGB consists of a combination of mantle peridotite, ultramafic and lower crustal rocks, bound in a heterogeneous structure. These new findings, based on vS tomography, corroborate the recent interpretation for which the Balmuccia peridotite outcrops are continuously linked to the IGB structure beneath. The new outcomes contribute to a multidisciplinary framework for the interpretation of the forthcoming results of the scientific drilling project DIVE. DIVE aims at probing the lower continental crust and its transition to the mantle, with two ongoing and one future boreholes (down to 4 km depth) in the IVZ area, providing new, complementary information on rock structure and composition across scales. In this framework, we constrain the upper crustal IGB geometries and lithology based on new evidence for vS, connecting prior crustal knowledge to recent active seismic investigations.

Description: Published

Description: 1089–1105

Description: OST1 Alla ricerca dei Motori Geodinamici

Description: JCR Journal

Description: On the 9th of January 2020, an Mw 6.4 strike-slip earthquake took place north of the Asian margin of the Bering Sea. The earthquake occurred within the known reverse-right-lateral active fault zone, called Khatyrka–Vyvenka, which transverses the Koryak Highland from SE to NW and is thought to be a surface manifestation of the Asian portion of either the Bering plate boundary or the northern edge of the Alaskan stream. No other strong earthquake has ever been recorded in this remote uninhabited area and the few existing seismic stations provide poor quality earthquake locations.We adopt SAR interferometry (InSAR) technique to define an improved location of the Koryak 2020 earthquake and constrain the seismic source. The analysis of the 2020 event revealed a previously unknown active fault of left-lateral kinematics that is possibly hidden and strikes NWtransversely to the Khatyrka–Vyvenka fault zone. Although several mechanisms could account for left-lateral kinematics of this fault, we propose that the structure is part of a more extended NW fault structure, that formed in pre-neotectonic times and has played a role of a pre-existing rheological discontinuity. This revived NW structure together with a similar structure located easterly, so far aseismic, make the plate/stream boundary segmented, step-like in plan view. The step-like boundary geometry may be the result of internal transform deformation of a rigid plate, but it is better explained by deflections of the Alaskan stream edge at local crustal asperities, which are pre-Cenozoic terrains.

Description: Published

Description: 1412–1421

Description: OST2 Deformazione e Hazard sismico e da maremoto

Description: JCR Journal

Description: This article has been accepted for publication in Geophysical Journal International ©:The Author(s) 2023. Published by Oxford University Press on behalf of the Royal Astronomical Society. All rights reserved.Uploaded in accordance with the publisher's self-archiving policy. All rights reserved.

Description: Estimation of local seismic response plays a key role in assessing local seismic hazard and particularly in the design of shaking scenarios. Modelling local seismic response involves knowing of the shear wave velocity (Vs) and quality factor (Qs) profiles for the site in question. The many techniques that have been developed to assess Vs in surface deposits produce reliable measurements of Vs , but these rarely correspond to direct measurements of Qs . The latter is often considered through damping measures from laboratory tests on small-scale soil samples, which can provide information primarily on intrinsic attenuation, neglecting the contribution of scattering effects. In this paper, using seismic recordings obtained at the surface and in boreholes at 100 m depth, we estimate an average value of Qs of some characteristic alluvial deposits of the Po Plain (northern Italy). Data come from a microseismic network which sampled an almost uniform lithology in the central Po Plain and consisted of three surface and four borehole stations with an interstation distance of about 2 km. The average value of Qs of the shallowest 100 m of the sedimentary strata, Qs100, is estimated by considering: (1) the high-frequency attenuation of seismic waves due to propagation through the corresponding stratigraphy and (2) the interference between incident and surface-reflected waves observed at borehole stations. We parametrize the first through k0_100, the difference between the values of the spectral decay parameter kappa (k) estimated at the surface and at the boreholes depth, respectively. We use the second in order to compute Vs100, the time-averaged Vs referred to the uppermost 100 m stratigraphy. We obtain: k0_100 = (11 ± 3) ms, Vs100 = (309 ± 11) m s −1 and Qs100 = 31 ± 10. At the surface, the estimated values of the site-specific kappa, k0, are found to range from 75 to 79 ms. As expected, these results are in good agreement with studies performed in other sites characterized by sandy or clayey lithologies, and can be usefully used in site response analysis at sites where the rigidity is mainly controlled by lithostatic pressure.

Description: Comune di Minerbio (grant: “Sperimentazione ILG Minerbio”; grant number: 0913.010).

Description: Published

Description: 2075–2094

Description: OST2 Deformazione e Hazard sismico e da maremoto

Description: JCR Journal

Description: Accurate quantification of seismic activity in volcanic regions is an important asset for im- proving hazard and risk assessment. This is especially true for densely populated areas, as in the case of Etna volcano (Southern Italy). There, the volcanic hazard is amplified by the seismic risk of acti ve faults, especiall y on the eastern flank of the volcano. In such a context, it is common to rely on moment magnitude ( M W ) to characterize seismicity and monitor the energy released during an eruption. In this study, we calculate the moment-based magnitude ( M W ) for selected seismic data sets, using different approaches in distinct magnitude ranges to cover the widest possible range of magnitude that characterizes Etna’s seismicity . Specifically , we computed the M W from a data set of moment tensor solutions of earthquakes that occurred in the magnitude range 3.4 ≤M L ≤4.8 during 2005–2020; we created a data set of seismic moment and associated M W for earthquakes 1.0 ≤M L 〈 3.4 obtained by analysing source spectra; we fine-tuned two relationships, for shallow and deep earthquakes, to obtain M W from response spectra. Finally, we calibrated a specific relationship between M W and M L for the Etna area earthquakes in the range 1.0 ≤M L ≤4.8. All the empirical relationships obtained in this study can be applied in real-time analysis of the seismicity to provide fast and robust information on the released seismic energy.

Description: INGV-DPC 2012- 2021 agreement; B2 DPC-INGV 2019-2021 project; IMPACT Department strategic project ; ‘Project PE0000005–RETURN (NRRP)

Description: Published

Description: 2520-2534

Description: OST2 Deformazione e Hazard sismico e da maremoto

Description: JCR Journal

Description: This study presents a new robust statistical framework, in which to measure relative differences, or deviations from a hypothetical reference value, of Gutenberg-Richter b-value. Moreover, it applies this method to recent seismicity in Italy, to find possible changes of earthquake magnitude distribution in time and space. The method uses bootstrap techniques, which have no prior assumptions about the distribution of data, keeping their basic features. Excluding Central Italy, no significative b-value variation is found, revealing that the frequency-magnitude distribution exponent is substantially stable or that data are not able to reveal hidden variations. Considering the small size of examined magnitude samples, we cannot definitively decide if the higher b-values in Central Italy, consistently founded by all applied tests, have a physical origin or result from a statistical bias. In any case, they indicate short-lived excursions which have a temporary nature and, therefore, cannot be associated solely to spatial variations in tectonic framework. Both the methodological issues and the results of the application to seismicity in Italy show that a correct assessing of b-value changes requests appropriate statistics, that accurately quantify the low accuracy and precision of b-value estimation for small magnitude samples.

Description: Published

Description: 729–740

Description: OST4 Descrizione in tempo reale del terremoto, del maremoto, loro predicibilità e impatto

Description: JCR Journal

Description: This article has been accepted for publication in Geophysical Journal International ©:The Author(s) 2023. Published by Oxford University Press on behalf of the Royal Astronomical Society. All rights reserved.Uploaded in accordance with the publisher's self-archiving policy. All rights reserved.

Description: We present the results from a fully unconstrained moment tensor inversion of induced seismic events in a complex and high seismic hazard region (Val d’Agri basin, Southern Italy). The study area hosts two well-documented cases of induced microseismicity linked to (i) a wastewater injection well of a giant oilfield (the largest in onshore Europe), and (ii) severe seasonal level changes of an artificial lake. In order to gather information on the non-doublecouple components of the source and to better understand the rupture mechanisms, we analyse seismic events recorded during daily injection tests in the disposal well. The computed moment tensors have significant non-double-couple components that correlate with the well-head injection pressure. The injection parameters strongly influence the rupture mechanism that can be interpreted as due to the opening/closing of a fracture network inside a fault zone of a pre-existing thrust fault. For the case of the reservoir-induced seismicity, no direct correlations are observed with the loading/unloading of the reservoir.

Description: Published

Description: 1617–1627

Description: OST3 Vicino alla faglia

Description: JCR Journal

Description: Serpentinites are polymineralic rocks distributed almost ubiquitously across the globe in active tectonic regions. Magnetite-rich serpentinites are found in the low-strain domains of serpen- tinite shear zones, which act as potential sites of nucleation of unstable slip. To assess the potential of earthquake nucleation in these materials, we investigate the link between me- chanical properties and fabric of these rocks through a suite of laboratory shear experiments. Our experiments were done at room temperature and cover a range of normal stress and slip velocity from 25 to 100 MPa and 0.3 to 300 μm s −1 , respecti vel y. We show that magnetite-rich serpentinites are ideal materials since they display strong sensitivity to the loading rate and are susceptible to nucleation of unstable slip, especially at low forcing slip velocities. We also aim at the integration of mechanical and microstructural results to describe the underlying mechanisms that produce the macroscopic behaviour. We show that mineralogical composi- tion and mineral structure dictates the coexistence of two deformation mechanisms leading to stable and unstable slip. The weakness of phyllosilicates allows for creep during the interseis- mic phase of the laboratory seismic cycle while favouring the restoration of a load-bearing granular framework, responsible of the nucleation of unstable events. During dynamic slip, fault zone shear fabric determines the mode of slip, producing either asymmetric or Gaussian slip time functions for either fast or slow events. We report rate/state friction parameters and integrate our mechanical data with microstructural observations to shed light on the mech- anisms dictating the complexity of laborator y ear thquakes. We show that mineralogical and fabric heterogeneities control fault slip behaviour.

Description: Published

Description: 1778–1797

Description: OST3 Vicino alla faglia

Description: JCR Journal

Description: This article has been accepted for publication in Geophysical Journal International ©:The Author(s) 2023. Published by Oxford University Press on behalf of the Royal Astronomical Society. All rights reserved.Uploaded in accordance with the publisher's self-archiving policy. All rights reserved.

Description: We report on about 20 yr of relative gravity measurements, acquired on Mt. Somma–Vesuvius volcano in order to investigate the hydrological and volcano-tectonic processes controlling the present-day activity of the volcano. The retrieved long-term field of time gravity change (2003–2022) shows a pattern essentially related to the subsidence, which have affected the central part of the volcano, as detected by the permanent GNSS network and InSAR data. After reducing the observations for the effect of vertical deformation, no significant residuals are found, indicating no significant mass accumulation or loss within the volcanic system. In the north-western sector of the study area, at the border of the volcano edifice, however, significant residual positive gravity changes are detected which are associated to ground-water rebound after years of intense exploitation of the aquifers. On the seasonal timescale, we find that stations within the caldera rim are affected by the seasonal hydrological effects, while the gravity stations at the base of the Vesuvius show a less clear correlation. Furthermore, within the caldera rim a multiyear gravity transient is detected with an increase phase lasting about 4 yr followed by a slower decrease phase. Analysis of rain data seem to exclude a hydrological origin, hence, we hypothesize a deeper source related to the geothermal activity, which can be present even if the volcano is in a quiescent state. We infer the depth and volume of the source by inverting the spatial pattern of the gravity field at the peak of the transient. A volume of fluids of 9.5 × 107 m3 with density of 1000 kg m−3 at 2.3 km depth is capable to fit reasonably well the observations. To explain the gravity transient, simple synthetic models are produced, that simulate the ascent of fluids from a deep reservoir up to the depth of 2.3 km and a successive diffusion within the carbonate aquifer hosting the geothermal system. The whole process appears to not significantly affect the seismicity rate and the deformation of the volcano. This study demonstrates the importance of a 4-D gravity monitoring of a volcano to understand its complex gravity signals that cover different spatial and temporal scales. Discriminating the different contributions that mix up in the observed gravity changes, in particular those due to hydrologic/anthropogenic activities form those due to the geothermal dynamics, is fundamental for a complete and reliable evaluation of the volcano state.

Description: Published

Description: 1565–1580

Description: OSV2: Complessità dei processi vulcanici: approcci multidisciplinari e multiparametrici

Description: JCR Journal

Description: This article has been accepted for publication in Geophysical Journal International ©:The Author(s) 2023. Published by Oxford University Press on behalf of the Royal Astronomical Society. All rights reserved.Uploaded in accordance with the publisher's self-archiving policy. All rights reserved.

Description: In a recent work, we applied the every earthquake a precursor according to scale (EEPAS) probabilistic model to the pseudo-prospective forecasting of shallow earthquakes with magni- tude M 5.0 in the Italian region. We compared the forecasting performance of EEPAS with that of the epidemic type aftershock sequences (ETAS) forecasting model, using the most recent consistency tests developed within the collaboratory for the study of earthquake predictabil- ity (CSEP). The application of such models for the forecasting of Italian target earthquakes seems to show peculiar characteristics for each of them. In particular, the ETAS model showed higher performance for short-term forecasting, in contrast, the EEPAS model showed higher forecasting performance for the medium/long-term. In this work, we compare the performance of EEPAS and ETAS models with that obtained by a deterministic model based on the occur- rence of strong foreshocks (FORE model) using an alarm-based approach. We apply the two rate-based models (ETAS and EEPAS) estimating the best probability threshold above which we issue an alarm. The model parameters and probability thresholds for issuing the alarms are calibrated on a learning data set from 1990 to 2011 during which 27 target earthquakes have occurred within the analysis region. The pseudo-prospective forecasting performance is as- sessed on a validation data set from 2012 to 2021, which also comprises 27 target earthquakes. Tests to assess the forecasting capability demonstrate that, even if all models outperform a purely random method, which trivially forecast earthquake proportionally to the space–time occupied by alarms, the EEPAS model exhibits lower forecasting performance than ETAS and FORE models. In addition, the relative performance comparison of the three models demonstrates that the forecasting capability of the FORE model appears slightly better than ETAS, but the difference is not statistically significant as it remains within the uncertainty level. However, truly prospective tests are necessary to validate such results, ideally using new testing procedures allowing the analysis of alarm-based models, not yet available within the CSEP.

Description: Published

Description: 1541–1551

Description: OST4 Descrizione in tempo reale del terremoto, del maremoto, loro predicibilità e impatto

Description: JCR Journal

Description: To understand the seismic hazard of a subduction zone, it is necessary to know the geometry, location and mechanical characteristics of the interplate boundary below which an oceanic plate is thrust downward. By considering the azimuthal dependence of converted P-to-S (Ps) amplitudes in receiver functions, we have detected the interplate boundary in the Makran subduction zone, revealing significant seismic anisotropy at the base of the accretionary wedge above the slab before it bends down beneath the Jaz Murian basin. This anisotropic feature aligns with a zone of reduced seismic velocity and a high primary/secondary wave velocity ratio (Vp/Vs), as documented in previous studies. The presence of this low-velocity highly anisotropic layer at the base of the accretionary wedge, likely representing a low-strength shear zone, could possibly explain the unusually wide accretionary wedge in Makran. Additionally, it may impact the location and width of the locked zone along the interplate boundary.

Description: Iranian National Science Foundation (INSF)

Description: Published

Description: 64-74

Description: OST1 Alla ricerca dei Motori Geodinamici

Description: JCR Journal

Description: This article has been accepted for publication in Geophysical Journal International ©:The Author(s) 2023. Published by Oxford University Press on behalf of the Royal Astronomical Society. All rights reserved.Uploaded in accordance with the publisher's self-archiving policy. All rights reserved.

Description: The Every Earthquake a Precursor According to Scale (EEPAS) forecasting model is a space– time point-process model based on the precursory scale increase (ψ ) phenomenon and associated predictive scaling relations. It has been previously applied to New Zealand, Cal- ifornia and Japan earthquakes with target magnitude thresholds varying from about 5–7. In all previous application, computations were done using the computer code implemented in Fortran language by the model authors. In this work, we applied it to Italy using a suite of computing codes completely rewritten in Matlab. We first compared the two software codes to ensure the convergence and adequate coincidence between the estimated model parameters for a simple region capable of being analysed by both software codes. Then, using the rewritten codes, we optimized the parameters for a different and more complex polygon of analysis using the Homogenized Instrumental Seismic Catalogue data from 1990 to 2011. We then perform a pseudo-prospective forecasting experiment of Italian earthquakes from 2012 to 2021 with Mw ≥ 5.0 and compare the forecasting skill of EEPAS with those obtained by other time in- dependent (Spatially Uniform Poisson, Spatially Variable Poisson and PPE: Proximity to Past Earthquakes) and time dependent [Epidemic Type Aftershock Sequence (ETAS)] forecasting models using the information gain per active cell. The preference goes to the ETAS model for short time intervals (3 months) and to the EEPAS model for longer time intervals (6 months to 10 yr).

Description: Published

Description: 1681–1700

Description: OST4 Descrizione in tempo reale del terremoto, del maremoto, loro predicibilità e impatto

Description: JCR Journal

Description: Thermobarometry provides a critical means of assessing locations of magma storage and dynamics in the lead-up to volcanic eruptions and crustal growth. A common approach is to utilise minerals that have compositions sensitive to changes in pressure and/or temperature, such as clinopyroxene, which is ubiquitous in mafic to intermediate magmas. However, clinopyroxene thermobarometry may carry significant uncertainty and require an appropriate equilibrium melt composition. In addition, the degree of magma undercooling (ΔT) affects clinopyroxene composition and zoning, with common sector zoning potentially obfuscating thermobarometry results. Here, we use a set of crystallisation experiments on a primitive trachybasalt from Mt. Etna (Italy) at ΔT = 25–233 °C, P = 400–800 MPa, H2O = 0–4 wt % and fO2 = NNO + 2, with clinopyroxene crystals defined by Al-rich zones (prisms and skeletons) and Al-poor zones (hourglass and overgrowths) to assess common equilibrium models and thermobarometric approaches. Under the studied conditions, our data suggest that the commonly applied Fe–Mg exchange (cpx-meltKdFe–Mg) is insensitive to increasing ΔT and may not be a reliable indicator of equilibrium. The combined use of DiHd (CaMgSi2O6 + CaFeSi2O6) and EnFs (Mg2Si2O6 + Fe2Si2O6) models indicate the attainment of equilibrium in both Al-rich and Al-poor zones for almost all investigated ΔT. In contrast, CaTs (CaAl2SiO6) and CaTi (CaTiAl2O6) models reveal substantial deviations from equilibrium with increasing ΔT, particularly in Al-rich zones. We postulate that this reflects slower diffusion of Al and Ti in the melt compared with Ca and Mg and recommend the concurrent application of these four models to evaluate equilibrium between clinopyroxene and melt, particularly for sector-zoned crystals. Thermobarometers calibrated with only isothermal–isobaric experiments closely reproduce experimental P–T at low ΔT, equivalent to natural phenocrysts cores and sector-zoned mantles. Models that also consider decompression experiments are most accurate at high ΔT and are therefore suitable for outermost phenocryst rims and groundmass microlites. Recent machine learning approaches reproduce P–T conditions across all ΔT conditions. Applying our experimental constraints to sector-zoned microphenocrysts and groundmass microlites erupted during the 1974 eccentric eruption at Mt. Etna, we highlight that both hourglass and prism sectors are suitable for thermobarometry, given that equilibrium is sufficiently tested for. The combination of DiHd, EnFs, CaTs and CaTi models identifies compositions closest to equilibrium with the bulk melt composition, and results in smaller differences in P–T calculated for hourglass and prism sectors compared with applying only DiHd and EnFs equilibrium models. This provides a framework to assess crystallisation conditions recorded by sector-zoned clinopyroxene crystals in mafic alkaline settings.

Description: Published

Description: egad074

Description: OSV2: Complessità dei processi vulcanici: approcci multidisciplinari e multiparametrici

Description: JCR Journal

Description: The joint ESA/NASA Mass-change And Geosciences International Constellation (MAGIC) has the objective to extend time-series from previous gravity missions, including an improvement of accuracy and spatio-temporal resolution. The long-term monitoring of Earth’s gravity field carries information on mass change induced by water cycle, climate change and mass transport processes between atmosphere, cryosphere, oceans and solid Earth. MAGIC will be composed of two satellite pairs flying in different orbit planes. The NASA/DLR-led first pair (P1) is expected to be in a near-polar orbit around 500 km of altitude; while the second ESA-led pair (P2) is expected to be in an inclined orbit of 65°–70° at approximately 400 km altitude. The ESA-led pair P2 Next Generation Gravity Mission shall be launched after P1 in a staggered manner to form the MAGIC constellation. The addition of an inclined pair shall lead to reduction of temporal aliasing effects and consequently of reliance on de-aliasing models and post-processing. The main novelty of the MAGIC constellation is the delivery of mass-change products at higher spatial resolution, temporal (i.e. subweekly) resolution, shorter latency and higher accuracy than the Gravity Recovery and Climate Experiment (GRACE) and Gravity Recovery and Climate Experiment Follow-On (GRACE-FO). This will pave the way to new science applications and operational services. In this paper, an overview of various fields of science and service applications for hydrology, cryosphere, oceanography, solid Earth, climate change and geodesy is provided. These thematic fields and newly enabled applications and services were analysed in the frame of the initial ESA Science Support activities for MAGIC. The analyses of MAGIC scenarios for different application areas in the field of geosciences confirmed that the double-pair configuration will significantly enlarge the number of observable mass-change phenomena by resolving smaller spatial scales with an uncertainty that satisfies evolved user requirements expressed by international bodies such as IUGG. The required uncertainty levels of dedicated thematic fields met by MAGIC unfiltered Level-2 products will benefit hydrological applications by recovering more than 90 per cent of the major river basins worldwide at 260 km spatial resolution, cryosphere applications by enabling mass change signal separation in the interior of Greenland from those in the coastal zones and by resolving small-scale mass variability in challenging regions such as the Antarctic Peninsula, oceanography applications by monitoring meridional overturning circulation changes on timescales of years and decades, climate applications by detecting amplitude and phase changes of Terrestrial Water Storage after 30 yr in 64 and 56 per cent of the global land areas and solid Earth applications by lowering the Earthquake detection threshold from magnitude 8.8 to magnitude 7.4 with spatial resolution increased to 333 km.

Description: Published

Description: 1288–1308

Description: JCR Journal

Description: In near-Earth space, a large population of high-energy electrons are trapped by Earth’s magnetic field. These energetic electrons are trapped in the regions called Earth’s ring current and radiation belts. They are very dynamic and show a very strong dependence on solar wind and geomagnetic conditions. These energetic electrons can be dangerous to satellites in the near-Earth space. Therefore, it is very important to understand the mechanisms which drive the dynamics of these energetic electrons. Wave-particle interaction is one of the most important mechanisms. Among the waves that can be encountered by the energetic electrons when they move around our Earth, whistler mode chorus waves can cause both acceleration and the loss of energetic electrons in the Earth's radiation belts and ring current. Using more than 5 years of wave measurements from NASA’s Van Allen Probe mission, Wang et al (2019) developed chorus wave models which depend on magnetic local time (MLT), Magnetic Latitude (MLat), L-shell, and geomagnetic condition index Kp. To quantify the effect of chorus waves on energetic electrons, we calculated the bounce-averaged quasi-linear diffusion coefficients using the chorus wave model developed by Wang et al (2019) and extended to higher latitudes according to Wang and Shprits (2019). Using these diffusion coefficients, we calculated the lifetime of the electrons with an energy range from 1 keV to 2 MeV. In each MLT, we calculate the lifetime for each energy and L-shell using two different methods according to Shprits et al (2007) and Albert and Shprits (2009). We make the calculated electron lifetime database available here. Please notice that the chorus wave model by Wang et al (2019) is valid when Kp 〈= 6. If the user wants to use this lifetime database for Kp 〉6, please be careful and contact the authors.

Description: This publication provides the codes produced for the article "Temporally dynamic carbon dioxide and methane emission factors for rewetted peatlands. Nature Communications Earth and Environment" by Aram Kalhori, Christian Wille, Pia Gottschalk, Zhan Li, Josh Hashemi, Karl Kemper, and Torsten Sachs (https://doi.org/10.1038/s43247-024-01226-9). In the article, the authors estimate the cumulative GHG emissions of a rewetted peatland in Germany using the long-term ecosystem flux measurements. They observe a source-to-sink transition of annual carbon dioxide (CO2) fluxes and decreasing trend of methane (CH4) emissions. This software is written in R and MATLAB. Running the codes ([R files and .m files](Code)) and loading the data files ([CSV files and .mat files](Data)) requires the pre-installation of [R and RStudio] (https://posit.co/downloads/) and ([MATLAB]. The RStudio 2022.07.2 Build 576 version has been used for the R scripts. The land cover classification work was performed in QGIS, v.3.16.11-Hannover. Data were analyzed in both MATLAB and R and plots created with R (R Core Development Team 2020) in RStudio®.

Description: The data provided here is an exemplary dataset for the flux site Zarnekow from one year (2018). The complete dataset that is needed to run the codes for all the years can be obtained from the European Fluxes Database Cluster under site ID DE-Zrk (Sachs et al., 2016) or provided upon request. This repository is intended to provide the necessary MATLAB and R code to reproduce the results by Kalhori et al. (2024). The data are provided as zip folder containing (1) a csv file with associated definition of variables and units (file: 2023-004_Kalhori-et-al_README_2018_units.txt), (2) a shapefile (file: 2023-004_Kalhori-et-al_2018_LAiV_DOP.shp) and (3) a Geotiff (file: 2023-004_Kalhori-et-al_2018_LAiV_DOP.tiff).

Description: This data set contains the results from a 2023 GFZ Innovative Research Expedition project to explore for natural hydrogen gas (H2) occurrences in the NW Pyrenean foreland, near the town of Biarritz in France. The data represent in-situ measurements of soil and spring water gas, as well as in-situ spring water property measurements, complemented with laboratory analysis results of gas contents and noble gas isotopic compositions of gas and spring water samples collected during the expedition. This GFZ Innovative Research Expedition was inspired by previous exploration efforts in the region by Lefeuvre et al. (2021, 2022). These authors detected elevated concentrations of natural H2 gas in the soil and interpreted this natural H2 to be derived from serpentinizing mantle rocks below the Pyrenees. The main aims of this expedition were the following: (1) in-situ measuring soil gas contents and taking soil gas samples for laboratory analysis at a site near the town of Peyrehorade in the NW of the general study area of Lefeuvre et al. (2021), thus improving the soil gas data coverage along the NW end of the North Pyrenean Frontal Thrust (NPFT); (2) taking gas samples from degassing springs (or water samples from non-degassing springs to be degassed in the lab) in the general Lefeuvre et al. (2021) study area for additional laboratory analysis of gas contents and noble gas isotopic compositions, which may be indicative of (deep) gas origins; and (3) performing a detailed soil gas analysis by means of a portable mass spectrometer at Sauveterre-de-Béarn, a site along the NPFT where Lefeuvre et al. (2022) measured elevated concentrations of natural H2 in the soil. Furthermore, we also measured the properties of the visited springs (temperature, pH, conductivity) while on site, and performed additional in-situ soil gas measurements from manual drillholes. Details on the measurement and sampling methods, on the laboratory analyses, as well as the results of these measurements and analyses are provided in the data description file The expedition involved six field days in July 2023, during which a total of 26 sites were visited. These sites were selected for their vicinity near a major geological contact or fault zone that could have facilitated upward circulation of gas or (thermal) water from the (deep) subsurface (i.e., potentially from the mantle).

Description: This dataset comprises event peak flows, representing extreme floods at 516 stations in Germany. The data generation process involves several key steps. Initially, observed rainfall events associated with 10 historical flood disasters from 1950 to 2021 are undergone spatial shifts. These shifts involve three distances (20, 50, and 100 km) and eight directions (North, Northeast, East, Southeast, South, Southwest, West, Northwest), resulting in 24 counterfactual precipitation events. Including the factual (no shift) event, a total of 25 distinct shifting events are considered. Subsequently, these shifted fields are used as atmospheric forcing for a mesoscale hydrological model (mHM) set up and calibrated for the entire Germany. The model produces daily stream flows across its domain, from which the event peak flows are derived. This dataset is expected to provide a valuable resource for analyzing and modeling the dynamics extreme flood events in Germany.

Description: The dataset is the basis for describing a 60-year-long evolution of groundwater dynamics and thermal field in the North German Basin beneath the Federal State of Brandenburg (NE Germany), covering the period between 1953 and 2014 with monthly increments. It was produced by one-way coupling of a near-surface distributed hydrologic model to a 3D basin-scale thermohydraulic groundwater model with the goal of investigating feedbacks between climate-driven forcing (in terms of time- and space-varying recharge and temperature), basin-scale geology, and topographic gradients. Modeled pressure and temperature distributions are validated against published groundwater level and temperature time series from observation wells. Our results indicate the spatio-temporal extent of the groundwater system subjected to nonlinear interactions between local geological variability and climate conditions. The dataset comprises of input files and scripts required to run the groundwater model in GOLEM and output files from the transient thermo-hydraulic simulations in EXODUS format. The input and output data is organized as separate archived folders (*.gz format).

Description: The data set was collected to identify hydrological processes and their evolution over it time. It consists of several individual files in tabstop delimeted text format. The data set contains the data obtained from deuterium and brilliant blue tracer experiments at two chronosequence studies in the glacier forefield of the Stone Glacier and the Griessfirn in the central Alps, Switzerland. Each chronosequence consisted of four moraines of different ages (from 30 to 13500 years). At each forefield sprinkling experiments with deuterium and dye tracer experiments with blue dye (Brilliant Blue) were conducted on three plots per moraine. The moraines at the forefield of the Stone Glacier developed from siliceous parent material and at the forefield of the Griessfirn from calcareous parent material. Data from the siliceous forefield are marked with (S) and data from the calcareous forefield are marked with (C). The data set consist of soil moisture time series and soil water isotope profiles of the sprinkling experiments with deuterium, as well as trinary images of stained vertical subsurface flow paths from the dye tracer experiment. The individual plots per moraine are distinguished via their position relative to one another on the moraine (left, middle, and right, looking upslope). The plots used for the sprinkling experiments were located in close vicinity to the plots used for the dye tracer experiments. For the sprinkling experiments with deuterium each plot (4m x 6m) per age class was equipped with 6 soil moisture sensors. Three of these sensors were installed as a sensor profile at one side of the plot about one meter downslope from the upper plot boundary. The sensors were installed at 10, 30, and 50 cm soil depth. On the other side of the plot, two sensors were placed in 10 cm depth, one opposite to the sensor profile and the second sensor one meter upslope from the lower plot boundary. The sixth sensor was placed at 10 cm depth in the center of the plot. The plots were irrigated on three consecutive days with three different irrigation intensities and deuterium concentrations. Per forefield, the soil moisture data are listed in one file per age class. The file contains for each plot, the time stamp and the soil moisture values of the 6 sensors.

Description: This data repository for the Southern Caribbean and NW South America contains a 3D thermal model computed down to 75 km depth, the modelled hypocentral temperatures and geothermal gradients at the locations of crustal earthquakes, and the crustal seismogenic depths calculated from earthquake statistics, as well as the associated modelled temperatures. We used the uppermost 75 km of the gravity-constrained structural and density model of Gómez-García et al. (2020, 2021) to derive the 3D thermal configuration of the study area (5°-15° N, 63°-82° W). A steady-state approach was followed, in which upper and lower boundary conditions were set to run the thermal calculations using the software GOLEM (Cacace & Jacquey, 2017; Jacquey & Cacace, 2017). A catalogue of earthquakes occurred within the study area and surroundings was compiled from public sources. In the database archived here, we provide data of the best located crustal earthquakes within the boundaries of this area, from January 1980 to June 2021. Earthquakes below the magnitude of completeness, or with poorly determined depths, were disregarded. Earthquakes were deemed crustal if their hypocentres were located between the topo-bathymetry from the GEBCO relief (Weatherall et al., 2015) and the Moho depth from the GEMMA model (Reguzzoni & Sampietro, 2015). We computed the crustal seismogenic depth as the 90th and 95th percentiles (D90 and D95), respectively, of the crustal hypocentral depths. These percentiles were mapped on a latitude-longitude grid, using for each grid node at least the 20 closest earthquakes as sample. The hypocentral temperatures, the geothermal gradient at the earthquake locations, and the temperatures at the D90 and D95 surfaces were calculated from the lithospheric-scale thermal model. For more details about the modelling approach and interpretation of the results, we kindly ask the reader to refer to the main publication: Gomez-Garcia et al. (2024).

Description: This data set includes the results of high-resolution digital elevation models (DEM) and digital image correlation (DIC) analysis applied to analogue modelling experiments. Twenty generic analogue models are extended on top of a rubber sheet. Two benchmark experiments are also reported. Detailed descriptions of the experiments can be found in Liu et al. (submitted) to which this data set is supplement. The data presented here are visualized as topography and the horizontal cumulative surface strain (principal strain and slip rake).

Description: The EU funded project CRM-geothermal aims to establish an overview of the potential for critical raw materials (CRM) in geothermal fluids across the EU and third countries (Ref). Within this framework, the geothermal sites of Tuzla, Seferihisar and Dikili in eastern Turkey have been visited in March 2023. To estimate the potential of CRM at these sites, a comprehensive sampling program was performed. Rock samples (drill gravel) of the production borehole and scaling from gas-water separators were obtained. Furthermore, sampling of geothermal fluids (gas and brine) and precipitates (salt) along the production line was performed. Here, the results of the geochemical analyses of solid sample materials (drill gravel, scales and salt) are presented. All analyses were performed in the ElMiE-Lab (Elements and Minerals of the Earth Laboratory) at German Research Centre for Geosciences Potsdam, Germany (https://labinfrastructure.geo-x.net/laboratories/8). For their major and minor element compositions, bulk samples of drill gravel and scales were analyzed with XRF and ICP-MS, respectively. Salt precipitates were analyzed for dry loss and mineral composition using XRD

Description: In autumn 2022, an expedition to Tanzania was undertaken within the framework of the research project “CRM-geothermal” and Scintific Priority Program (SPP) 2238 “Dynamics of Ore Metal Enrichment”. Within „CRM-geothermal“ we are looking for an environmentally friendly co-production of critical raw materials together with the provision of geothermal energy. In the EARS, high levels of rare earth elements (REE), Sr, Ba and Mg are expected in waters and solids in areas with alkaline volcanic rocks, while other critical elements, including helium, have been sought in other localities. In particular, the eastern branch is the most juvenile sector and has increased geothermal potential related to hot fluids migrating along permeable faults. Tanzania was crossed from north to south, along the eastern arm of the EARS, to collect gas, water, rock and sediment samples associated with natural hot springs, lakes and vents. On site, physical and chemical parameters were measured in-situ and documented together with the geology, infrastructure and the domestic use of the hot site. In the south, existing drill sites and geothermal development areas were visited and gas and water samples were taken from boreholes and rocks sampled from drill cores. The survey covered 13 sites, from Lake Natron in the north to Lake Malawi in the south (see map).

Description: This dataset is the result of an experimental series that was carried out in September/October 2022 at GFZ German Research Centre for Geosciences, Potsdam, Germany to observe biosorption of lead under extreme conditions. Synthetic solutions, simulating the geothermal fluids from the Heemskerk geothermal power plant were were prepared in 30 ml glass vials (Rotalibo screw neck ND24 EPA). To prepare the stock solutions, sodium chloride (NaCl, 99.8 %, Cellpure, Merck, DE) was added at 265 g/L and Pb(II), in form of lead nitrate (Pb(NO3 )2 , Merck, DE), at 1 g/L to ultrapure water. To assess the impact of acetic acid on lead biosorption, two treatments were done: one without acetic acid and one where acetic acid (100 %, Merck, DE) was added at 60 mg/L. Finally, dead biomass of the fungus Penicillium citrinum was added in the samples at a concentration of 4 g/L (Wahab et al., 2017). The samples were incubated in an autoclave at a pressure of 8 bars on a rotative shaker. The temperature was set at 25 °C, 60 °C or 98 °C with three contact times (1, 2 and 3 h). All treatments were performed in triplicates. For each treatment, two controls without biomass were done. Control samples without the addition of NaCl were done in duplicate, at 25 °C and for 2 h. After incubation, samples were filtered through a 0.22 µm nitrocellulose filter (Sartorius Stedim Biotech, FR) to separate the biomass from the liquid. The biomass on the filters was dried for 24 h at 45 °C before being scraped from the filter and kept in a Falcon tube at room temperature.

Description: The Turkey heat flow database includes several research articles obtained from the catalogue of The Global Heat Flow Data Assessment Project conducted by the International Heat Flow Commission (IHFC; www.ihfc-iugg.org). The presented database contains 725 heat-flow determinations compiled from 9 different publications generated between 1991-2023 reported within Turkey. For the reporting and sorting of the database, the structure documented by Fuchs et al. (2023) is followed. Within this dataset, 98% of the entries represent continental heat-flow data (onshore), while the remaining 2% correspond to marine data (offshore). 88% of the reported heat flow values were obtained via direct temperature measurements, while the remaining data (12%) were estimated from indirect Curie depth temperature calculations.

Description: In diesem Bericht werden die durch das GFZ Potsdam am 29. und 30. November 2023 durchgeführte bohrlochgeophysikalische Messungen in den Bohrungen Gt Khn 1/88 und Gt Khn 2/87 in Karlshagen (Mecklenburg-Vorpommern) dokumen-tiert. Die Messungen wurden mit dem Ziel der Gewinnung hochaufgelöster und un-gestörter Temperatur-Tiefen-Profile durchgeführt. Die Stillstandszeiten seit Erstel-lung liegen bei mehreren Jahrzehnten; jene seit letzter Befahrung bei fünfzehn Jahren, weshalb von ungestörten Gebirgstemperaturen ausgegangen werden kann. In der Bohrung Gt Khn 2/87 wurde bei 1786,5 m Teufe eine Temperatur von 57,8 °C, welches einem mittleren Temperaturgradienten von 27,8 °C/km entspricht, ge-messen. Die Bohrung Gt Khn 1/88 konnte bis zu einer Teufe von 325,1 m befahren werden, die gemessene Temperatur betrug 16,2 °C, der entsprechende mittlere ge-othermische Gradient beträgt ca. 23,6 °C/km. This report documents the borehole geophysical logging performed by GFZ Potsdam in the Gt Khn 1/88 and Gt Khn 2/87 boreholes in Karlshagen (Mecklenburg-Western Pomerania) on the 29th and 30th of November 2023. The measurements were conducted to achieve high-resolution and undisturbed temperature-depth pro-files. The shut-in times since the boreholes were drilled are several decades; the shut-in time since last activities in the boreholes are in the order of 15 years. There-fore, undisturbed formation temperatures can be expected in the boreholes. In the Gt Khn 2/87 borehole, a temperature of 57.8 °C was measured at a depth of 1786.5 m, which corresponds to an average temperature gradient of 27.8 °C/km. The Gt Khn 1/88 borehole could be logged to a depth of 325.1 m and the measured temperature at this depth was 16.2 °C, corresponding to an average geothermal gradient of approx. 23.6 °C/km.

Description: The data publication contains the compilation of global heat-flow data by the International Heat Flow Commission (IHFC; www.ihfc-iugg.org) of the International Association of Seismology and Physics of the Earth's Interior (IASPEI). The presented data update release 2024 contains data generated between 1939 and 2024 and constitutes the second intermediate update benefiting from the global collaborative assessment and quality control of the Global Heat Flow Database running since May 2021 (http://assessment.ihfc-iugg.org). The data release comprises new original heat-flow data published since April 2023 (the update 2023). It contains 91,182 heat-flow data from 1,586 publications. 57% of the reported heat-flow values are from the continental domain (n ~ 54,553), while the remaining 43% are located in the oceanic domain (n ~ 36,692).

Description: The dataset contains source parameters of acoustic emission (AE) events recorded during triaxial friction (stick-slip) experiments performed on the Westerly Granite sample WgN05. In addition we provide raw waveform data of AE events recorded in triggered mode with a network of 16 AE sensors. Basic seismic catalog associated with the stick-slip experiment contains origin time, hypocentral location in local Cartesian coordinate system of the sample (with associated uncertainties), and AE-derived magnitude. In addition, for a subset of AEs we provide full moment tensors. This catalog include information on fault parameters (strike, dip and rake of the two nodal planes), percentage of isotropic, compensated linear vector dipole and double-couple components of the full moment tensor, P, T, B axes orientations in the coordinate system of the sample, uncertainty assessment, as well as the six independent moment tensor components. Finally, we provide a time series of axial stress values as presented in the Kwiatek et al. (2023) as well as the coordinates of the AE sensors. The catalog and parametric data is supplemented with the raw waveform recordings stored in HDF5 format from 16 acoustic emission sensors placed on the surface of the sample.

Description: Teleseismic back-projection imaging has emerged as a powerful tool for understanding the rupture propagation of large earthquakes. However, its application often suffers from artifacts related to the receiver array geometry. We developed a teleseismic back-projection technique that can accommodate data from multiple arrays. Combined processing of P and pP waveforms may further improve the resolution. The method is suitable for defining arrays ad-hoc to achieve a good azimuthal distribution for most earthquakes. We present a catalog of short-period rupture histories (0.5-2.0 Hz) for all earthquakes from 2010 to 2022 with Mw {greater than or equal to} 7.5 and depth less than 200 km (56 events). The method provides automatic estimates of rupture length, directivity, speed, and aspect ratio, a proxy for rupture complexity. We obtained short-period rupture length scaling relations that are in good agreement with previously published relations based on estimates of total slip. Rupture speeds were consistently in the sub-Rayleigh regime for thrust and normal earthquakes, whereas a tenth of strike-slip events propagated at supershear speeds. Many rupture histories exhibited complex behaviors, e.g., rupture on conjugate faults, bilateral propagation, and dynamic triggering by a P wave. For megathrust earthquakes, ruptures encircling asperities were frequently observed, with down-dip, up-dip, and balanced patterns. Although there is a preference for short-period emissions to emanate from central and down-dip parts of the megathrust, emissions up-dip of the main asperity are more frequent than suggested by earlier results. The data are presented as follows (and described in detail in the associated README): SUPPORTING DATA SET S1 (2024-001_Vera-et-al_Supporting-Data-S1.zip) This Data Set (S1) consists of *.bp files containing (1) short-period earthquake rupture patterns, (2) energy radiated maps, and (3) source time functions derived from back-projections (0.5-2.0 Hz). The Data Set S1 includes 56 folders, representing 56 processed earthquakes between 2010 and 2022 with a moment magnitude (Mw) greater than or equal to 7.5 and a depth less than 200 km. These folders are labeled in the format YYYYMMDDhhmm_EVENT_NAME_REGION (UTC) in *.bp format. SUPPORTING DATA SET S2 (2024-001_Vera-et-al_Supporting-Data-S2.csv) This Data Set (S2) comprises a *.csv file containing earthquake source information used in the back-projection and the resulting ruptur...

Description: To enhance the EU's economic autonomy, feasible options for local sourcing of critical raw materials that would allow for shorter supply routes along with ethical and responsible value chains are under contemplation. Social acceptance of mining in Europe is, however, low, and the establishment of new mining sites faces strong public opposition. Therefore, innovative solutions for the production of primary raw materials need to be developed. A new idea for raw material extraction is the extraction of essential elements from geothermal fluids. Deep geothermal fluids, increasingly used for energy production, often contain high concen-trations of dissolved ions and gases in commercially interesting concentrations. The EU-funded project CRM-geothermal aims to develop new technologies to extract these highly relevant elements, including helium, during geothermal production cycles. In this way, an environmentally friendly and socially acceptable exploration and exploitation method could be deployed. One aim of the CRM-geothermal project is to gain an overview of the actual quantities of critical raw materials in various geothermal fluids in Europe by taking and analyzing fluid samples. In Turkey for instance, classical high enthalpy (volcanic) systems exist, which are representative for many geothermal areas worldwide. The sites are located at the edges of tectonic plates and close to areas undergoing volcanic activity. The brines are mixed with seawater and circulate in the deeper crust. The data publication contains analyses results of three gas samples from Tuzla, two samples from Seferihisar geothermal power plant and one sample from the Dikili geothermal field in Turkey, taken in 2023 as part of the CRM-geothermal project.

Description: This data publication provides a European assessment of building exposure, organized country-by-country. The dataset provides information about the number of buildings; the number of occupants; structural information and structural costs of buildings per geographical area. The main purpose of this data collection is risk assessment for natural hazards, however it can be used by anyone in need of a building exposure dataset. The data holds information about single buildings, with global estimates of built-up area on 10m x 10m pixels and exposure information per district. All OpenStreetMap (OSM) buildings existing in an OSM excerpt from 1 July 2023, 00:00 UTC (OpenStreetMap contributors, 2023), all buildings from the Global ML Building Footprint (GMLBF, Microsoft, 2023) dataset have been processed and for each building the occupancy type and number of stories have been identified based on data in OSM, such as land use and points of interest. The Global Human Settlement Built-up Characteristics 2022A Layer has been used as initial distribution of built area (Pesaresi, 2022). Aggregated exposure information, including the structural information and the number of occupants, stems the ESRM20 (Crowley et al., 2020). The resulting dataset is distributed per country as an SQLite/SpatiaLite database. Each database contains three tables and one view. The database is organized around three key concepts, that each have their own table. An Entity is a geographical unit that contains exposure. In this dataset, the entities are tiles in a multi-resolution grid, according to the Quad tree structure (Finkel & Bentley, 1974), with the tiles projected using the Web Mercator projection (EPSG:3857). The zoom-level of the Quadkeys inside the grid varies from level-15 to level-18, depending on the number of buildings inside each tile to preserve privacy-sensitive information. Practically, the size of the tiles varies between around 100m x 100m and 1km x 1km. Each entity consists of one or more Assets, defining the number of buildings of a particular structural type and their population and structural value. The structural type is described using a taxonomy string, describing for example structural properties, occupancy type and the expected number of stories. The exact definition of a taxonomy that is used in this dataset is described in the GEM Building Taxonomy v2.0 (Brzev et al., 2013). On top of the tables, one key view has been defined too. A view is essentially a query on the table that give some insights into the data. The `key_values_per_tile` provides the total number of buildings, total number of occupants at night and total structural costs summed over all assets in one tile entity.

Description: This data publication provides the results of the investigations and measurements of thermal rock properties conducted on site in the Tournemire field laboratory and at the Thermal Petrophysics Lab at GFZ. The thermal characterization of the clayey Jurassic (Upper Toarcian, ca. 180 My old) is contributing to the site characterization of the Tournemire Underground Research Lab (URL), located in Southern France. This URL is installed in a former railway tun-nel to better understand the physical processes resulting from thermal and hydrau-lic loading in a small fault zone in a highly consolidated shale formation (Bonnelye et al., 2023). At the Tournemire site, faults and fractures of different sizes extend from the surface (sedimentary cover) to the crystalline basement. At one specific gallery (Gallery East 03) installed in the former tunnel, thermally controlled in-situ fluid injection experiments are scheduled on a strike-slip fault zone outcropping at the URL (Bonnelye et al., 2023). In 2022, we visited the URL for baseline characteri-zation of thermal properties and to study the heterogeneity of the clay-dominated formation. Therefore, we took the chance to collect data and samples for a laborato-ry measurement campaign and to measure thermal conductivity in-situ in the tun-nel wall of Gallery East 03. The thermal data shall provide the baseline for the pa-rameterization of future numerical 3D models to better understand the thermal-hydraulic processes related to the experiment. This data publication provides the results of the investigations and measurements conducted on-site in the field la-boratory and at the Thermal Petrophysics Lab at GFZ.

Description: This article has been accepted for publication in Geophysical Journal International ©:The Author(s) 2023. Published by Oxford University Press on behalf of the Royal Astronomical Society. All rights reserved.Uploaded in accordance with the publisher's self-archiving policy. All rights reserved.

Description: The south-eastern sector of the Mount Etna, Italy, is characterized by numerous active faults, in particular the Belpasso–Ognina lineament, the Tremestieri–San Gregorio–Acitrezza fault, the Trecastagni fault and the Fiandaca–Nizzeti fault including the Timpe Fault System. Their activity is the result of both volcanism and tectonics. Here, we analyse the ground deformation occurred from 2016 to 2019 across those active faults by using the GNSS data acquired at 22 permanent stations and 35 campaign points observed by the Etna Observatory (INGV) and by the University of Catania. We also use the time-series of line of sight displacement of permanent scatterers SENTINEL-1 A-DInSAR obtained by using the P-SBAS tool of the ESA GEP-TEP (Geohazards Thematic Exploitation Platform) service. We discriminate the contributions of the regional tectonic strain, the inflations, the deflations of the volcano and the gravitational sliding in order to analyse the deformation along the faults of the south-eastern flank of Etna. The shallow and destructive Mw = 4.9 earthquake of 2018 December 26 occurred within the studied area two days after a dyke intrusion, that propagated beneath the centre of the volcano accompanied by a short eruption. Both GNSS and InSAR time-series document well those events and allow to investigate the post-seismic sliding across the faults of south-eastern flank. We analyse the slow slip events (SSE) that are observed in the GNSS and InSAR time-series in the vicinity of the Acitrezza fault. We quantify and discuss the tectonic origin of the Belpasso–Ognina lineament that we interpreted as a tear fault.

Description: Published

Description: 664–682

Description: OST3 Vicino alla faglia

Description: JCR Journal

Description: Ground shaking, whether it is due to natural or induced earthquakes, has always been a matter of concern since it correlates with structural/non-structural damage and can culminate in human anxiety. Industrial activities such as water injection, gas sequestration and waste fluid disposals, promote induced seismicity and consequent ground shaking that can hinder ongoing activities. Therefore, keeping in mind the importance of timely evaluation of a seismic hazard and its mitigation for societal benefits, the present study proposes specifically designed ground-motion prediction equations (GMPEs) from induced earthquakes in the St. Gallen geothermal area, Switzerland. The data analysed in this study consist of 343 earthquakes with magnitude −1.17 ≤ ML, corr ≤ 3.5 and hypocentral distance between 4 and 15 km. The proposed study is one of the first to incorporate ground motions from negative magnitude earthquakes for the development of GMPEs. The GMPEs are inferred with a two-phase approach. In the first phase, a reference model is obtained by considering the effect of source and medium properties on the ground motion. In the second phase the final model is obtained by including a site/station effect. The comparison between the GMPEs obtained in the present study with GMPEs developed for the other induced seismicity environments highlights a mismatch that is ascribed to differences in regional seismic environment and local site conditions of the respective regions. This suggests that, when dealing with induced earthquakes, GMPEs specific for the study should be inferred and used for both monitoring purposes and seismic hazard analyses.

Description: Published

Description: 820–832

Description: 5T. Sismologia, geofisica e geologia per l'ingegneria sismica

Description: JCR Journal

Description: © The Author(s), 2022. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Hudson, A. R., Peters, D. P. C., Blair, J. M., Childers, D. L., Doran, P. T., Geil, K., Gooseff, M., Gross, K. L., Haddad, N. M., Pastore, M. A., Rudgers, J. A., Sala, O., Seabloom, E. W., & Shaver, G. Cross-site comparisons of dryland ecosystem response to climate change in the US long-term ecological research network. Bioscience, 72(9), (2022): 889–907, https://doi.org/10.1093/biosci/biab134.

Description: Long-term observations and experiments in diverse drylands reveal how ecosystems and services are responding to climate change. To develop generalities about climate change impacts at dryland sites, we compared broadscale patterns in climate and synthesized primary production responses among the eight terrestrial, nonforested sites of the United States Long-Term Ecological Research (US LTER) Network located in temperate (Southwest and Midwest) and polar (Arctic and Antarctic) regions. All sites experienced warming in recent decades, whereas drought varied regionally with multidecadal phases. Multiple years of wet or dry conditions had larger effects than single years on primary production. Droughts, floods, and wildfires altered resource availability and restructured plant communities, with greater impacts on primary production than warming alone. During severe regional droughts, air pollution from wildfire and dust events peaked. Studies at US LTER drylands over more than 40 years demonstrate reciprocal links and feedbacks among dryland ecosystems, climate-driven disturbance events, and climate change.

Description: Funding was provided by the USDA-ARS SCINet Big Data Project (grant no. 0500–00093–001–00-D), and the National Science Foundation US LTER Program to New Mexico State University for the Jornada Basin (grant no. DEB 20–25166), Kansas State University for the Konza Prairie (grant no. DEB 2025849), the Kellogg Biological Station (grant no. DEB 1832042), Cedar Creek Ecosystem Science Reserve (grants no. DEB-1234162 and no. DEB-1831944), ARC (grant no. DEB-1637459), MCM (grant no. OPP-1637708), CAP (grant no. DEB-1832016), and SEV (grant no. DEB-1655499). Support was also provided by the Minnesota Supercomputer Institute and the University of Minnesota, Michigan State University AgBioResearch.

Description: Author Posting. © The Author(s), 2022. This is the author's version of the work. It is posted here by permission of Oxford University Press for personal use, not for redistribution. The definitive version was published in Geophysical Journal International 231(2),(2022): 1434–1445, https://doi.org/10.1093/gji/ggac257.

Description: Makran subduction zone is very active with ∼38 mm yr−1 convergence rate and has experienced great earthquakes in the past. The latest great earthquake of 1945 Mw 8.1 event also triggered a large tsunami and led to ∼4000 casualties. However, due to incomplete historical seismicity records and poor modern instrumentation, earthquake mechanism, co-seismic slip and tsunami characteristics in Makran remain unclear. On 2017 February 17, an Mw 6.3 earthquake rattled offshore Pasni of Pakistan in the eastern Makran, marking the largest event after the 1945 Mw 8.1 earthquake with good geodetic and geophysical data coverage. We use a combination of seismicity, multibeam bathymetry, seismic profile, InSAR measurements and tide-gauge observation to investigate the seismogenic structure, co-seismic deformation, tsunami characteristics of this event and its implication for future major earthquakes. Our results indicate that (1) the earthquake occurred on the shallow-dipping (3°–4°) megathrust; (2) the megathrust co-seismically slipped 15 cm and caused ∼2–4 cm ground subsidence and uplift at Pasni; (3) our tsunami modelling reproduces the observed 5-cm-high small tsunami waveforms. The Pasni earthquake rupture largely overlaps the 1945 slip patch and disturbs the west and east megathrust segments that have not ruptured yet at least since 1765. With such stress perturbation and possible stress evolution effect from the 1945 earthquake, the unruptured patches may fail in the future. This study calls for more preparedness in mitigating earthquake and associated hazards in the eastern Makran.

Description: his study is financially supported by the National Natural Science Foundation of China (Nos. 42076059, 41890813, 41976066 and 41976064), the Key Special Project for Introduced Talents Team of Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou) (No. GML2019ZD0205), Chinese Academy of Sciences (Nos. Y4SL021001, QYZDY-SSW-DQC005, 131551KYSB20200021, ISEE2021PY03, 133244KYSB20180029 and E1SL3C02), Guangdong Basic and Applied Basic Research Foundation (No. 2021B1515020098) and China–Pakistan Joint Research Centre on Earth Sciences.

Description: The Doldrums Megatransform System (~7–8°N, Mid-Atlantic Ridge) shows a complex architecture including four intra-transform ridge segments bounded by five active transform faults. Lower crustal rocks are exposed along the Doldrums and Vernadsky transform walls that bound the northernmost intra-transform ridge segment. The recovered gabbros are characterized by variably evolved chemical compositions, ranging from olivine gabbros to gabbronorites and oxide gabbros, and lack the most primitive gabbroic endmembers (troctolites, dunites). Notably, the numerous recovered gabbronorites show up to 20 vol. % of coarse-grained orthopyroxene. Although covariations in mineral and bulk-rock chemical compositions of the olivine and oxide gabbros define trends of crystallization from a common parental melt, the gabbronorites show elevated light over heavy rare earth elements (LREE/HREE) ratios in both bulk-rock and mineral compositions. These features are not consistent with a petrological evolution driven solely by fractional crystallization, which cannot produce the preferential enrichments in highly incompatible elements documented in the orthopyroxene-bearing lithologies. We suggest that gabbronorites crystallized from evolved melts percolating and partly assimilating a pre-existing olivine gabbro matrix. Saturation in orthopyroxene and selective enrichments in LREE relative to M-HREE are both triggered by an increase in assimilated crystal mass, which ranges from negligible in the oxide-gabbros to abundant in the gabbronorites. This melt–rock reaction process has been related to lateral melt migration beneath ridge-transform intersections, where variably evolved melts injected from the peripheral parts of the melting region towards the transform zone may interact with a gabbroic crystal mush to form abundant oxide-bearing gabbronoritic associations.

Description: Published

Description: egac086

Description: 3A. Geofisica marina e osservazioni multiparametriche a fondo mare

Description: JCR Journal

Description: The REHEATFUNQ Python package helps to work with the (residual) scatter of surface heat flow even in small regions. REHEATFUNQ uses a stochastic model for regional aggregate heat flow distributions (RAHFD), that is, the collected set of heat flow measurements within a region marginalized to the heat flow dimension. The stochastic model is used in a Bayesian analysis that (1) yields a posterior estimate of the RAHFD which captures the range of heat flow within the analysis region, and (2) quantifies the magnitude of a surface heat flow anomaly within the region, for instance through the generating frictional power. The stochastic model underlying REHEATFUNQ views heat flow data, uniformly sampled across the region of interest, as a random variable. A gamma distribution is used as a model for this random variable and information from the global data set of Lucazeau (2019) is introduced by means of a conjugate prior (Miller, 1980). The detailed science behind the model is described in Ziebarth et al. (202X). The analysis by Ziebarth et al. (202X) can be reproduced through the Jupyter notebooks contained in the subdirectory “jupyter/REHEATFUNQ/”. The location specified in the map below covers the region to which REHEAFUNQ is applied in this analysis. REHEATFUNQ is a Python package that uses a compiled Cython/C++ backend. Compiling REHEATFUNQ requires the Meson build system and a number of scientific libraries and Python packages (and their dependencies) that are listed in the documentation. A Docker image “reheatfunq” is provided as an alternative means of installation. The Docker image comes in two flavors, specified in “Dockerfile” and “Dockerfile-stable”. The former is based on the current “python:slim” image and downloads further dependencies through the Debian package manager, leading to a short image generation time. The latter bootstraps the REHEATFUNQ dependencies from source, aiming to create a reproducible model. To do so, “Dockerfile-stable” depends on the sources contained in “vendor-1.3.3.tar.xz”. If you plan to build the stable image, download both “REHEATFUNQ-1.3.3.tar.gz” and “vendor-1.3.3.tar.xz”, and see the README contained in the latter. Later versions of the “REHEATFUNQ” archive are compatible with the latest “vendor” archive. A quickstart introduction and the API documentation can be found in the linked documentation.

Description: This dataset includes 1-hour GNSS coordinate product processed by GFZ. The observations are from the two GNSS station installed by BKG on the small offshore island of Heligoland in the North Sea. These products are hourly position time series (North, East and Vertical). The 30-second daily RINEX files since 2020 are downloaded from BKG. Together with 5 IGS stations in Europe, the collected RINEX data are processed with the Earth Parameter and Orbit System (EPOS) software from GFZ. The EPOS software uses un-difference carrier phase and pseudo-range observables from GPS and GLONASS L1 and L2 frequencies. They formed an ionosphere-free linear combination to remove the first-order ionosphere effect in the observation. The phase center variation (PCV/PCO) of the satellite and ground station antenna are corrected by IGb14. The station deformation caused by ocean tide loading is modeled by the FES2004 model. Apriori zenith hydro-static/non-hydro-static delay is obtained using the Global Pressure and Temperature model (GPT2) and Vienna mapping functions (VMF) in a 6-hour grid file database. To ensure consistency in the GNSS data analysis, we took the GNSS precise satellite orbits as well as clock products from the 2nd reprocessed (before 2014) and routine (since 2015) yield by EPOS software. The same station parameters are set up as used for the GNSS orbit and clock estimation. All the GNSS data were processed in units of 24 hours periods. The estimated parameters are (i) the receiver clock error for every epoch as white noise, (ii) the hourly station coordinates, (iii) daily tropospheric gradients, (iv) the daily inter-system clock bias for GLONASS, and (v) 2-hour tropospheric wet zenith delays with random-walk constrain.

Description: The Community Stress Drop Validation Study has been organized as a technical activity group (TAG) of SCEC (Southern California Earthquake Center) with the aim of investigating the source parameters of the 2019 Ridgecrest seismic sequence in California. Information about the stress drop TAG are available trough the benchmark web-page (https://www.scec.org/research/stress-drop-validation). Several groups applied different techniques to a shared data set with the objective of extracting source parameters (e.g. seismic moment and corner frequency) and in turn to estimate the stress drop. We applied a spectral decomposition approach known as generalized inversion technique (GIT) and the overall analyses are presented in a series of two articles (Bindi et al 2023a; Bindi et al 2023b). Results in the form of files, figures, and tables are disseminated through this archive.

Description: The network consists of a vertical borehole array equipped with 3C sensors (geophones) for the analysis of swarm earthquakes in the Western Bohemia / Vogtland area located in the German/Czech border region. A surface array is completing the 3D observation of the wave field with 3C sensors (geophones). Waveform data is available from the GEOFON data centre, under network code 6A, and is embargoed until FEB 2035.

Description: An annually resolved chronologies of carbon isotopes from five living oak (Quercus robur) trees have been measured from tree ring cellulose covering up to the last 180 years (1836CE – 2020CE). This tree-ring stable isotope data set was established within the ‘Terrestrial Environmental Observatories’ (TERENO) of the Helmholtz Association. The site “Lake Tiefer See” is subject to the TERENO monitoring activities at the Northeast German Lowland Observatory coordinated by the GFZ German Research Centre for Geosciences in Potsdam. The data set comprises the δ13C records with respect to the international VPDB standard. Lake Tiefer See (53°350 N, 12°320 E) is located 90 km NNW of Berlin in the morainic terrain of the NE-German Polish Basin. It is part of in the N–S trending Klocksin Lake Chain. The sampled trees are growing at the southern shore of the lake. Fifteen co-dominant Quercus robur tree individuals were cored at about 1.3m above ground from two opposite positions using an increment corer of 5 mm diameter (Suunto, Finland or Mora, Sweden).

Description: An annually resolved chronologies of oxygen isotopes from five living oak (Quercus robur) trees have been measured from tree ring cellulose covering up to the last 180 years (1836CE – 2020CE). This tree-ring stable isotope data set was established within the ‘Terrestrial Environmental Observatories’ (TERENO) of the Helmholtz Association. The site “Lake Tiefer See” is subject to the TERENO monitoring activities at the Northeast German Lowland Observatory coordinated by the GFZ German Research Centre for Geosciences in Potsdam. The data set comprises the δ18O records with respect to the international VSMOW standard. Lake Tiefer See (53°350 N, 12°320 E) is located 90 km NNW of Berlin in the morainic terrain of the NE-German Polish Basin. It is part of in the N–S trending Klocksin Lake Chain. The sampled trees are growing at the southern shore of the lake. Fifteen co-dominant Quercus robur tree individuals were cored at about 1.3m above ground from two opposite positions using an increment corer of 5 mm diameter (Suunto, Finland or Mora, Sweden).

Description: Discordant iron-rich ultramafic pegmatites (IRUPs) intersect the UG2 chromitite at many places in the Bushveld Complex but the effects on ore grade and mineralogy have rarely been studied in detail. We investigated a drillcore profile through the UG2 layer affected by IRUP intrusions at Tha-ba mine, NW Bushveld. The work involved quantitative assessment of the ore mineralogy from mineral liberation analysis (MLA), chemical analyses of the main silicate and oxide minerals by elec-tron microprobe (EPMA) and micro-X-ray fluorescence (micro-XRF) element maps of the UG2-IRUP contacts. The data reveal features of UG2-IRUP interaction at different scales. The micro-XRF study of the contacts shows that a thin layer of Fe-Ti-Cr spinel and ilmenite formed on the IRUP side, while chromite on the UG2 side underwent grain coarsening with loss of interstitial space, and developed chemical gradients in Cr, Al, Fe and Ti. On a larger scale, the MLA data document changes in ore and gangue mineral assemblages through most of the meter-thick UG2 layer. This includes formation of secondary hydrous silicates, replacement of PGE- sulfides by PGE alloys and PGE- As-Sb-Bi-Te-Pb phases, and formation of secondary Ni-Cu-Fe sulfides after pentlandite and chalcopyrite.

Description: SAM ("Simplified Analytical Model") is a MatLab-based software that allows for fast and flexible simulations of three-dimensional dyke pathways in an elastic medium. The model was first introduced in "Mechanical modeling of pre-eruptive magma propagation scenarios at calderas" (Mantiloni, L. et al. 2023). In SAM, dykes are modelled as penny-shaped cracks of fixed radius, opening against the local direction of the least-compressive principal stress. The direction of propagation is determined by the gradient of the external stress normal to the crack's plane and the buoyancy force of the magma filling the dyke, calculated at a set of observation points along the crack's tipline. The model can also include a uniform internal pressure within the dyke and compute the stress intensity factor along the crack's tipline, comparing it to the fracture toughness of the host rock to determine if the dyke will advance. SAM needs a model for the stress field of the host rock as input, as well as magma and rock densities, rock elastic properties, the dyke's radius and the number of observation points. The model may be applied to simulate dyke pathways in realistic volcanic settings with different stress sources, and can perform large numbers of simulations in little time. The model does not, however, account for any viscous flow of magma within the dyke, nor the velocity of dyke propagation. Dykes cannot change shape or area during the propagation, and are always bound to be oriented normally to the local least-compressive principal stress axis. This repository also includes data and parameters of the synthetic scenarios discussed in "Mechanical modeling of pre-eruptive magma propagation scenarios at calderas".

Description: As the negative impacts of hydrological extremes increase in large parts of the world, a better understanding of the drivers of change in risk and impacts is essential for effective flood and drought risk management and climate adaptation. However, there is a lack of comprehensive, empirical data about the processes, interactions and feedbacks in complex human-water systems leading to flood and drought impacts. To fill this gap, we present an IAHS Panta Rhei benchmark dataset containing socio-hydrological data of paired events, i.e. two floods or two droughts that occurred in the same area (Kreibich et al. 2017, 2019). The contained 45 paired events occurred in 42 different study areas (in three study areas we have data on two paired events), which cover different socioeconomic and hydroclimatic contexts across all continents. The dataset is unique in covering floods and droughts, in the number of cases assessed and in the amount of qualitative and quantitative socio-hydrological data contained. References to the data sources are provided in 2023-001_Kreibich-et-al_Key_data_table.xlsx where possible. Based on templates, we collected detailed, review-style reports describing the event characteristics and processes in the case study areas, as well as various semi-quantitative data, categorised into management, hazard, exposure, vulnerability and impacts. Sources of the data were classified as follows: scientific study (peer-reviewed paper and PhD thesis), report (by governments, administrations, NGOs, research organisations, projects), own analysis by authors, based on a database (e.g. official statistics, monitoring data such as weather, discharge data, etc.), newspaper article, and expert judgement. The campaign to collect the information and data on paired events started at the EGU General Assembly in April 2019 in Vienna and was continued with talks promoting the paired event data collection at various conferences. Communication with the Panta Rhei community and other flood and drought experts identified through snowballing techniques was important. Thus, data on paired events were provided by professionals with excellent local knowledge of the events and risk management practices.

Description: The Illgraben is a 10 km² steep side valley located in Switzerland. This active debris flow catchment supplies 5-15% of the total sediment load of the Rhône River upstream of Lake Geneva. The 30-80° steep catchment slopes host frequent rock falls and slides. From 2012 to 2014, a network of up to ten Nanometrics Trillium Compact 120s broadband seismometers, sampled by Digos DataCube³ext loggers at 200 Hz (and later by centaur), was deployed in and around the catchment to monitor distributed geomorphic activity.

Description: 24 European annually resolved stable isotope chronologies have been constructed from tree ring cellulose for the last 400 years (1600CE – 2003CE) for carbon and oxygen and for the last 100 years for hydrogen. Data was produced within the ISONET project (400 Years of Annual Reconstructions of European Climate Variability Using a Highly Resolved Isotopic Network,) to initiate an extensive spatiotemporal tree-ring stable isotope network across Europe funded as part of the fifth EC Framework Programme “Energy, Environment and Sustainable Development”. This data set comprises the ISONET δ18O records.

Description: The objective of the Lomonosov satellite development concerns the studies of the ultra-high energy cosmic rays and fast processes in optic, X-ray and gamma-ranges which occur in the upper layers of the Earth's atmosphere and in the Universe. This dataset consists of calibrated electron particle detector (EPD) data of the ELFIN-L instrument which has been described in Shprits et al. (2018). The data rate is 2 measurements per second on 8 physical electron detectors with 12 sub-channels from 21 keV to 4.7 MeV. The data is available from May to November 2016. Due to commissioning-phase and changes in the detector configuration only data from August to November 2016 is usable. Some of the electron detector channels do not provide valid or other than noise measurements. The valid channels are for 21 keV, 30 keV, 44 keV, 1.006 MeV and 1.600 MeV.

Description: 24 European annually resolved stable isotope chronologies have been constructed from tree ring cellulose for the last 400 years (1600CE – 2003CE) for carbon and oxygen and for the last 100 years for hydrogen. Data was produced within the ISONET project (400 Years of Annual Reconstructions of European Climate Variability Using a Highly Resolved Isotopic Network,) to initiate an extensive spatiotemporal tree-ring stable isotope network across Europe funded as part of the fifth EC Framework Programme “Energy, Environment and Sustainable Development”. This data set comprises the ISONET δ13C records.

Description: The here referenced dataset provides eventbased Distributed Acoustic Sensing (DAS) recordings made with an approximately 22 km long dark telecommunication fiber lying in urban Potsdam and surroundings. For each of 164 Mgt;=5 earthquakes occurring in February 2023 and listed by the USGS, one hour of data is provided starting with the event's origin time. Additionally, the whole day of February 14 is provided in hourly files. The data was recorded in the frame of the global DAS month, an initiative to collaboratively record and share simultaneously recorded DAS data from all over the world (https://www.norsar.no/in-focus/global-das-monitoring-month-february-2023). DAS is an emerging technology increasingly used by seismologists to convert kilometer long optical fibers into seismic sensors.

Description: This data set is digital image correlation data, including surface displacement and strain data from laboratory subduction megathrust earthquake cycles. The data consists of grids of surface strain (elastic and permanent), trench-normal surface displacement, vorticity and divergence maps over analog seismic cycles, and time series of surface displacement. The data have been derived using a stereo camera setup and processed with LaVision Davis 10 software. Detailed descriptions of the experiments and results regarding the surface pattern of the strain can be found in Kosari et al. (2023), to which this data set is supplementary. We use three configurations to mimic the along-strike heterogeneous spatiotemporal distribution of frictional locking (Rosenau et al., 2019; Kosari et al., 2022b). A central patch separates two stick-slip zones as an aseismic barrier in all configurations. The frictional properties of the central patch vary as a velocity-strengthening (VS configuration), a velocity-neutral (VN), and a velocity-weakening (VW configuration). The VW zone generates smaller slip events with a higher frequency (i.e., recurrence interval) than the stick-slip zones. Four frictionally different materials have been emplaced on the interface: The sticky-rice as velocity-weakening material (a-blt;0) resulting in stick-slip cycles simulating earthquake cycles, fine-grained sugar and rubber-sand mixture as velocity-strengthening (a-bgt;0) and velocity-neutral (a-b=0) material, and fine-grained salt as velocity-weakening material (a-blt;0) (Kosari et al., 2023).

Description: This data set includes videos depicting the surface evolution (time-lapse photography, digital image correlation [DIC] analysis, and topography analysis), and internal evolution (X-ray CT-imagery and DIC analysis) of four laboratory experiments (analogue models) simulating lithospheric-scale rifting. All experiments were performed at the Tectonic Modelling Laboratory of the University of Bern (UB). Detailed descriptions of the model set-up and results, as well as the monitoring techniques can be found in Zwaan amp; Schreurs (2023a and b).

Description: The present dataset is a comprehensive earthquake catalogue for the Northern Chile subduction zone forearc covering the period 2007-2021, determined from IPOC seismic station data (GFZ and CNRS-INSU 2006; https://doi.org/10.14470/pk615318) plus some auxiliary stations (IPOC = Integrated Plate Boundary Observatory Chile; http://www.ipoc-network.org). The method of automatized earthquake catalogue retrieval, the different relocation steps as well as the different earthquake class labels, and the structures outlined by the seismicity are described in detail in Sippl et al. (2023). The catalogue builds on the one from Sippl et al. (2018; https://doi.org/10.5880/GFZ.4.1.2018.001), but uses a slightly deviating parameter set and a new event category. The columns of the data files are: year, month, day, hour, minute, second, latitude [dec. degrees], longitude [dec. degrees], depth [km], magnitude [ML], identifier The identifier term provides a first-order spatial classification of the seismicity, an explanation is given in Sippl et al. (2023).

Description: The Global Gravity-based Groundwater Product (G3P) provides groundwater storage anomalies (GWSA) from a cross-cutting combination of GRACE/GRACE-FO-based terrestrial water storage (TWS) and storage compartments of the water cycle (WSCs) that are part of the Copernicus portfolio. The data set comprises gridded anomalies of groundwater, TWS, and the WSCs glacier, snow, soil moisture and surface water bodies plus layers containing uncertainty information for the individual data products. All WSCs are spatially filtered with a Gaussian filter to be compatible with TWS. Spatial coverage is global, except Greenland and Antarctica, with 0.5-degree resolution. Temporal coverage is from April 2002 to December 2020 with monthly temporal resolution. Gridded data sets are available as NetCDF files containing variables for the parameter value as anomaly in mm equivalent water height and the parameter’s uncertainty as mm equivalent water height. The latest version of the data is visualized at the GravIS portal: http://gravis.gfz-potsdam.de/gws. From GravIS, the data is also available as area averages for several large river basins and aquifers, as well as for climatically similar regions. G3P was funded by the EU Horizon 2020 programme in response to the call LC-SPACE-04-EO-2019-2020 “Copernicus evolution – Research activities in support of cross-cutting applications between Copernicus services” under grant agreement No. 870353.

Description: GRACE monthly gravity field solutions starting from April 2002 to June 2017 up to degree and order 90 computed with the Celestial Mechanics Approach at AIUB. The time series is an updated of AIUB-RL02 GRACE monthly gravity field time series using Level-1B GRACE data and updated background models. The dataset is created within the framework of the G3P - Global Gravity-based Groundwater Product project (https://www.g3p.eu/), this project has received funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement No 870353.

Description: This dataset contains measurements of viscous and viscoelastic materials that are used for analogue modelling. Proper density and viscosity scaling of ductile layers in the crust and lithosphere, requires materials like Polydimethylsiloxane (PDMS), to be mixed with fillers and low viscoity silicone oils. Changing the filler content and filler material, the density, viscosity and power-law coefficient can be tuned according to the requirements. All materials contain a large amount of PDMS and all but one a small amount of an additional silicone oil. Adding plasticine or barium sulfate lead to shear thinning rheologies with power-law exponents of p〈0.95. Adding corundum powder only has a minor effect on the power-law exponent. Some mixtures also have an apparent yield point but all are in the liquid state in the tested range. In general, the rheologies of the materials are very complex and in some cases strongly temperature dependent. However, in the narrow range of relevant strain rates, the behaviour is well defined by a power-law relation and thus found suitable for simulating ductile layers in crust and lithosphere.

Description: “This ocean-bottom seismometer deployment is part of the SEAMSTRESS project examining tectonic stress effects on Arctic methane seepage. The project is led by PI Andreia Plaza-Faverola at the Centre for Arctic Gas Hydrates, University of Tromsö, Norway. A total of 10 ocean bottom seismometers (OBS) were deployed on Vestnesa Ridge, a sediment drift body just north Knipovich Ridge at its intersection with the Molloy Transform fault (cruise CAGE-20-5). The aim of the experiment was to look for stress release along faults that control seepage sites on Vestnesa Ridge. The network consisted of 8 Lobster type broadband OBS from the German Instrument Pool for Amphibian Seismology (DEPAS) and 2 3C geophones provided by the University of Tromsö. Instruments were free-fall deployed and spaced by about 10 km. They recorded continuously at 100 Hz for 11 months between August 2020 and July 2021.Short, intersecting refraction profiles were shot across all OBS stations, such that OBS positions at the seafloor could be determined within 10 m (cruise CAGE-21-3). Clock drift in this experiment was nonlinear and skew values were only obtained for 6 of the stations. Skew-corrected station VSN01 served as reference station to obtain the clock drift of all other stations using noise cross-correlation and subsequently correct also for the thus determined nonlinearity of time drift. Waveform data are available from the GEOFON data centre, under network code Y9 and are embargoed until July 2025.

Description: The USGS has suspended the distribution of the widely used whole rock reference materials BHVO, BCR and BIR. The goal of this work is to identify a material as similar as possible to the original BIR Islandic basalt. This material can then undergo an ISO-compliant certification of the whole rock powder major and trace element contents. The sampling quarry east of Reykjavik has multiple basalt flows and it is not known which one was originally sampled in 1980 for production of above mentioned reference materials. In this study, three samples were tested to see which is most similar to what was published by Flanigan (1984). Here, the results of this exploratory sample collection are presented, but note that these data are not part of the certification process or represent certified results.

Description: The objective of the Lomonosov satellite development concerns the studies of the ultra-high energy cosmic rays and fast processes in optic, X-ray and gamma-ranges which occur in the upper layers of the Earth's atmosphere and in the Universe. This dataset consists of electron particle detector (EPD) counts data of the ELFIN-L instrument which has been described in Shprits et al. (2018). The data rate is 2 measurements per second on 8 physical electron detectors with 12 sub-channels from 21 keV to 4.7 MeV. The data is available from May to November 2016. Due to commissioning-phase and changes in the detector configuration only data from August to November 2016 is usable. Some of the electron detector channels do not provide valid or other than noise measurements. The valid channels are for 21 keV, 30 keV, 44 keV, 1.006 MeV and 1.600 MeV.

Description: This data set includes overviews and videos depicting the surface evolution (time-lapse photographs, topography data and digital image correlation [DIC] analysis) of 6 analogue models simulating rotational rift tectonics. In these experiments we examined the links between rotational rifting and different distributions of lithospheric weaknesses, and the evolution of the East African Rift System. All experiments were performed at the Tectonic Modelling Laboratory of the University of Bern (UB). Detailed descriptions of the model set-up and results, as well as the monitoring techniques can be found in Zwaan & Schreurs (2023).

Description: The Early-Warning and Rapid Impact Assessment with real-time GNSS in the Mediterranean (EWRICA) is a federal Ministry of Education and Research funded project (funding period: 2020-2023) that aims to develop fast kinematic and point source inversion and modeling tools combining GNSS-based near field data with traditional broadband ground velocity and accelerometer data. Fast and robust estimates of seismic source parameters are essential for reliable hazard estimates, e.g. in the frame of tsunami early warning. Hence, EWRICA aims for the development and testing of new real time seismic source inversion techniques based on local surface displacements. The resulting methods shall be applied for tsunami early warning purposes in the Mediterranean area. In this framework, this repository is a suite of four packages that can be used and combined in different ways and are ewricacore, ewricasiria, ewricagm and ewricawebapp. These four packages can be deployed in a docker container (see instructions below) to demonstrate a possible output of Early-Warning and Rapid Impact Assessment. In the Docker, a probabilistic earthquake source inversion report (ewricasiria) and a Neural network based Shake map (ewricagm) are generated for two past earthquakes whose data (event and waveform) is continuously served by GEOFON servers at regualr intervals to produce and test a real case scenario. The whole workflow is managed by ewricacore, a central unit of work that first fetches the waveform data via the seedlink protocol and event data via event bus or FDSN web service, then collects and cuts waveforms segments according to a custom configuration, and eventually triggers custom processing (ewricasiria and ewricagm in the docker, but any processing can be implemented) whenever configurable conditions are met. The final package, ewricawebapp is a web-based graphical user interface that can be opened in your local browser or deployed on your web server in order to visualize and check all output produced by the docker workflow in form of HTML pges, images and data in various formats (e.g., JSON, log text files).

Description: We present a Python application to download events and data from FDSN webservices (https://www.fdsn.org/webservices/) and compute the events energy Magnitude (Me), producing outputs in several formats (QuakeML, HDF, CSV, HTML). This software has been used to compile a seismic catalogue including Me estimated form P-waves recorded at teleseismic distances in the range 20° ≤ ∆ ≤ 98°, available at GFZ Data Services (Bindi et al., 2023; https://doi.org/10.5880/GFZ.2.6.2023.010). The software complete pipeline (download and energy magnitude computation) can be deployed locally via terminal commands or chained and scheduled on a server to compute the energy magnitude in semi-realtime (e.g. daily or weekly).

Description: The Southern Permian Basin in Central Europe (in Germany and Poland) hosts several sediment-hosted Cu deposits (see Borg et al., 2012). The Cu- and Zn-Pb sulfide mineralization is preserved in the coarse-grained continental siliciclastics of the uppermost Rotliegend (S1), organic matter- and carbonate-rich marine mudstones of the Kupferschiefer (T1) and dolomitic Zechstein Limestone (Ca1). In these datasets, we provide quantitative mineralogical and geochemical data of drill core samples from the Saale Basin in East Germany. The samples include the uppermost Rotliegend sandstone (S1), Kupferschiefer (T1) and lowermost Zechstein Limestone (Ca1), referred as the Kupferschiefer system, from three drill cores (Sangerhausen, Allstedt and Wallendorf). This data publication includes quantitative mineralogy (X-ray diffraction), bulk rock major, minor and trace element geochemistry (X-ray fluorescence and inductively coupled mass spectrometry) and total organic carbon (elemental analyzer).

Description: In support of the Environmental Mapping & Analysis Program (EnMAP) mission [1], the acquisition of accurate and comparable spectroradiometric in-situ measurements is crucial for vicarious validation of the official EnMAP data products [2]. This document provides a guide on properly conducting spectroradiometric field measurements within the scope of EnMAP. It is a summary, of the detailed technical handbook developed by the Commonwealth Scientific and Industrial Research Organisation (CSIRO) [3], the approach established by the Remote Sensing Laboratories (RSL, University of Zurich) [4], on the bases of „Progress in field spectroscopy“ [5], “Field and airborne spectroscopy cross validation - Some considerations” [6] and the experience gained throughout numerous validation efforts for air- and spaceborne sensors by the Remote Sensing and Geoinformatics section at the GFZ Potsdam that have been specially adapted for EnMAP purposes. The following procedure should be used when conducting in-situ measurements of terrestrial surfaces to obtain consistent measurements by applying a repeatable approach throughout the validation phase of the EnMAP mission.

Description: This data publication represents the main outcomes of WP1.200 of Individual Project IP1 and Deliverable D1.1 of the research unit NEROGRAV. The goal of WP1.200 was the realistic representation of modern ocean tide model uncertainties in the form of empirical Variance-Covariance Matrices (VCMs) for the utilization in satellite gravimetric dealiasing. In the following, we describe the data set generation and format. A more detailed description of the processing strategy of the data set can be found in Abrykosov et al. (2021).

Description: The three datasets presented here are high-resolution catalogs containing origin time of seismic events for the same region and time range that have derived using AI-based techniques and a matched filter search. The corresponding standard catalogs from the agencies AFAD and KOERI are available under https://tdvms.afad.gov.tr/ (last accessed 28/07/2022) and http://www.koeri.boun.edu.tr/sismo/2/earthquake-catalog/ (last accessed 28/07/2022), respectively, when searching in the bulletin for longitude 28.80-29.10, latitude 40.4-40.625, and from November 1st 2018 to January 31th, 2019. Specifications for the three catalogs are. (i) Catalog derived utilizing AI-based techniques. We applied the PhaseNet deep learning method (Zhu & Beroza, 2019) to detect and pick the P-and S- waves of seismic events embedded in continuous seismic recordings from 16 stations surrounding the region of interest resampled at 100 Hz. The method was trained on a dataset from Northern California, but has been shown to generalize well to other tectonic settings. The picks were associated into seismic events using the GaMMA association method (Zhu et al., 2022). Manual check of the waveforms from all detections led to 516 seismic events with clear waveforms retained for further processing. (ii) Template matching catalog A. We applied the matched filter algorithm EQcorrscan (Chamberlain et al., 2017) to the two nearby seismic stations with the largest data recovery during the period of interest, ARMT and MDNY. We utilized 14 manually picked template events with M 〉 2 that occurred in the region of interest during the analyzed time period, which were recorded in both stations. As a first criteria to remove false detections, we retained only detections exhibiting a Median Absolute Deviation (MAD) larger than eight. We required detections from different templates to be at least 1.5 seconds apart. To remove duplicate detections (e.g., detections of the same event by different templates), we retained the detections with the highest average correlation if multiple detections occurred within 2.5 seconds. As a second criteria, we calculated cross-correlation derived phase-picks. A pick was declared if the maximum normalized correlation between the signal of the template event and of the detection exceeds 0.7. We correlated the signals in a short window of ±0.3 seconds around the assumed pick time based on a time-shifted version of the template phase-pick. We retained the S-pick exhibiting the higher cross-correlation value with respect to the template. Following this step, we considered only detections with ≥ 2 picks. In case of events with only two picks we ensured that that were from the same station to have control on the ts-tp and therefore the distance of the event from the detecting station. This catalog contains 2,462 seismic events (all manually reviewed) with magnitudes MW in the range [-2.4, 4.5]. Since we were not able to locate the events from this catalog, we considered as “origin time” the time of the first arrival. (iii) Template matching catalog B. We derived a second template matching catalog utilizing twelve of the closest seismic stations displaying high seismic data recovery during the analyzed time period. An initial list of detections was generated following the same steps as for the Template Matching Catalog A, with the additional requirement that all detections must contain at least one picks from one of the two closest stations, ARMT and MDNY. All detections from this catalog were also manually reviewed. The full description of the data processing and creation of the catalog is provided in the article “Stress changes can trigger earthquake sequences in a hydrothermal region south of Istanbul” by Martínez-Garzón et al., currently under review in Geophysical Research Letters.

Description: The GFZ-Landsvirkjun Theistareykir Fibre array is located in the Theytareykir geothermal area, in North Iceland. It is collocated with arrays of broadband seismometers and gravity meters (see e.g., https://doi.org/10.1186/s40517-021-00208-w). The geometry of the fibre array is following the telecom network in the area, and was chosen to test the seismological capabilities of telecom cables in this geothermal environment. We connected an iDAS V2 interrogator from Silixa. The interrogator location is lat=65.898041, lon=-16.966274. The array starts N-S and after 1.5 km, turns towards the East, up to a local transmission antenna station for mobile phones. The length of the path is ~5 km. The length of the cable is actually more than 15 km, as other fibre instance is connected at the transmission antenna station.. Jumps were performed along the cable to geo-locate the channels. The exact location of the fibre can unfortunately not be disclosed. Original recordings at 1000 Hz were downsampled to 200 Hz using a software from INGV-OE (michele.prestifilippo@ingv.it) and are provided in an h5 format. We provide here the first fibre instance (5 km long). The data contain 1 h long recording intervals framing M>5 teleseismic earthquakes recorded in the frame of the global DAS month, an initiative to collaboratively record and share simultaneously recorded DAS data from all over the world (https://www.norsar.no/in-focus/global-das-monitoring-month-february-2023). DAS is an emerging technology increasingly used by seismologists to convert kilometer long optical fibers into seismic sensors.

Description: The dataset is an extended and updated version of the homogenized regional earthquake catalogue of the Marmara region, north-western Turkey, presented in Bohnhof et al. (2017) and Wollin et al. (2018). It is built on the regional Turkish seismicity catalogues provided by AFAD (Disaster and Emergency Management Presidency of Turkey) and KOERI (Kandilli Observatory and Earthquake Research Institute) and spans the time interval 2006-2020. All events available in these two catalogues in the wider Marmara region were combined and dublicate events removed. A total of 13812 events having at least 6 P- and/or S-picks were located using the NLLoc software (Lomax et al., 2000, 2009) in Octtree mode utilizing automatic picks (see Wollin et al., 2018 for details) for all available waveforms. The magnitude range is between M0.3 and M5.7 with time-variable magnitude of completeness and covers the area 39.70S-41.50S and 26.0E-30.65E. The full description of the data and methods is provided in the data description file.

Description: This dataset provides point-shapefiles and geotiffs, related to the figures presented in (Frick et al., 2022a, 2022b). It covers most of northern Germany, with the boundaries defined by the extent of the North German Basin, which is part of the Central European Basin System. The files contain information on the depth (m.b.s. = meter below surface), thickness, temperature, heat in place and heat storage potential of selected geological units and the formations therein. These data are an addendum to the data presented in (Frick et al., 2022a, 2022b), resolving 5 geological units and 9 formations. The data are presented as regularly spaced point-shapefiles, with a spacing of 1000 m. The data were produced as part of the Helmholtz Climate Initiative (HICAM), which focuses on Net Zero 2050 (mitigation) and Adapting to Extreme Events (adaptation). As part of this initiative, estimates of the heat in place and heat storage potential of the subsurface play an important part for mitigation of fossil fuel bound emissions as they pose a promising alternative (geothermal energy). The data presented here, therefore give an overview of areas which might be suited for geothermal applications in the different geothermal target units and formations. We integrated the recently published TUNB Model (BGR et al., 2021) as well as available borehole data, data from the Sandsteinfazies and GeoPoNDD projects (Franz et al., 2018, 2015) and temperature data from two models (Agemar et al., 2014; Frick et al., 2021) the process of which will be described in the following.

Description: IGMAS+ is a software combining 3-D forward and inverse modeling, interactive visualization and interdisciplinary interpretation of potential fields and their applications under geophysical and geological data constrains. The software has a long history starting 1988 and has seen continuous improvement since then with input by many contributors. Since 2019, IGMAS+ is maintained and developed at The Helmholtz Centre Potsdam - GFZ German Research Centre for Geosciences by the staff of Section 4.5 – Basin Modelling and Section 5.2 – eScience Centre with strong ongoing support by H.-J. Götze and S. Schmidt from CAU Kiel. The official webpage of IGMAS+ is available at https://www.gfz-potsdam.de/igmas. Each major version of IGMAS+ is assigned with a DOI. Intermediate releases including changelog can be found at https://git.gfz-potsdam.de/igmas/igmas-releases/-/releases/. This is a collection DOI referring to all versions of IGMAS+. Links to each published version are redundantly available via the "Files" section and the Related Work section ("includes").

Description: The GFZ Potsdam HART (Hazard and Risk Team) in cooperation with the DFG research training group 2043 NatRiskChange at Potsdam University has enabled the acquisition of Airborne Laser Scanning (ALS) and high-resolution optical data which were acquired between 22 September 2021 and 24 October 2021 by the Milan Geoservice company, Spremberg, Germany. This data acquisition took place in the Eifel regions of North Rhine-Westphalia (NRW) and Rhineland-Palatinate (RLP), which were hit by the 14 July 2021 precipitation event leading to widespread severe inundations, flash floods and caused around 185 victims and massive damage to settlements, river geometry and other geomorphic features. The high-resolution ALS and optical data acquisitions aimed at the documentation and quantification of the extent of flood related changes and destructions as well as their reappraisal before diffusion erases traces. Thus, the generated data are valuable for forensic event analysis and future attempts on flood forecasting and warning in the context of scientific and practical purposes.

Description: The classical way to model the stress state in a rock volume is to estimate displacement boundary conditions that minimize the deviation of the modelled stress state with respect to model-independent stress information such as stress magnitude data. However, these data records are usually subject to significant uncertainties and measurement errors. Hence, it has to be expected that not all stress magnitude data records are representative and can be used in a model. In order to identify unreliable stress data records, the stress state that is based on individual data records is solved and compared with observations at a few discrete locations. While this method works, it is not efficient in that most of the solved model scenarios will be discarded. The solving of the entire model consumes immense amount of computation time for a high-resolution model. Yet, the stress state is required at only a very limited number of locations. For linear geomechanical models it is sufficient to estimate the stress state from three model scenarios with arbitrary, but different displacement boundary conditions. These three results can be used to estimate analytically using a linear regression at discrete points stress states based on user-defined boundary conditions. The tool Fast Automatic Stress Tensor Estimation (FAST Estimation) is a Python function that automatizes this approach. FAST Estimation provides very efficiently the stress states at pre-defined locations for all possible boundary conditions. It does not provide the continuous stress field as provided by a solved geomechanical model. Instead, it is a cost-efficient solution for the rapid assessment of stress states at a limited number of discrete locations based on pre-defined boundary conditions.

Description: We present a new, consistently processed seismicity catalogue for the Eastern and Southern Alps, based on the temporary dense Swath-D monitoring network. The final catalogue includes 6,053 earthquakes for the time period 2017-2019 and has a magnitude of completeness of −1.0ML. The smallest detected and located events have a magnitude of −1.7ML. Aimed at the low to moderate seismicity in the study region, we generated a multi-level, mostly automatic workflow which combines a priori information from local catalogues and waveform-based event detection, subsequent efficient GPU-based event search by template matching, P & S arrival time pick refinement and location in a regional 3-D velocity model. The resulting seismicity distribution generally confirms the previously identified main seismically active domains, but provides increased resolution of the fault activity at depth. In particular, the high number of small events additionally detected by the template search contributes to a more dense catalogue, providing an important basis for future geological and tectonic studies in this complex part of the Alpine orogen.

Description: Here, we present model files and example scripts for the Neural network-based model of Electron density in the Topside ionosphere (NET). The model is based on radio occultation data from Gravity Recovery And Climate Experiment (GRACE), Challenging Minisatellite Payload (CHAMP) and Constellation Observing System for Meteorology Ionosphere and Climate (COSMIC-1) missions from 2001 until 2019. The NET model is based on alpha-Chapman functions with a linear decay of scale height with altitude, and consists of 4 sub-models (2 parameters of the F2-peak and 2 parameters of the linear scale height decay). The model uses geographic and magnetic latitude and longitude, magnetic local time, day of year, altitude, solar flux index P10.7, geomagnetic activity index Kp, storm-time SYM-H index as inputs. An example data frame to run the model is provided, as well as the Jupyter notebook to perform an example run.

Description: The Atmosphere and Ocean non-tidal De-aliasing Level-1B (AOD1B) product is widely used in satellite gravimetry to correct for transient effects of atmosphere-ocean mass variability that would otherwise alias into monthly-mean global gravity fields. The most recent release is based on the global ERA5 reanalysis and ECMWF operational data together with simulations from the general ocean circulation model MPIOM consistently forced with fields of the same atmospheric data-set. As background models are inevitably imperfect, residual errors due to aliasing remain. Accounting for the uncertainties of the background model data has, however, proven to be a useful approach to mitigate the impact of residual aliasing. In light of the changes made in the new release of AOD1B, previous uncertainty assessments are deemed too pessimistic and have been revised in the new time-series of true errors: AOe07. One possible way to include the uncertainty information of background models in gravity field estimation or simulation studies is through the computation and application of a variance-covariance matrix that describes the spatio-temporal error characteristics of the background model. The AOe07 variance-covariance-matrix provides this information through (1) a fully populated matrix up to degree and order 40 as well as (2) a diagonal matrix up to degree and order 180.

Description: Because of the multi-stepped pathways of sediment comprising the foreland fold-thrust belt (FFTB), detrital quartz grains that recycle from the FFTB sources contain cosmogenic radionuclides (CRN), such as 10Be and 26Al, accumulated during previous exposure, resulting in inheritance and, hence, anomalously low erosion rates. This inhibits the straightforward use of 10Be as tracers for modern erosion rates and sediment discharge from the FFTB, prevalent at the external edges of collisional orogens such as the Himalaya. We present a novel approach for quantifying the erosion rates of FFTB by comparing measured and modeled CRN concentrations in fluvial sediments. We apply this approach to the Mohand Range, an emergent fault-related fold in the frontal part of the northwestern Himalaya (see the location map below). The 10Be and 26Al datasets presented here were used to calibrate our model, which we used to quantify the erosion rates in and sediment flux from the Mohand Range. Datasets provided here include a summary of the location and depositional age of 33 fluvial sediments and two sandstone samples collected from the Mohand Range, 10Be analysis results of 23 of these fluvial sediments and two bedrock samples, and 26Al-10Be pair analysis results of the remaining ten fluvial sediment samples (Dataset S1). Moreover, the data include the depositional age map of uplifted older foreland sediments across the western Mohand Range (Dataset 2) and the map of best-fit 10Be concentration inherited from Himalayan paleoerosion (Dataset 3) and sediment burial in the foreland (Dataset 4). We also include a map of the best-fit 10Be concentration produced during modern erosion of the Mohand Range (Dataset 5) and a map of the best-fit uplift/erosion rates across the western Mohand Range (Dataset 6). For more information (e.g., sampling method, analytical procedure, and data processing), please refer to the main article (Mandal et al., 2023).

Description: This data publication is supplementary to a study on the climatic controls on leaf wax hydrogen isotopes, by Gaviria-Lugo et al. (submitted). The dataset contains hydrogen isotope ratios from leaf wax n-alkanes (δ2Hwax) taken from soils, river sediments and marine surface sediments along a climatic gradient from hyperarid to humid in Chile. In addition, for each sampling site the hydrogen isotope ratios from precipitation (δ2Hpre) from the grids produced by the Online Isotopes in Precipitation Calculator (OIPC) (Bowen and Revenaugh, 2003). Furthermore, for each sampling site we report mean annual data of precipitation, actual evapotranspiration, relative humidity, and soil moisture, all derived from TerraClimate (Abatzoglou et al., 2018). Also provide data of mean annual temperature and the annual average of maximum daily temperature derived from WorldClim (Fick and Hijmans, 2017). As a final climatic parameter, we also derived data of aridity index from the Consultative Group of the International Agricultural Research Consortium for Spatial Information (CGIARCSI) (Trabucco and Zomer, 2022). In addition to climatic variables, for each site we include land cover fractions of trees, shrubs, grasses, crops, and barren land. These land cover fractions were obtained from Collection 2 of the Copernicus Global Land Cover layers (Buchhorn et al., 2020) via Google Earth Engine. For further comparison here we provide δ2Hwax compiled from 26 publications (see references) that reported both the n-C29 and n-C31 n-alkanes homologues from soils and lake sediments. For each sampling site of the global compilation, we provide δ2Hpre and the same climatic and land cover parameters as for the Chilean data (i.e., precipitation, actual evapotranspiration, relative humidity, soil moisture, aridity index, temperature, fraction of trees, fraction of grasses, etc.), using the same sources. The data is provided here as one single .xlsx file containing 9 data sheets, but also as 9 individual .csv files, to be accessed using the file format of preference. Additionally, 5 supplementary figures that accompany the publication Gaviria-Lugo et al. (submitted) are provided in one single .pdf file. The samples taken for this study were assigned International Geo Sample Numbers (IGSNs), which are included in the provided tables S4, S5 and S6.

Description: The Dec 22nd 2018 flank collapse and tsunami at Anak Krakatau in Indonesia is a key event in geosciences as little is known about the lead-up processes and deformation changes prior to flank failure. We processed Sentinel-1 satellite radar data in both ascending (orbit 171) and descending (orbit 47) acquisition using multi-temporal InSAR with the Small BAseline (SB) method during the 4 years prior to the collapse. The data shows that the flank was already moving for years prior to collapse, demonstrating that developing instability in volcano can be monitored long before a collapse. The southwest flank movement rates averaged approx. 27 cm/yr, but underwent intermittent accelerations coinciding with distinct intrusion events in Jan/Feb 2017 and in Jun 2018. The data archived here supplements the material detailed in Zorn et al. (202X, https://doi.org/XXXXX).

Description: Porphyry copper deposits provide most of the world’s, half its molybdenum reserves and are resources for Zn, Pb, Au, and Ag. The porphyry mineralization is inferred to form on time scales between 50 and 100kyrs whereby the mineralization forming magma chamber is generally built up by multiple intrusive events. The overall source magmatic system can be active for several millions of years. We used the Complex System Modeling Platform (CSMP++) to simultaneously model sill injection, heat transfer, the release of metal-bearing magmatic fluids, the multi-phase flow of saline hydrothermal fluids, and dynamic permeability variations with a continuum porous medium approach. Our modeling studies the volumetric injection rate and its impact on the growth of the magma chamber and the Cu-ore shell but also investigates the influence of hydrothermal convection and fluid release. The setup of each modeling run is changed slightly, either by changing the influx rate, changing the geometry of the magma chamber, or changing the location of fluid release. CSMP was modified to produce vtk and vtu files every 100 years which were read into the Paraview 4.3.1 software to perform the post-processing (including the calculation of the copper enrichment factor and the pore fluid factor). Paraview was then used to produce the displayed videos.

Description: In 2020 and 2021 the STIMTEC-X hydraulic stimulation experiment was performed at ca.~130 m below surface at the Reiche Zeche underground research laboratory in Freiberg, Saxony/Germany. The project temporally followed the STIMTEC experiment at the same site and aimed at understanding the stress heterogeneity of the anisotropic and metamorphic gneiss rock mass. The STIMTEC-X experiment applied the hydraulic stimulation technique in several boreholes at the mine-scale. Complementary to the stimulations, there were active seismic ultrasonic transmission data acquired before the stimulations. We use a seismic monitoring network consisting of six single-component acoustic emission (AE) sensors (sensitivity 1-60 kHz), six hydrophone-like AE sensors (sensitivity 1-40 kHz) and four to twelve single-component Wilcoxon accelerometers (sensitivity 50 Hz-25 kHz). The AE sensors and remained stationary in sub-horizontal and upwards reaching boreholes, the accelerometers were mostly installed along the tunnel walls with one accelerometer in a shallow borehole in each tunnel, and the hydrophone-like AE sensors were installed in the down-going water filled boreholes, but repositioned for each measurement campaign (Figure 1). This data set of 120 active ultrasonic transmission (UT) measurements is supplementary to Boese et al. (2022, in review), which introduces some of the active measurement campaigns of the STIMTEC-X experiment in detail. The whole data set togetter with the “Ultrasonic transmission measurements from six boreholes from the STIMTEC experiment, Reiche Zeche Mine, Freiberg (Saxony, Germany)” [https://doi.org/10.5880/GFZ.4.2.2021.002] was used to evaluate performance measures such as sensitivity and frequency bandwith, coupling, placement and polarity of the hydrophone-like AE sensor compared to AE sensors. The active seismic data provided here are from seven boreholes (BH01, BH05, BH06, BH10, BH14, BH18, BH19) as shown in Figure 1. There are nine tables provided as metadata of which seven contain the STIMTEC-X sensor coordinates for each measurement campaign, the event information of all the 120 UT measurements and the UT picks. The UT measurements were recorded with a sampling rate of 1 MHz and results from an automatic stack of 1024 UT pulses generated by the ultrasonic transmitter and recorded by the STIMTEC-X sensors. The UT measurements are saved in binary file format (fsf file format). Fsf-files can be processed with FOCI software: https://www.induced.pl/software/foci. Each fsf file contains 32768 samples, which corresponds to 0.032768 seconds. All UT event files were manual inspected and phase arrivals identified. These are stored in the fsf-file header as well as in the table STIMTECX_UT_picks.csv.

Description: The Mallik Anticline is a geologic structure in the Mackenzie Delta in the Canadian Arctic. Tectonics throughout the Cenozoic, with compressional phases in the early Eocene to the late Miocene, formed this large, domed structure that is today an important source of hydrocarbons. Gas hydrates occur in the clastic sedimentary rocks of the Oligocene to Pleistocene Kugmallite, Mackenzie Bay, and Iperk sequences, which were essentially formed by deltaic processes. The presence of hydrocarbon gases within the permafrost zone in the Canadian Arctic has led to extensive exploration and production activities in the region since the mid-1960s, and the investigations by geologists and geophysicists have already been published in numerous scientific articles to date. The associated report (Chabab and Kempka, 2023) describes the implementation of the first field-scale 3D static geologic model of the Mallik site, which was created using data from well logs and 2D seismic reflection profiles. The dataset presented here provides elevation depths and thickness data of the three distinct sequence boundaries Kugmallit-Richards, Mackenzie Bay-Kugmallit and Iperk-Mackenzie Bay as well as fault data from the Mallik site.

Description: For the integration of this dataset, several research articles were collected from the catalog of The Global Heat Flow Data Assessment Project. Specially, this data publication encloses all heat-flow data of onshore India. The resulting updated database contains 617 determinations of heat-flow from 36 publications. The data are presented according to the standards defined by the World Heat Flow Database Project and the International Heat Flow Commission (Fuchs et al., 2023)

Description: At valley glaciers, rockwall erosion supplies debris to glacier surfaces. Once deposited on the ice, rockwall debris is passively entrained and becomes part of the glacial system, e.g., forming medial moraines as downglacier transport continues. Where debris occurs supraglacial, it modifies ice ablation and, thus, changes in rockwall erosion and debris supply rates modify glacial debris cover and mass balance and may affect glacier retreat in response to climate change. Yet, estimates on rockwall erosion rates close to glacier surfaces are few and quantifying spatiotemporal supply patterns is not trivial. This data publication is supplementary to the study on rockwall erosion rates at five Swiss valley glaciers around Pigne d’Arolla, by Wetterauer & Scherler (2023). We temporally and spatially assess rockwall erosion by measuring in situ-produced cosmogenic 10Be concentrations ('[10Be]measured') in medial moraine debris, which we systematically sampled along downglacier-profiles, and by comparing records from various medial moraines, which are supplied by rockwalls differing in exposure and morphology. However, as '[10Be]measured' within supraglacial debris is the sum of '[10Be]rockwall', accumulated during rockwall erosion, and '[10Be]transport', accumulated during post-depositional downglacier transport, medial moraine '[10Be]measured' should be corrected for '[10Be]transport'. If glacier velocities through time are known, '[10Be]transport' can be estimated by downglacier debris trajectory modelling. Providing our 10Be dataset and ~40-year records of glacier surface velocities from four of the five valley glaciers (Glacier du Brenay, Glacier de Cheilon, Glacier de Pièce, Glacier de Tsijiore Nouve) is the main objective of this data publication. The dataset of the fifth glacier (Glacier d’Otemma) has already been published as case study by Wetterauer et al. (2022a,b).

Description: In October 2021, GFZ installed together with INGV Catania, Iraci and ASIR Ltd (Advances Seismic Instrumentation & Research) the very first seismic borehole broadband seismometers at two selected sites at Mt. Etna, Sicily (see Fig. 1). The installation was completed under the EU-funded project ‘SiC nano for PicoGeo’ (http://www.picogeo.eu/). Site one is located next to the Astrophysical Observatory at Serra La Nave (SLN) and site two is located in the city of Mascalucia (MAS). At each site one borehole broadband seismometer was permanently installed (cemented) at approximately 70 m depth. In approx. 1-2m distance, a second ground-level borehole 4.5 Hz Geophone was temporarily installed (sanded) at 1 m depth until July 2022 (see Fig. 2). The ground-level geophones served as a local surface reference sensor to better evaluate the increase of signal quality from surface to depth. Test data were evaluated between October 2021 and July 2022. Sensor settings were adjusted during this time period to obtain the best possible data resolution at both test sites. This data publication compiles a segment of waveform recordings utilized for the assessment of data quality from the two installed broadband borehole seismometers, along with noise plots (Fig. 3-5) illustrating the enhancements in the data quality of frequency ranges compared to surface sensors at Mt. Etna.

Description: The data publication contains all heat-flow data of offshore in the Guaymas basin. The data release contains data generated between 1959 and 2019 and constitutes a substantial update and extension compared to the last compilation provided by Becker & Fisher (1991). The data set comprises new heat-flow determinations published after 1991 as well as data from before 1991, which were not included in the Becker & Fisher (1991). The resulting updated database contains 487 determinations of heat-flow at 464 locations from 17 publications. 95% of the reported heat-flow values are determined from marine probe sensing technique and 5% in boreholes.

Description: Physical samples (or specimen or artefacts) represent the origin of research results in many scientific disciplines. Assigning persistent identifier (PID) to samples is a fundamental step to make them discoverable and traceable in unambiguous way over the Web. The International Generic Sample Number (IGSN) is a PID for physical samples and connecting these with their online description following a dedicated metadata schema. Sample descriptions of samples are available in various formats and detail. In order to publish them in a standardized manner and to automate and standardize the preparation and processing, the software product SAMIRA (Sample IGSN Registration Automation) was created as part of the Project FAIR WISH, funded by the Helmholtz Metadata Collaboration (HMC). SAMIRA aims to automate the generation of Metadata XML-Files for the Registration of PIDs from different input sources (e.g. the FAIR Samples Template, Wiezcorek et al., 2023). This first version of SAMIRA implements the creation of IGSN metadata and Datacite metadata and the respective registration.

Description: This dataset provides friction data from ring-shear tests on corundum sand “NKF120” used in analogue modelling of tectonic processes as a rock analogue for “strong” or “high density” layers in the earth’s upper crust (e.g. Klinkmüller et al., 2016) or as an additive to PDMS silicone oil to increase its density and non-linearity (Zwaan et al., 2018). According to our analysis the material shows a Mohr-Coulomb behaviour characterized by a linear failure envelope. Peak, dynamic and reactivation friction coefficients of corundum sand are µP = 0.75, µD = 0.57, and µR = 0.62, respectively (Table 5). Cohesion of the material ranges between 100-150 Pa. The material shows a minor rate-weakening of ~1% per ten-fold change in shear velocity v and a stick-slip behaviour at low shear velocities. The tested bulk material consists of corundum sand with grain size of 90-120 µm (Table 1). Corundum sand is produced as industrial abrasive materials and sold e.g. by the company Nico Bosse Strahlmittel Berlin. The data presented here are derived by ring shear testing using a SCHULZE RST-01.pc (Schulze, 1994, 2003, 2008) at HelTec, the Laboratory for experimental tectonics at the Helmholtz Center Potsdam – GFZ German Research Centre for Geosciences in Potsdam, Germany.