Mineral Carbonation of CO2 in Mafic Plutonic Rocks, II—Laboratory Experiments on Early-Phase Supercritical CO2-Brine-Rock Interactions

The potential for mineral carbonation of CO2 in plutonic mafic rocks is addressed through a set of laboratory experiments on cumulate gabbro and gabbro-diorite specimens from the Sines Massif (Portugal). The experiments were conducted in an autoclave, for a maximum of 64 days, using a CO2 supersaturated brine under pressure and temperature conditions similar to those expected around an injection well during early-phase CO2 injection. Multiple techniques for mineralogical and geochemical characterization were applied anteand post-carbonation experiments. New mineralogical phases (smectite, halite and gypsum), roughness increase and material loss were observed after exposure to the CO2 supersaturated brine. The chemical analysis shows consistent changes in the brine and rock specimens: (i) increases in iron (Fe) and magnesium (Mg) in the aqueous phase and decreases in Fe2O3 and MgO in the specimens; (ii) a decrease in aqueous calcium (Ca) and an increase in CaO in the cumulate gabbro, whereas in the gabbro-diorite aqueous Ca increased and afterwards remained constant, whereas CaO decreased. The geochemical model using the CrunchFlow code was able to reproduce the experimental observations and simulate the chemical behavior for longer times. Overall, the study indicates that the early-stage CO2 injection conditions adopted induce mainly a dissolution phase with mineralogical/textural readjustments on the external area of the samples studied.


Introduction
The International Energy Agency, in its flagship report "Energy Technology Perspectives", has demonstrated that CO 2 capture, utilization and storage (CCUS) is a key technology for CO 2 emissions reduction, essential to achieve the targets set in the Paris Agreement [1].
In general, the parameters that affect the rate of carbonate minerals' precipitation are brine composition, temperature, pressure and, principally, pH [24]. Mineral carbonation is favored over a higher pH-for instance, above 6.5 pH for Ca-Mg carbonates [25] or above 9.0 pH for CaCO 3 [9].
In the CarbFix project, the water in the Hellisheidi carbon injection site has a temperature ranging from 15 to 35 • C and the in situ pH ranges from 8.4 to 9.8 [26]. The injected water (with dissolved CO 2 ) has a temperature of 25 • C and a pH of 3.7 to 4.0 [27]. According to. [28], during the injection phase, the water with CO 2 will create porosity in the near vicinity of the injection by dissolving primary and secondary minerals. Furthermore, away from the injection well, secondary minerals will precipitate due to the reaction with the Ca-Mg-Fe-rich reservoir rocks.
Other authors [14,29,30] observed a decrease in the permeability of the carbonated rock. Inter-and intragranular pores (10 µm pore throats) were filled with magnesite and characterized using µRaman and SEM-TEM images. SEM-TEM observations [14] showed a magnesite layer separated from the olivine by submicron siderite grains, and poorly crystallized phyllosilicates. Iron oxide and amorphous silica are also observed as secondary phases in batch and flow-through experiments [14,22,29]. Phyllosilicates and chalcedony affect both the porosity and permeability, causing a decrease in the carbonation rate due to the creation of an exfoliated passivation layer rich in silica.
Carbonation experiments have been performed by [31] on continental flood basalt (CFB) samples (10.3 MPa and 90 • C) from eastern Washington, and by Schaef and McGrail [32] in basalt samples representing formations from North America, India and Africa.
Both studies noted that calcite was the first mineral to precipitate. For instance, Schaef and McGrail [32] observed small calcite nodules precipitating after 86 days. Post-reacted samples from Deccan basalts [32] displayed larger and opaque precipitates (after 280 days) and more reddish-brown grains containing a large calcite component (66.0-82.0 wt. %) and minor magnesite (9.1-22.0 wt. %) or siderite (3.1-15.0 wt. %) components. The basalt grains representing Southern Africa (Karoo) were insignificantly reactive, displaying very few carbonate nodules with high Mg (31.0 wt. %) and Fe (29.0 wt. %) contents. More recent basalt carbonation experiments have been performed by [33] on recent Auckland basalt (0.3 Ma). After 140-day experiments (100 • C and 5.5 MPa), they observed ankerite and aluminosilicates as secondary precipitation and an increase in both porosity and permeability. This observation contrasts with the major reported carbonation studies [29,30], where the permeability was observed to decrease due to carbonate growth. Theoretical and experimental studies of rock-CO 2 interactions in wet conditions and low temperature (25-90 • C) and first phases (0.5-100 days) indicate the presence of chemical reactions and, consequently, textural-mineralogical-chemical changes [25]. However, according to [25], in order to optimize low-temperature mineral carbonations, an "equilibrium" between initial CO 2 (acid supply), rock to water ratio and temperature needs to be adjusted in order to effectively mineralize CO 2 within a reasonable time scale.

Methodology
To evaluate the potential for mineral carbonation in plutonic mafic rocks, a six-step methodology was followed (Figure 1):

1.
Selection of representative samples for study and definition of conceptual conditions (rock-brine-CO 2 ) to be studied (Figures 1 and 2a).

2.
Mineralogical, textural and chemical characterization of the specimens before exposure to brine and supercritical CO 2 (SC CO 2 ) (Figures 1 and 2b).

