7.5 Evaporite geochemistry

Formation of evaporites

Evaporites are mineral deposits that form when water evaporates faster than it's replenished, concentrating dissolved salts until they precipitate out of solution. These deposits record information about past climates, ocean chemistry, and basin dynamics, making them valuable archives in sedimentary geochemistry.

This section covers the precipitation sequence, the factors that control brine concentration, and the specific solubility thresholds that govern which minerals form and when.

Evaporation sequence

As water evaporates from a basin, dissolved ions become progressively more concentrated. Minerals precipitate in a predictable order governed by their solubility: the least soluble minerals drop out first, and the most soluble minerals precipitate last.

• Carbonates (calcite, aragonite) precipitate at the earliest stages
• Sulfates (gypsum, anhydrite) form next as concentration continues
• Halite (NaCl) crystallizes from increasingly concentrated brine
• Bittern salts (potassium and magnesium salts like sylvite and carnallite) precipitate in the final stages from the most concentrated brines

This sequence is sometimes called the "Usiglio sequence" after the chemist who first documented it experimentally by evaporating Mediterranean seawater.

Brine concentration factors

Several variables control how quickly and completely a brine concentrates:

• Solar radiation intensity drives the rate of evaporation
• Wind patterns enhance surface evaporation and influence crystal nucleation
• Humidity controls the vapor pressure gradient between the brine surface and atmosphere; lower humidity means faster evaporation
• Basin morphology determines water depth, surface-area-to-volume ratio, and circulation patterns. Shallow, restricted basins concentrate brines more efficiently
• Inflow-outflow balance is the single most important control on long-term salinity evolution. A basin must lose more water to evaporation than it gains from rivers, groundwater, and marine inflow

Mineral precipitation order

The precipitation order is governed by the solubility product ($K_{sp}$) of each mineral phase. Minerals with lower $K_{sp}$ values precipitate at lower brine concentrations.

Approximate concentration factors relative to normal seawater:

• Calcite ($CaCO_3$): precipitates at ~3.5–4× concentration
• Gypsum ($CaSO_4 \cdot 2H_2O$): ~4.5–5× concentration (note: some sources cite ~11× for massive gypsum beds, but initial nucleation begins earlier)
• Halite ($NaCl$): ~10–12× concentration
• Sylvite ($KCl$) and carnallite ($KMgCl_3 \cdot 6H_2O$): ~60–70× concentration
• Bischofite ($MgCl_2 \cdot 6H_2O$): ~90× concentration, representing the most extreme evaporation

In practice, most natural evaporite basins never reach the bittern salt stage because it requires removing over 99% of the original water volume.

Major evaporite minerals

The mineralogy of an evaporite deposit reflects the composition of the parent brine and the degree of evaporation it reached. Each mineral carries distinct physical and chemical properties relevant to both geochemical interpretation and economic use.

Halite characteristics

Halite ($NaCl$) is the most volumetrically abundant evaporite mineral.

• Forms cubic crystals with perfect cleavage in three directions
• Highly soluble in water, with solubility increasing modestly with temperature
• Often contains fluid inclusions, tiny pockets of trapped brine that preserve snapshots of the water chemistry at the time of crystallization. These are among the most valuable geochemical tools in evaporite studies
• Occurs as massive beds, hopper crystals (where edges grow faster than faces), or interbedded with other evaporite and clastic minerals

Gypsum and anhydrite

These two calcium sulfate minerals are closely related but form under different conditions.

• Gypsum ($CaSO_4 \cdot 2H_2O$) contains structural water and forms monoclinic crystals. Common crystal habits include selenite (transparent, tabular), satin spar (fibrous), and desert roses (rosette-shaped aggregates). Swallowtail twinning is a diagnostic feature
• Anhydrite ($CaSO_4$) is the dehydrated equivalent, forming orthorhombic crystals. It's denser and harder than gypsum
• Gypsum converts to anhydrite at burial depths greater than ~600 m to 1 km, depending on the local geothermal gradient and pore fluid chemistry. The reaction is: $CaSO_4 \cdot 2H_2O \rightarrow CaSO_4 + 2H_2O$
• Both minerals form nodules, laminations, or massive beds in evaporite sequences

Potash minerals

Potash minerals precipitate from the most concentrated brines and are economically significant.

• Sylvite ($KCl$) often occurs intergrown with halite, forming the rock type called sylvinite
• Carnallite ($KMgCl_3 \cdot 6H_2O$) is highly hygroscopic and can actually dissolve itself by absorbing atmospheric moisture (deliquescence) in humid conditions
• Polyhalite ($K_2Ca_2Mg(SO_4)_4 \cdot 2H_2O$) is gaining attention as a multi-nutrient fertilizer source because it supplies K, Ca, Mg, and S simultaneously
• These minerals are the primary raw materials for potassium fertilizer (potash) production worldwide

