Mineral Carbonation Potential (MCP) of Mine Waste Material: Derivation of an MCP Parameter

: The heterogenous mineralogy of ultrama ﬁ c deposits hosting mining operations makes it challenging to accurately determine the waste rock’s mineral carbonation potential (MCP). Additionally, the signi ﬁ cantly higher carbonation capabilities of olivine than serpentine add to the di ﬃ - culty. To address this issue, in this work, a new and unique tool called the MCP calculator was developed as a Microsoft ExcelTM spreadsheet to accurately determine the amount of anthropogenic CO 2 that a given rock mass can sequester through mineral carbonation. The program estimates the modal mineral abundance of ultrama ﬁ c rocks to aid in MCP estimation. This tool is designed to be cost-e ﬀ ective and tailored for use by the mining industry, utilising abundant lithogeochemical data to evaluate their deposit as a potential substrate for industrial mineral carbonation operations. The paper introduces the MCP calculator, outlines a framework for developing the MCP parameter, and presents an example of its application. The calculator is speci ﬁ c to the mineral assemblage investigated at the Turnagain ultrama ﬁ c complex in northern British Columbia but can be adjusted to study comparable deposits. The paper acknowledges that using waste rock in a mineral carbonation operation requires economic and practical decisions beyond the scope of the research.


Introduction
Mineral carbonation is a process by which carbon dioxide (CO2) reacts with minerals in rocks to form solid carbonates [1]. The reaction can naturally occur over geological timescales or can be artificially accelerated for carbon-capture and storage purposes [2]. Mineral carbonation has the potential to significantly mitigate greenhouse gas emissions by removing CO2 from the atmosphere and permanently storing it in solid form [3,4]. The process can utilise waste materials, such as mine tailings, as a source of minerals, providing an additional benefit of reducing waste while mitigating emissions [5,6]. Mineral carbonation has been the subject of numerous studies in the fields of geology [7][8][9], mineralogy [10], engineering [11][12][13][14], and environmental science [15,16]. These studies have focused on several key aspects of the process, including, but not limited to, the reaction kinetics [11,17,18], mineral reactivity [17], rock types [9], and life cycle analysis [4,15,16]. Different rocks have been studied to understand their potential for mineral carbonation and the effectiveness of different carbonation methods for specific rock types [19][20][21]. The mineralogy of the rock and the presence of reactive minerals, such as olivine, serpentine, and brucite, play a critical role in determining the suitability of a rock for mineral carbonation [8,22]. In practice, the capability of a rock to sequester anthropogenic CO2 through mineral carbonation is determined through a combination of laboratory experiments and field studies, which provide information on the kinetics and thermodynamics of mineral carbonation reactions [23]. The typical procedure starts with a complete characterisation of the rock mineralogy and composition, followed by reactivity tests at the laboratory and pilot scales to measure the rate and extent of the carbonation reaction for the specific rock [24,25]. In the case of field studies, once carbonation has been initiated, the process is monitored to determine the extent of carbonation over time and ensure the stability of the carbonates in the long-term [7,26,27].
Although following these steps provides valuable information for designing and implementing mineral carbonation, conducting experiments on every available rock is impractical. Therefore, developing mathematical models is necessary to predict the performance of mineral carbonation for different rock types and conditions. Researchers worked extensively on kinetics and thermodynamic models that simulate the mineral carbonation reaction rate and extent over time and under conditions, such as the temperature, pressure, and pH. They also predict the type and stability of carbonates and the reaction conditions required for mineral carbonation [28,29].
The possibility of using mine tailings from ultramafic formations for ex situ mineral carbonation and additional mineral recovery has been the subject of several studies conducted by mining companies and research institutes [30][31][32]. Incorporating waste rock and tailings in a mineral carbonation operation involves several practical and economic considerations. These include developing a viable industrial-scale mineral carbonation process, determining the economic feasibility of mining enough substrate rock to supply a mineral carbonation plant, and, crucially, establishing a financial incentive to pursue mineral carbonation operations in the first place [33]. Diamond drill core sampling and analysis commonly provide lithogeochemical data in mining operations. Hence, there is a significant benefit in developing a method to utilise this data to determine a deposit's mineral carbonation potential and inherent internal variability.
The geological mapping of ultramafic deposits can provide an initial estimate of the volume of ultramafic rocks at a deposit, whether carried out by a national geological survey or a mining company. By analysing bulk rock samples taken during surface mapping and assaying them for principal cation oxides, the RCO2 of that rock mass can be derived, which can be used to provide a preliminary estimate of the capacity of that rock mass to sequester CO2. RCO2 is a widely used parameter to estimate the mineral carbonation potential of an ultramafic deposit. This measure was initially proposed by Lackner, Wendt [34]. RCO2 is the mass of ore required to convert a unit mass of CO2 into carbonate. The initial measure was based on the molar concentration of magnesium and calcium, as described by Lackner, Wendt [34,35]. However, the calculation was later revised [36] to include ferric iron as shown below: where ΣCa 2+ + Fe 2+ + Mg 2+ are the sum of the molar concentrations for the specified cations, and MWCO2 is the molecular weight of CO2. Theoretical minimum RCO2 values have been calculated as 2.1, 1.6, and 2.7 for serpentine, olivine and wollastonite, respectively [37]. This initial approach is sufficient to determine if the deposit justifies additional exploration as a potential material source for a mineral carbonation operation. To develop a more informed estimation, a more detailed investigation is required.
In this research, we aimed to create the mineral carbonation potential (MCP) parameter, which can be generated from existing lithogeochemical data in three-dimensional space. The parameter is exclusively intended for use in the mining industry due to data availability, typically only found in deposits with economic mineral importance. The MCP parameter measures the potential for a rock mass at a specific location to store carbon dioxide via mineral carbonation. Due to varying deposit types, mineralogy, and the range of mineral carbonation methods present at a research stage, generating an applicable parameter to suit all these factors is challenging. In this article, we first propose a framework for developing an MCP parameter. Then, this approach is illustrated using lithogeochemical data from the Turnagain ultramafic complex in northern British Columbia, Canada, as an example. Although the calculator is tailored to the specific mineral assemblage studied at Turnagain, the estimation process can be modified to analyse similar deposits.

