meteoric ablation is the source of metals such as calcium , sodium , and iron in the mesosphere and lower thermosphere ( mlt ) [ plane et al . ,
these metals occur as neutral atoms in layers between 80 and 100  km and atomic ions which become dominant above 95  km .
lidar observations of the neutral metal layers have shown that ca is significantly depleted with respect to na and fe : the ca / na and ca / fe column abundance ratios are ~0.01 and 0.004 [ plane et al . ,
2015 ] , compared to their chondritic ratios ( in ci type carbonaceous chondrites ) of 1.05 and 0.069 , respectively [ asplund et al . , 2009 ] .
the most likely reason for this depletion is differential ablation , where the highly refractory ca only ablates from meteoroids which reach temperatures above 2200  k , compared to the relatively volatile elements na and fe which largely ablate below 2000  k [ carrillosanchez et al . ,
however , in contrast to the neutral ratios , the observed atmospheric ca / na and ca / fe ratios are ~0.4 and 0.03 [ kopp , 1997 ] , which are much closer to their respective ci ratios .
this result implies that there are important differences in the ionmolecule chemistries of these metals . in this paper
we focus on calcium and iron chemistries , because ca and fe are chemically reactive ions and therefore behave more alike , compared with the chemically inert na ion . in the mlt , these ions are largely produced by charge transfer reactions with no and , to a lesser extent , o2 
  ( photoionization is a much less important ionization process [ plane et al . ,
2015 ] ) . for example [ rutherford et al . , 1972 ] , 
 ( r1)ca+no+ca++no 
  above 110  km ,
the slow dielectronic / radiative recombination of the metallic ions with electrons becomes the dominant neutralization pathway [ bones et al .
, 2016 ] . however , between 90 and 100  km the dominant pathway is reaction with o3 to produce the metal oxide ion [ broadley et al . , 2007 ] , 
 ( r2)ca++o3cao++o2followed by dissociative recombination ( dr ) with electrons 
 ( r3)cao++eca+o 
  since dr reactions generally have rate coefficients > 10  cm  molecule  s [ florescumitchell and mitchell , 2006 ] , if reactions r2 and r3 were the only chemistry occurring , the lifetime of ca would be ~10  s , with consequently much lower concentrations than measured by lidar [ gerding et al . , 2000
; granier et al . , 1985 ; raizada et al . , 2011 ] and rocketborne mass spectrometry [ kopp , 1997 ] .
in fact , the reaction of cao with atomic o 
 ( r4)cao++oca++o2converts the metal oxide ion back to the atomic ion , thus slowing the overall rate of neutralization [ broadley and plane , 2010 ] . the rate coefficients k 
 1 [ rutherford et al . , 1972 ] , k 
 2 [ broadley et al . , 2007 ] , and k 
 4 [ broadley and plane , 2010 ] have been measured previously , as have the rate coefficients for the analogous reactions of fe [ rollason and plane , 1998 ; rutherford and vroom , 1972 ; woodcock et al . ,
we have recently reported the first measurement of the dr of feo [ bones et al . , 2016 ] .
the present paper describes a study of reaction r3 , whose rate coefficient does not appear to have been reported previously .
we then use the measured value of k 
 3 to understand why the ca / fe ratio between 90 and 100  km is nearly an order of magnitude larger than the ca / fe ratio .
finally , we examine the role of ionmolecule chemistry in the phenomenon of sporadic ca layers ( cas ) .
these are very thin layers of ca which appear suddenly , generally at altitudes around 95  km , and can have peak concentrations more than an order of magnitude above the background average ca atom layer [ gerding et al . ,
simultaneous commonvolume measurements of ca and ca by lidar provide a unique tool for determining whether cas are produced by the neutralization of ca in descending sporadic e layers , as proposed originally to explain sporadic na layer formation from na [ cox and plane , 1998 ] .
the measurements were made in a flowing afterglow system with a pulsed laser source of metallic ions .
the system is described in detail elsewhere [ bones et al . , 2016 ] .
briefly , the system consists of a stainless steel flow tube connected to a differentially pumped quadrupole mass spectrometer ( hiden hpr60 ) running in positive ion mode with pulse counting .
an ar / electron plasma is generated by flowing he through a microwave discharge cavity ( mccarroll cavity , 145  w at 2450  mhz ) and then introducing a flow of ar into the helium / electron plasma , whereupon the charge is transferred to the ar largely through penning ionization .
the electron density was measured 12.5  cm upstream of the target by a langmuir probe ( hiden espion ) .
all experiments were conducted at pressures of 1  torr and a temperature of 295  k. cao ions ( and neutral cao molecules ) were introduced into this plasma by the pulsed laser ablation of a solid target .
a neodymium : yttrium / aluminum / garnet laser ( 8  mj  pulse at 532  nm , 10  hz repetition rate ) was directed at a sample of ca metal with a heavily oxidized surface , mounted on a rotating rod .
isotopic analyses of the ablated ions showed negligible contamination of the signal by fe ( which is isobaric with cao at m / z  =  56 ) .
for some experimental runs , the laser pulse was loosely focused on the target to change the ratio of cao to cao generated in the ablation pulse .
the relative concentration of cao was measured by the mass spectrometer as a function of flow rate ( 35  slm ) , which is inversely proportional to the time available for the cao ions to react with the electrons in the plasma before they enter the mass spectrometer ( 2.13.5  ms , corresponding to flow velocities of 69113  m  s ) .
the skimmer cone between the flow tube and the first stage of the mass spectrometer had a 4.0  mm orifice biased at 19  v , and the skimmer cone between the first and second stages of the mass spectrometer had a 1.8  mm orifice biased at 89  v. each cao measurement was the sum of the m / z  =  56 signal from three repeat accumulations of 500 laser shots each .
these repeat measurements had an average standard deviation of 15% , although for measurements at low flow velocities the standard deviation could be as high as 30% . during the course of an experiment , over typically 23  h , the cao ablation efficiency varied slowly ( usually decreasing ) , and this was corrected for by taking reference cao signals periodically ( at a median flow velocity of 91  m  s ) and fitting a normalization function .
materials : he ( 99.995% , boc gases ) was purified by passing through a trap of dry molecular sieve ( 4   , 12  mm , alfa aesar ) cooled to 77  k. ar ( 99.9999% , boc gases ) was used without further purification .
the calcium shot ( sigmaaldrich ) was at least a decade old and extensively oxidized .
the ablation of the oxidized ca target produced pulses of cao and cao ( figure  s1 in the supporting information ) .
the cao then undergoes charge transfer with ar ions in the plasma : 
 ( r5)cao+ar+cao++ar 
  the upper limit for the rate coefficient k 
 5 can be calculated from modified langevin theory [ su and chesnavich , 1982 ; troe , 1985 ] to be 8.2    10  cm  molecule  s , using a theoretical dipole moment of 8.8  d and polarizability of 6.6    10  cm for cao [ broadley and plane , 2010 ] .
the interactions of cao , cao , ar , and e in the flow tube are then described by the following coupled differential equations : 
 ( e1)dcao+dt = k5ar+caokdiffcao++k3ecao+ 
 ( e2)dcaodt=kdiffcao+k5ar+cao 
 ( e3)dar+dt=kdiffar++k5caoar+ 
 ( e4)dedt=kdiffe+k3cao+e 
  the average electron diffusion rate k 
 diff 
  was 948    96  s at 1.0  torr .
this was determined by measuring the absolute electron density with the langmuir probe and varying the flow velocity at constant pressure ( figure  s2 ) . the electron density immediately after formation of the ar plasma
was then estimated to be ( 3.2    0.6 )    10  cm . in the diffusion rate of ar , k 
 diff 
 , the dominant positive charge carrier was assumed to match the electron diffusion rate under the conditions of ambipolar diffusion which apply because the electrons and ions rapidly equilibrate to the bath gas temperature at the relatively high total pressure of 1  torr [ smirnov , 2001 ] ) .
k 
 diff 
  was estimated to be 3.7% smaller than k 
 diff 
  in the bath gas where the ar / he ratio was 1:10 [ bones et al . ,
the diffusion coefficients of neutral cao are calculated ( using the method described in self and plane ) to be 464 and 95  torr  cm  s in he and ar , respectively , which corresponds to a wall loss rate k 
 diff 
  of 651  s. relations e1e4 were solved numerically to predict the relative cao concentration at the downstream end of the flow tube .
there are two unknowns : k 
 3 and the initial neutral cao concentration ( [ cao]0 ) .
each data set was then normalized to the cao produced by pulsed ablation at the fastest flow velocity , in order to remove the daytoday variation in the ablation efficiency from the ca target .
the data were then averaged in six reaction time bins with an average bin width of 1.4    10  s , and a neldermead iterative simplex algorithm [ nelder and mead , 1965 ] used to find the best fit of the kinetic model to the six experimental points , by varying k 
 3 and [ cao]0 .
the fitting was weighted to the number of data in each time bin , which varied between 10 and 32 .
note that in order to model the timeresolved cao concentration , the model actually requires the product k 
 5 [ cao ] in relations e1 and e2 ; since we fix k 
 5 to the calculated langevin limit ( see above ) , [ cao]0 is then the other unknown fitted parameter .
figure  1a is a contour plot of  
  as a function of k 
 3 and [ cao]0 , where 
 ( e5)2=i=1ncaoexpt , i+caomodel , i+caoexpt , i+2 
  ( a ) plot of the  
  surface as a function of k 
 3 and initial neutral cao concentration produced by pulsed ablation of the ca target .
the white cross indicates the position of the minimum , and the white contour the 1 envelope of uncertainty .
( b ) abundance of cao as a function of time in the flow tube between the ablation source and detection at the mass spectrometer : measurements ( solid circles with 1 error bars ) , model fit ( solid line ) , and model run with the upper and lower limits of k 
 3 ( dashed lines ) .
the minimum on the  
  surface yields k 
 3(295  k )  =  ( 3.0    1.0 )    10  cm  molecule  s , where the uncertainty estimate is discussed below .
the corresponding value for [ cao]0 is 1.9    10  molecule  cm .
the ratio of [ cao]0/[cao]0 produced directly by ablation is estimated to be 17    5 .
this ratio decreased by 45% when the ablation laser was loosely focused , as expected because a greater fraction of ions would have been produced .
the uncertainty in k 
 3 was estimated by using monte carlo selection of the theoretically calculated and experimentally measured parameters used in the kinetic model , within their uncertainties listed in table  1 ( assuming a tophat probability distribution ) .
the model was then rerun to compare with the measurements : if the simplex algorithm found a minimum in the  
  surface , and the resulting value of k 
 3 was within the range of 1    101    10  cm  molecule  s expected for a dr reaction [ florescumitchell and mitchell , 2006 ] , then the value of k 
 3 was accepted .
the mean and standard deviation of k 
 3 from 2500 model fits then provided the final value of k 
 3 and the estimate of its uncertainty .
although the determination of k 
 3 is sensitive to the flow velocity and electron diffusional wall loss , these parameters are measured with high precision .
in contrast , k 
 3 is much less sensitive to the estimated [ cao]0/[cao]0 ratio and estimate of k 
 5 using langevin theory , but these parameters have a greater assigned uncertainty ( table  1 ) and so make a similar contribution to the overall uncertainty of k 
 3 ( see the supporting information for further details ) .
the white contour line in figure  1a depicts the envelope of uncertainty around the minimum , which is 44% below the contour line .
figure  1b shows a plot of the measured abundance of cao ( with 1 standard deviation ) as a function of reaction time .
the dashed lines show the model output with the upper and lower bounds of k 
 3 .
note that the overall uncertainty is larger than the uncertainty of the measurements because the monte carlo analysis included additional uncertainties ( table  1 ) when determining the uncertainty of k 
 3 .
in order to study the ionmolecule chemistry of calcium in the mlt , we assign a t 
  temperature dependence to k 
 3 , which is typical for dr reactions [ florescumitchell and mitchell , 2006 ] .
reaction r3 is 45% slower than the dr reaction of feo at room temperature [ bones et al . ,
we now examine the consequences of this for explaining the very different ca / fe and ca / fe ratios in the 90100  km region of the atmosphere .
the dashed black line shows the [ ca]/[fe ] ratio profile measured by rocketborne mass spectrometry .
the mean and standard deviation were calculated from seven flights covering a range of latitudes and local times ( note that the geometric , rather than arithmetic , mean is used because of the small data set ; the geometric mean is defined as 
 x=x1x2xnn ) .
the dashed blue line shows the mean [ ca]/[fe ] ratio computed from simultaneous and common volume lidar measurements of ca and fe at khlungsborn , germany ( 54n ) .
data from 22 nights ( each with 29  h of measurements ) were used ; the total data set is ~120  h of measurements . the measurements were performed between april 1997 and march 1999 .
vertical profiles ( 90100  km ) of the measured ca / fe ratio from a series of middle and highlatitude rocket flights ( dashed black line ) , and the ca / fe ratio measured by lidar at 54n ( dashed blue line ) , from which the ca / fe ratio is estimated by assuming chemical steady state ( black solid line ) .
the ci ratio of the metals ( red dashed line ) and the modeled meteor ablation ratio ( green dashed line ) are shown for comparison .
these two profiles should be compared with the chondritic ci ratio ( red dashed line ) [ asplund et al . ,
2009 ] , and the ratio of the meteoric ablation fluxes of ca and fe , calculated by applying the chemical ablation model to an astronomical model of cosmic dust [ carrillosanchez et al . , 2015 ] .
note that the ablation ratio is a factor of 7 times smaller than the ci ratio because of the differential ablation of the two elements .
the [ ca]/[fe ] ratio falls between the ci and the ablation ratios , and the [ ca]/[fe ] ratio is significantly below even the ablation ratio . since the rates of reactions r1r4 in the mlt are reasonably fast between 90 and 100  km , and the dominant forms of these metals in this region of the atmosphere are neutral and ionized atoms [ plane et al . ,
2015 ] , the kinetic steady state approximation can be applied to ca and cao : 
 ( e6)dca+dt = k1no+cak2o3ca++k4ocao+0 
 ( e7)dcao+dt = k2o3ca+k3ecao+k4ocao+0 
  hence , 
 ( e8)ca+ca = k1no+k3e+k4ok2k3o3e 
  since k 
 4[o ] is generally more than 3 orders of magnitude larger than k 
 3 [ e ] between 90 and 100  km [ plane et al . ,
2015 ] , 
 ( e9)ca+cak1no+k4ok2k3o3e 
  dividing equation ( e9 ) by the analogous equation for [ fe]/[fe ] , the atmospheric constituents [ no ] , [ o3 ] , [ o ] , and [ e ] cancel . rearranging , and using the rate coefficients listed in table  2 , [ ca]/[fe ]  =  8.3[ca]/[fe ] .
when the lidar measurement of the [ ca]/[fe ] ratio is multiplied by this factor , the estimated [ ca]/[fe ] ratio ( solid black line in figure  2 ) is seen to be in good agreement with the ion ratio measured by mass spectrometry .
inspection of table  2 shows that there are two reasons for this large increase in the ratio of the ions compared with the neutrals .
first , charge transfer between no and ca is ~5 times faster than for fe .
rate coefficients required to estimate [ ca]/[fe ] from [ ca]/[fe ] we now examine whether the mechanism for sporadic metal layer formation in a descending sporadic e layer [ cox and plane , 1998 ] can explain the appearance of cas below 100  km .
figures  3a and 3b show the lidar observations of ca and ca on 9 august 1995 at juliusruh on the island of rgen , germany ( 55n )
. in figure  3a , a pronounced layer of ca ions , confirmed by ionosonde measurements to be part of a sporadic e layer ( shown with white crosses in the figure ) , descends from 100  km to 93  km at a rate of ~4  km  h. figure  3b illustrates the simultaneous neutral ca measurements , with the background ca layer continuously present with a peak around 87  km .
note that a cas layer appears at 0140 lt once the ca layer has descended to 94  km .
in contrast , the patches of relatively high concentrations of ca above 95  km that occur around 0100 and 0200 lt are not accompanied by significant concentrations of ca . lidar observations of ( a ) ca and ( b ) ca on 9 august 1995 at juliusruh , germany . the white crosses in figure 
( c ) modeled production of ca from the ca ions in a sporadic e layer whose descent is indicated by the dashed white line .
we now use the model of cox and plane   to simulate the production of ca in this descending sporadic e layer .
the sporadic e layer is initialized with the measured ca density of 810  cm and fe , mg , na , and si densities which are scaled to this using their relative abundance factors measured by kopp .
the total metal ion ( and hence electron ) density at 105  km is then 8.9    10  cm , consistent with the ionosonde electron density measurements which ranged from ( 4.5 to 10.5 )    10  cm between 2300 and 0000 lt .
the layer is set to have a full width at half maximum of 1.5  km .
note that the sporadic layer is assumed to have a large horizontal homogeneous extent so that advection does not play a role while the layer is descending .
the ionmolecule chemistry of ca from the present study and the other four meteoric metals [ plane et al . ,
2015 ] are included in the model , along with background atmospheric densities of o3 , o , total density , and temperature from the 1d general mesospheric model mesomod [ murray and plane , 2005 ] .
figure  3c shows the modeled production of ca in the descending sporadic e layer ( the background ca layer is not modeled in this simulation ) .
the layer , whose position is indicated by the dashed white line , descends at a constant rate to 93  km where the ions are then dumped , as indicated by the ca observations in figure  3a .
the model successfully reproduces the height and absolute concentration of the resulting cas , which is an important test of the reaction kinetics , since inspection of figure  3 shows that only ~10% of the available ca is converted to ca during the observational period .
moreover , very little neutral ca is produced until the ion layer descends below 95  km , consistent with figure  3b , because the rate of neutralization of ca depends on [ o3 ] and inversely on [ o ] ( equation e9 ) . since [ o3 ]
increases ~100fold between 105 and 93  km , and [ o ] is approximately constant , these ion neutralization rates increase by roughly 2 orders of magnitude .
this causes the production of ca to occur at the dump altitude , but not higher up 
so long as the rate of descent of the sporadic e layer is sufficiently fast .
in order to study the ionmolecule chemistry of calcium in the mlt , we assign a t 
  temperature dependence to k 
 3 , which is typical for dr reactions [ florescumitchell and mitchell , 2006 ] .
reaction r3 is 45% slower than the dr reaction of feo at room temperature [ bones et al . ,
we now examine the consequences of this for explaining the very different ca / fe and ca / fe ratios in the 90100  km region of the atmosphere .
the dashed black line shows the [ ca]/[fe ] ratio profile measured by rocketborne mass spectrometry .
the mean and standard deviation were calculated from seven flights covering a range of latitudes and local times ( note that the geometric , rather than arithmetic , mean is used because of the small data set ; the geometric mean is defined as 
 x=x1x2xnn ) .
the dashed blue line shows the mean [ ca]/[fe ] ratio computed from simultaneous and common volume lidar measurements of ca and fe at khlungsborn , germany ( 54n ) .
data from 22 nights ( each with 29  h of measurements ) were used ; the total data set is ~120  h of measurements . the measurements were performed between april 1997 and march 1999 .
vertical profiles ( 90100  km ) of the measured ca / fe ratio from a series of middle and highlatitude rocket flights ( dashed black line ) , and the ca / fe ratio measured by lidar at 54n ( dashed blue line ) , from which the ca / fe ratio is estimated by assuming chemical steady state ( black solid line ) .
the ci ratio of the metals ( red dashed line ) and the modeled meteor ablation ratio ( green dashed line ) are shown for comparison .
these two profiles should be compared with the chondritic ci ratio ( red dashed line ) [ asplund et al . ,
2009 ] , and the ratio of the meteoric ablation fluxes of ca and fe , calculated by applying the chemical ablation model to an astronomical model of cosmic dust [ carrillosanchez et al . , 2015 ] .
note that the ablation ratio is a factor of 7 times smaller than the ci ratio because of the differential ablation of the two elements .
the [ ca]/[fe ] ratio falls between the ci and the ablation ratios , and the [ ca]/[fe ] ratio is significantly below even the ablation ratio . since the rates of reactions r1r4 in the mlt are reasonably fast between 90 and 100  km , and the dominant forms of these metals in this region of the atmosphere are neutral and ionized atoms [ plane et al . ,
2015 ] , the kinetic steady state approximation can be applied to ca and cao : 
 ( e6)dca+dt = k1no+cak2o3ca++k4ocao+0 
 ( e7)dcao+dt = k2o3ca+k3ecao+k4ocao+0 
  hence , 
 ( e8)ca+ca = k1no+k3e+k4ok2k3o3e 
  since k 
 4[o ] is generally more than 3 orders of magnitude larger than k 
 3 [ e ] between 90 and 100  km [ plane et al . ,
2015 ] , 
 ( e9)ca+cak1no+k4ok2k3o3e 
  dividing equation ( e9 ) by the analogous equation for [ fe]/[fe ] , the atmospheric constituents [ no ] , [ o3 ] , [ o ] , and [ e ] cancel . rearranging , and using the rate coefficients listed in table  2 , [ ca]/[fe ]  =  8.3[ca]/[fe ] .
when the lidar measurement of the [ ca]/[fe ] ratio is multiplied by this factor , the estimated [ ca]/[fe ] ratio ( solid black line in figure  2 ) is seen to be in good agreement with the ion ratio measured by mass spectrometry .
inspection of table  2 shows that there are two reasons for this large increase in the ratio of the ions compared with the neutrals .
first , charge transfer between no and ca is ~5 times faster than for fe .
we now examine whether the mechanism for sporadic metal layer formation in a descending sporadic e layer [ cox and plane , 1998 ] can explain the appearance of cas below 100  km .
figures  3a and 3b show the lidar observations of ca and ca on 9 august 1995 at juliusruh on the island of rgen , germany ( 55n ) . in figure  3a ,
a pronounced layer of ca ions , confirmed by ionosonde measurements to be part of a sporadic e layer ( shown with white crosses in the figure ) , descends from 100  km to 93  km at a rate of ~4  km  h. figure  3b illustrates the simultaneous neutral ca measurements , with the background ca layer continuously present with a peak around 87  km . note that a cas layer appears at 0140 lt once the ca layer has descended to 94  km .
in contrast , the patches of relatively high concentrations of ca above 95  km that occur around 0100 and 0200 lt are not accompanied by significant concentrations of ca . lidar observations of ( a ) ca and ( b ) ca on 9 august 1995 at juliusruh , germany . the white crosses in figure 
( c ) modeled production of ca from the ca ions in a sporadic e layer whose descent is indicated by the dashed white line .
we now use the model of cox and plane   to simulate the production of ca in this descending sporadic e layer .
the sporadic e layer is initialized with the measured ca density of 810  cm and fe , mg , na , and si densities which are scaled to this using their relative abundance factors measured by kopp .
the total metal ion ( and hence electron ) density at 105  km is then 8.9    10  cm , consistent with the ionosonde electron density measurements which ranged from ( 4.5 to 10.5 )    10  cm between 2300 and 0000 lt .
the layer is set to have a full width at half maximum of 1.5  km .
note that the sporadic layer is assumed to have a large horizontal homogeneous extent so that advection does not play a role while the layer is descending .
the ionmolecule chemistry of ca from the present study and the other four meteoric metals [ plane et al . ,
2015 ] are included in the model , along with background atmospheric densities of o3 , o , total density , and temperature from the 1d general mesospheric model mesomod [ murray and plane , 2005 ] .
figure  3c shows the modeled production of ca in the descending sporadic e layer ( the background ca layer is not modeled in this simulation ) .
the layer , whose position is indicated by the dashed white line , descends at a constant rate to 93  km where the ions are then dumped , as indicated by the ca observations in figure  3a .
the model successfully reproduces the height and absolute concentration of the resulting cas , which is an important test of the reaction kinetics , since inspection of figure  3 shows that only ~10% of the available ca is converted to ca during the observational period .
moreover , very little neutral ca is produced until the ion layer descends below 95  km , consistent with figure  3b , because the rate of neutralization of ca depends on [ o3 ] and inversely on [ o ] ( equation e9 ) . since [ o3 ] increases ~100fold between 105 and 93  km , and [ o ] is approximately constant , these ion neutralization rates increase by roughly 2 orders of magnitude .
this causes the production of ca to occur at the dump altitude , but not higher up 
so long as the rate of descent of the sporadic e layer is sufficiently fast .
the rate coefficient for the dissociative recombination of cao has been measured for the first time and is found to be ( 3.0    1.0 )    10  cm  molecule  s at 295  k. this result has two atmospheric implications .
first , the surprising observation that the ca / fe ratio is ~8 times larger than ca / fe between 90 and 100  km can now be explained by the known ionmolecule chemistry of these two metals .
second , the rate of neutralization of ca ions in a descending sporadic e layer is fast enough to explain the often explosive growth of sporadic ca layers .
 abstractthe dissociative recombination of cao+ ions with electrons has been studied in a flowing afterglow reactor . 
 cao+ was generated by the pulsed laser ablation of a ca target , followed by entrainment in an ar+ ion / electron plasma . 
 a kinetic model describing the gasphase chemistry and diffusion to the reactor walls was fitted to the experimental data , yielding a rate coefficient of ( 3.0    1.0 )    107  cm3  molecule1  s1 at 295  k. this result has two atmospheric implications . 
 first , the surprising observation that the ca+/fe+ ratio is ~8 times larger than ca / fe between 90 and 100  km in the atmosphere can now be explained quantitatively by the known ionmolecule chemistry of these two metals . 
 second , the rate of neutralization of ca+ ions in a descending sporadic e layer is fast enough to explain the often explosive growth of sporadic neutral ca layers . 
Introduction
Experimental Technique
Experimental Results
Atmospheric Implications
The Background Ca
Sporadic Ca Layers
Conclusions
Supporting information