Geothermal fields are prone to temperature fluctuations from natural hydrothermal activity, anthropogenic drilling practices, and magmatic intrusions. These fluctuations may elicit a response from the rocks in terms of their mineralogical, physical (i.e., porosity and permeability), and mechanical properties. Hyaloclastites are a highly variable volcaniclastic rock predominantly formed of glass clasts that are produced during nonexplosive quench-induced fragmentation, in both subaqueous and subglacial eruptive environments. They are common in high-latitude geothermal fields as both weak, highly permeable reservoir rocks and compacted impermeable cap rocks. Basaltic glass is altered through interactions with external water into a clay-dominated matrix, termed palagonite, which acts to cement the bulk rock. The abundant, hydrous phyllosilicate minerals within the palagonite can dehydrate at elevated temperatures, potentially resulting in thermal liability of the bulk rock. Using surficial samples collected from Krafla, northeast Iceland, and a range of petrographic, mineralogical, and mechanical analyses, we find that smectite dehydration occurs at temperatures commonly experienced within geothermal fields. Dehydration events at 130, 185, and 600°C result in progressive mass loss and contraction. This evolution results in a positive correlation between treatment temperature, porosity gain, and permeability increase. Gas permeability measured at 1 MPa confining pressure shows a 3-fold increase following thermal treatment at 600°C. Furthermore, strength measurements show that brittle failure is dependent on porosity and therefore the degree of thermal treatment. Following thermal treatment at 600°C, the indirect tensile strength, uniaxial compressive strength, and triaxial compressive strength (at 5 MPa confining pressure) decrease by up to 68% (1.1 MPa), 63% (7.3 MPa), and 25% (7.9 MPa), respectively. These results are compared with hyaloclastite taken from several depths within the Krafla reservoir, through which the palagonite transitions from smectite- to chlorite-dominated. We discuss how temperature-induced changes to the geomechanical properties of hyaloclastite may impact fluid flow in hydrothermal reservoirs and consider the potential implications for hyaloclastite-hosted intrusions. Ultimately, we show that phyllosilicate-bearing rocks are susceptible to temperature fluctuations in geothermal fields.

1. Introduction

Reservoir rocks in geothermal fields are exposed to thermal fluctuations from natural [1] and anthropogenic temperature sources [2]. These fluctuations range from ~250°C of cooling during thermal stimulation practices [3] to heating of up to 400°C during flow testing [4] and up to ~1200°C during basaltic magma intrusions [5]. Importantly, for the sustainability of hydrothermal systems, temperature variations cause volumetric changes that may impart damage [612] and have the potential to trigger mineral reactions, prompting precipitation or breakdown [13] in altered reservoir rocks [14]. These reactions can affect key reservoir rock properties, such as porosity, permeability, and strength [15], which may influence the capacity for fluid circulation [16, 17] and thus dictate energy production potential [18]. Understanding the lithology-specific development of these properties in response to temperature is important for improved fluid flow modelling in geothermal fields.

Hyaloclastite is a rock type prevalent in subaqueous and high-latitude, glaciated regions such as Iceland, where it forms a major reservoir constituent in several geothermal fields [1921]. It is often highly porous and permeable, such that it is frequently targeted for geothermal production, shallow freshwater aquifers, or carbon reinjection and mineralisation [2224], but it is also weak and can collapse to form a mechanically sealed, impermeable caprock [25]. Hyaloclastites are highly variable, altered volcaniclastic breccias that form explicitly by severe quench-induced fragmentation of magma interacting with a large volume of external water or ice [26]. However, the term is often applied ambiguously in the literature to any lava fragmented by interaction with water [27]; consequently, it is commonly identified in a wide range of water-rich environments, such as mid-ocean ridges [28], seamounts [29, 30], submarine volcanic flanks [31, 32], subglacial tuyas [33], tindar ridges [34, 35], nearshore waters [36], and emergent islands [37, 38]. Therefore, hyaloclastite is more appropriately considered a non-genetic term [39] that is more accurately defined by the descriptive lithological criteria, used herein, of quench-fragmented sideromelane (i.e., basaltic glass) supported by a palagonite matrix [40]. This is in contrast to the largely cohesionless perlite, which is produced following the hydration-induced alteration of the dacitic and rhyolitic glasses [41].

The abundance of palagonite in hyaloclastites arises from the inherent metastable nature of basaltic volcanic glass exposed to fluids [42], especially at moderate to high temperatures [43, 44]. Palagonite is considered the first alteration product of mafic glass; it forms initially as an amorphous phase during the complex, concurrent processes of glass hydration and devitrification [45, 46]. The crystallographic character of palagonite is time-dependent, transitioning from a clear amorphous phase displayed in concentric bands of gel-palagonite to a highly variable assortment of crystal habits (e.g., fibrous, lath-like, or granular structure) termed fibro-palagonite. Fibro-palagonite is commonly dominated by clays, namely smectite, and zeolites [47, 48].

Rocks composed of clay and zeolite minerals tend to be highly sensitive to moderate thermal fluctuations on the order of a few hundred degrees Celsius [49]. In particular, high temperatures may trigger devolatilisation reactions that prompt the breakdown of these minerals. This occurrence leads to a loss of material that results in the creation of porosity, thereby affecting the strength and permeability of the rock [4951]. This is also true for incomplete reactions in rocks experiencing short excursions to high temperature [52]. In the case of palagonite undergoing a temperature increase, dehydration can begin below 200°C and result in greater mass loss than in many other clays [53]. This suggests hyaloclastite may be particularly susceptible to thermally induced devolatilisation reactions that are likely to affect geomechanical properties [49].

Here, we investigate the impact of thermal treatment on hyaloclastite, constraining the mineralogical, mechanical, and physical evolution at a range of temperatures up to 600°C.

2. Material and Methods

2.1. Overview

Samples were initially characterised using a suite of thermal analysis equipment in order to select a series of treatment temperatures at which the hyaloclastite would be analysed in further detail. After fully characterising the mineralogical, physical, and mechanical properties of the as-collected material, cores were dwelled at the selected treatment temperature and cooled to room temperature; petrographic observations were made; and the impact on mineralogy, mass, porosity, permeability, and strength was determined.

2.2. Materials