Data Collection and Mineralogical Investigation for the Estimation of MCP
A well-structured data collection program is imperative to assess ultramafic deposits for their potential in mineral carbonation. Such a program should include a sampling protocol explicitly designed for a mineral carbonation investigation, utilising various geological, mineralogical, and geochemical analytical techniques. This section provides an overview of the data types commonly available during mining exploration or operation stages. In the following section, methodology is proposed for manipulating existing data to estimate the MCP of ultramafic deposits, using the available data from the Turnagain ultramafic deposit as a case study.
The data preparation required for the MCP calculator is relatively simple. All lithogeochemical assay results obtained from drill core samples at the ultramafic complex can be incorporated into the calculator. However, it is essential to store each drill core assay interval in a database and link it to the corresponding drill hole. A unique sample identification should be assigned to identify the depth interval represented by each sample. It is crucial to have a geoscientist classify the rock by logging the drill core sample intervals based on their geological characteristics. A representative sample from each interval must also be assayed for major cation quantities.
A mineralogical investigation is crucial to estimate the mineral carbonation potential (MCP) of ultramafic deposits. Mining or exploration companies usually conduct a detailed mineralogical study of the central rock units in targeted ultramafic deposits. This investigation involves optical microscopy of thin sections cut from representative samples of each lithology and an X-Ray diffraction analysis (XRD) to determine the principal mineral assemblage of each lithology. Estimating the proportions of magnesium and calcium silicates within each rock unit is essential due to the variable success of using differing magnesium and calcium silicates in experimental mineral carbonation. Identification of the major mineral phases for each lithology should be determined through optical microscopy and XRD, which can provide a preliminary assessment of the abundance of each mineral and the degree of alteration present for each lithology. Quantitative methods, such as thinsection point counting and XRD with Rietveld refinement, can provide a more precise determination of the principal mineral quantities. Mineralogical investigation can also aid any modal mineral estimation undertaken, such as electron microprobe analysis of individual minerals, which can aid in determining their precise geochemical makeup. For instance, determining the Mg 2+ /Fe 2+ ratio (magnesium number) of olivine minerals is crucial. It can vary significantly between rock types and provide a misrepresentation of the magnesium cation quantity available for mineral carbonation from olivine. A better understanding of the mineralogy of ultramafic deposits can improve the estimation of their MCP.

