ALBERT

Description: Spherical harmonic coefficients representing an estimate of Earth's mean gravity field during the specified timespan derived from GRACE-FO mission measurements. These coefficients represent the full magnitude of land hydrology, ice, and solid Earth processes. Further, they represent atmospheric and oceanic processes not captured in the accompanying GAC product.

Description: In May 2018 a volcano-seismic sequence accompanied the upward migration of a magmatic intrusion from Moho depth to the seafloor led to the drainage of the deep magmatic reservoir and to the birth of a submarine volcano offshore the island of Mayotte, Comoro Islands. This process of magma transport was accompanied by an intense seismic swarm and peculiar long-duration very long period signals. Between 1 January 2018 and 1 May 2019 we detected 407 sources of very long period signals and 6990 volcano-tectonic earthquakes. This report collects detection, location and source parameters catalogs for these two sets of earthquake sources. This data publication provides the catalogues of very long period (VLP) signals and volcano-tectonic (VT) earthquakes, as discussed in Cesca et al. (2019). Here, methods and data used to create the different catalogues are only briefly discussed; a more accurate description is given in Cesca et al. (2019), which furthermore discusses the different processes of dike migration, undersea eruption, deep reservoir drainage and overburden sagging which are responsible for the seismic activity.

Description: This dataset provides friction data from ring-shear tests (RST) for a quartz sand used in analogue experiments at the Institute of Geophysics of the Czech Academy of Science (IGCAS) (Kratinová et al., 2006; Zavada et al., 2009; Lehmann et al., 2017; Krýza et al., 2019). It is characterized by means of internal friction coefficients µ and cohesion C. According to our analysis the materials show a Mohr-Coulomb behaviour characterized by a linear failure envelope. Peak friction coefficients µP of the tested material is ~0.75, dynamic friction coeffi-cients µD is ~0.60 and reactivation friction coefficients µR is ~0.64. Cohesions of the material range between 90 and 130 Pa. The material shows a minor rate-weakening of 〈1% per ten-fold change in shear velocity v.

Description: QuantumFrontiers_HLSST_SLR_COMB2019s is a series of monthly gravity field models based on high-low satellite-to-satellite (HLSST) tracking and satellite laser ranging (SLR) data up to degree and order 60. The combination of HLSST and SLR data is done on the normal equation level using Variance Component Estimation. The series spans from 2003 to 2018 and thus covers the entire period between GRACE and GRACE Follow-On. It is therefore a prime candidate to bridge the data gap between these two satellite mission considering long-wavelength features on a global scale. The model has been developed with data contributions from the Astronomical Institute, University Bern (AIUB), the Institute of Geodesy, Theoretical Geodesy and Satellite Geodesy, Graz University of Technology, the Institute for Geodesy, Leibniz University Hannover and the European Space Agency. More details on the processing can be found in "Time-Variable Gravity Signal in Greenland Revealed by High-Low Satellite-to-Satellite Tracking" (Weigelt et al, 2013, https://doi.org/10.1002/jgrb.50283) Funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy – EXC-2123 QuantumFrontiers – 390837967.

Description: GRACE/GRACE-FO Level-3 product based on GFZ RL06 Level-2B products (Dahle & Murböck, 2019) representing Ocean Bottom Pressure (OBP) variations provided at 1° latitude-longitude grids as defined over ocean areas. The OBP grids are provided in NetCDF format divided into yearly batches. The files each contain seven different variables: 1) 'barslv': gravity-based barystatic sea-level pressure 2) 'std_barslv': gravity-based barystatic sea-level pressure uncertainties 3) 'resobp': gravity-based residual ocean circulation pressure resobp 4) 'std_resobp': gravity-based residual ocean circulation pressure uncertainties 5) 'leakage': apparent gravity-based bottom pressure due to continental leakage 6) 'model_ocean': background-model ocean circulation pressure 7) 'model_atmosphere': background-model atmospheric surface pressure These Level-3 products are visualized at GFZ's web portal GravIS (http://gravis.gfz-potsdam.de). Link to data products: ftp://isdcftp.gfz-potsdam.de/grace/GravIS/GFZ/Level-3/OBP

Description: We propose to investigate the structure and evolution of the Main Pamir Thrust (MPT) with a high-density seismological array. The MPT, with its surface expression along the east-west trending Alai Valley, marks the northern boundary of the Pamir. The Alai Valley, separating the Pamir and the Tien Shan, constitutes the last vestige of a formerly continuous basin that linked the Tarim and the Tajik Basins. The MPT manifests itself as a place of high seismic activity with frequently occurred disastrous earthquakes. The array is about 50 km long, consisted of 90 three-component geophones (stations G?? and C??) and 10 Trillium-Compact seismometers (stations T??), and equipped with 100 CUBE dataloggers. We will construct a high-resolution receiver function profile to image the MPT and accurately locate the local earthquakes associated with the MPT. Funded by BMBF, within the framework of CaTeNA project – Climatic and Tectonic Natural Hazards in Central Asia. Waveform data are available from the GEOFON data centre, under network code 7A and are embargoed until Jan 2024.

Description: The CNES/GRGS RL04 Earth gravity models are a set of gravity field solutions based on GRACE and SLR data, provided at different time samplings: (A) CNES/GRGS RL04 time series (A/1) A monthly GRACE+SLR time series of gravity field models (A/2) A 10-day GRACE+SLR time series of gravity field models (B) A mean gravity model EIGEN-GRGS.RL04.MEAN-FIELD, computed from the monthly RL04 GRACE+SLR time series and from GOCE data. (A) CNES/GRGS RL04 time series DATA: The data from the Star Camera Assembly (SCA), Accelerometer (ACC), K-Band Ranging (KBR) and GPS receiver are used. The KBR data is processed in the form of the relative velocity between the spacecrafts: K-Band Range-Rate (KBRR). In addition to the data from GRACE, the data from 5 SLR satellites are also used (Lageos, Lageos-2, Starlette, Stella and Ajisai), in order to provide an accurate and consistent description of the very low degrees of the gravity field (mainly degrees 1 and 2). The version of the GRACE data used for RL04 is L1B-v2 for the ACC and GPS data, L1B-v3 for the SCA and KBR data. INVERSION METHOD: By contrast with the GRACE solutions in spherical harmonics provided by other groups, the CNES/GRGS solutions are not obtained by a simple Cholesky inversion. The normal matrices are first diagonalized, ordered by decreasing order of the Eigen values and only the best defined sets of linear combinations of the spherical harmonics are solved. More details can be found here: https://grace.obs-mip.fr/variable-models-grace-lageos/grace-solutions-release-04/rl04-products-description/ (B) EIGEN-GRGS.RL04.MEAN-FIELD mean model EIGEN-GRGS.RL04.MEAN-FIELD is a mean model of Earth's gravity field spherical harmonics coefficients, based on the RL04 version of the CNES/GRGS time series of monthly gravity field determinations from GRACE & SLR data. EIGEN-GRGS.RL04.MEAN-FIELD is complete to degree and order 300. Between degrees 1 and 90, it contains time-variable gravity (TVG) coefficients ; above degree 90, it is a static field. EIGEN-GRGS.RL04.MEAN-FIELD is based on GOCE-DIR5 for the part between degree 91 and 300. The TVG coefficients between degrees 1 and 90 are obtained from a regression on the GRGS-RL04-v1 monthly time series of solutions (2002/09 – 2016/06). For degrees 1 and 2 this TVG part is temporally extended to 1993/01-2019/02 through the use of a GRGS SLR-only solution based on the data of 5 SLR satellites (Lageos, Lageos-2, Starlette, Stella, Ajisai). Outside of the measurements period (1993/01-2019/02 for degrees 1 and 2, 2002/09-2016/06 for degrees 3 to 90), the gravity field is extrapolated in the following way: - for degrees 1 and 2, before 1993/01 : average slope based on historical SLR data, mean annual and semi-annual periodic signals based on their average value between 1993 and 2019 - for degrees 1 and 2, after 2019/02 : average slope & mean annual and semi-annual periodic signals (based on their average value between 1993 and 2019) - for degrees 3 to 90, before 2002/09 : zero-slope extrapolation, mean annual and semi-annual periodic signals based on their average value between 2002 and 2016 - for degrees 3 to 90, after 2016/06 : average slope & mean annual and semi-annual periodic signals (based on their average value between 1993 and 2019) More details can be found here: https://grace.obs-mip.fr/variable-models-grace-lageos/mean-fields/release-04/

