Impact of Standing Column Well Operation on Carbonate Scaling

: Standing column well constitutes a recent promising solution to provide heating or cooling and to reduce greenhouse gases emissions in urban areas. Nevertheless, scaling issues can emerge in presence of carbonates and impact their efﬁciency. Even though a thermo-hydro-geochemical model demonstrated the impact of the water temperature on carbonate concentration, this conclusion has not been yet demonstrated by ﬁeld investigations. To do so, an experimental ground source heat pump system connected to a standing column well was operated under various conditions to collect 50 groundwater samples over a period of 267 days. These ﬁeld samples were used for mineral analysis and laboratory batch experiments. The results were analyzed with multivariate regression and geochemical simulations and conﬁrmed a clear relationship between the calcium concentrations measured in the well, the temperature and the calcite equilibrium constant. It was also found that operating a ground source heat pump system in conjunction with a small groundwater treatment system allows reduction of calcium concentration in the well, while shutting down the system leads to a quite rapid increase at a level consistent with the regional calcium concentration. Although no major clogging or biofouling problem was observed after two years of operation, mineral scales made of carbonates precipitated on a ﬂowmeter and hindered its operation. The paper provides insight on the impact of standing column well on groundwater quality and suggests some mitigation measures


Introduction
Nowadays, it is mostly accepted that greenhouse gases are responsible for global warming [1].Ground source heat pumps (GSHPs) have the capacity to reduce greenhouse gas emissions for heating and cooling buildings [2], which represented 61% of the total energy consumption in the commercial and institutional sector in Canada in 2016 [3].The GSHPs are sustainable technologies [4].Due to their renewable aspect and mostly their high energy efficiency, the number of GSHPs has been rising worldwide over the past 15 years [5].Among the different types of ground heat exchangers, standing column wells (SCWs) have proven to be a promising solution because of their high heat exchange rate and lower first costs [6,7].Another significant advantage is that SCWs can be installed more easily in urban areas as they require a smaller footprint than closed-loop systems, and less groundwater than open-loop systems [8].
Basically, an SCW is a coaxial uncased well (see Figure 1), in which the groundwater is continuously recirculated.This brings the groundwater into direct contact with the surrounding rock and allows heat transfer by conduction along the borehole wall [6,8].During peak periods, part of the circulation flow rate is discharged outside the SCW.This process creates a cone of depression around the SCW and induces a convergent groundwater flow, increasing advective heat transfer that enhances the thermal performances of the SCW [8,9].However, since the groundwater is directly used as the heat carrier fluid, clogging processes can affect the SCW [10] and the various mechanical devices [11] connected to it.Indeed, precipitation can increase resistance to groundwater flow inside the pipes, obstruct the heat exchanger and the submersible pump.The latter can increase pump fatigue, which could lead to its premature failure [11].Besides, precipitation within the heat exchanger can also reduce the energy supplied to the GSHPs [11,12] and lead to a lower coefficient of performance.For instance, a deposit of 0.8 mm on a heat exchanger surface can increase the GSHP consumption by 19% [13].Several elements are indicators of active clogging processes within GSHP systems.These processes can have a physical [14,15], microbiological [10,16,17] or chemical [18][19][20] source.However, in most cases clogging appears by a combination of more than one of those sources [21,22].

Injection well
Many factors can promote clogging processes in carbonate environments.Indeed, the higher partial pressure of carbon dioxide in groundwater with respect to the atmosphere promotes CO 2 degassing, which in turn increases calcite nucleation and precipitation.In a system using groundwater as the heat carrier fluid, these processes are active at the free surface since pumping operations induce a water level variation [23] and a turbulent flow inside the well [20].The same processes are also active in the injection well [24] or in the hydraulic circuit if zones of low pressure are present [25].The dynamic operation of a system has also an impact on water quality.For instance, at system startup, the particle load increases in the well [15] while stagnant conditions occurring during downtime phases promote oxygenation, that in turn favors microbial abundance and diversity in the well [26].
The operation of a GSHP modifies the groundwater temperature (hotter during cooling mode, cooler during heating mode).Recently, the operation of a SCW has been simulated with a coupled thermo-hydro-geochemical model [8,27,28].These works clearly showed that during heating mode, the continuous recirculation of cold groundwater in a SCW promotes calcite dissolution and enrichment of Ca 2+ .On the opposite, it was observed numerically that the cooling mode promotes precipitation mostly in the above-ground equipment.This advanced model was further exploited and revealed that sealing the SCW to prevent CO 2 degassing decreased by 33.3% the precipitation rate of calcite [28], and that a constant bleed rate of 10% reduced the calcite accumulation on borehole walls due to smaller temperature variations in the SCW [27].
The previous numerical results were obtained considering the equilibrium reactions and kinetics of calcite [29].However, for the operating temperatures of a GSHP, laboratory tests indicated that calcite precipitation is inhibited by the presence in groundwater of orthophosphates [30], magnesium and organic acids [31,32].Thus, simple geochemical models cannot be used alone to predict calcite precipitation in natural systems and should be complemented by field investigations.Study of biofouling occurring in an SCW has already been study for hydrous ferric oxides [10,33].However, even though calcium carbonates are the most common scale deposit in open-loop residential systems [34], carbonates precipitation has yet not been studied during the operation of real SCW systems.This work aims to fulfill this gap of knowledge by using an GSHP system connected to an experimental SCW operated during 267 days under various cooling and heating conditions.

