Data for NERC grant NE/L000660/1. This is the data supporting Fig. 4 of the publication: Ebigbo, A., Lang, P. S., Paluszny, A., and Zimmerman, R. W. (2016). Inclusion-based effective medium models for the permeability of a 3D fractured rock mass. Transport In Porous Media, DOI: 10.1007/s11242-016-0685-z. It contains numerically computed permeabilities for various realisations of fracture networks. There are six different cases (as explained in the paper).

Elastic constants were calculated by using the stress-strain method and density functional theory for crystals of different symmetry. A toolkit was developed to facilitate the input preparation and output processing for elastic calculations using the Vienna Ab Initio Simulation Packge (VASP). The details of the calculation. methods, and the toolkit will be published elsewhere and linked to this deposit. The deposit contains example folders for the monoclinic, orthorhombic, trigonal, tetragonal, hexagonal and cubic lattices. The purpose of this deposit is to provide examples for the toolkit users.

This dataset comprises ECLIPSE input decks for a 3D reservoir simulation of the CO2 plume at the Sleipner CO2 injection site. This whole reservoir model is an attempt to history match the growth of the plume observed on seismic data. A seismic velocity and density model derived from the 3D reservoir simulation is also included, together with a series of Seismic Unix scripts to create a synthetic seismic section through the Sleipner reservoir model, for comparison with released time-lapse seismic data.

Open source modeling code, with which all data were generated: https://github.com/kuangdai/AxiSEM-3D This code was primarily developed within the NERC-funded project, and used for a at least 10 publications over the past two years: [1] Wolf, Long, Leng, Nissen-Meyer. Sensitivity of SK(K)S and ScS phases to heterogeneous anisotropy in the lowermost mantle from global wavefield simulations, 2021. GJI, 228, 366–386, https://doi.org/10.1093/gji/ggab347 [2] Krier, Thorne, Leng, Nissen-Meyer: A compositional component to the Samoa ultralow-velocity zone revealed through 2- and 3-D waveform modeling of SKS and SKKS differential travel-times and amplitudes, Journal of Geophysical Research. doi:10.1029/2021JB021897 [3] Thorne, M. S., Leng, K., Pachhai, S., Rost, S., Wicks, J., & Nissen-Meyer, T. (2021). The most parsimonious ultralow-velocity zone distribution from highly anomalous SPdKS waveforms. Geochemistry, Geophysics, Geosystems, 22, e2020GC009467. https://doi.org/10.1029/2020GC009467 [4] Haindl, Leng, Nissen-Meyer, 2021. A 3D Complexity-Adaptive Approach to Explore Sparsity in Visco-Elastic Wave Propagation, Geophysics, doi.org/10.1190/geo2020-0490.1 [5] Tesoniero, Leng, Long, Nissen-Meyer. Full wave sensitivity of SK(K)S phases to arbitrary anisotropy in the upper and lower mantle, Geophysical Journal International, 222, 412–435, https://doi.org/10.1093/gji/ggaa171 [6] Thorne, M.S.; Pachhai, S.; Leng, K.; Wicks, J.K.; Nissen-Meyer, T, 2020. New Candidate Ultralow-Velocity Zone Locations from Highly Anomalous SPdKS Waveforms. Minerals 2020, 10, 211. [7] Fernando, Leng, Nissen-Meyer, 2020. Oceanic high-frequency global seismic wave propagation with realistic bathymetry, Geophysical Journal International, 222, 1178–1194, https://doi.org/10.1093/gji/ggaa248 [8] Leng, Korenaga, Nissen-Meyer, 2020. Three-dimensional scattering of elastic waves by small-scale heterogeneities in the Earth’s mantle, Geophysical Journal International, 223, 1, 502–525, https://doi.org/10.1093/gji/ggaa331 [9] Szenicer, Leng, Nissen-Meyer, 2020. A complexity-driven framework for waveform tomography with discrete adjoints, Geophysical Journal International, https://doi.org/10.1093/gji/ggaa349 [10] Leng, Nissen-Meyer, van Driel, Hosseini, Al-Attar, 2019. AxiSEM3D: broad-band seismic wavefields in 3-D global earth models with undulating discontinuities, Geophysical J Int., 217, 2125–2146 Each of publications is based on the code mentioned above, and metadata for running the simulations of the papers are given therein, in a reproducible manner.