Krafla caldera, located in northeast Iceland (Figure 1), hosts a well-developed geothermal field operated by Landsvirkjun, the national power company of Iceland. The caldera infill is dominated by hyaloclastite and basaltic intrusions to a depth of 1300 m [21], which comprise the reservoir rock hosting hydrothermal fluids harnessed for heating and energy production [54]. The surficial sample block was collected from the southeastern caldera edge (65°N 41.067; -16°W 43.089) in August 2015, where the outcrops protrude from the surface and are isolated from the hydrothermal reservoir, thus undergoing trivial debilitation from exposure to the high-temperature reservoir fluids. However, the youngest and shallowest hyaloclastites at Krafla formed during the last glacial period [55] and have possibly been buried and then subsequently exhumed due to glacial erosion [56]. To complement these 'fresh' samples, a limited suite of subsurface samples was retrieved during the coring of boreholes KH-4 (70 m) and KH-6 (556 m and 732 m) by Landsvirkjun. The in situ sample temperatures were recorded during and after drilling of the wells; the 70 m sample from KH-4 was measured at 41°C during drilling, whilst after leaving borehole KH-6 to thermally equilibrate for one week, the 556 m and 732 m depths were measured at 145 and 125°C, respectively [57]. Note that the sample collected from 556 m depth may have interacted with a basaltic dyke, located approximately 1 m below. The surficial sample was chosen to ensure it was texturally representative of the local geology and similar to the subsurface samples in terms of clast size and abundance. Upon visual inspection, the selected sample was relatively homogenous and lacked large features such as fractures or clasts greater than 1 cm in diameter. For all the mechanical and permeability measurements, heterogeneities were limited to less than 10% of the sample size of 26 mm by 52 mm for cores and 26 mm by 13 mm for discs. Following sample preparation and prior to all testing, samples were oven-dried at 70°C for 4 hours and subsequently stored in a desiccator at room temperature. All samples were prepared, characterised, and tested at the University of Liverpool.

2.3. Simultaneous Thermal Analysis (STA)

Mass loss and heat capacity were measured against temperature concurrently in a Netzsch STA 449 F1 Jupiter using a simultaneous thermal analyser (STA), which combines thermogravimetric analysis and differential scanning calorimetry (TGA-DSC). Samples were cored to 6 mm by 1 mm discs and loaded into a platinum crucible. Prior to each test run, the sample chamber was purged with argon gas and exhausted to vacuum conditions. Each test was initially conducted with an empty crucible to provide a correction baseline, a second time containing a sapphire standard with the same sample dimensions, to correlate the DSC data, and a final time containing the hyaloclastite sample. All temperature profiles were initially heated to 50°C at 2°C/min for a 10-minute isothermal period. For surficial samples, constant-rate temperature profiles, heated to 700°C at 5, 10, and 20°C/min, as well as stepped isothermal profiles, consecutively heated to each treatment temperature at 5°C/min and dwelled for 12 hours, were conducted. For subsurface samples, constant-rate temperature profiles were heated to 700°C at a rate of 10°C/min only. Data was collected from 50°C at a 0.1-minute resolution, with mass accurate to 0.025 μg and heat capacity to ±2%.

2.4. Thermomechanical Analysis (TMA)

Absolute thermal expansion coefficients were recorded using a Netzsch TMA 402 F1 Hyperion. Samples were cored in to cylinders measuring 6 mm in diameter by 5 mm in height. A core sample was loaded between the piston of the TMA, and the furnace was sealed shut around the sample assembly. Prior to each test, the sample chamber was purged with argon gas and exhausted to vacuum conditions. For each test, the temperature and load profiles were completed twice; initially, a ceramic cylinder, matching the sample proportions to within 5%, was used to supply a correction baseline that was removed from the second test run containing the hyaloclastite sample; this allowed an accurate determination of the sample length changes. All temperature profiles were initially heated to 40°C at 2°C/min for a 10-minute isothermal period. Constant-rate temperature profiles, heated to 700°C at 5, 10, and 20°C/min, were completed using a 1 N load. Subsurface samples were heated to 700°C at a constant rate of 5°C/min. Data was collected at a 0.01-minute resolution, with height change accurate to 1.25 nm and load to 0.01 mN.

2.5. Thermal Treatment

Thermal treatment temperatures of 130°C, 185°C, 400°C, and 600°C were selected following the TGA-DSC and TMA measurements. Sample cores, 26 mm in diameter by 52 mm in height, and Brazilian discs, 26 mm in diameter by 13 mm in height, were thermally treated in a Carbolite CWF 1300 box furnace, to be compared against the as-collected samples. Samples were heated at 5°C/min and, once at target temperature, were left to dwell for 12 hours before being cooled to room temperature at a rate of 5°C/min.

2.6. Mineralogical Analysis

The mineralogical and textural changes at each treatment temperature were investigated using a combination of X-ray diffraction (XRD) analysis, optical microscopy, and quantitative evaluation of minerals by scanning electron microscopy (QEMSCAN). Adjacent samples were cored and subjected to the temperature profile for each of the treatment temperatures discussed above in order to prepare samples for mineralogical analysis; parts of these rocks were cut for thin section preparation, and the rest was crushed to powder (see below). Adjacent samples were selected to minimise potential mineralogical differences from within the hyaloclastite, which are petrogenetically heterogeneous.

XRD analysis was completed on the dried, glycolated, and clay-separated samples using a Panalytical X’Pert Pro MPD diffractometer fitted with an X’Celerator detector. Texturally representative rocks were crushed, in distilled water, to a powder <10 μm using an agate McCrone micronising mill, and subsequently dried at 60°C before being further crushed into a loose powder using an agate pestle and mortar. Clay-separated samples were partially crushed in a ceramic pestle and mortar and then prepared with ultrasonication in distilled water, with the equivalent spherical particle size selected for by standard centrifugation methods. Clay-separated samples were then dried at 60°C and recrushed into a light random powder. For glycolated samples, saturation of a random powder was achieved by using ethylene glycol by vapour pressure at 60°C, for 24 hours. Samples were back-loaded into cavity holders as random powders. A copper X-ray tube was used, with a Ni filter to select for Cu K-α radiation. Scans covered the 2 theta range of 4-70°. Data was analysed using the Relative Intensity Ratio (RIR) method within the HighScore Plus® software, alongside reference patterns from the International Centre for Diffraction Data, Powder Diffraction File 2 Release 2008.

Mineral distribution was imaged at 20 μm resolution for each thin section and at 4 μm resolution for selected sites, using Scanning Electron Microscope Energy-Dispersive X-ray Spectroscopy (SEM-EDS) on uncovered, carbon-coated thin sections imbedded with luminescent dye. Using a QEMSCAN developed by FEI, elemental chemistry was mapped and quantified using two Bruker EDS detectors and matched to known compositions of minerals and glasses. Crystallographic features are not recorded, preventing the differentiation of polymorphs and mineral dissociation (in the case where chemical transport is limited). Mineral and glass abundance was quantified at 20 μm resolution by comparing the relative proportions of pixels in each image, normalised against the pore space.