Site Description
The groundwater samples used to perform this study were collected through the use of a mobile geothermal laboratory that operates a SCW on a large scale.The following section presents the main components of this laboratory, the geochemical composition of the geological units crossed by the boreholes and the chemistry of the groundwater prior to the operation of the laboratory.

Geothermal Mobile Laboratory
The geothermal mobile laboratory (see Figure 1) is located in the city of Varennes, near Montreal in Canada.This infrastructure was designed to acquire information on the thermo-hydro-geochemical processes occurring during the operation of any ground heat exchanger such as open-loop, closed-loop or SCW [35,36].The laboratory is composed of four water-to-air heat pumps of a total capacity of 56 kW.This allows emulation of the cooling and heating loads of a small building.
To prevent clogging in the GSHPs, a plate heat exchanger separates the groundwater loop from a propylene glycol loop that supplies the GSHPs.All the equipment is held in a marine container to easily move it to different sites, as shown in Figure 2a and Table 1.As shown in Figure 1, the groundwater loop includes a submersible pump, two flowmeters (F-GW and F-Bleed), two temperature sensors (EWT and LWT) and four sample points (SP-1 to SP-4).The recording equipment and its accuracy are presented in Table 2.Note that a groundwater treatment unit is also installed just before the plate heat exchanger to reduce clogging and scaling in the various parts of the system.Initially, the goal was to divert 20% of the groundwater to the treatment unit withdrawn by the submersible pump.Due to higher head loss than expected, only 3% to 7% of the total flow rate is actually treated.Since the groundwater is continually recirculated in the SCW, this allows treating all the groundwater of the SCW after a few cycles.The unit is composed of a sediment filter and two water softeners to protect the cation-ion exchange resin, which is regenerated with a brine.A disinfection system completes the groundwater treatment unit.The disinfection system is composed of one micro-filter, two low capacity ozone generators and one UV lamp.The characteristics of these devices are presented in Table 3 and shown in Figure 2b.Thus, sample point 1 (SP-1) is used to collect groundwater at the outlet of the SCW.Sample point 2 (SP-2) allows collection of the treated water while sample point 3 (SP-3) collects the mix of treated water and groundwater just before the plate heat exchanger.After heat transfer occurs at the plate heat exchanger, the water is reinjected into the aquifer, either by the SCW or injection well.The reinjected water is sampled from the sample point 4 (SP-4).

