Next Article in Journal
One-Point Statistics Matter in Extended Cosmologies
Previous Article in Journal
Entropy of Artificial Intelligence
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Constraints on the Production of Phosphine by Venusian Volcanoes

1
Earth, Atmospheric and Planetary Sciences, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
2
School of Physics & Astronomy, Cardiff University, 4 The Parade, Cardiff CF24 3AA, UK
3
Department of Earth Sciences, University of Cambridge, Downing Street, Cambridge CB2 3EQ, UK
4
Institute of Astronomy, University of Cambridge, Madingley Rd., Cambridge CB3 0HA, UK
5
Department of Astronomy & Astrophysics, Northwestern University, Evanston, IL 60201, USA
6
Center for Interdisciplinary Exploration and Research in Astrophysics, Northwestern University, Evanston, IL 60201, USA
7
Blue Marble Space Institute of Science, Seattle, WA 98154, USA
8
Cavendish Laboratory, University of Cambridge, JJ Thomson Ave., Cambridge CB3 0HE, UK
9
Department of Physics, Massachusetts Institute of Technology, 77 Mass. Ave., Cambridge, MA 02139, USA
10
Department of Aeronautics and Astronautics, Massachusetts Institute of Technology, 77 Mass. Ave., Cambridge, MA 02139, USA
*
Author to whom correspondence should be addressed.
Universe 2022, 8(1), 54; https://doi.org/10.3390/universe8010054
Submission received: 30 November 2021 / Revised: 10 January 2022 / Accepted: 13 January 2022 / Published: 17 January 2022
(This article belongs to the Section Planetary Sciences)

Abstract

:
The initial reports of the presence of phosphine in the cloud decks of Venus have led to the suggestion that volcanism is the source of phosphine, through volcanic phosphides ejected into the clouds. Here, we examine the idea that mantle plume volcanism, bringing material from the deep mantle to the surface, could generate observed amounts of phosphine through the interaction of explosively erupted phosphide with sulfuric acid clouds. The direct eruption of deep mantle phosphide is unphysical, but a shallower material could contain traces of phosphide, and could be erupted to the surface. The explosive eruption that efficiently transports material to the clouds would require ocean:magma interactions or the subduction of a hydrated oceanic crust, neither of which occur on modern Venus. The transport of the erupted material to altitudes coinciding with the observations of phosphine is consequently very inefficient. Using the model proposed by Truong and Lunine as a base case, we estimate that an eruption volume of at least 21,600 km3/year would be required to explain the presence of 1 ppb phosphine in the clouds. This is greater than any historical terrestrial eruption rate, and would have several detectable consequences for remote and in situ observations to confirm. More realistic lithospheric mineralogy, volcano mechanics or atmospheric photochemistry require even more volcanism.

1. Introduction

Greaves et al. [1] reported detecting ~20 ± 10 ppb phosphine (PH3) in the atmosphere of Venus, using the James Clerk Maxwell Telescope (JCMT) and Atacama Large Millimeter/submillimeter Array (ALMA) radio telescopes. Subsequent to the publication of their paper, a recalibration of ALMA data and reprocessing the JCMT data [2] led to a revision of the abundance to 5–7 ppb with 5–6 α confidence. Infrared observation made by the SOIR (solar occultation in the infrared) instrument on Venus Express suggests a much lower upper limit abundance of <0.5 ppb above the clouds, at 60 km [3]. Other groups challenged the detection on technical grounds, suggesting that the detected line was artefactual [4,5] or was contamination by SO2 [6], both of which were rebutted in [7,8]. A reanalysis of Pioneer Venus mass spectrometry data supports the presence of PH3 in the clouds [9], with masses detected which are not consistent with any other species. Thus, two independent techniques that were applied in different decades suggest, but do not conclusively prove, the presence of PH3.
If PH3 is present in Venus’ atmosphere, it is worth exploring why it might be there. A wide range of potential sources have been suggested, both in the literature [10] and in informal discussions. Bains et al. [11,12] discuss these at length, emphasizing that potential sources must be shown to quantitatively account for the presence of at least 1 ppb PH3 at cloud level in Venus’ atmosphere to explain the detections. They show that, within the limits of the knowledge of the chemistry and photochemistry of the atmosphere of Venus, there is no known abiotic atmospheric source for PH3.
The chemistry and mineralogy of the surface and subsurface of the planet are less well known, but [11] show that thermodynamics make a surface volcanic source unlikely. They explored which gases could erupt on Venus with different temperatures, water abundances, and rock reduction regimes. Specifically, they show that, even if the surface was as reduced as the QIF redox buffer (implying the presence of metallic iron in erupted material), only ~1 part in 106 of phosphorus in the erupted material would be present as phosphine. They also state that plume volcanism is unlikely, given a semi-quantitative argument. They comment that “on Earth, mantle plume magma is estimated to rise on a timescale of 106–107 years at temperatures in excess of 3000 K, during which time phosphorus species would reach thermodynamic equilibrium relevant to the temperature and pressure of the upper mantle”. This temperature and this pressure strongly disfavor phosphine.
Ref. [11] shows that a biological source in the clouds is not ruled out by thermodynamics, but comment that Venus’ clouds are an extremely challenging environment due to high acidity and very low water activity, some of which have been partially addressed in follow-up work [13,14,15].
Subsequently, Truong and Lunine [16] revisited the concept that plume volcanism could be a source of PH3. Specifically, they propose a model in which plume volcanism could bring the P(-3) species (presumed to be metal phosphides) present in the deep mantle to the surface; explosive volcanism could lift them to the cloud layer, and the sulfuric acid in the clouds could then react with phosphides to make PH3.
In this paper, we examine the volcanic supply of phosphine scenario and determine under what conditions it could apply. We conclude that it would require a highly unexpected set of conditions for the model to be a valid source of PH3. The only possible scenario would require Venus to be entering a global resurfacing epoch, such as that proposed to have occurred ~300 Mya [17], with highly unexpected lithospheric chemistry, volcano dynamics, and atmospheric processes. Evidence of this could be found in a range of past and planned observations of Venus.

2. Materials and Methods

2.1. Phosphorus Partition between Iron and Silicate

The partition coefficient of phosphorus between iron and silicate in mixed iron/silicate melts (Dp) was calculated following [18] using Equation (1)
D p = a · ln [ f ( O 2 ) ] + b T + c · T P + d · n b o t + f
where a, b, c, d and f are constants, such that a = 0.404, b = 17455, c = 579, d = −0.725, f = −3.15, f(O2) is the oxygen fugacity, T is the absolute temperature, P is pressure in gigapascals, and nbo/t is the ratio of non-bridging oxygens to tetrahedral oxygens, which is related to the degree of polymerization of the silicate glass structure; in [18] nbo/t ranged from 0.88 to 2.57, and was taken to be 1.5 here. Pressure as a function of depth was taken from the Preliminary Reference Earth model [19], corrected for Venus gravity (0.88 at the surface). The temperature of the Venusian mantle is poorly constrained. We have taken the average temperature profile for Earth from [20] as a plausible estimate. We note that the solubility of phosphorus in melts is an inverse function of temperature [18], and so if mantle plumes are hotter than the surrounding rocks, which is likely [21], this model will overestimate the partition of phosphorus into iron.

2.2. Phosphorus Thermodynamics in Silicate Melts

The thermodynamics of phosphorus speciation between iron phosphide and mineral phosphate in silicate melts were calculated as described in [11]. In brief, Gibbs free energy values were obtained from [22,23,24,25,26], and the density of minerals was taken from The Minerals Project [27]. The free energy of a reaction of reactants to give products was calculated using Equation (2)
Δ G T , P r = n = 1 , p ( S n · Δ G n T ) n = 1 , r ( S n · Δ G n T ) + n = 1 , p ( S n · V n · P 10 ) n = 1 , r ( S n · V n · P 10 )
where Δ G T , P r is the free energy of the reaction at temperature T (in Kelvin) and pressure P (in bar), Sn is the stoichiometric number of moles of reagent n, Δ G n T is the free energy of formation of reagent n from its elements at temperature T, Vn is the molar volume of reagent n, p are products, and r are reactants. The physical state of the reactants is taken as a solid or a liquid depending on temperature, except for O2, which is a gas throughout. SiO2 was assumed to be quartz and metals were assumed to be in their lower oxidation state where relevant (i.e., Fe(II), not Fe(III) etc.). The free energy of reaction Δ G T , P r per mole of oxygen is converted to an oxygen fugacity (i.e., the activity of O2  ( { O 2 } ) at which the reaction is at equilibrium when the activity of oxidized and reduced iron species are the same) according to Equation (3)
f ( O 2 ) = L O G 10 ( { O 2 } ) ~ Δ G T , P r R · T · 2.303
The oxygen fugacity of QIF, IM and FMQ buffers, as a function of temperature and pressure was taken from [28]. The activity of free phosphorus in solution in iron at <10% phosphorus is negligible [24], and so, all the phosphorus can be considered to be present as iron phosphide, or as phosphate. Details of the reactions and of the calculation of the free energy of reaction and oxygen fugacity are provided in the Supplementary Materials.

2.3. Phosphorus and Oxygen Diffusion through Solids

The diffusivity of an element in a solid can be approximated by an Arrhenius Equation (4)
D o = A · e E a R · T
where Do is the diffusion constant (in this paper in units of cm2/day), R is the gas constant, T is the absolute temperature, and constants A and Ea are as listed in Table 1.
We note that the values for phosphorus diffusion through iron were only measured between 1050 °C and 1200 °C, and so, values outside this range represent an extrapolation. The rate of diffusion of phosphorus out of an iron sphere assumes a zero concentration outside the sphere, and the diffusion of oxygen into the sphere assumes a constant external concentration and no oxygen inside the sphere at the start. These were estimated by numerical integration; the result can be summarized by Equation (5)
t ½ = r 2 K · D o
where K is a constant = 0.376 for diffusion of phosphorus out, and 0.355 for diffusion of oxygen into the sphere.

