How To Calculate Fractional Crystallization Triangle Diagram

Calculate initial and evolved liquid compositions, then plot them on a ternary-style triangle projection for fast petrologic interpretation.

How to Calculate a Fractional Crystallization Triangle Diagram: Expert Step-by-Step Guide

Fractional crystallization is one of the most important controls on magma evolution. In practical petrology, geochemistry, and volcanic hazard interpretation, you often need to answer a simple but technically rich question: if crystals are removed from a melt as cooling proceeds, how does the liquid composition move through compositional space? A triangle diagram (also called a ternary diagram) gives a fast visual answer when composition is expressed in three components that sum to 100%.

This guide shows exactly how to calculate fractional crystallization numerically and then convert those numbers into a triangle plot. The calculator above automates the arithmetic, but knowing the equations and assumptions is critical if you want defensible scientific interpretation. We will cover the logic, equations, worked sequence, data interpretation, common errors, and how to validate your result using real geochemical behavior.

1) What a fractional crystallization triangle diagram represents

A triangle diagram represents three compositional components, commonly major components (for example, Qz-Fs-Mafic proxies in simplified teaching systems) or selected geochemical groups that sum to 100%. Each point inside the triangle corresponds to one normalized composition. If a melt evolves from an initial composition to a more fractionated composition, the point migrates through the triangle.

• Initial liquid point: normalized starting composition of the melt.
• Evolved liquid point: composition after some crystals have been removed.
• Optional crystal/cumulate point: bulk composition of solids extracted during the interval.

The path shape depends on the mineral assemblage and partitioning. Compatible components (high bulk D) tend to be removed into solids, making the residual melt depleted. Incompatible components (low D) remain in the liquid and become enriched as F decreases.

2) Core equations used in the calculator

You first normalize the initial components so they sum to 1. Then you apply a crystallization model. The two most common are Rayleigh fractional crystallization and batch (equilibrium) crystallization:

1. Rayleigh fractional crystallization: C_l = C_0 * F^(D-1)
2. Batch crystallization: C_l = C_0 / (F + D*(1-F))

Where:

• C_0 is initial concentration of a component
• C_l is concentration in the residual liquid
• F is liquid fraction remaining (1.0 means no crystallization, 0.3 means 70% solid removed)
• D is bulk partition coefficient for the component

After computing C_l for A, B, and C, renormalize the three values to 100% for plotting in a triangle. The calculator also estimates extracted solids by mass balance: C_s = (C_0 - F*C_l)/(1-F) when F < 1.

3) Step-by-step manual workflow

1. Choose three components that are meaningful for your petrologic question.
2. Enter measured or modeled initial composition values.
3. Assign realistic bulk partition coefficients for each component.
4. Set F (liquid fraction remaining) based on your crystallization scenario.
5. Compute residual liquid with Rayleigh or batch equation.
6. Renormalize to 100% so the composition is triangle-ready.
7. Convert to XY coordinates for plotting in a triangular frame.
8. Interpret direction and magnitude of movement relative to phase compatibility.

4) How to choose good D values

The largest source of modeling error is unrealistic partition coefficients. Bulk D depends on mineral assemblage proportions and mineral-melt partitioning for each phase. In mafic systems dominated by olivine + clinopyroxene + plagioclase, components that are hosted in mafic silicates often have higher D than highly incompatible lithophile components.

Practical tip: run sensitivity tests by varying each D by plus/minus 25% to see if your triangle trend is robust. If trajectory changes strongly, your conclusion should be presented as conditional, not absolute.

5) Comparison table: typical partition behavior and enrichment at F = 0.40 (Rayleigh)

Element or Component Proxy	Typical Bulk D Range in Basaltic Systems	Midpoint D Used	Rayleigh Factor F^(D-1) at F=0.40	Residual Liquid Trend
Rb (incompatible)	0.01-0.20	0.10	2.28	Strong enrichment
Zr (incompatible to mildly compatible)	0.10-0.30	0.20	2.08	Enrichment
Sr (often plagioclase-sensitive)	1.00-3.00	2.00	0.40	Depletion
Ni (olivine-compatible)	5.00-20.00	10.00	0.00026	Very strong depletion
Cr (spinel or pyroxene-compatible)	10.00-50.00	20.00	0.000000027	Extreme depletion

These magnitudes illustrate why incompatible components dominate evolved liquids while compatible components collapse rapidly. Your triangle path should reflect this asymmetry. If your model predicts the opposite without a mineralogical reason, revisit D and normalization steps.

6) Real-world compositional context table

A useful reality check is comparing modeled trends to common basalt families. While natural magmas also reflect source heterogeneity, assimilation, and mixing, fractional crystallization often explains first-order movement from primitive to evolved compositions.

Basalt Type (Typical)	SiO2 wt%	MgO wt%	FeO* wt%	Na2O + K2O wt%	Typical Interpretation
MORB (mid-ocean ridge basalt)	49-51	7-10	8-11	2-4	Relatively depleted mantle source, moderate fractionation
OIB (ocean island basalt)	46-51	6-12	8-13	3-6	Enriched source signatures plus fractionation effects
Arc basalt (subduction settings)	50-54	4-8	7-10	3-7	Hydrous processes, plagioclase and pyroxene control, variable oxidation

If your modeled triangle path is intended to represent basalt differentiation, check whether trends are coherent with these broad ranges. For example, strong compatible-element depletion and incompatible-element enrichment at lower F should emerge naturally under Rayleigh behavior.

7) Interpreting the triangle path correctly

• Toward one vertex: means relative enrichment of that component in residual melt after normalization.
• Away from one vertex: means relative depletion due to crystal uptake.
• Short movement: can imply mild crystallization, similar D values, or high F.
• Long movement: often indicates strong incompatible-compatible contrast and low F.

Remember: triangle plots are compositional and normalized. Absolute concentration can rise for one component while still appearing visually constrained if another component enriches faster. Pair triangle interpretation with concentration-vs-F plots when precision matters.

8) Common mistakes and how to avoid them

1. Not normalizing inputs: triangle math requires sum = 100% or 1.0.
2. Using impossible F values: keep 0 < F ≤ 1.
3. Applying trace-element D to major oxides without justification: major elements may need phase-equilibrium constraints.
4. Confusing batch and Rayleigh models: they can diverge strongly at low F.
5. Ignoring mineral assemblage changes: D is not always constant through the whole path.

9) Validation strategy for professional use

For research or advanced project reports, validate your triangle model with at least three checks:

• Mass-balance consistency between initial liquid, residual liquid, and cumulate.
• Reasonable D values based on phase assemblage and literature constraints.
• Agreement with independent indicators such as Mg#, Ni, Cr, Sr anomalies, or petrography.

You can also simulate multiple crystallization steps with evolving D values, which better matches natural systems where crystallizing mineralogy changes with temperature and pressure.

10) Authoritative references for further study

For deeper geochemical context and high-quality educational resources, use these references:

• USGS Geochemistry Program (.gov)
• MIT OpenCourseWare Trace Element Geochemistry (.edu)
• Carleton SERC Fractional Crystallization Teaching Resource (.edu)

Conclusion

Calculating a fractional crystallization triangle diagram is straightforward when you follow a disciplined sequence: normalize initial composition, choose a physically defensible model, apply the proper equation with realistic D and F, renormalize outputs, and plot in ternary space. The calculator on this page handles the numerics and plotting, but the interpretation depends on your geological context. Use the trend direction, magnitude, and compatibility logic together, and always test sensitivity to model assumptions.