Geology and Geochemistry
The GSHP system located in the laboratory is connected to a 300 m deep SCW (Ø 165 mm) and to a 150 m deep injection well (Ø 165 mm) located at a distance of 10 m.However, due to a delay to install the various pipes in the SCW, the bottom walls collapsed, resulting in an effective depth of 215 m.Notice that the SCW and injection well are not sealed and can therefore exchange gases with the atmosphere.A stratigraphic column of the geologic materials observed during drilling was established by Beaudry et al. [35,36] from the macroscopic observation of 54 samples taken every 6.1 m (see Figure 3).This analysis indicated the presence of an overburden composed of clayed silt from the surface to a depth of 3.1 m.For depths of 3.1 to 239.3 m, the bedrock is composed mainly of gray mudstone intersected locally by beds of siltstone, beige-pink sandstone, beige limestone and igneous rocks.This material is attributed to the Nicolet Formation of the Lorraine Group [37].For the depths of 239.3 m to 305.3 m, a black shale limestone attributed to the Utica Group [37] was observed.Figure 3. Geological stratigraphy and results of a geochemical analysis of the bedrock around the experimental standing column wells (SCWs).The geological units were identified from macroscopic observation of 54 samples taken every 6.1 m by Beaudry et al. [35,36] while the geochemical analysis relies on three composite random samples.
Due to the impact of the geology on groundwater quality, a geochemical analysis of the bedrock was performed on three composite rock samples taken randomly along the borehole (see Figure 3).To analyze each sample, the lithium tetraborate fusion method was used to determine the percentage of SiO 2 , Al 2 O 3 , CaO, MgO, FeO, Na 2 O, K 2 O, MnO.A LECO CS-744 analyzer set at 1200 • C, was used to measure the concentration of S and of total carbon by infared spectrophotometry method.The CO 2 was determined by subtracting the total carbon to the carbon provided by the LECO CS-744 after a hydrochloric acid treatment (HCl at 50%) during 1 hour, according to the norm MA. 310-CS 1.0 [38].The Loss On Ignition (LOI) was set at 1000 • C.An average was calculated from all the samples and results are presented in Figure 3 (right).
The results of the geochemical analyses indicated that the bedrock was mainly composed of silicates, with high concentrations of aluminum, calcium, magnesium and iron oxides.A significant concentration of sulfur was also observed.The percentage of calcium carbonate (CaCO 3 ) in the bedrock is calculated according to Equation (1) [39].
where M CaCO 3 , M Ca are the molar mass.The resulting fraction of calcium carbonate in the samples is 12.78%.Note the Equation ( 1) is initially developed for Netherlands geology.The calcium carbonate percentage indicates that the bedrock is chemically reactive with the groundwater circulating in the SCW or within the fractures [40].It should be noted that a part of the calcium observed in the samples can be present in the form of calcium sulfates (CaSO 4 ), such as gypsum.In addition, due to the high percentage of magnesium oxides and CO 2 , magnesium carbonate or dolomite (CaMg(CO 3 ) 2 ) can also be present in the sample.

Initial Chemical Composition of Groundwater
An initial sampling of the groundwater present in the SCW was performed on 11 November 2016.The results obtained are shown in Table 4 and indicated a total hardness of 134 ppm CaCO 3 and a pH of 8.1.These two values were above the values of 80 ppm and 7 recommended to avoid carbonate precipitation in GSHP systems [13].Thus, carbonate scaling was expected to occur in various parts of the system.Note that it was not possible to perform more that one sampling before starting the operations of the laboratory and thus characterize the seasonal fluctuations of the groundwater chemistry.

Material and Method
First of all, this section presents the experience performed on the SCW.The groundwater samples were collected between January 2018 and October 2018 on a weekly basis (50 samples) under various conditions.These conditions are presented.A mineral scale sampled in January 2019 is presented.This sampling procedure and analysis method are presented.To interpret the groundwater calcium concentration during this operation, batch experiments of bedrock dissolution were performed in Polytechnique Montréal.The setup of batch experiments and the analysis method are presenting in this section.Finally, The chemical equations and software PHREEQC to calculated the kinetics and equilibrium for the batch experiments and the on-off sequence are presented.

Field Sampling and Groundwater Analysis
The test conditions were decomposed into five identifiable sequences.The first one was a heating mode, composed of demand peaks during the day.The second one was an On-Off sequence.The third one was composed of constant heating demand during the day.The fourth one was a period of cooling and recirculating, and the last sequence was a cooling mode.These conditions are summarized in Table 5.
To collect representative groundwater samples, the following experimental protocol was established and upheld.Samples were collected after the submersible pump had been running for at least 5 min (even if a system downtime occurred during the previous days).Three flushes were performed for each sample point (SP-1 to SP-4) to eliminate any dead volume [41] before sampling.The measurements of temperature, pH, electrical conductivity and oxidation-reduction potential were then carried out with a multiparameter probe (YSI Pro Plus and Pro Series).The water was filtered with a 0.45 µm hydrophilic polyvinylidene fluoride membrane (Millex-HV of Millipore).The samples dedicated to the analysis of iron and manganese were immediately acidified with a solution of HNO 3 having a pH lower than 2. All samples were stored at 4 • C for preservation until they were analyzed.
Analysis of the total soluble concentrations of Fe, Mn, Ca and Mg was carried out by atomic absorption with an Analyst 200 Spectrometer.The concentrations of major anions F − , Cl − , NO − 2 , Br − , NO − 3 , SO 2− 4 , PO 3− 4 were determined by ionic chromatography with a Dionex ICS-5000 system.The detection limits of these two analysers were respectively 0.02 mg/L and 0.5 mg/L.A potentiometric titrator HI 901 measured the alkalinity.