2.7. Porosity Determination

Prior to and following thermal treatment, the skeletal volume of each sample (i.e., the volume of solid rock, including isolated pore space) was measured using a Micromeritics AccuPyc II 1340 Helium Pycnometer, accurate to ±0.1% of the measured volume. The connected porosity (), into which fluids are able to flow, was then determined by the ratio of the connected pore volume (equivalent to the difference between the core volume () and the measured skeletal volume ()) to the core volume (), such that:

2.8. Permeability Measurements

Gas permeability was measured for each sample core, prior to and following thermal treatment, using a Vinci Technologies gas permeameter with nitrogen gas. The 26 mm by 52 mm sample cores were inserted into a compressible Viton jacket and loaded to 1 MPa confining pressure using a manual valve. Gas flow was automatically set to a constant rate through the sample, which increased by a factor of 2 until a differential pressure of >0.5 psi was achieved between the inlet and the outlet. The flow rate values ranged from 33 to 190 cm3/min. The sample permeability () was calculated using Darcy’s law: where is the flow rate; is the cross-sectional surface area of the sample; is the pressure differential, accurate to 1% full scale, measured across the sample; is the viscosity of the liquid (in this case nitrogen); and is the sample length.

For selected cores, steady-state water permeability was also measured in a hydrostatic loading cell from Sanchez Technologies. Within the load cell, the 26 mm by 52 mm sample cores were inserted into an impermeable Viton jacket and loaded to the desired effective pressure (confining pressure − pore pressure). Confining pressure was applied using low-viscosity silicone oil, and pore pressure was applied using demineralised water. Permeability was measured by maintaining 2 MPa of pore pressure at one side of the sample and 1 MPa of pore pressure at the other, such that a constant 1 MPa pressure differential was upheld with an average pore pressure of 1.5 MPa, and the flow rate was measured, accurate to 0.01 ml/min. For each sample, the permeability was measured at 5 MPa confining pressure increments consecutively, up to 30 MPa, to simulate depths up to and in excess of the deepest hyaloclastite units at Krafla. To ensure that permeability was measured under steady-state conditions and that no gas slippage was occurring, the need for the Klinkenberg [58] and Forchheimer [59] corrections was assessed for each sample and found to not be required for any of the samples in both the gas and water permeameters.

Porosity change associated with fracture closure during loading was also recorded in the hydrostatic cell, prior to each permeability measurement [8], using the monitored volume of water that was expelled from the sample. However, absolute porosity variations are not reported as volume change is not monitored during the first loading step; instead, the porosity reduction (in %) is calculated from the initial porosity measured by helium pycnometry. Note that due to potential damage inflicted during loading, samples measured in the hydrostatic load cell were not used for subsequent strength testing.

2.9. Strength Measurements

The uniaxial and triaxial compressive strength tests and indirect tensile strength tests were performed on the as-collected and thermally treated (TT) samples at dry, ambient room temperature conditions.

The uniaxial compressive strength (UCS) of the rocks was determined using an 8800 Instron uniaxial press. Cylindrical rock cores, for which porosity and gas permeability had been measured, were loaded at a constant strain rate of 10-5 s-1 until failure. Load was recorded with a resolution of 0.1 s at 10 Hz.

The compressive triaxial strength (TXL) of cores subjected to various confining pressures was measured using a Sanchez Technologies TRIAX100 press. Here again, cylindrical rock cores for which porosity and permeability had been determined were placed between the pistons and jacketed using an impermeable Viton sleeve. The sample assembly was subjected to confining pressure () by introducing argon gas in the pressure vessel; 5 MPa confining pressure was applied to each sample whilst a constant 1 MPa differential stress was maintained by controlling the axial load. Note that no pore pressure was applied, so that the confining pressure is equivalent to the effective pressure. Upon reaching the confining pressure, the sample was axially deformed at a strain rate of 10-5 s-1 until rupture, denoted by a stress drop in the mechanical data, or until a stress plateau was reached. The confining pressure, pore pressure, axial stress, and sample deformation were recorded at 1 Hz.

Young’s modulus, calculated by dividing stress over strain, was derived using the gradient of the manually defined, elastic linear loading section of each uniaxial and triaxial strength curve.

The indirect tensile strength (UTS) of the samples was measured by employing the Brazilian disc method using a 5969 Instron uniaxial press. Here, cylindrical discs, 26 mm diameter by 13 mm thickness, were radially loaded at a constant deformation rate of 26 μm/s until a stress drop was recorded, associated with failure. Sample deformation, accurate to ±0.1 μm, and load were recorded with a resolution of 0.05 s, and the tensile strength, , was calculated following the ASTM 2008 standard: where is the maximum applied load (in N) and and are the thickness and diameter of the specimen (in m), respectively.

Note that the mechanical data was corrected for compliance of the loading frame of each press in accordance with the ASTM D7070-16 standard procedure.

3. Results

3.1. Simultaneous Thermal Analysis

Thermogravimetric analysis (TGA) reveals progressive, nonlinear mass loss upon heating above 50°C (Figure 2(a)). For the surficial samples (referred to as 0 m), mass loss is initially negatively correlated with the heating rate, whereby mass is lost more rapidly at low heating rates; however, beyond 250°C, there is no correlation and the mass loss at 700°C ranges from 7 to 10%, reflecting sample heterogeneity. In contrast, the subsurface samples (referred to by their sampling depths of 70 m, 556 m, and 732 m) have distinct multistep mass loss shoulders around 100-200°C and 475-600°C; yet, the magnitude of mass loss in each of these temperature ranges varies with depth. Ultimately, the fraction of mass lost is higher in surficial samples than in subsurface samples.

Differential scanning calorimetry measurements associated with the above thermogravimetric analysis reveal further details associated with the mass loss events (Figure 2(b)). The surface samples show double-shouldered endothermic peaks at 130°C and 185°C (Figure 2(b)), consistent with the 100-200°C mass loss event recorded by the TGA (Figure 2(a)). With the increased heating rate, the endothermic peaks are pushed to slightly higher temperatures (up to ~150°C and 220°C at 20°C/min). At a higher temperature, a faint shoulder develops around 300°C and a broad, low-magnitude peak is evident around 650°C; yet, neither corresponds to distinct mass loss events in Figure 2(a). In contrast, the subsurface samples exhibit a less substantial double-shouldered peak at 130°C and 185°C than the surficial samples and faint, broad peaks around 300°C; however, they display strong endothermic peaks at high temperatures: for the 70 m sample, a wide peak develops between 550°C and 660°C; for the sample from 556 m, the data show a peak at 525°C; and for the sample from 732 m, the data show a peak at 560°C.

