Description: The formation of mud volcanoes in the Gulf of Cadiz is closely linked to diapirism in the deep subsurface. The Mercator mud volcano (MMV) is a rare example where diapiric emplacement, in addition to being key for upward fluid migration, is also an important zone for fluid and mineral diagenesis. The most intriguing findings in the near-surface muds of the MMV are extremely high salinities of up to 5.2 M of NaCl from diapiric and evaporitic halite dissolution and the occurrence of authigenic gypsum and anhydrite crystals, both of which have not been observed to date in the Gulf of Cadiz. Employing a thermodynamic model we elucidate how the interplay of temperature pulses, strong salinity gradients, and fluid flow dynamically drive mineral dissolution and re-formation. The strong increase in salinity in the pore fluids has important implications for thermodynamic equilibria by significantly lowering the activity of water, thereby raising the gypsum–anhydrite transition zone from 〉1 km to about 400 m sediment depth at the MMV. This transition is further shifted to immediately below the seafloor during intervals of active mud and fluid expulsion when the MV surface temperature is heated up to at least 30 °C. As a consequence, precipitation of authigenic gypsum near the sediment surface (1–2 mbsf) has been linked to the dissolution of evaporites below the MMV. More precisely, the mechanisms generating supersaturation in the ascending gypsum-saturated MMV fluids are (1) the slow and constant cooling of these fluids along the geothermal gradient during their ascent leading to formation of ubiquitous micro-crystals and (2) the more rapid cooling after a heat pulse or transport from greater and warmer depth during an active mud volcano phase leading to the precipitation of cm-scale gypsum crystals or even fist-size concretions. The MMV fluids approaching the salt diapir from farther below have experienced a genesis similar to those of other mud volcanoes in the Gulf of Cadiz located above deep-rooted faults. These processes include clay mineral dewatering, thermogenic degradation of organic matter and deep high-temperature leaching of terrigenous sediments or continental crust.

Description: We present a transport-reaction model (TRACTION) specifically designed to account for non-ideal transport effects in the presence of thermodynamic (e.g. salinity or temperature) gradients. The model relies on the most fundamental concept of solute diffusion, which states that the chemical potential gradient (Maxwell’s model) rather than the concentration gradient (Fick’s law) is the driving force for diffusion. In turn, this requires accounting for species interactions by applying Pitzer’s method to derive species chemical potentials and Onsager coefficients instead of using the classical diffusion coefficients. Electrical imbalances arising from varying diffusive fluxes in multicomponent systems, like seawater, are avoided by applying an electrostatic gradient as an additional transport contribution. We apply the model to pore water data derived from the seawater mixing zone at the submarine Mercator mud volcano in the Gulf of Cadiz. Two features are particularly striking at this site: (i) Ascending halite-saturated fluids create strong salinity (NaCl) gradients in the seawater mixing zone that result in marked chemical activity, and thus chemical potential gradients. The model predicts strong transport-driven deviations from the mixing profile derived from the commonly used Fick’s diffusion model, and is capable of matching well with the profile shapes observed in the pore water concentration data. Even better agreement to the observed data is achieved when ion pairs are transported separately. (ii) The formation of authigenic gypsum (several wt%) occurs in the surface sediments, which is typically restricted to evaporitic surface processes. Very little is known about the gypsum paragenesis in the subseafloor and we first present possible controls on gypsum solubility, such as pressure, temperature, and salinity (pTS), as well as the common ion and ion pairing effects. Due to leaching of deep diapiric salt, rising fluids of the MMV are saturated with respect to gypsum (as well as celestite and barite). Several processes that could drive these fluids towards gypsum supersaturation and hence precipitation were postulated and numerically quantified. In line with the varied morphology of the observed gypsum crystals, gypsum paragenesis at the MMV is likely a combination of two temperature-related processes. Gypsum solubility increases with increasing temperature, especially in strong electrolyte solutions and the first mechanism involves the cooling of saturated fluids along the geothermal gradient during their ascent. Secondly, local temperature changes, i.e. cooling during the transition from MMV activity towards dormancy results in the cyclic build-up of gypsum. The model showed that the interpretation of field data can be majorly misguided when ignoring non-ideal effects in extreme diagenetic settings. While at first glance the pore water profiles at the Mercator mud volcano would indicate strong reactive influences in the seawater mixing zone, our model shows that the observed species distributions are in fact primarily transport-controlled. The model results for SO4 are particularly intriguing, as SO4 is shown to diffuse into the sediment along its increasing (!) concentration gradient. Also, a pronounced gypsum saturation peak can be observed in the seawater mixing zone. This peak is not related to the dissolution of gypsum but is simply a result of the non-ideal transport forces acting on the activity profile of SO4 and Ca profiles.

Description: Gypsum makes up about one fifth of giant salt deposits formed by evaporation of seawater throughout Earth’s history. Although thermodynamic calculations and precipitation experiments predict that gypsum precipitates when the salinity of evaporating seawater attains about 110 g kg-1, gypsum deposits of the Mediterranean Salt Giant often bear the geochemical signature of precipitation from less saline water masses. Addressing this geochemical riddle is important because marine gypsum deposition and continental gypsum erosion affect the global carbon cycle. We investigated gypsum deposits formed in the marginal basins of the Mediterranean Sea during the Messinian Salinity Crisis (about 6 million years ago). These often bear low-salinity fluid inclusions and isotopically light crystallization water, confirming previous published reports that the Mediterranean Salt Giant harbors low-salinity gypsum deposits. A geochemical model constrained by fluid inclusion salinity and isotope (87Sr/86Sr, δ34SSO4, δ18OH2O, δDH2O) measurements excludes that Ca2+- and SO42--enriched continental runoff alone provides the trigger for gypsum precipitation at low salinity. We propose that, concurrent with the prevalent evaporative conditions and with Ca2+- and SO42--bearing runoff, the biogeochemical sulfur cycle is capable of producing a spatially-restricted and temporally-transient increase of Ca2+ and SO42- within benthic microbial mats, creating local chemical conditions conductive to gypsum precipitation. This hypothesis is supported by the presence of dense packages of fossils of colorless sulfur bacteria within gypsum in several Mediterranean marginal basins, together with independent geochemical and petrographic evidence for an active biogeochemical sulfur cycle in the same basins. Should this scenario be confirmed, it would expand the range of environments that promote marine gypsum deposition; it would also imply that an additional, biological coupling between the calcium, sulfur and carbon cycles exists.