Collection and Analysis of Mineral Scales
Installed in 2016, the flow meter F-GW began to malfunction between 10 October 2018 and 3 January 2019, until eventually it was decided to be disassembled.Observation of the flow meter piping indicated that mineral scales 0.8 mm thick were covering the inner piping (see Figure 4).The scales were tested with nitric acid (pH < 2) and reacted proving the presence of carbonates.Some mineral chips were then collected and stored in tubes for 2 days before a microscopic analysis was performed at the Centre for Characterization and Microscopy of Materials in Polytechnique Montréal.The analysis was carried out using three different technologies: a scanning electron microscope (SEM), an energy dispersive spectroscopy (EDS) and an X-ray diffraction.It is worth noting that no other major problems associated to scaling or biofouling were observed so far.

Batch Experiments
A batch experiment was designed to study the calcium equilibrium and dissolution kinetics of the rock samples gathered during the drilling of the experimental wells.To that end, rock samples with grain sizes between 1.25 mm and 0.63 mm (mean 0.66 mm) were selected.The samples were washed with demineralized water, followed by a bath of hydrochloric acid (0.012N) of 10 s to remove all impurities that were fixed to the surface during the drilling process [29].Finally, the samples were rinsed with demineralized water to stop the reaction with hydrochloric acid.The rock samples were then dried at 65 • C to prevent burning any organic matter.The samples relative density was measured at 2.76 following the standard ASTM C128-15 [42].Samples with a surface of approximately 8.25 cm 2 were immerged in 80 mL of demineralized water (manufacture by Mili-DI) at 11 • C in microbiological bottles, as shown in Figure 5.The microbiological bottles were sealed with rubber bushings to perform the dissolution in a closed system.Trapped air was replaced by injecting nitrogen N 2 .An orbital shaker then was used to stir the rock samples and the demineralized water together, without opening the solution to the atmosphere.Previous works indicated the rate of rotation has a pronounced impact at low pH (<5) while at higher pH the grain surface has the largest influence [29].For this reason, the stirring speed was set to 200 rpm, speed at which the grain was suspended in the solution without mechanical shocks resulting of the impact with the microbiological bottles.During the batch experiment, the concentration of calcium and manganese was measured twice per day during one week by atomic absorption.

Geochemical Simulation
The simulation software PHREEQC [43] was used to perform two geochemical simulations.These simulations were used to analyze the dissolution kinetics and equilibrium of the batch experiments and the impact of GSHP downtimes on calcium concentration.The third simulation allowed computation of the saturation index (SI) of calcite in groundwater.The thermodynamic database WATEQ4F was adopted because it accounts for partial carbon dioxide pressure P CO 2 and includes temperature of the solution [43].This database has been used in related past studies [18,24] for these specific reasons.To simulate the dissolution kinetics of carbonate rock samples, PHREEQC uses the temperature dependent kinetics constants of calcite as identified by Plummer et al. [29].
To interpret the equilibrium of a mineral, PHREEQC calculates the saturation index in water by [40]: where I AP represents the ionic activity product calculated from the activities of free ion species.
In Equation (2), K(T) is the temperature dependent equilibrium constant, with the temperature T in Kelvin.For a SI = 0 the mineral is at equilibrium; for SI < 0 the mineral is undersaturated, leading to a dominant dissolution state and if SI > 0 the mineral is oversaturated, leading to the precipitation of the mineral [40].
In this work, two temperature dependent equilibrium constants were used.The first equilibrium constant was for pure calcite and has been proposed by Plummer and Busenberg [44]: An important feature of Equation ( 3) is that equilibrium constant decreases with temperature, leading to precipitation increase and calcium concentration reduction.
The second equilibrium constant considered the fact that magnesium inhibits calcite precipitation and that magnesian calcite (Ca 1−x Mg x CO 3 ) is more stable than pure calcite.To model this situation, the concept of stoichiometric saturation proposed by Thorstenson and Plummer [45] was used.The dissolution reaction for magnesian calcite is then: where x is the molar concentration ratio between [Mg 2+ ]/[Ca 2+ ] in the aqueous and solid solution.The equilibrium constant for magnesian calcite is: where K Mg-Calcite is the equilibrium constant, a Ca 2+ , a Mg 2+ and a CO 2− 3 are the activity of species in aqueous solution.The right side of Equation ( 5) is the I AP used in Equation (2).The activities were calculated by PHREEQC using the final concentrations measured in the batch experiments described previously.At the end of dissolution, an equilibrium state was considered reached.The same kinetics of pure calcite developed by Plummer et al. [29] were used.