3. Results

We describe several constraints on the production of phosphine by Venusian volcanoes, and determine quantitatively under what conditions Venusian volcanism could be a valid source of PH3. We build on the baseline case of [16] using the approach described in [11]. We first discuss why the mineralogy (chemical species) of erupted material is necessarily determined by the upper mantle and lithospheric minerology (Section 3.1), which leads to the detailed modelling of the abundance of phosphides in the mantle and lithosphere (Section 3.2). We then discuss why explosive volcanism is likely to be rare on Venus, (Section 3.3), why all volcanism is less likely to produce large ash clouds on Venus (Section 3.4), and some other limits on the production of phosphine from phosphide (Section 3.5), before reviewing the rate at which phosphide must be delivered to the clouds (Section 3.6). Section 4 integrates these calculations to show that the minimum rate of eruption that explains 1 ppb phosphine in the clouds is 21,600 km3/year, not 0.03 km3/year as suggested by [16], even given the most optimistic assumptions.

3.1. Eruptive Melt Minerology Is Determined in the Upper Mantle and Lithosphere

Before addressing the likely minerology of an erupting volcano, we briefly discuss why the minerology of that volcano is determined by the prevailing conditions of pressure, temperature and oxygen fugacity in the lithosphere and upper mantle, not the deep mantle.
We note that the minerology of a rock (the chemical species present, such as phosphate or phosphide) is distinct from the elemental abundance of a rock (the relative number of atoms of the elements in a rock). An assemblage of phosphorus, hydrogen and oxygen atoms could be present as phosphate and H2, or as phosphide and H2O. The elemental abundance depends solely on the source of the materials that made up the rock. The minerology—the chemical forms that those elements take—depends on the elemental abundance and on the pressure and temperature of the rock. In this paper, we are primarily concerned with the minerology of materials, not their elemental composition, which is assumed to be similar to Earth.
It is important to emphasize that oxygen fugacity is both directly a function of temperature, pressure, and rock composition, and indirectly affected by the effect of temperature and pressure on the mineralogy that is stable under a set of conditions (e.g., [31]). Oxygen fugacity is a property of the minerology of a rock, because it is the equilibrium position of a hypothetical reaction
MO ↔ M + O2
where MO is an oxidized mineral, M is a reduced mineral, and O2 is molecular oxygen. The oxygen budget of the rock will of course be the same at low pressure as at high pressure, but the f(O2) that the abundance of O translates to will be dominated by the mineralogy attained by the sample at a specific temperature and pressure. The oxygen fugacity of a system therefore changes with changes in temperature and pressure [28], i.e., with depth in a planet, independently of any change in the composition of the system. Thus, the oxygen fugacity of a rock of fixed bulk composition at the base of the lithosphere will be substantially different from the oxygen fugacity of the same rock in the deep mantle, even though they have the same intrinsic oxidizing capacity and the same elemental composition. Therefore, for the oxygen fugacity of the deep mantle to dictate the composition of volcanism, which necessarily occurs at low pressure, the relaxation of a system to low pressure thermochemical equilibrium (and its associated higher f(O2)) must be suppressed.
The deep mantle has been suggested, primarily on the basis of diamond inclusion chemistry, to be reducing [32]. Mantle convection, including mantle plumes, lifts this material to the base of the lithosphere, where it can pool and form melts that are subsequently erupted. A mantle plume is a buoyant upwelling of abnormally hot rock within the Earth’s mantle [33]. The ‘head’ of a plume forms a mushroom-shaped bolus which establishes the channel behind it by pushing aside the cooler, denser, more viscous material of the mantle [34,35]. Such plumes have been suggested as dredging xenoliths from the 250 km-deep cratonic keels of continents [36].
While deep mantle minerology can be reflected in inclusions in highly refractile xenoliths, the xenoliths are not representative of the bulk erupted material from plume volcanism, which is overwhelmingly basaltic, for two key reasons. First, although mantle material may have upwelled from a great depth, its melting typically begins in the upper mantle, and certainly its volumetrically dominant melting occurs there. Therefore, at the point at which melting has occurred, the material will have an upper mantle mineral assemblage, which produces more intrinsically oxidizing conditions (for the same bulk oxygen content of the rock), as noted above. Second, even melts that are produced at great depth must transit through the overlying mantle and lithosphere. On Earth, melt transport out of the mantle dominantly occurs from the porous flow along grain boundaries, which brings the melt into close contact with progressively lower pressure mineral assemblages. The geometry of melt flow out of the mantle, and the high temperatures at which this is occurring, both promote the equilibration of melts surrounding rock. For this reason, even the largest outpourings of magma on Earth, which have been supplied by deep upwellings of mantle rock, evidence depths of last equilibration in the upper mantle.
The rate at which the geometry of melt transport brings the rock and fluid into equilibrium can be illustrated as follows. Iron is expected to be solid at any depth in the mantle [37], and as deep mantle phosphides are expected to be almost exclusively dissolved in the iron phase, we will use iron as a model material. We will ask how big a sphere of iron must be before diffusion will allow any phosphorus in it to escape, or for external oxygen atoms in silicate rock to diffuse in and turn the phosphide to phosphate, if that reaction is thermodynamically favored. We address the thermodynamics in Section 3.2 below. We assume that the iron is moving at 10 m/year, which is typical of mantle plumes [38]; other mantle convection will be substantially slower, so this is a conservative lower limit on the size of the iron bodies required.
The results of this calculation are shown in Figure 1.
Rock grains on the mantle are typically less than 1 cm [41,42], and specifically, iron grains in silicate melts are found experimentally to be typical of mm size [18]. The simple model in Figure 1 confirms that such grains will be ‘frozen’ in the phosphorus chemical state present at 100–200 km depth. The minerology below that depth will be erased by the traverse of the mantle. Thus, while the elemental and isotopic composition of plume volcanic melts can reflect that of the deep mantle (e.g., [43]), their chemical state will necessarily reflect the minerology of the shallow mantle or lithosphere.
There are therefore four mechanisms that could bring deep mantle material to the surface. The first, and most common, is solid state convection, during which deep mantle minerology would equilibrate to shallow minerology, erasing the deep mantle redox signature. The second is fluid flow along the surface of solid grains, which, as summarized above, would equilibrate the fluid with the redox state of the rocks over which it flowed. The third would be the rapid bulk flow of fluid through large fractures. This is the mechanism that brings rare xenoliths from the shallow mantle through the lithosphere (e.g., [32,44,45,46,47,48,49]). Fractures that propagate from the lithosphere through the deep mantle are physically unrealistic. Fractures are initiated by an overpressure of fluid in the underlying rocks. The overpressure needed to generate and sustain a fracture which is a thousand kilometers long from the deep mantle to the surface is unrealistic, and such fracturing has never been observed on Earth. Lastly, the fourth mechanism would be fluid flow from the deep mantle, which could bring meter-size fragments of phosphide mineral to the surface. This is also unrealistic as the channels between grains are a millimeter or smaller in size.
The formation of metallic phosphides by the very rapid freezing of the mantle material has happened in the fragmentation of planetisimals that formed the iron meteorites [50]. However, this event resulted in the partial or complete disruption of the parent body, which clearly has not happened to Venus.
The minerology of bulk erupted volcanic material is therefore defined by the minerology of the upper 100–200 km of the planet, which may include mantle material depending on how thick the lithosphere is. This does not preclude that some erupted material is highly reduced (e.g., [51]). On Earth, specifics of the mantle and lithosphere minerology result in localized ore bodies of reduced rock (e.g., [52]), and these (and rarely other rocks) can contain highly reduced xenoliths as inclusions. However, with the exception of rare xenoliths, these are the product of localized lithospheric minerology, not lower mantle mineralogy.

3.2. Phosphorus Redox State in the Lithosphere and Upper Mantle

