This information details the method of calculating dilatancy from pore volumometry measurements. In the velocity step tests, an initial shear enhanced compaction phase drastically reduces the sample pore volume until the sample yields. After the sample yields, the pore volume continues to decrease but at a lower rate of decrease. The imposed velocity steps cause compaction or dilation of the sample material that is superimposed on this overall compaction trend. Pore pressure is held at a set point in all tests and any volume changes in the control system are assumed to correspond directly to changes in pore volume in the sample. The method here is aimed at producing quantitative, reproducible values for dilatancy from experimental data. The script fits a polynomial function to all the volume data to give the overall trend of the shear enhanced compaction. The data position of the start of the velocity step of interest is entered manually into the function. When dilatancy occurs on a step change in velocity, this is not immediately recorded using volumometry, as the permeability of the sample will produce a transient response as pore fluid pressure equilibrates between the sample and the pore fluid pressure system. To discount time effects, every velocity step was processed with a ‘time_dep’ phase for the first 100µm of displacement after the imposed velocity step change. Using the values for ‘vel_step’ and ‘time_dep’, the volumometry data are split into separate matrices incorporating time and volume preceding and following the velocity step change. A linear regression model fits a polynomial curve to the pore volumometry data and returns the coefficient of determination (R2) for the fit of the model. The shear enhanced compaction phase prior to yield is included in the fit. The code incrementally adds the value entered for ‘step’ to the data in the velocity step of interest. This ‘step’ value is a positive or negative value depending on whether the velocity has increased or decreased, respectively. A new linear regression model is then fitted to the whole dataset and if the R2 value has increased, the code will continue to loop to add the value of ‘step’ to the pore volume data. It concludes when the R2 value reaches a peak and begins to decrease, as the fit is no longer improving. We assume at this point that the effect of the dilatancy due to the velocity step has been removed, and the cumulative sum of the ‘step’ values is equivalent to the dilatancy. As the loop goes one iteration past the optimum R2 value, the code reverts to the previous set of values with the best R2 value. In experiments with multiple velocity step changes, the code needs to retain the previous corrections of the data. The function ‘pf_correct’ is used to correct the velocity steps that have been previously processed. The values for ‘vel_step’, ‘time_dep’ and the returned value of ‘offset’ need to be given in the inputs for ‘pf_correct’.

HystLab (Hysteresis Loop analysis box), is MATLAB based software for the advanced processing and analysis of magnetic hysteresis data. Hysteresis loops are one of the most ubiquitous rock magnetic measurements and with the growing need for high resolution analyses of ever larger datasets, there is a need to rapidly, consistently, and accurately process and analyze these data. HystLab is an easy to use graphical interface that is compatible with a wide range of software platforms. The software can read a wide range of data formats and rapidly process the data. It includes functionality to re-center loops, correction for drift, and perform a range of slope saturation corrections.

The dataset describes the results of high pressure experimental measurements of three contrasting 'tight' rocks; a Bowland Shale, a Haynesville shale and Pennant sandstone. The results are tabulated as a csv file, listing experimental parameters, confining pressure, argon gas pore pressure and permeability. complementary measurements of key petrophysical data are provided - bulk modulus of compressibility, porosity TOC and density.

CO2 equation of state software from UKCCSRC project: Tractable equations of state for CO2 mixtures in CCS: Algorithms for automated generation and optimisation, tailored to end-users and tutorial presentation to support equation of state software. Grant number: UKCCSRC-C1-22.

Numerical models of mass flows and tsunamis that they generated with their entrance in the sea. The mass flows propagate in the Sciara del Fuoco of Stromboli. The mass flows are not real events, but are rather used as a sensitivity analysis to examine tsunamigenic potential of mass flows of landslides and pyroclastic flows of different durations, volumes and coherence. The data was generated with the two fluid version of the Volcflow model. The material includes tsunami height measurements (plotted as well as raw data) for each run recorded by virtual gauges located around the island (map of gauges included), maximum wave height data in the area around the island and at the shores, a video visualisation of the mass flows and resultant tsunami waves, and a figure of the final deposit from each modelled mass flow. The numerical model simulations were carried out by Symeon Makris. The Volcflow code is not included in the submitted material but it is open source and can be downloaded here: https://lmv.uca.fr/volcflow/