Results and Discussion
To identify the factors impacting carbonate scaling in a SCW, the geothermal laboratory and its groundwater treatment system were used between 16 January and 10 October 2018 under various heating and cooling conditions.During this period, the on-board data acquisition system recorded the operation parameters every minute, while groundwater sampling were performed on a weekly basis.

Geochemical Mapping
The analysis results of the samples collected at sample points 1 to 4 are presented in Figure 6.To compare the results coming from a large data set, whisker boxplots were used.The latter present the median, the 25% and 75% percentiles, as well as the minimum, maximum and outliers' values.To observe the impact of the temperature on the water chemistry, the results were separated into two sets for each sample point based on the temperature measured at SP-1 during the groundwater sampling.The first set, shown in blue, is associated to temperatures equal to or lower than the initial groundwater temperature of 11 • C. The second set, shown in red, is associated to higher temperatures.This separation allowed us to study the impact of the cooling and heating mode.
The temperature at the sample points is shown in Figure 6a.Note that the number of samples n were uniformly distributed between the two sets of temperatures.The temperatures ranged between 1.3 • C at SP-4 and 36.2 • C at SP-1, and covered the typical temperatures observed during the operation of a SCW system in heating and cooling mode.The median temperatures were almost the same for a given set.However, the median temperature at the outlet of the plate heat exchanger (SP-4) was 2.3 • C, which was significantly lower than the other medians for the cold dataset.This difference is due to the heat extraction at the plate heat exchanger occurring in heating mode, confirming the separation as being representative of the heating mode.
As shown in Figure 6b, it is clear that the pH evolved as a function of the temperature since the pH was 0.3 units higher in heating mode than in cooling mode.The same pattern was also observed for the calcium and magnesium concentrations shown in Figure 6c,d.Indeed, the median concentrations were systematically lower at higher temperatures, which indicates that calcite precipitation probably occurred when the GSHP system was in cooling mode.These experimental observations are consistent with the results of Eppner et al. [27,28] who used a coupled geochemical model and observed numerically that recirculation of cold groundwater in a SCW promotes calcite dissolution and calcium enrichment, while recirculation of warmer groundwater during cooling mode favors calcite precipitation and calcium depletion.The fact that magnesium concentrations were also lower in cooling mode is a good indicator that dissolution and precipitation of dolomite or magnesian calcite were also active in the GSHP system.This is supported by the energy-dispersive spectroscopy test and the X-ray diffraction test carried out on the mineral scales shown in Figure 4. Indeed, the microscopic analysis pointed out mostly calcite and some magnesian calcite in the scales, confirming that precipitation processes were active.
The effect of the groundwater treatment system was clearly visible on the pH and on the calcium and magnesium concentrations coming from SP-2.Indeed, the median calcium concentration of the two temperature sets dropped by 80.8% and 88.9% at SP-2 with respect to SP-1.A similar result was observed with the magnesium, which dropped by 88.4% and 83.4%.These significant concentration reductions were partly attributed to the cation exchange resin and water softener installed in the groundwater treatment unit.The percentage of groundwater treated being small (3% to 7%), the impact of the treatment system on the concentrations measured at SP-3 and SP-4 was however smaller.
Nevertheless, the median concentrations of calcium and magnesium were significantly below their initial concentrations of 32.4 and 13.0 ppm, indicating that some processes were actively removing calcium from the groundwater.Based on the concentrations measured at the outlet of the groundwater treatment system (SP-2), it is clear that the continuous operation of the treatment unit could explain at least a part of these reductions.A second explanation lies in the probable impact of temperature on the equilibrium constants, which induces calcite precipitation and pH drop when the SCW is operated in cooling mode [27,28].Degassing of CO 2 can also explain a part of the calcium reduction [28].Based on these sole results, it is however impossible to distinguish which process controls mostly precipitation.