Geochemical indicators

Evaporites preserve geochemical signals that allow reconstruction of ancient brine chemistry, water sources, and environmental conditions. The three main approaches are isotopic analysis, trace element geochemistry, and fluid inclusion studies.

Stable isotope signatures

Different isotope systems record different aspects of the evaporite environment:

• $\delta^{18}O$ and $\delta D$ in hydration water (e.g., in gypsum) reflect the isotopic composition and evaporation history of the parent brine
• $\delta^{34}S$ in gypsum and anhydrite tracks sulfate sources and redox cycling. Marine sulfate has a characteristic $\delta^{34}S$ that has varied through geologic time, so this isotope also helps with age correlation
• $\delta^{13}C$ in associated carbonates records carbon sources and biological productivity
• $^{87}Sr/^{86}Sr$ ratios distinguish marine brines (which track the seawater Sr curve) from continental brines (which reflect local bedrock weathering)
• $\delta^{11}B$ can reconstruct paleo-pH of ancient seawater, since boron isotope fractionation is pH-dependent

Trace element composition

Trace elements substitute into evaporite crystal lattices in amounts that depend on brine chemistry:

• Bromine in halite is one of the most widely used indicators. Br substitutes for Cl in the halite structure, and its concentration increases as evaporation progresses. Primary marine halite typically has Br concentrations of 60–200 ppm, while recycled (dissolved and reprecipitated) halite has much lower Br (<40 ppm)
• Sr in gypsum/anhydrite reflects brine evolution and can distinguish marine from non-marine sources
• Mg/Ca ratios in associated carbonates track brine chemistry changes
• REE patterns indicate terrigenous input and redox conditions in the depositional environment
• Elevated trace metals (Cu, Zn, Pb) may signal hydrothermal contributions to the brine

Fluid inclusion analysis

Fluid inclusions are microscopic pockets of ancient brine trapped within crystals as they grew. They're essentially tiny time capsules.

• Homogenization temperatures (measured by heating inclusions until the vapor bubble disappears) provide constraints on formation temperature
• Microthermometry determines salinity and major ion ratios by measuring freezing-point depression and eutectic temperatures
• Chemical analysis of opened inclusions (using techniques like crush-leach or laser ablation ICP-MS) reveals the actual ionic composition of ancient brines
• Gas compositions within inclusions can indicate the presence of organic matter or microbial activity in the depositional environment

Evaporite depositional environments

The setting in which evaporites form strongly influences their mineralogy, geochemistry, and preservation potential. The key distinctions are between coastal vs. inland settings and marine vs. continental brine sources.

Sabkha vs. playa settings

• Sabkhas are low-relief coastal flats with marine influence and high groundwater tables. Evaporite minerals (typically gypsum, anhydrite, and halite) precipitate within the sediment from capillary evaporation of shallow groundwater. Algal/microbial mats are commonly interbedded with the evaporites
• Playas are inland closed basins with no marine connection and fluctuating water levels. They tend to produce more mineralogically diverse deposits, potentially including borates, nitrates, and lithium-bearing brines in addition to standard evaporite minerals
• Both settings experience capillary evaporation, where water is drawn upward through sediment pores and evaporates at or near the surface, forming efflorescent salt crusts

Marine vs. continental basins

• Marine basins connected to the ocean produce evaporites with relatively predictable compositions because seawater chemistry is globally uniform. These deposits often form thick sequences during sea-level lowstands when basins become partially or fully restricted
• Continental basins have brine chemistries that vary widely depending on local bedrock geology and weathering inputs. They can produce unusual mineral assemblages including alkaline lake deposits (trona, nahcolite) and zeolites
• Brine evolution paths differ fundamentally between the two settings. Marine brines follow the well-known seawater evaporation pathway, while continental brines can evolve along multiple chemical divides depending on initial composition

Ancient vs. modern evaporites

• Ancient evaporites can form sequences hundreds of meters thick, preserved in sedimentary basins. The Messinian Salinity Crisis (~5.96–5.33 Ma) produced up to 1–3 km of evaporites in the Mediterranean when the basin became isolated from the Atlantic
• Modern evaporites serve as analogs for interpreting ancient deposits. The Great Salt Lake, Dead Sea, and coastal sabkhas of the Persian Gulf are actively studied for this purpose
• Comparing modern process observations with ancient rock records is essential for accurate paleoenvironmental reconstruction