Supplemental files (plotting scripts, tables) in support of Greene et al., 2019. Early Cenozoic Decoupling of Climate and Carbonate Compensation Depth Trends. PP2019.Greeneetal.SupplementaryTables.xlsx Tables S1-S3. Raw data for CCD time-slice reconstruction for NP8, NP10-11, and NP12-13, respectively. Raw data for each CCD time-slice reconstruction (program, site number, current and paleo- latitude/longitude, underlying basement age, current water depth, sediment cover, reconstructed paleodepth (see Methods), wt% CaCO3 mean, and individual wt% CaCO3 measurements. Tables S4-S18. Model output: time series of pCO2 and Ca2+ weathering flux for each ensemble #1 experiment. Tables S19-S22. Final year model output for each experiment in each ensembles 1-4. Final year model output for each experiment ensembles 1-4 (outgassing rate modification factor relative to x3 pre-industrial pCO2, pCO2, mean ocean [DIC], mean ocean [ALK], POC export, Ca2+ weathering flux, mean ocean temperature, mean land surface air temperature, overturning stream function min/max, and final year CSH and CCD (following Goodwin and Ridgwell [2010]). For ensembles 2-4 bolded columns indicate model variables roughly fixed across all experiments within the ensemble. PP2019_Greeneetal_Subsidence.m Matlab script for calculating paleodepth from current water depth and seafloor age. This script plots a subsidence curve for a single location and displays the paleodepth for that location at a user-specified point in the past using the subsidence equations from Cramer et al., 2009 'Ocean overturning since the Late Cretaceous: Inferences from a new benthic foraminiferal isotope compliation', Paleoceanography 24(4), PA4216 with a simplified sediment unloading term. PP2019_Greeneetal_plot_CCDcontour.m Matlab script for contouring depth/wt% carbonate data into contoured CCD snapshot. In support of Published Paper: Greene, S.E., Ridgwell, A., Kirtland Turner, S., Schmidt, D.N., Pälike, H., Thomas, E., Greene, L.K. and Hoogakker, B.A.A. (2019), Early Cenozoic Decoupling of Climate and Carbonate Compensation Depth Trends. Paleoceanography and Paleoclimatology, 34: 930-945. https://doi.org/10.1029/2019PA003601

These data consist of spatial and temporal datasets for 7 different small-scale laboratory experiments of fluid-driven fractures, described in the paper The hidden internal flow dynamics of shear-thinning magma in dikes (Kavanagh et al., 2025, accepted in AGU Advances, March 2025). These experiments, conducted at the University of Liverpool, are analogue models of magma transport via flux-driven dykes. The 7 experiments are named HEC1, HEC2, HEC3, XG1, XG2, W1, W2. Experiments HEC and XG involved the injection of a shear-thinning fluid (a hydroxyethyl cellulose polymer (HEC) and xanthan gum solution (XG)), whilst experiments W1 and W2 involved Newtonian water injections. Experiments HEC1, HEC3, XG1 and W1 were ‘seeded fluid experiments’ or ‘PIV experiments’, whilst experiments HEC2, XG2 and W2 were ‘seeded gelatine experiments’ or ‘dyke-thickness experiments’. We provide the raw experimental data along with the Matlab scripts used to process and plot the data. Further information is provided in the containing README documents.

Rheometric data for suspensions of bubbles and/or particles in a Newtonian suspending liquid. Experimental suspensions are intended as analogues for multiphase magma and lava (containing volatile bubbles and solid crystals). The dataset comprises data collected using a Anton Paar MCR702 rheometer. A Newtonian working fluid (sugar syrup) suspends variable fraction of gas bubbles and solid particles. Data relating to suspensions of bubbles only are found in the folder Bubble_Suspensions, and data relating to suspensions of bubbles and particles are found in the folder 3Phase_Suspensions. For all suspensions, rheometric data were collected in a wide-gap concentric cylinder geometry using both rotational rheometry (flow curves) and oscillatory rheometry (frequency sweeps). Bubble fraction and particle fraction are varied systematically. In each folder a matlab script is provided which facilitates extraction and analysis of the raw data, which are held in .csv files

These data represent a series of analyses exploring the seismic behaviours of low-cohesion volcanic sediments – in this case the Neapolitan Yellow Tuff - under varying strain rates. The data include deformation logs from triaxial compression experiments, and the accompanying 12-channel acoustic emission recordings at 10 MHz. These are paired with X-Ray Computed Tomography images of one of the cores from both before and after deformation, to examine damage behaviour. These data include: Deformation logs captured from the triaxial press Acoustic emission event data Processed acoustic emission sonograms for selected events Matlab code for processing of sonograms Matlab code for statistical analysis of the acoustic emission data Before and after X-Ray Computed tomography data for a core which underwent 2% strain at a rate of 4x10-6 s-1. These data relate to Rowley et al - Deformation controlled Long-Period seismicity in low cohesion volcanic sediments https://doi.org/10.31223/osf.io/7rkzv

The excel table presents the high-resolution LA-MC-ICP-MS Pb isotope data (± 2S%) produced for the Paleocene interval of Fe-Mn crust core sample 085_004 recovered in 2016 during the JC142 expedition to Tropic Seamount for the MarineE-Tech project. A detailed description of the sample can be found in Josso et al., 2019 and 2020. The Scripts data are Matlab and R readable file presenting all scripts and built-in functions used during the processing of the data to demonstrate the presence of astronomical imprint on the cyclicity of the Pb isotope ratios.