3.
Exposure of the specimens to CO 2 supersaturated brine at selected conditions (supercritical CO 2 : 8 MPa and 40 • C) in the autoclave (Figures 1 and 2c Upon conclusion of the laboratory experiments in Step 3, mineralogical, textural and chemical characterization of specimens and brine chemical analysis were conducted. 5.
Geochemical modelling of the mineral carbonation experiments using CrunchFlow. 6.
Interpretation of results and correlation of experimental and modelling data. where the permeability was observed to decrease due to carbonate growth. Theoretical and experimental studies of rock-CO2 interactions in wet conditions and low temperature (25-90 °C) and first phases (0.5-100 days) indicate the presence of chemical reactions and, consequently, texturalmineralogical-chemical changes [25]. However, according to [25], in order to optimize lowtemperature mineral carbonations, an "equilibrium" between initial CO2 (acid supply), rock to water ratio and temperature needs to be adjusted in order to effectively mineralize CO2 within a reasonable time scale.

Methodology
To evaluate the potential for mineral carbonation in plutonic mafic rocks, a six-step methodology was followed (Figure 1): 1. Selection of representative samples for study and definition of conceptual conditions (rockbrine-CO2) to be studied (Figures 1 and Figure 2a). 2. Mineralogical, textural and chemical characterization of the specimens before exposure to brine and supercritical CO2 (SC CO2) (Figures 1 and 2b).     [34] and [35]. (b) Sample preparation for mineralogical and geochemical analyses before and after SC CO2 exposition. (c) Layout of the experimental setup. Reactor system used for the pressurized CO2 injection (modified from. [36]).

Materials
Two different lithologies of igneous rock from the Sines Massif, Portugal were sampled and used in the experiments: (i) Experiment 1 with a cumulate gabbro (CG) from a cliff near Praia do Norte and (ii) Experiment 2 with a gabbro-diorite (GD) sampled at quarry Monte Chãos [11].
The coarse-grained CG sample was cut into 40 subsample cubes of 27 cm 3 each and four of 1 cm 3 . They were divided in four run sets and one reference set (0-day set). For Experiment 1, each set of specimens had seven cubes of 27 cm 3 , two parallelepipeds of 27/2 cm 3 and one of 1 cm 3 (Table 1). The coarse-grained gabbro-diorite sample was cut into 80 subsamples cubes of 27 cm 3 each and four of 1 cm 3 . For Experiment 2, each set of specimens had 15 cubes of 27 cm 3 , two parallelepipeds of 27/2 cm 3 and one of 1 cm 3 (Table 1).

Figure 2.
(a) Conceptual diagram of the reactive zones (Z1, Z2, Z3, Z4 and Z5) around the injection well according to [34] and [35]. (b) Sample preparation for mineralogical and geochemical analyses before and after SC CO 2 exposition. (c) Layout of the experimental setup. Reactor system used for the pressurized CO 2 injection (modified from. [36]).

Materials
Two different lithologies of igneous rock from the Sines Massif, Portugal were sampled and used in the experiments: (i) Experiment 1 with a cumulate gabbro (CG) from a cliff near Praia do Norte and (ii) Experiment 2 with a gabbro-diorite (GD) sampled at quarry Monte Chãos [11].
The coarse-grained CG sample was cut into 40 subsample cubes of 27 cm 3 each and four of 1 cm 3 . They were divided in four run sets and one reference set (0-day set). For Experiment 1, each set of specimens had seven cubes of 27 cm 3 , two parallelepipeds of 27/2 cm 3 and one of 1 cm 3 ( Table 1). The coarse-grained gabbro-diorite sample was cut into 80 subsamples cubes of 27 cm 3 each and four of 1 cm 3 . For Experiment 2, each set of specimens had 15 cubes of 27 cm 3 , two parallelepipeds of 27/2 cm 3 and one of 1 cm 3 (Table 1). The brine used in the experiments is a natural brine sampled from an old borehole in a saline aquifer (see Section 4.2).

Experimental Procedure (Autoclave)
The experimental setup of the autoclave employed in this experiment (Figure 2c) is based on similar systems described by [37] and [38]. Specific initial conditions in the autoclave were considered due to the planned target: sample material (rock-type and representative sample size), geological environment (pressure, temperature and salinity) and technical equipment (materials for chamber, software, pumps, etc.) for the final arrangement of the experimental device and run conditions. The autoclave (Figure 2c) [36,39] has two CO 2 cylinders (standard industrial CO 2 at 4.5 MPa) that are linked to the other elements of the system by steel connectors (diameter: 5 mm). The first CO 2 cylinder is directly connected with the chamber. The second CO 2 cylinder is connected to a piston pump that operates with a flow of 0.01 g/s. In case of gas leakage in the chamber during the experiment, this pump maintains the experimental pressure defined for the test. The inside of the chamber has a capacity of 2 dm 3 . This is coated with polytetrafluoroethylene (PTFE) to protect the material against corrosion. At the bottom of the chamber, a thermostat controls the internal temperature. The calorimeter and pump are linked to the chamber with pressure and temperature sensors and are connected to a computer.
In detail, the experiment consisted of exposure of mafic rocks to CO 2 -supersaturated brine (i.e., SC CO 2 -rich brine) in the autoclave to a pressure (P) of 8 MPa and to a temperature (T) of 40 • C without flow. The P and T conditions were selected to exceed the CO 2 supercritical (SC CO 2 ) point [40,41] and to simulate the conditions of injection and storage of CO 2 [2,42]. These conditions are representative to a depth of approx. 800 m. The selected exposure time (1, 4, 16 and 64 days) was chosen to identify possible changes in the rock (dissolutions and/or precipitation) during the first injection phases.
The experiments began with the immersion of specimen rocks within natural brine in the chamber. Once the rock is introduced and fully immersed in the brine, CO 2 is injected ( Table 1).
The experimental runs comprised: (a) a pressurized CO 2 injection (3 h, from 4.5 MPa and 20 • C conditions to the SC condition); (b) a pressurized stabilization (period of test, no CO 2 flow inside the chamber) and (c) CO 2 pressure release (3 h, from supercritical conditions to ambient conditions). The final volume of CO 2 in chamber to 80 bar and 40 • C was (i) 1184 cm 3 in experiment with CG and (ii) 218 cm 3 in experiment with gabbro-diorite. The times of filling and emptying the chamber with SC CO 2 were the same (3 h, from ambient conditions to supercritical conditions and supercritical conditions to ambient conditions, respectively), following the chamber manufacturer's recommendations. This is the time required to reach the target pressure and temperature values from the initial ambient conditions.

Material Characterization
To obtain a precise characterization of the rock specimens and to evaluate the changes after SC CO 2 exposure, in mineralogy, texture and chemistry, a set of complementary analytical techniques was repeatedly applied.
The X-ray powder and in situ diffractograms were produced using a Bruker AXS-D8 Advance (Bruker Corp., Billerica, MA, USA), with Cu-Kα radiation (λ = 0.1540598 nm), under the following conditions: scanning between 3 • and 75 • (2θ), scanning velocity of 0.05 • 2θ/s, accelerating voltage of 40 kV, and current of 40 mA. In order to evaluate the mineralogical composition of the specimen surface, in situ grazing incidence geometry experiments were conducted, with incidence of 1.5 • and 2θ scanning from 8 to 60 • .
A Hitachi S-3700N SEM (Hitachi High Technologies, Berlin, Germany), coupled with a Bruker XFlash 5010 SDD detector (Bruker Corp, Billerica, MA, USA), was used for the surface sample chemical analysis. The analysis was performed under a low vacuum at 40 Pa, with a current of 20 kV. An infrared spectrometer Bruker Hyperion 3000 equipped with a single-point MCT detector cooled with liquid nitrogen and a 20 × ATR objective with a Ge crystal of 80 µm diameter was used.
An infrared spectrometer Brüker Hyperion 3000 equipped with a single-point MCT detector cooled with liquid nitrogen and a 20 × an attenuated total reflectance (ATR) objective with a Ge crystal of 100 µm diameter was used. The infrared spectra were acquired with a spectral resolution of 4 cm −1 , 32 scans, in the 4000-650 cm −1 region. In order to ensure the representativeness of the data, each cube was analysed in nine different spots, screening the most exposed surface, and each spot was analysed three times.
The whole-rock geochemistry analysis was provided by XRF, which allows for the quantification of major oxides (SiO 2 , TiO 2 , Al 2 O 3 , Na 2 O, K 2 O, CaO, MgO, MnO, FeO, P 2 O 5 ), sulfur and some minor elements (Rb, Sr, Y, Zr, Nb, Th, Cr, Co, Ni, Cu, Zn, Ga, As, Pb, Sn, V, U, Cl). Analyses were performed with an S2 Puma energy-dispersive X-ray spectrometer (Bruker), using a methodology similar to that adopted by [43]. A description of the standard reference materials (SRM) utilized in the calibration method can be found elsewhere [44]. After the determination of loss on ignition (LOI), samples were fused on a Claisse LeNeo heating chamber, using a flux (Li-tetraborate) to prepare fused beads (ratio sample/flux = 1/10). The software utilized for acquisition and data processing was Spectra Elements 2.0, which reported the final oxide/element concentrations and the instrumental statistical error (Stat. error) associated to the measurement.
Linear roughness was measured by a homemade microprofilometer consisting of 55 parallel needles put in contact. The spacing between the needles' center-line is 500 µm for a total measurable length of 3.5 mm. After the profile measure, a photo of the needles' position was taken. Profiles of the samples were detected by drawing a line joining all the extremities of the needles, with the help of Adobe Photoshop CC 2018. The profile line in JPG format was uploaded on online software (WebPlotDigitizer, https://apps.automeris.io/wpd/) that allowed us to identify the X and Y values of the line according to the pixels composing the JPG photo. Pixel values were later converted to microns according to a scale.
The brine analyses were performed at IGME (Madrid, Spain) by ion chromatography (Dionex 600 de Vertex) and ICP-OES (Varian Vista MPX) before and after the experiment on each run. Iron and magnesium brine content determination, from Experiment 1, was performed at HERCULES Laboratory (Évora, Portugal) with an Agilent 8800 ICP Triple Quad (ICP-QQQ), operating with an RF power of 1550 W, RF matching of 1.7 V, a sample depth of 10 mm, carrier gas (Ar) of 1.1 L/min and plasma gas (Ar) of 15 L/min. Prior to the analysis, the equipment was calibrated with a tuning solution from Agilent, and the sensitivity and resolution were optimized and the doubly charged ions (< 1.84%) and oxides (< 1.10%) were minimized.

Code Descriptions and Capabilities
The simulations were performed using CrunchFlow [45,46], a multicomponent reactive flow and transport code for studying fluid rock interactions in porous media. The reactive transport code has many advantages such as a complete equilibrium thermodynamic treatment, flexible kinetic rate law formulations for each mineral depending on the reaction mechanism, and a range of 3D flow capabilities. The reactive transport code [47] numerically solves the mass balance of solutes, as shown in Equation (7): where φ is the updated porosity, C j is the concentration of component j (mol m −3 ), q is the Darcy velocity (m s −1 ), R j is the total reaction rate affecting component j (mol m −3 s −1 ) and D is the combined dispersion-diffusion coefficient (m 2 s −1 ). The carbonation experiments were conducted in closed batch conditions, that is, without flow, so that Equation (7) can be simplified as Equation (8): The reaction rate is only described as function of time. The mineral dissolution to aqueous phases and precipitation of secondary phases are kinetically controlled, thus a reaction rate is given in term of primary species. The reaction kinetic is treated based on the rate law types including the transition state theory (TST), irreversible, monod, dissolution only and precipitation only [46]. In this work, a common theoretical framework provided by the TST rate law [48][49][50] is adopted. The dissolution and precipitation are treated as reversible at equilibrium, by explicitly including a dependence on Gibbs energy or saturation state (Equation (9)): where R m is the rate of precipitation (rate > 0) or dissolution (rate < 0) of mineral m in mol L −1 s −1 , A m is the reactive surface area and k m the kinetic constant of mineral m (in mol.m −2 .s −1 ). The sign of the saturation index Ω = log Q m K eq determines the sign of the reaction rate. Negative means dissolution and positive means precipitation. The kinetic constant at given temperature T (K) is calculated from Equation (10): where k 25 is the kinetic constant at 25 • C, E a is the apparent activation energy (KJ mol −1 ) and R is the gas constant (J mol −1 K −1 ). The change in initial porosity (φ i ) and mineral bulk surface area (A m ) owing to dissolution is computed from Equation (11): and the change due to precipitation is computed from Equation (12):

Input Conditions for Rock and Fluid Composition
The numerical simulations presented in this article only describe the CG experiments; modelling of the gabbro-diorite experiments is ongoing and will be described elsewhere. The initial mineral composition of the rock was derived from Canilho [51] and adjusted with petrography and mineral chemistry, obtained through SEM-EDS [11]. For the numerical simulations, it was assumed that CG is mainly composed of clinopyroxene and olivine, associated with minor amounts of amphibole and calcic plagioclase. Iron oxides and small retrogradation mineral phases were not considered. As in the autoclave experiments, the modelling considered a closed system in which the solid rock takes 20% of the volume and the remaining space is taken up by the brine. As a requirement of the code database, the CG mineral phases are described with the endmembers of the solid solutions, as represented in Table 2. The reactive surface area for each mineral was calculated based on a random variable of the specific surface areas (SSA), defined in the literature as ranging from 10 to 250 cm 2 /g (e.g., applying Brunauer-Emmett-Teller (BET) measurements) for mafic and ultramafic rock [25,[52][53][54][55][56] from Equation (13): The term SA (geometric surface area = 0.67 cm 2 /g) is the ratio between the rock area and the rock volume = 1.68 cm 2 /cm 3 . Because the dissolution reactions depend on the minerals' surface areas, we adjusted the SA with the first term in Equation (13), in order to reproduce the ideal conditions for CO 2 mineralization. In fact, mineral surface area is the most uncertain and complex kinetic parameter and the evolution (especially for multimineral systems) is not quantitatively understood at present [53]. The definition of this parameter requires random values as the measured values using BET overestimates the reaction rates [53,55,57].
The kinetic constants for the reacting minerals were taken from the literature [52,53]. The fit of the model to the experimental data (aqueous calcium, silica, magnesium and iron) concentrations and pH) was performed based on Table 2. The initial concentration of dissolved CO 2 was calculated to be 5.5 × 10 −1 mol kgw −1 according to the Duan and Sun [16] model, with an imposed total pressure of 80 bars. From the thermodynamic and kinetic point of view, 73 aqueous species were considered in the simulations. The equilibrium constants were taken from the EQ3/6 thermodynamic database [58] included in CrunchFlow. The activity coefficients were calculated using the extended Debye-Hückel formulation (b-dot model) [59], with parameters from the same database. The brine composition in the simulations is similar to that of the brine used in the experiments, with pH = 6.85 (before CO 2 addition). The brine is enriched in Ca and Mg (5.12 × 10 −2 and 2.62 × 10 −2 mol.kgw −1 ) and contains important amounts of sulfur and sodium chloride (Table 3).

Petrographic and Chemical Characterisation of Rocks before and after SC CO 2 Exposure
The two lithologies (CG and GD) used in experiments and tested as described above show similar textural and mineralogical results. After interaction of the rock specimens with the SC CO 2 supersaturated brine in the autoclave chamber, an increase in the dissolution of specimens can be perceived; after the one-day run there is a noticeable surface roughness, which increases for the longer runtimes and eventually leads to material fragmentation for 16-day and 64-day runs ( Figure 3). The mean linear roughness (Ra) is higher on CG from 14.80 µm to 21.89 µm for 16 and 64 days, respectively, whereas for GD the Ra increases from 1.01 µm to 1.85 µm for 16 and 64 days, respectively.

Petrographic and Chemical Characterisation of Rocks before and after SC CO2 Exposure
The two lithologies (CG and GD) used in experiments and tested as described above show similar textural and mineralogical results. After interaction of the rock specimens with the SC CO2 supersaturated brine in the autoclave chamber, an increase in the dissolution of specimens can be perceived; after the one-day run there is a noticeable surface roughness, which increases for the longer runtimes and eventually leads to material fragmentation for 16-day and 64-day runs ( Figure  3). The mean linear roughness (Ra) is higher on CG from 14.80 μm to 21.89 μm for 16 and 64 days, respectively, whereas for GD the Ra increases from 1.01 μm to 1.85 μm for 16 and 64 days, respectively. The petrographic analyses, through transmitted light microscopy, did not show significant mineralogical or textural differences after the experiments. It is important to note that even if there were any changes, they should be located on the surface and the thin section production process could have eliminated them. The powder-XRD performed on global fraction before immersion [11] reflects the modal igneous composition given by clinopyroxene, olivine, amphibole, plagioclase and magnetite for CG and plagioclase, clinopyroxene, olivine and mica for the GD (Table 4). Clinochlore has been identified on both rock samples and probably derived from the alteration of the mafic mineralogical phases. After the 1-, 4-, 16-and 64-day runs, the obtained powder diffractograms (Table  4) shows besides the igneous mineralogy the new presence of talc, vermiculite, and halite on both CG and GD in most specimens. Moreover, XRD performed on the surface of CG and GD 64-days specimens, through grazing incidence (Figure 4), additionally reveals the presence of smectite and gypsum, which were not detected in 0-day specimens. The petrographic analyses, through transmitted light microscopy, did not show significant mineralogical or textural differences after the experiments. It is important to note that even if there were any changes, they should be located on the surface and the thin section production process could have eliminated them. The powder-XRD performed on global fraction before immersion [11] reflects the modal igneous composition given by clinopyroxene, olivine, amphibole, plagioclase and magnetite for CG and plagioclase, clinopyroxene, olivine and mica for the GD (Table 4). Clinochlore has been identified on both rock samples and probably derived from the alteration of the mafic mineralogical phases. After the 1-, 4-, 16-and 64-day runs, the obtained powder diffractograms (Table 4) shows besides the igneous mineralogy the new presence of talc, vermiculite, and halite on both CG and GD in most specimens. Moreover, XRD performed on the surface of CG and GD 64-days specimens, through grazing incidence (Figure 4), additionally reveals the presence of smectite and gypsum, which were not detected in 0-day specimens.  As described for the CG [60], the SEM images obtained before and after each experiment on GD (Figure 5a,b) essentially show an increase in ferromagnesian minerals' dissolution, reflected in the higher surface roughness, as reported by other authors [61]. The added phases show different reflectance on backscattered electron (BSE) images and, according to EDS elemental map distribution (Figure 5c,d), correspond to significant enrichment in chlorine (Cl), sodium (Na), sulfur (S) and carbon (C). The precipitation of salts in the form of efflorescence covers a large part of the surface of the specimens and the cavities generated by the dissolution. The carbon enrichment observed on the edge of 64-day specimens (CG and GD) and associated with salts (Figure 5d and Figure 6a) corresponds to organic material, which accumulates due to runoff during the drying process at 40 °C. The origin of this material remains to be identified, but is probably related to the presence of organic material (hydrocarbon) within the original brine.  The carbon enrichment observed on the edge of 64-day specimens (CG and GD) and associated with salts (Figures 5d and 6a) corresponds to organic material, which accumulates due to runoff during the drying process at 40 • C. The origin of this material remains to be identified, but is probably related to the presence of organic material (hydrocarbon) within the original brine.
In fact, the carbon particles observed at the surface (Figure 6a) are not related to the crystallization of carbonates or any other mineral phase [60]. The ATR-FTIR spectra of these carbon-enriched areas present infrared fingerprints pointing to the presence of triglyceride-enriched areas whose origin is still to be determined [58]. In addition, the analysis by ATR-FTIR of CG and GD after 16 days (Figure 6b) of interaction with CO 2 revealed a decrease in the intensity of the absorption band Whole-rock geochemistry data for CG and GD before and after 64 days within a SC CO 2 -brine solution are shown in Table 5. For CG the most significant variations after a six-day run were an increase in CaO (+0.4 wt. %) and a decrease in The carbon enrichment observed on the edge of 64-day specimens (CG and GD) and associated with salts (Figure 5d and Figure 6a) corresponds to organic material, which accumulates due to runoff during the drying process at 40 °C. The origin of this material remains to be identified, but is probably related to the presence of organic material (hydrocarbon) within the original brine. In fact, the carbon particles observed at the surface (Figure 6a) are not related to the crystallization of carbonates or any other mineral phase [60]. The ATR-FTIR spectra of these carbonenriched areas present infrared fingerprints pointing to the presence of triglyceride-enriched areas

Brine Evolution
Before starting the carbonation experiments, the composition of the brine was analysed using a combination of techniques (see Tables 6 and 7, "Pure Brine" column). The brine used in each run of experiments was recovered immediately after the end of the run for a comparative chemical analysis. Table 6. Chemical analysis of pure brine and brine taken from the reaction chambers (one, four, 16 and 64 days) for Experiment 1 with CG ( Figure S1 as supplementary material).  The variation in the brine chemistry with time for both experiments-Experiment 1 with CG and Experiment 2 with GD-is shown in Tables 6 and 7. The pH, measured just after concluding each run, is more acidic for Experiment 1 (with CG). In fact, the initial pH of 6.85 dropped to 4.5, 4.87, 5.1 and 5.5, after one, four, 16 and 64 days, respectively, in the brine-rock-CO 2 system (Table 1) whereas for Experiment 2 (with GD) the pH decreased from 7 (brine without CO 2 ) to 5.81, 5.61, 6.38 and 6.61. In Experiment 2, the number of specimens was doubled (higher rock/brine ratio), but there was not a clear increase in the concentration of ions in solution, as suggested by similar values of conductivity for both experiments.
For Experiment 1 (CG), in general, the behavior of cations in solution was variable, despite the tendency of Ca 2+ decreasing, whereas K + , Mg 2+ , SiO 2 and total Fe increased.
For Experiment 2 (GD), besides Na + , all other analysed cations had lower concentrations in solution when compared to Experiment 1. Furthermore, the analysed elements did not always show the same behavior as described for Experiment 1. After an increase until the four-day run, Ca 2+ tended to stabilize. K + showed an increase, whereas Mg 2+ concentrations tended to stay the same or increase slightly. SiO 2 clearly increased, but the 64-day run saw a decrease in silica.
In both experiments, SO 4 2− had similar behavior: a smooth enrichment, followed by a decrease for the 64-day run; Na + and Cl − had more unpredictable behavior, displaying either enrichment or impoverishment trends along the four runs.

Outlet Solution Composition
The results of the experimental and simulated variations and tendencies in the output concentration over time are shown in Figure 7. The measured and simulated pH in brine supersaturated with CO 2 are presented as a function of the reactive surface area (Figure 7a) and compared to the variation of dissolved CO 2 (Figure 7b). For Experiment 2 (GD), besides Na + , all other analysed cations had lower concentrations in solution when compared to Experiment 1. Furthermore, the analysed elements did not always show the same behavior as described for Experiment 1. After an increase until the four-day run, Ca 2+ tended to stabilize. K + showed an increase, whereas Mg 2+ concentrations tended to stay the same or increase slightly. SiO2 clearly increased, but the 64-day run saw a decrease in silica.
In both experiments, SO4 2-had similar behavior: a smooth enrichment, followed by a decrease for the 64-day run; Na + and Clhad more unpredictable behavior, displaying either enrichment or impoverishment trends along the four runs.

Outlet Solution Composition
The results of the experimental and simulated variations and tendencies in the output concentration over time are shown in Figure 7. The measured and simulated pH in brine supersaturated with CO2 are presented as a function of the reactive surface area (Figure 7a) and compared to the variation of dissolved CO2 (Figure 7b). Initially, the fluid pH at the onset of the modelling was 4.0 and increased rapidly to 4.5, and then gradually to 7.0, in the range of the measured values. Clearly, over the simulated time (0 to 200 days), pH ranged from 4.0 to 5.1 for the low SSA, increasing up to 7 for high SSA (170-220 cm 2 .g −1 ). The interaction between the rock and acidified fluid is a function of the dissolved CO2 and the reactive surface area, since the decrease in HCO3 -is inversely correlated to the increase in pH. Hence, the increase in pH above 5.5 was correlated to the consumption of the dissolved CO2. As a result, at t < 50 days for SSA < 75 cm 2 .g −1 , and t < 10 days, for SSA = 170-220 cm 2 .g −1 , the pH remains lower than 5.0 and the dissolved CO2 concentration remained constant. The constant concentration of HCO3 -at a low pH denotes a dominant dissolution phase, where the concentration of the primary aqueous species (Ca, Mg, Fe, SiO2) increases with time ( Figure 8). However, at moderate and high pH values, CO2 consumption is inferred from the decline in HCO3 -concentration (below 0.48 mol.kgw −1 ), triggering CO2 mineralization into carbonates. Initially, the fluid pH at the onset of the modelling was 4.0 and increased rapidly to 4.5, and then gradually to 7.0, in the range of the measured values. Clearly, over the simulated time (0 to 200 days), pH ranged from 4.0 to 5.1 for the low SSA, increasing up to 7 for high SSA (170-220 cm 2 .g −1 ). The interaction between the rock and acidified fluid is a function of the dissolved CO 2 and the reactive surface area, since the decrease in HCO 3 − is inversely correlated to the increase in pH. Hence, the increase in pH above 5.5 was correlated to the consumption of the dissolved CO 2. As a result, at t < 50 days for SSA < 75 cm 2 .g −1 , and t < 10 days, for SSA = 170-220 cm 2 .g −1 , the pH remains lower than 5.0 and the dissolved CO 2 concentration remained constant. The constant concentration of HCO 3 − at a low pH denotes a dominant dissolution phase, where the concentration of the primary aqueous species (Ca, Mg, Fe, SiO 2 ) increases with time ( Figure 8). However, at moderate and high pH values, CO 2 consumption is inferred from the decline in HCO 3 − concentration (below 0.48 mol.kgw −1 ), triggering CO 2 mineralization into carbonates. Simulated aqueous silica, Ca, Mg and Fe were observed to initially increase and then level off or decrease over increasing time. For instance, a minor increase in Ca concentration with respect to the initial concentration (∆C Ca = 0.14 mol·kgw −1 ) was observed for t < 40 days before a continuous decrease began. The pattern of increase and decrease match perfectly with the evolution of dissolved CO 2 , and are in agreement with the observations in previous studies [14,25,29,55]. The simulated magnesium concentration doubled before levelling off (after 100 days). Moreover, the SiO 2 (aq) concentration stabilized shortly before a sharp and progressive decrease, indicating quantitatively important secondary silicate phases' precipitation. A steady increase in iron concentration was observed; the simulation indicates levelling off occurred later, after 150 days, and was attributed to fayalite dissolution. Aluminum values, initially 1.0 × 10 −4 mol/kg, dropped sharply by five orders of magnitude, following the albite, Ca zeolite and kaolinite precipitation observed in the first simulation.
Because of the assumptions made in the initial mineralogy, the simulation could not reproduce the potassium and SO4 − behaviors.

Minerals' Dissolution and Precipitation
Simulated primary and secondary minerals saturation indices (SI) are shown in Figure 9a,b. All primary minerals (except albite) are seen to dissolve (negative saturation index and rate, Table 3). SI of albite remains positive over time, implying that the modest decrease in Na + (from 0.371 to 3.69 mol.kgw −1 ) is a consequence of albite volume increase. The dissolution rates (Table 8) indicate that olivine, plagioclase and diopsidic minerals are the most reactive phases.  Two mineral assemblages resulting from the CG alteration are predicted by the geochemical model. A Ca and Fe-Mg carbonate assemblage and another composed of Ca-Mg-Fe silicates (and aluminum) such as calcium-rich zeolites, clay minerals and amorphous silica. Ca-zeolites (mesotile and mordenite) and kaolinite are the first secondary minerals to form simultaneously with albite volume increase, followed by the amorphous silica (Figure 9c). The Si-Al-rich assemblage formation results in a sudden decrease in aluminum and silicon from the outlet solution. The carbonates'  Two mineral assemblages resulting from the CG alteration are predicted by the geochemical model. A Ca and Fe-Mg carbonate assemblage and another composed of Ca-Mg-Fe silicates (and aluminum) such as calcium-rich zeolites, clay minerals and amorphous silica. Ca-zeolites (mesotile and mordenite) and kaolinite are the first secondary minerals to form simultaneously with albite volume increase, followed by the amorphous silica (Figure 9c). The Si-Al-rich assemblage formation results in a sudden decrease in aluminum and silicon from the outlet solution. The carbonates' mineral phases correspond to calcite, dolomite, magnesite, siderite and ankerite. Calcite and dolomite appear in the simulations after 16 days, followed by siderite and magnesite (after 40 days). Ankerite was the last carbonate phase to precipitate, concomitant with Fe-Mg-smectite.

Experimental Results
The petrographic study of the two samples, with SEM-EDS, did not reveal any carbonate phase precipitation. However, an increase in textural roughness in the samples was detected with the number of days of exposure into brine SC CO 2 , as also mentioned in [61]. The diffractograms on global powdered fractions for the two experimental specimens are not considerably different between 0 and 64 days. Nevertheless, they reflect trace amounts of talc and vermiculite, as alteration products, after most experimental runs of mafic silicates. Also present is halite. The grazing incidence diffraction of CG and GD specimens, which helps to detect the mineralogy overlaid on the face specimens, additionally indicated the presence of gypsum and smectite ( Table 9). Considering that talc and vermiculite were not detected by this technique within 64 days, those mineral phases should already be present before the experiments, resulting from the replacement of iron-magnesium igneous phases. Their trace amounts may justify the fact that they are not systematically detected in powder XRD. The absence of significant differences in the global mineralogical and chemical compositions of samples from the two experiments allows us to conclude that the reactions between minerals and fluids were not significant. On the other hand, the mineralogical and physical changes are limited to external areas of the CG and GD specimens exposed to CO 2 -rich brine. The changes could be due to local precipitations/dissolutions at the specimen surface and could represent the early effects of influence of the CO 2 -rich brine on the rock. Table 9. Comparative summary of modelled and experimental mineral phases from cumulate gabbro (GXRD: grazing incidence geometry XRD; phases expressed as vol %).

Primary Phases Secondary Phases
Modelling Clinopyroxene (52) Orthopyroxene (1.2) Olivine (10) Plagioclase (11) behavior of Mg 2+ and total Fe, what is generally described is an increase in its concentration. This was also observed with the GD experiments. • The Ca 2+ (aqueous) decrease can be explained by the high content of this ion in the initial brine. Under the experimental conditions, Ca 2+ should be consumed when forming secondary minerals (Ca-Al-Si; calcium aluminum silicates), as predicted by the modelling. The hypothesis that precipitation of gypsum within the autoclave could account for the Ca 2+ decrease should be discarded, since the sulphate ion concentration is not sufficiently high. Modelling results indicate that for gypsum to precipitate the sulphate concentration should be at least double.

•
The 0.4 wt. % increase in CaO in the total rock composition after a 64-day run seems to reflect not only its incorporation in the structure of Al-silicates but also the crystallization of gypsum during the drying process, after concluding the run. Regarding the experiment with GD, the variation of Ca 2+ and Mg 2+ concentrations in brine shows a very small increase, or a tendency to remain constant. The increase is in line with what is described in the literature [25]. Such an increase in solution is reflected in the decreasing concentrations of CaO (0.42 wt. %) and MgO (0.21 wt. %) in total rock for the 64-day run specimens. FeO shows a similar trend, decreasing by 0.3 wt. % in the same specimens.

Numerical Simulation and Correlation with Experimental Results
Observations on fluid chemistry and pH changes over time indicate that, under similar conditions of P CO 2 and temperature, the simulation with SSA = 120 cm 2 /g shows a relative agreement with the experimental data (Figures 7 and 8). The considered SSA is in the range of the specific area of olivine, plagioclase and pyroxene used in [53]. The predicted mobility of calcium is higher than for other components, with a slight increase in concentration for the first 24 days, before being incorporated within the secondary phases. As mentioned above, this increase is not reflected in Experiment 1.
The model also shows that Ca 2+ has been consumed in Ca-zeolites and/or carbonates; for instance, mesolite and mordenite are the first phases to consume calcium at a low pH < 5.5. Moreover, the gradual increase in Mg concentration is justified by the progressive dissolution of diopside and forsterite, accompanied by a moderate precipitation of dolomite and later of magnesite. The results of the newly formed mineralogy (Table 9) indicate the presence of early neoformation phases (e.g., smectite) that can capture calcium within its structure. Nevertheless, no carbonates were observed or detected during the experiment. In CG specimens, most likely, the secondary phases might be local and covered by the considerable amount of salts lining the rock surface.
The simulation also indicates that the ∆(Mg/Si) and ∆(Ca/Si) ratios are 1.45 and 0.35, respectively, and suggests that olivine, diopside and Ca-plagioclase (the main CG minerals) dissolve stoichiometrically [25]. Therefore, the high ratio between silicon in the simulation and the experiment (Si/Si* = 10) testify to the nonstoichiometric dissolution of CG in the experiment. Furthermore, (∆(Mg*/Si*) = 4) is significantly higher and denotes that magnesium dissolves faster than silica (incongruent dissolution). At pH > 5.5, Ca-Mg carbonates (namely, calcite and dolomite) start precipitating and continue up to a neutral/alkaline pH, consistent with the observations from other authors [25,55]. Calcite remains the dominant carbonate phase at a high pH, indicating its stability under alkaline conditions. Iron displayed a similar trend to Mg 2+ , showing a high release rate. Consequently, a tiny fraction of siderite and Fe-Mg-rich smectite was precipitated, concomitant with the dissolution of the primary Fe-rich minerals (fayalite) and followed by ankerite (Ca (Fe, Mg) (CO 3 ) 2 ). The calculated silica and iron concentrations are significantly higher compared to Experiment 1.
Besides the incongruent dissolution observed in the experiments, CG rocks are less reactive than basaltic rocks. As a matter of comparison, Gysi and Stefánsson [25] measured elevated P CO 2 , ion and Si values in glass basalts, in the same range of the simulations (2 mmol/kgw). Note that the brine used in Experiment 1 has a CO 2 concentration 1.5 times higher than that used previously by Gysi and Stefánsson [55], and thus, according to the those authors' results, further secondary phases are expected, since dissolved CO 2 is available during the time of the experiment. The limitation in secondary phases (in CG) is due to the rock reactivity. The experiment's duration is another factor influencing the secondary phases' precipitation.
The total volume of secondary minerals predicted is 0.8 vol %, dominated mostly by zeolites (0.61% in rock volume), carbonates, especially calcite (0.172% in rock volume), and insignificant amount of clays. The small amounts of secondary phases render their identification in SEM very difficult, as reported in previous studies [62]. After 50 days, the ∆(Ca) < 0 corresponds to a stage where zeolites and carbonates competed with a dominant volume of Ca-zeolites. The Mg/Si and Fe/Si ratios increased over time, resulting from the continuous dissolution of clinopyroxene and olivine and precipitating trace amounts of Mg-Fe-rich carbonates, as well as clay minerals.
SEM and XRD observations showed a significant amount of salt crystals lining the samples' surfaces. However, small amounts of Na 2 O were measured in the rock, contrasting with the noticeable quantity of crystals observed in SEM images. The Na and Cl concentrations also showed a cyclic decrease and increase in concentration, not consistent with an unlikely continuous evaporation of the brine during the experiments. We believe that halite and gypsum crystallized after opening the chamber and degassing, or during drying.
As discussed in Section 2, other authors with similar experimental approaches have observed mineral carbonation occurring in time frames comparable to those applied in our experiments. Still, the lack of noticeable precipitation of carbonate minerals during our experiments was not unexpected, since the experiment was designed to look at early-phase dissolution conditions. The proportion of SC CO 2 /brine/rock, was defined to ensure that the brine would always be supersaturated CO 2 during the experiments, as is expected in early-phase injections when dissolution is the dominant process. Such proportions imply that, during the whole duration of the experiments, the pH remained very low in the autoclave, and carbonates could not precipitate. Exposure of the samples to the CO 2 -supersaturated brine for longer periods would probably result in consumption of the CO 2 , pH increase and precipitation of carbonate minerals. Given the time constraints of the project, for the second run of experiments, with the gabbro-diorite sample, the duration of the runs was kept constant, but the CO 2 /brine/rock proportion was changed. The volume of CO 2 was decreased and the mass of rock increased, to see if mineral carbonation would be observed during at least the 64-day runs. The results of this second round of experiments have not yet been fully explored as geochemical modelling is still being done, but they will be the basis for defining the physical conditions and duration of the next set of experiments.

Conclusions
The results obtained reveal, in the first analysis, the complexity of the systems considered: not only the mineralogy observed in the CG and GD but also the composition of the highly saline brine used in the experiments. Thus, a full fit between modelling and experimental results was not always verified.
The preliminary results of the numerical modelling indicate that carbonate precipitation is possible in a batch reactor at a low temperature and pressure (40 • C and 80 kbar). However, the amount of carbonated minerals is very small. In all the tests conducted, Ca-Mg-Fe types of carbonates (calcite, dolomite, magnesite, siderite and ankerite) are predicted, mostly after 64 days (about 0.18% in volume). At this time, carbonates are coupled with the secondary silicates (0.62 vol %) namely zeolites, clays and amorphous silica. Zeolites and amorphous silica were the first predicted phases, forming just after two days.
After 64 days, the experimental data are in line with the modelling as far as the Ca-Al-Si mineral phases. The presence of smectite was verified by grazing XRD. However, calcite, or any other carbonate also predicted by the model in small quantities, was not recorded by any of the analysis techniques (SEM-EDS, XRD, FTIR). Both experiments detected a large amount of salt (gypsum and halite) lining on specimens' surfaces that was not reproduced by the model. Those should have precipitated after opening the chamber/during drying. Some vestigial secondary phases not detected by XRD can be hidden under the salts and not observed by SEM-EDS. Experiment 2 with GD was not modelled, but the experimental results are in line with the observations from Experiment 1: textural features, given enhanced dissolution of ferromagnesian minerals and increased roughness with run time, and secondary mineralogy precipitation (smectite, halite and gypsum).
The measured complex variations in brine composition and total rock compositions reflect the dissolution under SC CO 2 conditions but also the contribution of the high-salinity brine.
The sensitive analysis conducted on the reactive surface area demonstrated that, for the cumulate gabbros, a specific surface area of 120 cm 2 /g is required to fit the geochemical model instead of the 250 cm 2 /g reported for basalt glass. Alternatively, the BET method can be applied. This preliminary model's results will provide interesting benchmarking, and other alternatives are being considered to adequately match the laboratory results.
The experimental and modelling results for 64 days essentially reflect a dissolution mechanism. For higher pH values, and for 64-day runs, the prevailing mechanism changes and significant carbonation occurs. Although in small amounts, the carbonate precipitation would be significant if upscaled to a reservoir scale. Subsequent experiments within the INCARBON project will attempt to replicate the conditions under which mineral precipitation will be the prevailing process, in order to establish a clear picture of the mineral carbonation potential in the target rocks.