We estimate the abundance of phosphide in volcanic ejecta using thermodynamics, extending the analysis done by [11].
The chemical state of an element in rock can be estimated from the overall redox state of elements in the rock, the temperature, and the pressure. The redox state of elements of the rock is dominated by the state of its most abundant redox active element(s), which for mantle minerology is usually taken to be iron. Standard redox buffers are used as references; the three cited here are quartz–iron–fayalite (QIF), which is Fe(II)/Fe(0), iron/magnetite (IM), which is Fe3O4/Fe(0), and fayalite–magnetite–quartz (FMQ), which is Fe3O4/Fe(II). The redox state of a system may be represented by its oxygen fugacity f(O2) [28], but we note that the f(O2) of a mineral assemblage changes with pressure and temperature, even if the elemental abundance in that assemblage does not change (see [28] for review).
The redox state of the Venusian lithosphere is not known, but is likely to be FMQ for two reasons. Firstly, all in situ measurements of the elemental composition of Venus’ surface are consistent with it being similar to terrestrial basalt [53,54,55], which itself typically varies between WM and FMQ oxygen fugacity. Secondly, SO2 is expected to react with surface rocks on Venus, and so, the presence of 30–300 ppm SO2 below the Venusian clouds [14] requires the replenishment of the SO2 by volcanism [56,57]. The <15 ppb of H2S [58] in the atmosphere suggests that the source of sulfur gases is oxidized relative to the H2S/SO2 couple, which requires the O2 fugacity of FMQ or WM depending on assumptions about water content and temperature [11]. Under Venus surface pressure (~90 bar), H2S is slightly more abundant for any given temperature and water content than it would be on Earth (Ref. [11] and the S.I. Section 2.4.2 of Ref. [11]), which reinforces the idea that to achieve dominance of SO2, the Venusian surface must be oxidized. Again, this does not preclude local pockets of extremely reduced (or extremely oxidized) rocks, but it does argue strongly that, as a global average, volcanic outgassing is from rocks in the WM—FMQ range of oxygen fugacity.
In principle, phosphorus could be present in iron or silicate phases of the upper mantle, and hence, in erupted magmas. We therefore calculated the partitioning of phosphorus between metal and silicate, which depends on pressure and temperature. The result in Figure 2 shows that down to 600 km depth phosphorus is overwhelmingly partitioned into the silicate phase of melts. We therefore focus on the thermodynamics of phosphorus in the silicate phase.
We calculated the oxygen fugacity of the phosphate/phosphide equilibrium for all couples for iron and for all more electropositive metals than iron for phosphate and phosphide, for which thermodynamic data were available for the temperature range in the lithosphere and upper mantle. Other metals such as nickel and copper that are more electronegative than ison were not considered. If copper phosphide was present in an oxidizing environment (FMQ), it would be oxidized to copper phosphate; if present in a reduced environment (QIF), it would be displaced by more electropositive iron to form iron phosphide and copper metal. This is illustrated by the observation that meteoritic phosphide is found overwhelmingly as Fe3P (Schreibersite) and only a rare, trace mineral as Ni3P (Rhabdite). We emphasize that this does not mean that copper phosphide or other metal phosphides could not occur, just that they are unlikely to be the major phosphides in the overall shallow mantle minerology, and hence are unlikely to play a dominant role in the overall minerology of volcanism on Earth or on Venus.
The results of the thermodynamic analysis are shown in Figure 3.
The calculation on Figure 3 clearly shows that all of the phosphates for which thermodynamic data are available are preferred over phosphides at depths shallower than 200 km. If the mantle is reduced at QIF/IM levels, then iron phosphides become preferred over phosphates at a depth of ~850 km, where the temperature is ~2000 K, which is above the melting temperature of iron phosphides [59]. This is consistent with iron phosphides being dredged from the continental ‘keels’ of continental cratons on Earth [60] as xenoliths. However, such xenoliths represent milligrams of material in kilograms of rock, not the bulk melt [52,61,62,63,64,65,66] (hence the name xenolith—“Foreign rock”). At lower mantle depths (Figure 3A), several phosphides are strongly favored over phosphate, consistent with such minerals being found as inclusions in mantle diamonds.
However, no phosphide is preferred over phosphate in the upper 200 km of the lithosphere. This is true whether we assume a melt temperature of 1100 K, or of 660 K (the ambient temperature at the top of the highest mountain on Venus [67]). This supports our previous calculations, which showed that it would require extraordinary circumstances for phosphorus to be present as phosphine or elemental phosphorus in near-surface rocks [11]. For an 1100 K melt and a QIF surface rock redox, and if Fe3P were a dominant phosphide, then 8·10−6 of the phosphorus in the rock would be present as phosphide. (Barringerite—Fe2P—has only been found in meteorites and in pyrometamorphic rocks [68,69], and is an improbable phosphide to form in large scale volcanic melts, as iron is 40 times as abundant in the lithosphere of Earth as phosphorus, and a similar abundance ratio is expected on Venus). If the lithosphere is at FMQ redox state, and Fe3P is the most likely phosphide, then phosphide will constitute only 5·10−11 of the erupted phosphorus. We emphasize that it is the minerology of these near surface rocks, not the minerology of the deep mantle, that defines the minerology of volcanic erupta. We can therefore say that, even if the crust has an oxygen fugacity of QIF/IM (which, as mentioned above, is unlikely), phosphides will comprise only a very small fraction of erupted phosphorus. Therefore, to produce phosphine from plume volcanism, a substantially larger amount of rock would have to erupt than is assumed by [16]. We present a quantitative calculation of how much magma would have to erupt in Section 4 below.
In summary, phosphides might be present as <10−5 of the phosphorus in surface (<200 km depth) melts if the melt is at the QIF redox buffer, or if they are transported from ~200 km depth as rare xenoliths. In melts at the FMQ buffer redox, which is more characteristic of terrestrial melts, phosphide would be no more than 10−15 of the phosphorus content.
We conclude that the optimistic scenarios suggest that melts erupting on the surface could contain 10−5–10−17 of their phosphorus content as phosphide.

3.3. Mechanisms of Explosive Volcanism on Venus

Having reached the surface, phosphide minerals have to interact with acidic liquids to generate phosphine efficiently. This means that they must be fragmented into fine ash or dust and then lifted to heights of >35 km. Fragmentation is needed for two reasons; to allow the material to be lifted to high altitude and to allow it to react with sulfuric acid efficiently. Truong and Lunine’s model [16] postulates that explosive volcanism lifts erupted material tens of kilometers into the Venusian atmosphere so that it can react with sulfuric acid to make PH3. The most obvious site of the reaction is the sulfuric acid droplets in the clouds, but in principle phosphides, could also react with sulfuric acid vapor below the clouds, down to 35 km, below which H2SO4 is predominantly dissociated from H2O and SO3 [70]. Truong and Lunine propose that ash is lifted by explosive volcanism, and use the example of the 1883 eruption of Krakatoa, which lifted several cubic miles of ash and dust into the Earth’s stratosphere, as a model for such an explosive eruption. This type of volcanism contrasts with effusive eruptions, where the product is primarily lava flows which remain on the ground. We therefore next discuss whether explosive volcanism of the Krakatoa type is a realistic model for Venusian volcanism. We note that, even on Earth, explosive volcanism rarely generates ash columns that reach above 10 km. However, explosive volcanism is not likely to be common on Venus.
Explosive volcanism in this case is not just volcanism that provides a gas column that can rise to stratospheric heights, but one that also fragments a substantial fraction of the erupted magma into fine dust and ash. The observation of tropospheric SO2 by Pioneer Venus [71,72] and the requirement for a source of SO2 to balance reaction with crustal rocks [56,57] suggest that current volcanism is erupting SO2 into Venus’ atmosphere (although this itself is controversial), but is not evidence that those eruptions produce an ash column.
Explosive volcanism can only occur if erupted magma has a high content of a volatile, which drives the fragmentation of the magma, as well as the physical explosive lifting of the material into the ash cloud. On Earth, by far the most common volatile driver of explosive volcanism is water. Systems that provide explosive eruptions and ash columns are overwhelmingly seen at subduction zones or where surface water is able to interact with ascending magma (e.g., in Iceland, where volcanoes may erupt beneath the ice). This is because the subduction of the oceanic seafloor traps water, sulfates, and organics, which provide the high volatile content of resulting magmas. The 1992 Cerro Negro eruption cited by [16] is a classic example of a basaltic eruption that is nevertheless part of a subduction system [73]; in this case, the Central American Volcanic Arc resulting from the subduction zone at the western edge of the Caribbean Plate [74]. The Tarawera eruptions are also the product of subduction arc volcanism. Tarawera is part of the Taupo volcanic zone, which is subduction zone volcanism [75,76,77]. The 1883 eruption of Krakatoa, an unusually violent event (described as “the loudest sound heard in history” [78]) was not only the product of a subduction zone volcano but also of ocean–magma interactions which drove the final, catastrophic explosion. Similarly, the explosive eruption of the Reunion plume volcano was due to water: ocean interaction [79]. Other historically destructive explosions, such as the Minoan explosion of Thera/Santorini [80], also probably had a component of ocean–magma interactions to drive their unusually powerful eruptions and the consequent injection of ash into the high atmosphere.
Melt viscosity plays a major role in explosive volcanism, and most viscous melts on Earth, e.g., rhyolites, come from processing subducted crust [81]; the Tarawera eruptions were viscous and rhyolite-rich. However, whatever the geology, without a high content of volatiles, explosive eruption is not possible. This poses two problems for postulating explosive volcanism on Venus.
First, Venus is notoriously water-poor. Venus probably does not have Earth-like subduction zones, and certainly does not have surface oceans which could provide subducted water or interact with erupting magma. The net rate of loss of H from Venus is 1 to 2 orders of magnitude less than that from Earth [82], so the volcanic input of H into the atmosphere (as H2O or any other gas) in total must be no more than 10% of Earth’s, and water-driven explosive volcanism correspondingly less frequent. Modelling suggests that explosive volcanism on Venus might occur if the water content of the melts was between 2% and 5%, depending on altitude, conduit geometry, and melt minerology [67,83]. Such water contents are characteristic of oceanic subduction zones on Earth [84,85,86]; the hotspot volcanism, which might be a better analogue for Venusian magmatism, typically has water content <0.5% [87,88]. Thus, explosive volcanism on Venus driven by water is unlikely in average geology, and hence will be rare.
The second issue with explosive volcanism as a mechanism for generating phosphide-containing ash is the consistency of having phosphide and water present at the same time in the magma. If explosive volcanism is driven by water degassing from the melt, then water must be present in the melt. Although the reaction of phosphide with water is 1000-fold slower than the reaction with acid [89] at room temperature, it is expected to be 1030 times faster at 1300 K than at 300 K, assuming Q10 = 2. Any phosphide in the melt would therefore react with water in the melt to form PH3, which is thermally unstable at melt temperatures [90]. Again, the primary driver of explosive volcanism on Earth is unlikely to deliver phosphides in large amounts.
Explosive volcanism therefore would have to be driven by a volatile other than water. Explosive volcanism could be driven by CO2 (e.g., potentially the 1085 eruption of Sunset Crater, Flagstaff [91]) or SO2 (e.g., the 1982 El Chichon eruption [92]—although this was a subduction volcano). This is, however, also inconsistent with the presence of phosphide. Either the bulk melt is sufficiently reduced to preserve P as phosphide, in which case C would be present as carbide or methane and S as sulfide or H2S, or the melt is sufficiently oxidized for C to be present as CO2 and S as SO2; in which case, the P (which is more electropositive than C and S, i.e., more readily oxidized) would be present as P(+5) oxides or anions. We also note that SO2 is less likely to drive explosive volcanism on Venus than on Earth, as the higher surface pressure reduces outgassing of SO2 (as well as H2S) (See Figure 3 in ref. [83]).
(These considerations do not affect xenolith eruption, as xenoliths are chemically isolated from the melt on the melt transport timescale, but as noted above, xenoliths form an extremely small fraction of a magma.)
Chemical consistency therefore requires a melt that is very low in water but contains at least 5% of another volatile that is not CO2 or SO2. While such chemistry could be imagined, it would be unprecedented.
This is confirmed by the radar imaging of Venus’ surface. Imaging shows very little evidence of explosive volcanism on Venus, such as pyroclastic deposits, and probable volcanic features are consistent with an overwhelmingly effusive mode of volcanism [93,94,95]. The same studies suggest that Venus has not undergone exceptional volcanic activity in the recent past. We conclude that explosive volcanism is both expected to be rare, and is observed to be rare, on Venus.
The alternative is that the ash clouds generated by effusive volcanism are sufficiently large to drive the required phosphide load to 35 km. Effusive volcanism does generate ash and dust, but several orders of magnitude less than explosive eruption. Studies of a number of eruptions of mantle hot-spot volcanoes show that they can produce ash clouds of several kilometers high, but that these contain 0.1–1% of the total erupted volume at most, usually much less; up to 10% of the erupted mass is produced as non-effusive material, but this reaches only tens or hundreds of meters into the air (the ‘lava fountains’ that characterize such eruptions) and only 2–3% or less of the ash ejecta reaches heights of 1 km or more (e.g., [96,97,98]). It is plausible to suggest that a much larger volcano could loft some of the plume-derived material to stratospheric heights through entrainment in the rising gas column, a possibility we address in Section 3.4 below.
We conclude that volcanism on Venus is likely to primarily be effusive, which is no more than 1–10% as efficient at delivering fine ash to high altitudes as explosive eruptions, and probably would deliver only trivial amounts.