Case Study: Turnagain Ultramafic Deposit, British Columbia, Canada
The Turnagain ultramafic deposit was evaluated as a potential site for mineral carbonation based on preliminary geological mapping that estimated a volume of approximately 10.8 km 3 of ultramafic rocks. This deposit has been explored and drilled by Hard Creek Nickel Corporation to develop a resource estimate for a potential nickel sulphide mining operation.

Preliminary Data Analysis
This research performed X-ray fluorescence (XRF) on 16 samples to estimate their major oxide content. The loss on ignition (LOI) value, a key indicator targeted for the estimation of serpentine in the rocks, was determined. Optical mineralogy was also undertaken by geoscientists at the Turnagain ultramafic complex [38,39]. The major cumulate rock units at Turnagain comprise dunite, wehrlite, hornblendite, and minor dioritic phases. More details on the predominant mineral assemblage of each of these significant rock types were published in a Master's thesis by Scheel [39]. This provided the basis for estimating the mineral assemblage using the MCP calculator in this research.
Five specific rock types were targeted at the deposit based on their potential to sequester CO2 via direct aqueous mineral carbonation. These rock types were green dunite, dunite, wehrlite, olivine clinopyroxenite, and serpentinite. The rock identification of each interval was extracted from the Turnagain exploration database, and the drill hole collar and survey information, the logged geology codes, and relevant ICP-ES geochemical assay results were added to the database. The rock types were selected based on the recorded high magnesium contents and expected high olivine contents.
The green dunite lithology was identified as separate from dunite due to the olivine crystals' higher Mg 2+ /Fe 2+ composition. The MCP calculator requires the magnesium number to be identified for each lithology to determine the proportions of magnesium and iron cation quantity that should be estimated into olivine during the mineral estimation process of the MCP calculator. The magnesium number is calculated using the following formula: Magnesium number= Mg The assigned rock codes of each lithology studied and their average magnesium number can be found in Table 1. Rocks identified as serpentinite were also analysed to test the effectiveness of the MCP calculator in estimating olivine and serpentine abundances.