Impact of Temperature and GSHP Operation
It was shown in Section 4.1 that groundwater temperatures can impact calcium concentrations at a detectable level.However, this analysis did not allow to identify the impact of more complex GSHP operations.One can now see in Figure 7 the temporal evolution of calcium concentration, temperature and flow rate for various GSHP operation modes.First, note how calcium concentrations varied in time as a function of the GSHP operation mode in a quite complex manner.Indeed, concentrations of calcium ranged between nearly 0 to 62.80 ppm, with sharp increases during system' downtimes.From Figure 7a, it is also clear that the water treatment system reduced systematically calcium concentrations at SP-2.Although less striking, the concentrations measured at the system outlet (SP-4) were also lower than the ones measured at the inlet (SP-1) for most samples.When the groundwater treatment system was active, the average concentration difference between SP-4 and SP-1 was 2.37 ppm.By comparison, this value was 83.1% smaller at 0.40 ppm when the groundwater was not treated.Using the concentration difference between SP-4 and SP-1 and the flow rate shown in Figure 7d, the cumulative calcite mass potentially removed by the groundwater treatment system between 16 January and 10 October 2018 (267 days) was calculated to 33.5 kg (see Figure 7b).This suggests that by continuously removing small amounts of calcium, the treatment system allowed bringing down the concentrations of calcium and could help preventing calcite precipitation even if only a small fraction of the total flow rate was treated.Notice in Figure 7b how EWT and LWT differed during downtime from the temperature measured by the optic fiber installed in the SCW and how the temperature in the SCW rapidly returned to its initial value of 11.0 • C after each GSHP shutdown.Conversely, during a downtime the circulation flow rate being nil, EWT and LWT stabilized at the laboratory air temperature and were not representative of the groundwater temperature.Nevertheless, this indicates that temperatures of 20-25 • C (and even 35 • C in June 2018) developed within the pipes of the GSHP system.Such high temperatures could therefore promote carbonate precipitation and biofouling in the pipes' stagnant water.

SP
Since analysis of the temporal series of Figure 7 did not reveal clear links between the calcium concentrations and the various operational parameters, a stepwise multiple linear regression was performed to identify the variables controlling the measured calcium concentrations.To identify the parameters having the highest impact on observations, the standardized regression coefficients were obtained by converting the explanatory variables to z-scores.Variables having a p-value (F-statistic) greater than 0.05 were deemed insignificant and were not included in the regression model.Alternatively, some variables were automatically combined and included if deemed significant.To capture linear and non linear relationships between the dependent variable (Ca) and the operational parameters of the GSHP system, the latter were transformed, combined and/or averaged as shown in Table 6.For example, the average of EWT and LWT was used to represent the temperature within the SCW.This temperature signal was then used to evaluate the logarithm of equilibrium constant log(K Calcite ).To measure short and long term effects and cyclic patterns, backward moving averages and moving standard deviations were applied on windows of 3 min and 24 h.The heating, recirculation and cooling modes were coded with integer values of −1, 0 and 1, respectively.Because of the calcium spikes observed during the downtime, a variable named Downtime was used to represent the time since system shutdown.A value of 0 was used if the system was operating.Finally, a variable describing the relative decrease in calcium due to the treatment system was used to represent the impact of the water treatment system on calcium concentrations.
Results of the stepwise regression are summarized in Table 6 while Figure 8 illustrates the retained regression model (R 2 = 0.88).Note that several different variables were tested, leading to regression models having R 2 ranging between 0.81 and 0.96.These models were almost always consistent with the results presented in this work.Based on the rank of the absolute value of the standardized coefficients, it is clear that the groundwater temperature explained most of the observed calcium concentrations.The standardized coefficients associated to T (24 h) being negative, this indicates that calcium concentration was inversely proportional to the temperature within the experimental SCW.The second highest standardized coefficient was also negative and associated to the equilibrium constant of calcite K Calcite .Note how these results are consistent with each other and indicate that lower calcium concentrations were expected at higher temperature due to a higher precipitation rate.This statistical result is also consistent with the kinetic reactions of carbonates.
A second important set of variables regroups the downtime and the operation mode.The standardized coefficients were all positives, indicating that calcium concentration was proportional to downtime.Recall that in Figure 7a calcium spikes were clearly present during downtime.As shown in Figure 8, the regression model integrated very well these events that range from approximately 7 to 60 ppm.The next section will show that this quite strong relationship with downtime was linked to the kinetic reactions of magnesian calcite and the time required to reach a chemical equilibrium.The standardized coefficient associated to the operation mode (−1, 0 and 1) was positive.This is not surprising since this variable was coded with −1 for heating (lower temperature) and 1 at higher temperatures and was basically a surrogate of T (24 h).
The third set of variables considered in the regression model gathered the bleed flow rate (F-Bleed) and the pumping flow rate (F-GW), as well as the variables described previously.Analysis of the standardized coefficients indicated that F-Bleed was proportional to the calcium concentrations.The most plausible explanation is that the increased bleed flow rate promoted a flow of groundwater toward the SCW at a higher calcium concentration.The results of the batch tests presented in Section 4.3 supported this explanation.Treatment e f f iciency Finally, the standard deviation of T and F-GW was considered significant and included in the regression model.The fact that the standard deviation of these variables was so often included in the regression is surprising and is probably connected to the lower concentrations measured during periods of high variability.This relationship could be due to the cyclic pattern of the temperature (see Figure 7c) and its impact on kinetic reactions involving calcium.Similarly, changes of the water level in the SCW due to variations of the flow rate might impact CO 2 degassing at the well head and the precipitation rates in the SCW.Additional work will however be necessary to prove these hypotheses.
It is worth noting that several variables were deemed insignificant by the stepwise regression and were not retained to explain the variations of calcium concentrations.For instance, the efficiency of the groundwater treatment system was not included, probably because the small flow rate of the treatment unit (3% to 7%) was too small to induce noticeable impacts.The fact that most of the variables calculated with a 3-min window were not included in the regression indicates that the processes affecting calcium concentration were more a matter of several hours than a few minutes.

