ALBERT

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: 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: 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 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: 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: This data set contains measurements of an underground hydraulic fracture experiment at Äspö Hard Rock Laboratory in May and June 2015. The experiment tested various injection schemes for rock fracture stimulation and monitored the resulting seismicity. The primary purpose of the experiment is to identify injection schemes that provide rock fracturing while reducing seismicity or at least mitigate larger seismic events. In total, six tests with three different injection schemes were performed in various igneous rock types. Both the injection process and the accompanied seismicity were monitored. For injection monitoring, the water flow and pressure are provided and additional tests for rock permeability. The seismicity was monitored in both triggered and continuous mode during the tests by high-resolution acoustic emission sensors, accelerometers and broadband seismometers. Both waveform data and seismicity catalogs are provided.

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.

Description: This dataset contains predictions of Earth orientation parameters (EOP) submitted during the Second Earth Orientation Parameters Prediction Comparison Campaign (2nd EOP PCC). The 2nd EOP PCC has been carried out by Centrum Badań Kosmicznych Polskiej Akademii Nauk CBK PAN in Warsaw in cooperation with the GFZ German Research Centre for Geosciences in Potsdam (Germany) and under the auspices of the International Earth Rotation and Reference Systems Service (IERS) within the IERS Working Group on the 2nd EOP PCC. The purpose of the campaign was to re-assess the current capabilities of EOP forecasting and to find most reliable prediction approaches. The operational part of the campaign lasted between September 1, 2021 and December 28, 2022. Throughout the duration of the 2nd EOP PCC, registered campaign participants submitted forecasts for all EOP parameters, including dX, dY, dPsi, dEps (components of celestial pole offsets), polar motion, differences between universal time and coordinated universal time, and its time-derivative length-of-day change. These submissions were made to the EOP PCC Office every Wednesday before the 20:00 UTC deadline. The predictions were then evaluated once the geodetic final EOP observations from the forecasted period became available. Each participant could register more than one method, and each registered method was assigned an individual ID, which was used, e.g., for file naming. The dataset contains text files with predicted parameters as submitted by campaign participants and MATLAB file which is a database with all correct predictions from each participant loaded into a structure.

Description: This dataset provides friction data from ring-shear tests on garnet sand 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). According to our analysis the material shows a Mohr-Coulomb behaviour characterized by a linear failure envelope. Peak, dynamic and reactivation friction coefficients of garnet sand are µP = 0.83, µD = 0.61, and µR = 0.73, respectively (Table 5). Cohesion of the material ranges between 20-120 Pa. The material shows a no significant rate-dependency. 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. The RST is specially designed to measure friction coefficients µ and cohesions C in loose granular material accurately at low confining pressures (〈20 kPa) and shear velocities (〈1 mm/sec) similar to sandbox experiments. In this tester, a granular bulk material layer is sheared internally at constant normal stress σN and shear velocity v while shear force and lid displacement (corresponding to density and volume change ΔV) are measured continuously. For more details see Klinkmüller et al. (2016).

Description: In the southern Central Andes (~32°S), subduction of the Nazca oceanic plate beneath the South American continental plate becomes horizontal. The growth of the Altiplano-Puna Plateau is covalently related to the southward migration of the flat subduction, but the role of subduction geometry and the plate strength on current and long-term deformation of the Andes remains poorly explored. This study takes a data-driven approach of integrating the previous structural and thermal model of the lithosphere of the southern central Andes into a 3D geodynamic model to explore the different parameters contributing to the localization of deformation. We simulate visco-plastic deformation using the geodynamic code ASPECT. The repository includes parameter files and input files for the reference model (S1) and the following alternative simulations: a series of models with variation in friction at the subduction interface (S2a-d), a series of models with variation in sedimentary strength (S3a-d), a series that studies the effect of topography (S4), and a series that studies the effect of plate velocities. In addition, a readme file gives all the instructions to run them.

Description: The Northeast Atlantic (NEA) region has long been a subject of interest due to its complex geological history, particularly regarding the interaction between the Iceland plume and the lithospheric plates. In this data publication, we present a comprehensive three-dimensional structural and density model of the NEA crust and uppermost mantle, consolidating and integrating a wide range of previously fragmented data sets. Our model highlights the influence of the Iceland plume on the region's geological evolution, shedding light on the mechanisms that facilitated the continental breakup between Europe and Laurentia during the earliest Eocene period. The whole workflow and methods are described in Gomez Dacal et al. (2023) and its Supplementary Information.