The MCP Calculator
The MCP calculator was developed to provide a more accurate estimate of the mineral carbonation potential of a rock mass by considering the variability in silicate mineralogy throughout a rock mass, as different silicates have variable mineral carbonation capabilities. The calculator provides the volume of potential carbonate minerals that can be formed for a discrete volume of substrate rock by converting the major cation weight percentage composition into major oxides, estimating their abundance based on a linear best fit stoichiometric conversion into a mineral assemblage suite defined by the user, and estimating the potential stoichiometric magnesium silicate-to-magnesium carbonate conversion achievable by that discrete volume of rock based on the conversion extent for a given mineral through experimental mineral carbonation.
The MCP calculator requires adaptation on a deposit-by-deposit basis to account for the specific mineral assemblage local to the deposit. The instructions for data input into each worksheet and potential modifications that can be made to customise the program to a specific ultramafic deposit are provided in the Ph.D. thesis by Jacobs [33]. The program consists of four worksheets: the Sample Preparation Worksheet, the Mineral Composition Calculator, the MCP Calculator Worksheet, and the Multiple Sample Worksheet.
To facilitate a better understanding of the program's structure, visual representations of these worksheets are displayed in Figure 1. Subsequent sections will provide a more comprehensive explanation of each worksheet. The program is provided as a macro-enabled worksheet allowing the user to process up to 30 samples simultaneously.  The Sample Preparation Worksheet is designed to allow users to enter rock types, rock codes, and their average magnesium numbers. Additionally, it provides an option to enter the measured major cation values for each sample in either weight percent (wt.%) or parts per million (PPM). The sulphur (S) content can be converted to FeS or FeS2. The estimated component abundance will be automatically generated and can be used in the MCP Calculator Worksheet.
Based on the field Herrmann and Berry [40] method, the Mineral Composition Calculator Worksheet is a tool designed to determine the chemical composition of individual minerals to estimate their assemblage. The worksheet is divided into two main sections: (1) the mineral calculator and (2) the mineral library. The mineral calculator allows users to calculate the major oxide values of each mineral. The mineral calculator for oxide conversion includes the major cations commonly found within minerals from ultramafic rocks. The user must enter the name of the mineral being assessed and input the number of each respective cations for the mineral formula. The molecular weight of the mineral is automatically updated by summing the totals of the number of each cation multiplied by the molecular weight of each cation converted into their respective oxides. The weight percentage of each chemical component making up the mineral is generated and sums up to 100%. The general formula used to calculate the weight percentage of each chemical component is as follows: The mineral library in this worksheet provides a place where the results for each mineral can be pasted for storage before being inputted into the MCP calculator.
The MCP Calculator Worksheet consists of four sections, which provide data sources for the calculator and display its results. The working part of the worksheet contains a series of data grids that require no input from the user unless the estimated mineral assemblage differs from the test data. These grids comprise a series of equations for modal mineral estimation following a general linear stoichiometric estimation procedure through numbered stages. The first section contains the imported modal profile of the entire sample data analysed, providing the data source for the iterative working data grids in the same worksheet. The second section is designed to store any ideal modal mineral stoichiometric compositions to be estimated. Although these cells are inactive, they can be incorporated within the working data grids for mineral estimation. Section 3 of the worksheet displays the estimated mineral compositions in weight percentage, which can be customised to match the mineral assemblage identified at the studied deposit. This section is populated from the working data grid labelled by a series of stages. Each stage represents a calculated section of the linear stoichiometric estimation process. The MCP calculator can be customised for different deposits by altering the calculation formulas within the worksheet.
The stoichiometric linear mineral estimation stages created by the author for use at Turnagain are summarised in Figure 2. Stage one within the calculator is the normalisation of the modal profile to 100%, which assumes that the remaining weight percent of each sample tested comprises SiO2 and H2O. Stages two, three, and four assign the modal weight percentages from the sample data preparation worksheet for the sulphides (FeS or FeS2), base metal quantity, and all other accessory components not identified as a component of the major mineral assemblage. The TiO2 and Cr2O3 component of the major oxide profile is identified from quantitative XRD based on the Turnagain samples as ilmenite (FeTiO3) and chromite (FeCr2O4), respectively. These mineral abundances are estimated using all the respective oxides available in each sample, and the iron content is renormalised as FeO for further estimation. Suppose minerals other than ilmenite and chromite are the major mineral phase of their respective oxides. In that case, the associated formulas within the data grid cells can be adjusted as needed. The analysis of major cations in multiple samples from Turnagain revealed high levels of aluminium (up to 20 wt.%), concerning the generally low mean aluminium contents of ultramafic rocks. To consider this elevated aluminium content, the weight percent anorthite (CaAl2Si2O8) content is estimated in stage seven. Plagioclase is present in a significant proportion of the samples, and anorthite can be substituted for a generic plagioclase mineral formula for estimation or other plagioclase end-member mineral formulas. Enstatite was identified as another plagioclase mineral at Turnagain, but estimating both enstatite (MgSiO3) and forsterite (Mg2SiO4) needed to be more practical. Hence, forsterite being the dominant mineral, it was selected over enstatite for the estimation.
The low concentrations of K2O, MnO, Na2O, and P2O5 measured in the samples analysed indicate that these major oxides are incorporated into the accessory minerals in stage three. They mainly form minor components of the silica minerals already estimated. The presence of the minerals from XRD or similar analytical techniques should define the choice of plagioclase mineral estimation. This is a deposit-specific variable.   • Renormalisation to 100 % modal mineral balance The magnesium-bearing minerals are the principal controlling factor based on estimated MCP values. In stage eight, the magnesite (MgCO3) content is estimated if the CO2 content is accurately determined. If no accurate CO2 content is available, the estimated rock sample is considered to have no carbonates. The working data grid can be adapted to include calcite (CaCO3) if required. The principal mineral source of calcium at Turnagain was identified in XRD analysis as diopside (CaMgSi2O6) and estimated in stage nine. Diopside is particularly abundant in wehrlite and olivine clinopyroxenite. Other pyroxene minerals could be substituted in the calculator at this stage. The diopside estimate uses all the CaO content of each sample. The mineral estimate is extracted from the associated stoichiometric MgO and SiO2 contents. The stage one value and the MgO weight abundances are renormalised for use in later stages of the estimation.
Olivine, serpentine, magnetite/hematite, quartz, and periclase are estimated in stages ten to fourteen. These are estimated through a series of iterations that stoichiometrically best estimate the modal weight percentages based on user-defined factors. The estimation process for Turnagain ultramafic rocks starts with olivine since it is the most abundant mineral in this deposit. The calculator then assumes that any remaining unaccounted magnesium oxide (MgO) from the earlier estimation stages is present in this mineral, resulting in a complete estimation procedure. The magnesium number assigned to each sample within the modal profile determines the proportions of MgO, FeO, and SiO2 for the olivine estimation. The olivine conversion matrix in the worksheet controls the active cells for use in the olivine estimation. This is referred to within the series of iterations and can estimate olivine with any assigned magnesium number between 70 and 99.
The Initial olivine quantity is estimated by first calculating the SiO2 and FeO content required to form olivine with the assigned magnesium number to fix all the MgO into olivine. The SiO2 + H2O content from stage one and the FeO content is subsequently recalculated. If the SiO2 + H2O or FeO content is in a deficit, the MgO concentration is renormalised to account for this deficit, as such, to return zero values for SiO2 + H2O and/or FeO. If the MgO content is renormalised for a SiO2 deficit, the olivine content is determined. If no serpentine is estimated, the remaining FeO is assigned to magnetite/hematite. The MgO content is assigned as periclase in stage 14. In the final stage, the sample is renormalised to balance the modal mineral composition to 100% stoichiometrically, completing the modal mineral profile for the sample.
The serpentine content is estimated through a series of iterations in stage 11. Where necessary, the olivine abundance is adjusted in stage ten to appropriately balance the FeO, MgO, SiO2 + H2O contents within the estimated mineral assemblage. The general serpentine formula used for estimation was Mg3Si2O5(OH)4. The active cells for internal calculations are determined from the serpentine conversion matrix to proportion the mineral's MgO and SiO2 + H2O components. When estimating olivine, serpentine abundance is first estimated to mineralogically fix any MgO released from an FeO deficit. Enough SiO2 + H2O must be available to fix the surplus MgO.
Suppose there is a surplus SiO2 + H2O content above a user-defined quantity (a 1 wt.% tolerance is defined at Turnagain based on XRD analyses), the required MgO to estimate the amount of serpentine is taken from the provisional olivine weight percent. This releases FeO from the estimated olivine quantity and assigns it to stage 12. The olivine weight percent is renormalised to account for the removed MgO, FeO, and SiO2 content. This process is repeated through three further sequences to balance the MgO, FeO, and SiO2 iteratively. This normalises the olivine, serpentine, magnetite/hematite, and quartz ± water weight abundance following the iteration. The minimum olivine content is set to 10% within the equation formulas at Turnagain, if each sample contains at least 10 wt.% olivine. This is to avoid generating negative olivine weight percentages. Further iterations could be added if samples have an extremely high serpentine content resulting in a large residual SiO2 + H2O estimated content remaining after three iterations.
In stage 12, all the FeO content is assigned to magnetite/hematite. The MCP calculator has no method of determining the ratio of magnetite versus hematite in the samples being examined, and as such, the FeO content is grouped as both. The mineral profile is normalised to 100% by assigning the remaining weight percent to silica ± water. The final estimated modal mineral composition profile is added to Section 3 of the worksheet. If the reaction extents are not known for each mineral, the values produced in this section could be used to model the mineral distribution within the deposit.
The final MCP values for each sample estimated are determined in the fourth section. This is automatically completed, dependent on the values reported in the previous section and those defined by the user in the mineral carbonation extents table. The mineral carbonation extents table defines either the maximum or average mineral carbonation reaction extent of the selected mineral carbonation process. The table can contain up to five minerals capable of supplying cations that can be bound to carbon to form carbonates. In the example provided by the Turnagain ultramafic deposit, olivine, serpentine, and diopside are potential silicate minerals that can provide magnesium and/or calcium to be bound with carbon dioxide. At Turnagain, olivine is targeted as the principal contributor to the MCP value; however, this could easily be adapted to focus on serpentine.
The MCP values are determined using Equation (4): MCP= S min a ×R x + S min b ×R x +…+ S min n ×R x ×CO 2 wt.% 100 (4) where MCP is the volume weight percentage of magnesite that can be formed from a given volume of rock. Smin is the weight percentage of substrate mineral available for carbonation (a, b,…, n), calculated using Equation (5), Rx is the average or maximum extent of silicateto-carbonate conversion achieved through experimental mineral carbonation for the given mineral, as in Jacobs [33]. CO2 wt.% is the volume weight percent carbon dioxide required to sequester the available MgO as magnesite.
S min = wt.% Mineral×wt.% MgO in mineral 100 The mass of CO2 that a given rock mass can sequester is calculated as Mass CO 2 seq= Mass of rock × MCP 100 × wt.% CO 2 in magnesite 100 (6) where wt.% CO 2 in magnesite 100 =0.522

