Process and rate of dedolomitization: Mass transfer and C14 dating in a regional carbonate aquifer

Regional dedolomitization is the major process that controls the chemical character of water in the Mississippian Pahasapa Limestone (Madison equivalent) surrounding the Black Hills, South Dakota and Wyoming. The process of dedolomitization consists of dolomite dissolution and concurrent precipitation of calcite; it is driven by dissolution of gypsum.

Deuterium and oxygen isotopic data from the ground water, coupled with regional potentiometric maps, show that recharge occurs on the western slope of the Black Hills and that the water flows northward and westward toward the Powder River Basin. A significant part flows around the southern end of the Black Hills to replenish the aquifer to the east of the Hills. Depth of flow was inferred from interpretation of the silica geothermometer based on the temperature-dependent solubilities of quartz and chalcedony in water. Chemical effects of warm water in the Pahasapa Limestone include changes in the solubility products of minerals, conversion of gypsum to anhydrite, solution and precipitation of minerals, and increases in the tendency for outgassing of carbon dioxide. Where sulfate reduction is not important, sulfur isotope data show that (1) in the Mississippian aquifer, most of the sulfate is from dissolution of gypsum and (2) some wells and springs have a hydrologic connection with overlying Permian and Pennsylvanian evaporites. Sulfate ion concentration, a progress variable, shows a strong correlation with pH as a result of the combined effects of the dedolomitization reactions.

Mass-balance and mass-transfer calculations were used to adjust ¹⁴C values to determine a range of ground-water flow velocities between 2 and 20 m/yr. These velocities are characteristic of carbonate aquifers. The average rates of dolomite and gypsum dissolution are 1.7 × 10⁻⁴ and 3.4 × 10⁻⁴ mmol/kg of H₂O/yr, respectively. The precipitation of calcite is occurring at the rate of 3.4 × 10⁻⁴ mmol/kg of H₂O/yr. The close agreement among the model results demonstrates that dedolomitization is controlling water-rock interactions in this regional carbonate aquifer system.