3.4. Ash Plume Generation on Venus

Several factors mitigate against any volcanic eruption on Venus being as efficient at delivering ash to >35 km altitudes as it would be on Earth, regardless of its minerology. These are discussed in detail in [95]. Specifically, ash generation is less efficient, and plumes are likely to be smaller.
To generate fine ash, magma must be fragmented from bulk material into very small fragments. Regardless of magma chemistry, high surface pressure will suppress ash generation. At higher surface barometric pressure, gases will exsolve less and into smaller volumes, which will suppress rapid bubble expansion, leading to lower degrees of fragmentation.
Eruption itself only projects magma a few hundred meters into the atmosphere. To rise further, ash must be entrained in a rising column of hot gas. The lift driving volcanic plume ascent into planetary atmospheres comes from heat exchange from the fragmented magma to the atmosphere, heating the gas, and thereby producing buoyancy. The efficiency of this heat exchange is a function of magma fragmentation, which, as noted above, is likely to be lower on Venus than on Earth. The resulting buoyant driving force will depend in part on the temperature lapse rate of the atmosphere, which is less on Venus than on Earth [72]. Head and Wilson [95] estimates that plumes on Venus will be ~1/3 the height those generated by similar eruptions on Earth. Extremely high volatile fractions (>12% SO2 for example) are needed to create a 50 km height ash cloud under Venus conditions [95].
Thus, although we adopt a simple scaling that effusive volcanism produces 10–100 times less ash than explosive volcanism, the above considerations point to this being highly conservative in the case of Venus, where the real decrease in fine ash production may be more like >1000 times less than for a terrestrial explosive eruption.

3.5. Other Criteria for a Volcanic Source

We note that, in the literature, data on the hydrolysis of phosphides cited by Truong and Lunine relate to aqueous solutions of sulfuric acid, not concentrated acid. Concentrated sulfuric acid is an oxidizing agent, not a hydrolyzing reagent; noted in [11]. We also believe that Truong and Lunine overestimate the mantle phosphorus content. We present isotope record data in [11] that suggests that, outside very limited regions, the mantle phosphorus abundance is likely to be ~0.2 wt.%. Further analysis of >2000 rock samples from Earthchem [99] suggests an average phosphorus abundance of 0.074 wt.% in terrestrial igneous rock samples.

3.6. Required Rate of Phosphide Volcanic Eruption

We will briefly address the amount of mineral that is required to produce the flux of phosphine necessary to maintain 1 ppb gas at cloud level in the atmosphere. The required PH3 production rate has been dealt with in detail in [11], so we will only summarize the argument here.
Ref. [11] provides a detailed photochemical analysis of the rate of PH3 production necessary to sustain a 1 ppb PH3 level in the clouds. Ref. [16] simplifies this, and estimates the phosphine production rate required to explain the putative phosphine signal by
R P = M P H 3 t P H 3 ,
where M P H 3 is the mass of phosphine in the 53–61 km detection layer reported by [1] and t P H 3 is the photochemical lifetime of PH3 at 53–61 km. Ref. [16] further estimates t P H 3 by the transport timescale for a molecule to cross ~10 km, which is in the order of magnitude of the distance a molecule, would need to traverse to reach the bottom or top of the Venusian clouds, where its lifetime is short due to photochemical destruction [1].
The ref. [16] calculation underestimates the required production flux of phosphine, for three reasons:
  • It assumes that phosphine is only present from 53–61 km. However, the very same diffusive processes invoked to transport phosphine to the upper atmospheres mean that phosphine must be present at altitudes beyond 53–61 km as well. In other words, this approach underestimates M P H 3 .
  • It overestimates the lifetime of high-altitude phosphine. The eddy diffusion coefficient increases with altitude in the Venusian atmosphere, meaning that the use of 103 cm2 s−1 underestimates the phosphine destruction rate. Therefore, it overestimates t P H 3 for PH3 at z > 63 km, which must exist due to the diffusive transport.
  • It assumes that phosphine is only destroyed after transport. In fact, phosphine is continually destroyed throughout the Venusian atmosphere due to radical sinks [1].
When we attempt to account for these factors by running a full 1D photochemical model, we find a required PH3 production rate of 26 kg/s, compared to 1 kg/s for [16] (to produce 1 ppb PH3). We note that this probably under-estimates the necessary PH3 production rate, as PH3 is oxidized efficiently by liquid concentrated sulfuric acid [13], a sink not accounted for in any models for lack of kinetic data.
This is not to discount the usefulness of such calculation; indeed, we executed a variant of it ourselves ([11] and Section 3.2.4.1 in ref. [11]). This calculation is valuable because it involves few parameters and is therefore less vulnerable to the considerable photochemical uncertainties. Further, such calculation has the philosophical advantage that it seeks to explain solely the detected phosphine, and thus, in a sense, provides the most conservative approach (easiest for known mechanisms to explain). However, it has the disadvantage that, in order to be true, it requires transport in the Venusian atmosphere to behave in ways atmospheric transport is not known to behave. Especially, it requires PH3 not to exist outside of the ~53–61 km layer, while the same diffusion which transports the PH3 to be destroyed should also result in PH3 existing outside the 53–61 km layer. For the [16] explanation to function, Venus must host a unique atmospheric vertical transport pattern which sequesters PH3 in the cloud deck.
For the purposes of this paper, we adopt a range of fluxes; a lower limit of the 1 kg/s adopted by Truong and Lunine [16] as a conservative value, and a higher limit of the [11] value of 26 kg/s.

4. Discussion

Truong and Lunine [16] revisited the idea of mantle plume volcanism as a source of the PH3 tentatively detected on Venus. We have built on this hypothesis quantitatively in terms of atmospheric chemistry, mantle minerology, and eruption physics, and chemistry. As a base case, [16] estimate that 1 kg/sec of phosphides needs to be delivered to clouds, and this requires 0.03–0.15 km3/year of material containing 1% phosphorus, 5% of this is phosphides, and a t PH 3 = 109 s. In Table 2, we summarize the factors needed to correct for a more complete model of phosphorus geochemistry, assuming that the lithosphere is reducing (QIF oxygen fugacity). We conclude that the actual amount of magma that needs to erupt is ≥21,600 km3/year. Even the Deccan and Siberian trap emplacements (which might have been up to 10 km3/year at their peak [16,100,101]), fall short of the required rate of phosphide injection by at least three orders of magnitude. This therefore represents an unprecedented level of volcanic activity.
21,600 cubic kilometers is quite a lot of magma; enough to build a volcano the size of Olympus Mons every 138 years [102], or to resurface the planet to a depth of 1 km in 22,000 years. It clearly cannot be from one volcanic plume; if this were a single volcanic conduit through which magma was flowing at 10 m/sec (an exceptionally high velocity for terrestrial volcanoes [103,104]), then such a flux would require a ‘dike’ with an area of ~68,500 square meters, or 12 American Football fields, an unprecedented event. The result summarized in Table 2 assumes an extremely conservative model of surface geochemistry, volcanic mechanics, and of phosphine photochemistry. If we instead adopt the more physically realistic photochemical model of [1,11], then the rate of phosphide delivery to the clouds is 26 times that assumed here, and the volcanic eruption rate is similarly higher at ~560,000 km3/year. If we assume that the lithosphere has an oxygen fugacity nearer FMQ than QIF, then the required eruption rate is several orders of magnitude higher still. If we adopt a more realistic rate of ash cloud generation, then the required eruption rate is again increased.
We conclude that large-scale volcanism is a plausible source for PH3 on Venus if
(a)
The atmospheric structure and gas transport on Venus are such that PH3 (and only PH3) is retained solely in the lower cloud layer, and does not diffuse upwards or downwards, and
(b)
The Venusian lithosphere is substantially more reduced than is expected from lander data and atmospheric chemistry, and
(c)
Volcanism is currently erupting tens of thousands of cubic kilometers of magma onto the surface per year, and
(d)
Venusian volcanism is unexpectedly efficient at generating high altitude clouds of fine magma ash.
None of these requirements is unphysical, but all are unexpected and can be tested by observation. We note that the unexpected atmospheric physics can only apply to PH3, and not to other minor gases (which are adequately modelled using conventional mixing processes).
Very high volumes of volcanic eruption should have altered the surface over hundreds of square kilometers since radar observations of Venus began in the 1970s, so systematic searches for such changes could test for the reality of such a very large-scale volcanism. Additionally, the thermal imprint of such large-scale volcanism would be expected to be visible in the cloud layer as a hot spot whose source co-rotated with the planet and not with the clouds.
There remains the possibility that very high volumes of magma happened to be erupting in 1978 (when Pioneer Venus sampled the atmosphere), and in 2017–2019, when the Greaves et al. [1] observations were made, but eruptions have been quiescent for most of the intervening 40 years. This would require only one or two volcanoes the size of the Hawaiian island complex [105] to be accumulated since Venus’ surface observations began under our most optimistic assumptions, but would still point to Venus entering a new phase of volcanism unlike that of the recent past.

