Evolution of porosity in the saltwater–freshwater mixing zone of coastal carbonate aquifers: An alternative modelling approach

Dissolution of calcium carbonate in the saltwater–freshwater mixing zone in carbonate islands and coastal aquifers up to now has been treated by coupling geochemical equilibrium programmes to a reactive-transport model. The result is a complex nonlinearly coupled set of differential transport-advection equations, which need high computational efforts. If dissolution rates of calcite are sufficiently fast, such that one can assume the solution to be in equilibrium with respect to calcite a highly simplified modelling approach can be used. By use of the concept of Phillips in his book: Flow and Reactions in Permeable Rock, Cambridge University Press, New York (1991) and its more general formulation by De Simoni et al. [De Simoni, M., Carrera, J., Sánchez-Vila, X., Guadagnini, A., 2005. A procedure for the solution of multicomponent reactive transport problems, Water Resource Research, 41, W11410.], it is possible to decouple the complex set of equations. To calculate initial changes of porosity in the rock matrix one needs solely to solve the advection-transport equation for salinity s. Current codes on density driven flow such as SEAWAT can be used. To obtain the dissolution capacity of the mixed saltwater–freshwater solution the calcium equilibrium concentration ceq(s) is obtained as a function of salinity by PHREEQC-2. Initial porosity changes can then be calculated by a simple analytical expression of the gradient of the spatial distribution s(x, y) of salinity, the distribution of flow fluxes q(x, y) and the second derivative of the calcium equilibrium concentration ceq(s) with respect to salinity s. This approach does no longer require large dispersivities on the order of meters; but reasonable results are obtained with dispersivity on the scale of the pore size of the porous rock, which correctly describes mixing on the molecular scale.