Volcanology Commons^™

Pre-Eruptive Timescales Via Olivine Diffusion Chronometry For Lava And Tephra Of The Schreibers Meadow Cinder Cone Eruption, Mount Baker (Kulshan), Wa, Sabrina Ansari

All Master's Theses

Mount Baker is a very-high threat Cascade volcano and is located near densely populated areas in northern Washington, putting many at risk if an eruption were to occur. Schreibers Meadow cinder cone is located on the southern flank of Mount Baker and produced the Sulphur Creek (SC) tephra and lava 9.8 ka. Previous work on the SC eruption used petrology and geochemistry to reconstruct the pre-eruptive magmatic architecture beneath Mount Baker Volcanic Field, including olivine-hosted melt inclusions to determine crystallization pressures, as well as plagioclase- and clinopyroxene-derived diffusion timescales to determine magma storage timescales of several months and recharge-to-eruption timescales …

Eruption Initiation Mechanisms And Timescales At A Very-High-Threat Washington Volcano: Clues From Crystal Cargo In Lavas From Mount Baker (Koma Kulshan), Emily Yoder

All Master's Theses

Understanding the timescales of pre-eruptive processes is key for improving forecast models of future eruptions. Mount Baker (Koma Kulshan), WA, is part of the northern Cascade volcanic arc and classified as a very high threat volcano. However, we lack an understanding of the duration of time between the final magmatic process before eruption (potentially the process that initiated the eruption) and the eruption itself. This study focuses on (1) estimating eruption initiation timescales via diffusion chronometry and (2) determining eruption initiation mechanisms for three andesitic lava flows from Mount Baker: Dobbs Creek (~119 ka), Dobbs Cleaver (~105 ka), and Swift …

Timescales Of Magma Storage And The Pre-Eruptive History For The Most Recent Lava Flow At Mount Baker (Koma Kulshan), Wa, Desiree Cunningham

All Master's Theses

Mount Baker, ranked as a very-high threat volcano in the United States, is in a well-populated area in northern WA, putting many at risk in the event of a future eruption. Previous work on Mount Baker focused on understanding the magmatic architecture beneath the volcano and has identified multiple magmatic components that contributed to different eruptions throughout its history. The youngest lava flow in the Mount Baker Volcanic Field is the 9.8 ka Sulphur Creek lava flow. The last study conducted on this flow by Garvey (2022) determined the most recent magma mush configuration that was tapped for this eruption. …

Serial Interaction Of Primitive Magmas With Felsic And Mafic Crust Recorded By Gabbroic Dikes From The Antarctic Extension Of The Karoo Large Igneous Province, Jussi S. Heinonen, Arto V. Luttinen, Frank J. Spera, Saku K. Vuori, Wendy A. Bohrson

Geological Sciences Faculty Scholarship

Two subvertical gabbroic dikes with widths of ~ 350 m (East-Muren) and ≥ 500 m (West-Muren) crosscut continental flood basalts in the Antarctic extension of the ~ 180 Ma Karoo large igneous province (LIP) in Vestfjella, western Dronning Maud Land. The dikes exhibit unusual geochemical profiles; most significantly, initial (at 180 Ma) ε_Nd values increase from the dike interiors towards the hornfelsed wallrock basalts (from − 15.3 to − 7.8 in East-Muren and more gradually from − 9.0 to − 5.5 in West-Muren). In this study, we utilize models of partial melting and energy-constrained assimilation‒fractional crystallization in deciphering the …

The Role Of Cumulate Plagioclase Entrainment In The Steens Basalt, Se Oregon: Insights From The Magma Chamber Simulator, Valerie Strasser

All Master's Theses

The Steens Basalt, located in SE Oregon, is the oldest and most mafic member of the Columbia River Basalt Group (CRBG). Because the Steens Basalt contains large (up to 6 cm) and abundant (up to 40%) plagioclase crystals, it is classified as a Giant Plagioclase basalt (GPB). Understanding the origin of these giant plagioclase yields insight into magmatic processes like recharge, crustal assimilation, cumulate entrainment, and fractional crystallization. Because large and abundant plagioclase are found in other flood basalt provinces worldwide, these results could also enhance understanding of magma evolution and storage conditions for GPB worldwide. The Steens basalt is …