5. Conclusions

The presence of PH3 in Venus’ atmosphere is unexpected and invites explanation. Volcanism is only an adequate explanation if our understanding of basic diffusion processes in Venus’ atmosphere, the minerology of Venus’ surface, and of the mechanisms of ash cloud formation are incorrect, and if Venus is currently undergoing a major eruption event of a magnitude greater than any seen on Earth in the Phanerozoic, and different from the volcanism regime in Venus’ recent past. It has been suggested that Venus has undergone such major resurfacing as global resurfacing ~300 Mya, in possibly a single catastrophic event that lasted less than 10 My [17]. If Venus were starting such an event now, it could explain the presence of phosphine in the clouds. Such an exciting possibility should be explored.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/universe8010054/s1: Spreadsheet containing the following sheets: Buffers (table of standard buffer fugacity as a function of depth), Fugacity Calculation (table of free energy of formation of minerals and calculation of oxygen fugacity for selected phosphide/phosphate buffers), References (literature references for the thermodynamic data), Read Me (details of how the calculations are performed in this spreadsheet).

Author Contributions

Conceptualization, W.B. and O.S.; methodology, W.B., O.S. and S.R.; writing—original draft preparation, W.B.; writing—review and editing, W.B., O.S., S.R., P.B.R., J.J.P., J.S.G. and S.S.; funding acquisition, S.S., P.B.R. and O.S. All authors have read and agreed to the published version of the manuscript.

Funding

We thank the Change Happens Foundation for funding. S.R. acknowledges the funding from the Simons Foundation (495062). P.B.R. thanks the Simons Foundation for funding (SCOL awards 599634).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All relevant data is provided in Supplementary Materials to this article.

Acknowledgments

