Icarus 212 (2011) 24–41
Contents lists available at ScienceDirect
Icarus
journal homepage: www.elsevier.com/locate/icarus
On the chemistry of mantle and magmatic volatiles on Mercury ⇑ Mikhail Yu. Zolotov
School of Earth and Space Exploration, Arizona State University, Tempe, AZ 85287-1404, USA article info abstract
Article history: The surface of Mercury contains ancient volcanic features and signs of pyroclastic activity related to Received 21 November 2009 unknown magmatic volatiles. Here, the nature of possible magmatic volatiles (H, S, C, Cl, and N) is dis- Revised 8 November 2010 cussed in the contexts of formation and evolution of the planet, composition and redox state of its mantle, Accepted 19 December 2010 solubility in silicate melts, chemical mechanisms of magma degassing, and thermochemical equilibria in Available online 25 December 2010 magma and volcanic gases. The low-FeO content in surface materials (<6 wt%) evaluated with remote
sensing methods corresponds to less than 2.3 fO2 log units below the iron–wüstite buffer. These redox Keywords: conditions imply restricted involvement of hydrous species in nebular and accretion processes, and a lim- Mercury, Interior ited loss of S, Cl, and N during formation and evolution of the planet. Reduced conditions correspond to Mercury, Surface Volcanism high solubilities of these elements in magma and do not favor degassing. Major degassing and pyroclastic Cosmochemistry activity would require oxidation of melts in near-surface conditions. Low-pressure oxidation of graphite in moderately oxidized magmas causes formation of low-solubility CO. Decompression of reduced N-sat-
urated melts involves oxidation of nitride melt complexes and could cause N2 degassing. Putative assim-
ilation of oxide (FeO, TiO2, and SiO2) rich rocks in magma chambers could have caused major degassing through oxidation of graphite and S-, Cl- and N-bearing melt complexes. However, crustal rocks may never have been oxidized, and sulfides, graphite, chlorides, and nitrides could remain in crystallized
rocks. Chemical equilibrium models show that N2, CO, S2,CS2,S2Cl, Cl, Cl2, and COS could be among the most abundant volcanic gases on Mercury. Though, these speciation models are restricted by uncer- tain redox conditions, unknown volatile content in magma, and the adopted chemical degassing mechanism. Ó 2010 Elsevier Inc. All rights reserved.
1. Introduction: Volcanic features and composition of surface (Robinson and Lucey, 1997; Murchie et al., 2008; Robinson et al., materials 2008; Ernst et al., 2010). The ground mid-infrared emission spectra
suggest silicate surface materials with spatially variable SiO2 con- Morphological and color analyses of Mariner 10 images suggest tent and dominant mafic to ultramafic compositions (Boynton a volcanic origin of some surface features on Mercury (Strom et al., et al., 2007). These interpretations are consistent with high heights 1975; Basaltic Volcanism, 1981; Spudis and Guest, 1988; Robinson of lobate flow fronts indicating high-volume eruptions of mafic and Lucey, 1997). The presence of volcanic rocks is confirmed by lava similar to terrestrial flood-basalts (Wilson and Head, 2008). the initial MESSENGER results (Head et al., 2008, 2009a; Robinson Some regions (e.g., the inner part of Caloris basin) could be repre- et al., 2008; Murchie et al., 2008; Prockter et al., 2010). The volca- sented by evolved igneous rocks with intermediate compositions, nic features demonstrate mostly effusive lunar-like activity with- as suggested from mid-infrared spectra (Sprague et al., 2002, out stratovolcanoes or large sheet structures (Head et al., 2008, 2009; Boynton et al., 2007). Intermediate compositions may repre- 2009a). However, several features could be interpreted as pyroclas- sent rocks formed through upward fractional crystallization of a tic deposits (Head et al., 2008, 2009a; Murchie et al., 2008; Robin- putative magma ocean (Brown and Elkins-Tanton, 2009)or son et al., 2008; Kerber et al., 2009; Prockter et al., 2010). through differentiation of impact melts and basaltic magma in Spectroscopic data in visible, near- and mid-infrared ranges chambers. Although magma chambers are suggested from some indicate the compositional diversity of the Mercury’s surface images (Head et al., 2009b), there are no morphological signs of (e.g., Sprague et al., 2009; Denevi et al., 2009; Vernazza et al., volcanic activity associated with evolved magmas (Head et al., 2010). Suggested volcanic plains and pyroclastic deposits are 2008, 2009a). brighter and redder than underlying rocks and impact ejectae Visible, near- and mid-infrared spectra of Mercury’s surface are consistent with the presence of Mg-rich olivine and pyroxene, Mg- ⇑ Fax: +1 480 965 8102. and Ca-clinopyroxene, plagioclase feldspars, silicate glasses, and E-mail address: [email protected] sulfides (Sprague et al., 1995, 2002, 2009; Warell et al., 2006,
0019-1035/$ - see front matter Ó 2010 Elsevier Inc. All rights reserved. doi:10.1016/j.icarus.2010.12.014 M.Yu. Zolotov / Icarus 212 (2011) 24–41 25
2010; Boynton et al., 2007; Vernazza et al., 2010), though interpre- The identity of the volatiles that caused the pyroclastic erup- tations are model-dependent and surface mineralogy remains tions is enigmatic. Kerber et al. (2009) showed that abundant mag- uncertain. The lack of spectral features in visible and near-infrared matic gases (up to 1.3 wt% depending on molecular mass) are ranges can be accounted for by the presence of low-albedo opaque needed to form the large pyroclastic deposit at Caloris basin. They phase(s), which could represent sulfides and/or Ti-bearing oxides suggested Earth-like volatiles (H2O, CO2, CO, H2S) along with abun- (Rava and Hapke, 1987; Denevi and Robinson, 2008; Denevi dant S gases to explain the eruptions. This paper discusses the ele- et al., 2009; Vernazza et al., 2010). The detection of Mg in the exo- mental abundances of H, S, C, Cl, and N in Mercury’s silicate sphere of Mercury (McClintock et al., 2009) agrees with the identi- interior, the composition of possible magmatic volatiles, and fication of Mg-silicates by ground-based spectroscopy, although chemical aspects of degassing. Then, thermochemical calculations Mg sulfides could be present as well (Sprague et al., 1995). The are used to estimate volcanic gas speciation at appropriate bulk presence of Ca and alkalis in the exosphere (e.g., McClintock compositions, redox conditions (fO2), temperatures (T), and pres- et al., 2009) also agrees with the identification of Ca-pyroxenes sures (P). and feldspars in surface materials. Surface spectra of Mercury in near- and mid-infrared are distinct from lunar spectra, probably 2. Volatiles in the interior indicating a fundamental difference in abundance and composition of ferrous silicates (Vernazza et al., 2010). Observational data on Mercury’s volatiles are limited to exo- Although the abundance and speciation of iron in surface mate- spheric Na, K, and possibly, S- and H-bearing ionized molecules rials remains uncertain, the majority of remote-sensing data indi- formed through sputtering of surface materials (Potter and Mor- cate a low-FeO content. The lack of 1 lm absorption in most of gan, 1985, 1986; Zurbuchen et al., 2008). The abundance and spe- Mercury’s spectra, the color of the surface, and microwave emissiv- ciation of volatiles in the interior would reflect the composition ity data indicate an FeO-depleted composition (<3–6 wt% FeO) and redox state of parent planetesimals, and specifics of accretion, (e.g., McCord and Clark, 1979; Vilas, 1988; Blewett et al., 1997; impact re-processing, differentiation, and initial degassing. The Robinson and Taylor, 2001; Boynton et al., 2007). Mid-infrared pathways of formation and evolution of Mercury remain poorly ground emission spectra do not reveal FeO-rich compositions constrained, although diverse scenarios have been suggested, as (Sprague et al., 2002, 2009; Boynton et al., 2007; Vernazza et al., reviewed by Solomon (2003), Taylor and Scott (2003), Benz et al. 2010). These ground observations are consistent with the interpre- (2008), and Breuer et al. (2008). Most of these models are aimed tations of Mariner 10 color data (Robinson and Lucey, 1997; Denevi at explaining a high core/mantle mass ratio and the observed and Robinson, 2008). Initial MESSENGER results also do not indi- FeO deficiency in surface materials. In this section, we discuss vol- cate a presence of abundant surface ferrosilicates (Solomon et al., atiles in the context of formation and evolution pathways. 2008; McClintock et al., 2008; Robinson et al., 2008). The interpretation of MESSENGER neutron spectrometer data obtained during three flybys indicates elevated (lunar-mare type) 2.1. Bulk composition in the framework of formation and evolution concentrations of Fe and Ti, though phase composition cannot be scenarios obtained (Lawrence et al., 2010). These authors suggest an ilmenite content of 7–18 wt% in surface materials. Lawrence et al. (2010) The high density of Mercury (5.43 g/cm3) implies a high Fe/Si also proposed early crystallization of ilmenite in oxidized lunar- ratio and a significant fractionation of materials that originated type magmas (cf., Usselman and Lofgren, 1976) that may lead to with solar Fe/Si ratio. The fractionation could have been associated late crystallization of FeO-poor pyroxenes. Alternatively, a pres- with nebular condensation under specific conditions, chemical gra- ence of Fe in surface sulfides (e.g., FeS, Sprague et al., 2009; Vern- dients in the nebula, differential re-condensation related to chon- azza et al., 2010) could also be consistent with the low-FeO content drule-forming and impact events, aerodynamic separation of in silicates (Blewett et al., 2009). In the latter case, Ti could be pres- particles, and preferential evaporation and/or impact removal of ent in less oxidized minerals than ilmenite. silicates. The latter two processes could have involved km-size The presence of Fe-bearing oxides (magnetite, ulvospinel, planetesimals, Moon- to Mars-sized planetary embryos, and Mer- ilmenite) in Mercury’s rocks seems unlikely because their forma- cury itself. tion requires oxidized conditions that are inconsistent with FeO- Thermodynamic models for condensation of the solar nebular poor silicates (Riner et al., 2010, this paper). In reduced (FeO-poor) gas (e.g., Grossman, 1972; Grossman and Larimer, 1974; Ebel, magmatic conditions, Ti could be present in Mg-rich geikielite– 2006) demonstrate that high temperature condensates ( 1300– ilmenite solid solution (MgTiO3–FeTiO3)(Riner et al., 2010), rutile 1450 K) have an elevated metal/silicate ratio compared to low-T (TiO2), perovskite (CaTiO3), or even sulfides. Although ilmenite condensates. Preferential accumulation these condensates in has been suggested as an opaque spectrally neutral phase to ac- high-T inner nebula could have accounted for a high Fe/Si ratio count for low reflectance in visible and near-infrared ranges (Rava of planetesimals that formed Mercury. High gas (H2) pressure in and Hapke, 1987; Robinson et al., 2008; Denevi et al., 2009; Vern- the innermost nebula would also facilitate early high-T condensa- azza et al., 2010), it does not provide a good match with near-infra- tion of Fe-metal compared to Mg-silicates (Grossman, 1972; Lewis, red emission spectra (Sprague et al., 2009). In contrast, rutile could 1972). Subsequent chondrule-forming events caused by gas shock be present in surface materials based on interpretation of mid- waves (Desch and Connoly, 2002) could have favored metal con- infrared emission spectra (Sprague et al., 2009) and near-infrared densation at high gas densities over formation of silicates (includ- reflectance spectra (Vernazza et al., 2010). Titanium-rich sulfides ing Fe-silicates) through re-condensation. This process could have are present in enstatite meteorites (Keil, 1969) that are depleted been more efficient in the innermost nebula because of high pres- in FeO-bearing minerals (Keil, 1968). The presence of spectral fea- sure. High-T removal of silicates through evaporation (cf. Fegley tures of enstatite in Mercury’s mid-infrared spectra (e.g., Sprague and Cameron, 1987) from separate chondrule melts and/or nebular et al., 2009) does not exclude the occurrence of Ti sulfides in dust could have also contributed to a high Fe/silicate ratio on Mer- FeO-depleted surface materials. This paper (Section 5) shows that cury. Extremely reduced conditions, which may reflect elevated C/ rutile, Ti2O3, and some Ti sulfides could form in FeO-poor mag- O and H2/H2O ratios in the inner nebula (e.g., Ciesla and Cuzzi, matic conditions. It is possible that a majority of Fe and Ti is in re- 2006), also support preferential formation of Feo-rich condensates duced forms that reflect redox conditions during magmatic (Larimer and Batholomay, 1979), and account for the low-FeO con- crystallization. tent in accreted materials. Note that condensation of nebular gas 26 M.Yu. Zolotov / Icarus 212 (2011) 24–41 with solar abundances also leads to reduced condensates with 1989; Keil et al., 1994) demonstrate that extremely reduced plane- minute FeO content in silicates (e.g., Ebel, 2006; Fedkin and Gross- tesimals were also involved in these processes. Independent of the man, 2006; Grossman et al., 2008). origin of metal-rich planetesimals, an aerodynamic separation of High T–P conditions in the inner solar nebula may imply a low denser and/or larger planetesimals could have led to an elevated amount of condensed volatiles (at least H, S, and alkalis). However, Feo/silicate ratio in the high-P inner nebula and may have contrib- elevated C/O and H2/H2O ratios in the inner nebula could have led uted to the formation of metal-rich Mercury (Weidenschilling, to accumulation of S, C, and N in reduced solids (e.g., Larimer and 1978). Batholomay, 1979). If Mercury’s primary materials were formed at Our working hypothesis is that Mercury formed from reduced, C/O > 0.8, the volatiles could have been condensed in phases high-T, high-P, anhydrous nebular condensates that were affected which occur in enstatite chondrites (Mg, Ca, Cr and Ti sulfides, C by removal of silicates in major collisions occurring before and/or species, and nitrides; Keil, 1968, 1989; Brearley and Jones, 1998) after formation of the planet. Although enstatite chondrites and (e.g., Larimer, 1975; Larimer and Batholomay, 1979; Grossman CB-like chondrites provide insights about reduced nebular solids et al., 2008). This notion and the low FeO in enstatite meteorites and metal–silicate fractionation through impacts, it is possible that led to the suggestion of significant contribution from enstatite none of these groups of meteorites represent the composition of chondrite (EC) material to the bulk composition of Mercury (Was- Mercury. In particular, FeO-poor and Feo-rich condensates of the son, 1988; Sprague et al., 1995). Enstatite chondrites are isotopi- solar gas composition could be among Mercury’s protoliths, though cally identical to Earth and could be similar to the dominant they do not contain minerals found in enstatite chondrites. Like- constituent of our planet (Javoy et al., 2010). A deficiency of H2O wise, the presence of hydrated carbonaceous clasts in CB chon- gas in the inner nebula caused by diffusion of water vapor toward drites (Bonal et al., 2010) and their oxygen isotopic composition the water condensation front (Stevenson and Lunine, 1988; Cyr (Clayton and Mayeda, 1999; Rubin et al., 2003; Ivanova et al., et al., 1998; Ciesla and Cuzzi, 2006) and/or evaporation of con- 2008) indicate their formation in the main asteroid belt. densed C-bearing materials that migrated toward the Sun (Larimer and Batholomay, 1979; Lodders and Fegley, 1993) could account 2.2. Hydrogen and water for even greater contribution of the EC-like material to the compo- ÀÁ þ þ sition of Mercury. Condensation models show that a 85% depletion Hydrogen-bearing ionized molecules H2O ; H2O; OH and + of nebular H2O gas content would lead to formation of minerals H2S have been tentatively identified in Mercury’s exosphere by found in enstatite chondrites (Hutson and Ruzicka, 2000; Pasek MESSENGER (Zurbuchen et al., 2008). These molecules could have et al., 2005). Mn–Cr isotopic systematics (Shukolyukov and Lugm- resulted from sputtering of surface materials, which may include air, 2004) and xenon isotopes (Jee-Yon et al., 2009) in enstatite water ice (Zurbuchen et al., 2008) and/or protons implanted with meteorites also indicate their formation in the hot inner solar neb- the solar wind. Water ice (Slade et al., 1992) may be associated ula sunward of the asteroid belt. However, the Fe/Si ratio in EC is with radar-bright craters in the north polar region (Harmon and similar to solar abundances and a metal–silicate separation is re- Slade, 1992; Harmon et al., 2001) and could represent molecules quired to form Mercury from reduced condensates, which may or trapped after cometary/meteoritic impacts (Moses et al., 1999). may not contain minerals found in enstatite chondrites. Although endogenic sources of Mercury’s H are not excluded,
Disruptive collisions of planetesimals and differentiated plane- the interior of Mercury is likely to be depleted in H and H2O. The tary embryos (e.g., Wetherill, 1986; Asphaug et al., 2006) would low-FeO content in surface silicates and suggested FeO-poor man- provide an efficient mechanism for metal–silicate separation. tle (Robinson and Taylor, 2001) are indicators of water scarcity that Numerical simulations of Benz et al. (1988, 2007) demonstrate that also implies reduced conditions. Water is the major oxidizing agent the density of Mercury could be accounted for by stripping of its in the solar nebula and inside asteroids and planetary bodies (Le- silicate mantle by major impact(s). Dynamical models show that wis and Prinn, 1984; Dreibus and Wänke, 1989; McSween and Lab- major collisions were more frequent within 2 AU from the Sun, otka, 1993; Arculus and Delano, 1987; Kasting et al., 1993; Ciesla and that iron meteorites could have formed through impact dis- and Cuzzi, 2006), and reduced rocks often reflect processes that oc- ruptions of differentiated bodies in the terrestrial planet region curred under deficiency of water. The low-FeO content in surface (Bottke et al., 2006). High relative velocities of planetesimals and silicates suggests a small (if any) contribution from processes such planetary embryos at low AU should have led to strong collisions as aqueous alteration of parent planetesimals, accretion of icy and o (e.g., Wetherill, 1986). Rapid accretion of planetesimals at low hydrated planetesimals, and Fe –H2O interaction in the interior. AU also favored their fast differentiation (through decay of 26Al, Outward diffusion of nebular water vapor and a limited pene-
Grimm and McSween, 1993) and made metal–silicate separation tration of H2O-bearing planetesimals into the inner solar nebula more efficient. In addition to silicates, some volatiles could have could have created H2O-deficient conditions (e.g., Stevenson and been lost through decompression and boiling of silicate fluids in Lunine, 1988; Ciesla and Cuzzi, 2006). Together with high temper- collisions (Melosh et al., 2004; Asphaug et al., 2006) followed by ature, these conditions would prevent gas–solid type hydration rapid quenching before re-condensation of volatiles. (Lewis, 1972, 1988). A metal–silicate separation through major col- The metal-rich (60–70 vol.%) CB carbonaceous chondrites (e.g., lisions (Benz et al., 1988, 2007) implies a loss of H and water Weisberg et al., 2001; Krot et al., 2003) could have formed through (Asphaug et al., 2006). Even if some water had accreted on Mer- collisions of planetary embryos (Campbell et al., 2002; Rubin et al., cury, it would have been consumed in oxidation reactions with 2003; Krot et al., 2005). CB chondrites are strongly depleted in Feo-metal (e.g., Dreibus and Wänke, 1989), elemental C (Holloway moderately volatile lithophile elements (Na, K, Mn, Zn, Se, etc.) and Jakobsson, 1986; Kadik et al., 2006) and other reduced solids and have high abundances of lithophile refractory elements (Weis- (Fe3P, MgS, etc.) leading to the formation and escape of H2 and berg et al., 2001; Krot et al., 2003). They have low-FeO contents in CO. If Mercury’s mantle formed through upward crystallization of silicates (<4 mol.%), and the high Feo/silicate ratio makes them a magma ocean (Brown and Elkins-Tanton, 2009), water might plausible analogs for Mercury’s bulk composition (Taylor and Scott, have been expelled toward the surface and partially degassed at 2003). Several other metal-rich meteorites (CH chondrites, palla- low-P conditions. Low solubility of water in Feo-metal saturated sites, ureilites, mesosiderites, and IVA iron meteorites) with signs magmas (Kadik et al., 2004 and refs. therein) would also favored of rapid cooling could have formed through major collisions as well water loss. Finally, possible impact removal of the upper part of (Keil et al., 1994; Yang et al., 2007; Scott, 2007). Signs of a major, the mantle (Benz et al., 1988, 2007) would have led to major loss early impact in the Shallowater enstatite achondrite (aubrite) (Keil, of water and H. M.Yu. Zolotov / Icarus 212 (2011) 24–41 27
The discovery of H-bearing glasses (4–46 ppm H O; Saal et al., 2 -8 o 2008) and phosphates in lunar basalts (Boyce et al., 2010) implies IW: 2Fe0.947O = 1.894Fe + O2 IRI: 2FeTiO = 2TiO + 2Feo + O IW that some water either survived a Moon-forming collision or was 3 2 2 -10 o trapped after the event. The estimated pre-eruptive water content ITOI: 2FeTiO3 = Ti2O3 + 2Fe + 1.5O2 RTO: 4TiO = 2Ti O + O in lunar glasses ( 700 ppm; Saal et al., 2008) could be considered 2 2 3 2 -12 M-A as an upper limit for water content in Mercury’s rocks because IRI , bar eO melt is enriched in water compared to mantle source rocks, and be- 2 F wt% CB/CH O 6 cause of the possible water-bearing nature of proto-Earth. f -14 eO % F log t 2 w EC ITOI 2.3. Mantle oxidation state -3 -16 IW RTO Oxygen fugacity (fO2) of lunar mare basalts (18–21 wt% FeO; -6 IW Basaltic Volcanism, 1981)is0to 2 log units below the iron–wüs- -18 tite buffer (IW, Eq. (1) in Table 1)(Wadhwa, 2008). The lesser FeO 1200 1300 1400 1500 1600 1700 1800 content on Mercury compared to lunar values implies more re- Temperature, K duced conditions. For 2–6 wt% FeO in silicate rocks, Malavergne Fig. 1. Oxygen fugacity at reduced conditions that corresponds to some mineral et al. (2010) estimated Mercury’s mantle fO2 as 3–6 log units below buffers and suggested analogs for Mercury’s mantle melts. The dotted curves refer IW (IW 3toIW 6) (Fig. 1). They assume ideal melt in which FeO to conditions of 3 and 6 log fO units below the iron–wüstite buffer (IW, Eq. (1) in o 2 activity is equal to concentration in equilibrium with Fe -metal. Table 1). IRI, ilmenite–rutile–iron buffer (Eq. (2)); ITOI, ilmenite–titanium oxide– Consideration of the non-ideality of silicate melts saturated iron buffer (Eq. (3)); RTO, rutile–titanium oxide buffer (Eq. (4)). The gray curves with Feo-metal using empirical equations from Ariskin et al. represent melts with 2 wt% and 6 wt% FeO based on the Ariskin et al. (1993) melt (1993) leads to more definitive fO evaluations. These equations model for 1500–1700 K. The curves marked CB/CH and M-A correspond to partial 2 melts of CB/CH chondrites and Mercury’s mantle composition from Morgan and describe relationships between fO2, FeO molar content, and con- Anders (1980), respectively. The circle symbol refers to a partial melt of an enstatite centrations of other major oxides in melts. For plausible magma chondrite. formed through partial melting of Mercury’s mantle compositions
Table 1 Chemical equilibria and their equilibrium constants.
Reaction Log K at 1 bar 298 K 1500 K 1700 K
o (1) 2Fe0.947O (wüstite) = 1.894Fe +O2(g); IW buffer 85.8 11.6 9.53 o (2) 2FeTiO3 (ilmenite) = 2TiO2 (rutile) + 2Fe +O2(g); IRI buffer 93.6 12.9 10.7 o (3) 4/3FeTiO3 = 2/3Ti2O3 + 4/3Fe +O2(g); ITOI buffer 103 14.1 11.6
(4) 4TiO2 = 2Ti2O3 +O2(g); RTO buffer 121 16.4 13.2
(5) 2MgS + O2(g) = 2MgO + S2(g) 65.8 11.9 10.3
(6) 2CaS + O2(g) = 2CaO + S2(g) 32.6 5.90 5.1
(7) 2MnS + O2(g) = 2MnO + S2(g) 36.6 6.42 –
(8) TiS2 +O2(g) = TiO2 +S2(g) 72.5 14.1 –
(9) TiS2 + 0.75O2(g) = 0.5Ti2O3 +S2(g) 42.2 10.0 –
(10) 2TiS + 1.5O2(g) = Ti2O3 +S2(g) 139 23.6 –
(11) 2MgS + 2SiO2 +O2(g) = 2MgSiO3 +S2(g) 77.3 14.1 12.4
(12) 2MgS + SiO2 +O2(g) = Mg2SiO4 +S2(g) 76.2 13.8 12.0
(13) 2CaS + 2SiO2 +O2(g) = 2CaSiO3 +S2(g) 62.9 12.1 10.6
(14) 2CaS + SiO2 +O2(g) = Ca2SiO4 +S2(g) 55.4 10.8 9.57
(15) 2MnS + 2SiO2 +O2(g) = 2MnSiO3 +S2(g) 45.5 8.03 6.93
(16) 2MnS + SiO2 +O2(g) = Mn2SiO4 +S2(g) 45.1 7.18 6.66
(17) 2TiS2 = 2TiS + S2(g) 54.3 3.69 – o a (18) 2FeS = 2Fe +S2(g) 50.0 4.8 – o (19) MgS + FeO = MgO + Fe + 0.5S2(g) 11.1 0.08 0.26 o (20) MgS + FeTiO3 = MgO + TiO2 +Fe + 0.5S2(g) 13.9 0.52 0.23 o (21) TiS2 + 2FeTiO3 = 3TiO2 + 2Fe +S2(g) 21.1 1.14 – o (22) 1.5MgS + FeTiO3 = 0.5Ti2O3 + 1.5MgO + Fe + 0.75S2(g) 27.7 1.65 0.96
(23) MgS + 2TiO2 =Ti2O3 + MgO + 0.5S2(g) 27.7 2.73 1.46
(24) Mg2SiO4 + SiO2 = 2MgSiO3 1.16 0.31 0.36 o (25) SiO2 =Si +O2(g) 150 22.4 18.7 (26) Co + FeO = Feo + CO(g) 20.0 2.45 3.09 o (27) C + 0.5O2(g) = CO(g) 24.0 8.47 7.99 o o (28) C + FeTiO3 =Fe + TiO2 + CO(g) 22.8 2.0 2.62 o o (29) C + 0.667FeTiO3 = 0.667Fe + 0.333Ti2O3 +CO 27.4 1.43 2.21 o (30) C + 2TiO2 =Ti2O3 +CO 39.6 0.27 1.39 o o (31) C + SiO2 =Si + 0.5O2(g) + CO(g) 126 14.0 10.7
(32) TiN + O2(g) = TiO2 + 0.5N2(g) 102 16.6 14.2
(33) 2TiN + 1.5O2(g) = Ti2O3 +N2(g) 143 25.1 21.7
(34) 2NH3(g) = N2(g) + 3H2(g) 5.75 8.41 8.87
(35) H2O(g) = 0.5O2(g) + H2(g) 40.0 5.73 4.67
Notes: The equilibrium constants are calculated using thermodynamic properties of corresponding solids, liquids, and gases mainly from Gurvich et al. (1989–1994), Pankratz
(1987), Robie and Hemingway (1995), and Barin (2004). Although O2 gas is not present in magmas, it is used to formally quantify fO2 at chemical equilibria. In reality, some of o o the species (FeO, MgS, SiO2, etc.) could be components of magma, a phase or a solid solution. In all equations, C ,Fe , TiO2, FeTiO3, and MgSiO3 represent graphite, iron metal, rutile, ilmenite, and enstatite, respectively. a Extrapolated from T < 1400 K. 28 M.Yu. Zolotov / Icarus 212 (2011) 24–41
(Taylor and Scott, 2003), concentrations of 2–6 wt% FeO result in the redox conditions of IW 0.5 imply a significant oxidation of fO2 of IW 2.3 to IW 3.4 at 1600–1700 K (Figs. 1 and 2). Specifi- the innermost solar nebula through incorporation of water-bearing cally, the redox state of a partial melt of the Indarch (EH4) enstatite planetesimals migrated from larger heliocentric distances (Ciesla chondrite at 1700 K with 0.25 wt% FeO (McCoy et al., 1999; Taylor and Cuzzi, 2006; Grossman et al., 2008). Although this migration and Scott, 2003) corresponds to IW 5.3 (log fO2 = 14.9). Melt de- could explain elevated oxidation states of some materials (e.g., rived from a mantle composition consisting of the silicate compo- Type II chondrules with high FeO content) formed in the asteroid nents of metal-rich CB/CH-like chondrites (Krot et al., 2001) with belt, these planetesimals are unlikely to have survived high-T con- 3.2 wt% FeO (Taylor and Scott, 2003) corresponds to IW 3 ditions in the innermost nebula. Considering limited contribution
(log fO2 = 13.5) at 1600 K, and IW 3.1 (log fO2 = 12.6) at from water-bearing bodies to the composition of Mercury (O’Brien 1700 K. Similar results are obtained for melt with 3.7 wt.% FeO et al., 2006), a reduced mantle seems likely. Metal–silicate frac- (Taylor and Scott, 2003) derived from the mercurian mantle com- tionation through collisions, gas drag, and igneous differentiation position of Morgan and Anders (1980):IW 2.9 (log fO2 = 13.4) at (see Section 1) would not have changed the reduced nature of 1600 K, and IW 3 (log fO2 = 12.5) at 1700 K. materials that formed Mercury. These estimates show that fO2 in Mercury’s magmas could be close to the conditions determined by mineral equilibria that in- 2.4. Sulfur, carbon, chlorine, nitrogen, and alkalis clude Ti(IV) and Ti(III) oxides (ITOI and RTO buffers, Figs. 1 and 2). The conditions of the IRI buffer (Eq. (2) in Table 1) correspond 2.4.1. Sulfur to high FeO melt contents that are not consistent with remote Estimations of Mercury’s bulk S content strongly depend on the sensing observations of surface materials (Fig. 2) (except, possibly, formation scenario. Accretion of Mercury from EC-like materials Lawrence et al., 2010, interpretation). In other words, the presence would lead to solar abundances of S (Table 2). In addition to FeS, of ilmenite contradicts low FeO concentration in silicates, and pure Mg-, Ca- and Cr-sulfides are the likely products of nebular pro- ilmenite may not be present in FeO-poor Mercury’s rocks (Riner cesses in H2O-depleted inner nebula (Cyr et al., 1999; Pasek et al., 2010). et al., 2005; Grossman et al., 2008). Conversely, formation of Mer- The possibility of crystallization of ilmenite before pyroxene at cury from metal-rich products of planetary embryo collisions could IW 0.5 (Usselman and Lofgren, 1976) could be reconciled with lead to an S-depleted planet, considering the low S content in CB the low-FeO content in surface silicates, and the lunar-type oxi- chondrites (Table 2). However, independent of the formation sce- dized mantle cannot be excluded (Lawrence et al., 2010). However, nario and redox state of accreted solids, Mercury’s silicate mantle abFeO wt % 12 5 6 7 -1 IRI, 1700 K IRI, 1500 K K 2 0 - 600 K -12 170 1 IW -2 ITOI, 1700 K K 1500 ITOI, 1500 K
K 3 0 - -3 150 W K
, bar -14 I 1700 2 4 - RTO, 1700 K O W
I f
relative to IW -4 2 O log -16 f RTO, 1500 K -5 log
-18 -6 0 5 10 15 20 0123456 FeO mole % FeO mole %
Fig. 2. Oxygen fugacity as a function FeO concentration in silicate melts. The data are obtained with empirical equations for silicate systems saturated withFeo-metal (Ariskin et al., 1993). The curves refer to magmas formed through partial melting of plausible Mercury’s mantle compositions from Taylor and Scott (2003). Solid, dashed, and dotted curves represent partial melts of an enstatite chondrite (McCoy et al., 1999), averaged silicate components of metal-rich CB/CH like chondrites (Krot et al., 2001), and a
Mercury’s mantle (Morgan and Anders, 1980), respectively. In (a), thin dashed curves show fO2 values relative the IW buffer at specific temperatures and melt compositions.
In (b), the data are presented relative to the IW buffer for low-FeO contents. The dotted lines refer to fO2 buffers depicted in Fig. 1. The upper horizontal axis corresponds to wt% FeO in the melt derived from metal-rich chondrites.
Table 2 Abundances of elements in some chondrites and aubrites.
S, wt% C, wt% Cl, ppm N, ppm CI carbonaceous chondrites 5.41 3.45 700 3180 H ordinary chondrites 2.0 0.21 140 48 L ordinary chondrites 2.2 0.25 270 43 LL ordinary chondrites 2.1 0.31 200 70 EH enstatite chondrites 5.6 0.39 570 420 EL enstatite chondrites 3.1 0.43 230 240 Aubrites (enstatite achondrites) 0.24–0.86a 0.02–0.14b 1–31c 10–100b CB chondrites 0.52e 0.37e – 25–106d
Notes: Chondritic data are from Lodders and Fegley (1998). a Gibson et al. (1985). b Grady et al. (1986). c Garrison et al. (2000), without Shallowater aubrite which contains 400 ppm Cl. d Krot et al. (2003) and Ivanova et al. (2008). e Rubin et al. (2003) and Ivanova et al. (2008). M.Yu. Zolotov / Icarus 212 (2011) 24–41 29 may only contain accessory sulfides because they would likely a 10 accumulate in the core (cf., Riner et al., 2008). Besides FeS, this accumulation could have affected Mg-, and Ca-sulfides, if they EC gas 1 were present (McCoy et al., 1999). The effect of subsequent major collisions (Benz et al., 1988, 2007) on concentration of S in the mantle is not clear. 0.1 Mercury's mantle? Sun 0.01 EH L LL CI 2.4.2. Carbon Cl/S mole ratio EL H C/S Mole ratio Kilauea The solar abundance of C is an order of magnitude higher than 0.001 the abundance in CI carbonaceous chondrites and is about two or- 0.1 1 10 ders of magnitude higher than in ordinary, enstatite, and CB chon- b 10 drites (Table 2; Asplund et al., 2005; Lodders and Fegley, 1998). EC gas This corresponds to a decrease in C content in the main asteroid 1 Sun belt where carbonaceous bodies become less abundant toward the Sun (Gradie et al., 1989). The limited formation of C-bearing CI 0.1 Mercury's mantle? solids in the inner solar nebula could be related to formation of EH EL CB CO from condensed organic species transported from higher helio- 0.01 N/S mole ratio LL centric distances (AU > 5) (Cuzzi et al., 2003). Reduced nebula H L conditions have not cased accumulation of abundant graphite in 0.001 enstatite chondrites, though these meteorites are slightly enriched 0.1 1 10 20 in C compared to ordinary chondrites (Table 2). Nebular condensa- C/S mole ratio tion models of Lewis (1988) do not predict C solids formed at high T–P conditions at the heliocentric distance of Mercury. It follows Fig. 3. Normalized elemental abundances in chondritic materials, gases released that solids condensed in the vicinity of the Sun could be C depleted. from heated enstatite chondrites (EC gas, Muenow et al., 1992), and Kilauea 1918 Nebular condensates (and Mercury) could be C-rich only if con- summit eruption (Symonds et al., 1994). Data for carbonaceous (CI, CB), ordinary (H, L, LL) and enstatite chondrites (EH, EL) are from Table 2, and solar photosphere data densed at higher C/O ratios compared to enstatite chondrites, that are from (Asplund et al., 2005). The gray ovals show possible elemental ratios in is at C/O > 0.83 when high-T condensates become rich in graph- Mercury’s mantle and magmas. ite, TiC, and SiC (e.g., Larimer and Batholomay, 1979; Grossman et al., 2008). Although such enrichment could be caused by mas- sive infall of C-rich solids into the inner solar nebula (Lodders ula could have led to formation of metal halides, as exemplified by and Fegley, 1993) no meteorite rich in these phases is found. Dy- reactions namic models show that only a limited fraction of carbonaceous 2HClðgÞþNa2Oðin gas or solidsÞ!2NaCl þ H2OðgÞ planetesimals formed at AU > 2 could have accreted on Mercury (O’Brien et al., 2006). Major collisions of planetary embryos may HClðgÞþNaAlSi O ðsolidÞ!NaCl þ 0:5Al O þ 3SiO have no major effect on C abundance, as implied from similar C 3 8 2 3 2 abundance in CB and ordinary chondrites (Table 2). Although C/S þ 0:5H2OðgÞ ratio in Mercury could vary from that in enstatite and carbona- These interactions of HCl(g) could have led to initial Cl abundances ceous chondrites (Fig. 3), lower ratios similar to EC seem probable. on Mercury as high as in CI chondrites and the Sun. During differentiation, accreted organic carbon could have been Highly metamorphosed EL6 chondrites are depleted in Cl com- lost as CO and hydrocarbon gases formed through thermal process- pared to EH4 chondrites (Rubin and Choi, 2009). One can suggest ing, oxidation, and graphitization of organics. However, in anhy- some loss of Cl during thermal/impact metamorphism after accre- drous high-P igneous systems, the majority of C could have not tion of Mercury and/or planetary embryos. Chlorine and alkali ele- degassed because of precipitation of low-solubility graphite (see ments could have been accumulated in the upper horizons of Section 4). Independent of the origin of Mercury and bulk C con- Mercury through magmatic differentiation possibly involving up- tent, a major portion of accreted C could be sequestered in the me- ward crystallization of a magma ocean (Brown and Elkins-Tanton, tal core. Although reducing conditions in the mantle may imply a 2009). Subsequently, these elements could have been removed stable existence of native C species, Mercury’s silicate rocks and through impact stripping. Despite the uncertainties, a Cl/S mole ra- magmas may not be enriched in C compared to other terrestrial tio of 0.01 characterizes many solar system objects (Fig. 3a) and planets. could be used a proxy for Mercury’s mantle rocks and magmas.
2.4.3. Chlorine 2.4.4. Nitrogen The abundance of Cl in Mercury’s mantle rocks and magmas is Chondrites are strongly depleted in N compared to cometary hard to constrain because of the unknown bulk planetary content materials (Jessberger, 1999; Bockelée-Morvan et al., 2004) and so- and limited solubility data for dry reduced melts. On the one hand, lar abundances (Asplund et al., 2005). This indicates a major devol- a low Cl content in ordinary chondrites compared to CI chondrites atilization of N-bearing compounds (organics, CHON particles, and
(Table 2) and solar abundances implies Cl loss in high-T processes NH3 hydrate) in the inner solar nebula. The position of Mercury in in the inner nebula and thermally processed planetesimals. The the Solar System implies a limited delivery of N in organic-bearing secondary nature of all Cl phases in ordinary and carbonaceous planetesimals and planetary embryos (O’Brien et al., 2006). How- chondrites also suggests a low-T (<200 K) condensation of Cl (as ever, Mercury’s rocks may not be depleted in N compared to other HCl hydrates, Zolotov and Mironenko, 2007) which may account terrestrial planets because of possible accretion of refractory for Cl depletion in the inner Solar System. On the other hand, re- nitrides (TiN, Si3N4) and N-bearing metal alloys condensed from duced conditions in the inner nebula could have favored high-T reduced nebular and/or impact-generated gases. In enstatite chon- formation of halides consistent with elevated Cl abundances in drites, the presence of reduced N solids (Keil, 1968, 1989; Rubin EH enstatite chondrites (Table 2; Rubin and Choi, 2009). A defi- and Choi, 2009) accounts for significantly larger N content com- ciency of water vapor and elevated P–T conditions in the inner neb- pared to ordinary and some carbonaceous (CB, CV, CO) chondrites 30 M.Yu. Zolotov / Icarus 212 (2011) 24–41
(Table 2; Lodders and Fegley, 1998). The association of N with the have high FeO/MgO molar ratios ( 0.3–2, Mayne et al., 2009; McS- metal phase and its high-T ( 1000 °C) release from a CB/CH mete- ween et al., 2010) and may not be considered as close analogs for orite (Ivanova et al., 2008) may indicate trapping of degassed N Mercury’s igneous melts and cumulates. through condensation of metal from impact gases (Sugiura et al., Low-FeO content in Mercury’s surface silicates also implies a 1999). It follows that major collisions involving differentiated FeO-poor mantle (Robinson and Taylor, 2001) that would produce planetary embryos and Mercury may not have resulted in a major high-T melts. A mantle rich in refractory oxides is consistent with N loss. Abundance and speciation of N in enstatite and CB/CH chon- the non-global appearance of volcanic features and spectral signs drites could be used as a proxy for Mercury’s initial materials (Ta- of Mg- and Ca-rich silicates in surface materials (see Section 1). ble 2 and Fig. 3b). The refractory mineralogy may represent accreted Fe-depleted sil- The high stability of nitrides during metamorphism at extre- icates (enstatite, forsterite, diopside, and feldspars) and remained mely reduced conditions (Petaev and Khodakovsky, 1986; Fogel in the mantle after segregation of Feo–FeS melts into the core. et al., 1989) could have prevented major loss of N before the igne- Alternatively, the mantle could have formed through upward crys- ous differentiation of Mercury and planetary embryos. The high tallization of a magma ocean leading to a FeO-rich upper layer stability of nitride complexes in reduced, H-free and high-P mag- though fractional crystallization (Brown and Elkins-Tanton, mas (see Section 4) could also have limited N degassing from a 2009). Later on, this layer could have been stripped by impact(s) putative magma ocean, as suggested for the terrestrial counterpart (Benz et al., 1988) or even sunk into the lower mantle (Brown (Libourel et al., 2003). The presence of N in meteoritic Feo alloy also and Elkins-Tanton, 2009). The putative evaporation of an early implies some sequestering of N into the metal cores on planetary mantle also results in refractory rocks (Cameron, 1985; Fegley embryos and Mercury (cf., Fegley, 1983). and Cameron, 1987). The interpretation of MESSENGER neutron spectrometer data 2.4.5. Alkalis suggest Mercury’s rocks similar to lunar mare basalts from Luna Alkalis could have not been abundant in high-T nebular conden- 24 and Luna 16 samples (Lawrence et al., 2010). These data are sates and alkali oxides may compose only a few tenths of weight inconsistent with the melt models for FeO-poor protoliths (Taylor percent of materials accreted on Mercury (Lewis, 1988). Alkalis and Scott, 2003; Burbine et al., 2002). Although neutron spectrom- could have been depleted through high-T evaporation in the inner eter data would revolutionize our views on Mercury’s magmatism, nebula (through chondrule formation and/or collisions), after they need to be confirmed with independent orbital or in situ accretion (Fegley and Cameron, 1987), or stripped by impacts to- measurements. gether with a feldspathic crust (Benz et al., 1988; Lewis, 1988). Although the presence of Na and K in the exosphere (Potter and 4. Solubility of S, C, Cl, and N in reduced anhydrous magmas Morgan, 1985, 1986; Domingue et al., 2007; McClintock et al., 2009) is consistent with feldspar-bearing surface rocks, it does The specifics of degassing during eruptions and the abundances not imply an alkali-rich mantle. Enhanced K emission over the Cal- of volatiles in magma are affected by their solubility, which is non- oris basin (Sprague et al., 1990) may indicate a low degree of melt- linearly dependent on T, P, and fO2. Reduced melts could be en- ing rather than alkali enrichment in the mantle. Note that the riched in S, Cl, and N because of the formation of non-volatile melt apparent lack of crustal subduction implies limited recycling of complexes without oxygen. However, some volatiles could be alkalis and other incompatible elements in igneous processes. A undersaturated because of limited abundances in mantle source deficiency of alkalis in the mantle also implies elevated liquidus regions. temperatures. Petrological and material science experiments performed in sil-
icate melts at fO2 below IW 2(Fincham and Richardson, 1954; Fo- 3. Temperature and composition of mantle magmas gel et al., 1996; McCoy et al., 1999; Siebert et al., 2004; Berthet et al., 2009) reveal S solubility over an order of magnitude higher Temperature and composition of magmas mainly depend on the than at more oxidizing conditions (Mavrogenes and O’Neill, uncertain mineralogy and volatile content in Mercury’s mantle. 1999; Holzheid and Grove, 2002)(Fig. 4a). The highest S solubility Petrologic experiments demonstrate that dry mantle rocks melt has been reported by Malavergne et al. (2007) (up to 10 wt% at at temperatures up to several hundred degrees higher than H2O- 10 kbar and IW 7) and Fogel et al. (1996) (7.8 wt% in CaS-satu- rich mantle substrates (e.g., Green, 1975; Righter et al., 2006). rated enstatite chondrite melts at 1 bar and 1673 K). In contrast Dry lavas on the Moon and Io are hotter than Earth’s counterparts to more oxidized systems, S solubility is higher in FeO-poor melts (Basaltic Volcanism, 1981; McEwen et al., 1998; Keszthelyi et al., and could be caused by formation of metal (Mg, Ca, Mn, Ti, etc.)
2007) and a deficiency of H2O on Mercury would contribute to ele- sulfide melt complexes (Fogel et al., 1996; McCoy et al., 1999; Fo- vated melt temperatures affecting the composition of silicate melt. gel, 2005). The anomalous S solubility in reduced melts is illus- Numerical modeling of partial melting using plausible mantle trated by the presence of 0.8–2.5 wt% S in basaltic glasses in compositions of Mercury demonstrates the formation of mafic to aubrites, which are the most S-rich samples among achondrites ultramafic melts with low-FeO contents (Taylor and Scott, 2003; and silicate planetary rocks (Fogel, 2005). These data show that S Section 2.3). Experimental melting of the Indarch (EH4) enstatite content in Mercury’s igneous materials would indicate redox con- chondrite (McCoy et al., 1999) suggests magma temperatures of ditions in magma and mantle source regions. 1670–1720 K at 20–50% partial melting. The high melt tempera- In anhydrous reduced magmas, the solubility of C is mainly con- ture reflects low abundances of FeO, alkalis, and water. The basaltic trolled by the equilibria of C-bearing solids (e.g., graphite) with CO normative composition of those melts ( 34 vol% plagioclase, 42% (Eqs. (26) and (27)). Low fO2 and temperature and high pressure fa- enstatite, 15% diopside, and 6% sulfides) has been suggested as vor the stability of graphite. The low solubilities of CO, CO2 (Pan a proxy for Mercury’s lavas (Burbine et al., 2002). Note that the use et al., 1991) and graphite, which precipitates, account for the low of enstatite achondrites (aubrites, Watters and Prinz, 1979)as solubility of C in magma (<