Based on the simultaneous thermal analysis (Figures 2(a) and 2(b)), four temperatures were selected as thermal treatment targets: 130, 185, 400, and 600°C. Results from the 12-hour isothermal TGA measurements at these temperatures show that mass reaches a new stable value over long time periods, resulting in mass loss of 5%, 5.8%, 7.1%, and 7.5%, respectively (Figure 2(c)).

3.2. Thermomechanical Analysis (TMA)

Thermomechanical analysis was employed to constrain the length changes associated with the mass loss events observed in Figure 2(a). Surface samples initially exhibit limited thermal expansion, followed by minor contraction between 110°C and 290°C (Figure 2(d)). A return to limited thermal expansion ends at 500°C, after which significant contraction occurs. The temperature of the maximum extent of contraction between 110°C and 290°C and the rate of collapse beyond 500°C are positively correlated with the heating rate. Subsurface samples from 70 m depth show moderate expansion up to 475°C, at which point expansion ceases, before rapidly accelerating beyond 640°C. Note that a visual inspection of the two 70 m samples, following the TMA analysis, indicated the occurrence of vesiculation of the glass phase; we do not study this further as no surficial samples underwent such a process. The 556 m sample is entirely dominated by thermal expansion up to 700°C, whilst the 732 m sample expands rapidly up to 185°C and collapses beyond ~575°C.

3.3. Mineralogical Analysis

A range of mineralogical and petrographic techniques were applied to further investigate the crystallographic and textural impact of suspected reactions. XRD analysis (Table 1) shows that the surficial hyaloclastite is primarily composed of the common mafic igneous minerals, anorthite, augite, and quartz, along with secondary minerals chabazite (i.e., zeolite group) and smectite (i.e., clay group). The analysis indicates that, other than small fluctuations in the anorthite and augite proportions, the mineralogical assemblage remains essentially stable up to 130°C. The fraction of smectite identified in the samples decreases by ~6% following treatment to 185°C and is mostly absent beyond 400°C (Table 1). Analysis of the X-ray diffractograms show that the d(060) peak for unaltered surficial material is at 1.535 Å, whist the d(001) peak expands from 15 Å to 16.2 Å upon glycolation and collapses to 10 Å after dwelling at 400°C for 12 hours, typical of saponite (Supplementary Figure 1). XRD analysis of the subsurface hyaloclastite from Lévy et al. [60] highlights that the smectite transitions from saponite in the surficial and shallow hyaloclastites to interlayered saponite-chlorite in the deeper (556 m and 732 m) hyaloclastites.

Petrographic analysis provides a textural counterpart to the quantitative mineralogical XRD. Thin section analysis, using both optical microscopy and scanning electron microscopy, reveals the highly variable nature of hyaloclastite (Figure 3). In particular, we note significant textural heterogeneity in terms of glass distribution and geometry, crystal fraction and size, and pore space distribution (Figure 3). The pore space is shown to be predominantly hosted in the matrix in both the as-collected and TT samples (Figures 3(b) and 3(r)), with additional, isolated porosity hosted within the highly variable vesicular glass clasts (Figure 3(j)). Textural examination of the thin sections for the as-collected material shows a low quantity of intra-mineral microfractures within the phenocrysts and glass (Figure 3(b)). A qualitative inspection of the TT sample thin sections found a comparable number of intra-mineral microfractures. Given the material heterogeneity, we do not quantify this further, finding that material characterisation which examines porosity and permeability before and after TT on the same sample to be a more robust measure of sample evolution with TT. Thin section analysis also offers no evidence of reaction rims around crystals or of sintering textures, such as particle necking during agglutination occurring during TT (Figure 3). Changes in the abundance of constituent phases as a result of TT are revealed by QEMSCAN (Table 2, Figure 3). Changes in crystallinity are not recorded in the QEMSCAN maps, and so the dissociation of smectite is not quantified; however, the textural evolution highlights that the porosity increase is dominated by changes in the matrix, which we associate with a reduction in smectite revealed by XRD (from ~19% to trace content at , Table 1). Loss of smectite with temperature results in an irregular desiccated texture, distinct from a fractured surface due to its smoother edges and globular distribution (Figures 3(b) and 3(r)), causing the matrix to transition from a continuous, pervasive network (Figures 3(c) and 3(g)) to discrete patches with increasing TT temperature, particularly in the 400°C and 600°C samples (Figures 3(o) and 3(s)).

3.4. Thermally Induced Changes in Porosity and Permeability

The evolution of mass during thermal treatment is reported for the 26 mm by 52 mm samples in order to explore the impact of material heterogeneity, discussed in Section 3.3, on mass loss (indicated by thermal analysis in Section 3.1). Sample mass change in the surface cores is positively correlated with the TT (Figure 4(a)). The samples cored from the as-collected material show a large degree of heterogeneity in the material properties despite being sourced from a single block ~40 cm in length. This hyaloclastite heterogeneity corresponds to substantial scatter in the porosity-gas permeability relationship (Figure 4(b)). However, the effect of this scatter can be accounted for and removed from the TT data by normalising the pore volume fraction change against the initial skeletal volume fraction of the corresponding core prior to TT, : where is the post-thermal treatment porosity; similarly, the resultant permeability changes are normalised by dividing the post-thermal treatment sample permeability with the original sample permeability (Figure 4(c)). The normalised plot shows that the permeability and porosity changes of all samples increase with TT (Figure 4(c)).

The water permeability of the hyaloclastite measured in the hydrostatic cell shows a negative correlation with effective pressure (Figure 5(a); Supplementary Table 1) as the permeable pathways are constricted by pore space closure (Figure 5(b)). Again, the samples exhibit a large degree of scatter in permeability for all TTs, which is independent of the mass change induced by TT. The as-collected samples as well as those TT to 130°C, and to an extent 185°C, show a moderate-to-large drop in permeability above 17.5 MPa that is not observed in the 400 and 600°C TT samples. These effective pressure thresholds are also observed in the porosity evolution plot (Figure 5(b)), which shows a transition from low to high densification with effective pressure across these values. The samples that have been TT to higher temperatures exhibit contrasting evolution of permeability and porosity with effective pressure.

In these cases, the permeability of samples does not change significantly with effective pressure (Figure 5(a)), whilst the porosity shuts linearly with effective pressure (Figure 5(b)); the exception to this behaviour is for one sample (thermally treated to 400°C) which shows an increase in permeability between 7.5 and 17.5 MPa (Figure 5(a)); interestingly, this sample was accompanied by almost no porosity decrease between 17.5 and 23.5 MPa (Figure 5(b)).

3.5. Thermally Induced Impacts on Strength