We are very grateful for the constructive comments provided in review.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Greaves, J.S.; Richards, A.M.; Bains, W.; Rimmer, P.B.; Sagawa, H.; Clements, D.L.; Seager, S.; Petkowski, J.J.; Sousa-Silva, C.; Ranjan, S. Phosphine gas in the cloud decks of Venus. Nat. Astron. 2020, 5, 655–664. [Google Scholar] [CrossRef]
  2. Greaves, J.S.; Richards, A.M.S.; Bains, W.; Rimmer, P.B.; Sagawa, H.; Clements, D.L.; Seager, S.; Petkowski, J.J.; Sousa-Silva, C.; Ranjan, S.; et al. Addendum: Phosphine gas in the cloud deck of Venus. Nat. Astron. 2021, 5, 726–728. [Google Scholar] [CrossRef]
  3. Trompet, L.; Robert, S.; Mahieux, A.; Schmidt, F.; Erwin, J.; Vandaele, A. Phosphine in Venus’ atmosphere: Detection attempts and upper limits above the cloud top assessed from the SOIR/VEx spectra. Astron. Astrophys. 2021, 645, L4. [Google Scholar] [CrossRef]
  4. Villanueva, G.L.; Cordiner, M.; Irwin, P.G.J.; de Pater, I.; Butler, B.; Gurwell, M.; Milam, S.N.; Nixon, C.A.; Luszcz-Cook, S.H.; Wilson, C.F.; et al. No evidence of phosphine in the atmosphere of Venus from independent analyses. Nat. Astron. 2021, 5, 631–635. [Google Scholar] [CrossRef]
  5. Akins, A.B.; Lincowski, A.P.; Meadows, V.S.; Steffes, P.G. Complications in the ALMA detection of phosphine at Venus. Astrophys. J. 2021, 907, L27. [Google Scholar] [CrossRef]
  6. Lincowski, A.P.; Meadows, V.S.; Crisp, D.; Akins, A.B.; Schwieterman, E.W.; Arney, G.N.; Wong, M.L.; Steffes, P.G.; Parenteau, M.N.; Domagal-Goldman, S. Claimed detection of PH3 in the clouds of Venus is consistent with mesospheric SO2. Astrophys. J. Lett. 2021, 908, L44. [Google Scholar] [CrossRef]
  7. Greaves, J.S.; Rimmer, P.B.; Richards, A.; Petkowski, J.J.; Bains, W.; Ranjan, S.; Seager, S.; Clements, D.L.; Silva, C.S.; Fraser, H.J. Low Levels of sulphur dioxide contamination of phosphine spectra from Venus’ atmosphere. arXiv 2021, arXiv:2108.08393. [Google Scholar]
  8. Greaves, J.S.; Richards, A.M.; Bains, W.; Rimmer, P.B.; Clements, D.L.; Seager, S.; Petkowski, J.J.; Sousa-Silva, C.; Ranjan, S.; Fraser, H.J. Reply to: No evidence of phosphine in the atmosphere of Venus from independent analyses. Nat. Astron. 2021, 5, 636–639. [Google Scholar] [CrossRef]
  9. Mogul, R.; Limaye, S.S.; Way, M.J.; Cordova, J.A. Venus’ mass spectra show signs of disequilibria in the middle clouds. Geophys. Res. Lett. 2021, 48, e2020GL091327. [Google Scholar] [CrossRef]
  10. Omran, A.; Oze, C.; Jackson, B.; Mehta, C.; Barge, L.M.; Bada, J.; Pasek, M.A. Phosphine generation pathways on rocky planets. Astrobiology 2021, 21, 1264–1276. [Google Scholar] [CrossRef]
  11. Bains, W.; Petkowski, J.J.; Seager, S.; Ranjan, S.; Sousa-Silva, C.; Rimmer, P.B.; Zhan, Z.; Greaves, J.S.; Richards, A.M. Phosphine on Venus cannot be explained by conventional processes. Astrobiology 2021, 21, 1277–1304. [Google Scholar] [CrossRef] [PubMed]
  12. Bains, W.; Petkowski, J.J.; Seager, S.; Ranjan, S.; Sousa-Silva, C.; Rimmer, P.B.; Zhan, Z.; Greaves, J.S.; Richards, A.M. Venusian phosphine: A ‘wow!’signal in chemistry? In Phosphorus, Sulfur, and Silicon and the Related Elements; Taylor & Francis: London, UK, 2021; pp. 1–6. [Google Scholar]
  13. Bains, W.; Petkowski, J.J.; Zhan, Z.; Seager, S. Evaluating alternatives to water as solvents for life: The example of sulfuric acid. Life 2021, 11, 400. [Google Scholar] [CrossRef] [PubMed]
  14. Rimmer, P.B.; Jordan, S.; Constantinou, T.; Woitke, P.; Shorttle, O.; Hobbs, R.; Paschodimas, A. Hydroxide salts in the clouds of Venus: Their effect on the sulfur cycle and cloud droplet pH. Planet. Sci. J. 2021, 2, 133. [Google Scholar] [CrossRef]
  15. Bains, W.; Petkowski, J.J.; Rimmer, P.B.; Seager, S. Production of ammonia makes venusian clouds habitable and explains observed cloud-level chemical anomalies. Proc. Natl. Acad. Sci. USA 2021, 118, e2110889118. [Google Scholar] [CrossRef]
  16. Truong, N.; Lunine, J. Volcanically extruded phosphides as an abiotic source of Venusian phosphine. Proc. Natl. Acad. Sci. USA 2021, 118, e2021689118. [Google Scholar] [CrossRef]
  17. Strom, R.G.; Schaber, G.G.; Dawson, D.D. The global resurfacing of Venus. J. Geophys. Res. Planets 1994, 99, 10899–10926. [Google Scholar] [CrossRef]
  18. Gu, T.; Stagno, V.; Fei, Y. Partition coefficient of phosphorus between liquid metal and silicate melt with implications for the Martian magma ocean. Phys. Earth Planet. Inter. 2019, 295, 106298. [Google Scholar] [CrossRef]
  19. Dziewonski, A.M.; Anderson, D.L. Preliminary reference Earth model. Phys. Earth Planet. Inter. 1981, 25, 297–356. [Google Scholar] [CrossRef]
  20. Spiliopoulos, S.; Stacey, F.D. The earth’s thermal profile: Is there a mid-mantle thermal boundary layer? J. Geodyn. 1984, 1, 61–77. [Google Scholar] [CrossRef]
  21. Dannberg, J.; Sobolev, S.V. Low-buoyancy thermochemical plumes resolve controversy of classical mantle plume concept. Nat. Commun. 2015, 6, 6960. [Google Scholar] [CrossRef] [Green Version]
  22. Chase, M.W.J. NIST-JANAF thermochemical tables. J. Phys. Chem. Ref. Data 1998, 9. [Google Scholar] [CrossRef]
  23. Robie, R.A.; Hemingway, B.S. Thermodynamic Properties of Minerals and Related Substances at 298.15 K and 1 bar (105 Pascals) Pressure and at Higher Temperatures; US Government Printing Office: Washington, DC, USA, 1995.
  24. Zaitsev, A.; Dobrokhotova, Z.V.; Litvina, A.; Mogutnov, B. Thermodynamic properties and phase equilibria in the Fe–P system. J. Chem. Soc. Faraday Trans. 1995, 91, 703–712. [Google Scholar] [CrossRef]
  25. Barin, I. (Ed.) Thermochemical Data of Pure Substances, 3rd ed.; VCH Verlagsgesellschaft mbH: Berlin, Germany, 1995. [Google Scholar]
  26. Schlesinger, M.E. The Thermodynamic Properties of Phosphorus and Solid Binary Phosphides. Chem. Rev. 2002, 102, 4267–4302. [Google Scholar] [CrossRef] [PubMed]
  27. Jain, A.; Ong, S.P.; Hautier, G.; Chen, W.; Richards, W.D.; Dacek, S.; Cholia, S.; Gunter, D.; Skinner, D.; Ceder, G. Commentary: The Materials Project: A materials genome approach to accelerating materials innovation. APL Mater. 2013, 1, 011002. [Google Scholar] [CrossRef] [Green Version]
  28. Frost, B.R. Introduction to oxygen fugacity and its petrologic significance. In Reviews in Minerology Volume 25. Oxide Minerals: Petrologic and Magnetic Significance; Lindsleys, D.H., Banerjee, S.K., Eds.; Minerological Society of America: Washington, DC, USA, 1991; pp. 1–10. [Google Scholar]
  29. Uemura, K. Diffusion of elements in solid iron and steel. Tetsu-to-Hagane 1943, 29, 329–335. [Google Scholar] [CrossRef] [Green Version]
  30. Wang, X.; Posselt, M.; Faßbender, J. Influence of substitutional atoms on the diffusion of oxygen in dilute iron alloys. Phys. Rev. B 2018, 98, 064103. [Google Scholar] [CrossRef]
  31. Stolper, E.M.; Shorttle, O.; Antoshechkina, P.M.; Asimow, P.D. The effects of solid-solid phase equilibria on the oxygen fugacity of the upper mantle. Am. Mineral. J. Earth Planet. Mater. 2020, 105, 1445–1471. [Google Scholar] [CrossRef]
  32. Griffin, W.; Gain, S.; Cámara, F.; Bindi, L.; Shaw, J.; Alard, O.; Saunders, M.; Huang, J.-X.; Toledo, V.; O’Reilly, S. Extreme reduction: Mantle-derived oxide xenoliths from a hydrogen-rich environment. Lithos 2020, 358, 105404. [Google Scholar] [CrossRef]
  33. Choudhuri, M.; Nemčok, M. Plumes and Hotspots. In Mantle Plumes and Their Effects; Choudhuri, M., Nemčok, M., Eds.; Springer International Publishing: Cham, Switzerland, 2017; pp. 19–42. [Google Scholar]
  34. Sleep, N.H. Hotspot volcanism and mantle plumes. Ann. Rev. Earth Planet. Sci. 1992, 20, 19–43. [Google Scholar] [CrossRef]
  35. Loper, D.E. Mantle plumes and their effect on the Earth’s surface: A review and synthesis. Dyn. Atmos. Ocean. 1998, 27, 35–54. [Google Scholar] [CrossRef]
  36. Sleep, N.H.; Ebinger, C.J.; Kendall, J.-M. Deflection of mantle plume material by cratonic keels. Geol. Soc. Lond. Spec. Publ. 2002, 199, 135–150. [Google Scholar] [CrossRef]
  37. Williams, Q.; Jeanloz, R.; Bass, J.; Svendsen, B.; Ahrens, T.J. The melting curve of iron to 250 gigapascals: A constraint on the temperature at earth’s center. Science 1987, 236, 181–182. [Google Scholar] [CrossRef] [PubMed]
  38. d’Acremont, E.; Leroy, S.; Burov, E.B. Numerical modelling of a mantle plume: The plume head–lithosphere interaction in the formation of an oceanic large igneous province. Earth Planet. Sci. Lett. 2003, 206, 379–396. [Google Scholar] [CrossRef]
  39. Smrekar, S.; Auerbach, V.; Ostberg, C.; O’Rourke, J. How thick is the lithosphere of Venus? LPI Contrib. 2019, 2193, 8035. [Google Scholar]
  40. Uppalapati, S.; Rolf, T.; Crameri, F.; Werner, S.C. Dynamics of lithospheric overturns and implications for Venus’s surface. J. Geophys. Res. Planets 2020, 125, e2019JE006258. [Google Scholar] [CrossRef]
  41. Behn, M.D.; Hirth, G.; Elsenbeck, J.R. Implications of grain size evolution on the seismic structure of the oceanic upper mantle. Earth Planet. Sci. Lett. 2009, 282, 178–189. [Google Scholar] [CrossRef] [Green Version]
  42. Dannberg, J.; Eilon, Z.; Faul, U.; Gassmöller, R.; Moulik, P.; Myhill, R. The importance of grain size to mantle dynamics and seismological observations. Geochem. Geophys. Geosyst. 2017, 18, 3034–3061. [Google Scholar] [CrossRef]
  43. Mazza, S.E.; Gazel, E.; Bizimis, M.; Moucha, R.; Béguelin, P.; Johnson, E.A.; McAleer, R.J.; Sobolev, A.V. Sampling the volatile-rich transition zone beneath Bermuda. Nature 2019, 569, 398–403. [Google Scholar] [CrossRef]
  44. Beyer, C.; Frost, D.J. The depth of sub-lithospheric diamond formation and the redistribution of carbon in the deep mantle. Earth Planet. Sci. Lett. 2017, 461, 30–39. [Google Scholar] [CrossRef]
  45. Kaminsky, F.V.; Wirth, R.; Schreiber, A. Carbonatitic inclusions in deep mantle diamond from Juina, Brazil: New minerals in the carbonate-halide association. Can. Mineral. 2013, 51, 669–688. [Google Scholar] [CrossRef]
  46. Palot, M.; Pearson, D.G.; Stern, R.A.; Stachel, T.; Harris, J.W. Isotopic constraints on the nature and circulation of deep mantle C–H–O–N fluids: Carbon and nitrogen systematics within ultra-deep diamonds from Kankan (Guinea). Geochim. Cosmochim. Acta 2014, 139, 26–46. [Google Scholar] [CrossRef]
  47. Pearson, D.; Brenker, F.; Nestola, F.; McNeill, J.; Nasdala, L.; Hutchison, M.; Matveev, S.; Mather, K.; Silversmit, G.; Schmitz, S. Hydrous mantle transition zone indicated by ringwoodite included within diamond. Nature 2014, 507, 221–224. [Google Scholar] [CrossRef]
  48. Smith, E.M.; Shirey, S.B.; Nestola, F.; Bullock, E.S.; Wang, J.; Richardson, S.H.; Wang, W. Large gem diamonds from metallic liquid in Earth’s deep mantle. Science 2016, 354, 1403–1405. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Tschauner, O.; Huang, S.; Greenberg, E.; Prakapenka, V.; Ma, C.; Rossman, G.; Shen, A.; Zhang, D.; Newville, M.; Lanzirotti, A. Ice-VII inclusions in diamonds: Evidence for aqueous fluid in Earth’s deep mantle. Science 2018, 359, 1136–1139. [Google Scholar] [CrossRef] [Green Version]
  50. Scott, E.R. Iron meteorites: Composition, age, and origin. In Oxford Research Encyclopedia of Planetary Science; Oxford University Press: Oxford, UK, 2020. [Google Scholar]
  51. Ulff-Moller, F. Solidification history of the Kitdlit Lens: Immiscible metal and sulphide liquids from a basaltic dyke on Disko, central West Greenland. J. Petrol. 1985, 26, 64–91. [Google Scholar] [CrossRef]
  52. Xiong, Q.; Griffin, W.L.; Huang, J.-X.; Gain, S.E.; Toledo, V.; Pearson, N.J.; O’Reilly, S.Y. Super-reduced mineral assemblages in “ophiolitic” chromitites and peridotites: The view from Mount Carmel. Eur. J. Mineral. 2017, 29, 557–570. [Google Scholar] [CrossRef]
  53. Surkov, Y.A.; Barsukov, V.L.; Moskalyeva, L.P.; Kharyukova, V.P.; Kemurdzhian, A.L. New data on the composition, structure, and properties of Venus rock obtained by Venera 13 and Venera 14. J. Geophys. Res. Solid Earth 1984, 89, B393–B402. [Google Scholar] [CrossRef]
  54. Surkov, Y.A.; Moskalyova, L.P.; Kharyukova, V.P.; Dudin, A.D.; Smirnov, G.G.; Zaitseva, S.Y. Venus rock composition at the Vega 2 Landing Site. J. Geophys. Res. Solid Earth 1986, 91, E215–E218. [Google Scholar] [CrossRef]
  55. Gilmore, M.S.; Mueller, N.; Helbert, J. VIRTIS emissivity of Alpha Regio, Venus, with implications for tessera composition. Icarus 2015, 254, 350–361. [Google Scholar] [CrossRef]
  56. Fegley, B. 2.7–Venus. In Treatise on Geochemistry, 2nd ed.; Holland, H.D., Turekian, K.K., Eds.; Elsevier: Oxford, UK, 2014; pp. 127–148. [Google Scholar]
  57. Fegley, B.; Prinn, R.G. Estimation of the rate of volcanism on Venus from reaction rate measurements. Nature 1989, 337, 55–58. [Google Scholar] [CrossRef]
  58. Krasnopolsky, V.A. High-resolution spectroscopy of Venus: Detection of OCS, upper limit to H2S, and latitudinal variations of CO and HF in the upper cloud layer. Icarus 2008, 197, 377–385. [Google Scholar] [CrossRef]
  59. Yin, Y.; Zhai, K.; Zhang, B.; Zhai, S. Electrical resistivity of iron phosphides at high-pressure and high-temperature conditions with implications for lunar core’s thermal conductivity. J. Geophys. Res. Solid Earth 2019, 124, 5544–5556. [Google Scholar] [CrossRef]
  60. Pearson, D.G.; Scott, J.M.; Liu, J.; Schaeffer, A.; Wang, L.H.; van Hunen, J.; Szilas, K.; Chacko, T.; Kelemen, P.B. Deep continental roots and cratons. Nature 2021, 596, 199–210. [Google Scholar] [CrossRef]
  61. Leung, I.; Guo, W.; Friedman, I.; Gleason, J. Natural occurrence of silicon carbide in a diamondiferous kimberlite from Fuxian. Nature 1990, 346, 352–354. [Google Scholar] [CrossRef]
  62. Robinson, P.T.; Bai, W.-J.; Malpas, J.; Yang, J.-S.; Zhou, M.-F.; Fang, Q.-S.; Hu, X.-F.; Cameron, S.; Staudigel, H. Ultra-high pressure minerals in the Luobusa Ophiolite, Tibet, and their tectonic implications. Geol. Soc. Lond. Spec. Publ. 2004, 226, 247–271. [Google Scholar] [CrossRef]
  63. Pierro, S.D.; Gnos, E.; Grobéty, B.; Armbruster, T.M.; Bernasconi, S.; Ulmer, P. Rock-forming moissanite (natural ∂-silicon carbide). Am. Mineral. 2003, 88, 1817–1821. [Google Scholar] [CrossRef]
  64. Yang, J.; Meng, F.; Xu, X.; Robinson, P.T.; Dilek, Y.; Makeyev, A.B.; Wirth, R.; Wiedenbeck, M.; Cliff, J. Diamonds, native elements and metal alloys from chromitites of the Ray-Iz ophiolite of the Polar Urals. Gondwana Res. 2015, 27, 459–485. [Google Scholar] [CrossRef]
  65. Shiryaev, A.A.; Griffin, W.L.; Stoyanov, E. Moissanite (SiC) from kimberlites: Polytypes, trace elements, inclusions and speculations on origin. Lithos 2011, 122, 152–164. [Google Scholar] [CrossRef]
  66. Xu, S.; Wu, W.; Xiao, W.; Yang, J.; Chen, J.; Ji, S.; Yican Liu, U. Moissanite in serpentinite from the Dabie Mountains in China. Mineral. Mag. 2018, 72, 899–908. [Google Scholar] [CrossRef]
  67. Airey, M.W.; Mather, T.; Pyle, D.; Glaze, L.; Ghail, R.; Wilson, C. Explosive volcanic activity on Venus: The roles of volatile contribution, degassing, and external environment. Planet. Space Sci. 2015, 113, 33–48. [Google Scholar] [CrossRef] [Green Version]
  68. Dera, P.; Lavina, B.; Borkowski, L.A.; Prakapenka, V.B.; Sutton, S.R.; Rivers, M.L.; Downs, R.T.; Boctor, N.Z.; Prewitt, C.T. High-pressure polymorphism of Fe2P and its implications for meteorites and Earth’s core. Geophys. Res. Lett. 2008, 35, L10301. [Google Scholar] [CrossRef] [Green Version]
  69. Britvin, S.; Murashko, M.; Vapnik, E.; Polekhovsky, Y.S.; Krivovichev, S. Barringerite Fe2P from pyrometamorphic rocks of the Hatrurim Formation, Israel. Geol. Ore Depos. 2017, 59, 619–625. [Google Scholar] [CrossRef]
  70. Krasnopolsky, V.A. Chemical kinetic model for the lower atmosphere of Venus. Icarus 2007, 191, 25–37. [Google Scholar] [CrossRef]
  71. Glaze, L.S. Transport of SO2 by explosive volcanism on Venus. J. Geophys. Res. Planets 1999, 104, 18899–18906. [Google Scholar] [CrossRef]
  72. Esposito, L.W. Sulfur dioxide: Episodic injection shows evidence for active Venus volcanism. Science 1984, 223, 1072–1074. [Google Scholar] [CrossRef]
  73. La Femina, P.C.; Connor, C.B.; Hill, B.E.; Strauch, W.; Saballos, J.A. Magma–tectonic interactions in Nicaragua: The 1999 seismic swarm and eruption of Cerro Negro volcano. J. Volcanol. Geotherm. Res. 2004, 137, 187–199. [Google Scholar] [CrossRef]
  74. Carr, M.; Feigenson, M.; Patino, L.; Walker, J. Volcanism and geochemistry in Central America: Progress and problems. Geophys. Monogr.-Am. Geophys. Union 2003, 138, 153–174. [Google Scholar]
  75. Cole, J.W. Structure, petrology, and genesis of Cenozoic volcanism, Taupo Volcanic Zone, New Zealand—A review. N. Zeal. J. Geol. Geophys. 1979, 22, 631–657. [Google Scholar] [CrossRef] [Green Version]
  76. Hogg, A.G.; Higham, T.F.; Lowe, D.J.; Palmer, J.G.; Reimer, P.J.; Newnham, R.M. A wiggle-match date for Polynesian settlement of New Zealand. Antiquity 2003, 77, 116–125. [Google Scholar] [CrossRef] [Green Version]
  77. Sahetapy-Engel, S.; Self, S.; Carey, R.J.; Nairn, I.A. Deposition and generation of multiple widespread fall units from the c. AD 1314 Kaharoa rhyolitic eruption, Tarawera, New Zealand. Bull. Volcanol. 2014, 76, 836. [Google Scholar] [CrossRef]
  78. Colgate, S.A.; Sigurgeisson, T. Dynamic mixing of water and lava. Nature 1973, 244, 552–555. [Google Scholar] [CrossRef]
  79. Ort, M.H.; Di Muro, A.; Michon, L.; Bachèlery, P. Explosive eruptions from the interaction of magmatic and hydrothermal systems during flank extension: The Bellecombe Tephra of Piton de La Fournaise (La Réunion Island). Bull. Volcanol. 2016, 78, 5. [Google Scholar] [CrossRef] [Green Version]
  80. Antonopoulos, J. The great Minoan eruption of Thera volcano and the ensuing tsunami in the Greek Archipelago. Nat. Hazards 1992, 5, 153–168. [Google Scholar] [CrossRef]
  81. Bachmann, O.; Bergantz, G.W. Rhyolites and their Source Mushes across Tectonic Settings. J. Petrol. 2009, 49, 2277–2285. [Google Scholar] [CrossRef] [Green Version]
  82. Tian, F.; Chassefière, E.; Leblanc, F.; Brain, D.A. Atmosphere escape and climate evolution of terrestrial planets. In Comparative Climatology of Terrestrial Planets; University of Arizona Press: Tucson, AZ, USA, 2013; pp. 567–581. [Google Scholar]
  83. Gaillard, F.; Scaillet, B. A theoretical framework for volcanic degassing chemistry in a comparative planetology perspective and implications for planetary atmospheres. Earth Planet. Sci. Lett. 2014, 403, 307–316. [Google Scholar] [CrossRef] [Green Version]
  84. Kelley, K.A.; Plank, T.; Newman, S.; Stolper, E.M.; Grove, T.L.; Parman, S.; Hauri, E.H. Mantle melting as a function of water content beneath the Mariana Arc. J. Petrol. 2010, 51, 1711–1738. [Google Scholar] [CrossRef] [Green Version]
  85. Eggler, D.H. Water-saturated and undersaturated melting relations in a Paricutin andesite and an estimate of water content in the natural magma. Contrib. Mineral. Petrol. 1972, 34, 261–271. [Google Scholar] [CrossRef]
  86. Baker, M.B.; Grove, T.L.; Price, R. Primitive basalts and andesites from the Mt. Shasta region, N. California: Products of varying melt fraction and water content. Contrib. Mineral. Petrol. 1994, 118, 111–129. [Google Scholar] [CrossRef]
  87. Jamtveit, B.; Brooker, R.; Brooks, K.; Larsen, L.M.; Pedersen, T. The water content of olivines from the North Atlantic Volcanic Province. Earth Planet. Sci. Lett. 2001, 186, 401–415. [Google Scholar] [CrossRef]
  88. Moore, J.G. Water Content of Basalt Erupted on the ocean floor. Contrib. Mineral. Petrol. 1970, 28, 272–279. [Google Scholar] [CrossRef]
  89. Geng, J.; Zhang, R.; Wang, X. Chemical Origin of Phosphine in Nature. In Proceedings of the 2010 4th International Conference on Bioinformatics and Biomedical Engineering, Chengdu, China, 18–20 June 2010; pp. 1–12. [Google Scholar]
  90. Cardelino, B.; Moore, C.; Cardelino, C.; McCall, S.; Frazier, D.; Bachmann, K. Semiclassical calculation of reaction rate constants for homolytical dissociation reactions of interest in organometallic vapor-phase epitaxy (OMVPE). J. Phys. Chem. A 2003, 107, 3708–3718. [Google Scholar] [CrossRef]
  91. Allison, C.M.; Roggensack, K.; Clarke, A.B. Highly explosive basaltic eruptions driven by CO2 exsolution. Nat. Commun. 2021, 12, 217. [Google Scholar] [CrossRef]
  92. Krueger, A.J. Sighting of El Chichon sulfur dioxide clouds with the Nimbus 7 total ozone mapping spectrometer. Science 1983, 220, 1377–1379. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  93. Byrne, P.K. A comparison of inner Solar System volcanism. Nat. Astron. 2020, 4, 321–327. [Google Scholar] [CrossRef]
  94. Head, J.W.; Crumpler, L.; Aubele, J.C.; Guest, J.E.; Saunders, R.S. Venus volcanism: Classification of volcanic features and structures, associations, and global distribution from Magellan data. J. Geophys. Res. Planets 1992, 97, 13153–13197. [Google Scholar] [CrossRef] [Green Version]
  95. Head, J.W., III; Wilson, L. Volcanic processes and landforms on Venus: Theory, predictions, and observations. J. Geophys. Res. Solid Earth 1986, 91, 9407–9446. [Google Scholar] [CrossRef]
  96. Aramaki, S.; Hayakawa, Y.; Fujii, T.; Nakamura, K.; Fukuoka, T. The October 1983 eruption of Miyakejima volcano. J. Volcanol. Geotherm. Res. 1986, 29, 203–229. [Google Scholar] [CrossRef]
  97. Parfitt, E.A.; Wilson, L. A Plinian treatment of fallout from Hawaiian lava fountains. J. Volcanol. Geotherm. Res. 1999, 88, 67–75. [Google Scholar] [CrossRef]
  98. Mastin, L.G.; Guffanti, M.; Servranckx, R.; Webley, P.; Barsotti, S.; Dean, K.; Durant, A.; Ewert, J.W.; Neri, A.; Rose, W.I.; et al. A multidisciplinary effort to assign realistic source parameters to models of volcanic ash-cloud transport and dispersion during eruptions. J. Volcanol. Geotherm. Res. 2009, 186, 10–21. [Google Scholar] [CrossRef]
  99. Earthchem Collaboration. Earthchem. 2021, Volume 2021. Available online: https://www.earthchem.org/ (accessed on 8 August 2021).
  100. Sen, G. Generation of Deccan trap magmas. J. Earth Syst. Sci. 2001, 110, 409–431. [Google Scholar] [CrossRef] [Green Version]
  101. Renne, P.R.; Basu, A.R. Rapid eruption of the Siberian traps flood basalts at the permo-triassic boundary. Science 1991, 253, 176–179. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  102. Blasius, K.R.; Cutts, J.A. Topography of Martian central volcanoes. Icarus 1981, 45, 87–112. [Google Scholar] [CrossRef]
  103. Hardee, H.C.; Larson, D.W. The extraction of heat from magmas based on heat transfer mechanisms. J. Volcanol. Geotherm. Res. 1977, 2, 113–144. [Google Scholar] [CrossRef]
  104. Sachs, P.M.; Stange, S. Fast assimilation of xenoliths in magmas. J. Geophys. Res. Solid Earth 1993, 98, 19741–19754. [Google Scholar] [CrossRef]
  105. Robinson, J.E.; Eakins, B.W. Calculated volumes of individual shield volcanoes at the young end of the Hawaiian Ridge. J. Volcanol. Geotherm. Res. 2006, 151, 309–317. [Google Scholar] [CrossRef]