Impact of Downtime
The previous section revealed that the duration of system shutdown influenced the concentration of calcium observed in the experimental SCW.To understand the link between these two variables, the results of the batch tests described in Section 3.3 were used to simulate during a shutdown the evolution of calcium within a SCW.The geochemical simulations were performed with the equilibrium constant of calcite (K Calcite ) calculated according to Equation (3), and with the equilibrium constant for magnesian calcite (K Mg-Calcite ).Note that K Mg-Calcite was based on the activities corresponding to the concentrations of calcium, magnesium and alkalinity at the end of the batch test as shown in Table 7.The simulations calculated the evolution of calcium concentrations over 9 days, starting from an initial state in a specific environmental condition and considering a contact with a mineral surface.The initial state of the water corresponded to the initial groundwater quality for experimental SCW as shown in the Table 8.The minerals were either pure calcite or magnesian calcite.The x was calculated at 0.124 with calcium and magnesium concentrations presented in Table 7.The ratio between specific surface and volume were determined by the batch experiments.The same ratio was used for downtime periods.The simulations were performed using a temperature of 11 • C, which corresponded to the initial groundwater temperature and with initial conditions representing either the batch experiments or the experimental SCW (see Table 8).First, observe in Figure 9a how the calcium concentrations obtained with the equilibrium constant of magnesian calcite had a better fit to the measurements than the concentrations simulated with K Calcite .The relative difference between the concentrations reached approximately 30% in the plateau.At a temperature of 11 • C, three days were required to reach the equilibrium.As indicated in Section 3.1, the submersible pump was activated during 5 min before sampling the groundwater in the geothermal laboratory.Thus, the concentrations measured during a downtime and shown in Figure 9b were deemed representative of the concentrations in the SCW at a temperature of approximately 11 • C. Both simulations showed an overall kinetics equivalent to that of the batch tests and in relative adequacy with the concentrations measured during downtime.The calcium concentration increased mainly during the first two days following a shutdown of the GSHP system.Note how the stoichiometric saturation model successfully represented the long term equilibrium condition for calcium concentration, approaching the initial groundwater concentration upon stabilization.These conditions could be observed for a shutdown of 8 days or more.For shorter downtime it seemed that calcium concentrations obtained with K Calcite fitted better with the experimental values, stabilizing at 49.47 ppm.Analysis of the temporal series of Figure 7 showed that a fairly rapid return to the initial temperature of 11 • C was observed after a shutdown of the GSHP system.At this temperature, geochemical simulations showed that an equilibrium concentration was reached in 3 days at levels compatible with calcium concentrations observed initially on site or during downtime.Regression model also indicated that calcium concentrations were proportional to shutdown duration and bleed flow rate.We concluded that the downtime promoted a return to the equilibrium in the SCW, while the bleed rate fostered groundwater inflow to the SCW at an equilibrium concentration higher than what was usually observed.These conclusions are, for the moment, limited to the specific test site studied in this work.

