Reactive Transport Modelling of Elevated Dissolved Sulphide Concentrations in Sedimentary Basin Rocks

Groundwater with total dissolved sulphide concentrations in excess of $1.0\times {10}^{-4}\text{mol}{\text{L}}^{-1}\left(3\text{mg}{\text{L}}^{-1}\right)$ is relatively common at intermediate depths in sedimentary basins. However, the mechanisms responsible for the formation and spatial distribution of these sulphidic waters in sedimentary basins, which have been affected by periods of glaciation and deglaciation, are not fully understood. Sulphate reduction rates depend on many factors including redox conditions, salinity, temperature, and the presence and abundance of sulphate, organic matter, and sulphate-reducing bacteria. Two-dimensional reactive transport modelling was undertaken to provide potential explanations for the presence and distribution of sulphidic waters in sedimentary basins, partially constrained by field data from the Michigan Basin underlying Southern Ontario, Canada. Simulations were able to generally reproduce the observed depth-dependent distribution of sulphide. Sulphate reduction was most significant at intermediate depths due to anoxic conditions and elevated sulphate concentrations in the presence of organic matter in waters with relatively low salinity. The simulations indicate that glaciation-deglaciation periods increase mixing of waters at this interfacial zone, thereby enhancing rates of sulphate reduction and the formation of sulphide. In addition, the simulations indicate that glaciation-deglaciation cycles do not significantly affect sulphide concentrations in low permeability units, even at shallow depths (e.g., 25 m), while concentrations in permeable units remain stable below depths of 500 m.