Description: he joint ESA/NASA Mass-change And Geosciences International Constellation (MAGIC) mission 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. The MAGIC mission 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 degrees at approximately 400 km altitude. The ESA--led pair P2 Next Generation Gravity Mission (NGGM) 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. sub--weekly) resolution, shorter latency, and higher accuracy than GRACE and GRACE-FO. This will pave the way to new science applications and operational services. The performances of different MAGIC mission scenarios for different application areas in the field of geosciences were analysed in the frame of the initial ESA Science Support activities for MAGIC. The data sets provided here are the Level-2a simulated gravity field solutions of MAGIC scenarios and the related reference signal that were used for these analyses. The .gfc files in the folders monthly (31-day solutions) and weekly (7-day solutions) contain the estimated (HIS) coefficients (Cnm, Snm) as well as the formal errors (SigCnm, SigSnm) of the different MAGIC scenarios. In order to compute the coefficient errors, the reference/true HIS coefficients contained in the folder HIS_reference_fields need to be subtracted from the estimated HIS coefficients. The data sets provided here comprise the Level-2a simulated gravity field solutions of MAGIC scenarios and the related reference signal (based on Dobslaw et al. 2014; 2015) that were used for the above analyses.

Description: The Kupferschiefer districts in Central Europe contain some of the world’s highest-grade sediment-hosted stratiform Cu (SSC) deposits (see Borg et al., 2012). The high-grade sulfide mineralization in the organic matter-rich marine mudstones of the Kupferschiefer (T1), and also in the underlying continental sandstones of the uppermost Rotliegend (S1) and overlying Zechstein Limestone (Ca1), in the Saale subbasin (Eastern Germany) are dominantly formed as a replacement of calcite cement (Mohammedyasin et al., 2023). We provide carbonate major element chemistry, carbon isotope composition of organic matter, and calcite carbon and oxygen isotope microanalysis datasets of drill core samples from the Saale subbasin in Eastern Germany. The samples include the uppermost Rotliegend sandstone (S1), Kupferschiefer (T1) mudstones and lowermost Zechstein Limestone (Ca1), referred as the Kupferschiefer system, from three drill cores (Sangerhausen, Allstedt and Wallendorf).

Description: This dataset comprises 47 fluid samples from 11 geothermal sites (Germany, Austria, Iceland, Turkey, Netherlands, Belgium, French West Indies). The samples were collected within the REFLECT project (Redefining geothermal properties at extreme conditions to optimize future geothermal energy extraction). The focus with these analyses were on the organic compound composition of the fluids, since they are rarely included in the analyses of fluids taken from geothermal power plants. Understanding the organic compound composition of geothermal fluids might help to better understand chemical reactions within the fluids and might help to mitigate problems that arise with the operation of a geothermal power plant such as mineral precipitation (scaling) and corrosion of the casing and pipes.

Description: We provide present-day glacial isostatic adjustment (GIA) gravity changes simulated with the numerical model VILMA. The effects of Earth and ocean pole tide due to rotational deformation (considered in VILMA) were removed. The dataset contains the solutions for 56 GIA model ensemble members including 54 3D models and 2 1D models. The results are provided as Stokes coefficients with a resolution of degree/order 170.

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/.

Description: The dataset presented here encompasses the results of the geochemical analyses of water and recent carbonate samples collected in the El Peinado basin located in the Southern Puna Plateau in Catamarca, Argentina. This system formed by the hypersaline lake Laguna del Peinado, numerous hydrothermal springs, and the small hypersaline lake Laguna Turquesa, provides a natural laboratory to study carbonate formation and the mechanisms that control the incorporation of various elements and isotopes into their structure under a broad range of geochemical conditions. Geochemical analyses include data on the physicochemical parameters, elemental, and isotopic (δ18O, δ2H, δ11B) composition of the waters, and data on the elemental and isotopic (δ18O, δ13C, δ11B) composition of the carbonates. These data allowed us to calculate element partition coefficients and isotopic fractionation between coupled water-carbonate samples from this natural setting, which are also included here. This dataset also includes the results of water modelling using the software PHREEQC, which contains data on the chemical speciation of carbon and boron, the species contributing to total alkalinity, and mineral saturation indices. This information is useful for all those dealing with geochemistry of hypersaline lakes, geochemistry of continental carbonates, as well as paleoenvironmental and paleoclimatic studies using lake carbonates as archives. These data correspond to the research article “On the origin and processes controlling the elemental and isotopic composition of carbonates in hypersaline Andean lakes”. The full description of the data is provided in the data description file.