Reliability and Limitations of the MCP Calculator
This section discusses the limitations of the MCP calculator, providing insight into the sources of variation in the calculated MCP values when using this method instead of other methods.
An incomplete modal profile without information on SiO2 and H2O contents only partially estimates mineral proportions and CO2 sequestration capacity in the MCP calculator. Major cation data input creates an internally consistent profile through conversion to major oxides using fixed conversion factors.
A high sulphide content in rock samples can result in lower FeO values in the mineral estimation stage when Fe is incorporated as FeS or FeS2. If assigning sulphur content as iron sulphides is problematic, the calculator can estimate this based on the inputted sulphur content and a user-defined sulphur mineral formula within stage two of the MCP calculator.
The MCP calculator's estimation process is linear, which imposes certain limitations on the number of minerals that can be included in the mineral profile. It is recommended to use generic mineral formulas for accessory minerals other than the principal magnesium and calcium-bearing silicates, which are the primary focus of the calculator. In the example provided in this work, anorthite and diopside are used as principal sources of aluminium and calcium, respectively, due to their familiar presence in rocks and their contribution to the modal profile. Realistically estimating the remaining SiO2 content is essential for estimating the magnesium silicate mineral, as the magnesium-to-silica ratio controls the olivine and serpentine proportioning in the mineral profile.
The iterative estimation of olivine in stage ten of the MCP calculator assumes SiO2 as the remaining weight percentage based on assumptions about the rock type and magnesium number. However, variations in magnesium numbers within a rock type can affect the accuracy of the estimation and may require internal renormalisation of the iteration equations.
After estimating olivine, any unaccounted oxide percentages are assumed to be SiO2 + H2O, with MgO being removed from the estimated olivine weight percentage to satisfy the stoichiometric serpentine formula. The SiO2 content removed from the olivine matrix is assumed as SiO2 and H2O to balance the olivine and serpentine proportioning, resulting in an internally consistent iteration.
Using the data grids presented in the MCP calculator, the mineral estimation consistently provides a modal mineral profile totalling 100%. MCP value calculations utilising the user-defined reaction extents are consistent for all samples. These can be adapted to account for improvements in the mineral carbonation extent for a given mineral. The MCP value calculation relies only on minerals with an assigned reaction extent. It does not consider external factors within the modal profile that may influence mineral carbonation on a sample-by-sample basis.

