Dolomitization Models

Dolomitization is the process by which existing calcium carbonate minerals (calcite or aragonite) are replaced by dolomite [CaMg(CO₃)₂], expressed by the replacement reaction: 2CaCO₃ (solid) + Mg²⁺ (aq) = CaMg(CO₃)₂ (solid) + Ca²⁺ (aq). ^[1] Direct precipitation from aqueous solution — the other chemical pathway — requires temperatures far in excess of normal surface conditions, so most geologists believe that dolomitization by replacement was the dominant mechanism for forming ancient dolomite bodies.

Kinetic Requirements

Theoretical considerations suggest that dolomite formation is favored kinetically by high Mg²⁺/Ca²⁺ ratios, low Ca²⁺/CO₃^2– ratios, and low salinity. ^[1] Higher temperatures are also favorable, and at temperatures exceeding about 100°C most kinetic inhibitors — including Mg²⁺ hydration — become ineffective. ^[1] These requirements explain why researchers have focused heavily on finding natural settings where elevated Mg/Ca ratios and fluid movement can drive the replacement reaction at near-surface conditions.

The Three Principal Models

Three principal models have been proposed to explain how dolomitization occurs in natural environments near the Earth’s surface. ^[1] Models 1 and 3 invoke seawater — concentrated by evaporation in the case of Model 1 — as the dolomitizing fluid, while Model 2 requires mixing of seawater and fresh (meteoric) water. ^[1]

1. Hypersaline (Sabkha / Evaporation-Reflux) Model

In hypersaline environments such as the sabkhas of the Persian Gulf — coastal plains characterised by the presence of evaporites — intense evaporation concentrates Mg in brines, raising the Mg²⁺/Ca²⁺ ratio to levels that favor dolomitization, commonly in association with gypsum and anhydrite formation. ^[1] The mechanism involves evaporative concentration of seawater in a supratidal or lagoonal setting; as CaSO₄ minerals (gypsum, anhydrite) precipitate and remove Ca from solution, the Mg/Ca ratio rises. Dense hypersaline brine then sinks and refluxes downward and seaward through underlying carbonate sediments, driving dolomitization in the subsurface. This explains the frequent spatial association of ancient dolomites with evaporite sequences in the stratigraphic record.

2. Mixed-Water (Mixing-Zone) Model

The mixing-zone model proposes that where meteoric (fresh) groundwater and marine porewaters meet in the subsurface, the resulting mixed fluid is undersaturated with respect to calcium carbonate but supersaturated with respect to dolomite. ^[1] This zone of mixing — typically beneath coastal carbonate platforms — provides the chemical conditions (lower Ca²⁺/CO₃^2– ratios, lower effective salinity) that reduce kinetic inhibitors and allow dolomite to nucleate and replace precursor carbonates. The position of the mixing zone shifts with sea level, so dolomitization can sweep through large volumes of rock as the zone migrates over time.

3. Seawater (Shallow-Subtidal) Model

The seawater model proposes that normal or slightly modified seawater, pumped or flushed through carbonate sediments by tidal or thermal circulation, can serve as the dolomitizing fluid. ^[1] Although modern seawater has a Mg/Ca ratio too high for calcite to precipitate freely but not high enough to overcome all dolomite kinetic barriers, the enormous volumes of seawater that can flow through a permeable carbonate platform over geologic time provide a vast reservoir of Mg ions. Given sufficient time and fluid flux, even normal-salinity seawater may drive dolomitization in shallow-subtidal settings where temperatures are relatively elevated.

Scale of the Problem

The three models each address some of the kinetic requirements for dolomitization, but none of them has been universally accepted as a complete explanation, because none fully accounts for the immense volumes of ancient dolomite preserved in the stratigraphic record. ^[1] The mismatch between the small volumes of modern dolomite being formed today and the thick ancient sequences raises the question of whether early-formed near-surface processes, or burial diagenetic replacement, or some combination of both, were responsible for the bulk of the world’s dolomite.