Description: This data publication contains the datasets generated in a study aiming at reconstructing paleoclimatic conditions during the late Holocene in northern Philippines. The data come from samples taken from sediment lakes retrieved from Bulusan Lake on the Luzon Island, Philippines. On these samples we measured the stable-hydrogen-isotopic composition of terrestrial-lipid biomarkers to reconstruct ENSO dynamics and past hydrological conditions, pollen data to reconstruct past vegetation, and magnetic susceptibility measurements of the sediment cores to reconstruct past erosion rates. This is complemented with isoGSM2 data to constrain modern hydrological conditions. The data was generated between 2013-04 and 2020-9. The data files are provided in Excel and tab-delimited text versions.

Description: This data publication presents quantitative DNA data obtained through fluorometric detection of genomic DNA and the estimation of 16S rRNA gene copies using quantitative Polymerase Chain Reaction (qPCR). The data encompasses various soil and rock samples collected across a climate gradient. The DNA was extracted using a protocol enabling the separate analysis of intracellular DNA (iDNA) and extracellular DNA (eDNA) from the same sample. The primary objective of this study was to enhance a previously established method developed by Alawi et al. (2014) for analyzing terrestrial samples by introducing modifications to the extraction buffer. Phosphate buffers at two different concentrations (120 mM and 300 mM), EDTA (300 mM), and a high-concentration phosphate buffer in combination with EDTA (300 mM each) were tested in conjunction with a detergent mix (detailed in Medina et al., 2023; submitted). Thorough tests, including spiked DNA experiments and cell counts, were conducted on one low biomass sample to validate the extraction setups. The two most effective extraction protocols were then applied to all samples from the four designated sites and compared with the phosphate buffer described by Alawi et al. (2014), resulting in the calculation of improvement factors. The resulting dataset provides valuable quantitative DNA information and estimates of 16S rRNA gene copies across diverse soil and rock samples along a climate gradient. The modifications made to the extraction buffer demonstrated improved efficiency in extracting especially iDNA compared to the original method. These findings contribute to the refinement and optimization of DNA extraction protocols for terrestrial samples, enabling more accurate and comprehensive analyses of microbial communities in different environments.