Uniaxial testing was used to constrain the behaviour of samples under compressive loading (Figure 6(a)). The data show typical stress-strain relationships, such as documented in Heap et al. [61], transitioning from an initial concave-upward segment to being quasilinear, followed by a concave-downward segment. In terms of absolute strain, the beginning and end of the quasilinear segments are not correlated with treatment temperature; however, the average strain of the linear portion decreases from 0.0034-0.0037 in the as-collected 130 and 185°C samples to 0.0021-0.0023 in the 400 and 600°C samples. As the 400 and 600°C samples experience more strain prior to failure than the as-collected and lower temperature samples, the shortening of the linear portion also corresponds to a significantly smaller proportion of the total stress-strain curve. The average gradient of the linear segment, which is associated with rock elasticity and is used to calculate Young’s modulus (Figure 7(d)), also decreases with temperature beyond 400°C, from between ~1.3 GPa and 1.56 GPa for the as-collected 130 and 185°C samples, to between ~0.7 GPa and 1.07 GPa for the 400 and 600°C samples. The data show that the as-collected, surficial hyaloclastites exhibit a range of compressive strength from 8 to 10 MPa (Figure 6(a)). Upon thermal treatment to 130 and 185°C, we observe no systematic changes in the sample strength, although the data show wider scatter (Figure 6(a)), whilst samples thermally treated to 400 and 600°C weaken to between 4.2 and 6.8 MPa (Figure 6(a)).

Triaxial testing was used to constrain the behaviour of the hyaloclastite at a confining pressure of 5 MPa (Figure 6(b)), which is representative of ~300 m depth in the geothermal reservoir, assuming nominal rock and water densities of 2500 kg m-3 and 800 kg m-3, respectively [62]. Under such a confinement, the compressive strength of the as-collected hyaloclastite is ~30 MPa, exhibiting a small stress drop upon rupture. Following thermal treatment at 130 and 185°C, hyaloclastite exhibits a moderate strength decrease of up to 3.3 and 7.2 MPa, respectively; yet, the stress-strain curves show systematic reductions in the stress drop upon rupture. Hyaloclastites subjected to higher temperatures of 400 or 600°C reveal contrasting behaviour in which samples tend to yield at lower differential stress, beginning to flow upon strain (Figure 6(b)); the one exception to this is the strongest sample, treated at 400°C, which was texturally comparable to the other TT samples and yet shows mechanical behaviour similar to samples subjected to lower temperatures. Thus, thermal debilitation slightly weakens hyaloclastite and promotes a shift towards a ductile regime (Figure 6(b); Supplementary Table 2).

Brazilian tests were used to quantify the UTS of hyaloclastite (Figure 6(c); Supplementary Table 3). The as-collected hyaloclastite shows a near-linear stress loading curve and a UTS of 1.5 MPa. Upon thermal treatment to 130, 185, and 400°C, we note a systematic weakening trend down to 0.6 MPa. Hyaloclastites subjected to the higher 600°C temperature do not weaken further but rather undergo more compaction prior to failure.

The peak strengths from the UCS, TXL, and UTS curves show a decreasing trend with porosity and TT (Figure 7(a)). Analysis of the UCS (Figure 7(b)), UTS (Figure 7(c)), and Young’s modulus data (Figure 7(d)) shows they follow trends published for igneous rocks. Young’s moduli obtained from the uniaxial (Figure 6(a)) and triaxial (Figure 6(b)) strength tests reveal a negative correlation with the connected porosity of samples, imparted by thermal treatment. Young’s moduli obtained during uniaxial tests are generally lower than the equivalent triaxial results for a given porosity (Figure 7(d)).

The data presented here highlights the substantial impact of thermal treatment on the mineralogical, physical, and mechanical properties of altered reservoir rocks as present in active geothermal systems. Notably, increasing temperature results in lower smectite contents, higher porosity and permeability, and lower compressive and tensile strengths.

4. Interpretation and Discussion

4.1. Application of Laboratory Results to the Field

The investigation on the thermal stability of hyaloclastite present up to 1300 m depth in the shallow geothermal reservoir within the Krafla caldera provides important constraints on the properties of a common reservoir rock. However, due to the physically and mineralogically heterogeneous nature of hyaloclastite, studies at a laboratory scale [8, 48] are challenging to apply to field scales. The variations noted in the mineralogical assemblage of the surficial hyaloclastite studied here (Tables 1 and 2) are likely due to the physically and mineralogically variable source volcaniclastic components, despite being prepared from a single 40 cm block. Influence of this heterogeneity continuously transpires in the physical and mechanical characterisation of both the as-collected and TT samples, particularly highlighted by the contrasting permeability behaviour in Figure 5(a). Heterogeneity at the field scale may be markedly greater due to sorting and componentry [35] and the degree of alteration, particularly between hyaloclastites exposed to different fluid chemistry [19]. Therefore, the data and trends of thermal treatment are most applicable to smectite-bearing hyaloclastite exposed to meteoric-sourced fluids, such as at Krafla. Here, we interpret the influence of thermal treatment on the physical and mechanical properties of surficial hyaloclastite and compare the resultant properties with those obtained from cores sampled from depth.

4.2. Impact of Temperature on the Mineralogical Assemblage of Hyaloclastite

Thermally treating hyaloclastite results in mass loss (Figure 2(a)), which can be used as a proxy for water loss in hydrous clays [53, 64, 65]. However, it is important to note that phyllosilicates, such as smectite, comprise a group of highly variable minerals that have distinct dehydration and dehydroxylation temperatures [66]. The reduction in the d(001) spacing in the XRD data (Supplementary Figure 1) is linked to the dissociation of the palagonite matrix resulting from smectite dehydration (Table 1). In detail, gylcolation firstly causes the d(001) peak to increase to 16.2 Å, as ethylene glycol replaces water in the interlayer space [67]; subsequent heating to 400°C for 12 hours drives ethylene glycol and water out of the interlayer space, decreasing the d(001) peak to 10 Å [68]. Therefore, the reduction in smectite quantity with temperature is more accurately a case of progressive alteration (dehydration), which ultimately damages the crystal structure and reduces the effectiveness of the diffraction process. This produces smaller peaks in the diffractogram, which are interpreted as a lower mineral quantity; a similar reduction in smectite is not observed in the QEMSCAN data (Table 2) as only chemistry is monitored, not crystal structure. The low-temperature (<200°C) devolatilisation reactions identified in the TGA-DSC data (Figures 2(a) and 2(b)) correlate well with the clay-separated XRD results (Supplementary Figure 1). The endothermic DSC peaks at 130°C and 185°C (Figure 2(b)) are encompassed within the 100–200°C dehydration window expected for saponite [69] and are likely due to the loss of sorbed water and interlayer water [68]. Saponite commonly precipitates from hydrothermal fluids found at depth [70, 71] and has been previously identified in cores from the Krafla geothermal reservoir [19]. The subsurface samples have less prominent saponite dehydration peaks due to a combination of lower saponite-chlorite ratios and the higher in situ temperatures potentially causing partial dehydration prior to measurement [60]. The ability of smectite to rehydrate is dependent upon temperature, cation size, and saturation state [72]. Under many conditions, saponite is able to fully or partially rehydrate [73], suggesting the impacts discussed here may be partially reversible.

