These spatial and temporal image datasets are of mixing fluids in a dyke-like (slot) geometry investigated using low-temperature laboratory experiments which are described in the paper Analogue experiments to investigate magma mixing within dykes (Havard et al., in review in Bulletin of Volcanology, January 2025). The experiments are low-temperature analogue models of magma mixing within dykes and are dynamically scaled to the natural system. Seven unique experiments using seven different miscible fluid pairs, representing two magmas of differing composition, were conducted at the University of Liverpool. The miscible fluid pairs were (and the shortened experiment name is): water and diluted glycerol (ExpW1); water and pure glycerol (ExpW2); water and glycerol mixed with golden syrup (ExpW3); water and pure golden syrup (ExpW4); diluted glycerol and pure glycerol (ExpD1); pure glycerol and diluted golden syrup (ExpG1); pure glycerol and pure golden syrup (ExpG2). The experiment setup consisted of the lower density fluid placed above the higher density fluid in the tank apparatus, then the tank was inverted to initiate the interaction between the fluids. The raw data is a temporal series of photographs (JPGs) collected over the course of each experiment located within the subdirectories “Raw”. The raw data is processed using the scripts accompanying the dataset to show how density and mixing evolves spatially and temporally within the slot. Further details are available in the directory metadata file.

Grant: SCAPE (UKCCSRC Call 2). This is a Matlab code that calculates the change in internal concentration of CO2 in a building as a cloud of CO2 engulfs the building. The CO2 is assumed to enter through any openings in the building. This Matlab code take a series of inputs including wind speed, building geometry, geometry of external temperature and external CO2 concentration (all inputs are listed in the headers of the spreadsheet ‘inputs.xlsx’) and calculates how the internal CO2 concentration and temperature changes over time and the toxic dose of CO2 received by individuals inside the building. Full details will be given in publication https://doi.org/10.1016/j.ijggc.2019.102849. File formats: .xlsx and .m

A number of processes, both natural and anthropogenic, involve the fracture of rocks subjected to tensile stress, including vein growth and mineralization, and the extraction of hydrocarbons through hydraulic fracturing. In each case, the fundamental material property of mode-I fracture toughness must be overcome in order for a tensile fracture to nucleate and propagate. Whilst measuring this parameter is straightforward at ambient pressure, estimating the fracture toughness of rocks at depth, where they experience a confining pressure, is technically challenging. Here, we report a new analysis combining results from standard thick-walled cylinder burst tests with quantitative acoustic emission to estimate the mode-I fracture toughness (KIc) of Nash Point Shale at elevated confining pressure, simulating in-situ conditions to approximately 1km. In the most favorable orientation, the pressure required to fracture the rock shell (injection pressure, Pinj) increases from 6.1 MPa at 2.2 MPa confining pressure (PC), to 34 MPa at 20 MPa confining pressure. When fractures cross the shale bedding, the required injection pressures are 30.3 MPa (at Pc = 4.5 MPa) and 58 MPa (at Pc = 20 MPa), respectively. Applying the model of Abou-Sayed (1978) to estimate initial flaw size, we calculate that this pressure increase equates to an increase in KIc from 0.6 MPa.m1/2 at 3.2 MPa differential pressure (Pinj - PC) to 4.1 MPa.m1/2 at 22 MPa differential pressure. We conclude that the increasing pressure due to depth in the Earth will have a significant influence on fracture toughness, which is also a function of the inherent anisotropy.

The supporting data for C. Harris et al., 2021, 'The impact of heterogeneity on the capillary trapping of CO2 in the Captain Sandstone', International Journal of Greenhouse Gas Control. We supply experimental and numerical simulation data used in the paper. The supplied codes reproduce each figure. The codes are split into 2 folders, descriptions of each of the folders are given below: 0 - README. This contains detailed instructions on using the supplied files. 1 - Main simulations. This contains the code to produce the main CMG (Computer Modelling Group) simulations outlined in the paper, with various input variable files. 2 - Other figures. This contains the code to produce other figures within the paper which do not rely on numerical simulations, including the experimental data.