Hypersaline Model — Detailed Mechanism

The sabkha environments of the Persian Gulf, and the supratidal zones of arid climates more generally, host many of the best-documented modern and Holocene dolomites, and these settings are the main evidence for the hypersaline model. ^[1] Under strongly evaporative conditions — where rates of evaporation exceed rates of precipitation — seawater beneath the sediment surface is concentrated by evaporation, and this concentration process causes aragonite and gypsum to precipitate, preferentially removing Ca²⁺ from the water and raising the Mg/Ca ratio. ^[1] The Mg/Ca ratio in normal seawater is about 5:1; when it rises possibly in excess of 10:1, dolomite is thought to form. ^[1]

Two mechanisms concentrate the brine to these high ratios. One involves evaporative pumping: capillary water in the sediments of the sabkha evaporates at the surface, and upward flow from the saturated groundwater zone continuously replenishes and concentrates the pore fluid. ^[1] The other involves surface evaporation of brine in ponds or bays, which raises the brine density above that of normal seawater so the concentrated brine sinks and flushes downward through underlying carbonate sediment in a process called seepage refluxion, potentially driving dolomitization in the subsurface. ^[1] Despite the clarity of the mechanism in principle, the overall volume of dolomite forming in modern sabkha environments is believed to be relatively small, and controversy persists about whether the dolomite forms by replacement of aragonite or high-magnesian calcite, or by primary precipitation of disordered protodolomite that subsequently achieves better cation ordering. ^[1]

Mixing-Zone Model — Mechanism and Critique

The mixing-zone model holds that brackish groundwaters produced by mixing of seawater with meteoric water can be supersaturated with respect to dolomite at Mg²⁺/Ca²⁺ ratios far lower than those required in hypersaline conditions. ^[1] In coastal areas where meteoric waters meet seawater in the subsurface, mixing is suggested to lower salinities sufficiently so that dolomites can form at Mg²⁺/Ca²⁺ ratios ranging from the normal seawater value of about 5:1 to as low as 1:1 — presumably because less competition from other ions in the less saline water reduces the kinetic barrier to dolomite nucleation. ^[1] The model has also been termed the Dorag model and the schizohaline model by its earlier proponents. ^[1]

The model has attracted serious criticism, however. The original calculations for the Dorag model used solubility values for well-ordered dolomite when the relevant mineral under surface conditions is actually a less ordered, Ca-rich dolomite with higher solubility, making the model’s thermodynamic case less secure than originally presented. ^[1] Furthermore, dolomite does not form in most modern freshwater/seawater mixing zones, and where it does form the volume is small; the true significance of the mixing zone may lie instead in its role of driving fluid movements in the marine groundwater below. ^[1]

Seawater (Shallow-Subtidal) Model — Tidal Pumping

The seawater model proposes that dolomitization can occur in normal, unmodified seawater if a sufficient volume is forced through the sediment so that each pore volume is constantly renewed, continuously supplying new Mg²⁺ while flushing away Ca²⁺ ions and other ions that might poison the dolomite crystal structure. ^[1] A documented example is Sugarloaf Key, Florida, where seawater is forced upward and downward through Holocene carbonate mud by the rise and fall of spring tides — a process called tidal pumping — and dolomite is actively forming in the sediment even though little or no evaporation has occurred. ^[1] Dolomitization by this mechanism involves both precipitation as a cement and later replacement of preexisting crystallites, and similar flushing may also occur during sea-level rise as marine porewaters move landward within a carbonate platform. ^[1]

Sulfate Inhibition of Dolomitization