The 525–560°C endothermic peaks identified in the 556 m and 732 m samples (Figure 2(b)) are indicative of chlorite dehydroxylation from interstratified saponite-chlorite [68], resulting in a second higher temperature mass loss response (Figure 2(a)). The endothermic peak and associated mass loss are greater within the 732 m sample, suggesting it has a lower saponite-chlorite ratio than the 556 m sample. The combined saponite-chlorite content at Krafla increases erratically from 19% at the surface to ~50% at depth, not including glass content [60], suggesting thermally induced devolatilisation reactions at depth may be more impactful. However, note that all of the TGA-DSC measurements were performed under near-atmospheric pressure conditions and that increased pressure from the subsurface may suppress dehydration and dehydroxylation to higher temperatures [74], such that the temperatures for reactions presented here are potentially conservative compared to those occurring in situ. Following the devolatilisation of the phyllosilicates, the bulk rock density increases due to the residual minerals emitting light elements; upon compaction, this may impact the physical and mechanical rock properties.

4.3. Impact of Temperature on the Thermal Expansivity of Hyaloclastite

In addition to mass loss, the thermal treatment of hyaloclastite also results in direct volume change as the palagonite phyllosilicate phases are prone to swelling at low temperatures and subsequently collapse following thermally induced dehydration and dehydroxylation (e.g., >500°C; Figure 2(d)). However, identifying the swelling properties of individual components in a porous polymineralic rock is challenging due to the potential for overlapping responses and the possibility of accommodating expansion into the pore space [75]. In the surface samples, the impact of saponite dehydration in the thermomechanical data (Figure 2(d)) is subtle, causing a minor compaction of the bulk rock, followed by significant compaction correlating with the chlorite dehydroxylation temperature. This suggests that hyaloclastite retains structural integrity throughout saponite dehydration whilst a higher temperature reaction correlated to chlorite dehydroxylation causes a partial loss of cohesion between the palagonite matrix and the supported clasts. The similar, more distinct collapse in the relatively chlorite-rich 732 m sample at 575°C correlates well with the endothermic peak at 560°C and the mass loss identified between 515°C and 615°C (Figures 2(a) and 2(b)). However, the 556 m subsurface sample has lower saponite-chlorite ratios and does not exhibit this collapse; instead, the response shows limited thermal expansion throughout the heating profile. Thus, we find that hyaloclastites with thermally stable mineralogical assemblage expand with temperature, whereas unstable assemblages may breakdown and cause contraction, as exemplified by the thermal treatment of surficial hyaloclastites (Figure 2(d)). The thermal expansion responses measured under near-atmospheric pressure conditions, and any damage attributed to them, may be pushed to higher temperatures under the increased pressure conditions associated with burial, as dehydration is suppressed [74]. However, without a change in temperature, an increase in pressure alone does not result in smectite dehydration [76].

The occurrence of substantial acceleration in thermal expansion above 640°C from the 70 m samples (Figure 2(d)) coincided with vesiculation of glass clasts in the experimental products; this was not reproduced in samples from other depths, despite the abundance of glass in all samples (Figure 3). In detail, thermal expansivity started to accelerate around 600°C, following a moderate endothermic peak (Figure 2(b)), associated with crossing of the glass transition for hydrated basaltic glass [77, 78], expected in shallow, water-rich environments [43]. Beyond this temperature, structural relaxation of the stress enabled water exsolution and increased vesiculation that sustained the heat absorption leading to a prolonged endothermic plateau (Figure 2(d)), as the samples rapidly expanded (Figure 2(b)).

4.4. Impact of Temperature on the Porosity and Permeability of Hyaloclastite

The mass loss associated with the dehydration/dehydroxylation of the phyllosilicate phases during thermal treatment of hyaloclastite is also responsible for the increase in porosity and permeability displayed in Figure 4. A temperature increase can strongly impact the physical attributes for fluid storage capacity and flow, as the basal spacing of palagonite decreases upon devolatilisation, causing densification and generating intrinsically connected micropores [79]. In detail, the amount and ratio of saponite-chlorite available for reaction provides a first-order control on the maximum extent of porosity possibly created by excursion to high temperature. As such, the time-dependent crystallisation of sideromelane to palagonite [47] suggests that older hyaloclastite units, which can contain more phyllosilicates, are more susceptible to temperature increases and the resultant enhancements to the porous permeable network. Thus, the lithology-specific initial glass chemistry, hydrothermal fluid chemistry, and pressure-temperature conditions will therefore impact the potential for the temperature-induced mineralogical and physical alterations that regulate the evolution of the storage capacity and permeability of a reservoir subjected to thermal fluctuations. In addition, the reservoir fluid properties will impact the hydration state of the phyllosilicate minerals; the increased scatter from water permeability (Figure 5(a)) to gas permeability (Figure 4(b)) is likely a result of interaction between the pore fluid and the minerals (especially clays) lining the permeable network.

At the sample scale, the lack of additional fractures in the TT samples suggests that dehydration/dehydroxylation is the primary mechanism responsible for porosity and permeability gain. Whilst cracking due to thermal stresses during heating and cooling has been shown to impart changes in the physical properties of some rocks, we advance that thermally generated cracks may not necessarily accumulate in hyaloclastite. Thermal cracking has been ascribed to the generation of stresses resulting from (1) a mismatch in thermal expansion or contraction between minerals, (2) thermal expansion anisotropy within a mineral, and (3) thermal gradients across a sample [80, 81]. Hyaloclastites are fragmental rocks with variable degrees of cohesion arising from the presence of palagonite, so thermal stressing is expected to cause minimal damage in hyaloclastite as mineral expansion is in part buffered by palagonite, which is weaker than typical volcanic minerals, and due to the abundant and ubiquitous pore space that accommodates expansion without stress development. Eggertsson et al. [8] showed that thermal stimulation may avoid fracture genesis and result in negligible permeability change, independent of the cooling rate, if the thermal stress associated with expansion is alleviated by preferential closure of the existing porous network, such as present in these samples (Figure 2(b)). We also note no evidence of intra-mineral thermal cracking induced by anisotropic expansion. However, it is possible that inter-phenocryst/glass fragment microfracturing occurred, as noted in hyaloclastite quenched from 350°C [82], but that textural evidence was subsequently overprinted by phyllosilicate dissolution. Thermal microfractures are often non-pervasive in nature [6, 49], suggesting the limited contribution they provide to permeability is unable to account for the permeability evolution revealed in this study. Deformation-induced macrofractures may induce a more substantial permeability increase [83, 84], yet such damage is not imparted by thermal treatment alone. However, the impact of thermal fracturing on the broad range of volcanic lithologies is yet to be well constrained, as discussed in Heap et al. [61], and very few studies focus on highly altered material and the impact of pressure on thermal stimulation [52].