This work presents a detailed three-dimensional finite element based model for wave propagation, combined with a postprocessing procedure to determine the fracture intensity caused by blasting. The data generated during this project includes output files of all simulations with detailed fields, geometries and meshes. The model incorporates the Johnson-Holmquist-2 constitutive model, which is designed for brittle materials undergoing high strain rates and high pressures and fracturing, and a tensile failure model. Material heterogeneity is introduced into the model through variation of the material properties at the element level, ensuring jumps in strain. The algorithm for the combined Johnson-Holmquist-2 and tensile failure model is presented and is demonstrated to be energy-conserving, with an open-source MATLABTM implementation of the model. A range of sub-scale numerical experiments are performed to validate the modelling and postprocessing procedures, and a range of materials, explosive waves and geometries are considered to demonstrate the model's predictive capability quantitatively and qualitatively for fracture intensity. Fracture intensities on 2D planes and 3D volumes are presented. The mesh dependence of the method is explored, demonstrating that mesh density changes maintain similar results and improve with increasing mesh quality. Damage patterns in simulations are self-organising, forming thin, planar, fracture-like structures that closely match the observed fractures in the experiments. The presented model is an advancement in realism for continuum modelling of blasts as it enables fully three-dimensional wave interaction, handles damage due to both compression and tension, and relies only on measurable material properties. The uploaded data are the specific simulation outputs for four explosion models occurring on two different rock types, and the specific fracture patterns generated.

Supporting data for the paper 'Small-scale capillary heterogeneity linked to rapid plume migration during CO2 storage'. We supply experimental, analytical and numerical simulation data used in the paper. The supplied zipped folders follow the same order as the main paper, with codes to reproduce each figure (and those in the supporting information PDF). There are also video files (in the 5_Field_scale_simulation zipped folder) showing the final CO2 plume evolution from the static images in the main paper Figure 4. Descriptions of each of the folders are given below: 0 - README. This contains detailed instructions on the data and using the supplied files. 1 - Scaling analysis. This contains the scaling analysis analytical methods, with figure generation for Figure 1 in the main paper. 2 - Petrophysics. This contains all the petrophysical experimental data, analysis files and core flood simulation files. This is used to produce Figure 2 in the main paper. 3 - Fine_resolution_simulations. This contains the simulation files, Matlab post processing files and figure generation for the fine resolution simulations, presented in Figure 3 in the main paper. 4 - MIP_upscaling. This contains simulations files, Matlab post processing files and figure generation for the macroscopic invasion percolation scheme. The results of this are presented in the supporting information document. 5 - Field_scale_simulation. This contains the simulations files, Matlab post processing files and figure generation for the final field scale simulations in the main manuscript Figure 4 and in the supporting information. In each folder are separate READMEs containing specific information relevant for the included files.

The mechanics of olivine deformation play a key role in long-term planetary processes, including the response of the lithosphere to tectonic loading or the response of the solid Earth to tidal forces, and in short-term processes, such as post-seismic creep within the upper mantle. Previous studies have emphasized the importance of grain-size effects in the deformation of olivine. Most of our understanding of the role of grain boundaries in the deformation of olivine is inferred from comparison of experiments on single crystals to experiments on polycrystalline samples, as there are no direct studies of the mechanical properties of individual grain boundaries in olivine. In this study, we use high-precision mechanical testing of synthetic forsterite bicrystals with well characterized interfaces to directly observe and quantify the mechanical properties of olivine grain boundaries. We conduct in-situ micropillar compression tests at high-temperature (700°C) on bicrystals containing low-angle (4• tilt about [100] on (014)) and high-angle (60• tilt about [100] on (011)) boundaries. During the in-situ tests, we observe differences in deformation style between the pillars containing the grain boundary and the pillars in the crystal interior. In the pillars containing the grain boundary, the interface is oriented at ∼ 45° to the loading direction to promote shear. In-situ observations and analysis of the mechanical data indicate that pillars containing the grain boundary consistently support elastic loading to higher stresses than the pillars without a grain boundary. Moreover, the pillars without the grain boundary sustain larger plastic strain. Post-deformation microstructural characterization confirms that under the conditions of these deformation experiments, sliding did not occur along the grain boundary. These observations support the hypothesis that grain boundaries are stronger relative to the crystal interior at these conditions. This data set is associated with the pre-print manuscript with the DOI: 10.22541/essoar.167979601.17867144/v1