Strontium incorporation into calcite generated by bacterial ureolysis was investigated as part of an assessment of a proposed remediation approach for
90Sr contamination in groundwater. Urea ...hydrolysis produces ammonium and carbonate and elevates pH, resulting in the promotion of calcium carbonate precipitation. Urea hydrolysis by the bacterium
Bacillus pasteurii in a medium designed to mimic the chemistry of the Snake River Plain Aquifer in Idaho resulted in a pH rise from 7.5 to 9.1. Measured average distribution coefficients (D
EX) for Sr in the calcite produced by ureolysis (0.5) were up to an order of magnitude higher than values reported in the literature for natural and synthetic calcites (0.02–0.4). They were also higher than values for calcite produced abiotically by ammonium carbonate addition (0.3). The precipitation of calcite in these experiments was verified by X-ray diffraction. Time-of-flight secondary ion mass spectrometry (ToF SIMS) depth profiling (up to 350 nm) suggested that the Sr was not merely sorbed on the surface, but was present at depth within the particles. X-ray absorption near edge spectra showed that Sr was present in the calcite samples as a solid solution. The extent of Sr incorporation appeared to be driven primarily by the overall rate of calcite precipitation, where faster precipitation was associated with greater Sr uptake into the solid. The presence of bacterial surfaces as potential nucleation sites in the ammonium carbonate precipitation treatment did not enhance overall precipitation or the Sr distribution coefficient. Because bacterial ureolysis can generate high rates of calcite precipitation, the application of this approach is promising for remediation of
90Sr contamination in environments where calcite is stable over the long term.
Precipitation reactions influence transport properties in porous media and can be coupled to advective and dispersive transport. For example, in subsurface environments, mixing of groundwater and ...injected solutions can induce mineral supersaturation of constituents and drive precipitation reactions. Magnetic resonance imaging (MRI) and microcomputed tomography (μ-CT) were employed as complementary techniques to evaluate advection, dispersion, and formation of precipitate in a 3D porous media flow cell. Two parallel fluids were flowed concentrically through packed glass beads under two relative flow rates with Na2CO3 and CaCl2 in the inner and outer fluids, respectively. CaCO3 became supersaturated and formed a precipitate at the mixing interface between the two solutions. Spatial maps of changing local velocity fields and dispersion in the flow cell were generated from MRI, while high resolution μ-CT imaging visualized the precipitate formed in the porous media. Formation of a precipitate minimized dispersive and advective transport between the two fluids and the shape of the precipitation front was influenced by the relative flow rates. This work demonstrates that the combined use of MRI and μ-CT can be highly complementary in the study of reactive transport processes in porous media.
Background
A proposed strategy for immobilizing trace metals in the subsurface is to stimulate calcium carbonate precipitation and incorporate contaminants by co-precipitation. Such an approach will ...require injecting chemical amendments into the subsurface to generate supersaturated conditions that promote mineral precipitation. However, the formation of reactant mixing zones will create gradients in both the saturation state and ion activity ratios (i.e.,
a
C
O
3
2
-
/
a
C
a
2
+
). To better understand the effect of ion activity ratios on CaCO
3
precipitation kinetics and Sr
2+
co-precipitation, experiments were conducted under constant composition conditions where the supersaturation state (Ω) for calcite was held constant at 9.4, but the ion activity ratio
(
r
=
a
C
O
3
2
-
/
a
C
a
2
+
)
was varied between 0.0032 and 4.15.
Results
Calcite was the only phase observed, by XRD, at the end of the experiments. Precipitation rates increased from 41.3 ± 3.4 μmol m
-2
min
-1
at
r =
0.0315 to a maximum rate of 74.5 ± 4.8 μmol m
-2
min
-1
at
r =
0.306 followed by a decrease to 46.3 ± 9.6 μmol m
-2
min
-1
at
r
= 1.822. The trend was simulated using a simple mass transfer model for solute uptake at the calcite surface. However, precipitation rates at fixed saturation states also evolved with time. Precipitation rates accelerated for low
r
values but slowed for high
r
values. These trends may be related to changes in effective reactive surface area. The
a
C
O
3
2
-
/
a
C
a
2
+
ratios did not affect the distribution coefficient for Sr in calcite (D
P
Sr
2+
), apart from the indirect effect associated with the established positive correlation between D
P
Sr
2+
and calcite precipitation rate.
Conclusion
At a constant supersaturation state (Ω = 9.4), varying the ion activity ratio affects the calcite precipitation rate. This behavior is not predicted by affinity-based rate models. Furthermore, at the highest ion ratio tested, no precipitation was observed, while at the lowest ion ratio precipitation occurred immediately and valid rate measurements could not be made. The maximum measured precipitation rate was 2-fold greater than the minima, and occurred at a carbonate to calcium ion activity ratio of 0.306. These findings have implications for predicting the progress and cost of remediation operations involving enhanced calcite precipitation where mineral precipitation rates, and the spatial/temporal distribution of those rates, can have significant impacts on the mobility of contaminants.
The spectral induced polarization (SIP) technique is a promising approach for delineating subsurface physical and chemical property changes in a minimally invasive manner. To facilitate the ...understanding of position and chemical properties of reaction fronts that involve mineral precipitation in porous media, we investigated spatiotemporal variations in complex conductivity during evolution of urea hydrolysis and calcite precipitation reaction fronts within a silica gel column. The real and imaginary parts of complex conductivity were shown to be sensitive to changes in both solution chemistry and calcium carbonate precipitation. Distinct changes in imaginary conductivity coincided with increased hydroxide ion concentration during urea hydrolysis. In a separate experiment focused on the effect of hydroxide concentration on interfacial polarization of silica gel and well-sorted sand, we found a significant dependence of the polarization response on pH changes of the solution. We propose a conceptual model describing hydroxide ion adsorption behavior in silica gel and its control on interfacial polarizability. Our results demonstrate the utility of SIP for noninvasive monitoring of reaction fronts, and indicate its potential for quantifying geochemical processes that control the polarization responses of porous media at larger spatial scales in the natural environment.
Ureolytically driven calcite precipitation is a promising approach for inducing subsurface mineral precipitation, but engineered application requires the ability to control and predict precipitate ...distribution. To study the coupling between reactant transport and precipitate distribution, columns with defined zones of immobilized urease were used to examine the distribution of calcium carbonate precipitation along the flow path, at two different initial flow rates. As expected, with slower flow precipitate was concentrated toward the upstream end of the enzyme zone and with higher flow the solid was more uniformly distributed over the enzyme zone. Under constant hydraulic head conditions the flow rate decreased as precipitates decreased porosity and permeability. The hydrolysis/precipitation zone was expected to become compressed in the upstream direction. However, apparent reductions in the urea hydrolysis rate and changes in the distribution of enzyme activity, possibly due to CaCO3 precipitate hindering urea transport to the enzyme, or enzyme mobilization, mitigated reaction zone compression. Co-injected strontium was expected to be sequestered by coprecipitation with CaCO3, but the results suggested that coprecipitation was not an effective sequestration mechanism in this system. In addition, spectral induced polarization (SIP) was used to monitor the spatial and temporal evolution of the reaction zone.
•A hybrid pore-continuum multiscale flow and reactive transport model is proposed.•The multiscale model approach is demonstrated for a mixing-controlled reaction.•The new model provides improved ...predictive capability over single-scale models.
Continuum-scale models, which employ a porous medium conceptualization to represent properties and processes averaged over a large number of solid grains and pore spaces, are widely used to study subsurface flow and reactive transport. Recently, pore-scale models, which explicitly resolve individual soil grains and pores, have been developed to more accurately model and study pore-scale phenomena, such as mineral precipitation and dissolution reactions, microbially-mediated surface reactions, and other complex processes. However, these highly-resolved models are prohibitively expensive for modeling domains of sizes relevant to practical problems. To broaden the utility of pore-scale models for larger domains, we developed a hybrid multiscale model that initially simulates the full domain at the continuum scale and applies a pore-scale model only to areas of high reactivity. Since the location and number of pore-scale model regions in the model varies as the reactions proceed, an adaptive script defines the number and location of pore regions within each continuum iteration and initializes pore-scale simulations from macroscale information. Another script communicates information from the pore-scale simulation results back to the continuum scale. These components provide loose coupling between the pore- and continuum-scale codes into a single hybrid multiscale model implemented within the SWIFT workflow environment. In this paper, we consider an irreversible homogeneous bimolecular reaction (two solutes reacting to form a third solute) in a 2D test problem. This paper is focused on the approach used for multiscale coupling between pore- and continuum-scale models, application to a realistic test problem, and implications of the results for predictive simulation of mixing-controlled reactions in porous media. Our results and analysis demonstrate that the hybrid multiscale method provides a feasible approach for increasing the accuracy of subsurface reactive transport simulations.
•Double diffusion reactant mixing can control mineral precipitation in porous media.•Two chemical systems studied were calcium carbonate and calcium phosphate.•Complex coupling of physical and ...chemical processes controls precipitation outcomes.•Simulations involving reaction kinetics with porosity changes requires high spatial resolution.•Properties of precipitation zones will affect potential engineering applications.
Multi-component mineral precipitation in porous, subsurface environments is challenging to simulate or engineer when in situ reactant mixing is controlled by diffusion. In contrast to well-mixed systems, the conditions that favor mineral precipitation in porous media are distributed along chemical gradients, which evolve spatially due to concurrent mineral precipitation and modification of solute transport in the media. The resulting physical and chemical characteristics of a mixing/precipitation zone are a consequence of coupling between transport and chemical processes, and the distinctive properties of individual chemical systems. We examined the spatial distribution of precipitates formed in “double diffusion” columns for two chemical systems, calcium carbonate and calcium phosphate. Polyacrylamide hydrogel was used as a low permeability, high porosity medium to maximize diffusive mixing and minimize pressure- and density-driven flow between reactant solutions. In the calcium phosphate system, multiple, visually dense and narrow bands of precipitates were observed that were reminiscent of previously reported Liesegang patterns. In the calcium carbonate system, wider precipitation zones characterized by more sparse distributions of precipitates and a more open channel structure were observed. In both cases, formation of precipitates inhibited, but did not necessarily eliminate, continued transport and mixing of the reactants. A reactive transport model with fully implicit coupling between diffusion, chemical speciation and precipitation kinetics, but where explicit details of nucleation processes were neglected, was able to qualitatively simulate properties of the precipitation zones. The results help to illustrate how changes in the physical properties of a precipitation zone depend on coupling between diffusion-controlled reactant mixing and chemistry-specific details of precipitation kinetics.
► We solve reactive transport PDAEs in a parallel fully-coupled implicit manner. ► Physics-based blocking preconditioning strategy is used with JFNK solution method. ► PBP significantly reduce linear ...iteration number for highly nonlinear problems. ► Solution efficiency and parallel scalability of the approach are demonstrated.
Modeling large multicomponent reactive transport systems in porous media is particularly challenging when the governing partial differential algebraic equations (PDAEs) are highly nonlinear and tightly coupled due to complex nonlinear reactions and strong solution-media interactions. Here we present a preconditioned Jacobian-Free Newton-Krylov (JFNK) solution approach to solve the governing PDAEs in a fully coupled and fully implicit manner. A well-known advantage of the JFNK method is that it does not require explicitly computing and storing the Jacobian matrix during Newton nonlinear iterations. Our approach further enhances the JFNK method by utilizing physics-based, block preconditioning and a multigrid algorithm for efficient inversion of the preconditioner. This preconditioning strategy accounts for self- and optionally, cross-coupling between primary variables using diagonal and off-diagonal blocks of an approximate Jacobian, respectively. Numerical results are presented demonstrating the efficiency and massive scalability of the solution strategy for reactive transport problems involving strong solution-mineral interactions and fast kinetics. We found that the physics-based, block preconditioner significantly decreases the number of linear iterations, directly reducing computational cost; and the strongly scalable algebraic multigrid algorithm for approximate inversion of the preconditioner leads to excellent parallel scaling performance.