Results and Evaluation of the MCP Calculator
The MCP calculator was evaluated using lithogeochemical data from the Turnagain ultramafic complex. Sixteen samples were chosen for validation due to a full suite of geochemical analyses being undertaken. The mean stoichiometric olivine to magnesite achieved from experimental carbonation based on the samples was used as the average carbonation extent for olivine at the deposit. The effectiveness of the calculator was evaluated using the XRD analysis of mineral carbonation products. For XRD analysis, all samples were reduced into fine powder to the optimum grain-size range (<10 µm).
Step-scan X-ray diffraction data were collected over a range 3-80 °2θ with CoKa radiation and a step size of 0.04 °2θ. The X-ray diffractograms were analysed using the International Centre for Diffraction Database (ICDD) PDF-4 and Search-Match software by Siemens (Bruker, Billerica, MA, USA). X-ray diffraction data of the samples were refined with Rietveld program Topas 4.2 (Bruker AXS).

Major Cation-to-Oxide Conversion
To assess the accuracy of the MCP calculator, all 16 samples from the Turnagain ultramafic deposit were analysed via ICP-ES at Acme Analytical Laboratories in Vancouver, British Columbia. The ICP-ES data were added to the sample data preparation worksheet of the MCP calculator program, and a modal profile of the major oxides for each sample was generated. The resulting major oxide values were compared to those obtained using XRF from the same samples.. The results were validated using the sample data preparation worksheet. Table 2 shows that the total absolute residual values for the estimated values were all below 5 wt.%, which is acceptable. The residual values are the calculated difference between estimated and measured values. The total residual values sum the absolute differences observed between the estimated and measured components. The most significant difference between estimated and measured values was for MgO, which is unsurprising given that it makes up almost 50% of each sample. However, the difference between estimated and measured values was still under 3 wt.%. Figure 3 compares measured MgO, Fe2O3, and CaO oxide values from XRF analysis with estimated oxide values from the MCP calculator program for all 16 samples from the Turnagain ultramafic complex. The residual values are the calculated difference between estimated and measured values. The total residual values are the sum of the squared differences observed between the estimated and measured components. * Value not measured or estimated.

