Barrovian Metamorphism and Index Mineral Zones

Barrovian metamorphism is the canonical medium-pressure, medium-temperature (MP–MT) regional metamorphic sequence developed within thickened continental crust during collisional orogenesis, operating along apparent geothermal gradients of 20 °C to 30 °C per kilometer (pressures of 0.3 to 1.0 GPa and temperatures of 300 °C to 750 °C). Pioneered by British geologist George Barrow in his landmark 1893 mapping of the Neoproterozoic Dalradian Supergroup metapelites across the southeastern Scottish Highlands (Glen Esk, Angus) and subsequently refined by C. E. Tilley (1925), the Barrovian sequence demonstrated that progressive regional metamorphism of a chemically uniform pelitic protolith (aluminous marine shale or mudstone) produces a predictable spatial succession of six index minerals with increasing thermal grade:

Chlorite → Biotite → Garnet (Almandine) → Staurolite → Kyanite → Sillimanite

In the field, the first appearance of each index mineral across the terrain defines a mapped isograd—a line on a geological map (representing the intersection of a three-dimensional surface of constant metamorphic grade with the Earth's topography) along which a specific prograde dehydration reaction was first crossed. As a pelitic rock progresses upward through the Barrovian zones across the greenschist and amphibolite metamorphic facies, its structural fabric coarsens systematically from fine-grained cleaved slate (chlorite zone) and silky phyllite (chlorite–biotite zones), through porphyroblastic garnet-muscovite schist and staurolite-mica schist, into coarse-grained kyanite–sillimanite schists and migmatitic gneiss.

Lower-Grade Pelitic Zones: Chlorite, Biotite, and Garnet

Pelitic shales enter the Barrovian sequence rich in detrital clay minerals (illite, smectite, kaolinite), detrital quartz (SiO₂), and sodic plagioclase, holding up to 8–12 wt% pore and interlayer H₂O. Across the sub-greenschist to greenschist transition (~300 °C, 0.2–0.4 GPa), clays recrystallize into a fine-grained phyllosilicate matrix of chlorite + phengitic muscovite [KAl₂(AlSi₃O₁₀)(OH)₂] + quartz + albite, defining the Chlorite Zone (slates and lower phyllites).

1. Biotite Isograd (~400–425 °C): As temperature increases, celadonite-rich phengitic muscovite reacts with chlorite (or detrital K-feldspar reacts with chlorite) to nucleate brown pleochroic porphyroblasts of biotite [K(Mg,Fe)₃(AlSi₃O₁₀)(OH)₂]:

Phengitic Muscovite + Chlorite → Biotite + Al-rich Muscovite + Quartz + H₂O

1. Garnet Isograd (~450–500 °C): Near the upper boundary of the greenschist facies, iron-rich chlorite reacts with muscovite and quartz to form dodecahedral porphyroblasts of red almandine-rich garnet [Fe₃Al₂(SiO₄)₃], often stabilized at slightly lower temperatures (~450 °C) if manganese (spessartine component, Mn₃Al₂Si₃O₁₂) is present:

Fe-Chlorite + Muscovite + Quartz → Almandine Garnet + Biotite + H₂O

Rocks in this zone form classic garnet-muscovite schist, in which garnet porphyroblasts frequently preserve curved internal inclusion trails (S_i) of quartz and ilmenite ("snowball garnets") recording syntectonic rotation during ductile shearing.

Medium- to High-Grade Zones: Staurolite, Kyanite, and Sillimanite

Crossing ~520–550 °C at 0.4–0.7 GPa moves the metapelite into the amphibolite facies, where progressive dehydration consumes remaining chlorite and nucleates dense aluminosilicate phases:

1. Staurolite Isograd (~520–560 °C): Matrix chlorite and muscovite (or chloritoid in high-Al, low-K pelites) react to form prismatic, penetration-twinned porphyroblasts of monoclinic staurolite [Fe₂Al₉Si₄O₂₂(OH)₂], yielding staurolite-mica schist:

Chlorite + Muscovite + Almandine → Staurolite + Biotite + Quartz + H₂O

1. Kyanite Isograd (~580–620 °C at P > 0.5–0.6 GPa): Once matrix chlorite is exhausted in normal pelites, or via the terminal breakdown of staurolite + muscovite + quartz, triclinic bladed kyanite (the high-pressure polymorph of Al₂SiO₅) crystallizes:

Staurolite + Muscovite + Quartz → Kyanite + Biotite + Almandine Garnet + H₂O

1. Sillimanite Isograd (~640–700 °C): At the highest Barrovian grade, kyanite crosses the univariant polymorphic phase boundary (d*P*/d*T* ≈ +1.8 MPa/K) into fibrolitic or prismatic sillimanite (Al₂SiO₅). At still higher temperatures (~680–730 °C, the "Second Sillimanite Isograd"), muscovite undergoes dehydration melting or breakdown in the presence of quartz:

Muscovite + Quartz → Sillimanite + K-Feldspar (KAlSi₃O₈) + H₂O (or Granitic Melt) converting the rock into a banded sillimanite–K-feldspar gneiss or stromatic migmatite.

Because pelitic schists comprise at least six major chemical oxide components—SiO₂, Al₂O₃, FeO, MgO, K₂O, and H₂O (the KFMASH system), plus minor Na₂O, CaO, MnO, and TiO₂—plotting phase equilibria directly on a compositional diagram requires reducing the dimensionality of the system. In 1957, James B. Thompson Jr. introduced the AFM projection, which projects KFMASH mineral compositions onto a ternary plane with apices A (Al₂O₃), F (FeO), and M (MgO) from the ubiquitous coexisting phases quartz (SiO₂), muscovite (KAl₃Si₃O₁₀(OH)₂), and a pure H₂O fluid. The vertical plotting coordinate A and horizontal coordinate M/(F + M) are defined in molar oxide proportions as:

A = (Al₂O₃ − 3 K₂O) / (Al₂O₃ − 3 K₂O + FeO + MgO)     and     X_Mg = MgO / (FeO + MgO)

On the Thompson AFM diagram, two-phase tie-lines and three-phase triangles reveal why two pelitic layers at the exact same temperature and pressure can display different index minerals: a highly aluminous, iron-rich shale plots above the garnet–chlorite tie-line and forms chloritoid and staurolite early, whereas a less aluminous, magnesium-rich greywacke-shale plots below the garnet–chlorite tie-line and remains a simple biotite–garnet–muscovite schist right up to the kyanite or sillimanite zone without ever nucleating staurolite. Prograde Barrovian isograds correspond to discontinuous topological tie-line flips (such as the breaking of the Almandine–Chlorite tie-line to establish the Staurolite–Biotite tie-line at ~540 °C).

England and Richardson (1977) Thermal Relaxation Model

The physical mechanism generating Barrovian P–T arrays was established quantitatively by Philip England and Stephen Richardson (1977). During continental collision (such as the Cenozoic India–Asia Himalayan collision or the Ordovician Grampian orogeny in Scotland), thrust stacking doubles the thickness of the continental crust from ~35 km to 60–70 km over 5 to 15 million years. Because tectonic burial occurs an order of magnitude faster than conductive heat transport (given a crustal thermal diffusivity κ ≈ 10⁻⁶ m²·s⁻¹, the thermal relaxation time over a thickness L = 35 km is τ ~ L²/κ ≈ 38 million years), buried rocks initially experience a rapid isothermal or low-T pressure increase (P_max), followed by slow radiogenic self-heating (from decay of ⁴⁰K, ²³²Th, ²³⁵U, and ²³⁸U in the thickened upper/middle crust) and basal conductive heat flux as erosion strips the overburden.

Consequently, every rock parcel traces a clockwise pressure–temperature–time (P–T–t) loop (with pressure plotted increasing upward on the vertical axis): peak pressure (P_max, kyanite stability field) is reached millions of years before peak temperature (T_max, sillimanite stability field) is attained during decompressional exhumation. The mapped Barrovian field gradient ("metamorphic piezothermic array") represents the locus of (P_*T*max, T_max) points frozen in as rocks from progressively deeper structural levels dehydrate and cool during exhumation.

In Akiho Miyashiro's (1961) framework of metamorphic facies series, the medium-P/T Barrovian series (~20–30 °C/km, kyanite → sillimanite) occupies the intermediate position between the low-P/high-T Buchan (or Abukuma) series (>35–50 °C/km, characterized by andalusite and cordierite rather than kyanite and staurolite, driven by mid-crustal granitic plutonism or back-arc extension) and the high-P/low-T Franciscan (or Sanbagawa) subduction series (<10–12 °C/km, characterized by glaucophane blueschist and eclogite).

Why did George Barrow map the sequence in pelitic shales rather than sandstones or limestones?

Hydrous marine shales (pelites) are chemically sensitive "petrological thermometers" because they combine high concentrations of Al₂O₃, K₂O, FeO, MgO, SiO₂, and H₂O, enabling a cascade of discontinuous dehydration reactions every 50–70 °C between 300 °C and 750 °C. Pure quartz sandstones (orthoquartzites) and pure calcitic limestones lack the chemical components needed to nucleate new silicate phases and simply coarsen texturally into quartzite and marble across the entire Barrovian P–T range.

What is the difference between Barrovian and Buchan metamorphism in the Scottish Highlands?

Immediately northeast of Barrow's classic Glen Esk district, in the Buchan region of Aberdeenshire, Dalradian metapelites were metamorphosed at lower pressures (0.2–0.35 GPa, <12 km depth) alongside widespread syn-orogenic gabbroic and granitic intrusions. Instead of the Barrovian kyanite → sillimanite progression, Buchan pelites cross the Al₂SiO₅ phase diagram below the aluminosilicate triple point (~0.42 GPa, 530 °C), producing the low-pressure index mineral sequence Chlorite → Biotite → Cordierite → Andalusite → Sillimanite.

How are "inverted" Barrovian metamorphic zones formed in the Himalaya?

Along the Main Central Thrust (MCT) of the Himalaya, Barrovian zones are structurally inverted: high-grade sillimanite and kyanite gneisses overlie lower-grade staurolite, garnet, and biotite schists upward through the structural section. This inverted Barrovian sequence results from a combination of (1) hot hanging-wall crystalline nappes thrust rapidly over cooler footwall metasediments ("hot-iron" conductive heating + shear heating), and (2) ductile channel-flow extrusion and kilometre-scale recumbent folding of pre-existing right-way-up isograds.

1. Barrow, G. (1893). "On an intrusion of muscovite-biotite gneiss in the south-eastern Highlands of Scotland, and its accompanying metamorphism." Quarterly Journal of the Geological Society of London, 49(1–4), 330–358. https://doi.org/10.1144/GSL.JGS.1893.049.01-04.52
2. Thompson, J. B. (1957). "The graphical analysis of mineral assemblages in pelitic schists." American Mineralogist, 42(11–12), 842–858. https://pubs.geoscienceworld.org/msa/ammin/article-abstract/42/11-12/842/540169/The-graphical-analysis-of-mineral-assemblages-in
3. Miyashiro, A. (1961). "Evolution of metamorphic belts." Journal of Petrology, 2(3), 277–311. https://doi.org/10.1093/petrology/2.3.277
4. England, P. C., & Richardson, S. W. (1977). "The influence of erosion upon the mineral facies of rocks from different metamorphic environments." Journal of the Geological Society, 134(2), 201–213. https://doi.org/10.1144/gsjgs.134.2.0201