Experimental work demonstrated that dissolved SO₄^2– inhibits the formation of dolomite, offering a chemical explanation for why dolomite is scarce in open-marine environments where sulfate is abundant. ^[1] Dissolved SO₄^2– ions can allegedly inhibit dolomitization of calcite at SO₄^2– concentrations as low as 5 percent of normal seawater values, meaning that any natural process reducing sulfate — such as bacterial sulfate reduction or precipitation of calcium sulfate (CaSO₄·2H₂O) — would favour dolomite formation. ^[1] Subsequent experiments confirmed that dissolved sulfate at a concentration of 0.005 M retards but does not prevent dolomitization of calcite at high temperatures, with the effect likely related to calcite dissolving faster in a sulfate-free environment due to its greater degree of undersaturation there. ^[1]

Bacterial Dolomite Precipitation

Bacteria can play a role in precipitating dolomite under some conditions, with documented evidence from Lagoa Vermelha — a shallow-water coastal lagoon near Rio de Janeiro, Brazil — where dolomite precipitates at normal Earth-surface temperatures in black, organic-rich sediments through the activities of sulfate-reducing anaerobic bacteria. ^[1] Precipitation occurs because sulfate reduction releases excess Mg along with other by-products, and saturation of Mg at the submicron scale in microenvironments around bacterial cell bodies creates conditions favourable for preferential dolomite precipitation; the precipitate begins as a Ca-rich dolomite that undergoes aging to increase its cation ordering over time. ^[1] Laboratory experiments using sulfate-reducing bacteria cultured from the same lagoon have confirmed the biological mechanism. ^[1]

Climate and Ocean Chemistry Effects on Dolomitization

Because Mg²⁺ concentrations in the ocean were higher during aragonite-sea intervals — when lower rates of seafloor spreading and lower sea levels reduced Mg²⁺ absorption onto hot seafloor basalts — dolomite precipitation may have been favoured during such times compared to calcite-sea intervals. ^[1] Ocean temperature also appears to influence dolomite formation: rapid seafloor spreading raises atmospheric CO₂ through outgassing, producing greenhouse (hothouse) conditions, while slow spreading produces icehouse conditions, and dolomite has been reported to be more common in rocks deposited during greenhouse times. ^[1] The proposed explanation is that higher temperatures increase terrestrial weathering, delivering more dissolved carbon — mainly as bicarbonate (HCO₃^–) — to the ocean; this increased ocean bicarbonate causes greater supersaturation with respect to dolomite than to calcite, favouring dolomite precipitation during warm greenhouse periods. ^[1]

Subsurface (Burial) Dolomitization

Much dolomite in the geological record carries relict textures showing it formed by replacement of a precursor limestone in the subsurface long after deposition — sometimes millions to hundreds of millions of years later — rather than through early penecontemporaneous processes. ^[1] A striking example is the Western Canada Sedimentary Basin, where massive replacement dolomites constitute 50–90 percent of all Devonian dolomites, with dolomitization having occurred in both the intermediate subsurface (500–1,500 m) and the deep subsurface (1,500–3,000 m). ^[1] At the elevated temperatures present in the subsurface, dolomitization can proceed readily in any buried limestone with sufficient porosity and permeability to allow large volumes of Mg-bearing water to circulate through it; the central challenge is identifying a mechanism to drive that circulation. ^[1]

Three fluid-circulation mechanisms have been proposed: (1) gravity-driven flow controlled by the elevation of the meteoric recharge area, creating a hydraulic head; (2) thermal convection driven by a geothermal heat source below the basin; and (3) buoyant circulation within freshwater lenses along the mixing zone with saline waters, where brackish-water discharge at the coast draws compensating saline inflow at depth. ^[1] Each mechanism is capable of delivering the enormous volumes of Mg-rich fluid that the scale of burial dolomitization demands, and in practice more than one may operate simultaneously within a basin over geologic time.

Master UPSC Geology Optional

Ex-ONGC Geologist & Rank Holder

Learn the exact analytical answer-writing patterns needed for UPSC Optional from an AIR 2 & AIR 25 holder.

1-on-1 Personalized Mentorship

Elite Batch (Strictly 10 Seats)

Targeted Strategy for AIR 1-100

Bilingual Conceptual Lectures

Join Us

Offline in Delhi