John M. C.
Plane
*
School of Chemistry, University of Leeds, Woodhouse Lane, Leeds LS2 9JT, U.K. E-mail: j.m.c.plane@leeds.ac.uk; Fax: 44 113 343 6401; Tel: 44 113 343 8044
First published on 7th June 2012
This review discusses the magnitude of the cosmic dust input into the earth's atmosphere, and the resulting impacts from around 100 km to the earth's surface. Zodiacal cloud observations and measurements made with a spaceborne dust detector indicate a daily mass input of interplanetary dust particles ranging from 100 to 300 tonnes, which is in agreement with the accumulation rates of cosmic-enriched elements (Ir, Pt, Os and super-paramagnetic Fe) in polar ice cores and deep-sea sediments. In contrast, measurements in the middle atmosphere – by radar, lidar, high-flying aircraft and satellite remote sensing – indicate that the input is between 5 and 50 tonnes per day. There are two reasons why this huge discrepancy matters. First, if the upper range of estimates is correct, then vertical transport in the middle atmosphere must be considerably faster than generally believed; whereas if the lower range is correct, then our understanding of dust evolution in the solar system, and transport from the middle atmosphere to the surface, will need substantial revision. Second, cosmic dust particles enter the atmosphere at high speeds and undergo significant ablation. The resulting metals injected into the atmosphere are involved in a diverse range of phenomena, including: the formation of layers of metal atoms and ions; the nucleation of noctilucent clouds, which are a sensitive marker of climate change; impacts on stratospheric aerosols and O3 chemistry, which need to be considered against the background of a cooling stratosphere and geo-engineering plans to increase sulphate aerosol; and fertilization of the ocean with bio-available Fe, which has potential climate feedbacks.
![]() John M. C. Plane | Following his PhD in physical chemistry at Cambridge University (1983), John Plane was a Research Fellow at St. John's College and then Associate Professor at the Rosenstiel School of Marine and Atmospheric Science (University of Miami, Florida). He moved to the University of East Anglia (1991), becoming Professor of Environmental Science in 1999. Since 2006 he has been Professor of Atmospheric Chemistry at the University of Leeds. |
This review addresses an apparently simple question: what is the magnitude of the cosmic dust input to the earth's atmosphere? Table 1 shows that even very recent estimates of the Interplanetary Dust Particle (IDP) input vary from 5 to 270 t d−1 (tonnes per day). Zodiacal cloud observations and spaceborne dust detection (dark blue shading in Table 1) indicate a daily input of 100–300 t d−1, which is mostly in agreement with the accumulation rates of cosmic elements in polar ice cores and deep-sea sediments (grey shading). In contrast, measurements in the middle atmosphere (light blue shading) – by radar, lidar, high-flying aircraft and satellite remote sensing – indicate that the input is only 5–50 tonnes.
![]() |
There are two reasons why this matters. First, if the upper range of estimates is correct, then vertical transport in the middle atmosphere must be considerably faster than is generally thought to be the case, so that meteoritic material is removed more rapidly from the atmosphere in order to sustain a higher rate of injection; whereas if the lower range is correct, then our understanding of dust evolution in the solar system, and transport mechanisms from the middle atmosphere to the earth's surface, will need substantial revision. Second, cosmic dust particles enter the atmosphere at high speeds (11–72 km s−1) and in most cases completely ablate.5 The resulting metals injected into the atmosphere are involved in a diverse range of impacts, including the formation of layers of metal atoms and ions,6 nucleation of noctilucent clouds,7 effects on stratospheric aerosols and O3 chemistry, and fertilization of the ocean with bio-available Fe.8 These impacts of meteoric ablation obviously depend on the magnitude of the IDP input.
Although the actual measurements in Table 1 appear to be sound, their interpretation to yield the estimated IDP input is potentially compromised by significant uncertainties (final column of Table 1) which are discussed in detail below. The enormous range of values in Table 1 implies a fundamental lack of understanding in at least some parts of the system. In fact, the range is so large that even if the true input is in the “middle” (∼20–50 t d−1), significant reassessment of many of these processes will be required.
In this review, the different ways of estimating the IDP input and assessing its impact are grouped into five Science Topics (STs), illustrated in Fig. 1. Moving downwards from the top of the atmosphere, these topics are: the source of IDPs in the inner solar system and the process of meteoric ablation (ST1); the layers of metal atoms and ions which result from ablation (ST2); the formation of meteoric smoke particles and their role as ice nuclei in the mesosphere (ST3); the impacts of meteoric smoke on aerosols and O3 chemistry in the stratosphere (ST4); and the deposition of smoke to the earth's surface (ST5). Fig. 1 illustrates the relationship between these STs. Each ST (grey shaded box) contains phenomena in the Earth System (red boxes) which provide an estimate of the IDP input through relevant observations (blue boxes). Note that a list of acronyms appears at the end of the review.
![]() | ||
Fig. 1 Structure of the review material, grouped into five Science Topics. Each topic (grey shaded box) contains phenomena in the Earth System (red boxes) which provide an estimate of the IDP input through relevant observations (blue boxes). |
The input flux of meteoroids into the atmosphere is so uncertain because no single technique can observe particles over the mass range from about 10−12 to 1 g which make up the bulk of the incoming material.3Fig. 2 shows that the particle mass can vary by 30 orders of magnitude, although the largest contribution of mass entering the atmosphere on a daily basis comes from particles around 10 μg. Assuming a meteoroid density of ∼2.8 g cm−3, these particles will have a diameter of ∼200 μm. There is a population of huge impactors with masses greater than 1010 g which make a significant contribution, but only on a geological timescale! Any single measurement technique will only sample a subset of this size distribution. For instance, optical camera networks which observe visible meteors detect particles larger than about 1 mg in mass, or 1 mm in radius. Larger particles (approaching 1 g in mass) are much rarer, so that counting statistics on a time scale of months start to matter.
![]() | ||
Fig. 2 Mass influx (per decade of mass) plotted against particle mass [data taken from Flynn21]. |
Meteor radars measure particles with masses between about 10−9 and 10−3 g, and therefore cover the most important mass range (Fig. 2). Meteor radar data was used to produce a much-quoted estimate of 44 t d−1 for the global input, although this involved artificially increasing the size distribution to match visual meteor observations.11 The evaporating atoms, particularly metals, ionize through hyperthermal collisions with air molecules.5 This creates a trail of electrons behind the meteoroid, which can be detected by radar. The mass and speed of the meteoroid then have to be estimated indirectly.22 Furthermore, the wavelength of the radar only samples a subset of the mass/velocity/altitude distribution of the meteoroids, so that some extrapolation is required to estimate the total mass input.11
In the past two decades, high-powered large aperture (HPLA) radars, such as the Arecibo Observatory and the EISCAT radars in the Arctic, have been able to detect by incoherent scatter the meteor head echo (i.e. the ball of plasma around the ablating particle as it descends through the atmosphere). This enables measurements of the direction of origin, velocity, deceleration and (indirectly) mass to be made.10,23–25 While initially the mean entry velocity seemed to be significantly higher, around 40–50 km s−1, than the velocity measured with conventional meteor radars, it has now been realised that there is a sampling bias towards high-speed meteors.26 Conventional meteor radars do not efficiently detect meteors which occur at higher altitudes (>100 km), because of the rapid diffusion of the ionized trails. Since faster meteors generally occur at higher altitudes, distributions measured by meteor radars are biased towards the lower speeds. In fact, it has now been shown27 that HPLA radars observe the same population of meteors as observed by meteor radars, and in addition detect a population of faster meteors that ablate at altitudes where trails are not efficiently detected. However, the magnitude of the head echo still depends on the meteoroid mass and velocity, and each HPLA radar is sensitive to a particular mass range.27 This implies that the velocity distribution of the smallest particles measured by an HPLA radar will be biased towards faster speeds: small and slow particles will not have sufficient kinetic energy to ablate, and hence will not produce sufficient electrons to be detected. The average entry speeds are now thought to be between 25 and 30 km s−1.28 It should be noted that particles which originate within the solar system must have entry velocities that range from 11.5 km s−1 for a particle in the same prograde orbit as the Earth (i.e. orbiting in the same direction), to 72 km s−1 for a particle in a retrograde orbit.2
The population of IDPs smaller than 10−9 g can only be measured by impact detectors on satellites. An important estimate of the IDP input was provided by the Long Duration Exposure Facility (LDEF), an orbital impact detector placed on a spacecraft for several years, which yielded an estimate of 110 t d−1.9 However, the LDEF experiment measured crater size, which was treated as a proxy for particle kinetic energy. Hence, the particle velocity distribution had to be assumed in order to determine the mass distribution. If the average velocity is higher (see above) than the value of only 18 km s−1 that was employed in the LDEF analysis, then the corresponding mass distribution would be shifted down by more than an order of magnitude.10
Because of their very high entry velocities, meteoroids undergo rapid frictional heating by collision with air molecules. If the particles reach melting point (∼1800 K), their constituent minerals will then rapidly vaporize – the process termed meteoric ablation. Ablation tends to occur where the atmospheric pressure is around 1 μbar. In the case of the Earth, the peak ablation rate is around 90 km, compared with 80 km on Mars, 115 km on Venus, and 500 km on Titan.29 The physical chemistry of ablation has been treated in detail by several investigators.30–32 The problem becomes manageable for particles smaller than about 250 μm in radius, because heat conductivity through the particle is then fast enough for the particle to be treated as isothermal.30 Assuming a density of 2.8 g cm−3 (typical of an ordinary chondrite), a 250 μm meteoroid has a mass of 180 μg (which, as Fig. 2 shows, is larger than the bulk of the IDP mass input).
The frictional heating of the meteoroid by collisions with air molecules is balanced by radiative loss and the consumption of heat energy through temperature increase, melting and vaporization.5 In order to calculate these terms, parameters such as the meteoroid shape, density, and composition are needed. The question of composition has been discussed in detail recently.33,34 There is some uncertainty here because of the great variability in composition of different types of meteorites.33 Furthermore, it may be that the composition of the meteoroids that ablate in the upper atmosphere is different from that of the meteorites that have survived transit through the atmosphere.33 Nevertheless, the current assumption is that most IDPs have the composition of ordinary chondrites, which is essentially olivine (FeMgSiO4). The elemental abundances of the major metallic constituents relative to Si (1.0) are then: Mg (1.07), Fe (0.90), Al (0.085), Ca (0.061), Na (0.057) and Ni (0.049).35
If the particle is roughly spherical (which should be the case once it has melted), then the energy balance equation is given by:36
![]() | (I) |
The deceleration of the meteoroid is given by
![]() | (II) |
Inspection of eqn (I) shows that for very small particles the heat capacity term will be much smaller than the radiative loss term, in which case the meteoroid will not become hot enough to ablate. This is the case for all particles smaller than 10−12 g. A 10−10 g particle must enter the atmosphere at over 40 km s−1 in order for its temperature to exceed 1800 K, the temperature at which volatile elements such as Na begin to evaporate rapidly. All meteoroids larger than 10−7 g should reach this ablation temperature.5
Interpreting radar measurements requires a model to predict the evaporation rate and subsequent ionization rate of the elemental constituents during ablation. The most recent example is the Chemical Ablation Model (CABMOD),5 which treats the physics and chemistry of ablation by including sputtering by inelastic collisions with air molecules before the meteoroid melts, evaporation of atoms and oxides from the molten particle, and impact ionization of the ablated fragments by hyperthermal collisions with air molecules. Evaporation is based on the assumption of thermodynamic equilibrium within the molten meteoroid and between the particle and the surrounding vapour phase. The thermodynamics are provided by the well-established MAGMA code,37 which is based on a fractionation model used by planetary scientists to account for mass loss from the planet Mercury. The loss rate of each element is then calculated by applying Langmuir evaporation, which assumes that the rate of evaporation into a vacuum is equal to the rate of evaporation needed to balance the rate of uptake of a species i in a closed system. The rate of mass release of species i with molecular weight μi is given by the Herz–Knudsen equation:
![]() | (III) |
Fig. 3 illustrates the elemental injection profiles calculated by CABMOD for a meteoroid entering the atmosphere with the most likely mass and velocity according to LDEF.5 The mass loss which occurs above 105 km, when the particle has not yet melted (Tm < 1800 K), is due to sputtering. Note that CABMOD predicts differential ablation, i.e. the most volatile elements – Na and K – ablate first, followed by the main constituents Fe, Mg and Si, and finally the most refractory elements such as Ca. CABMOD has been used successfully to model the variation of meteor head echoes with height, using measurements made at the Arecibo Observatory.38 However, there are two important caveats. First, much of the thermodynamic data in MAGMA has to be extrapolated to temperatures above 2000 K from measurements below 1700 K.37 Second, Langmuir evaporation represents an upper limit to the evaporation rate, since evaporation into a vacuum can be significantly slower than into a vapour at equilibrium (because diffusion from the bulk into the surface film can become rate-determining). For instance, the evaporation coefficient from a molten silicate into vacuum can be much less than 1 for elements such as Fe (0.25) and Ca (0.06),39 which affects the height in the atmosphere where these elements ablate (and also the heating rate of the meteoroid).
![]() | ||
Fig. 3 Elemental ablation profiles for a 5 μg meteoroid entering at 20 km s−1, as predicted by the Chemical Ablation Model (CABMOD). The particle temperature is shown on the top abscissa. |
Further work is need to remove these uncertainties in chemical ablation models, so that they can be used to correct for biases in the meteor mass/velocity distribution measured by radars.26 One important reason for doing this is to refine the ZCM described above, by constraining the predicted IDP orbits and velocities by using radar data. An initial attempt in this direction was published very recently, showing that the IDP input could be reduced to only 41 t d−1 by changing the initial orbital characteristics of IDPs ejected by Jupiter family comets.40
Lidar observations can be made with extremely good time and height resolution (typically 60 s and 40 m, respectively, for the Na layer), so that the metal layers can be used as tracers of atmospheric motions such as tides and gravity waves.43 In the case of Na, K and Fe, a narrow line-width laser can be used in the lidar transmitter to measure temperature and wind profiles in the upper mesosphere.42,44,45 This is achieved by scanning the laser across the resonance line to measure the degree of Doppler broadening, and hence the local temperature. By employing a laser with a tuning accuracy and frequency stability of at least 50 MHz, the temperature can be measured with an error of ∼2 K. The wind along the line-of-sight of the lidar can be determined from the net Doppler shift of the resonance line, typically with an error of less than 3 m s−1. By pointing the lidar sequentially in both zenith and off-zenith directions, the wind can then be resolved into the zonal, merional and even vertical components. Another way to measure the temperature is to use a two-colour lidar to measure simultaneously the relative populations of the 5D4 and 5D3 spin–orbit multiplets of ground-state Fe, from which the temperature can be derived assuming Boltzmann equilibrium.46
Metal resonance lidars are thus an extremely important tool for studying the chemistry and physics of the mesosphere/lower thermosphere (MLT), which is largely inaccessible to direct measurements. High altitude aircraft and research balloons reach altitudes of about 22 and 45 km, respectively, whereas satellites cannot operate below 150 km without atmospheric drag causing rapid re-entry. The only way to sample the MLT in situ is via rocket-borne payloads. However, because the payload traverses the MLT region at a velocity typically in excess of 1 km s−1, the instruments must have a rapid time response.47 There have been a number of measurements of the concentrations of positive metallic ions made by rocket-borne mass spectrometry.48,49 These flights have been motivated to establish whether there is a link between meteor showers and the abundance of metallic ions (there does not appear to be a significant correlation49), and the role of metallic ions in forming sporadic E layers and noctilucent clouds48 (see ST3).
Metallic ions such as Mg+ and Fe+ have also been observed by resonant scattering of sunlight, using spectrometers on space vehicles.50 More recently, satellite observations have also been made of the neutral Na and Mg layers.51–53 The vertical resolution achievable with limb-scanning spectrometers on satellites is only about 2 km (cf. the lidar vertical resolution of 40 m), but satellites in polar sun-synchronous orbits provide near-global coverage. Fig. 4 illustrates the Na column abundance (i.e. the concentration of the Na layer integrated over height), as a function of latitude and season. The satellite data-set51 has been supplemented with ground-based lidar data during the polar winter when the Na layer is not solar-illuminated. Fig. 4 illustrates that there is very little seasonal variation at low latitudes (less than a factor of 2), but the winter/summer ratio increases to nearly an order of magnitude at high latitudes. The primary reason for the large summertime depletion at high latitudes is the very low temperature during the summer because of the adiabatic cooling of upwelling air.6 A secondary reason is the efficient removal of metallic species on noctilucent cloud particles (see ST3).54,55
![]() | ||
Fig. 4 Column abundance of the Na layer (units: 109 atom cm−2) as a function of latitude and month. |
Understanding the characteristic features of the metallic layers in the upper atmosphere has required studying the reaction kinetics of neutral and ionized metallic species with atmospheric constituents such as O3, O2, O and H. Over the past 30 years, the two classical techniques of flash photolysis and the fast flow tube have provided a great deal of kinetic data on the pertinent reactions of the meteoric metals Fe, Mg, Na, Ca and K.6,41,56 The results have been used to construct atmospheric models which successfully explain the metal layers above 80 km.6,12,13,41,57Fig. 5 illustrates the author's current model of iron chemistry in the upper atmosphere. The rate coefficients for the reactions depicted with black arrows have all been measured. Above 95 km, Fe+ ions are produced by charge transfer with the lower E region ions NO+ and O2+. Neutralisation of Fe+ then occurs via reaction with O3 to form FeO+, or recombination with O2 to form FeO2+, producing molecular ions which can then undergo dissociative recombination (DR) with electrons to generate Fe (red arrows in Fig. 4).
![]() | ||
Fig. 5 Fe chemistry in the MLT. Black arrows: reactions with measured rate coefficients; red, blue, brown arrows indicate reactions which need laboratory study (see text). |
Metallic ions, particularly Fe+ and Mg+, are the major constituents of sporadic E layers, thin layers of concentrated plasma which occur in the lower thermosphere (95–130 km).48,49 Sporadic E layers are important for radio communications, both enabling over-the-horizon radio propagation and attenuating ground-to-space communications. The chemical lifetimes (against neutralisation) of these metallic ions are controlled by the ion-molecule chemistry illustrated for the case of Fe+ in Fig. 5. This shows that once molecular ions such as FeO+ or FeO2+ form, there is a competition between DR and conversion back to Fe+ by reactions with atomic O, which effectively slows down neutralisation. When a sporadic E layer descends below 100 km, the concentration of atomic O decreases markedly in relation to O3,O2 and N2.6 Thus, molecular ions form more quickly and are converted back to Fe+ more slowly. The combined effect is to reduce the lifetime of Fe+ from days above 100 km, to only minutes at 90 km.58
Recently, low-lying sporadic ion layers have been observed, using radio occultation from spacecraft, to occur around 90 km on Mars,59,60 120 km on Venus61 and 550 km on Titan.62 The CO2 atmospheres of Mars and Venus pose a particular challenge to the existence of metallic ions, because they should form CO2-clusters very rapidly and undergo DR. In fact, it turns out that sporadic layers in the Martian atmosphere are most likely Mg+ rather than Fe+, following a recent laboratory study63 which showed that atomic O reacts much more rapidly with molecular Mg-containing ions to form Mg+, compared with their Fe-containing analogues. This very effectively slows down the neutralisation of Mg+ in a CO2-rich atmosphere. Titan is interesting because ablation occurs over a much greater altitude range, as a result of its small atmospheric scale height.29 An important reason for carrying out comparative studies with these other atmospheres is that this provides a self-consistency check of the meteor input functions produced from the ZCM throughout the solar system.
Below 85 km in the terrestrial atmosphere, Fig. 5 shows that atomic Fe is oxidised in a series of reactions involving O3, O2, CO2 and H2O to form FeO3, Fe(OH)2 and FeOH.64 The latter may also be oxidized by O3 to form FeOOH (which is the building block of the mineral goethite). The reaction kinetics (blue arrows in Fig. 5) and photochemistry of species such as FeOH, FeOOH, FeO3 with H need to be studied experimentally, because they are key to understanding how quickly the metal reservoir species are permanently removed. This is thought to happen through the polymerisation of these compounds together with SiO2 vapour in the mesosphere over several days (brown arrows in Fig. 5), forming nanometre-sized meteoric smoke particles (MSPs) which provide a permanent sink for gas-phase metallic compounds (see ST3).65
The chemistry in Fig. 5 was recently incorporated into a 1-D model of the MLT to study the behavior of the Fe layer over Rothera (Antarctica), where a Fe lidar had been operated for two years.13Fig. 6 (upper panel) illustrates the observed seasonal variation of the layer height profile. As in the case of Na layer (Fig. 4), the minimum occurs during the summer months when the MLT is cold and noctilucent clouds are present. The maximum in the Fe layer occurs in late autumn rather than mid-winter, because the MLT temperature and the meteor input function both maximize then.13 The lower panel of Fig. 6 shows that the model is able to capture this seasonal variation, as well as the correct height of the peak of the layer.
![]() | ||
Fig. 6 Seasonal variation of the Fe layer vertical profile at Rothera (Antarctica). Top panel: lidar measurements. Bottom panel: model simulation with an Fe injection rate equivalent to an IDP input of 6 t d−1.13 |
Strikingly, an IDP input of only 6 t d−1 was required for the model to produce the observed absolute Fe concentration. However, this depends crucially on the rate of vertical transport of the ablated Fe through the MLT to below 80 km, where it is presumed to be permanently removed as MSPs. A previous modelling study of the Na layer12 showed that the size of the IDP input required to model the observed atomic Na layer correlates roughly linearly with the vertical eddy diffusion coefficient (Kzz). Since Kzz is itself a poorly known parameterisation of vertical transport, measurements of metal atom concentrations do not directly constrain the IDP input rate.
Four components of vertical transport in the MLT have been identified recently.66 These are: the residual mean circulation (downwards in winter, reverse in summer); turbulent (eddy) diffusion, produced by breaking gravity waves; downwards dynamical transport caused by dissipating gravity waves; and chemical transport, where wave action and irreversible chemical loss at a lower altitude (e.g. to form MSPs) produces a net flux. A high performance metal resonance lidar has been used to measure the Na atom density and vertical wind profiles simultaneously; the average of their product yields the vertical Na atom flux as a function of height.66 The annual-average downward Na flux measured at the Starfire Optical Range (New Mexico) corresponds to an IDP input of about 24 t d−1. Interestingly, this input is ∼4 times higher than that needed for the 1-D simulation in Fig. 6, which only includes vertical transport by eddy diffusion.13 Dynamical and chemical transport appear to be much more important than turbulent transport, so that until these transport mechanisms are included in models it will not be possible to use the measured metal atom concentrations to constrain the IDP input (or at least the ablated input) reliably. Unfortunately, both dynamical and chemical transport are driven by relatively short period/wavelength gravity waves which cannot yet be resolved explicitly in general circulation models.
Nevertheless, the relative concentrations of the different meteoric metal atoms provide valuable information. Now that the relevant ion-molecule and neutral reactions of Na-, Fe-, Mg- and Ca-containing species have been studied in detail in the laboratory,6,12,56,64,67–71 the atmospheric chemistries of these metals can be reasonably well quantified in a model. This means that more robust conclusions can be drawn regarding the relative meteoric ablation rates of Na, Fe, Mg and Ca that are required to reproduce lidar measurements of the metals. In the case of Na and Fe, modelling a set of lidar observations of both metals at South Pole showed that the Fe ablation rate, relative to that of Na, needed to be reduced by a factor of ∼4.72 That is, the Fe:
Na ablation rates need to be ∼4
:
1, compared to their chondritic ratios of 15
:
1. The Mg ablation rate needs to be decreased relative to that of Na by a factor of ∼5, while Ca ablation rate needs to be decreased by a factor of ∼40.57 These relative ablation factors therefore increase as the element becomes more refractory. As discussed in ST1, significant differential ablation will occur if most of the incoming IDs are small and/or slow, i.e. in close to prograde orbits. Alternatively, the evaporation coefficients39 for the refractory elements in the Herz–Knudsen equation (eqn (III)) could be much less than unity.
However, this creates a potential difficulty. If there is a population of meteoroids that lose all their volatile elements (Na and K), but do not then ablate completely, this implies that the particles must have melted (diffusion of Na and K through the solid particles would be too slow to allow them to escape into the gas phase). If a particle of mass greater than 10−9 g survives atmospheric entry, then it should sediment rapidly to the earth's surface (within a day). Hence, there should be a fairly homogeneous scattering of once-molten IDPs at the surface. These particles are termed cosmic spherules, and can be identified in polar ice because they are close to perfectly spherical and glassy, having melted during atmospheric entry.
The accretion rate of cosmic spherules has been measured by retrieving them from the bottom of an ice chamber used for a drinking water tank at the Scott–Amundsen base at South Pole.73 The flux and size distribution of 50–700 μm diameter particles, corresponding to particles in the mass range from 0.2 to 500 μg, were used in the study. This mass range covers the bulk of IDP particles (Fig. 2), so these cosmic spherules should provide a useful measure of the bulk of IDPs which underwent partial ablation. Furthermore, because the spherules sediment rapidly through the atmosphere, they are not subject to atmospheric circulation, as in the case of the much smaller MSPs (see ST 5). Thus, the flux measured at South Pole should be a good measure of the flux of unablated IDP material. This flux corresponds to a global input of 7.4 ± 0.9 t d−1.73
If the average meteoroid entry velocity peaks around 25 km s−1 based on radar measurements (see ST 1), then about 90% of the incoming mass should ablate.5,73 This would imply that the total IDP input was around 74 t d−1. However, that would not explain why a fraction of only ∼0.25 of the incoming Fe and Mg ablated relative to Na. Since these are the major meteoric constituents, the implication would be that the total IDP input is only around 7.4/(1–0.25) = 10 t d−1. That in turn would imply that the average entry velocity is less than 15 km s−1i.e., most of the IDPs are in a prograde orbit, in accord with the ZCM.1
An important reason for studying MSPs in the mesosphere is their relation to noctilucent clouds (NLCs). These clouds were first reported in 1885, and have been growing brighter and spreading to lower latitudes through much of the last century, so that they appear to be a clear signal of climate change.79 NLCs occur between 80 and 86 km, at high latitudes in the summer where the temperature falls below 150 K and H2O vapour, which is present at mixing ratios of only a few parts per million, is then able to form ice particles spontaneously.7 An important uncertainty in NLC research is the nature of the nuclei on which the ice particles grow. Understanding this is important because changes to the dominant meridional circulation in the mesosphere (which is driven by gravity waves from the lower atmosphere6) may alter the supply of nuclei, which would then affect both the occurrence frequency and brightness of the clouds. Furthermore, the increasingly consistent estimates of ice cloud particle numbers obtained from lidar, radar and satellite observations7 can be linked back to the MSP number density and hence to the IDP input.
Electronic structure calculations have recently been used to demonstrate that the smallest MSPs which should act as ice nuclei are the metal silicate molecules FeSiO3 and MgSiO3.80 This is because MgSiO3 and FeSiO3 have extremely large electric dipole moments of 12.2 and 9.5 Debye, respectively, so that H2O molecules bind to them with large negative free energies. The hydration thermodynamics indicate that ice nucleation should occur at a temperature around 140 K for a H2O mixing ratio of 4 ppm, typical of the polar summer mesosphere where NLCs form.
An important constraint on the IDP input is the optical extinction caused by MSPs between about 40 and 75 km. It has recently become possible to measure these small extinctions (as low as 10−8 km−1), using a visible/near-IR spectrometer on the Aeronomy of Ice in the Mesosphere (AIM) satellite.15 Extinction measurements indicate that the MSP composition is probably olivine (Mg2xFe2−2xSiO4, x = 0.4). However, the refractive indices used for this study were for bulk crystalline minerals, whereas MSPs are nm-size (and likely amorphous) particles; very small particles often have quite different refractive indices.86 The particles are also likely to be chemically weathered by H2O and H2SO4 during the months they spend descending into the stratosphere; this will probably lead to hydroxide and sulphate groups on the particle surfaces which again may change their optical properties significantly. Measured refractive index data on realistic particles are needed to relate the observed atmospheric extinctions to the MSP volume densities, and hence through a general circulation model to the IDP input.
Metal-rich MSPs should readily remove acidic species (e.g. H2SO4, HCl, HNO3) from the gas phase. This may explain the unexpected decrease of H2SO4 measured by mass spectrometry on balloons in the upper stratosphere.87–89 Indeed, a recent study83 using the UK Met Office's Unified Model has shown that the observed removal of H2SO4 could be explained if the uptake coefficient of H2SO4 on MSPs is greater than 10−2 (for an IDP input of around 20 t d−1). If the uptake coefficient of H2SO4 on MSP analogue particles could be determined in the future, this type of modelling exercise could be reversed to obtain the volumetric surface area of MSPs in the upper stratosphere, and hence the IDP input.
Airborne flights of an aerosol mass spectrometer in the mid-latitude lower stratosphere have shown that sulphate particles contain ∼0.75 wt% and ∼0.2 wt% of meteoric Fe and Mg, respectively.14 These fractions are even higher inside the winter polar vortices.84 A recent laboratory study83 showed that amorphous Fe-Mg-silicate particles dissolve in concentrated H2SO4 solutions at temperatures down to 230 K (typical lower stratospheric temperatures) on a time scale of less than a week. In fact, the very high concentrations of Fe and Mg measured by the airborne mass spectrometer14 imply that these metals are mostly in the form of solid sulphate particles in the droplets. Solid particles can act as efficient heterogeneous nuclei; this may explain an earlier laboratory study90 which found that concentrated binary H2SO4-H2O solutions containing FeSO4 and MgSO4 freeze to form sulphuric acid tetrahydrate (SAT) between 12 and 20 K higher than supercooled pure solutions.
Explaining nitric acid trihydrate (NAT) formation from the tertiary HNO3-H2SO4-H2O system in the winter polar stratosphere has been a long-standing problem. Homogeneous nucleation of NAT is too slow to account for observed NAT particles in the polar vortex,91 and MSPs have therefore been proposed as likely heterogeneous nuclei.92 Although experiments have been performed with silica particles93 and ground-up meteorites,94 neither of these surrogates is representative of the amorphous, fractal-like nature of MSPs,65,85 so that the nucleating ability of realistic MSP analogues should be investigated.
MSPs could therefore modify the O3 depletion in the lower stratosphere resulting from chlorine activation on PSCs and de-nitrification.92 A detailed understanding of interactions between MSPs and stratospheric aerosols will be important for accurate predictions as the stratosphere cools over the next century.95 Also, in the context of proposed geo-engineering plans to increase the sulphate aerosol by pumping SO2 into the stratosphere (in order to increase the amount of solar radiation scattered directly back to space, thus counteract greenhouse gas-driven warming in the troposphere),96 a quantitative assessment should be made of the possible effects caused by meteoric debris.
How can these very high fluxes be reconciled with the estimates from within the atmosphere, which seem to be consistent with a flux of less than 50 t d−1 (Table 1)? Interpreting the ice core flux measurements requires understanding the transport of MSPs into the troposphere, and their subsequent deposition mechanisms. For instance, in the study18 of Ir/Pt in the Greenland ice core it was postulated that MSPs descended from the mesosphere into the troposphere within the winter polar vortex, causing a localised concentration of MSPs at polar latitudes and hence the high deposition flux. However, a recent study83 of MSP transport using the Unified Model indicates that once MSPs reach the lower stratosphere in either polar vortex, they mostly enter the troposphere through mid-latitude tropopause folding. Dry deposition should then be reasonably uniform (within a factor of 2) over the entire earth's surface.
There are two problems with a uniform distribution of MSP deposition. First, it implies that the IDP flux really is much greater than 200 t d−1. Second, it does not explain why the deposition in Greenland is so much higher than Eastern Antarctica. The answer to the second point is that the snowfall rate in central Greenland is about 7 times greater than the Antarctic interior, implying that wet deposition is a more important removal mechanism for MSPs than dry deposition.97 Assuming that this is the case, the deposition map in Fig. 7 was generated by simply multiplying the precipitation rate by the MSP concentration in the lowest layer of the Unified Model, and then scaling to the Greenland ice core flux. Inspection of the figure shows that the deposition flux in central Greenland is about 7 times larger than in eastern Antarctica, in agreement with the measurements. However, the global IDP input is only reduced to 210 t d−1, so this serious discrepancy with the atmospheric estimates in Table 1 remains.
![]() | ||
Fig. 7 Annual-averaged flux of MSPs to the earth's surface by wet deposition predicted by the Unified Model (Units: 10−11 g cm−2 d−1). Plot provided by S. Dhomse (University of Leeds). |
Fig. 7 shows that the model predicts a large MSP deposition flux into the Southern Ocean around the coastal shelf of Antarctica (particularly in the Pacific sector), where the supply of bio-available iron to phytoplankton is limited.8 The estimated input from the Unified Model is ∼3 μmol Fe m−2 y−1, compared with an Aeolian dust input of ∼30 μmol Fe m−2 y−1.99 However, unlike continental mineral dust which has a low solubility (estimates vary from <1 to 10%), the MSP Fe should be in the form of highly soluble Fe2SO4 after processing in the stratospheric sulphate layer (see ST4 above). Thus, the input of bio-available Fe from IDPs is likely to at least as large as (and perhaps an order of magnitude greater than) the Aeolian dust input. This could have significant climate implications because phytoplankton both draw down CO2 and produce dimethyl sulphide, which evades into the atmosphere and contributes to the formation of ultra-fine aerosol which may grow large enough to act as cloud condensation nuclei.100
• 5–10 t d−1: this range is supported by: models of the metal layers in the MLT (although these models only include vertical transport by eddy diffusion); and the micrometeorite flux if the average entry velocity <15 km s−1, which then also explains the differential ablation of Ca and Fe relative to Na.
• 20–50 t d−1: this range is supported by: the micrometeorite flux if the average entry velocity >20 km s−1 (although then differential ablation is limited); MSP optical extinction in the middle atmosphere; the meteoritic content of stratospheric sulphate particles; and probably by MLT models if dynamical and chemical transport terms are included.
• 100–300 t d−1: this range is supported by spaceborne IDP detection (the LDEF experiment), and by the ZCM, although in both cases the estimates are reduced into the middle range if constrained by meteor radar measurements of entry velocity; the high range is also supported by surface deposition measurements, though there is uncertainty concerning transport into the troposphere and deposition mechanisms.
Finally, it is worth emphasising two points. First, the ablation of cosmic dust plays a number of significant roles throughout the atmosphere from the thermosphere to the surface, and some of these may be even more important if the IDP input is close to the upper end of current estimates. Second, even if the IDP input is close to the lower end of the range, this indicates that there are important gaps in understanding the evolution of dust in the solar system, and the transport processes coupling the different regions of the atmosphere.
CABMOD | Chemical Ablation Model |
DR | dissociative recombination (with electrons) |
HPLA | high performance large aperture (radar) |
IDP | interplanetary dust particle |
LDEF | Long Duration Exposure Facility |
MIF | meteoric input function |
MLT | mesosphere/lower thermosphere |
MSP | meteoric smoke particle |
NAT | nitric acid trihydrate |
NLC | noctilucent cloud |
PSC | polar stratospheric cloud |
SAT | sulphuric acid tetrahydrate |
ST | Science Topic |
UM | Unified Model |
ZCM | Zodiacal Cloud Model. |
Footnote |
† Part of the atmospheric chemistry themed issue. |
This journal is © The Royal Society of Chemistry 2012 |