A further measure of changing material properties as a response to TT can be the sensitivity of materials to effective pressure; water permeability shows a negative correlation with effective pressure as porosity is isolated (Figure 5). The as-collected samples as well as those TT to 130°C, and to 185°C, show a significant drop in porosity and permeability during confinement to effective pressures above 17.5 MPa that is not observed in the 400 and 600°C TT samples. The samples that have been TT to higher temperatures exhibit contrasting evolution of permeability and porosity with effective pressure, where porosity decreases linearly and more significantly than the lower TT samples, but permeability is not as sensitive to increasing effective pressure, suggesting that the porous network remains highly connected even as effective pressure is increased. There is one exception (thermally treated to 400°C) which shows an increase in permeability with increasing effective pressure, and a less significant reduction in porosity; this could be a result of irreversible compaction in the sample whereby connected pore space was generated, yet we observed no textural evidence for this.

4.5. Impact of Temperature on the Mechanical Properties of Hyaloclastite

The addition of thermal stress can also impact the resultant mechanical properties of hyaloclastite. In particular, a pressure increase causes a relatively abrupt reduction in pore space and permeability within the as-collected hyaloclastite and the samples thermally treated to 130°C (Figure 5). This transition from elastic to inelastic compaction upon loading is termed [85]. Beyond , compaction and grain crushing lead to a loss of pore space available for fluid flow [8, 50, 86]. Figure 5(b) indicates that shifts to lower effective pressure with thermal treatment. The samples thermally treated up to 185°C exhibit poorly defined changes in permeability associated with , and samples treated to higher temperatures develop no distinct changes with . This suggests that thermal treatment promotes a style of compaction that has limited impact on permeability, despite higher rates of porosity loss (Figure 5(b)).

The devolatilisation of the palagonite matrix results in a weakening of the rock, both in tension and in compression, where the UCS is 6–10 times greater than the UTS (Figure 7(a)). The UTS of hyaloclastite decreased with TT and, following treatment at 600°C, lost coherence and underwent more deformation before rupture (Figure 6(c)). The mode of deformation evolves similarly in compression, as noted by the onset of dilation occurring at lower stress and the accommodation of more substantial strain hardening prior to failure (Figure 6(a)), concordant with the lower Young’s modulus upon increasing treatment temperature (Figure 7(d)). Deformation under an effective pressure of 5 MPa accentuated the distinction between the low TT (≤185°C) and high TT (≥400°C) hyaloclastites as they macroscopically behaved in a brittle and ductile manner, respectively (Figure 6(b)).

The resultant strengths measured in all test types display a porosity control, irrespective of the stress field experienced (Figure 7(a)), which follows the common porosity-strength trend (Figures 7(b) and 7(c)) for a range of igneous rock types, regardless of TT [7, 11, 63, 87, 88]. The treatment temperature has a strong impact on porosity (Figure 4(c)), thereby further influencing mechanical compaction; however, comparison with strength and Young’s modulus data collected from subsurface samples [63] shows that temperature alone cannot explain the mechanical changes occurring within the reservoir, instilling the roles of compaction and alteration on strength changes. In part, this is due to the opposing influences of temperature and compaction on porosity; however, in concert, increased temperature will create additional pore space that may enable more complete compaction and densification to occur at depth [63]. The change in style, from a dominantly brittle to ductile failure, may be explained by a change in geometry of the pore space (e.g., Figures 3(b) and 3(r)), as dehydration increases the connectivity and irregularity of the desiccated pore network [89]. The influence of porosity in controlling sample strength can be further assessed using micromechanical modelling such as the pore-emanated crack model developed by Sammis and Ashby [90] and analytically modified by Zhu et al. [91] to derive the UCS, :

The strength is dependent upon the fracture toughness (or critical stress intensity factor), , the porosity (), and pore size (). This model was demonstrated to successfully approximate the UCS of limestone [91] and porous glass sintered in the laboratory [92] but show arguable efficiency to approximate the strength of heterogeneous, coherent volcanic rocks, owing to their common abundance of microfractures [7, 61]. Here, the UCS values decrease as porosity increases with TT (Figure 8), suggesting that would decrease from ~5 MPa down to ~3 MPa. Assuming remains constant within a single lithology (even upon smectite hydration), the analysis would suggest that the reduction in strength and may be the result of pore creation and widening, as suggested by Heap et al. [61] when evaluating the mechanisms underlying thermally induced rock weakening. However, the applicability of the pore-emanated crack model to hyaloclastite is questionable as the pore space is highly irregular and rock failure is promoted by increasingly more pervasive fracture architecture in thermally treated hyaloclastites.

4.6. Implications for Hyaloclastite-Bearing Geothermal Reservoirs

The temperatures investigated here cover the range of conditions hyaloclastite generally encounters in the shallow geothermal reservoir at Krafla [93] and encompass several important mineral transitions within the palagonite matrix. The transition from smectite to chlorite is common in hyaloclastite-bearing geothermal regions and is predominately controlled by temperature, assuming similar fluid chemistry [19]. Saponite is dominant below 200°C, transitions to saponite-chlorite interlayers between 200 and 240°C, and is fully replaced with chlorite beyond 240°C [19]. Therefore, these results provide a first-order approximation of the changes occurring in hyaloclastite during burial and regional heating. Well and core logging at Krafla suggests that the deepest hyaloclastite units at 1300 m naturally experience temperatures of approximately 320°C [94] and thus would have been subjected to the full range of transitions. However, as the temperature profile is spatially variable across the field, the depth of mineral transitions range over several hundred meters [94].