Diagenesis of evaporites

After deposition, evaporites undergo significant changes during burial and interaction with pore fluids. These diagenetic processes can dramatically alter the original mineralogy, texture, and geochemical signatures.

Dehydration reactions

The most important dehydration reaction in evaporite diagenesis is the gypsum-to-anhydrite conversion:

$CaSO_4 \cdot 2H_2O \rightarrow CaSO_4 + 2H_2O$

1. This reaction typically occurs at burial depths of ~600 m to 1 km, depending on the geothermal gradient and pore water activity
2. The reaction involves a ~38% volume reduction, which can generate collapse structures and breccias
3. Released water increases pore pressures and can drive fluid migration
4. Carnallite can also dehydrate, breaking down to sylvite and bischofite under certain conditions

Recrystallization processes

Primary depositional textures are frequently destroyed during burial:

• Halite recrystallizes readily, forming coarser crystals that obliterate original features like chevron growth bands and hopper crystal morphologies
• Anhydrite nodules can coalesce into massive beds, losing their original nodular fabric
• Recrystallization generates secondary fluid inclusions along healed fractures. These are distinct from primary inclusions and must be carefully distinguished during analysis, since they record later fluid compositions rather than original brine chemistry
• Fabric-destructive diagenesis can make it difficult to determine original depositional environments

Dissolution and reprecipitation

• Undersaturated fluids flowing through evaporite sequences can dissolve minerals, creating porosity and evaporite karst features
• Dissolved salts may reprecipitate elsewhere, forming secondary evaporite cements or vein fills (e.g., selenite veins in fractures)
• Under sustained pressure, halite deforms plastically and flows, forming salt domes and diapirs. These structures can pierce overlying strata and rise kilometers above their source beds
• Dissolution-reprecipitation cycles can concentrate economically valuable minerals by selectively removing more soluble phases

Economic importance

Evaporite deposits are among the most economically significant sedimentary resources, with applications spanning agriculture, chemical manufacturing, construction, and energy.

Salt deposits

• Halite is used for road de-icing, food preservation, and as feedstock for the chlor-alkali industry (producing chlorine and caustic soda)
• Solution mining of salt domes creates large underground caverns that serve as strategic storage for petroleum reserves and natural gas
• Trona (natural sodium carbonate) from evaporite deposits is the primary source of soda ash for glass manufacturing
• Salt domes create structural traps for hydrocarbons, making them important targets in oil and gas exploration. The cap rock above salt domes can also host native sulfur deposits

Potash resources

• Potassium-bearing evaporite minerals (sylvite, carnallite) are the world's primary source of potash fertilizer, critical for global food production
• Major deposits occur in Saskatchewan (Canada), the Perm Basin (Russia), and Belarus, which together dominate global production
• Extraction uses both conventional underground mining and solution mining techniques
• Polyhalite is gaining commercial interest as a lower-cost, multi-nutrient fertilizer alternative

Sulfate mineral extraction

• Gypsum is mined extensively for construction materials, particularly plasterboard (drywall) and plaster of Paris
• Anhydrite is used in cement production and as a soil conditioner
• Sodium sulfate minerals (mirabilite, thenardite) are extracted for use in detergents and glass manufacturing
• Celestite ($SrSO_4$) from evaporite deposits is the primary commercial source of strontium
• Barite ($BaSO_4$), sometimes associated with evaporites, is used as a weighting agent in drilling fluids

Environmental implications

Evaporite deposits and their dissolution products create distinctive environmental challenges, particularly in arid and semi-arid regions where these deposits are most commonly exposed at or near the surface.

Saline soil formation

• Evaporite dissolution combined with capillary rise of saline groundwater leads to soil salinization in arid regions
• Sodium accumulation degrades soil structure by dispersing clay particles, reducing permeability and making the soil hostile to most crops
• Remediation approaches include leaching with fresh water, chemical amendments (e.g., gypsum to displace sodium), and planting salt-tolerant crops (halophytes)
• Improper irrigation management in evaporite-bearing terrain can cause secondary salinization, turning previously productive land into salt flats

Groundwater salinity issues

• Dissolution of subsurface evaporites can raise total dissolved solids (TDS) in aquifers to levels that make water unsuitable for drinking or irrigation
• Upconing of deeper saline water into freshwater aquifers occurs when pumping rates exceed recharge, drawing the fresh-saline interface upward
• Evaporite karst creates preferential flow paths that can rapidly transmit saline water into otherwise fresh aquifer zones
• Careful monitoring and managed extraction rates are essential in regions underlain by evaporites