Description: The software package “ClassifyStorms” version 1.0.1 performs a classification of geomagnetic storms according to their interplanetary driving mechanisms based exclusively on magnetometer measurements from ground. In this version two such driver classes are considered for storms dating back to 1930. Class 0 contains storms driven by Corotating or Stream Interaction Regions (C/SIRs) and class 1 contains storms driven by Interplanetary Coronal Mass Ejections (ICMEs). The properties and geomagnetic responses of these two solar wind structures are reviewed, e.g., by Kilpua et al. (2017, http://doi.org/10.1007/s11214-017-0411-3). The classification task is executed by a supervised binary logistic regression model in the framework of python's scikit-learn library. The model is validated mathematically and physically by checking the driver occurrence statistics in dependence on the solar cycle phase and storm intensity. A detailed description of the classification model is given in Pick et al. (2019) to which this software is supplementary material. Under “Files” you can download ClassifyStorms-V1.0.1.zip, which contains the jupyter notebook “ClassifyStorms.ipynb” (https://jupyter.org/) and the python modules “Imports.py”, “Modules.py” and “Plots.py”. Check for an up-to-date release of the software on GitLab via https://gitext.gfz-potsdam.de/pick/public/ClassifyStorms (under Project, Releases). The “Readme.md” file provides all information needed to run or modify “ClassifyStorms” from the GitLab source. The software depends on the input data set “Input.nc”, an xarray Dataset (http://xarray.pydata.org/en/stable) saved in NetCDF format (https://www.unidata.ucar.edu/software/netcdf), which you can also download under “Files”. It contains 1. the HMC index: a three-hour running mean with weights [0.25,0.5,0.25] of the original Hourly Magnetospheric Currents index (HMC index, http://doi.org/10.5880/GFZ.2.3.2018.006). 2. the geomagnetic observatory data: vector geomagnetic disturbances from 34 mid-latitude observatories during 1900-2015 in the Cartesian Centered Dipole coordinate system. The original observatory data was downloaded from the WDC for Geomagnetism, Edinburgh (http://www.wdc.bgs.ac.uk/) and processed as described in section 2.1 of Pick et al. (2019). 3. the “reference” geomagnetic storms: universal time hours of 868 geomagnetic storm peaks together with their interplanetary drivers (class labels 0 or 1, see above) as described in section 2.2 of Pick et al., 2019. These events are taken from published lists (Jian et al., 2006a, 2006b, 2011; Shen et al., 2017; Turner et al., 2009), which are gathered in the separate ASCII file “ReferenceEvents.txt” (under “Files”) for a quick overview. 4. additional quantities for plotting: time series of Kp (since 1932) and Dst (since 1957) geomagnetic indices from the WDC for Geomagnetism, Kyoto (http://wdc.kugi.kyoto-u.ac.jp/wdc/Sec3.html) as well as the yearly mean total sunspot number from WDC-SILSO, Royal Observatory of Belgium, Brussels (http://sidc.be/silso/datafiles). The output of ClassifyStorms is "StormsClassified.csv" (under “Files”). This table lists the Date (Year-Month-Day) and Time (Hour:Minutes:Seconds) of 7546 classified geomagnetic storms together with the predicted interplanetary driver class label (0 or 1) and the corresponding probability (between 0 and 1). Version history: 20 Sep 2019: Version 1.0.1: Correction of plotting mistake in Figure m / Figure S4 (see gitlab repository for details)

Description: Surrogate playground is an automated machine learning approach written for rapidly screening a large number of different models to serve as surrogates for a slow running simulator. This code was written for a reactive transport application where a fluid flow model (hydrodynamics) is coupled to a geochemistry simulator (reactions in time and space) to simulate scenarios such as underground storage of CO2 or hydrogen storage for excess energy from wind farms. The challenge for such applications is that the geochemistry simulator is typically slow compared to fluid dynamics and constitutes the main bottleneck for producing highly detailed simulations of such application scenarios. This approach attempts to find machine learning models that can replace the slow running simulator when trained on input-output data from the geochemistry simulator. The code may be of more general interest as this prototype can be used to screen many different machine learning models for any regression problem in general. To illustrate this it also includes a demonstration example using the Boston housing standard data-set.

Description: Global spherical harmonic paleomagnetic field model LSMOD.2 describes the magnetic field evolution from 50 to 30 ka BP based on published paleomagnetic sediment records and volcanic data. It is an update of LSMOD.1, with the only difference being a correction to the geographic locations of one of the underlying datasets. The time interval includes the Laschamp (~41 ka BP) and Mono Lake (~34 ka BP) excursions. The model is given with Fortran source code to obtain spherical harmonic magnetic field coefficients for individual epochs and to obtain time series of magnetic declination, inclination and field intensity from 49.95 to 30 ka BP for any location on Earth. For details see M. Korte, M. Brown, S. Panovska and I. Wardinski (2019): Robust characteristics of the Laschamp and Mono lake geomagnetic excursions: results from global field models.

Description: This is a Level 3 data daily file product from various scientific and utility sensors on board of the `LEO' satellite 'CHAMP' with magnetic field data given by a time resolution of 1 Hz. Thise Level 3 data type is build to hold and merge finally corrected data, focusing on mature data calibration and corrections -- as well as internal consistency. This Level 3 data product is intended to supersede the various Level 2 versions with calibrated magnetic field readings from the CHAMP mission distributed hitherto and should be fitted for scientific use, assembling time series of scalar magnetic field values (but not directly readings from the scalar Overhauser sensor), vector magnetic field data from the boom-mounted Fluxgate 'FGM' sensors and attitude data from the ('ASC') boom-mounted Star Cameras. The vector data are given both in the satellite-bound sensor ('FGM') system and the Earth Centered Earth Fixed local 'NEC' (North-East-Center) system. The attitude time series, processed and cleaned, are represented by quaternions describing the satellite attitude related to the celestial system. The readings of the scalar OVM (Overhauser) absolute magnetometer at the top of the boom are not supplied directly, but were used during calibration of the vector magnetometer readings. The files with daily time coverage are in the (binary and self-describing) 'CDF' file format and accompanied, beside the generic 'CDF'-format timestamp, by the satellite's geocentric positions and utility information like quality flags.

Description: Post-processed GRACE/GRACE-FO spherical harmonic coefficients of GFZ RL06 Level-2 GSM products representing an estimate of Earth's gravity field variations during the specified timespan. Post-processing steps comprise: (1) subtraction of a long-term mean field; (2) optionally, decorrelation and smoothing with VDK filter (anisotropic filter taking the actual error covariance information of the underlying GSM coefficients into account, see Horvath et al. (2018)); (3) replacement of coefficient C20 and its uncertainty by values estimated from Satellite Laser Ranging (SLR); (4) subtraction of linear trend caused by Glacial Isostatic Adjustment (GIA) as provided by a numerical model; (5) insertion of coefficients of degree 1; and (6) removal of estimated signal with 161 days period. These coefficients represent signals caused by water mass redistribution over the continents and in the oceans. These post-processed GRACE/GRACE-FO GSM products are denoted as Level-2B products.

Description: These data are supplementary material to "3D Modelling of Vertical Gravity Gradients and the delimitation of tectonic boundaries: The Caribbean oceanic domain as a case study" (Gómez-García et al., 2019). This dataset contains information about the structure of the Caribbean oceanic crust, based on the modelling of the Vertical Gravity Gradients, which are gravity derivatives especially sensitive to density contrasts in the upper layers of the Earth.

Description: The distribution of data records for the maximum horizontal stress orientation S_Hmax in the Earth’s crust is sparse and very unequally. To analyse the stress pattern and its wavelength and to predict the mean S_Hmax orientation on regular grids, statistical interpolation as conducted e.g. by Coblentz and Richardson (1995), Müller et al. (2003), Heidbach and Höhne (2008), Heidbach et al. (2010) or Reiter et al. (2014) is necessary. Based on their work we wrote the Matlab® script Stress2Grid that provides several features to analyse the mean S_Hmax pattern. The script facilitates and speeds up this analysis and extends the functionality compared to the publications mentioned before. This script is the update of Stress2Grid v1.0 (Ziegler and Heidbach, 2017). It provides two different concepts to calculate the mean S_Hmax orientation on regular grids. The first is using a fixed search radius around the grid points and computes the mean S_Hmax orientation if sufficient data records are within the search radius. The larger the search radius the larger is the filtered wavelength of the stress pattern. The second approach is using variable search radii and determines the search radius for which the standard deviation of the mean S_Hmax orientation is below a given threshold. This approach delivers mean S_Hmax orientations with a user-defined degree of reliability. It resolves local stress perturbations and is not available in areas with conflicting information that result in a large standard deviation. Furthermore, the script can also estimate the deviation between plate motion direction and the mean S_Hmax orientation. The script is fully documented by the accompanying WSM Technical Report 19/02 (Ziegler and Heidbach, 2019) which includes a changelog in the beginning.

Description: The database presented here contains radiogenic neodymium and strontium isotope ratios measured on both terrestrial and marine sediments. It was compiled to help assessing sediment provenance and transport processes for various time intervals. This can be achieved by either mapping sediment isotopic signature and/or fingerprinting source areas using statistical tools (e.g. Blanchet, 2018b, 2018a). The database has been built by incorporating data from the literature and the SedDB database and harmonizing the metadata, especially units and geographical coordinates. The original data were processed in three steps. Firstly, a specific attention has been devoted to provide geographical coordinates to each sample in order to be able to map the data. When available, the original geographical coordinates from the reference (generally DMS coordinates, with different precision standard) were transferred into the decimal degrees system. When coordinates were not provided, an approximate location was derived from available information in the original publication. Secondly, all samples were assigned a set of standardized criteria that help splitting the dataset in specific categories. We defined categories associated with the sample location ("Region", "Sub-region", "Location", which relate to location at continental to city/river scale) or with the sample types (terrestrial samples – “aerosols”, “soil sediments”, “river sediments”, “rocks” - or marine samples –“marine sediment” or “trap sample”). Thirdly, samples were discriminated according to their deposition age, which allowed to compute average values for specific time intervals (see attached table "Age_determination_Sediment_Cores_V2.txt"). A first version of the database was published in September 2018 and presented data for the African sector. A second version was published in April 2019, in which the dataset has been extended to reach a global extent. The dataset will be further updated bi-annually to increase the geographical resolution and/or add other type of samples. This dataset consists of two tab separated tables: "Dataset_Nd_Sr_isotopes_V2.txt" and "Age_determination_Sediment_Cores_V2.txt". "Dataset_Nd_Sr_isotopes_V2.txt" contains the assembled dataset of marine and terrestrial Nd and/or Sr concentration and isotopes, together with sorting criteria and geographical locations. "Age_determination_Sediment_Cores_V2.txt" contains all background information concerning the determination of the isotopic signature of specific time intervals (depth interval, number of samples, mean and standard deviation). Column headers are explained in respective metadata comma-separated files. A full reference list is provided in the file “References_Database_Nd_Sr_isotopes_V2.rtf”. Finally, R code for mapping the data and running statistical analyses is also available for this dataset (Blanchet, 2018b, 2018a).

Description: Earth's magnetic field vector time series from `LEO' satellite 'CHAMP' for the 'CHAMP' mission period in high, unaveraged 50 Hz time resolution, using measurements from the FGM vector magnetometers and `ASC' Star Sensors on the mid-boom optical bench. The vector data are corrected and calibrated (by using the Overhauser scalar magnetometer as reference). The magnetic field vector data are given both in the satellite-bound sensor (`FGM') system and in the Earth Centered Earth Fixed local `NEC' (North-East-Center) system. For the latter the attitude time series (`ASC'), processed and cleaned, represented by quaternions describing the satellite attitude related to the celestial system, were used for the transformation. The files with daily time coverage are in the (binary and self-describing) `CDF' file format and accompanied, beside the `CDF'-format generic timestamp, by the satellite's geocentric positions and utility information like quality flags.

Description: Spherical harmonic coefficients that are zero over the continents, and provide the anomalous simulated ocean bottom pressure that includes non-tidal air and water contributions elsewhere during the specified timespan. These coefficients differ from GLO (or GAC) coefficients over the ocean domain by disregarding upper air density anomalies. The anomalous signals are relative to the mean field from 2003-2014.

Description: BayHunter is an open source Python tool to perform an McMC transdimensional Bayesian inversion of receiver functions and/ or surface wave dispersion. It is inverting for the velocity-depth structure, the number of layers and noise parameters (noise correlation and amplitude). The forward modeling codes are provided within the package, but are easily replaceable with own codes. It is also possible to add (completely different) data sets. The BayWatch module can be used to live-stream the inversion while it is running: this makes it easy to see how each chain is exploring the parameter space, how the data fits and models change and in which direction the inversion progresses.

Description: In Bindi et al. (2019) a harmonized local magnitude scale across Europe has been derived using data disseminated by network operators through the European Integrated Data Archive (EIDA). This data set contains regionalized non-parametric attenuation tables, attenuation corrections to the parametric model and station corrections for both non-parametric and parametric models for more than 2000 stations in Europe. Regionalization has been performed considering six different regions covering Europe and the polygons defining them are also provided. Data are subject to updates that can be triggered by the availability of new and substantial input data (reviewed earthquake catalogues and/or new waveforms). Each update will be released with a new version of the data. The data are provided in ASCII format (.csv).

Description: Extensive passive seismic monitoring was carried out between September 2017 and September 2018 over the Los Humeros geothermal field in Mexico. This experiment, in addition to several geophysical, geological, and geochemical surveys was conducted in the framework of the European H2020 and Mexican CONACyT-SENER project GEMex for a better understanding of the structures and behavior of the local geothermal system currently under exploitation, and for investigating future development areas. 25 broadband stations (22 Trillium C-120s and 3 Trillium C-20 PH) recording at 200 Hz, and 20 short period stations (Mark L-4C-3D) recording at 100 Hz comprised the network which is sub-divided into two sub-networks. An inner and denser (~1.6-2 km inter-station distance) pseudo-rhomboidal array (27 stations) was laid out to cover the producing zone and retrieve local seismicity mainly associated to injection and production operations, and to comply with beamforming of ambient noise and time reverse imaging techniques. An outer and sparser (~5 km minimum spacing) array was placed at around 30 km radius surrounding the inner network, and was mainly dedicated to larger scale imaging techniques, such as seismic ambient noise tomography, and regional earthquakes tomography. The GEMex project is supported by the European Union’s Horizon 2020 programme for Research and Innovation under grant agreement No 727550 and the Mexican Energy Sustainability Fund CONACYT-SENER, project 2015-04-68074. Waveform data are available from the GEOFON data centre, under network code 6G, and are embargoed until January 2023.

Description: We present SCOTER, an open-source Python programming package that is designed to relocate multiple seismic events by using direct P- and S-wave station correction terms. The package implements static and shrinking-box source-specific station terms techniques extended to regional and teleseimic distances and adopted for probabilistic, non-linear, global-search location for large-scale multiple-event location. This program provides robust relocation results for seismic event sequences over a wide range of spatial and temporal scales by applying empirical corrections for the biasing effects of 3-D velocity structure. Written in the Python programming language, SCOTER is run as a stand-alone command-line tool (requiring no knowledge of Python) and also provides a set of sub-commands to develop required input files (e.g. phase files, travel-time grid files, configuration) and export relocation results (such as hypocenter parameters, travel-time residuals) in different formats -- routine but non-trivial tasks that can consume much user time. This package can be used for relocating data sets in local, regional, and teleseimic scales.

Description: Time series of processed, cleaned attitude readings in quaternion format of the two boom-mounted 'ASC' star sensors of the 'LEO' satellite 'CHAMP', describing the satellite system attitude in respect to the celestial background. The nominal time resolution of the time series in the 'ASCII'-file listing is 1 Hz.

Description: Electron density and electron temperature time series from 'LEO' satellite 'CHAMP' for the CHAMP mission period at satellite position in low time resolution of 15 second and given in daily files. This are processed readings from the Planar Langmuir probe, which, in normal flight mode, was exposed in flight direction at the front of the `CHAMP' satellite body. The files are formatted as simple 'ASCII'-listings with white-space delimited columns.

Description: This publication contains the supplementary data set to Mikolaj et al. "Resolving geophysical signals by terrestrial gravimetry: a time domain assessment of the correction-induced uncertainty" (2019, JGR-Solid Earth). The aim of the article is to estimate the uncertainty of terrestrial gravity corrections applied to resolve small-scale gravity effects. The uncertainty of the gravity corrections is assessed using various models of the tidal effect, large-scale hydrology, non-tidal ocean loading, and atmosphere. Taken into account are widely recognized models with global spatial coverage, sufficient temporal resolution and coverage, and available to the public for research purposes. The uncertainty is expressed in terms of a root-mean-square and mean-absolute error of the deviations between all available models. The data set comprises models for 11 sites worldwide. The processing scripts are provided along with an explanatory file with all instructions for results reproduction and application of the uncertainty analysis for an arbitrary location.

Description: This data publication is part of the 'P³-Petrophysical Property Database' project, which has been developed within the EC funded project IMAGE (Integrated Methods for Advanced Geothermal Exploration, EU grant agreement No. 608553) and consists of a scientific paper, a full report on the database, the database as excel and .csv files and additional tables for a hierarchical classification of the petrography and stratigraphy of the investigated rock samples (see related references). This publication here provides a hierarchical interlinked petrographic classification according to standardized and internationally defined petrographic terms. The petrography or rock type classification scheme is structured based on a hierarchical subdivision with nine different ranks, where the rock description generally becomes more detailed with increasing rank of petrographic classification (based on the well database of the Geological Survey of Hessen, Germany: Hessisches Landesamt für Umwelt, Naturschutz, Umwelt und Geologie (HLNUG)). This hierarchical subdivision and the definitions of the petrographic terms are based on international conventions (e.g. Bates & Jackson 1987, Gillespie & Styles 1999, Robertson 1999, Hallsworth & Knox 1999, Bas & Streckeisen 1991, Schmid 1981, Fisher & Smith 1991). Furthermore, the classification corresponds to the subdivision provided by existing property data compilations such as e.g. Hantschel and Kauerauf (2009), Schön (2011), Rybach (1984) and Clauser and Huenges (1995). Petrographic classifications from rank 1 to rank 4 can usually be identified from macroscopic descriptions of well logs, cores and geological mapping. The petrographic classifications from rank 5 to rank 9 require additional information on the texture or grain size, the modal composition or the geochemistry etc., which can usually only be acquired by microscopic or comparable special investigations. Overall, the nine ranks cover a total of 1494 petrographic terms and thus goes well beyond other standardized catalogues (e.g. 'Simplified Lithology' in GeoSciML). The petrographic classification of a sample in P³ is based on the sample description within the original literature reference. A petrographic ID and a corresponding petrographic parental ID directly correlate the different classifications and their ranks.

Description: Spherical harmonic coefficients that represent anomalous contributions of the non-tidal atmosphere to the Earth's mean gravity field during the specified timespan. This includes the contribution of atmospheric surface pressure over the continents, the static contribution of atmospheric pressure to ocean bottom pressure elsewhere, and the contribution of upper-air density anomalies above both the continents and the oceans. The anomalous signals are relative to the mean field from 2003-2014.

Description: The data set contains mineral chemical analyses of 20 different copper bearing minerals and their corresponding hyperspectral spectra. The hyperspectral data were acquired with the HySpex system in a range of 400 – 2500 nm and are presented in a spectral library. Detailed information about the mineral specimen, sample area and geochemistry is presented in the data sheets and associated data description. The spectral library presented here is part of a bigger collection of spectral libraries.

Description: This report describes the passive seismic data acquired by the TOPASE network deployed over Rittershoffen geothermal field (Alsace, France). The monitoring period extends from March 2013 to November 2014, which includes the stimulation of the first well of the doublet, the drilling of the second well and well tests. These data were acquired using 31 Earth Data Loggers PR6-24 and MARK-SERCEL L-4C-3D 1 Hz seismometers of the Geophysical Instrument Pool Potsdam (GIPP), which were provided to the KIT-AGW-Geothermal research division.

Description: The present dataset is a high-resolution earthquake catalog for the region of the 2014 M8.1 Iquique earthquake sequence, in the Northern Chile subduction zone. Events in the catalog were obtained processing seismic waveforms from 〉100 permanent and temporary seismic stations. The list of waveform data sources used, the multistage automatic earthquake detection and location procedure implemented to build the catalog, as well as the patterns outlined by the seismicity are described in Soto et al. (2019).

Description: A seismic network was installed in the Helsinki capital area of Finland to monitor the response to a 6 km deep geothermal stimulation experiment in 2018. The Institute of Seismology, University of Helsinki (ISUH), installed these 100 geophones in addition to five surface broadband sensors and a 13-site borehole network deployed by the operating company. The stations operated for 106 days between 7 May and 20 August 2018 (day 127 to 232). The data set consists of raw CUBE-recorder data and converted MSEED data.

Description: The Alps are one of the best studied mountain ranges in the world, yet significant unknowns remain regarding their crustal structure and density distribution at depth. Previous published interpretations of crustal features within the orogen have been primarily based upon 2D seismic sections, and those that do integrate multiple geo-scientific datasets in 3D, have either focused on smaller sub-sections of the Alps or included the Alps, in low resolution, as part of a much larger study area. Therefore the generation of a 3D, crustal scale, gravity constrained, structural model of the Alps and their forelands at an appropriate resolution has been created here to more accurately describe crustal heterogeneity in the region. The study area of this work focuses on a region of 660 km x 620 km covering the vast majority of the Alps and their forelands are included, with the Central and Eastern Alps and the northern foreland being the best covered regions.

Description: An electric current is induced by the motion of electrical conducting seawater through the ambient geomagnetic field. The periodic oceanic tidal flow induces an electric current that emits periodical time-variable electromagnetic field signals. The radial component of the ocean tide induced magnetic field signals has successfully been extracted from magnetic field observations of the satellite missions CHAMP and Swarm. It is known that the amplitudes of these electromagnetic signals are modulated by, among other influences, variations of the electrical seawater conductivity distribution of the ocean. The electrical seawater conductivity in return depends on seawater temperature and salinity. In order to analyse the influence of variations in oceanic temperature and salinity, we modelled a complete set of monthly time slices of three dimensional global complex amplitudes of these electromagnetic field signals for the years 1990 to 2016. In order to analyse solely the influence of variations in the climate sensitive seawater temperature and salinity on the ocean tide induced magnetic field signals, the influences of the secular variation of the geomagnetic field and temporal variations in ocean tide transports have been neglected.

Description: We present a 3-D lithospheric-scale model covering the area of Germany that images the regional structural configuration. The model comprises 31 lithostratigraphic units: seawater, 14 sedimentary units, 14 crystalline crustal units and 2 lithospheric mantle units. The corresponding surfaces are integrated from previous studies of the Central European Basin System, the Upper Rhine Graben and the Molasse Basin, together with published geological and geophysical data. The model is a result of a combined workflow consisting of 3-D structural, gravity and thermal modelling applied to derive the 3-D thermal configuration. The top surface elevations and thicknesses of corresponding layers of the 3-D-D model are provided as ASCII files, one for each individual layer of the model. The columns in each file are identical: the Easting is given in the “X COORD (UTM Zone 32N)”, the Northing is in the “Y COORD (UTM Zone 32N)”, the top surface elevation of each layer is given as "TOP (m.a.s.l)", the thickness of each layer is given as "THICKNESS (m)".

Description: In Münchmeyer et al. 2019 magnitudes scales for Northern Chile have been derived with a focus on low uncertainties. The data set consists of three parts. First, a version of the IPOC catalog with the derived magnitude scales ML and MA and their uncertainties. Second, the attenuation functions for different waveform features. Third, the full matrix of features and the resulting single station magnitude predictions.The underlying IPOC catalog was obtained from Sippl et al. (2018). Detailed data description is provided in the README and in Münchmeyer et al. (2019) to which these data are supplementary material.

Description: EMMA – End Member Modelling Analysis of grain-size data is a technique to unmix multimodal grain-size data sets, i.e., to decompose the data into the underlying grain-size distributions (loadings) and their contributions to each sample (scores). The R package EMMAgeo contains a series of functions to perform EMMA based on eigenspace decomposition. The data are rescaled and transformed to receive results in meaningful units, i.e., volume percentage. EMMA can be performed in a deterministic and two robust ways, the latter taking into account incomplete knowledge about model parameters. The model outputs can be interpreted in terms of sediment sources, transport pathways and transport regimes (loadings) as well as their relative importance throughout the sample space (scores).

Description: Spherical harmonic coefficients that represent anomalous contributions of the non-tidal dynamic ocean to ocean bottom pressure during the specified timespan. The anomalous signals are relative to the mean field from 2003-2014.

Description: This data publication is part of the 'P³-Petrophysical Property Database' project, which was developed within the EC funded project IMAGE (Integrated Methods for Advanced Geothermal Exploration, EU grant agreement No. 608553) and consists of a scientific paper, a full report on the database, the database as excel and .csv files and additional tables for a hierarchical classification of the petrography and stratigraphy of the investigated rock samples (see related references). This publication here provides a hierarchical interlinked stratigraphic classification according to the chronostratigraphical units of the international chronostratigraphic chart of the IUGS v2016/04 (Cohen et al. 2013, updated) according to international standardisation. As addition to this IUGS chart, which is also documented in GeoSciML, stratigraphic IDs and parent IDs were included to define the direct relationships between the stratigraphic terms. The P³ database aims at providing easily accessible, peer-reviewed information on physical rock properties relevant for geothermal exploration and reservoir characterization in one single compilation. Collected data include hydraulic, thermophysical and mechanical properties and, in addition, electrical resistivity and magnetic susceptibility. Each measured value is complemented by relevant meta-information such as the corresponding sample location, petrographic description, chronostratigraphic age and, most important, original citation. The original stratigraphic and petrographic descriptions are transferred to standardized catalogues following a hierarchical structure ensuring intercomparability for statistical analysis, of which the stratigraphic catalogue is presented here. These chronostratigraphic units are compiled to ensure that formations of a certain age are connected to the corresponding stratigraphic epoch, period or erathem. Thus, the chronostratigraphic units are directly correlated to each other by their stratigraphic ID and stratigraphic parent ID and can thus be used for interlinked data assessment of the petrophysical properties of samples of an according stratigraphic unit.

Description: Increasing flood losses over the last decades emphasize the need towards significantly improved and more efficient flood risk management. One key requirement is reliable risk assessment in conjunction with consistent flood loss modeling. Current risk assessments and flood loss estimations for Europe are until now based on regional approaches using deterministic depth-damage function and do rarely report associated uncertainties. To reduce these shortcomings, we present the results of a novel, consistent approach based on the Bayesian Network Flood Loss Estimation MOdel for the private sector (BN-FLEMOps). The dataset is consistent in terms of the input data used to drive the model and because we use the same vulnerability model to derive the flood loss estimation. Essential inputs for any flood loss estimation are hazard (usually water depth), asset (value of objects at risk) and flood experience parameters. The hazard input was given by a European inundation scenario for a continent-wide flood with 100 years return period (Alfieri et al., 2014). Asset values were computed following the the approach by Huizinga et al. (2017) and the flood experience was derived using the database of the Dartmouth Flood Observatory (DFO) (Brakenridge, 2018).

Description: The data set contains mineral chemical analyses of 32 rare earth element (REE) -bearing minerals (REMin) and rare-earth oxides (REO) and their corresponding hyperspectral spectra. The hyperspectral data was acquired with the HySpex system in a range of 400 – 2500 nm and is presented in a spectral library. The two Rare Earth Element (REE) libraries consist of the spectra of 16 rare earth oxides powders (REO) and 14 REE-bearing minerals (REMin). In addition, it contains the spectra of niobium- and tantalum oxide, two elements technically not part of the REEs.

Description: We have installed 20 new Global Positioning System (GPS) markers in the West Pamir and the Tajik Depression and measured 25 markers once a year between 2013 and 2016 in survey mode. The stations are positioned along two dense NW-SE oriented profiles with an average spacing of 5-10 km. The profiles cross the Darvaz and the Vakhsh/Ilyak fault and thus monitor the recent slip of these two profiles, which are expected to accommodate the gravity-driven westward extrusion of the West Pamir into the Tajik Depression. Some of the stations include millimeter to centimeter offsets potentially caused by the 2015 Mw7.2 Sarez, Pamir, earthquake.

Description: The profile 9N was recorded in 1988 as part of the DEKORP project, the German deep seismic reflection program. The seismic survey of the ca. 92-km long line 9N was conducted to investigate the deep crustal structure of the northern Upper Rhine Graben with high-fold near-vertical incidence vibroseis acquisition. The objectives of the survey were to delineate the geometry of the major faults, which control the graben subsidence, to map the geometry of deep crustal reflection patterns and to reveal variations of the seismic signature of the lower crust in the context of rift formation. The first results were discussed by Wenzel et al. (1991), summarized by Brun et al. (1992) and supplemented by many other researches. Since the Eocene the Upper Rhine Graben has represented an active rift system. It obliquely intersects the Saxothuringian and the Moldanubian domains, which are separated by the NW vergent and dextrally sheared Lalaye‐Lubine‐Baden‐Baden fault. In the northern Vosges and Black Forest massifs the shear zone is characterized by low-grade Devonian metasediments. The profile starts in the crystalline Odenwald in the east, intersects the Tertiary and Quaternary fill of the Rhine Graben and ends in the late Palaeozoic sequences of the Saar-Nahe Basin in the west, where it crosses the Permian rhyolitic Donnersberg intrusion. The profile 1C creates a continuation of the survey to the west. The seismic section of 9N shows different crustal structures on both sides of the graben and some indications of dipping reflections in the mantle on the western side, which could refer to the genesis of the Upper Rhine Graben.

Description: Multi-resolution exposure model for seismic risk assessment in Turkmenistan. The model has been developed according to the methodology outlined in Pittore, Haas and Silva (2020) "Variable resolution probabilistic modeling of residential exposure and vulnerability for risk applications", Earthquake Spectra. The model is aggregated over a Central Voronoidal Tessellation (CVT) composed of geo-cells covering the territory of Turkmenistan (provided as a separate file). The model prior is based on user-elicited knowledge. The following specific modelling parameters have been employed: Two exposure models are provided, with prior strength pw 10 and 100. Both models have epsilon=0.001 (see publication indicated in the metadata for details on the modelling process) For each geo-cell the model includes the expected number of buildings , total occupancy and replacement cost for each of the 15 building types defined in the EMCA taxonomy (see Pittore et al, 2019b), plus the buildings that are belonging to other, non specified typologies (described by building type OTH). Each geo-cell also includes the area of the geo-cell itself in squared km. The data package contains three components: 1) exposure models in .csv 2) exposure models in .xml - the file is encoded in NRML 0.5 format and is compatible with the GEM openquake processing engine 3) shapefile of the tessellation that aggregates the exposure model. The field "cell_id" is the linkage with the exposure models

Description: The complete network consisted of eleven stations deployed on the island of Fogo, Cape Verde. Eight of the stations formed an array with an aperture of 700 m, deployed in the south of the island near the village of Achada Furna. Seven of the array stations were equipped with 3-component 4.5 Hz geophones, one with a Trillium Compact (broad-band) sensor. The remaing three stations were distributed across the island and equipped with Trillium Compact sensors. Data were recorded continuously from October 2015 to December 2016 with a sample rate of 200 Hz. Due to limited data storage, there are four recording gaps (20/12/2015-14/01/2016; 28/03/2016-04/04/2016; 17/06/2016-18/07/2016; 01/10/2016-18/10/2016). The network served as a pilot study for the more comprehensive study, FoMaPS, from 2017 to 2018 (FDSN code 9J), involving station deployments on Fogo and Brava.

Description: Petrophysical properties are key to populate numerical models of subsurface process simulations and for the interpretation of many geophysical exploration methods. They are characteristic for specific rock types and may vary considerably as a response to subsurface conditions (e.g. temperature and pressure). Hence, the quality of process simulations and geophysical data interpretation critically depend on the knowledge of in-situ physical properties that have been measured for a specific rock unit. Inquiries for rock property values for a specific site might become a very time-consuming challenge given that such data are (1) spread across diverse publications and compilations, (2) heterogeneous in quality and (3) continuously being acquired in different laboratories worldwide. One important quality factor for the usability of measured petrophysical properties is the availability of corresponding metadata such as the sample location, petrography, stratigraphy, or the measuring method, conditions and authorship. The open-access database presented here aims at providing easily accessible, peer-reviewed information on physical rock properties in one single compilation. As it has been developed within the scope of the EC funded project IMAGE (Integrated Methods for Advanced Geothermal Exploration, EU grant agreement No. 608553), the database mainly contains information relevant for geothermal exploration and reservoir characterization, namely hydraulic, thermophysical and mechanical properties and, in addition, electrical resistivity and magnetic susceptibility. The uniqueness of this database emerges from its coverage and metadata structure. Each measured value is complemented by the corresponding sample location, petrographic description, chronostratigraphic age and original citation. The original stratigraphic and petrographic descriptions are transferred to standardized catalogues following a hierarchical structure ensuring intercomparability for statistical analysis. In addition, information on the experimental set-up (methods) and the measurement conditions are given for quality control. Thus, rock properties can directly be related to in-situ conditions to derive specific parameters relevant for modelling the subsurface or interpreting geophysical data.

Description: The data set contains mineral chemical analyses of 37 different surface materials from the copper-gold-pyrite mine Apliki in the Republic of Cyprus and their corresponding hyperspectral spectra. The field samples were sampled in March 2018 in cooperation of the Cyprus Geological Survey Department of the Republic of Cyprus (GSD) and the German Research Centre for Geosciences (GFZ). The hyperspectral data was acquired with the HySpex system in a range of 400 – 2500nm and is presented in a spectral library. Detailed information about the mineral specimen, sample area and geochemistry is presented in the data sheets.

Description: Our understanding of the effects of ice on channel morphodynamics and bedload transport in northern rivers, frozen for several months, are hindered by the difficulties of ‘seeing’ through the ice. We use continuous seismic records of a small network at the Sävar River in northern Sweden to interpret processes and quantify water and sediment fluxes. We apply a seismic inversion approach to determine seasonal differences in hydraulics and bedload sediment transport under ice-covered vs. open-channel flow conditions and provide a first-order estimation of sediment flux in that Fennoscandian river. Analysis of seismic signals of ice-cracking support our visual interpretation of ice break-up timing and the main ice break-up mechanism as thermal rather than mechanical. Waveform data are available from the GEOFON data centre, under network code 8E.

Description: This dataset is related to the question of whether communities inside of certified "International Dark Sky Places" have different levels of lighting change in comparison to communities of similar size that are located further away. It is a supplement to Coesfeld et al. (2019), in which this question was examined for the time period 2012-2018. This dataset contains the boundaries of the analysis areas (i.e. community boundaries) in the directory "Coordinates and Graphs". These boundaries are stored as polygons in plain text format. Additional data related to the publication (e.g. Excel tables containing summary data of measured lighting trends for each community) are also included. Details of the data are available in the data description file.

Description: The interactive web page contains supplementary information for a publication by Hensch et al. 2019: "Deep low-frequency earthquakes reveal ongoing magmatic recharge beneath Laacher See Volcano (Eifel, Germany)". Details on the analysis of three tectonic and nine deep low-frequency earthquakes are given, including parameter results, error estimates, and figures. The analysis has been performed using the Grond software package (Heimann et. al 2018).

Description: This data set provides data from subduction zone earthquake experiments and analysis described in Rosenau et al. (2019). In the experiments analogue seismotectonic scale models of subduction zones characterized by two seismogenic asperities are used to study the interaction of asperities over multiple seismic cycles by means of static (Coulomb failure) stress transfer. Various asperity geometries (lateral/along-strike of the subduction zone distance and vertical/across-strike of the subduction zone offset) are tested on their effect on recurrence pattern of simulated great (M8+) earthquakes. The results demonstrate the role of stress coupling in the synchronization of asperities leading to multi-asperity M9+ events in nature. The data set contains time series of experimental surface velocities from which analogue earthquakes are detected and classified into synchronized events and solo events. The latter are subcategorized into main events and aftershocks and into normal and thrust events. An analogue earthquake catalogue lists all categorized events of the 12 experiments used for statistical analysis. Moreover, results from elastic dislocation modelling aimed ate quantifying the stress coupling between the asperities for the various geometries are summarized. Basic statistics of classified events (e.g. percentage of categorized events, coefficient of variation in size and recurrence time etc.) are documented. Matlab scripts are provided to visualize the data as in the paper.

Description: Radiance Light Trends is a GIS web application that is designed to quickly display information about radiance trends at a specific location (available online at https://lighttrends.lightpollutionmap.info). It uses data from two satellite systems, DMSP-OLS and VIIRS DNB, with data processing by NOAA. New VIIRS layers are added automatically as soon as NOAA makes them available to public.

Description: The scripts and workflow are supplementary material to "3D Modelling of Vertical Gravity Gradients and the delimitation of tectonic boundaries: The Caribbean oceanic domain as a case study" (Gómez-García et al., 2019). The codes include the calculation of the VGG response of a 3D lithospheric model, in spherical coordinates, using the software Tesseroids (Uieda, 2016). The "Readme_Workflow_2019_002.pdf" file provide very detail information about the structure of this repository, as well as the step-by-step for the scripts execution, and the list of the requiered software for the correct workflow performance.

Description: The determination of the global gravity field has gained momentum due to high accuracy satellite-derived observations and development of forward gravity modelling. Forward modelling computes the global gravitational field from mass distribution sources instead of actual gravity measurements and helps improving and complementing the medium to high frequency components of the global gravity field models. In this study, we approximate the global gravity potential of the Earth’s upper crust based on ellipsoidal approximation and a mass layer concept. Lateral density variations within a sequence of thin volumetric shells bounded by confocal lower and upper shell ellipsoids are used in the computation of the ellipsoidal harmonic coefficients which are then transformed into spherical harmonic coefficients on the Earth’s surface in the final step. The main outcome of this research is a spectral representation of the Earth’s upper crust’s gravitational potential, computed up to degree and order 3660 in terms spherical harmonic coefficients (ROLI_EllApprox_SphN_3660).

Description: Geysers are accessible sites of hot pots, springs and pools that regularly erupt. To investigate the frequency and dynamics of water eruptions we setup a local broadband seismic network at Strokkur geyser, Iceland. The experiment was running for 1 year, from June 2017 to June 2018. Four broadband seismic stations (Nanometrics Trillium Compact Posthole 20s) were buried 30-40 cm deep in the ground at a distance of 39 m (G4, SE), 47 m (G3, SW), 43m (G2, N) and 96 m (G1, NE) from the center of the Strokkur pool geyser. Regular visits and 2-month interval battery replacement allowed to power the stations without solar panels, therewith limiting visibility and site impact. From this data we picked a catalog of 73,466 eruptions, that are statistically further evaluated in Eibl et al. (2020), allowing to distinguish 50,135 single eruptions, and over 20,000 multiplet eruptions (i.e. several eruptions in close succession). The mean waiting time after an eruption at Strokkur linearly increased from 3.7 to 16.4 min for single and multiplets, respectively. This data publications releases the catalog of 73,466 eruptions.

Description: Since the early 1990s, the GFZ has operated a global GNSS station network with currently about 70 stations for precise satellite clock & orbit determination, realization of the terrestrial reference frame, radio occultation measurements or studies on crust dynamics. A subset of these stations contributes also to the tracking networks of the International GNSS Service (IGS) and the EUREF Permanent GNSS Network (EPN). Other stations contribute to GFZ observatories (IPOC, DESERVE, TERENO), to the GPS Atmosphere Sounding Project (GASP), to WMO Global Climate Observing System Reference Upper-Air Network (GRUAN) or to other external cooperations.

Description: P-phase arrival times automatically created by the SeisComP3 (https://www.seiscomp3.org/) software at the GFZ scanning all stations available in real-time at GEOFON Data Centre and listed in the contributors list. Data have been used in the publication by Steed et al 2018 to test the new CsLoc method, sent in relatime to EMSC using the HMB application (Heinloo, 2016, http://doi.org/10.5880/GFZ.2.4.2016.001). The data sets includes P-phases from 2016 and 2017.

Description: Due to poor global data coverage and large data uncertainties, Holocene magnetic field models are of limited resolution. The most commonly used modeling approach is based on an inversion strategy based on truncated spherical harmonics and does not provide uncertainty estimates. Focusing on snapshots of the magnetic field, the CORBASS (CORrelation Based Archeomagnetic SnapShot model) algorithm breaks with the tradition of truncated models by implementing a field model based on Gaussian processes (GP). This way regions covered well by data can be modeled with higher resolution, while areas of poor coverage remain uncertain. The GP approach provides a posterior covariance, together with the (mean) field model, which naturally serves as model uncertainty estimate. A full description of the model can be found in 〈TODO: reference to our upcoming paper〉. CORBASS is provided as a collection of python scripts together with an install file for a conda environment. This way CORBASS can be easily set up under various system dependencies. The gitlab page provides installation instructions, documentation and example notebooks and can be found under "Externe Ressourcen" in the left panel of this page.

Description: Geochemical models are used to seek answers about composition and evolution of groundwater, spill remediation, viability of geothermal resources and other important geoscientific applications. To understand these processes, it is useful to evaluate geochemical model response to different input parameter combinations. Running the model with varying input parameters creates a large amount of output data. It is a challenge to screen this data from the model to identify the significant relationships between input parameters and output variables. For addressing this problem we developed a Visual Analytics approach in an ongoing collaboration between Geoinformatics and Hydrogeology sections of GFZ German Research Centre for Geosciences. We implement our approach as an interactive data exploration tool called the GCex. GCex is a Visual Analytics approach and prototype that supports interactive exploration of geochemical models. It encodes many-to-many input/output relationships by the simple yet effective approach called Stacked Parameter Relation (SPR). GCex assists in the setup of simulations, model runs, data collection and result exploration, greatly enhancing the user experience in tasks such uncertainty and sensitivity analysis, inverse modeling and risk assessment. While in principle model-agnostic, the prototype currently supports and is tied to the popular geochemical code PHREEQC. Modification to support other models would not be complicated. GCex prototype was originally written by Janis Jatnieks at GFZ-Potsdam. It relies on Rphree (R-PHREEQC geochemical simulation model interface) written by Marco De Lucia at GFZ-Potsdam. A compatible version of Rphee is bundled with this installation.

Description: GRACE/GRACE-FO Level-3 product based on GFZ RL06 Level-2B products (Dahle & Murböck, 2019) representing Terrestrial Water Storage (TWS) anomalies provided at 1° latitude-longitude grids as defined over all continental regions except Greenland and Antarctica. The TWS anomaly grids are provided in NetCDF format divided into yearly batches. The files each contain four different variables: 1) 'tws': gravity-based TWS 2) 'std_tws': gravity-based TWS uncertainties 3) 'leakage': spatial leakage contained in TWS 4) 'model_atmosphere': background model atmospheric mass These Level-3 products are visualized at GFZ's web portal GravIS (http://gravis.gfz-potsdam.de). Link to data products: ftp://isdcftp.gfz-potsdam.de/grace/GravIS/GFZ/Level-3/TWS

Description: Operations such as time and coordinate conversions and data cleaning are routine tasks in geodesy and geophysics. Nevertheless, simple and efficient high-level functions to help those kinds of jobs are barely available, and has to be developed, again and again, by each student, engineer for each new project, and even by senior scientists. On another hand, Python became little by little within the last decade a well-used programming language in the academic world. Despite the fact that countless toolboxes already exist in Python for scientific purposes, none really exists for geodetic-oriented purposes. The geodeZYX toolbox aims to fill this gap. The objective of this toolbox, written in Python 3, is to provide a simple but useful and efficient set of functions to help geodesists and geophysicists to spend less time on the pre-processing steps and focus faster on their research, according to the KISS Principle.

Description: SMARTIE1 is a joint seismological experiment of the Karlsruhe Institute of Technology (KIT) and the Leipzig University. We installed in total 36 seismic stations as ring-like and profile-like measurements near to a single wind turbine (WT) at the Fraunhofer Institute for Chemical Technology (ICT) in Pfinztal, SW Germany, for 21 days. The main goals of this project are a better understanding of a single WT as a seismic source and the development of propagation models for the WT-induced seismic signals, depending on the geological properties. Waveform data are available from the GEOFON data centre, under network code X8, and are embargoed until end of 2019.

Description: The data available in this supplementary publication are: - A folder (2019-003_Corbi-et-al_Fig6.zip) containing: 1. top-view pictures (e.g. ‘lunayyr1_0025.JPG’) and displacement data obtained with MatPiv (e.g. ‘uun25.mat’ and ‘uvn25.mat’; dike parallel and orthogonal components; respectively) shown in figure 6 of Xu et al 2016. 2. a Matlab script (‘fig6_a_h.m’) that allows reproducing the same figure setup as in figure 6 panels a-h of Xu et al 2016. The thick red line highlights dike position. The background shading refers to dike orthogonal displacement. - A folder (2019-003_Corbi-et-al_PIV_data.zip) containing: 1. surface deformation data obtained with MatPiv. Each file (‘vel_fine_piv#.mat’) contains 4 elements (x, y, u, v) representing the coordinates and horizontal and vertical component of incremental velocity field organized in a 143 x 215 matrix; 2. the run_movie.m Matlab script. Running it the user can visualize the space-time evolution of cumulative surface displacement. The background shading refers to dike orthogonal component of displacement. The thick red line highlights dike position. - A folder (2019-003_Corbi-et-al_pictures.zip) containing the whole set of pictures from the experiment shown in Xu et al., 2016. - A movie (2019-003_Corbi-et-al_graben formation.mp4) obtained using the whole set of pictures (96 photos). The thick red line highlights dike position. The amount of dike opening is reported as header. - A movie (2019-003_Corbi-et-al_cum_displacement.mp4) showing the space-time evolution of cumulative surface displacement, where the background shading refers to dike orthogonal component of displacement. The thick red line highlights dike position.

Description: BEAT is an open-source software tool for the robust characterization of the temporal and spatial evolution of earthquake rupture processes. It uses kinematic rupture models that include low-parametric models like Moment Tensors but also complex high-parametric, finite-extent sources. In other words, BEAT allows studying earthquakes on a first-order level as points with location, size and mechanisms. In consecutive steps, the complexity of the source model may be increased by various details up to the potential to resolve rupture dimension, fault segmentation, slip-distribution and slip-history. The source model parameters and their uncertainties are estimated based on seismic waveforms, and/or geodetic observations like InSAR and GNSS data. Rapid forward modeling is enabled by using pre-computed Green's function databases, handled through the Pyrocko software library. Based on these, synthetic data are provided for arbitrary earthquake rupture models embedded in heterogeneous media. For an extensive exploration of the often high-dimensional model parameter space, BEAT offers a suite of sampling algorithms for high-standard Bayesian inference. The implementations of these sampling algorithms exploit the parallel architecture of modern computers for optimal performance. Finally, BEAT offers easy configuration and automatic visualization of relevant results. The software relies on functionality from PYROCKO (Heimann et al., 2017) and KITE (optionally, Isken et al., 2017).

Description: This data set includes 40 videos (+ 1 image) depicting the surface evolution of 39 experiments on crustal extension, as well as 4D CT imagery (figures and videos) of 6 of these experiments. The experiments examined the influence of the method for driving extension (foam base, rubber base, plate base or conveyor base) for localization of deformation in overlying layers of brittle-only and brittle-viscous materials representing the earth’s crust. All experiments were performed at the Tectonic Modelling Laboratory of the University of Bern. Detailed descriptions of the experiments and monitoring techniques can be found in Zwaan et al. (2019) to which these data are supplementary material. All experiments were monitored with top view photographs (SLR camera Nikon D-100 6.1 MPx). The photograph time steps depend on the applied extension velocity, but are generally 1 or 2 min. Six experiments were also monitored with an X-Ray computed tomography technique using a 64 slice Siemens Somatom Definition AS X-ray CT-scanner (Zwaan et al., 2016) with varying time intervals (5-30 min). CT-data was analyzed with the software OsiriX (Pixmeo SARL).

Description: Extensive passive seismic monitoring was carried out between September 2017 and September 2018 over the Los Humeros geothermal field in Mexico. This experiment, in addition to several geophysical, geological, and geochemical surveys was conducted in the framework of the European H2020 and Mexican CONACyT-SENER project GEMex for a better understanding of the structures and behavior of the local geothermal system currently under exploitation, and for investigating future development areas. 25 broadband stations (22 Trillium C-120s and 3 Trillium C-20 PH) recording at 200 Hz, and 20 short period stations (Mark L-4C-3D) recording at 100 Hz comprised the network which is sub-divided into two sub-networks. An inner and denser (~1.6-2 km inter-station distance) pseudo-rhomboidal array (27 stations) was laid out to cover the producing zone and retrieve local seismicity mainly associated to injection and production operations, and to comply with beamforming of ambient noise and time reverse imaging techniques. An outer and sparser (~5 km minimum spacing) array was placed at around 30 km radius surrounding the inner network, and was mainly dedicated to larger scale imaging techniques, such as seismic ambient noise tomography, and regional earthquakes tomography. The GEMex project is supported by the European Union’s Horizon 2020 programme for Research and Innovation under grant agreement No 727550 and the Mexican Energy Sustainability Fund CONACYT-SENER, project 2015-04-68074. Waveform data are available from the GEOFON data centre, under network code 6G, and are embargoed until January 2023.

Description: The Liquiñe-Ofqui fault system (LOFS) in south-central Chile provides a natural laboratory to assess the interplay between magma/hydrothermal fluid flow and crustal deformation. Understanding these processes is of paramount importance for geothermal energy exploration and seismic hazard assessment. We deployed a dense seismic network (Sielfeld et al., 2019) at the northern termination of the LOFS in south-central Chile (~38°S) between 2014 March and 2015 June. The main aim was to better understand the significance and implications of seismic activity in relation to geological information such as the complex fault-fracture network, volcanoes, and the stress field estimated from geological data. As a result, the network was designed to monitor the northern segment of the LOFS on a more regional scale rather than concentration on the activity of one volcano. The network covered a ~200‐km‐long section of the Southern Volcanic Zone, including several Holocene stratovolcanoes (Callaqui, Copahue, Caviahue Caldera, Tolhuaca, Lonquimay, Llaima, Sierra Nevada, Sollipulli, Villarrica, Quetrupillán, Lanín (La), and Mocho‐Choshuenco). Waveform data are open and available under network code 3H from the GEOFON data centre under license CC BY 4.0.

Description: The unrest of el Hierro Islands started in 2011 with a submarine volcanic eruption. In order to better characterize unrest of El Hierro Island 9 landstations were installed on the Island of ElHierro (Figure 1) between March 2015 and June 2016. General remarks about the seismological data: Data quality is very good due to the low anthropogenic noise level of the Island of ElHierro. All stations register with 200 samples per second. H09S has timing problems, most likely related fast-growing scrub and bushes around the station. The broadband stations H01B and H07B only registered the first 30 minutes of every hour. Recorders were EarthData PR 6-24 3-channel Short period seismometers were MARK L-4C-3D Broadband sensors were: Guralp-CMG-3ESP 3-channel seismometer Waveform data is available from the GEOFON data centre, under network code 2L, and is embargoed until January 2021.

Description: Definitive digital values of the Earth's mangetic field recorded during 2014 at INTERMAGNET observatories around the world. Data includes minute, hourly and daily vector values, along with observatory baseline values for quality control. Annual means are also included. All data is included on the single downloadable archive file (gzipped tar format) available from this landing page. This is the 24th annual publication in the series. Some national data institutions may have related DOIs that describe subsets of the data. These DOIs are shown under "Related DOIs to be quoted". For more information on the data formats used in this publication and the technical standards used to create the data, please refer to the INTERMAGNET Technical Manual and the Technical note TN6 "INTERMAGNET Definitive One-second Data Standard".

Description: A line of 6 broadband seismometers have been deployed across a ridge in the Hualien County (Eastern Taiwan). From March 2015 to June 2016 the network has been continuously recording waves incoming from the Taiwanese regional seismicity. During that period, more than 2000 earthquakes with magnitudes Ml〉3 and distant from less than 200km were recorded. The hill is well approximated by a triangular topography of 3600m in length by 900m in height. Waveform data are open and available from the GEOFON data centre, under network code 5K.