Figure 1. Solid pieces of iron can only be transported substantially unaltered to the surface from the deep mantle by a mantle plume if they are >1 cm radius. X-axis: depth (km). Y-axis: minimum radius for a sphere of iron to be transported from that depth to the surface, while losing no more than 50% of their phosphorus to diffusion. Plots show the temperature profile of the lithosphere, or that of a plume 300 °C hotter than the surrounding mantle [33], pooling below a 30 km [39] or a 150 km [40] thick lithosphere. Red lines show the required radius for the internal phosphorus concentration to be reduced by 50%; blue lines show the radius for the internal oxygen concentration to reach 50% of the external concentration.
Figure 1. Solid pieces of iron can only be transported substantially unaltered to the surface from the deep mantle by a mantle plume if they are >1 cm radius. X-axis: depth (km). Y-axis: minimum radius for a sphere of iron to be transported from that depth to the surface, while losing no more than 50% of their phosphorus to diffusion. Plots show the temperature profile of the lithosphere, or that of a plume 300 °C hotter than the surrounding mantle [33], pooling below a 30 km [39] or a 150 km [40] thick lithosphere. Red lines show the required radius for the internal phosphorus concentration to be reduced by 50%; blue lines show the radius for the internal oxygen concentration to reach 50% of the external concentration.
Universe 08 00054 g001
Figure 2. Phosphorus partitions into silicate above 1000 km depth. The partition coefficient of phosphorus Dp between metallic and silicate phases predicted for the mantle shows that below 600 km (for a mantle at FMQ reduction) or 1000 km (for a QIF mantle) is the majority of phosphorus in the metallic phase as iron phosphide. X-axis—depth below the surface (km). Left-hand Y-axis: Partition coefficient Dp of phosphorus between metallic and silicate phase expected at that depth. Right Y-axis: pressure (green) and temperature (red) at that depth.
Figure 2. Phosphorus partitions into silicate above 1000 km depth. The partition coefficient of phosphorus Dp between metallic and silicate phases predicted for the mantle shows that below 600 km (for a mantle at FMQ reduction) or 1000 km (for a QIF mantle) is the majority of phosphorus in the metallic phase as iron phosphide. X-axis—depth below the surface (km). Left-hand Y-axis: Partition coefficient Dp of phosphorus between metallic and silicate phase expected at that depth. Right Y-axis: pressure (green) and temperature (red) at that depth.
Universe 08 00054 g002
Figure 3. Thermodynamic analysis shows that phosphide is strongly disfavored over phosphate at lithospheric depths. X-axes: depth in kilometers. Y-axes: oxygen fugacity (log10[f(O2)]), calculated from free energy of reaction forming metal phosphide and oxygen from metal phosphate. (A) Full analysis to depth of 2000 km, for melt temperature of 1100 K. (B) Detail on the top 200 km, for melt temperature of 1100 K. (C) Top 200 km for melt temperature of 660 K. Iron phosphides are shown to be stable over phosphates at depths of 850 km or 1600 km if the mantle is QIF/IM or FMQ redox state respectively, but no phosphide is preferred over phosphate above 800 km in depth.
Figure 3. Thermodynamic analysis shows that phosphide is strongly disfavored over phosphate at lithospheric depths. X-axes: depth in kilometers. Y-axes: oxygen fugacity (log10[f(O2)]), calculated from free energy of reaction forming metal phosphide and oxygen from metal phosphate. (A) Full analysis to depth of 2000 km, for melt temperature of 1100 K. (B) Detail on the top 200 km, for melt temperature of 1100 K. (C) Top 200 km for melt temperature of 660 K. Iron phosphides are shown to be stable over phosphates at depths of 850 km or 1600 km if the mantle is QIF/IM or FMQ redox state respectively, but no phosphide is preferred over phosphate above 800 km in depth.
Universe 08 00054 g003
Table 1. Constants for Equation (4) to calculate diffusion coefficient.
Table 1. Constants for Equation (4) to calculate diffusion coefficient.
ElementSolidA (cm2/Day)Ea (kJ/mol)NotesReference
PIron9.84·104306Calculated from data in paper[29]
OIron2.396·10−1118Minimally affected by dissolved Si, Al or Mg[30]
Table 2. Summary of factors required for Truong and Lunine [16] base model.
Table 2. Summary of factors required for Truong and Lunine [16] base model.
FactorCommentSectionMultiplier
Phosphine lifetime t P H 3 could be as specified by Truong and Lunine or full photochemical model3.61–26
Abundance of erupted phosphide, assuming QIF f(O2)<10−5 of phosphorus likely to be phosphide within 100 km of surface3.1, 3.2>6000
Delivery of phosphide to cloudsEffusive volcanism at least 10-fold less efficient at producing fine ash as explosive volcanism3.3, 3.4101–102
Abundance of phosphorusMost likely 0.08%, not 1%3.5, [11]12
Overall scale-up of [16] estimate required to meet PH3 production requirement ≥7.2·105
Volume of volcanic eruption required, based on Truong and Lunine base case of 0.03–0.15 km3/year≥26,100 km3/year
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bains, W.; Shorttle, O.; Ranjan, S.; Rimmer, P.B.; Petkowski, J.J.; Greaves, J.S.; Seager, S. Constraints on the Production of Phosphine by Venusian Volcanoes. Universe 2022, 8, 54. https://doi.org/10.3390/universe8010054

AMA Style

Bains W, Shorttle O, Ranjan S, Rimmer PB, Petkowski JJ, Greaves JS, Seager S. Constraints on the Production of Phosphine by Venusian Volcanoes. Universe. 2022; 8(1):54. https://doi.org/10.3390/universe8010054

Chicago/Turabian Style

Bains, William, Oliver Shorttle, Sukrit Ranjan, Paul B. Rimmer, Janusz J. Petkowski, Jane S. Greaves, and Sara Seager. 2022. "Constraints on the Production of Phosphine by Venusian Volcanoes" Universe 8, no. 1: 54. https://doi.org/10.3390/universe8010054

APA Style

Bains, W., Shorttle, O., Ranjan, S., Rimmer, P. B., Petkowski, J. J., Greaves, J. S., & Seager, S. (2022). Constraints on the Production of Phosphine by Venusian Volcanoes. Universe, 8(1), 54. https://doi.org/10.3390/universe8010054

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop