Global Electric Circuit of Mars

Variations in Conductivity with Height

Above the ionization peak due to cosmic rays (on Earth, this peak occurs at approximately 10 km; at Mars, this peak occurs at the surface), the ion pair production rate drops exponentially. This is due to the fact that the ion production rate is directly proportional to the atmospheric number density. The pressure, and, therefore, the number density (if temperature is assumed to be approximately constant) of an atmosphere in hydrostatic equilibrium will decrease exponentially with altitude with an e-folding distance of one scale height, H. H = kT/mg where m is the average mass of the atmospheric molecules and g is the local acceleration due to gravity. The scale height in the middle atmosphere of Earth is about 7 km; at Mars the scale height is between 10 and 11 km. This difference is mainly due to the lower gravitational acceleration on Mars.

In Earth's middle atmosphere, above the ionization peak due to cosmic rays, the atmospheric conductivity nearly exponentially increases with altitude [Holzworth, 1995]. The main reason for this phenomenon is the presence of the collision frequency in the denominator of the expression for conductivity given above. Therefore, the mean free path of the ions increases as the neutral number density drops exponentially with altitude, and the collision frequency decreases, which more than makes up for the drop off in ionization rate with altitude.

As a simple example, consider an isothermal atmosphere in hydrostatic equilibrium. Typically, the change in ion number density, n, with time is given by

dn/dt = S - an2 - ßNAn

where S is the ion pair production rate, a is the recombination coefficient, ß is the attachment coefficient, and NA is the number density of aerosols. Ignoring the effect of aerosols and solving for a steady state solution,

n = SQRT(S/a)

Since S decreases with atmospheric number density above the ionization peak,

S ~ exp(-z/H)

and the ion number density

n ~ exp(-z/2H)

Now the collision frequency must be proportional to the ion number density and the neutral number density, so

vni ~ exp(-z/2H) x exp(-z/H) = exp(-3z/2H)

Finally,

õ = nq2/mvni ~ exp(-z/2H) / exp(-3/2H) = e(z/H)

So this result states that the atmospheric conductivity should exponentially increase with altitude with the same scale height, H, with which the atmospheric number density decreases. However, Gringel et al. [1986] showed that the mobility, k, is inversely proportional to the neutral number density, nn. Therefore,

õ = nq2/mvni = nqk ~ n/nn = exp(-z/2H) / exp(-z/H) = e(z/2H)

and conductivity increases more slowly. Regardless of its exact form, both methods result in an exponentially increasing atmospheric conductivity.

These estimates are in actuality upper limits on conductivity. In this simple treatment, we have neglected the effects of aerosols. Aerosols essentially decrease conductivity by vastly increasing the mass of the charge carriers. Typically aerosols such as dust are found in greater quantities near the surface, and, therefore, lower the atmospheric conductivity near the surface.

Briefly, we can estimate the total resistance of the Martian atmosphere if we make some simplifying assumptions. First, assume that the atmospheric conductivity at the surface is 10-11 S/m as given by Grard [1995] and exponentially increases with altitude only due to ion pair production by cosmic rays. Secondly, assume that the scale height of the atmospheric conductivity is the same as the neutral atmosphere scale height, that is, 10 km. Finally, assume that the conducting ionosphere occurs at an altitude of 100 km. Then

specific resistance = p = 1/õ

columnar resistance = P = § p dz = Hõ[exp(-zo/H) - exp(-z/H)] ~ 1015 Ohms m2

(Here § is not meant as a contour intergral but is used due to HTML restictions.)

resistance = R = P/4 pi rplanet = P/4 pi (3400 km)2 ~ 10 Ohms

Therefore, the global resistance of the Martian atmosphere appears to be on the order of 10 Ohms. Aerosols may increase the specific resistance of the lower atmosphere, and , thus, increase this value. This "back of the envelope" estimate is probably within an order of magnitude or so of the true value.

So far we have shown that the Martian atmosphere should have a finitely conducting atmosphere as a result of galactic cosmic rays, radioactive crustal material, and ionizing solar ultraviolet radiation which reaches the surface. The total atmospheric conductivity can vary on daily (due to solar ultraviolet radiation), yearly (due to global atmospheric pressure variations), and 11 year (due to the solar cycle) cycles. Also, the atmosphere is located between the highly conducting surface and ionosphere. Finally, the atmospheric conductivity increases exponentially with altitude.


Previous: Introduction and Sources of Atmospheric Conductivity
Next: Possible Current Sources

Back to Matt's Home Page


Global Electric Circuit of Mars
23 June 1998
matt at ssl dot berkeley dot edu