Description: The overarching goal of the Drilling Overdeepened Alpine Valleys (DOVE) project will be to date the age and extent of past glaciations. Formerly-glaciated areas are often characterized by deeply incised structures, often filled by Quaternary deposits. These buried troughs and valleys were formed by glacial overdeepening, likely caused by pressurized subglacial meltwater below warm-based glaciers. Results of this drilling campaign, supported by new dating technologies, will further provide critical data on 'how' and 'at which rate' glacial erosion affects such mountain ranges and their foreland. These processes are also of fundamental importance for evaluating the safety of radioactive waste disposal sites, which are planned in areas of former glaciations. Moreover, results of this project will fill gaps in the knowledge of paleoclimate and atmospheric circulation patterns during past glacial epochs and how these patterns affected ice build-up. The operational data sets include the drill core documentation from the mobile Drilling Information System (mDIS), full round core scans, MSCL data sets, a preliminary core description and the geophysical downhole logging data that were acquired during and subsequent to the drilling operations. All downhole logs and core depth were subject to depth correction to a common depth master (cf. operational report for detailed information). The data are described by two scientific reports, the Operational Report (https://doi.org/10.48440/ICDP.5068.001) and the Explanatory Remarks on the Operational Datasets (https://doi.org/10.48440/ICDP.5068.002).

Description: This is an Arctic-delta reduced-complexity model that can reproduce the 2-m ramp feature observed in most Arctic deltas. The model is built by first reconstructing from published descriptions of the DeltaRCM-Arctic model (Lauzon et al., GRL, 2019), which is, in turn, based on DeltaRCM by Liang et al. (Esurf, 2015). All the modifications and refinements leading to this model (ArcDelRCM.jl) are detailed in a manuscript submitted to Earth Surface Dynamics journal for publication (Chan et al., 2022: esurf-2022-25). Options are retained to run this model with the "DeltaRCM-Arctic" (reconstruction) setting. The code is written purely in Julia language.

Description: The Gravity Recovery and Climate Experiment-Follow-On (GRACE-FO) satellite mission, consisting of two satellites, each carry a magnetometer as part of its attitude orbit control system (AOCS). After careful calibration, the data acquired through them can be used for scientific purposes by removing artificial disturbances from other satellite payload systems. This dataset is based on the dataset provided by Michaelis et al. (2021, https://doi.org/10.5880/GFZ.2.3.2021.002) and uses a similar format. The platform magnetometer data has been calibrated against CHAOS-7 magnetic field model predictions for core, crustal and large-scale magnetospheric field (Finlay et al., 2020, https://doi.org/10.1186/s40623-020-01252-9) and is provided in the ‘chaos’ folder. The calibration results using a Machine Learning approach are provided in the ‘calcorr’ folder. Michaelis’ dataset can be used as an extension to this dataset for additional information, as they are connected using the same timestamps to match and relate the same data points. The exact approach based on Machine Learning is described in the referenced publication. Additionally, in the folder ‚fac’, field-aligned current derived from the magnetic field data are provided. There exists a similar dataset with calibrated magnetic data from the GOCE satellite mission under https://doi.org/10.5880/GFZ.2.3.2022.002 (Styp-Rekowski et al., 2022).

Description: The current global dataset of drainage system shapes has a relatively low spatial resolution. We obtained a new dataset (Basin90m) by calculating the drainage basins larger than 50 km2 globally using a 90-meter resolution Digital Elevation Model (DEM). The total number of drainage basins is 667629. For each drainage basin, we extracted the spatial distribution of the longest river channel and the sinuosity of the river. We computed fundamental geometric parameters for the drainage basins, such as area, length, width, aspect ratio, slope, and elevation. Basin90m consists of vector files (ESRI Shapefile format) containing global drainage basins and river channels. The file sizes for the basin and river data are 7.8 and 2.5 GB, respectively. All calculations were automated using a MATLAB script. For a more detailed description of Basin90m, please refer to our submitted data description article titled "A global dataset of the shape of drainage systems". The Basin90m dataset includes data in four sections. The first section comprises drainage basins globally with an area larger than 50 km². The data format is ESRI Shapefile. Eight morphometric indi-ces of the drainage system are stored in the attribute table of the basin shapefile. The "Basins" folder contains six subfolders, each representing a continent. Each continent's subfolder contains all the basins in that continent, categorized by different stream orders. For instance, the "South America" subfolder contains nine shapefile files corresponding to stream orders 1-9. The names of the shapefile files include their continent and stream order information. For example, "South_America_Basin_8.shp" represents all basins in South America with a stream order of 8. The second part of the Basin90m data consists of global main river channels. The longest river channel of each basin is stored in a folder named "Rivers". The internal structure of this folder is the same as the "Basins" folder. For instance, "South_America_River_8.shp" represents the main river channels in South America with a stream order of 8. The third part of Basin90m data is an Excel file named "Basin90m". This file contains eight morphometric parameters for all the basins. It includes both a globally merged sheet and sheets distinguishing different stream orders. The fourth part of Ba-sin90m data is a folder named "Matlab_code", which contains Matlab code for the automated ex-traction of drainage systems and their morphometric parameters.

Description: This data set provides a series of experiments from ring-shear tests (RST) on various materials that are used at several laboratories worldwide. The data contains the results of slide-hold-slide tests and the processed outputs of standardized ring shear tester data from related publications. Additionally, microscopy images of the materials under plain and polarized light are provided. The time dependent restrengthening of the materials is quantified using slide-hold-slide tests. This restrengthening has implications on the reactivation potential of granular shear zones in analogue models. With the provided software we first analyze the experimental data and then compare the angles and stresses needed to reactivate normal faults in the materials. We find that while healing rates are low, the majority of samples can not reactivate normal faults that are generated through extension of an analogue model.

Description: The META-WT project was designed to perform a 4-weeks seismic experiment in Germany with a dense seismic array of ~400 three-component geophones that covered (1) a 2.5km x 2.5km wind farm area in Brandenburg, Germany, with almost 200 wind turbines (WTs) and a well-studied subsurface structure and (2) a 20-km long radial line from the center of the wind farm with one geophone every half-kilometer. The objective was to capture the spatio-temporal seismic wave-field signature of the wind farm from continuous recordings of ambient noise. Due to the dense interstation distance and proposed geometry the experiment allowed for analyzing both small-scale wave field characteristics at an unprecedented spatial resolution and the longer distance radiation pattern of the wind farm. Waveform data is available from the GEOFON data centre, under network code XF, and is embargoed until Jan 2025.

Description: The Gt BTrKoe 1/2021 borehole was drilled in the framework of a research project called GeoFern, funded by the German Ministry for Economic Affairs and Climate Action under the grant number 03EE4007. The overall objective of this research project was to support the development of the geothermal heat utilization for urban areas. Therefore, the integration of reservoir utilization concepts into heat supply systems need to be studied. The GeoFern project aimed to contribute to the knowledge on geological structure and the lithological composition of the subsurface to minimize the explorational risks for future site developments in SE Berlin, Germany. It focused on the exploration of possible Mesozoic aquifers, suitable for aquifer thermal energy storage (ATES) in depths of up to 500 m. As stopping criteria for drilling, the presence of terrestrial (arid) clayey Keuper sediments of the Exter Formation (Upper German Triassic) were defined. In this data publication we provide the results of the investigations and measurements conducted on site in the field laboratory as well as the open-hole geophysical well-logging data of the Gt BTrKoe 1/2021 borehole acquired by a commercial contractor. In addition, a temperature log of the borehole, measured by the GFZ about two months after the end of drilling activities, is part of this data publication. The drilling of the Gt BTrKoe 1/2021 borehole started at the 15th of November 2021 with the setting of the conductor pipe and reached its final depth of 456 m in Triassic sediments on the 19th of December 2021. The drilling was conducted in two main sections using two different technologies. For the upper section, covering Cenozoic sediments and reaching a depth of 211 m, reverse drilling technology was used. This section comprises the Quaternary to Tertiary groundwater system and the Tertiary “Rupelton” (Oligocene, Rupelian). The latter represents an about 100 m thick clayey succession that do act as a regional aquitard, separating the deeper saline groundwater systems from the upper utilized (freshwater) groundwater levels. After setting and cementing of the casing, the borehole was further deepened by using conventional Rotary drilling technology. Due to the lack of knowledge on the geological situation of the pre-Cenozoic strata before the drilling, this section represents the most relevant part for answering the research goals of the project. In order to allow the most accurate description and characterization of the drilled strata, this section was completely cored using wireline coring equipment with 3-m core barrels. In total, 90 core runs were conducted and 197.4 m of cores retrieved, showing a core recovery factor of 81%. The core show a mean core diameter of about 100 mm. The drilling was stopped after encountering the multicolored terrestrial playa sediments of the Upper Triassic in the last core run. While the token cutting samples were not assigned with International Generic Sample Numbers (IGSN), the borehole (Norden, 2022) and all taken cores were registered with IGSNs.

Description: This data set is a description of a novel analogue modelling method used to run lithospheric-scale tectonic models, and to uniquely monitor these models through X-Ray CT-scanning techniques at the Tectonic Modelling Lab of the University of Bern (Switzerland). It includes information on the model set-up and model materials, and includes a step-by-step description of the general modelling procedure. A first application of this novel procedure, for the simulation of lithospheric scale rifting processes can be found in Zwaan & Schreurs (2023a) in Tectonics, with supplementary data publicly available via GFZ Data Services (Zwaan & Schreurs 2023b). The results of this work prove the feasibility of the method, and opens the door to a broad variety of new tectonic modelling studies

Description: This dataset provides friction data from drained ring-shear tests on a wet (water saturated) silica powder-glass beads-PVC powder mixture (40:40:20 wt.%) “CM2”, used in analogue modelling of tectonic and erosion processes as a rock analogue for the earth’s upper crust (e.g. Conrad et al., 2023, Reitano et al., 2020, 2022. 2023). According to our analysis the materials show a Mohr-Coulomb behaviour characterized by a linear failure envelope. Peak, dynamic and reactivation friction coefficients of CM2 are µP = 0.66, µD = 0.58, and µR = 0.61, respectively. Cohesion of the material ranges between 60-230 Pa. The tested bulk material CM2 consists of a mixture of 40 wt. % silica powder, 40 wt.% glass beads and 20 wt.% PVC powder which has been saturated with water (Table 1). Specification of silica powder is “Ventilated Quartz VR16” (https://www.valligranulati.it/products-granules-quartz-marble-sands-premixed/sheet-m/ventilated-quartz) by the company Valli Granulati S.r.l. (Italy). Ventilated quartz is obtained by micronisation of quartz sands with a high content of SiO2 (around 96%), and used e.g. in paints and abrasives. It should be handled with care to omit generation of dust and a half mask (filter class FFA1P2 RD) should be worn because it can harm the human respiratory tract with the potential of causing silicosis. Glass beads used here have a size (diameter) of 700-110 µm and their individual properties are described in detail Pohlenz et al. (2020). The commercial name for the PVC powder is “PVC K.57 Inovyn 257RF” by the company TPV Compound (Italy). PVC powder is mainly used for cleaning industrial structures (as abrasives) or for the production of PVC tubing, plastic sheets etc. The composition of this PVC powder is the same of the common Polyvinyl chloride. According to the regulation CE n.1272/2008 (CLP), this type of PVC powder is classified as not dangerous for the supply, also thanks to its low value of density and round shape.

Description: Under the coordination of the German Research Centre for Geosciences GFZ the first European onshore CO2 storage project was initiated in 2004 at Ketzin, approximately 25 km west of Berlin, Germany. About 67 kt of CO2 (purity 〉 99.9%) were injected into a saline aquifer from June 2008 until August 2013. All project stages were accompanied by a comprehensive monitoring and modelling program, focusing on the investigation of the processes involved and to assure leakage-free CO2 injection and save geological storage. Hence, methods from different geoscientific disciplines were applied, targeting the reservoir itself, the cap rock, the above-zone and the surface (Martens et al. 2015, Wipki et al. 2016). Here we report on the results of the long-term surface monitoring with continuous soil CO2 flux measurements. A profound and extensive database of measurements performed before injection started serves to interpret data during and after CO2 injection (Zimmer et al. 2011). As the CO2 flux measurements reflect the specific site conditions, which can vary locally and over time, trends must be interpreted carefully. After an exploration phase in 2004 and drilling of the first wells in 2007, CO2 was injected between 2008 and 2013 into Upper Triassic sandstones at a depth of 630 to 650 m. This reservoir is overlain by more than 165 m of shaley cap rocks. The site itself is located at the southern flank of the Roskow-Ketzin double-anticlinal structure (Martens et al. 2014, Förster et al. 2006, Förster et al. 2009) and the stored CO2 mainly migrated in northern to western direction (Wipki et al. 2016 ). Additionally, soil profiles of 70 cm depth were studied for their structure and carbon and nitrogen concentrations. The results helped to explain the spatial variations of the soil CO2 fluxes at the different locations (Zimmer et al. 2011). However, as most of the sampling positions are located next to agricultural roads and fields, an influence of used fertilizers and arable farming on the soil structure, chemical composition and the soil biology cannot be ruled out.

Description: The occurrence of exposed high-metamorphic rocks (granulites) in combination with various gravity anomalies aligned along the direction of Variscan strike characterize a special terrain (Saxothuringikum) which has been sandwiched between two major tectonic units during the Variscan orogeny. Near surface geological studies show evidence that the Saxothuringian zone represents extended crust. Therefore the model of a "metamorphic core complex" is often used to explain the exhumation of the "Saxonian granulites". The thickness of the crust, the geometry of the Moho, and the composition of middle and lower crust that underlie such" metamorphic core complexes" have remained largely unconstrained. Because these physical parameters are critical for understanding the extensional processes acting at depth, we have carried out a seismic refraction experiment in order to resolve the deeper structure of an exposed "granulite-complex". From May 6th to May 13th 1995 a seismic refraction - wide angle reflection experiment was carried out as part of the DFG-priority program: "Orogenic processes – their quantification and simulation at the example of the Variscides". Two lines, A and B, were completed in two deployments (see map in GRANU95_report.pdf). In total 12 shots were fired and over 4500 seismograms were collected using 130 instruments. Only two different types of instruments (Reftek and PDAS) have been used for recording the explosions. All instruments were equipped with a 3-component 1Hz seismometer. The 90 km long NW-SE line (deployment A, 74 instruments) from Leipzig to the Erzgebirge through the Saxonian Granulites was carried out on the 8th and 9th May 1995. Additionally 56 stations were placed symmetrically to shotpoint D along line B (perpendicular to line A). Shots were fired on locations A1, A3, A4, A2 (see map in GRANU95_report.pdf). The station spacing for this deployment was around 1.3 km. The 260 km long SW-NE line (deployment B, 93 instruments) from Dresden to Bamberg, also crossing the Saxonian Granulites was completed from 11th to 13th May 1995. Every second instrument from deployment A was moved to complete line B. Shots were fired on locations B, C, D, E, F, G, H and I (see map in GRANU95_report.pdf) and recorded along line B and line A (perpendicular to line B) at a receiver spacing of about 2.6 km.