Some geothermal regions, such as Krafla, are also host to magmatic intrusions that can transfer heat to their host rock, inducing much greater localised temperature changes [5]. Similarly, harnessing of very hot fluids, as experienced during flow testing of a 2100 m deep, near-magma borehole during IDDP-1 [21], resulted in fluids up to 550°C reaching the wellhead. Thus, a shallow magmatic intrusion or superheated/supercritical fluid ingress during flow testing can alter the mineralogical, physical, and mechanical properties of shallow hyaloclastite, common in conventional geothermal reservoirs. This may promote more efficient fluid flow in geothermal fields; yet, fluid flow in clay-bearing rocks is partially controlled by the chemistry of the fluids which can interact with the host [95]; thus, the impact of fluid chemistry on fluid flow is expected to decrease upon palagonite dissociation. Here, comparing the surficial hyaloclastite with the hyaloclastites sampled from the reservoir, we find that the thermally induced mass loss and mineralogical distribution of the deep hyaloclastites cannot be explained by simple heating events, as those performed in our experiments, but rather evoke the importance of pressure-temperature history and chemistry in the hydrothermal system. For instance, the onset of saponite dehydration may increase by several hundred degrees by increasing pressure to 30 MPa [74], such as through rock burial or glaciation, suggesting the physical and mechanical changes detailed here may be shifted to higher temperatures for buried samples. Additionally, the impact of pore fluid pressure may suppress dehydration reactions to higher temperatures than recorded in the dry measurements presented here [96]. As such, the interplay between temperature and pressure produces a complex relationship with the resultant mineralogical, physical, and mechanical properties. Geothermal regimes with low effective pressure and high temperature, such as the shallow subsurface, may exhibit a positive correlation between the temperature experienced and the resultant porosity and permeability, whilst in a deeper high-pressure environment, the breakdown from thermal fluctuations, if any, may promote compaction.

Anthropogenic-induced compaction may possibly occur due to the extraction of geothermal fluids inducing a reduction in the reservoir pore pressure [97]; this would increase the effective pressure and subsequently raise the depth threshold of compaction [98], locally altering the dominant permeability regime. Porosity change associated with temperature fluctuations may also influence the depth threshold of compaction. Similarly, a temperature reduction caused by fluid extraction may result in rehydration and swelling due to the potentially reversible nature of smectite dehydration, thereby clogging fluid pathways.

4.7. Implications for Magma Intrusions in Hyaloclastite

The intrusion of magma into hyaloclastite may result in complex intrusion-host rock interaction, evolution, and feedback, owing to variable degrees of devolatilisation and compaction. Initially, an intrusion would bake the margin, promoting phyllosilicate devolatilisation that improves the porous permeable network, allowing for efficient fluid flow; however, this may be accompanied by changes in strength and a transition from brittle to ductile deformation, which may instead favour compaction of the hyaloclastite along the magma boundary, causing a decrease in permeability [99] and compartmentalisation of fluid flow [100]. The construction versus destruction of a permeable porous network may have crucial impacts for the evolution of shallow magmatic systems. Firstly, the liberation of fluids from smectite and chlorite may, if trapped, locally increase pore pressure and generate induced seismicity akin to hydraulic fracturing [101] and, in extreme cases, cause phreatic explosions, as seen at Viti crater, Krafla [102]. Secondly, the volatiles liberated may promote magma hydration and increase the likelihood of explosive activity [103] or result in cooling and quenching of the magma. Thirdly, and similarly, margin compaction and shutting of the permeable network may limit magma outgassing, which again increases the likelihood of explosive volcanism. And finally, magma intrusion may result in viscous relaxation of sideromelane or melting of hyaloclastite, thus generating new magma, which may mix and interact with the intrusion. Thus, the evolution and feedback between magma and hyaloclastite are likely to be very complex depending on the nature (e.g., chemistry, temperature) and size of an intrusion and on the state and properties of hyaloclastites. We suggest that further research should be undertaken to better understand these complex feedback systems and the implications for the volcanic and geothermal processes.

5. Conclusion

In this study, we show that rock heterogeneity impacts the mineralogical, physical, and mechanical properties of hyaloclastite, which are susceptible to thermal fluctuations experienced within geothermal fields. This is due to the devolatisation of the palagonite matrix, composed of varying phyllosilicate minerals which are sourced from the tendency of metastable basaltic glass to crystallise and produce a smectite-dominant (i.e., saponite) palagonite matrix in cool settings (<200°C) and a chlorite-dominant matrix in hotter settings.

Thermal treatment of surficial hyaloclastite indicates that smectite dehydrates upon heating, with increased treatment temperature causing up to 10 wt.% mass loss at 700°C. This dehydration, quantified at 130, 185, 400, and 600°C, results in an enhanced porous permeable network, a decreased compressive and tensile strength, a decreased Young’s modulus, and a shift from brittle to ductile mode of deformation upon thermal treatment (≥400°C), even at a moderate effective pressure of 5 MPa (as experienced at shallow depths in a geothermal system). We assess the impact of temperature on hyaloclastite sampled from boreholes at 70, 556, and 732 m depth in the reservoir, finding that the glass in the shallow rock may have been hydrated whereas the deeper rocks contain abundant chlorite, which dehydroxylates at ~560°C.

The increased temperatures and pressures experienced in geothermal fields will reduce the hydration state of phyllosilicate minerals, causing weakening and potentially lowering the depth threshold for compaction. Such deformation in geothermal systems and areas surrounding magmatic bodies may result in the construction/destruction of fluid pathways and the compartmentalisation of reservoirs, possibly impacting the progression of volcanic activity.

Data Availability

The data used to support the findings of this study are included within the article and supplementary information.

Conflicts of Interest

The authors declare that they have no conflicts of interest


This study has been conducted as part of a PhD supported by the Natural Environment Research Council (NERC) EAO Doctoral Training Partnership that is funded by NERC (NE/L002469/1), whose support is gratefully acknowledged. This research was further supported by Landsvirkjun National Power Company of Iceland and by Landsvirkjun’s Energy Research fund, in addition to a Starting Grant from the European Research Council (ERC) on Strain Localisation in Magma (SLiM, no. 306488) and a Research Fellowship of the Leverhulme Trust (RF 2019-526\4) to Y. Lavallée and an Early Career Fellowship of the Leverhulme Trust (ECR 2016-325) granted to J.E. Kendrick.

Supplementary Materials

Supplementary Tables 1-3: physical and mechanical data for the surficial samples, including density, porosity, permeability, strength, and Young’s modulus change associated with thermal treatment. Supplementary Figure 1: X-ray diffractograms of the as-collected, glycolated, and unglycolated surficial samples. Supplementary Figures 2-6: QEMSCAN mineral distribution images, at 20 μm resolution, for the as-collected and thermally treated surficial samples. Supplementary Figure 7: annotated examples of the petrographic images presented in figure 3 of the main text. Figure S8: additional tensile strength curves for the as-collected and TT surficial materials. (Supplementary Materials)