Evaporite karst development

Evaporite karst develops through dissolution of gypsum and halite, and it progresses much faster than carbonate karst because these minerals are far more soluble.

• Halite dissolves roughly 7,500 times faster than limestone under comparable conditions
• Sinkholes, subsidence, and collapse features are common hazards in areas underlain by shallow evaporites
• These features pose serious challenges for infrastructure, including roads, buildings, and pipelines
• Geophysical methods (ground-penetrating radar, microgravity surveys, electrical resistivity) are used to detect subsurface voids before construction

Analytical techniques

Studying evaporites requires a combination of mineralogical, geochemical, and spatial analysis methods. Each technique addresses different questions about evaporite composition, formation conditions, and distribution.

X-ray diffraction methods

• Powder XRD is the standard technique for identifying and quantifying evaporite mineral phases, especially useful for fine-grained or mixed assemblages
• Rietveld refinement enables accurate quantitative mineral proportions from XRD patterns
• Synchrotron XRD provides high-resolution structural data for detailed crystallographic studies
• XRD of oriented mounts helps identify clay minerals interbedded with evaporites, which can indicate periods of clastic input during evaporite deposition

Geochemical modeling approaches

Modeling brine evolution requires specialized software because evaporite brines have very high ionic strengths, which means standard activity models (like Debye-Hückel) break down.

• PHREEQC and similar codes use Pitzer equations to accurately calculate activity coefficients in concentrated brines
• Reaction path modeling predicts the sequence and amounts of minerals that precipitate as a brine evaporates, allowing comparison with observed mineral assemblages
• Inverse modeling works backward from measured water chemistry to determine what water-rock interactions must have occurred
• Coupled reactive transport models simulate how fluid flow and chemical reactions interact during evaporite diagenesis

Remote sensing applications

• Satellite imagery maps surface evaporite deposits and monitors seasonal changes in salt pans and playas
• Hyperspectral remote sensing can distinguish specific evaporite minerals based on their diagnostic absorption features in the shortwave infrared
• Synthetic Aperture Radar (SAR) detects surface deformation (subsidence) related to subsurface evaporite dissolution
• Integration of remote sensing data with GIS supports regional-scale resource assessment and hazard mapping

Evaporites in Earth history

The distribution and composition of evaporite deposits through geologic time provide records of ocean chemistry, climate, and tectonic configuration. Because evaporites require specific conditions to form and preserve, their occurrence is unevenly distributed through the stratigraphic record.

Precambrian evaporite records

• Precambrian evaporites are rare, largely because they're highly soluble and have had billions of years to be dissolved and recycled
• Some of the oldest known evaporite evidence comes from the Neoarchean (~2.7 Ga) Steep Rock carbonate platform in Ontario, Canada
• Mesoproterozoic (1.6–1.0 Ga) deposits, including those in the Belt Supergroup of western North America, show improved preservation
• Neoproterozoic evaporites are associated with the breakup of the supercontinent Rodinia and provide clues about ocean chemistry during a time of dramatic environmental change (Snowball Earth events)

Phanerozoic evaporite cycles

Major evaporite deposition events correlate with specific tectonic and sea-level conditions:

• Cambrian–Ordovician: widespread evaporites formed in extensive epicontinental seas on low-relief continents
• Permian: the Zechstein Basin (northern Europe) and Castile Formation (west Texas) represent some of the thickest evaporite accumulations in the geologic record
• Triassic–Jurassic: evaporites associated with rifting of Pangea, preserved in basins like the Newark Supergroup of eastern North America
• Late Miocene: the Messinian Salinity Crisis (~5.96–5.33 Ma) produced massive Mediterranean evaporites when tectonic closure of the Strait of Gibraltar connection isolated the basin from the Atlantic

Paleoclimate interpretations

Evaporites contribute to paleoclimate reconstruction in several ways:

• Their geographic distribution maps ancient arid climate belts, typically located in subtropical high-pressure zones
• Stable isotope compositions ($\delta^{18}O$, $\delta^{34}S$) track changes in global hydrologic and sulfur cycles
• Fluid inclusions in halite provide direct measurements of past seawater temperatures and major ion chemistry
• Evaporite mineralogy (particularly whether $MgSO_4$ or $KCl$ salts dominate the bittern stage) reflects secular changes in ocean $Mg/Ca$ ratios and atmospheric $CO_2$ levels
• These data are most powerful when integrated with other paleoclimate proxies such as oxygen isotopes in foraminifera, paleobotanical records, and climate models