Reactant Modal Mineral Estimation
A representative sample from the Twin Sisters ultramafic complex and 16 samples from the Turnagain ultramafic complex were subjected to analysis using the MCP calculator. The accuracy of the MCP calculator was evaluated by comparing its results with those obtained through quantitative XRD with Rietveld refinement. An example of the Rietveld refinement plot of the sample 04-25-AJ6 carbonated product is presented in Figure 4. The Twin Sisters sample was analysed in the MCP calculator to evaluate olivine and serpentine estimation on what is known to be a relatively unaltered dunite sample with high olivine content. The results are provided in Table 3. The sample tested was assigned a magnesium number of 96 [41] in the calculator. The results show an excellent correlation between the measured (XRD) and estimated (MCP calculator) values.  The modal mineral abundances estimated using the MCP calculator based on samples tested from Turnagain are compared with the mineral abundances determined from the same drill core samples via XRD in Figure 5. Some minor phases were identified in  In Figure 5, the mineral abundances of olivine, serpentine, diopside, and magnetite estimated using the MCP calculator program are compared graphically with the mineral abundances determined via quantitative XRD. According to Scheel [39], the degree of serpentinisation in the Turnagain ultramafic deposit varies significantly from 0 to 100%, with an average of approximately 10% based on field and drill core analysis. As shown based on both XRD and MCP calculator analyses, the olivine and serpentine modal abundances fluctuate between 40% to 100% and 0% to 40% in all tested samples, respectively. The discrepancy between the calculated XRD modal abundances and the MCP calculator estimated abundances for olivine and serpentine could vary up to 20 wt.% from sample to sample.
When the estimated serpentine and olivine quantities are compared to XRD analyses, the mean absolute error for serpentine estimation is 11 wt.%, and for olivine estimation, it is 12 wt.%. Figure 6 displays a strong positive correlation between the absolute difference in olivine and serpentine estimated quantities between the calculated XRD and estimated MCP calculator result, indicating that an overestimation of olivine values based on the MCP calculator leads to an underestimation of serpentine and vice versa. This relationship has important implications for the MCP value estimation of the tested sample, particularly if the mineral carbonation method selected can only sequester one substrate (either olivine or serpentine) and not both, as the degree of error in mineral estimation is directly related to the absolute error.  The serpentine content estimated based on the MCP calculator is renormalised based on the degree of crystallisation determined via XRD analyses. Still, it should be noted that the XRD analyses only provide a partial estimate of the serpentine content. This estimate is derived by assuming that all amorphous material is serpentine during Rietveld refinement. The total MgO, SiO2, and H2O contents estimated in the MCP calculator's sample data preparation worksheet positively correlate with the same components measured via XRF [33]. The absolute error in olivine and serpentine modal abundances between the MCP calculator program-derived values and the calculated XRD values can be attributed to several factors, including the estimation error in the individual oxide components used for the olivine and serpentine estimation at the sample data preparation stage of the MCP calculator, error in proportioning these components into the respective mineral formulas, and/or inaccuracies in the serpentine and olivine quantities resulting from the partially quantitative XRD analyses being renormalised.
It is easier to accurately quantify the estimation error expected by the MCP calculator program with a precisely determined serpentine modal abundance for each sample tested. However, the guideline mean estimation errors of 11 wt.% and 12 wt.%, compared with semi-quantitative results of the same samples, seem reasonable. The MCP calculator program can determine the high serpentine contents in samples that are also measured with high serpentine contents via XRD with Rietveld refinement. This is demonstrated by analysing samples from the Twin Sisters ultramafic complex, where the calculator determined an extremely high olivine content with an acceptable estimation error. Although some degree of error is expected when estimating rock samples' olivine and serpentine contents, the MCP calculator program is an effective tool for estimating the degree of serpentinisation throughout an ultramafic deposit from lithogeochemical data. Table 4 compares the MCP values obtained from the analysis of sixteen samples from the Turnagain ultramafic complex using the MCP calculator program with the estimated MCP values projected from XRD analyses of reactant samples. The table also includes the measured magnesite contents obtained from quantitative XRD analyses of the experimental mineral carbonation products. The MCP calculation in Equation (4) only considered olivine as the mineral contributing MgO for mineral carbonation since it was the only mineral that contributed to mineral carbonation during experimental mineral carbonation. The experimental results are reported in the Ph.D. thesis by Jacobs [33] The MgO weight percentage available for mineral carbonation from olivine in Equation (4) depended on the magnesium number of the estimated olivine content defined in Equation (4). These magnesium numbers were applied in both the MCP calculator and the projected MCP values based on the quantitative XRD values of the reactant samples. The projected stoichiometric olivine-to-magnesite conversion extent was set at 37% based on the mean carbonation extent achieved using the experimental mineral carbonation procedure and conditions based on material from the Twin Sisters ultramafic deposit.
When comparing the average estimation difference in MCP values between those estimated using the MCP calculator program and those projected using the same reaction extent of 37% from quantitative XRD, there is little difference, at 5.15% to 4.89%, respectively. Both the MCP calculator values and those projected from reactant XRD are consistent throughout, typically ranging between 10% and 20%. The significant differences observed between estimated and measured values appear to be a function of the variability in the mineral carbonation extent. The use of the MCP parameter here is restricted by the industrial mineral carbonation processing regime. Until an established mineral carbonation process has been defined, the accurate use of MCP values is difficult to validate. Some samples, including 04-24-AJ3, 06-110-AJ2, 06-110-AJ17, 06-111-AJ1, 06-116-AJ9, and 06-116-AJ12, were estimated within 1% of the observed experimental mineral carbonation values, suggesting promising estimation results. The relative consistency of estimated MCP values and those projected from XRD indicate that the MCP calculator is an effective tool for predicting MCP behaviour, provided that the variable nature of experimental mineral carbonation is considered.

Conclusions
The MCP calculator and the MCP values it generates represent an initial attempt at maximising the use of lithogeochemical data available at mining operations to create a more informed estimate of their CO2 capacity based on the rocks available. In addition to the mining industry, this research can be beneficial to any carbon capture sequestration (CCS) project. This tool provides a cost-effective solution for mining companies to evaluate their deposit as a potential substrate for industrial mineral carbonation operations without additional geochemical testing. Although the MCP calculator has limitations, such as the use of oxide values converted from elemental measurements and sensitivity to the inherent variability of the carbonation process, in this work, it is shown that it can accurately determine samples with a high serpentine content, as demonstrated by its analysis of samples from the Twin Sisters ultramafic complex. Despite some margin of error, the program effectively estimated the degree of serpentinisation throughout the deposit from lithogeochemical data. The program also determined an extremely high olivine content with an acceptable estimation error, showing its effectiveness in determining the mineral contents of rock samples.