Data output from the numerical flow modelling in GRL manuscript ""Evidence for the top-down control of lava domes on magma ascent dynamics"", by Marsden, L., Neuberg, J. & Thomas, M., all of University of Leeds. The models were created using the Laminar Flow module in COMSOL Multiphysics v5.4 by L. Marsden. The following files are uploaded: Archive_Reference_Model.txt (Reference flow model: Gas loss function, Initial H2O content = 4.5 wt.% Excess pressure at depth = 10 MPa, Constant corresponding to crystal growth rate = 4e-6 s^-1 ) Archive_High_H2O.txt (Gas loss function, Initial H2O content = 10 wt.% Excess pressure at depth = 10 MPa, Constant corresponding to crystal growth rate = 4e-6 s^-1) Archive_No_Gas_Loss.txt (No gas loss, Initial H2O content = 4.5 wt.% Excess pressure at depth = 10 MPa, Constant corresponding to crystal growth rate = 4e-6 s^-1) Archive_Gamma_Low.txt (Gas loss function, Initial H2O content = 4.5 wt.% Excess pressure at depth = 10 MPa, Constant corresponding to crystal growth rate = 1e-6 s^-1) Archive_Excess_Pressure_0MPa.txt (Gas loss function, Initial H2O content = 4.5 wt.% Excess pressure at depth = 0 MPa, Constant corresponding to crystal growth rate = 4e-6 s^-1) Archive_Excess_Pressure_20MPa.txt (Gas loss function, Initial H2O content = 4.5 wt.% Excess pressure at depth = 20 MPa, Constant corresponding to crystal growth rate = 4e-6 s^-1) The files uploaded include the reference flow model and where a single key parameter has been changed in the flow modelling. We include data where the key parameter is at the upper or lower limit of the values tested. Data are not included where magma ascent is modelled to stall without the extrusion of a lava dome, as a time dependent model is not run in this case. A solution is provided using equilibrium modelling only. The following variables are output, at conduit centre unless specified: Depth (m), Time(s), Ascent velocity (m/s), Bulk Viscosity (Pa s), Crystal Content, Dome height (m), Gas Volume Fraction, Overpressure (Pa), Shear Stress at Conduit Wall (Pa)

The download .rar file contains a groundwater model of the coastal aquifer in Kwale County, Kenya (ModelMuse Text File) produced by Dr Nuria Ferrer and Dr Albert Folch at the Universitat Politècnica de Catalunya. The model can be used to explore future climate and groundwater abstraction scenarios to provide management recommendations. The download does not include proprietary abstraction data from industry project partners, thus running the model provided here will not reproduce published research findings. The file named”np67IH.bhd” are the initial heads file required to run the model.

This dataset contains digital models of three carbonate (Indiana limestone, Estaillades limestone and Edwards Brown Dolomite) and two sandstone (Bentheimer and Bunter) cores with varying degrees of heterogeneity types and length scales. The models define the spatial distribution of porosity and capillary heterogeneity, and can simulate core floods. To build these models, the workflow by Jackson et al. (2018) was followed. It combines observations from coreflood experiments with a numerical optimisation scheme to characterise the capillary heterogeneity within a rock sample. The data presented here contains the experimental X-ray CT dataset used for the paper "History Matching Core Floods in Heterogeneous Sandstone and Carbonate Rocks" by Wenck et al. (Under Review) along with CMG IMEX modelling files. Funding source: EPSRC iCase studentship with BP

Dataset contains 3D synthetic seismic waveforms for axisymmetric global Earth velocity models. The waveforms were calculated using the finite difference approach with the PSVAxi algorithm (Jahnke, G et al., 2008. doi:10.1111/j.1365-246X.2008.03744.x). The Earthmodels are 1D and use PREM parameters except close to the core-mantle boundary (CMB) where 3D ultra-low velocity zones (ULVZs) are added to the PREM background model. ULVZ are thin layers of strongly reduced seismic velocities located at the CMB that have been observed in several regions of the Earth. The dataset models interaction of the seismic wavefield with ULVZ structure with varying elastic parameters (P-wave, S-wave velocity, density), location (location at source or receiver side along the great circle path), ULVZ length, shape (box, Gaussian, trapezoid) and height. Detailed description of the approach and the model space are given in Vanacore et al, (2016). Data format is SAC (Seismic Analysis Code). Vanacore, E.A., Rost, S., Thorne, M.S., 2016. Ultralow-velocity zone geometries resolved by multidimensional waveform modelling. Geophys. J. Int. 206, 659–674. doi:10.1093/gji/ggw114

Experimental results used to parameterise and a test a mathematical model of uranium diffusion and reaction in soil. The exeperiments and model are described in Darmovzalova J., Boghi A., Otten W., Eades, L., Roose T. & Kirk G.J.D. (2019) Uranium diffusion and time-dependent adsorption-desorption in soil: a model and experimental testing of the model. Eur. J. Soil Sci., doi: 10.1111/ejss.12814. The research was funded by NERC, Radioactive Waste Management Ltd and the Environment Agency through the Radioactivity and the Environment (RATE) programme (Grant Ref NE/L000288/1, Long-lived Radionuclides in the Surface Environment (LO-RISE)).

General circulation model (HadCM3) output of the study by Matero et al. (2017) “The 8.2 ka cooling event caused by Laurentide ice saddle collapse. Data has been processed into netCDF4 - timeseries, and includes the following variables at model resolution: ocean temperature, ocean salinity, precipitation, air temperature at 2m height, depth of the oceanic mixed layer, sea ice concentration and meridional overturning circulation strength. The atmosphere component of the model has a horizontal resolution of 2.5° x 3.75° with 19 unevenly spaced vertical layers. The ocean component has a horizontal resolution of 1.25° x 1.25° with 20 unevenly spaced vertical layers. For more information see published paper, https://doi.org/10.1016/j.epsl.2017.06.011