Conclusions
In this work, an experimental ground source heat pump system connected to a standing column well was operated under various conditions and allowed collection of 50 groundwater samples over 267 days.Results of a multivariate statistical analysis and batch experiments pointed out the link between operations and carbonate equilibrium.Even though theoretical studies suggested an impact of operations on carbonate concentration, this was not yet demonstrated on full-scale standing column wells.As a result, the main contribution of the present paper is that temperature constitutes the most decisive parameter to explain the presence of calcium in the standing column wells and confirmed the dependence of calcium to reaction constants.A second important finding is the impact of system shutdowns.Indeed, this operation allows a quite rapid calcium increase in the well, and can promote precipitation and biofouling in the above ground pipes.For the specific test sites studied in this work, the regression model and batch experiments also highlighted a proportional relationship between bleed flow rate and calcium concentrations.
The results obtained in this study also showed that operating a groundwater treatment unit allowed a global decrease in calcium concentration and prevented precipitation of a significant amount of calcite.Although no major clogging and biofouling problems were observed during the first two operation years of the system, the chosen flowmeter was proved unsuitable for the groundwater conditions of the experimental site.Use of other flowmeter types is therefore recommended.To limit precipitation and biofouling during downtime, it is suggested to purge the above ground piping network with an automatic system triggered during a prolonged downtime period.Finally, although the link between CO 2 degassing and temperature variability was tiny, sealing the well could be an efficient way to limit CO 2 degassing and carbonate precipitation.

Figure 1 .
Figure1.Illustration of the experimental standing column well connected to the mobile geothermal laboratory and its groundwater treatment unit.The four sample points (SP-1 to SP-4), the temperature probes (EWT, LWT) and the flowmeters (F-GW, F-Bleed) are also shown.

Figure 2 .
Figure 2. Mobile geothermal laboratory used during this study, (a) general outside view, (b) groundwater treatment unit and (c) view of the SCW head.

Figure 4 .
Figure 4. Visual observation of mineral deposits on the tip of flow meter after two years of operation.(a) Mineral deposits on the flow meter tip, (b) chips used for microscopic analysis and (c) image taken by a scanning electron microscope.

Figure 5 .
Figure 5. Experimental setup for the batch test in a closed system with temperature control at 11 • C. A rock sample having a surface of 8.25 cm 2 is brought into contact with a 80 mL of demineralizing water with setting continuous rotation at 200 rpm.

Figure 6 .
Figure 6.Mapping of (a) temperature, (b) pH, (c) calcium and (d) magnesium concentration of a few parameters as a function of the sample point.The sample points are located at the outlet of the SCW (SP-1) and groundwater treatment system (SP-2), and at the inlet (SP-3) and outlet (SP-4) of the plate heat exchanger.The horizontal red dashed line corresponds to the initial values.The red cross represents the outliers' values.The blue bars are associated to samples having a temperature lower or equal to 11 • C, while the red bars are associated to higher temperatures.

Figure 7 .
Figure 7. Summary of the GSHP operation and evolution of calcium concentration between 16 January and 10 October 2018.Evolution of (a) calcium concentration at sample points 1 to 4, (b) cumulative calcite mass potentially removed by the treatment system computed by interpolating the concentrations every minute, (c) entering, leaving and SCW water temperatures and (d) total, bleed and treated water flow rates.The red shaded areas represent downtime periods caused either on purpose, by mechanical or control problems, or power outage.The blue shaded area represents a period without water treatment of the groundwater.

Heating- 1 :Figure 8 .
Figure 8. Prediction of the calcium concentrations measured at SP-1 using a multiple linear regression model.The operation phases indicated in the legend correspond to the phases identified in Figure 7.Note that the circled symbols highlight a downtime period.

Table 8 .
Initial conditions and parameters used for the geochemical simulations presented in Figure9.

Figure 9 .
Figure 9. Evolution of calcium concentration for a rock sample composed of pure calcite or magnesian calcite as simulated by PHREEQC and comparison with the experimental measurements obtained during the (a) batch tests described in Section 3.3 and (b) downtime of the experiment shown in Figure 7.

Table 1 .
Equipment installed in the mobile geothermal laboratory.

Table 2 .
[36]racy of the sensors installed in the mobile geothermal laboratory, from[36].

Table 3 .
Equipment of the groundwater treatment unit.

Table 4 .
Initial groundwater quality within the experimental SCW on 11 November 2016.

Table 5 .
Summary of test condition.

Table 6 .
Results of a stepwise multiple linear regression performed with explanatory variables converted to z-scores.The temperatures are the average between EWT and LWT.The durations of 3 min and 24 h corresponded to the time interval used to calculate the backward moving average or the standard deviation.

Table 7 .
Species activity at the end of the batch experiment.