Deep Open Storage And Shallow Closed Transport System For A Continental Flood Basalt Sequence Revealed With Magma Chamber Simulator, Jussi S. Heinonen, Arto V. Luttinen, Frank J. Spera, Wendy A. Bohrson

All Faculty Scholarship for the College of the Sciences

The Magma Chamber Simulator (MCS) quantitatively models the phase equilibria, mineral chemistry, major and trace elements, and radiogenic isotopes in a multicomponent–multiphase magma + wallrock + recharge system by minimization or maximization of the appropriate thermodynamic potential for the given process. In this study, we utilize MCS to decipher the differentiation history of a continental flood basalt sequence from the Antarctic portion of the ~ 180 Ma Karoo large igneous province. Typical of many flood basalts, this suite exhibits geochemical evidence (e.g., negative initial εNd) of interaction with crustal materials. We show that isobaric assimilation-fractional crystallization models fail to …

Quantifying Magma Recharge At Parinacota Volcano, Northern Chile (52 Ka - Recent): Insight From The Magma Chamber Simulator, Jennifer Mcleod

All Master's Theses

The Magma Chamber Simulator (MCS; Bohrson et al., 2014) was used to produce new thermodynamic models that constrain and quantify recharge and fractional crystallization (RFC) processes for the Old Cone (OC; 52-20 ka; 56-66 wt.% SiO₂) and Healing Flows (HF; >8 ka-recent; 58-64 wt.% SiO₂) eruptive stages of Parinacota volcano, northern Chile. Using the MCS, over 200 RFC models were produced that constrain host magma initial conditions, including pressure, fO₂ and H₂O, and host and recharge magma masses, temperatures and compositions. A Rhyolite Domes (RD; 47-40 ka; 74 wt.% SiO₂) …

Giant Plagioclase In The Steens Basalt, Se Oregon: Cumulate Entrainment Revealed By Textural And In Situ Chemical Analysis, Conner H. Toth

All Master's Theses

Many of the lava flows comprising the Steens Basalt in SE Oregon, the oldest and most mafic formation of the Columbia River Basalt Group (CRBG), contain exceptionally large (up to 5 cm) and abundant (up to 40 modal %) plagioclase crystals. This thesis presents a model for giant plagioclase formation in a complex system that integrates petrographic and in situ plagioclase compositional analyses with whole-rock data and computational modeling to explain systematic textural and chemical differences between plagioclase of the lower and upper Steens Basalt sub-sections. Lower Steens plagioclase are typically more chemically homogenous than upper Steens plagioclase and exhibit …

Quantifying Crustal Assimilation In Historical To Recent (1329-2005) Lavas At Mt. Etna, Italy: Insights From Thermodynamic Modeling, Marie Takach

All Master's Theses

The nearly continuous volcanic eruption record at Mt. Etna dating back approximately 700 years provides an excellent opportunity to investigate the geochemical evolution of a highly active volcano. Of particular interest is elucidating the cause of a selective enrichment in alkali elements (K, Rb, Cs) and ⁸⁷Sr/⁸⁶Sr. This alkali enrichment trend, which began in the 17^th century and accelerated after 1971, was accompanied by an increase in the volume, frequency, and explosivity of eruptions. To explain this enrichment, two major arguments are invoked: (1) crustal contributions (e.g., assimilation of the sedimentary basement), and (2) changes in …

Actively Forming Kuroko-Type Volcanic-Hosted Massive Sulfide (Vhms) Mineralization At Iheya North, Okinawa Trough, Japan, Christopher J. Yeats, Steven P. Hollis, Angela Halfpenny, Juan-Carlos Corona, Crystal Laflamme, Gordon Southam, Marco Fiorentini, Richard J. Herrington, John Spratt

All Faculty Scholarship for the College of the Sciences

Modern seafloor hydrothermal systems provide important insights into the formation and discovery of ancient volcanic-hosted massive sulfide (VHMS) deposits. In 2010, Integrated Ocean Drilling Program (IODP) Expedition 331 drilled five sites in the Iheya North hydrothermal field in the middle Okinawa Trough back-arc basin, Japan. Hydrothermal alteration and sulfide mineralization is hosted in a geologically complex, mixed sequence of coarse pumiceous volcaniclastic and fine hemipelagic sediments, overlying a dacitic to rhyolitic volcanic substrate. At site C0016, located adjacent to the foot of the actively venting North Big Chimney massive sulfide mound, massive sphalerite-(pyrite-chalcopyrite ± galena)-rich sulfides were intersected (to 30.2% …

Documenting Mantle And Crustal Contributions To Flood Basalt Magmatism Via Computational Modeling Of The Steens Basalt, Southeast Oregon, Sylvana J. Bendaña

All Master's Theses

Flood basalts are enormous volcanic events with volumes of volcanic cover and intrusive equivalents that are affected by and significantly affect the crust. Steens Basalt represents 31,800 km³ of flood basalt lavas that erupted in eastern Oregon ~16.8 Ma in less than 300,000 years. Analytical data of flows from a 1 km vertical exposure at Steens Mtn. documents time-transgressive changes in composition of two geochemically distinct units: (1) lower Steens, MgO-rich lavas with lower incompatible trace element concentrations and ⁸⁷Sr/⁸⁶Sr, and (2) upper Steens, MgO-poor, with higher incompatible trace element concentrations and ⁸⁷Sr/⁸⁶Sr. …

Evolution Of A Flood Basalt Crustal Magmatic System: In Situ Mineral Data And Computational Modeling Of The Steens Basalt, Megan Graubard

All Master's Theses

Flood basalts are gigantic basaltic eruptions that modify the mass of the crust via intrusion of mantle-derived magma, and change its composition through interaction between magma and crust. The Steens Basalt, located in southeast Oregon, erupted approximately 16.8 million years ago, and is the oldest member of the Columbia River Basalt Group, the youngest and best-exposed flood basalt province on Earth. The Steens Basalt has an eruption volume of approximately 31,800 km³ and the duration of eruption is estimated to be between ~50,000 and 300,000 years. Major- and trace element whole rock data from 111 stratigraphically controlled samples from …

Mesoamerican Sacred Geography And Archeological Landscapes: A Case Study In The Volcanic Highlands Of Michoacan, Mexico, Jarod Scott Drozdowski

All Master's Theses

US and Mexican teams have investigated the La Alberca Structure Complex with the community of Parangaricutiro. Their forest reserve, located near San Juan Nuevo in Michoacan, Mexico, contains several volcanic cones, including the famous Paricutin cinder cone. Sixteen earthen-rock structures have been mapped using aerial photography, GPS, GIS, and total station surveys. The structures range from 300 to 2,400 cubic meters in volume. The larger mound-like structures are traditionally called yacatas. Early radiocarbon dates for two of the three central structures (ca. 5200-2000 BCE) bracket dates for the Late Archaic burial in the La Alberca caldera. Arrangements of the larger …

Timing And Source Of Alkali Enrichment At Mt. Etna, Sicily: Constraints From Clinopyroxene Geobarometry And In Situ Sr Isotope Data, Kaitlyn Nelson

All Master's Theses

Since 1971, Mt. Etna, Europe’s largest and most active volcano, has exhibited increased eruption frequency and explosivity. Associated with this increased activity, researchers have documented higher abundances of alkali elements (K, Rb, Cs) as well as elevated ⁸⁷Sr/⁸⁶Sr in Etnean lavas. The source of this alkali-enrichment has been hotly debated, with end-member hypotheses involving mantle vs. crustal contributions. To further characterize the evolution of the subvolcanic magma storage and transport system, as well as the timing and source of alkali-enrichment, clinopyroxene from ten samples erupted between 1329 and 2004 was targeted for in situ textural, major element …

Thermodynamic Model For Energy-Constrained Open-System Evolution Of Crustal Magma Bodies Undergoing Simultaneous Recharge, Assimilation And Crystallization: The Magma Chamber Simulator, Wendy A. Bohrson, Frank J. Spera, Mark S. Ghiorso, Guy Brown, Jeffrey Creamer, Aaron Mayfield

All Faculty Scholarship for the College of the Sciences

The Magma Chamber Simulator quantifies the impact of simultaneous recharge, assimilation and crystallization through mass and enthalpy balance in a multicomponent–multiphase (melt + solids ± fluid) composite system. As a rigorous thermodynamic model, the Magma Chamber Simulator computes phase equilibria and geochemical evolution self-consistently in resident magma, recharge magma and wallrock, all of which are connected by specified thermodynamic boundaries, to model an evolving open-system magma body. In a simulation, magma cools from its liquidus temperature, and crystals ± fluid are incrementally fractionated to a separate cumulate reservoir. Enthalpy from cooling, crystallization, and possible magma recharge heats wallrock from its …

Structural Deformation Across The Southwest Mina Deflection, California-Nevada: Field Studies In The Huntoon Spring Quadrangle, Eliya R. Hogan

All Master's Theses

New geologic data provides insight into Pliocene strain rates across the southwest Mina deflection. ⁴⁰ Ar/³⁹ Ar geochronology of geologic units indicate cooling of Miocene ignimbrite at 11.86 ± 0.01 Ma, Pliocene basalt lavas at 3.7241 ± 0.0051 and 3.5262 ± 0.0172 Ma, and basaltic andesite and andesite lavas at 3.4715 ± 0.0107 Ma and 3.1078 ± 0.0088 Ma. Pliocene cinder cones and Quaternary sediments are also exposed in the Huntoon Spring quadrangle. Volcanic units are cut by a dominant set of NNE-, NE-, and ENE- striking sinistral faults and a lesser set of NNW- and NW- striking normal …

Bulk Rock Composition And Geochemistry Of Olivine-Hosted Melt Inclusions In The Grey Porri Tuff And Selected Lavas Of The Monte Dei Porri Volcano, Salina, Aeolian Islands, Southern Italy., Angela L. Doherty, Robert J. Bodnar, Benedetto De Vivo, Wendy A. Bohrson, Harvey E. Belkin, Antonia Messina, Robert J. Tracy

All Faculty Scholarship for the College of the Sciences

The Aeolian Islands are an arcuate chain of submarine seamounts and volcanic islands, lying just north of Sicily in southern Italy. The second largest of the islands, Salina, exhibits a wide range of compositional variation in its erupted products, from basaltic lavas to rhyolitic pumice. The Monte dei Porri eruptions occurred between 60 ka and 30 ka, following a period of approximately 60,000 years of repose. The bulk rock composition of the Monte dei Porri products range from basaltic-andesite scoria to andesitic pumice in the Grey Porri Tuff (GPT), with the Monte dei Porri lavas having basaltic-andesite compositions. The typical …

Documenting Magma Evolution Of The Fossa Delle Felci (Salina Island, South Tyrrhenian Sea) By Integrating In Situ Plagioclase Data With Quantitative Modeling, Aaron Mayfield

All Master's Theses

Compositional diversity occurs via processes dominated by recharge, assimilation and fractional crystallization (RAFC), quantification of which can provide insights into eruption catalysts, and therefore assist in hazard mitigation. This study evaluates the relative roles of RAFC processes in Fossa delle Felci (Salina Island, Italy) magmas, which span the basalt to dacite range (52-65 wt% SiO2). Integration of petrographic, whole rock, and in situ plagioclase data with quantitative magma chamber models suggests differentiation dominated by fractional crystallization and magma recharge/mixing in reservoirs located at different crustal levels. At ~12 km depth, magma evolved through crystallization of a pyroxene dominated assemblage; plagioclase …

Petrogenetic Relationship Of The Postcaldera Eruptions Of Mount Mazama, Crater Lake, Oregon; Evolution Of A Sub-Volcanic Magma Chamber Following A Large Silicic Eruption, Michelle Leanna Tebbe

All Master's Theses

Mount Mazama is the volcanic edifice that cataclysmically erupted ~503 km of relatively homogeneous rhyodacite lava ~ 7,700 years ago, forming the caldera known as Crater Lake. Within a few hundred years, andesitic eruptions built three distinct volcanic edifices on the floor of Crater Lake; ~ 3000 years later, rhyodacite eruptions formed a dome (Bacon et al., 2002). How magmatic systems evolve following a shallow, relatively large silicic eruption is the focus of this study. In situ geochemical analysis coupled with high-resolution textural images of plagioclase crystals in the four postcaldera volcanic edifices were used to identify distinct crystal populations …

Documenting Magnatic Processes At Filicudi Island, Aeolian Arc, Italy: Integrating Quantitative Modeling And Plagioclase Textural And In Situ Compositional Data, Michelle Harris

All Master's Theses

Documenting the physiochemical processes that influence magma composition is critical for forecasting eruption styles and managing volcanic hazards. Compositional diversity of magmas develops through recharge, assimilation, and fractional crystallization (RAFC) within subvolcanic magma reservoirs. Integration of MELTS modeling, whole rock, plagioclase textural and in situ elemental and isotopic data from Filicudi Island, Italy allow documentation of the roles and relative chronology that RAFC played in the magmatic evolution and elucidates aspects of the magma plumbing system structure.

Results indicate a polybaric magma plumbing system with deeper (3.5-4 kilobars) and shallower (0.5-1.2 kilobars) storage regions. Within the deeper system, FC acted …

Documenting The Origin Of Compositional Diversity Of Subduction Zone Magmatism, Alicudi, Aeolian Arc (Southern Tyrrhenian Sea, Italy) Using In Situ Plagioclase Data, Rachel Hunt

All Master's Theses

Processes that generate magma compositional diversity are important to document because composition affects degree of explosivity, which impacts hazard mitigation. Magnesium to silicon rich magmas are produced by processes such as magma recharge, assimilation, and fractional crystallization (RAFC) that occur in subvolcanic magma chamber(s). This study evaluates how magma chamber processes contribute to compositional diversity at Alicudi Volcano, Italy. Analytical and petrographic data from ten samples that span the subaerial history include whole rock major and trace elements and strontium/neodymium isotopes, and plagioclase textural types, major and trace elements, and strontium isotopes; numerical modeling was also conducted. Integration of these …

Crustal Assimilation And Magma Recharge In The Recent Mt. Etna Magma Plumbing System: Evidence From In Situ Plagioclase Textural And Compositional Data, Brad Pitcher

Undergraduate Honors Theses

Mt. Etna is Europe's largest and most active volcano, and as a result of its proximity to populated areas, understanding the structure of its magma plumbing system and the nature of its magmatic processes is essential for better predicting eruptive hazards. The aim of this study is to document core to rim textural, chemical, and isotopic variations in plagioclase, in order to investigate the structure of the subvolcanic magma system and the processes by which magmas change composition. Nomarski Differential Interference Contrast (NDIC) imaging was used to characterize the complex textures of plagioclase crystals in six trachybasaltic samples from eruption …

Lacustrine Sediment Record Of Multiple Quaternary Lava Dams On The Owyhee River, Southeastern Oregon, Caitlin Anne Orem

All Master's Theses

Multiple lava dams and correlating lakes impacted the Quaternary evolution of the Owyhee River. Sediment records from lava-dammed lakes were investigated to understand effects of the West Crater (WC) lava dam (~70 ka), the Saddle Butte 2 lava dam (~144 ka), and the Bogus Rim lava dam (~1.9 Ma). Evidence from the WC lava dam and related features indicates that dam duration consisted of five stages (1) dam and lake formation at ~70 ka; (2) dam overflow and lake sedimentation from ~70–46 ka; (3) removal of lava dam and lake termination from ~46 ka to at least 36 ka; (4) …

Paleomagnetic Behavior Of Volcanic Rocks From Isla Socorro, Mexico, Elise Sbarbori, Lisa Tauxe, Avto Goguitchaichvili, Jaime Urrutia-Fucugauchi, Wendy A. Bohrson

All Faculty Scholarship for the College of the Sciences

The direction and magnitude of the geomagnetic field vary both spatially and temporally and undergo significant departures from that of a geocentric axial dipole. In order to properly characterize persistent behaviors, time-averaged field models must be based on the highest quality data. Here we present full-vector paleomagnetic data for volcanic units exposed in the southeast quadrant of the island of Socorro, Mexico. We carried out a joint expedition between the Scripps Institution of Oceanography and the Universidad Nacional Autónoma México to Isla Socorro in January of 2005 during which we collected oriented paleomagnetic samples from 21 sites, representing as many …

Petrogenetic Processes Characterizing The Mount Bachelor, Oregon Magmatic System: Open- Versus Closed-System Processes, Sara Elizabeth Johnson

All Master's Theses

Mount Bachelor volcanic chain (MBVC), located in central Oregon, is one of the larger basalt and basaltic-andesite edifices in central Oregon. Preliminary studies have defined how eruptions have changed composition with time, but a detailed assessment of magma chamber processes has not been conducted. To gain a more thorough understanding of the magmatic processes that have contributed to the observed compositional evolution, this study focuses on one of four eruptive episodes, specifically episode III, the most voluminous episode of the four. Magmatic processes are assessed by focusing on both whole-rock and single crystal data.

The geochemistry and textural diversity of …

Energy-Constrained Recharge, Assimilation, And Fractional Crystallization (Ec-Raxfc): A Visual Basic Computer Code For Calculating Trace Element And Isotope Variations Of Opensystem Magmatic Systems, Wendy A. Bohrson, Frank J. Spera

All Faculty Scholarship for the College of the Sciences

Volcanic and plutonic rocks provide abundant evidence for complex processes that occur in magma storage and transport systems. The fingerprint of these processes, which include fractional crystallization, assimilation, and magma recharge, is captured in petrologic and geochemical characteristics of suites of cogenetic rocks. Quantitatively evaluating the relative contributions of each process requires integration of mass, species, and energy constraints, applied in a self-consistent way. The energy-constrained model Energy-Constrained Recharge, Assimilation, and Fractional Crystallization (EC-RaxFC) tracks the trace element and isotopic evolution of a magmatic system (melt + solids) undergoing simultaneous fractional crystallization, recharge, and assimilation. Mass, thermal, and compositional (trace …

Implications Of Climate Change And Tephra Accumulation On The Lacustrine Environment In A Volcanic Caldera In Michoacan, Mexico, Tabitha Trosper

Undergraduate Honors Theses

La Alberca is a volcanic caldera in the Michoacan-Guanajuato volcanic field of central Mexico. The oldest known burial ( ~5000 yrs BP) in the state of Michoacan was discovered beneath prehispanic cliff paintings on the caldera wall. The stratigraphy of fluvial and lacustrine sediments exposed in trenches in the caldera floor record environmental changes that have implications for the type of human use of the site. La Alberca is a closed system with an alluvial fan that extends from the northeastern wall into the lowest southern comer of the caldera floor. Paricutin eruption in 1943-52 and other earlier eruptions deposited …

Fluvial Response To Intra-Canyon Lava Flows, Owyhee River, Southeastern Oregon, Cooper Cooke Brossy

All Master's Theses

At least six lava flows have entered the Owyhee River Canyon north of Rome, Oregon, since the Pliocene and directly impacted the Owyhee River. The effects on the river of the two youngest lava flows, the West Crater (60–80 ka) and Saddle Butte (> 60–90 ka), are readily apparent. These two lava flows entered a paleo-Owyhee Canyon several kilometers wide via three different tributary drainages. The flows dammed the Owyhee River, created lakes, and effectively confined the river to the opposite side of the valley from the flows’ entrance. Lava from these flows filled a paleo-Owyhee Canyon to depths of …

Deciphering The Signature Of Magma Mixing: Examples From The Castle Creek Eruptive Period, Mount St. Helens, Washington, Seth Taylor Mattos

All Master's Theses

Mount St. Helens (MSH) volcano in southwestern Washington has intermittently erupted dacitic products for the last 40,000 years. On limited occasions, the volcano has produced andesite lava flows, and during one short-lived period, basaltic lava flows. This time interval has been termed the Castle Creek eruptive period and occurred between approximately 2500 and 1700 years B.P. The Castle Creek period erupted dacite, andesite and basalt within this short span of time. Andesite and dacite eruptions dominate the first approximately 700 years of the period, and all basaltic units were erupted in approximately the last 100 years of the period. This …

Energy‐Constrained Open‐System Magmatic Processes Iv: Geochemical, Thermal And Mass Consequences Of Energy‐Constrained Recharge, Assimilation And Fractional Crystallization (Ec‐Rafc), Wendy A. Bohrson, Frank J. Spera

All Faculty Scholarship for the College of the Sciences

A wealth of geochemical and petrological data provide evidence that the processes of fractional crystallization, assimilation, and magma recharge (replenishment) dominate the chemical signatures of many terrestrial igneous rocks. Previous work [Spera and Bohrson, 2001; Bohrson and Spera, 2001] has established the importance of integrating energy, species and mass conservation into simulations of complex magma chamber processes. An extended version of the energy‐constrained formulation, Energy‐Constrained Recharge, Assimilation, Fractional Crystallization (EC‐RAFC), tracks mass and compositional variations of melt, cumulates, and enclaves in a magma body undergoing simultaneous recharge, assimilation, and fractional crystallization [Spera and Bohrson, 2002]. …