← Back
The induction of apoptosis in SGC-7901 cells through the ROS-mediated mitochondrial dysfunction pathway by a Ir(III) complex.
Space Sci Rev (2016) 203:299–345
DOI 10.1007/s11214-016-0241-8
Applications of Electrified Dust and Dust Devil
Electrodynamics to Martian Atmospheric Electricity
R.G. Harrison1 · E. Barth2 · F. Esposito3 · J. Merrison4 · F. Montmessin5 · K.L. Aplin6 ·
C. Borlina7 · J.J. Berthelier5 · G. Déprez5 · W.M. Farrell8 · I.M.P. Houghton6 ·
N.O. Renno7 · K.A. Nicoll1 · S.N. Tripathi9 · M. Zimmerman10
Received: 28 August 2015 / Accepted: 15 February 2016 / Published online: 12 April 2016
© The Author(s) 2016. This article is published with open access at Springerlink.com
Abstract Atmospheric transport and suspension of dust frequently brings electrification,
which may be substantial. Electric fields of 10 kV m−1 to 100 kV m−1 have been observed
at the surface beneath suspended dust in the terrestrial atmosphere, and some electrification has been observed to persist in dust at levels to 5 km, as well as in volcanic plumes.
The interaction between individual particles which causes the electrification is incompletely
understood, and multiple processes are thought to be acting. A variation in particle charge
with particle size, and the effect of gravitational separation explains to, some extent, the
charge structures observed in terrestrial dust storms. More extensive flow-based modelling
demonstrates that bulk electric fields in excess of 10 kV m−1 can be obtained rapidly (in less
than 10 s) from rotating dust systems (dust devils) and that terrestrial breakdown fields can
be obtained. Modelled profiles of electrical conductivity in the Martian atmosphere suggest
the possibility of dust electrification, and dust devils have been suggested as a mechanism
of charge separation able to maintain current flow between one region of the atmosphere
and another, through a global circuit. Fundamental new understanding of Martian atmospheric electricity will result from the ExoMars mission, which carries the DREAMS (Dust
characterization, Risk Assessment, and Environment Analyser on the Martian Surface)—
B R.G. Harrison
r.g.harrison@reading.ac.uk
1
Department of Meteorology, University of Reading, Reading, UK
2
Southwest Research Institute, Boulder, CO, USA
3
INAF – Osservatorio Astronomico di Capodimonte, Naples, Italy
4
University of Aarhus, Aarhus, Denmark
5
Laboratoire Atmosphères, Milieux, Observations Spatiales (LATMOS), Guyancourt, France
6
Department of Physics, University of Oxford, OX1 3RH, Oxford, UK
7
Atmospheric, Oceanic and Space Sciences, University of Michigan, Ann Arbor, USA
8
NASA Goddard, Greenbelt, USA
9
Centre for Environmental Science and Engineering, Indian Institute of Technology, Kanpur, India
10
Johns Hopkins University, Baltimore, USA
300
R.G. Harrison et al.
MicroARES (Atmospheric Radiation and Electricity Sensor) instrumentation to Mars in
2016 for the first in situ electrical measurements.
Keywords Planetary electrostatics · Lightning discharge · Particle electrification · Global
circuit
1 Charge Separation in Dust
It has long been known that dust lofted or transported can become highly electrified. Strong
effects were reported on a gold leaf electrometer connected to a roof-level electrode which
was exposed to a dust storm in Lahore in 1847, together with the generation of sparks (Baddeley 1860). An accompanying characteristic of the Dust Bowl disaster in the US Great
Plains during the 1930s was severe static electricity, audible on domestic medium wave
radios and able to cause electric shocks, with earthing chains necessary on automobiles.
Under fair weather conditions, the magnitude of the vertical terrestrial atmospheric electric
field is about 100 V m−1 , sustained by distant thunderstorm and disturbed weather electrification. Dust storms are considerably electrified in comparison. Quantitatively, electric fields
exceeding 100 kV m−1 have been measured during blowing sand, dust storms and dust devils (Rudge 1913; Demon et al. 1953; Freier 1960; Crozier 1964; Harris 1967; Stow 1969;
Kamra 1972; Schmidt et al. 1998; Renno et al. 2004; Jackson and Farrell 2006; Kok and
Renno 2006, 2008a, 2008b; Williams et al. 2009). As an example, Fig. 1 shows the variation
in electric field during the passage of a dust devil in the Sahara desert, opportunistically
obtained during measurements undertaken for another purpose (Freier 1960).
As discussed later in this paper (Sect. 3), it is generally accepted that dust electrification is
due to contact and triboelectric charging between blowing particles during these phenomena.
The details of the charge transfer mechanism remain incomplete, but several experiments
and observations (Freier 1960; Inculet et al. 2006; Duff and Lacks 2008) suggest that on
average, during collisions, the smallest grains acquire negative charge with respect to larger
particles. In general, the smallest particles are brought into suspension and transported aloft
into the atmosphere by local turbulence, whereas the larger particles stay close to the surface.
This gravitational charge separation is consistent with increase of the atmospheric electric
field observed during dusty phenomena.
Laboratory and theoretical work (Kok and Renno 2006, 2008a, 2008b) suggest that intense electric fields can influence charged particle dynamics through changing their trajectories and reducing the threshold friction speed necessary for initiating their motion. This
facilitates the lifting of particles from the surface, so increasing the concentration of saltating grains at a given wind speed. Electric forces may therefore play a further role in the
evolution of dust events, including dust devils.
Observations of dust electrification are sparse in literature, both because this field of
study is relatively young and because of the difficulty in sampling stochastic dust activity.
There is renewed interest in dust electrification due to its consequences for human activities
such as the breakdown of power transmission lines, electrical alignment of atmospheric particles (Ulanowski et al. 2007), volcanic plume electrification (Mather and Harrison 2006)
and planetary exploration (Helling et al. 2016). In this latter case, the presence of electric discharges or electromagnetic noise can potentially affect communications or communications equipment, but substantial electric fields have also been argued to affect atmospheric chemistry and planet habitability (Atreya et al. 2006; Delory et al. 2006) and the
possible development of life (Miller 1953), in particular for Mars. New measurements from
Applications of Electrified Dust and Dust Devil Electrodynamics. . .
301
Fig. 1 Variation of electric field with time beneath the passage of a dust devil in the Sahara desert. (Reproduced from Freier 1960 with permission of Wiley)
the ExoMars 2016 space mission (Bettanini et al. 2014) using the DREAMS-MicroARES
(acronym formed from Dust characterization, Risk Assessment, and Environment Analyser
on the Martian Surface, and Atmospheric Radiation and Electricity Sensor) instrument will
provide the first direct measurements of the electric properties of Martian atmosphere. This
is likely to mark a hugely important step change in the knowledge of Martian atmospheric
electricity.
This paper therefore presents a pre-MicroARES view of the electrical properties of dusts
and dust devils, with particular reference to Mars. To do this, Sect. 2 first briefly summarises
knowledge on Martian atmospheric electricity. Section 3 then provides an overview of particle electrification mechanisms concerning dust. Section 4 reports the results of some of
the most recent observations of dust electrification during dust storms/devils in terrestrial
deserts both close to the surface and aloft. Section 5 provides an overview of modelling
of dust devils, and Sect. 6 describes the DREAMS-MicroARES sensor to be used on the
forthcoming ExoMars mission.
2 The Martian Electrical Environment
2.1 Atmospheric Electricity on Mars
Atmospheres electrify to varying extents, depending on the charged particles they contain
and whether winds, convection or other meteorological processes can actively separate local
charge (Harrison et al. 2008; Aplin 2006). The consequences of electrification include the
motion of charged particles and ions under the action of electric fields, and electrical discharges where substantial accumulations of charge occur. In the Martian atmosphere, cosmic
rays and ultra violet radiation generate molecular cluster ions which, together with free electrons, make the atmosphere electrically conductive. With active charge separation therefore,
current flow can potentially occur, and a planetary electrical system analogous to the terrestrial global atmospheric electric circuit has accordingly been considered (Fillingim 1998;
302
R.G. Harrison et al.
Farrell and Desch 2001; Aplin 2006). The abundance of dust in Mars’ atmosphere (and remote sensing of its dust devils), combined with knowledge that terrestrial dust electrification
is common has led to the expectation that Martian dust devils can become electrified.
The basis for electrical activity on Mars has been assumed to be similar to electrical activity in Earth’s atmosphere. This has been considered through the global circuit concept introduced for the Earth in the early 20th century by Wilson (1921), Harrison (2011). A global
circuit relies on the existence of several key characteristics: a global generator (lightning,
discharges, etc.), conductive surfaces (the ionosphere and the rocky surface) overall forming a spherical capacitor system (Rycroft et al. 2008, 2012). Lightning is a key indicator
of planetary atmospheric electricity which can be remotely sensed. For example, at Earth,
VLF signals propagating in the surface-ionosphere waveguide allow location of lightning
events, and high frequency “sferic” (>3 MHz) emissions can be detected by Earth-orbiting
spacecraft (Herman et al. 1973). Global circuits may prevail on a number of solar system
bodies including Mars (Aplin 2006; Aplin et al. 2008).
The case for Martian atmospheric electricity has taken a unique path motivated by a
strong astrobiological context. After the initial tribo-electricity experiments of Eden and
Vonnegut (1973), Mills (1977) mixed grains in a low pressure CO2 gas and created impressive electrical activity (glow and spark discharges) in order to gain insights into the
environment that may have created ambiguities in the Viking lander biological experiments.
It was concluded that this active dust-created electrical environment could make Martian
dust storms an effective ‘scavenger’ of organic material. To date though, there have been no
direct measurements of the atmospheric electrical environment on Mars. Indirect evidence
of electrification however exists from the apparent adhesion of dust to the wheels of the Mars
Pathfinder and Sojourner rovers, suggested to be electrostatic in origin (Farrell et al. 1999;
Ferguson et al. 1999) and from laboratory measurements using Martian analogue materials
(e.g. Krauss et al. 2003; Aplin et al. 2012).
Given the lack of in situ measurements, there are three strong but circumstantial arguments that dust storms in the low pressure CO2 atmosphere of Mars may be electrical in
nature:
(1) Laboratory experiments (Farrell et al. 2015) of the breakdown from mixing sand grains
all display measurable electrical effects. These not only include the early glow-creating
laboratory experiments by Eden and Vonnegut (1973) and Mills (1977), but also a later
set of dust mixing studies in a low pressure CO2 gas performed by Krauss et al. (2003,
2006) and pre-glow and spark discharge studies by Farrell et al. (2015).
(2) Modelling studies (see Sect. 5) suggest that an electron avalanche and collisional plasma
can be created when dust storm E-fields stress the low pressure CO2 gas.
(3) Measurements in terrestrial desert environments reveal that dust devils and dust features
can generate large electric fields (Freier 1960; Crozier 1964; Farrell et al. 2004; Renno
et al. 2004; Delory et al. 2006; Jackson and Farrell 2006; Seran et al. 2013; Esposito
et al. 2015, 2016). If this same electrical generator process occurred in a lower pressure
atmosphere (like that at Mars), there would be the initiation of enhanced electron impact
ionization and atmospheric breakdown.
In terms of arguments against Martian atmospheric electricity, Ruf et al. (2009) reported
the possible remote-sensed detection of lightning emission at 8 GHz from Mars using a
terrestrial-based radio telescope. Given the detection, the electric dipole moment to account
for the discharge was found to be relatively large. However, subsequent work (Anderson
et al. 2012) could not confirm the initial ∼8 GHz observation. Also, a sensitive 4–5 MHz
Applications of Electrified Dust and Dust Devil Electrodynamics. . .
303
radio system aboard Mars Express carried out an unsuccessful search for lightning RF discharges, in close proximity to the planet (Gurnett et al. 2010).
These remote sensing observations do not confirm impulsive lightning events from
Mars. However, instead of intense impulsive cloud-to-ground discharges known at Earth,
laboratory experiments suggest the effects may be more subtle, with mixing of particles creating dark Townsend (gas ionization) discharges, low current glows, and local weaker but numerous spark discharges. In other words, the Martian system may be
continually discharging (or ‘leaking’) at low current levels into the low-pressure atmosphere, thereby avoiding the excessively large accumulation of charge which leads to
impulsive lightning events. In situ measurement systems offer the possibility to resolve
this apparent fundamental difference in our understanding of the terrestrial and Martian
systems.
2.2 Ion Balance in Atmospheres
On Mars, the main permanent ionization source is the bombardment by cosmic rays and the
cascade of ions and electrons generated when they collide with neutral atoms and molecules.
The maximum ion production rate by cosmic rays is at the surface because the atmospheric
thickness of Mars is not sufficient to allow cosmic rays to deposit all their energy in the
atmosphere, hence most of the energy reaches the surface. During day-time photo-ionization
by solar EUV and soft X-rays can also contribute with an efficiency increasing with altitude
while, at night, energetic electron precipitation from the magnetic tail of the planet may also
provide a supplementary source of ionization.
Since the initial work of Whitten et al. (1971), progresses in laboratory measurements of
ion-neutral reactions as well as radio-electric observations of the electron density profiles
(Patzold et al. 2005) have allowed several authors to build more detailed models of ionized species in the Martian lower atmosphere. According to Molina-Cuberos et al. (2002)
the most abundant ions below ∼50 km are hydronium ions H3 O+ (H2 O)n while O+
2 domi−
−
,
NO
nates above ∼70 km. The main negative ion species are water clusters of CO−
3
2 , NO3
−3
reaching a density of about 4500 ions cm at the surface. At 35 km the electron density reaches a peak of ∼102 e cm−3 and after a decrease with altitude, a secondary peak of
∼70 e cm−3 appears at ∼5 km. This secondary peak may significantly impact the negative
conductivity at ground level. Although not indicated in the paper, this secondary peak close
to the surface might be due to photo-electron production from UV impacts on the surface
(Grard 1995).
An atmosphere’s electrical conductivity depends on the local concentrations of ions and
electrons and ion species it contains. The total conductivity σ is given by
σ = e(μ+ n+ + μ− n− + μe ne )
(1)
where n+ , n− , and ne are the number concentrations per unit volume of positive ions, negative ions and free electrons respectively, and μ+ , μ− , and μe are their associated electrical
mobility (the drift speed per unit electric field). The concentrations of ions and electrons
are determined by the production rates by ionization and photoelectron emission, and their
losses by recombination and attachment to other particles. In the presence of a dust-laden
atmosphere, the concentration of ions and electrons (and therefore the electrical conductivity) will be reduced. The ratio of ions to electrons depends on the concentration of atmospheric electrophilic species; on Mars both negative CO−
2 ions and electrons are thought to
be present (e.g. Aplin 2006). The ion and electron concentrations are governed by a set of
304
R.G. Harrison et al.
rate equations for ion production and loss, written separately for positive ions, negative ions
and electrons by Tripathi et al. (2008) as
dn+
βj+ Nj ,
= q − αn+ n− − αe n+ ne − n+
dt
j
(2)
dn−
= q − αn+ n− − n−
βj− Nj − F n− ,
dt
j
(3)
dne
= qe + F n− − αe n+ ne − ne
βje Nj ,
dt
j
(4)
assuming no transport of ionization. In Eqs. (2) to (4), the volumetric ion and electron production rates are q and qe respectively, the ion-ion and ion-electron recombination coefficients are α and αe , and the negative ion detachment rate to yield electrons is F . The positive
ion-particle attachment coefficient to particles of concentration N carrying j elementary
charges is βj+ , that for negative ion attachment to particles βj− , and that for electron attachment to particles βje . Theoretical formulations for recombination and attachment coefficients
exist (e.g. Harrison and Tammet 2008; Tripathi et al. 2008), and the electron detachment rate
can be found from thermodynamic considerations.
2.3 Martian Conductivity Profiles and Relaxation Timescales
Solution of the ion balance equations with time gives the variations in ion and electron
concentrations, which, when the production and loss rates are equal, become steady-state
concentrations. If the ion production rates are known or can be calculated as a function of
height (for example from information about the cosmic ray energies and the atmospheric
density profile), and assumptions are made about the vertical profile of particles, the steady
state ion and electron concentrations can be found. In turn, this allows calculation of the
electrical conductivity, σ . An important parameter which can be determined from σ is the
charge relaxation timescale, τ , given by
τ=
ε0
σ
(5)
where ε0 is the permittivity of free space. This defines the timescale for the decay of charge
on a particle in a medium with conductivity σ , as if the particle were considered to be of a
finite capacitance and connected to a distant fixed potential through the medium. It therefore
provides a characteristic timescale for charge to be sustained on a particle. In the Earth’s
lower atmosphere, timescales of 10 s to 100 s are typical, which allows electrification to
develop with moderately active electrical processes such as those in terrestrial dust storms,
thunderstorms or volcanic plumes.
Figure 2(a) shows calculations from Tripathi et al. (2008) for the vertical profile of
electrical conductivity in the Martian atmosphere, under different assumptions of aerosol
loading. It is immediately apparent that the electrical conductivity is usually two orders
of magnitude greater than that for Earth. Only when the Martian atmosphere is assumed
to contain a dust storm does its lower atmosphere have an electrical conductivity comparable to that of Earth. (This has practical consequences for testing Martian atmospheric
electrical instrumentation on Earth, see Sect. 6). The associated relaxation times are shown
Applications of Electrified Dust and Dust Devil Electrodynamics. . .
305
Fig. 2 Profiles of (a) atmospheric electrical conductivity σ and (b) the associated charge relaxation
timescale ε0 /σ , for Mars (from Tripathi et al. 2008) and Earth (from Nicoll 2012). In the Martian case,
three scenarios of conductivity profile are considered, for a clear atmosphere, an aerosol-laden atmosphere,
and a dust storm of opacity 5. In the terrestrial case, the total conductivity is considered to be twice the
positive conductivity, under the quiet solar conditions observed by Gringel (1978)
in Fig. 2(b). This indicates that, for charge to be sustained on particles in the Martian atmosphere for timescales appreciably greater than one second, a dust storm conductivity environment would be required. The combination presented by a dust devil generating particle
electrification in a low conductivity environment may therefore offer unique circumstances
on Mars for active charge generation.
The effect of including the aerosols in the electrical model of the lower Martian atmosphere very substantially reduces the free electron concentration at low altitude. It decreases by more than two orders of magnitude compared to the simple model of ion-neutral
chemistry of Molina-Cuberos et al. (2002) to be less than ∼0.5 e cm−3 at ground with the
∼100 e cm−3 peak at 35 km lifted to ∼45 km, and also reduced by two orders of magnitude
at ∼1 e cm−3 .
Michael et al. (2007, 2008) comprehensively modelled the role of aerosols on the charged
particle distribution and conductivity in the lower atmosphere using aerosol properties
(Chassefière et al. 1995). These studies were performed for both night-time and daytime
conditions when dust particles are photo-ionized by low energy (∼6 eV) solar UV that reach
the Martian surface and produce a photo-electron population whose properties were taken
from the work of Grard (1995). A role for dust particles, in particular close to the surface, is
clearly emphasized by this study.
At night aerosols become charged by attachment of positive or negative ions. More
aerosols become charged close to the surface than at higher altitudes, and the dust content from Phobos 2 spacecraft data indicates an average decrease of conductivity by a factor
of two over the whole range of altitudes (Chassefière et al. 1995).
2.3.1 Tribo-Electric Charging Processes
The build-up of intense electric fields in the Martian atmosphere is expected to result from
charge exchange between particles colliding in a turbulent environment (see Fig. 3).
In order to understand the differences between Earth and Mars charging processes, the
combined system must be considered, where the E-field is generated by the competing tribo-
306
R.G. Harrison et al.
Fig. 3 An artist’s impression of
an electrified dust devil at Mars
(from Farrell et al. 2004)
electric charging currents and atmospheric dissipation currents. This can be expressed as
−ε0
dE
= −Jtribo + Jdiss = −ndust qdust vdust + σ E + n0 exp(αT d)eμe E
dt
(6)
where n0 is the nominal near-surface electron content, αT is Townsend’s first coefficient representing the number of electron impact ionizations per unit distance, d is the extent of the
driving E-field in the dust devil, and μe is the electron mobility. The last two terms in square
brackets represent the atmospheric dissipation currents, the first of which is the nominal dissipation current and the second an enhanced current associated with the electron avalanche
process for a gas under electrical stress. In the following subsections, these different terms
are compared for both Earth and Mars.
2.3.2 Tribo-Electric Current Generation, Jtribo
The mixing of dust grains creates electricity via tribo-electric processes with the grain polarity being a function of the particle mass and composition. Forward et al. (2009a, 2009b)
found, for particles of identical composition, that smaller grains tend to charge negatively
and larger grains tend to charge positively. As discussed more extensively in Sect. 3, they
argued that the mass-size dependency involves the collisional transfer of electrons trapped
in defect-created meta-states upon contact or rubbing, with the net exchange of electrons
statistically favouring the smaller grains due to their increase likelihood of defect region
‘rubbing’ with larger grains (see Fig. 3 in Forward et al. 2009a, 2009b). In the centre of a
convective feature, the grains become stratified in the upward flow of warm air. Small negative grains tend to lift upward and large positive grains remain lower in such features, giving
rise to the large-scale dipole moment and E-field like that measured within terrestrial dust
devils (Freier 1960; Crozier 1964; Farrell et al. 2004; Renno et al. 2004; Delory et al. 2006;
Jackson and Farrell 2006).
The dust load for a nominal terrestrial dust devil has been estimated to be up to
∼10−3 kg m−3 (Metzger et al. 2011). Assuming a dust devil primarily contains particles
of ∼5 µm diameter carrying −5 fC per particle with an upward flux of ∼2 × 109 m2 s−1
(Farrell et al. 2004; 2006b), an estimate of the associated upward tribo-charging current
source would be ∼10−5 A m−2 . Similar mass-dependent tribo-charging and mass stratification process are anticipated for Mars. However, given that dust devils on Mars can be larger,
larger tribo-charging current sources can be anticipated.
Applications of Electrified Dust and Dust Devil Electrodynamics. . .
307
Fig. 4 Effective electron
conductivity measured in a
representative Martian
environment, as a function of an
electric field maintained between
two plates with a varying
separation (from left to right in
the figure) of 60, 40, 20, and
8 mm. The different current
regimes (nominal, Townsend, and
spark) are delineated. (Reprinted
from Farrell et al. 2015, with
permission of Elsevier)
2.3.3 Atmospheric Dissipation Currents, Jdiss
The large contrast between terrestrial and Martian dust devil electrification lies in the nature
of the response by the atmosphere. At Earth, for typical disturbed electric field values, the
atmosphere does not initiate the electron avalanche process, and the dissipation currents are
defined primarily by conduction (i.e. proportional to σ E). For the terrestrial atmosphere
with conductivity ∼50 fS m−1 only a modest removal of current will occur in a dust devil
having E ∼ 100 kV m−1 , with the competing dissipation current at Jdiss ∼ 10−8 A m−2 . This
current is far less than the driving tribo-electric current, and is not the process limiting the
development of dust devil electric field. It is likely that it is the efficiency of the dust lifting
process itself which limits the electric fields generated.
In contrast, in a 6 hPa CO2 atmosphere like that at Mars, the nominal value of
the atmospheric conductivity is predicted to be about ∼100 times greater than that at
Earth (Farrell et al. 2015), with σ ∼ 10−12 S m−1 (see Fig. 2). However, even at modest
E-field values of 10 kV m−1 , electron avalanche processes will develop, with the exponential increase of electron currents directed along the driving E-field (Delory et al. 2006;
Jackson et al. 2010). In the low-pressure Martian atmosphere, the mean free path for electrons under the driving E-field is larger, thereby creating more energetic electron-CO2 impacts compared to the case of the denser Earth-like pressure. This greater number of electron
impact ionizations will lead to the exponential increase in electron content.
In essence, the stress of even a modest E-field at Mars initiates the electron avalanche
process, creating a substantial dissipation current that increases itself exponentially with
increasing driving E-field. As found in recent laboratory work (Farrell et al. 2015; see
Fig. 4), the dissipation current can become comparable to the driving Jtribo , creating a situation where the E-field becomes limited by the atmospheric dissipation (i.e. dE/dt ∼ 0 in
Eq. (6)). These dissipation currents can take the form of the Townsend dark discharge or a
glow or even a spark if breakdown is reached. This visual effect was reported by Eden and
Vonnegut (1973), Mills (1977) and Krauss et al. (2003). The atmosphere at Mars therefore
creates a substantial and significant competing dissipation current that act to deplete the
development of large separated charge centres.
2.3.4 Lightning
Given the possibility of strong dissipation currents that limit the growth of storm charge
centres, it seems that large impulsive ‘dissipation’ currents via lightning would be unlikely.
However, Jackson et al. (2008) and Kok and Renno (2009) found that the same dust grains
responsible for the charging current, Jtribo , can also absorb electrons from the dissipation
308
R.G. Harrison et al.
currents. This electron absorption would then reduce dissipation currents and thus could
lead to anomalously large charge centres in the convective feature, thus requiring impulsive
breakdown as a means to remediate excessive charge. Clearly, resolving which of these
competing electrical processes dominates requires in situ measurements.
2.3.5 Implications for Dust and Climate
The existence of a planetary electric field with very intense local enhancements and possibly
significant breakdown currents may have a significant influence on the physics and chemistry of the surface material. A detailed description of the expected impacts of electrification
on processes relevant to the climate of Mars is presented in Kok and Renno (2008b).
Dust is the most important source of local heating in the Martian atmosphere and is
known to control in large part the thermal structure of the troposphere. Many studies in the
past have focused on identifying and further characterizing the processes at work in levitating dust from the surface of Mars. Well-established work (Bagnold 1962) indicated that
direct wind lifting on Mars typically occurs at wind speeds exceeding 25 ms−1 near the surface (Greeley et al. 1992), a condition that is difficult to reach, according to current general
circulation models. Instead, several other mechanisms have been proposed and explored in
laboratories that provide the necessary lifting strength without requiring winds as strong as
the standard theory predicts. The proposed mechanisms comprise sand blasting process and
levitation induced by electrostatic forces. Considering the potential electrical force exerted
on micron size dust particles by electric fields predicted inside dust devils or dust storms,
electrification may potentially contribute to the global transport of dust at a level comparable
as the drag force due to the wind (Berthelier et al. 2000). Studies have explored the more
specific impacts of electrification on the lifting processes of dust from the surface of Mars,
unveiling the unexpected role of electrification in the trajectories of saltating particles, an
effect first suggested by Schmidt et al. (1998) and Zheng et al. (2003). A major effect concerns the limitation imposed by electric forces on the height of the saltation layer. Contrary
to the standard theory, the observed height does not increase with increasing wind speed, as
it simultaneously leads to development of a downward-pointing electric force exerted over
the saltating particles (Kok and Renno 2008a, 2008b). A related effect is the emergence of a
positive lifting feedback loop as a charge gradient exists between the positively charged particles lying at the surface and the negatively charged saltating particles. This phenomenon
leads to a significant reduction of the shear velocity and thus aids hydrodynamic forces in
lifting particles from the surface (Kok and Renno 2008a, 2008b). Direct lifting by the electric force may also occur in parallel, and thus establishes a combination of effects that can
contribute significantly to dust lifting and potentially ultimately atmospheric dust loading
on Mars. This remains something which global modelling of the Martian atmosphere has
yet to consider quantitatively.
3 Contact and Tribo-Electrification of Environmental Dust and Sand
The importance of particle–particle interactions in the electrification of dust clouds and
structures such as dust devils highlights the understanding needed in contact and triboelectrification. This section now reviews understanding of the physical process(es) of contact
and tribo-electrification. The discussion is based principally on laboratory studies, but also
presents contemporary models and theories considering mineral dust and sand, because of
their relevance to dust devils.
Applications of Electrified Dust and Dust Devil Electrodynamics. . .
309
3.1 Definitions, History and Background
It is worth discussing the terminology regarding particulate electrification, since even this
can be a source of confusion. Contact electrification was originally used as a broad term
for electrostatic charge transfer resulting from contact, including contact modes such as
detachment, sliding, rolling, impact, etc. Tribo-electrification was later used to describe this
phenomenon resulting from ‘rubbing’. Currently both contact- and tribo-electrification are
widely and interchangeably used in the literature. Conventionally (though not formally)
contact electrification has since become associated with ‘contact and macroscopic separation
leading to charge transfer’ (e.g. McCarty and Whiteside 2008), whereas tribo-electrification
is typically used where emphasis is placed upon (empirical) dependence on composition
(i.e. the triboelectric series) and impact velocity, and it is also interchangeably used with
‘frictional electrification’.
Despite the importance of contact electrification and the amount of research in this field,
there is little agreement upon the core mechanism of charge transfer. Generally it is now
thought that there are several competing process occurring, rather than one single mechanism being responsible, and that this is dependent upon the nature of the interacting surfaces
and especially the surface chemistry. On the most fundamental level one may consider contact electrification to occur via charge exchange with the transfer of electrons and/or ions.
There are several models involving each of these processes which have been successfully
applied in respective experimental cases. A complete picture for all materials in all cases
still eludes researchers and has been a source of active debate for decades (McCarty and
Whiteside 2008; Harper 1998; Lowell and Rose-Innes 1980).
3.2 Experimental Studies of the Phenomenon
Numerous contact or tribo-electrification studies have used macroscopic contacting surfaces, often involving an insulator and a metal electrode (McCarty and Whiteside 2008).
Here, however, focus will be placed on granular materials (sand/dust). Typical methods
for quantifying electrification of coarse granular material (sand) involve removal and direct measurement of the electric charge, often using a combination of a Faraday cup and
electrometer. The specifics of transport and collection vary, including; cascades, fountains, fluidized beds, blow-away experiments, aerosolizers or single particle impact studies
(Matsusaka et al. 2010; Sickafoose et al. 2001; Poppe et al. 2000; Kok and Lacks 2009;
Merrison et al. 2012).
3.2.1 Surface Charge Density
Techniques involving direct electric charge extraction are ineffective when dealing with dust
sized (micrometer) particulates since they are typically well suspended and difficult to extract (e.g. into a Faraday cup). Suspended particles may also only be slightly charged, and,
in contact with a surface, they may not be electrically conductive enough to allow rapid
charge extraction. For fine suspended dust other techniques have typically been used to
quantify particulate electrification. They rely upon the application of an electric field and
the drift of electrified dust grains. The field induced drift velocity can be determined using
optical or laser systems to study particle trajectories (Kunkel 1950; Merrison et al. 2012;
Mazumder et al. 1991) or by extracting dust onto a surface (Merrison et al. 2012). The fieldinduced drift velocity will be proportional to the electrical force and therefore the electric
charge and polarity of the dust grains.
310
R.G. Harrison et al.
In determining order of magnitude limits for expected values of contact electrification
there has been considerable success using the ratio of charge to surface area. In a variety of
experimental techniques and over a broad range of grain sizes (µm–mm) the electrification
has been seen to be of order 0.1 mC m−2 (Poppe et al. 2000; Lowell 1986; Merrison et al.
2004; Nieh and Nguyen 1988). In many granular electrification studies the electrified surface
area is taken as the (total) surface area of the grain; since the contact area is not known, this
possibly explains the lower level of electrification seen in such cases e.g. ∼10−3 mC m−2
(Sickafoose et al. 2001; Gross et al. 2001).
3.2.2 Specific Charge
Many workers, however, still express electrification in terms of the specific charge (i.e. the
charge per unit mass, Q/m), due often to experimental convenience. There is less success
in establishing expected/limiting values. In aeolian sand transport studies values have been
quoted of order 60 µC kg−1 , (Schmidt et al. 1998; Zheng et al. 2003, 2004; Qu et al. 2004; Bo
et al. 2014). In gas-solid pipe flows continuous sand transport is performed while measuring
the current to metal pipe, here values of order 3 mC kg−1 have been measured (Matsusaka
et al. 2010). In metal oxides charges after blow-off exceeding ±200 mC kg−1 have been seen
(Oguchi and Tamatani 1993) (see the section on wind tunnel experiments).
3.2.3 Practical Experimental Limitations
Another common cause for observing reduced electrification in charge collection systems is
charge leakage via surface water, which is related to ambient humidity (Nieh and Nguyen
1988). In dust collection experiments electrostatic aggregation can also de-electrify suspended dust, especially when in high concentrations and after long suspension times (Merrison et al. 2012). In relatively recent studies using Atomic Force Microscopy (AFM) techniques and well prepared surfaces, significantly higher values were reported, although in
these no account was made for lateral spreading of charge (Horn et al. 1993).
In many cases electrical breakdown can limit the charge on grains during separation, e.g.
dielectric breakdown in terrestrial surface air is at an electric field of around 3 × 105 V m−1 .
(The lower surface pressure on Mars yields a smaller breakdown voltage there (Laughton
and Warne 2004)). This could provide a quantitative explanation for the experimentally observed upper surface charge concentration Q and has been supported experimentally in detailed investigation using polymer micro-spheres (50–500 µm radius r) where Q/r 2 is seen
to be constant, of order 0.1 mC m−2 . Electrical breakdown has also been directly observed
in some experiments (McCarty et al. 2007; Harper 1998; Matsuyama and Yamamoto 1997;
Horn et al. 1993; Matsusaka et al. 2010). The observed charging will therefore in many cases
be limited (by dielectric breakdown) to values less than 0.1 mC m−2 .
Even lacking a detailed physical understanding of the electrification process, there is
therefore general agreement on the order of magnitude of the electrification and this is sufficient in most cases to quantify the effect of electrification on, for example, entrainment and
transport of sand and dust.
3.3 Particulate Size Dependence
It is widely accepted that particulates of the same composition and of differing size will
show a tendency for larger particulates to electrify positively and smaller ones to electrify negatively, upon contact and separation. In experimental studies the electric fields
Applications of Electrified Dust and Dust Devil Electrodynamics. . .
311
generated by terrestrial dust devils support the idea that suspended dust becomes electrified negatively with respect to the sand/sand-bed (Schmidt et al. 1998; Zheng et al. 2003;
Qu et al. 2004). This size dependence has been best demonstrated in laboratory experiments involving sand cascading (Lacks et al. 2008; Kok and Lacks 2009; Lacks and
Levandovsky 2007, Forward et al. 2009a, 2009b; Bilici et al. 2014). However, not all experimental studies reproduce this behaviour (Trigwell et al. 2003; Sowinski et al. 2010;
Kunkel 1950). It should also be noted that in many of these electrification studies multiple
particle interactions are involved, including particle-wall interactions (Aplin et al. 2012),
and material purity (surface composition) is not well controlled. A more complex charge
exchange behaviour can therefore often not be ruled out.
In laboratory experiments of sand transport in the absence of dust, it has been seen
that sand may electrify either positively or negatively depending on the size distribution
within the sand bed (Bo et al. 2014; Zheng et al. 2003, Kok and Renno 2008a, 2008b).
Laboratory experiments have also shown that dust re-suspension in the absence of sand
shows little net (size dependent) dust electrification (Merrison et al. 2012), however the
role of dust aggregates (acting like large sand-sized particles) could possibly complicate this behaviour. This complexity makes describing/predicting contact electrification
for example within a Martian dust devil problematic. Modelling has been unsuccessful
in satisfactorily explaining this size dependence in contact electrification, although several promising models are being pursued. For example, one model is based upon electron
transfer through so-called high energy electron surface states (Lacks and Sankaran 2011;
Apodaca et al. 2010; Bo et al. 2014). It would be extremely useful here to experimentally
identify the precise charge/polarity dependence of contact electrification with grain size: this
has yet to be done.
3.4 Material and Humidity Dependence
Attempts to apply a single model (i.e. involving either electron or ion transfer) for all cases of
contact electrification have failed. Although the use of various so called tribo-electric series
have been empirically useful (and widely published) there are many cases where electrification does not follow such tribo-series, and examples exist of circular tribo-series (Harper
1998) wherein a repeating sequence of materials paradoxically generates increasing values
on the tribo-series. Also the tribo-series cannot explain electrification of like materials with
different sized particles (Lacks and Sankaran 2011). This suggests that there is not one single
electrification process occurring (McCarty and Whiteside 2008).
It is however apparent that material properties are crucial in applying models of
contact electrification. For contacting metals, electron transfer models have been extremely successful in reproducing observed electrification based upon differences in electron work function (e.g. Matsusaka et al. 2010). Similarly for insulators containing mobile ions, modelling involving ion exchange has provided a consistent picture and even
allowed control in contact electrification (Diaz and Felix-Navarro 2004; Law et al. 1995;
Mizes et al. 1990). Contact electrification in (non-ionic) insulating materials though, remains problematic. Currently most researchers assume that contacting (non-ionic) insulators
involve electron transfer (Lowell and Rose-Innes 1980; Grzybowski et al. 2005).
Experimentally, a typical problem is purity, especially the presence of surface contamination. Since the number concentration of charges at the surface involved in electrification is
low compared to the concentration of molecules/atoms, even low concentrations of impurities (of order ppm) can in principle dominate the electrification process. This makes studies
of contact electrification technically challenging.
312
R.G. Harrison et al.
Under ambient terrestrial conditions layers of water molecules are invariably present at
surfaces and can play role in electrification by enabling ion or electron transport. Since
humidity and surface ions such as H+ and OH− , are ubiquitous, models have been proposed to describe electrification of insulating materials (not containing ions) in which for
example OH− accumulate on surfaces with water layers. In experiments using polymers
and bound/unbound ions (Diaz and Felix-Navarro 2004) contact electrification was seen to
be eliminated at 0 % RH (relative humidity) and rise to a maximum at 30–40 % RH, supporting this ion exchange model. Again, however, there are experimental cases in which
contact electrification has been observed despite the absence of surface water (low humidity/vacuum) (Harper 1998; Lowell and Rose-Innes 1980) and also with well-prepared dry
surfaces (e.g. Horn et al. 1993; Gady et al. 1998). Hence, despite advances made (Diaz and
Felix-Navarro 2004; Law et al. 1995; Mizes et al. 1990) involving ion transfer (and the use
of charge control agents), this model still cannot be applied in all cases of contact electrification.
In the case of real planetary silicate minerals, ions will likely be present at the surfaces
(e.g. OH− , alkali metals, halide ions, etc.), and, terrestrially, water vapour (and therefore H+ ,
OH− ) will also be present and allow ion-transfer models to be applied. Generally, however,
most current modelling of aeolian (contact/tribo-) electrification involves electron transfer
(only) (Lacks and Sankaran 2011).
3.5 Models and Theories for Dust and Sand Electrification
For mineral dust and sand it is likely that both electron and/or ion transfer mechanisms
may be relevant. Here material transfer is considered to be a form of ion transfer (Tanoue
et al. 1999). For non-ionic insulators (e.g. polymers and possibly including mineral sand
or dust) most researchers attempt to employ electron exchange models (Lowell and RoseInnes 1980). However contact electrification does not correlate with surface or bulk electron
properties such as dielectric constant, atomic properties, ionization energy, electron affinity or electro-negativity (Wiles et al. 2003). Despite this, some success has been achieved
with electron transfer models involving electron donor/acceptors or so called high energy electron states in insulator-insulator contact, (also known as a molecular ion state
model) (Lowell and Rose-Innes 1980; Bailey 2001; Duke and Fabish 1978). Similar electron transfer models involving so-called high energy surface electron states (likely due to
impurities/contamination) are used in describing size dependent electrification (Lacks and
Sankaran 2011; Bo et al. 2014).
Specific electron transfer models have been developed for electrification in metalinsulator contact, for example; involving an effective potential difference (Davies 1969) or
performing quantum chemical (electron state) calculations (Shirakawa et al. 2008). Other
models have also studied multiple surface impacts; one involved a capacitance charging model including charge relaxation (electrical discharge) (Matsusyama and Yamamoto
2006).
Despite these advances in electron transfer models, similar success and progress has also
been achieved using ion exchange models specifically involving (proton exchange) which
not only works well for ionic or ion doped materials, but also where water layers may be
involved relating to chemical properties such as pH and/or zeta potential (Diaz and FelixNavarro 2004; Law et al. 1995; Mizes et al. 1990). It has even been argued that there is ‘never
electron transfer with insulators’ (Harper 1998), however this has not been demonstrated
experimentally and is not widely accepted.
Based on an ion transfer (entropy driven) model, an order of magnitude calculation can
be performed assuming a concentration of mobile ion groups (e.g. H+ , OH− , metal ions,
Applications of Electrified Dust and Dust Devil Electrodynamics. . .
313
halide ions, etc.) of, for example, around 0.1 nm2 generating around 3000 e µm−2 , this is
in reasonable agreement with the expected upper limit (Diaz and Felix-Navarro 2004). This
still corresponds to a surface concentration of <1 %. Experimental tests have yet to be
performed of the temperature dependence in electrification which could verify the entropy
driven model.
3.6 Wind Tunnel Experiments
On the scale of laboratory wind tunnels it has proved difficult to generate and therefore measure electric fields produced by sand/dust transport. It has been more practical to measure
and determine the electrification of single grains (and distributions) and then apply modelling to predict electric field generation.
The electrification of saltating sand grains has been studied using wind tunnels in an
extensive series of experiments, which also combined modelling with theory, involving determining the effects on grain trajectories of applied electric fields. In these experiments
sand grain electrification was measured in a conventional manner using an Faraday cup type
electrode to collect the charge of impacting grains and an electrometer to quantify the discharging current (Bo et al. 2014).
Grain electrification distribution was typically seen to be bipolar and broadly in the range
from −300 to +600 µC kg−1 . In the earliest wind tunnel studies the average electrification
was seen to be negative (Zheng et al. 2003; Qu et al. 2004; Zhang et al. 2004), however
in a more recent study the net charge is seen to be positive and also a higher degree of
electrification (agreeing with field experiments of Schmidt et al. 1998), Bo et al. (2014).
This discrepancy was interpreted as due to interaction (charge exchange) with finer grains
(non-saltating) within the wind tunnel sand bed. Notably a (broad) grain size distribution was
used, which was not representative of those seen in the field. From this work it seems that
there is a dependence of the measured electrification on the grain size of both the saltating
and non-saltating grains, as well as on wind speed and height (increasing strongly with
height).
Based on this sand electrification work electric fields have been predicted to be generally
less than 50 kV m−1 and typically a few kV m−1 (Bo et al. 2013a, 2013b), this is below
values at which significant effects would be expected on saltation and below (upper limit
for) value measured in the field (Farrell et al. 2004; Schmidt et al. 1998; Renno et al. 2004).
It should be noted however, that the models applied did not include dust entrainment, i.e.
fine, suspended fraction which may have enhanced the electric field.
For fine (suspended) dust grains the techniques based on measuring their electrical discharge current to an electrode is difficult and ineffective due to: the small charge degree, the
particles not being ballistic (having low impact velocity) and not liberating electrical charge
rapidly. An alternative technique has been applied in a series of wind tunnel studies of dust
electrification. Here an electric field is applied in order to drift dust grains out of suspension
and collect them on electrodes. The quantity of dust collected can then be quantified either
optically or by mass. By determining the amount of collected dust as a function of applied
electric field, the average electrification per grain can be established (Merrison et al. 2012;
Merrison et al. 2004).
It was found that the degree of electrification of re-suspended dust is typically of the order of 103 to 104 e/grain for µm scale particles (or around 0.1 mC m−2 ). In the absence of
sand the net electrification of the dust was seen to be close to zero (therefore not expected to
generate an electric field). However, dust aggregates (agglomerates) are also seen to become
electrified and may show the electrostatic behaviour of real dust clouds, even in the absence
314
R.G. Harrison et al.
of larger sand particles. Specifically in wind tunnel studies, it was found that electrification is affected by competing processes of aggregation (electrostatic self-assembly), causing
reduced charge and aggregate dispersion causing electrification.
Xie and Han (2012) performed experiments in a wind tunnel to investigate the effect
of relative humidity on aeolian E-field. The working section of the chamber was 20 m ×
1.3 m × 1.3 m. A sand bed of 8 m × 12 cm was deposited in the chamber. Sand particles
were collected in the Badain Jaran desert in China. The electric field was measured with
a KDY-IV field mill. At a given wind speed, they measured the aeolian electric fields at
different air relative humidity (RH) and sand relative moisture conditions. They observed
the electric field linearly increased with increasing relative humidity up to a critical value
and then exponentially decrease. This critical RH value was observed to increase with the
wind speed. For a wind speed of 14 m s−1 the critical RH was 32.7 %.
Wind tunnels have been more extensively used to investigate the influence of applied
electric fields upon sand transport (saltation). Several groups have reported reduction in the
threshold for saltation and enhanced transport rates of sand due to applied electric fields
(Kok and Renno 2006, Rasmussen et al. 2009), typically at electric fields above around
100 kV m−1 . This work has mainly been based upon the assumption of a conductive sand
bed (surface). However, recently it has been experimentally demonstrated that, for the case
in which the sand bed is insulating, the application of an electric field in fact significantly
increases the threshold shear stress for saltation (Holstein-Rathlou 2012). This is due to
dielectric attraction, entirely analogous in its operation to diamagnetic chain formation in
the presence of an applied magnetic field.
3.7 Charging and Radioactivity in Volcanic Plumes
Theoretical studies of triboelectric charging of single particle systems demonstrate that the
particle size distribution determines the magnitude and size of charging (Lacks and Levandovsky 2007). Laboratory studies of electrical charging of volcanic ash systems confirm that
the particle size distribution plays an important role in the electrical charging of terrestrial
volcanic plumes (Houghton et al. 2013; Cimarelli et al. 2014) and indicate that the composition of the particulate matter making up the plume also affects electrical charging. The factors that affect triboelectric charging in terrestrial volcanic plumes, namely the particle size
distribution, particle composition and relative humidity, are expected to have similar effects
in the aeolian Martian environment. While recent work has contributed significantly to the
understanding of triboelectric charging in terrestrial volcanic plumes, further investigation
is required to fully understand this process and the effects of the particle size distribution,
material composition and ambient conditions on charging.
Alpha and beta particles are emitted following the radioactive decay of the uranium, thorium and potassium radioisotopes, leaving a residual charge on the particle from which they
have been emitted (Clement and Harrison 1992). Studies of ash collected from the 2010
Eyjafjallajökull and 2011 Grimsvötn eruptions demonstrate that radioactive decay is only
seen in larger particles (Aplin et al. 2014). This size dependence means that for terrestrial
volcanic plumes, charging associated with radioactive decay is only likely to occur near the
vent itself, as the larger particles containing radioisotopes are rapidly removed by gravitational settling. The composition and size of particles in Martian dust storms are expected to
determine whether or not radioactive decay contributes to charging in Martian dust storms;
larger particles are expected to be lost first from the dust cloud.
During volcanic eruptions, magma fragmentation can contribute to plume charging near
the vent (Gilbert et al. 1991; James et al. 2008). Charged species are ejected from cracked
Applications of Electrified Dust and Dust Devil Electrodynamics. . .
315
surfaces following material fracture, and the loss of these charged species results in the
formation of charged fragments. However, particle fragmentation is not expected to occur
in Martian dust storms.
Clearly there is much more research needed in this area before we reach an understanding
of these disparate electrification processes. For wind tunnel investigations it would be informative to perform electrification studies which involve both sand and dust grains (i.e. using
an extremely broad size distribution from 1 µm–1 mm) in order to gain an understanding of
the interplay of electrification and field generation processes.
4 Atmospheric Observations of Dust Electrification
Atmospheric electrical parameters can be sensed in a variety of ways. The most commonly
measured quantity is the vertical electric field,1 but the air conductivity, charge density and
vertical conduction current have also been routinely observed at some measurement sites.
Many of these techniques have developed rapidly as electronic systems have improved, as
exceptionally good performance—very often ultra-low leakage current—is required from
the signal processing electronics across a wide range of conditions. A widely used fastresponse instrument for electric field measurement is the field mill, which operates by the
alternate exposure and screening of a sensing electrode using a rotating shutter followed by
phase-sensitive detection. A field mill can be orientated with its sensing surface upwards or
downwards for the vertical component of the field, although spherical geometries have also
been used in dust devil studies, to obtain the horizontal field component (Ravichandran and
Kamra 1999). Vertical electric fields can also be measured using passive (long horizontal
wire) antennas or other electrodes which come into electrical equilibrium with the potential
of the surrounding air. These can have a slow response because of the low conductivity of
atmospheric air: a radioactive coating can increase the air conductivity and the associated
time response. The position of a sensing instrument mounted above the surface acts to distort
the local electric field, and calibration to the electric field over a flat undisturbed surface is
needed. The negligible distortion of a long wire antenna provides one possible reference
technique to calibrate atmospheric electric field measurements, which has been used for
over a century (Harrison 1997, 2013; Bennett and Harrison 2006).
Another approach used in high field conditions is to measure the so-called point discharge current, which is the current flowing into a vertically facing sharpened vertical point
under the influence of the atmospheric electric field (Marlton et al. 2013). Techniques for
measurement of air conductivity and conduction current density are not discussed further
here, but summaries of the associated methods required are given in Harrison (2004), Aplin
and Harrison (2001) and Bennett and Harrison (2008).
4.1 Surface Measurements
As introduced in Sect. 1, Freier (1960) was one of the first researchers to measure electric
fields in dust devils. It was found that the electric field was upward pointing in dust devils.
Freier’s findings were later confirmed by Crozier (1964, 1970). These pioneering measurements of dust devil electric fields are consistent with results of more recent measurements
1 Note than a convention in fair weather atmospheric electricity is to report the Potential Gradient (PG) rather
than the vertical electric field Ez . In fair weather, the potential increases positively with height, and the PG
is typically +120 V m−1 near to the surface. (Ez in these circumstances would be −120 V m−1 , downward
pointing).
316
R.G. Harrison et al.
Fig. 5 ACI Project site at the
Owens Lake playa with electric
field sensors mounted at 0.5, 1.5,
2, and 3 m above the surface. The
results reported here are from
measurements of the electric field
at 3 m above the ground. The
PM2.5 aerosol concentration is
measured at approximately 3.5 m
above the ground, and saltation is
measured at 7 cm above the
ground
also indicating upward pointing electric fields in dust devils (e.g., Farrell et al. 2004; Renno
et al. 2004).
In dust storms the near surface electric field is stronger and in the opposite direction to
the downward pointing fair weather electric field (e.g., Latham 1964). This suggests that
dust particles become negatively charged after colliding with the larger sand particles and
the surface (e.g., Schonland 1953; Latham 1964; Kok and Renno 2008a, 2008b; Renno
et al. 2003). Schmidt et al. (1998) reported that the electric field is also upward pointing in
saltation alone (that is, without dust lifting). Indeed, Schmidt et al. (1998) reported upward
pointing electric fields of 160 kV m−1 at a few cm above the surface during saltation in
sand dunes. Electric fields of this magnitude could have significant effects on saltation and
therefore dust lifting.
4.1.1 Recent Measurements of Electric Fields in Dust-Devils
Except for the measurements of Schmidt et al. (1998), all previous measurements summarized above were made using grounded instruments. Since grounded instruments typically
measure the difference in potential between the ground and their sensing elements, rather
than the electric field at the sensor height, measurements with these grounded instruments
must be interpreted carefully. The presence of grounded instruments distorts the electric
fields being measured and the instrument configuration must be calibrated on site (e.g.
Sullivan 2013). In addition, the impact of charged dust or sand particles introduces noise
on grounded instruments by the effect of electric currents between the instruments and the
ground (Sullivan 2013).
Schmidt et al. (1998) used a cylindrical field mill isolated from the ground to measure the
electric fields in the saltation layer directly (Johnston et al. 1986). Since charge transfer by
the collision of charged particles with isolated field mills vanishes as soon as the field mills
reach the same potential as the charged particles, the errors produced by these collisions also
vanish. Renno et al. (2008) developed an isolated cylindrical field mill for making measurements of the electric fields in saltation, dust storms and dust devils. This new instrument is
capable of measuring the electric fields accurately even when the instrument is subject to
the impact of charged particles (Renno and Rogacki 2013).
Here we report results of measurements with this instrument in the Aerosols-Climate
Interaction (ACI) project site at the Owens Lake salty playa, in California (Halleaux and
Renno 2014). Figure 5 shows four electric field sensors installed at the ACI project site at
the Owens Lake playa. Figure 6 indicates that, in mid-March, a saltation event triggered an
Applications of Electrified Dust and Dust Devil Electrodynamics. . .
317
Fig. 6 PM2.5 aerosol concentration (upper panel), near surface electric field strength (middle panel), as
well as electric field direction and indication of saltation events in blue (lower panel) at ACI project site in
the Owens Lake in 2014. The measurements of the electric field strength and direction presented in Fig. 6
were made at 10 Hz. (Upward and downward pointing electric fields are defined as electric fields within
approximately 30◦ of the vertical. Saltation is measured with a Sensit sensor placed at about 7 cm above the
surface e.g. Halleaux and Renno (2014). Saltation events are defined as those time periods in when the impact
rate at the sensor is greater than 10−4 particles/second. Saltation intensity refers to the normalized number of
impacts on the Sensit sensor per unit time.)
order-of-magnitude increase in PM2.5 aerosol concentration. The hourly mean value of the
electric field exceeded 10 kV m−1 during this event. This value is two orders of magnitude
larger than the fair weather electric field.
Figure 6(c) indicates that, in general, the electric field is downward pointing, but that
when saltation and therefore dust lifting occurs the electric field points upward. The few
exceptions occur when saltation is weak and during stormy weather when charged clouds
are also present.
Figure 7 shows the results of measurements of PM2.5 aerosol concentration, electric field
and saltation during a day of intense saltation and dust lifting. The measurements indicate
that, in general, the electric field changes its polarity to upwards when saltation occurs. The
exceptions to this occurred during periods when saltation was weak or when charged clouds
were also present above the site.
4.1.2 Measurements of Electric Fields in Dust Storms
Measurements of electric fields during dust storms are sparse in the literature but show both
upward (Rudge 1913; Harris 1967; Stow 1969; Zhang et al. 2004) and downward (Demon
et al. 1953) electric field directions.
318
R.G. Harrison et al.
Fig. 7 PM2.5 aerosol concentration (upper panel), near surface electric field strength (middle panel), as well
as electric field direction and saltation events (lower panel) at the ACI project site in the Owens Lake on a
day of active saltation and dust lifting events (May 10, 2014)
Kamra (1972) performed measurements of the potential gradient between the surface and
at 1 m above, in several US locations using a 500 µCi Polonium probe as a potential sensing
probe. (The use of radioactivity enhances the time response of the sensor, without adding the
mechanical complexity of a field mill.) The potential of the radioactive probe was measured
by a Keithley model 600B electrometer, a laboratory grade device. During dust storms, this
instrument recorded both upward and downward E-fields depending on the nature of soil
in the measurement site. In particular, Kamra observed that dust storms dominated by clay
minerals produced only negative potential gradients, while, when silica dominated in the
dust storm, both polarities of potential gradients were produced.
During measurements in the Sahel with a Mission Instruments field mill, Williams
et al. (2009) observed E-fields directed both downward and upward. Very recently, Esposito et al. (2016) carried out several measurements in the West Sahara (Merzouga region, Morocco) during the dust storm season in 2013 and 2014. They performed simultaneous measurements of atmospheric electric field, using a Campbell CS110 field mill,
mounted at 2 m and directed toward the surface, atmospheric dust abundance and size distribution (in the range 0.265–34 µm), sand saltation rate (with the Sensit sensors), wind
speed and direction at six different heights, air and soil temperature, pressure, relative humidity, soil moisture content and solar irradiance. They observed predominantly negative
electric fields, i.e. with the same sign of fair weather E-field, with some sporadic and short
inversions in polarity (Fig. 8). This was observed both during dust storms and dust devils
(Fig. 9).
Applications of Electrified Dust and Dust Devil Electrodynamics. . .
319
Fig. 8 Typical dust storm
measured in the Sahara desert in
July 2013 (Esposito et al. 2016)
The relation between humidity and particle charging mentioned in Sect. 3.4 was observed
by Esposito et al. (2016) during dust storms monitored in the Sahara desert. By a statistical analysis of the data, they found a linear relationship between the concentration of dust
entered in suspension and atmospheric electric field. The slope of this linear trend was observed to increase when relative humidity reached a critical value. This value seemed to
depend also on the content of soil moisture and was around 20 % in the 2013 campaign and
30 % in 2014 where the soil was drier (Fig. 10).
The increase in the slope indicates that, for each value of dust concentration a lower
E-field is observed for relative humidity larger than the critical value. Taken together, all
these observations suggest there are complex mechanisms underlying natural dust electrification, and that different effects should be expected in different regions on Earth and Mars.
320
Fig. 9 Properties of a dust devil
observed in the Sahara desert in
July 2014. Upper panel: rate of
impact of particles. Second panel:
dust particle concentration. Third
panel: wind direction. Fourth
panel: wind speed. Fifth panel:
electric field. Bottom panel:
atmospheric pressure
R.G. Harrison et al.
Applications of Electrified Dust and Dust Devil Electrodynamics. . .
321
Fig. 10 Relationship between
atmospheric electric field and
concentration of dust emitted
during several dust storms
observed in the West Sahara in
2013 (upper) and 2014 (lower)
field campaigns by Esposito et al.
(2016). The slope of the linear
trend increases when relative
humidity exceeds a critical value.
This value depends also on the
moisture content of the soil
emitting dust grains
4.2 Measurements of Particle Charging Aloft
Although measurements of dust electrification in different circumstances have been made at
and near Earth’s surface, measurements made aloft are rare in comparison, due to the difficulties associated with measurement platforms and lack of suitable instrumentation. Of the
handful of airborne studies which do exist, all have focused on passive dust layers, which
are typically isolated from the surface, travelling many hundreds of km from their source
regions. As yet, no attempt has been made to measure electrification in dust storms or dust
devils above the surface, which is certainly something that should be attempted in future to
establish the in situ charge structure. Of the measurements that do exist, Gringel and Mühleisen (1978) were amongst the first to publish measurements of electrical conductivity, σ ,
322
R.G. Harrison et al.
Fig. 11 Vertical profile of total
aerosol particle concentration
(grey) and magnitude of space
charge density (black) measured
on a specially instrumented
radiosonde flight through
Saharan dust layers from Sal,
Cape Verde. Adapted from Nicoll
et al. (2011) “© IOP Publishing
Ltd. CC BY-NC-SA”
measured by a Gerdien condenser, made from a free balloon platform through a 2 km thick
layer of Saharan dust off the west coast of Africa. Their data show a clear decrease in the
conductivity (caused by the attachment of ions to the dust particles) within the layer, dropping by a factor of two compared to its clean air value at the same altitude. These findings
are to be expected from a passive dust layer where no active charge generation or charge
separation occurs. This is in accordance with Ohm’s law, which governs charge conservation in the fair weather atmosphere, and states that the vertical conduction current, Jc , is
equal to the product of the electrical conductivity and the local atmospheric electric field
Ez , (Jc = σ Ez ). In the case of no active charge generation, Jc through the dust layer should
be continuous, therefore any decrease in conductivity (caused by ion-aerosol attachment)
will lead to an increase in the local electric field within the dust layer. This hypothesis is
supported by measurements through aerosol layers (e.g. Markson 2007), which are similar
to dust layers in terms of the presence of large particles, and show clear increases in Ez
inside the layer, as a result of the conductivity decrease.
Following on from the work of Gringel and Mühleisen (1978), Nicoll et al. (2011) made
a series of instrumented balloon flights from the Cape Verde Islands, near the west coast
of Africa, which frequently experiences elevated layers of dust transported from the Sahara
in the so-called “Saharan Air Layer”, which can reach altitudes of 5–7 km. Simultaneous
measurements of dust particle concentration and space charge, ρ, were obtained by two
small sensors attached to standard balloon-carried meteorological radiosondes. The space
charge sensor measures the net charge (i.e. the difference between positive and negative
charge) carried by the dust particles and comprises a small spherical electrode connected to
a sensitive electrometer circuit, measuring the voltage on the electrode.
Figure 11 shows measurements from one of the instrumented balloon flights from Sal,
Cape Verde, which shows two distinct layers of dust particles at 2 km and 4 km (with particle diameters up to 2.6 µm). Measurements from the charge sensor show weak charge on the
dust particles in both layers (up to 10 pC m−3 ). Particle–particle charging (i.e. triboelectric
charging) of dust grains is likely to have originated near the surface during the initial lofting
of the dust, but this charge will decay quickly (within several minutes) due to recombination
with ions. This therefore suggests the presence of a continual charging (albeit weak) mechanism which is active in elevated dust layers, likely to be related to tribo-electrification of the
dust particles (Houghton et al. 2013). Such an effect has also been observed in elevated layers of volcanic ash from the eruption of Eyjafjallajökull in Iceland in 2010 (Harrison et al.
2010). The effect of such charge on long range transport of dust and ash is as yet unknown,
but may have implications for modification of deposition speed as well as particle–particle
agglomeration rates and particle wet removal by droplets.
Applications of Electrified Dust and Dust Devil Electrodynamics. . .
323
4.3 Particle Electrification Observed in Volcanic Plumes
A special case of particle electrification aloft is that associated with charge separation in
volcanic plumes. Observations have shown that charging of terrestrial volcanic plumes
can occur up to hundreds of kilometres from the eruption site (Harrison et al. 2010;
Hatakeyama 1949). The presence of charge in the Earth’s electrically conductive atmosphere
indicates that self-charging occurs in the plume and is attributed to triboelectric charging.
Observations of volcanic lightning during the Eyjafjallajökull eruption in 2010 demonstrate
that ambient atmospheric conditions also contribute to electrical charging of terrestrial volcanic plumes (Bennett et al. 2010), known as the ash-rich ice electrification system (ARIES)
mechanism (Aplin et al. 2016). Convective processes, as in the ARIES mechanism, are expected to contribute to charging of dust storms on Mars, but without the enhancement associated with the presence of ice.
5 Models of charge exchange between particles
Simulations of the particle charging mechanism derive from laboratory measurements and
have been developed in order to understand a number of different phenomena. The basic
studies are described below and applications to dust devils follow.
Melnik and Parrot (1998) were interested in the conditions required for electric discharge
in dust disturbances on Mars. Their charging mechanism assumed that the average charge on
a particle was proportional to the particle size. They adopted the relation that each collision
resulted in charge exchange between particles of q = 1 fC/µm, with the larger particle
losing electrons to the smaller particle.
Desch and Cuzzi (2000) proposed an alternate mechanism of charging where both particle size and composition influenced the amount of charge exchange between the colliding
particles. While their simulations were conducted in order to understand the mechanisms for
lightning generation in the solar nebula, the concepts are transferrable to particle collisions
within dust storms or dust devils. The important parameter in their study is the contact potential difference between the two colliding particles. The contact potential of a material is
the energy needed to remove an electron from that material. When two particles of differing
material collide, the electrons will migrate to the material with the higher contact potential.
For metals, this quantity is directly related to the work function. For insulators, the contact
potential energy is more difficult to determine from intrinsic properties of the material, but
potential differences between two materials can be found from the triboelectric series (Faraday 1855). A formula was derived to calculate the charge exchange, which is a function of
potential difference and the mutual capacitance of the two materials. See Sect. 3.2 of Desch
and Cuzzi (2000) for the full derivation, resulting in
q = f1 Φ − (1 − f2 )qtot
(7)
where f1 = (c12 c21 − c11 c22 )/(c11 + c12 + c21 + c22 ) and f2 = (c11 + c12 )/(c11 + c12 + c21 +
c22 ) are functions of the mutual capacitances, cij . (Φ is the potential.) When mixing large insulator particles with small metallic particles, f2 approaches unity and the charge exchange
can be rewritten as a function of potential difference and particle size as
q ∼ 2668(Φ/2 V)(rf /0.5 µm)e,
(8)
324
R.G. Harrison et al.
where rf = (rL−1 + rS−1 )−1 ∼ rS is the reduced radius, and e is the elementary charge. (The
approximate potential difference between a metal and an insulator is 2 V.)
Kok and Renno (2008a, 2008b) modelled charge exchange between sand particles and
the Earth’s surface employing the Desch and Cuzzi charging mechanism, with an additional
term for collisions between particle pairs of similar composition,
Φeff = (rL − rS )/(rL + rS ),
(9)
where S is a scaling factor determined by calibrating the model with field measurements
(S = 6 ± 4 V). Their study was applied to reconcile discrepancies between classical saltation
theory and terrestrial observations of the phenomenon. By including the electric force in the
equation of motion for particles in saltation, they found that the wind stress needed to lift
particles is reduced and the particle trajectories remain closer to the surface with a reduced
horizontal speed.
Kok and Renno (2009) applied their terrestrial saltation model to explore saltation with
charged particles on Mars. They concluded that the breakdown E-field in the Martian saltation layer is ∼43 kV m−1 , but their modelling could only achieve fields of ∼15–20 kV m−1 .
Therefore, electric discharges in the saltation layer seem unlikely, but this does not prohibit
discharges from occurring in larger scale phenomena such as dust storms and dust devils.
5.1 Applications of the Charge Exchange Mechanism to Dust Devils
The Melnik-Parrot and Desch-Cuzzi grain charging processes described above were applied
to a comparative analytical study of dust devils on Earth and Mars in the work of Farrell et al.
(2003). They created an electrodynamic model of a dust devil by applying the methodology
used for modelling terrestrial thunderstorms, e.g. see Mathpal et al. (1980) to the induction
electrification process between graupel and water/ice in terrestrial thunderstorms. The time
derivative of the electric field, E, is related to rate of charging on the larger grains, QL , by
E + σ E /ε0 = −nL vQL /ε0
(10)
where σ is the local atmospheric conductivity, nL is the particle number density, v is the
relative velocity between the large/small grains (v = vL − vS < 0), and ε0 is the permittivity of free space. Here QL is based on triboelectric processes (for dust devils) rather than
induction processes (for thunderstorms) and is proportional to the charge exchange per collision, q, as QL = vq, where v is the grain collision frequency. The model is applied to
both Earth and Mars by varying the atmospheric conductivity appropriately. Results from
their simulations are shown in Fig. 12.
A number of aspects become evident from Fig. 12. First, exponentially growing E-fields
develop in both the terrestrial and Martian dust devils due to the tribocharging and subsequent mass (and charge) separation within the convective cloud. The electric fields in the
Martian case are always lower than in the terrestrial case due to greater dissipation leakage current which removes charge from the stratified grains. Additionally, the Desch and
Cuzzi charging results in electric fields consistently lower than those produced using Melnik and Parrot. The two cases are approximately equal for the Desch and Cuzzi Φ = 2 V
which represents the compositional difference of a metal and an insulator (i.e., simulating
the collision between a metallic grain and a basaltic grain).
A follow-up study by Farrell et al. (2006a) added the acceleration term to Eq. (10), such
that the right hand size becomes: −nL (QL v + QL v )/ε0 . The QL v acceleration term
depends on the lifting process associated with the fluid, and the addition of this term links
Applications of Electrified Dust and Dust Devil Electrodynamics. . .
325
Fig. 12 Comparison of the size-only grain charging mechanism from (a) Melnik and Parrot with the (b)
Desch and Cuzzi grain charging mechanism which also includes the effects of grain composition (from Figs. 2
and 3 in Farrell et al. 2003)
the electrostatic effects to the fluid properties of the medium, particularly the wind speed.
The effects are seen only during the early dust devil formation period, but are evident as a
rapid rise in E-field value by as much as a factor of 10 over that shown in Fig. 12 in the first
ten seconds. As time progresses, the QL v term primarily drives the system and the E-field
asymptotes to approach results similar to the cases shown in Fig. 12.
5.2 Simulation Codes Integrating Particle Charging with Vortex Dynamics
Melnik & Parrot (1998) conducted a simulation where a number of particles of like composition but differing sizes were placed in a box and subjected to wind movement. These
particles were then placed in a larger simulation domain (50 m × 100 m) where the wind
vortex was halted and further particle movement was controlled by the forces of gravity, air
drag, and the electrostatic field. The redistribution of charge was tracked over time until an
electric field built up to 20 kV m−1 , the breakdown field for the Martian atmosphere. The
maximum electric field seen in the simulation reached ∼10 kV m−1 in the first 2 s; the breakdown field was reached in 6.6 s. A second set of simulations tracking the movement of the
electric field as particles roll up a hill (under the influence of a wind parallel to the surface)
were also conducted. The electric field at the top of the hill reached the Martian breakdown
value in about 13.5 s; growth of the electric field was exponential and much more regular
than in the wind vortex case. In both wind environments, the amount of time to reach the
breakdown field could be considerably lengthened (by greater than a factor of 10) by decreasing the number of particles involved, the amount of charge exchange upon collision,
or the wind velocity. Over these longer timescales other processes (e.g. recombination of
charges through subsequent collisions) are likely to occur, which could prevent the electric
field increasing to breakdown strength.
Huang et al. (2008) numerically simulated a dust devil in a cylindrical domain of radius
100 m and height 200 m, in a terrestrial environment. They modelled the entire process
iteratively beginning with development of the convective vortex from local surface heating,
lifting of sand grains and their subsequent movement due to gravity and vortex winds, and
applying a specified charge-to-mass ratio (Huang and Zheng 2001) on each grain in order to
calculate the resulting electric field from Coulomb’s law. The closest fit to the observations
of Farrell et al. (2004) employed charge-to-mass ratios of −120 µC kg−1 , 60 µC kg−1 , and
57 µC kg−1 for sand grains of diameter 0.15 mm, 0.2 mm, and 0.25 mm, respectively. The
maximum E-field within the dust devil was reached in about 80 seconds of simulation time.
326
R.G. Harrison et al.
Fig. 13 Illustration of the forces
on the grains in the Macroscopic
Triboelectric Simulation model.
Along the z-direction, grains are
subject to an imposed constant
upward wind and the downward
gravitational force. In the
horizontal plane, the prevailing
winds are defined by a pattern of
cyclostrophic balance between
the outward centripetal force and
inward pressure gradient force.
(The z and xy motions are not
coupled.) (From Farrell et al.
2006b)
Fig. 14 3-D perspective of the grain motion in an MTS simulation. The plots correspond to (a) 10 time steps,
(b) 300 time steps, and (c) 499 time steps, where each timestep is approximately 1 ms (adapted from Fig. 4
in Farrell et al. 2006b)
Motivated by the detection of ULF magnetic radiation emitted in terrestrial dust devils
(Houser et al. 2003; Farrell et al. 2004), Farrell et al. (2006b) developed the Macroscopic
Triboelectric Simulation (MTS). MTS is a 3-D particle code which quantifies charging associated with swirling, mixing dust grains. Grains of pre-defined sizes and compositions are
placed in simulation box (2 m × 2 m × 1 m) and allowed to move under the influence of
winds and gravity. The model tracks the movement of individual grains in the prevailing
winds, and charge exchange upon grain-grain collision. An illustration of the model and a
3-D perspective of the grains during a simulation are shown in Fig. 13 and Fig. 14.
The grains are initially placed in distinct clumps in close proximity to one another surrounding the box centre. Grain sizes are randomly generated but constrained such that the
radius ranges 0.05–10 µm with a median value near 0.1 µm. To explore the Desch and Cuzzi
compositional charging effects, the smaller grains (r < 2 µm) are given a metallic composition (conductors) and the less abundant larger grains are modelled as silicates (insulators),
with the potential difference between the two populations of Φ = 1.6 V. The motion of the
smaller grains traces the wind pattern, while gravitational effects on the larger grains force
them into a less circular wind pattern and they move towards the walls of the box as they
fall. The distribution of charge on each of the grains at the point in the simulation before the
grains reach the walls of the box (t = 0.5 s) is shown in Fig. 15(a).
The spatial distribution of the charged grains establishes a distinct current system involving the relative displacement of negative and positive charge in the dust devil; the transverse
part of this current is responsible for a time-variable magnetic field, and the vertical part
is responsible for a large-scale vertical electrostatic E-field. The assumption of solenoidal
Applications of Electrified Dust and Dust Devil Electrodynamics. . .
327
Fig. 15 (a) Distribution of grain charge as a function of grain radius from an MTS simulation. The smaller
metallic grains charge negatively and the large silicate grains charge positively as prescribed by the D&C
charging mechanism. Additionally, and particularly evident with the metallic grains is the increase in grain
capacity to hold charge with increasing particle radius. (b) The magnetic field generated as a function of time
over the course of the same simulation (from Farrell et al. 2006b)
Fig. 16 The MTS model returns charge values on individual dust grains. For each size-bin in MRAMS (Mars
Regional Atmospheric Modeling System) the charge values were binned and fit to a log-normal function
(from Barth et al. 2015, with permission of Elsevier)
current flow has allowed the magnetic effects of dust devils to be modelled (Kurgansky et al.
2007), which can explain some aspects of the ULF magnetic observations.
The MTS work described above used prescribed winds and the results were scaled up to
the size of a real (terrestrial size) dust devil. Barth et al. (2015) applied the MTS simulations to Martian dust devils by incorporating the MTS-charged dust grains into their Mars
Regional Atmospheric Modeling System (MRAMS; see Chap. 7 for additional description).
The charged grains shown in Fig. 15 were partitioned into the MRAMS radius bins and a
log-normal function was fitted to the charge values in each radius bin as shown in Fig. 16.
Then the first and second moments of the charge distribution were added as tracers so that
the charge on the dust particle tracers could be reconstructed at any MRAMS grid point.
As dust devils developed in the MRAMS Large Eddy Simulations (LES), the temporal/spatial change in charge distribution within the simulated Martian dust devil could be
tracked and resulting E-fields were calculated. They found that the magnitude of the resulting E-field had strong dependence on particle size, amount of charge, and amount of dust
lifted; the range of E-field values (mV m−1 to kV m−1 ) was in general, more modest than
in previous studies (an example is shown in Fig. 17). As this study was a first look at the
charge environment, a number of processes were not yet included. However, with this code
one can examine the dust devil E-field in three dimensions, make estimates of the overall
328
R.G. Harrison et al.
Fig. 17 Time and vertical
variation of the electric field
(V/m) an observer would see if
traveling within the dust devil at
the minimum pressure point, as
generated by the MRAMS-MTS
model (from Barth et al. 2015,
with permission of Elsevier)
dipole moment, and get new insights on the development of higher order moments that form
in regions of inhomogeneous fluid flows.
5.3 Breakdown Effects
Large E-fields like those measured in situ in terrestrial dust devils have been found which
exceed 50 kV m−1 (Farrell et al. 2004; Renno et al. 2004; Jackson and Farrell 2006). These
same electrical convective features placed in the low pressure Martian atmosphere are expected to initiate atmospheric breakdown. Specifically, as the E-fields grow in the Martian
dust devil, the atmospheric currents transition from nominal flows described by J = σ E to
the creation of a current-enhancing electron avalanche from Eq. (6) as
Ja = n0 exp αT (E)d eve E
(11)
where n0 is the initial electron density, αT is Townsend’s first coefficient defining the number of ionizing events per unit length, and ve is the electron drift velocity in the gas. In
the electron avalanche, those electrons with energy greater than 14.5 eV can ionize a CO2
molecule via electron impact ionization to thus create extra electrons that get accelerated by
the driving E-field. The population of electrons above 14.5 eV and their overall drift speed
increases with electric field (Nighan 1970).
Farrell et al. (2015) noted that the electron avalanche represents a substantial modification of the atmospheric conductivity. At low E-fields, the atmospheric conductivity is represented as a bulk isotropic quantity, σ , where electron and ion mobility are assumed to be
nearly equal. In contrast, within the electron avalanche itself, the electron conductivity is
exponentially greater than that from the positive ion component, and is highly directional.
As such, electrical dissipation of any charged object has to be examined based on the flow
of fast beam-like electrons vs slow ions reaching the surface. Negatively charged dust grains
immersed within the electron avalanche will not quickly dissipate their surface charge but
instead reach a negative charge equilibrium state consistent with current balance to the grain
surfaces. The interruption of dissipation of the negative dust grains is likely to act to reinforce the overall dust charge build-up in the convective feature.
Applications of Electrified Dust and Dust Devil Electrodynamics. . .
329
Fig. 18 Electron chemistry
associated with (a) CO2 and
(b) H2 O, in the vicinity of an
electrified dust devil or storm
(after Delory et al. 2006)
Delory et al. (2006) modelled the electron avalanche process and growth in electron density in a dust feature and found that in dipole electric fields greater than 15 kV m−1 the mean
free path for electron impact ionization of less than a metre leads to substantial ionization.
These energized electrons can then interact with CO2 and H2 O, via dissociative attachment
creating CO, O− , OH, and H− (see Fig. 18). The dust devil-created CO and OH interactions are then suspected to generate hydrogen peroxide at concentrations well above those
predicted via photochemistry (Atreya et al. 2006). These excesses of hydrogen peroxide become adsorbed onto aerosols and ultimately end up in the regolith, giving the H2 O2 a longer
lifetime than the chemically formed H2 O2 , which is photochemically dissociated during the
day with a mean lifetime of 6 hours. Implications of this longer lifetime are that H2 O2 could
be an effective oxidiser (Encrenaz et al. 2012). The energetic electrons in the avalanche also
can possibly create methane loss via dissociative attachment (Farrell et al. 2006c).
5.3.1 Atmospheric Chemistry Effects and Breakdown
Delory et al. (2006) were the first to suggest that catalytic production of hydrogen peroxide
can be triggered by electrostatic fields generated in the Martian dust devils and dust storms
as well as near the surface within the saltation layer. The origin of this hydrogen peroxide
enhancement lies in the possibility for electrons, after acceleration by large electric fields, to
ionize CO2 and H2 O molecules, dissociating them and thereby releasing hydroxyl radicals
as oxidizers of the Martian atmosphere. Atreya et al. (2006) predict that the HOx enhancement factor can be up to 100 in electric fields reaching the breakdown limit (∼20 kV m−1 ).
However, no dedicated laboratory measurement to date has been performed in support of
these assertions and the chemical implications for such catalytic production of oxidizer remain difficult to reconcile with our present understanding of Martian chemistry (Lefèvre
and Forget 2009). The electron avalanche process including dust/electron absorption losses
and the active conversion to O− via electron dissociative attachment have been the focus
of several later modelling studies (Jackson et al. 2008, 2010; Kok and Renno 2009). All of
these models indicate that electron avalanche processes could develop, leading to electron
impact excitation and new chemical products. However, each varied in the intensity of the
electron impact ionization and chemical production. Delory et al. (2006) and Atreya et al.
(2006) predicted an electron avalanche and plasma content having values below 1 part in
1014 . Kok and Renno (2009) suggested that the activity could be of lower intensity when
including grain charge dissipation effects. Jackson et al. (2010) found that the final equilibrium values for the plasma (electrons, O− , CO+
2 ) depend upon the assumed saturation
process, with the most optimistic case to be 1 part in 104 , and in the least ideal case to be
1 part in 1013 . These models all uniformly suggest the electron avalanche and ionization
process is very mild, producing ions at concentrations 1 % of atmospheric gas densities.
While electron dissociative processes are considered a loss process for methane (Farrell
et al. 2006c), laboratory work has suggested that electrical activity associated with CO2 and
water could stimulate methane formation (Robledo-Martinez et al. 2012). Specifically, the
330
R.G. Harrison et al.
electrons dissociate CO2 and H2 O into CO/O− and OH/H− , with the CO and H2 then create methane and water. Farrell et al. (2015) concluded that the added energetic electrons in
the atmosphere associated with relatively modest dust-generated E-fields act to create new
chemistry, and in the case of the methane, can act to both destroy and create this critical
biomarker. An obvious remaining question regarding methane is which process dominates
at any given time, and how do these competing processes act to affect the overall global
methane stability. To date, models have described the development of a bulk E-field, and its
effect on storm-size chemistry. However, one can imagine the E-fields near individual dust
particles to be amplified via the converging field near a local point source. A 5 micron particle charged to ∼ −2.5 fC will have a local surface potential of −5 V (Farrell et al. 2015).
However, in its immediate vicinity the local E-field values can become very large. For example, within 30 µm the E-fields exceed about 27 kV m−1 , and within 10 µm in excess of
250 kV m−1 . Thus very active local chemistry in the vicinity of these particles could be expected. Specifically, the electrons will initially collect or be removed on an individual small
lofted particle depending upon the particle charge state and current balance to the particle
in the electron avalanche flow. However, as these particles reach their negative equilibrium
charge state, the electrons in the flow will be locally accelerated (i.e., redirected) about the
particle. This effect will stimulate local chemistry in the near-grain environment. In equilibrium, the negative charged particles will also have a cloud of positive CO+
2 ions about it,
but will tend to reject the O− ions. The O− ions should then congregate near the positively
charged surface and in the vicinity of the larger positively charged grains near the bottom
of the dust devil. The O− ion is highly reactive and should modify near-surface chemistry.
These avenues have yet to be fully examined.
6 Measurements of Martian Atmospheric Electricity
To provide a thorough picture of the atmospheric electricity on Mars, three principal parameters need to be measured in situ. These are the electrical conductivity, the DC atmospheric
electric field and the AC electric fields associated with electromagnetic wave emissions.
The atmospheric electric fields are directly related to atmospheric charging and discharging processes, and are characterized by large amplitudes, extending from DC to frequencies
representative of the dynamics of charging and discharging processes.
The electrical conductivity of an atmosphere is an essential parameter of its global electric circuit. It also provides information on the ionization processes that control the electrical
state of the atmosphere and on the major charged particles species, free electrons, positive
and negative ions that are present in the medium. The conductivities associated with either
positively or negatively charged species have to be measured separately to gain the relevant
information. From model calculations (Sect. 2.3), the anticipated conductivities range from
10−12 S m−1 to 10−10 S m−1 .
The quasi-DC electric fields provide direct evidence of the charging mechanisms and of
their temporal variations. Based on the properties of the Martian atmosphere, measurements
of the quasi DC electric fields should ideally cover the wide range from a few mV m−1 to
∼10 kV m−1 in order to encompass all situations of interest, from the faint electrification due
to photo-electrons charging to potential electrical breakdown. Impact charging processes are
directly associated with winds and thus their temporal variations at a fixed location should
arise from e.g. the motion of dust devils and/or turbulence in the wind flow with a frequency
of the order of 10 Hz. On the other hand, following observations in thunderstorms, discharge
processes and the resulting redistribution of charges give rise to large amplitude variations
Applications of Electrified Dust and Dust Devil Electrodynamics. . .
331
as rapid as 0.5 to 1 ms. Ideally, therefore, in situ measurements should cover a frequency
range up to ∼2 kHz.
Measuring the small amplitude AC electric fields will provide information on the electromagnetic wave emissions and the impact of charged dust grains on the lander. Electromagnetic waves can originate from atmospheric processes, such as the EM waves emitted
by filamentary discharge currents in case of atmospheric breakdown or the Schumann resonances that result from the trapping of these waves in the surface-ionosphere resonator.
They may also arise from totally different phenomena in the distant ionized environment of
Mars through interactions with the solar wind. To achieve these observations AC electric
fields have to be measured in the frequency range from a few Hz to more than 3 kHz, with
a sensitivity of better than µV m−1 Hz−1/2 . This upper limit of 3 kHz is estimated from the
greatest frequency signals expected, those of transverse resonances at ∼1500 Hz (Simões
et al. 2009), and consideration of sampling requirements.
Several proposals have been made for an atmospheric electricity package on the surface
of Mars (Berthelier et al. 2000; Farrell et al. 2000, 2004; Renno et al. 2004). The upcoming
ExoMars Entry Descent Module (EDM) will carry the DREAMS package to the surface,
which includes the MicroARES sensor that will provide an opportunity for the first-ever
in situ measurements of Martian atmospheric electricity. These observations are likely to
provide great insights into the Mars atmospheric electricity-chemistry-astrobiology interconnection. MicroARES is motivated by two fundamental science questions which are well
addressed by in situ measurements: (1) is there an electric field on Mars?, and (2) are Martian
dust devils electrified?
Ideally, an array of field mills could be deployed to answer these questions, operating indefinitely to sample a wide range of conditions. However, the very tight constraints imposed
on the dimensions, mass and power of the sensors on the EDM surface payload restrict the
electric field measurements to a short time window of only a few days, during which there is
some expectation of dust storms. As on Earth, the planetary large scale electric field of atmospheric origin is essentially vertical and this component is of primary interest. MicroARES
is based on a simplified version of the standard double probe technique used to measure
electric field in space. Rather than deploying a boom with two sensors, MicroARES was
designed to measure the potential with respect to the lander ground, using a single spherical electrode positioned 30 cm above the lander’s upper surface. Previous balloon flights of
double-probe instruments together with numerical models have shown that the intended MicroARES configuration provides reliable measurements of the local atmospheric potential
at the position of the electrode and higher frequency variations of electrostatic or electromagnetic origin.
6.1 MicroARES Electric Sensor
The MicroARES electric field sensor integrated in the ExoMars/DREAMS experiment (Esposito et al. 2014, 2015; Bettanini et al. 2014) will measure electrical properties of the
Martian atmosphere. A modified version of this sensor has also been developed and tested
in desert conditions, demonstrating its capability to measure electric fields found in terrestrial conditions, and under a simulated Martian conductivity environment. Because of the
fundamental scientific results in planetary atmospheric electricity anticipated from the ExoMars/MicroARES measurements, the sensor to be deployed is described in some detail
here.
In its final configuration, the sensor consists of a 27 cm mast supporting a 3 cm spherical
electrode (see Fig. 19). The electrode is free to float electrically to a local potential, from
332
R.G. Harrison et al.
Fig. 19 Photograph of electrode
Fig. 20 Concept of the electrode—atmosphere interaction used in the MicroARES measurement system.
The sheath around the electrode is represented by the Resistance Capacitance combination of Rs and Cs , and
the input parameters of the electronic pre-amplifier by Ri and Ci . RL (Relay Low range: 25 G ) and RH
(Relay High range: 10 G ) are range setting resistances, which can be switched in successively (RL first,
then RH) as necessary as the electric field increases
which the local electric field can be derived using known electrical characteristics of the
input stage. For weight and power consumption reasons, the reference electrode is the lander
chassis. The spherical sensing electrode system has been designed upon the presumption that
the Martian atmospheric conductivity lies between 10−10 and 10−13 S m−1 . Measurements
of the conductivity itself can also be made, from determining the charging time to restore the
equilibrium potential of the electrode after bringing it to a fixed potential (e.g. Aplin 2005;
Bennett and Harrison 2006).
The configuration used for the electrode-atmosphere interaction and associated electronic
system is summarised in Figure 20. The interaction of the sensing electrode with the atmosphere can be modelled by a very simple RC circuit characterized by the sheath impedance
(Rs and Cs ). Accurate modelling of the sheath around the electrode is therefore important in order to derive the conductivity from the measurement of Rs alone and vice-versa.
Two methods can be used to achieve this, firstly a mathematical model approach (MolinaCuberos et al. 2010, §2.2; Berthelier et al. 2000, §3.2), where the spherical electrode is
considered immersed in a collisional and weakly ionized medium, and secondly an empirical approach of modelling of the equipotential lines around the electrode in order to retrieve
the coefficients ke and kφ which link the sheath resistance to the atmospheric conductivity
(Seran et al. 2013).
Applications of Electrified Dust and Dust Devil Electrodynamics. . .
333
For the MicroARES sensor, the first approach is used, using
R s · Cs =
ε0
σ
with
Cs = 4πre ε0
and
Rs =
1
4πσ re
(12)
where re is the electrode radius, ε0 is the permittivity of free space and σ is the local atmospheric conductivity. Further simulations of the electrode coupling with its environment
will enable a comparison of the two approaches. With the Martian atmospheric conductivity
assumptions, Rs will be between 1010 and 1013 .
The input of the preamplifier is also modelled using a RC circuit. The Ri value is around
1014
(based on the input impedance of the LMC6041 operational amplifier used combined with some leakage across the board), and the circuit has been optimised to maximize
this parameter since Rs and Ri form a potential divider at the input. To allow measurement
of the DC potential of the electrode regardless of the atmospheric conductivity, Rs has to
be negligible compared to Ri . For large electric fields, dividers formed by the RL and RH
grounding resistors can be activated using controllable high impedance relays, similar to
the method used for terrestrial measurements by Harrison and Aplin (2000), using the same
model of operational amplifier. The RL and RH values of 25 G and 10 G respectively
have been chosen to suit the Martian atmospheric conductivity assumption and the corresponding sheath resistance Rs .
The range of measurable frequencies extends from DC to 3.2 kHz for sensor potential
values ranging across ±100 V, corresponding to vertical electric fields of ±300 V m−1 . At
low frequencies, the electrode is coupled to the undisturbed plasma through the resistance of
its surrounding sheath, to measure the DC electric fields. Above a cut-off frequency, which
depends on the sheath resistance and electrode capacitance, the preamplifier is capacitance
coupled to the ambient atmosphere (see Fig. 20). The high sheath resistance makes it necessary to use high input impedance (>1014 ), low leakage current (<10−14 A) electronic
preamplifiers.
To cope with possible large electric fields exceeding the normal range of voltages applied
to the electrode, the high impedance relay can switch in a parallel input resistance to the
preamplifier (Fig. 20). This resistance is significantly smaller than the anticipated value of
the electrode sheath resistance, which will therefore produce a potential divider, lowering
the potential to lie within the input voltage range of the preamplifier. In practice there are
two different shunt resistances that can be switched independently, 10 G and 25 G . The
measurements performed with these two shunt resistances allows the sheath resistance to be
determined, which in disturbed conditions, may differ from that obtained in quiet conditions.
A further measurement possible is that of the atmospheric conductivity through the relaxation technique: the sensor potential is displaced from the local floating potential and
then allowed to recover its equilibrium potential. The time constant for this charging (or
discharging) process τ = ε0 /σ readily provides the local electric conductivity σ , provided
that the conductive area of the lander is much larger than that of the sensor. This condition
ensures that the return current collected by the lander to counterbalance the current collected
by the polarized sensor does not entail any significant variation of the lander potential with
respect to the atmosphere, and ensures that the lander ground stays constant with respect to
the ambient medium.
334
R.G. Harrison et al.
Fig. 21 Terrestrial adaptation of
MicroARES antenna as used in
the Sahara desert
6.2 MicroARES Tests in the Terrestrial Atmosphere
MicroARES can also be used for the measurement of the electric properties of Earth’s atmosphere, and this provides a test opportunity. A terrestrial version of MicroARES (Fig. 21)
was developed and operated in the Sahara desert (at Merzouga, Morocco) for a week in July
2014 during the local dust storm season. It operated during the field campaign described in
Sect. 4.1 and fully reported in Esposito et al. (2016). As Fig. 2 shows, however, the atmospheric conductivity profiles are very different between Earth and Mars, with the surface
atmospheric conductivity at Mars comparable with that of the terrestrial stratosphere. Although space instrumentation can be tested using stratospheric balloons (Fulchignoni et al.
2004), use at this altitude would not be practical, not least as there would be no prospect
of obtaining a ground connection. Since the Earth’s atmospheric conductivity2 at sea-level
is about 100 times smaller than that near the surface on Mars, the antenna-atmosphere coupling will generate a greater sheath resistance of up to 1015 , which will cause the signal to
be almost fully attenuated and the electrode will not be able to maintain itself at the floating
local potential because of the small bias current flow to the input amplifier.
Since the standard MicroARES configuration will not properly measure the Earth’s DC
electric field, possible solutions for a test system are to either change the relay resistances
and the pre-amplifier input resistance or to change the antenna. Given the fact that the input resistance is an intrinsic property of the design, increasing it is clearly not possible.
Therefore, the antenna needs to be adapted to the Earth atmospheric properties instead.
For atmospheric conductivity roughly 100 times smaller than that on Mars, an electrode
roughly 100 times larger is needed in order to keep Rs at the same order of magnitude.
For the tests undertaken in the Sahara, for which there were numerous practical issues,
a 60 cm side copper cube constructed from 0.8 mm thickness sheets was used. Its 2.16 m2
surface area was equivalent to a 41.6 cm radius sphere, around 30 times larger than the
original antenna. Thus, by calculation, the associated Cs becomes around 46 pF, and Rs (with
assumed atmospheric conductivity of 50 fS m−1 ) will be around 3.5 T instead of 100 T .
With this modified sheath resistance, the input relays are expected to provide switching
of the input stage as originally planned, and the pre-amplifier input resistance should not
2 The conductivity considered here as the mean conductivity is the mobility weighted average of the bipolar
conductivities.
Applications of Electrified Dust and Dust Devil Electrodynamics. . .
335
Fig. 22 Time series of the
MicroARES unprocessed signal.
Measurements were acquired in
the Sahara desert, Morocco in
July 2014
significantly affect the measurements. The mast height was also increased to 50 cm, which
lowers the parasitic capacity to a few pF. A further consequence is that the centre of the cube
is now located at 80 cm from the ground, which will lead to a greater measured potential
(Fig. 21).
As the cube antenna and Earth atmosphere have different electrical properties from
the anticipated Martian input configuration, some different post-processing is required.
The response time of the instrument input is inversely proportional to the conductivity
(τ = ε0 /σ = 1/Rp Cp ), hence the expected response time on Earth with this cube antenna
will be ∼200 s. The two main consequences of this are that the signal is highly deformed,
which can be corrected numerically, and that, in the conductivity measurement mode,
a switching sequence as rapid as 4 s is no longer possible. Since the post-processing correction depends on σ , its determination is performed with the signal correction. The instrument
can only record 20 min worth of data, so it is restarted, with the signal correction applied
every 20 min.
To reconstruct the input response, the instrument input stage parameters are needed accurately. The model employed for the electronics is that presented in
Figure 20, with Ri = 10 T , Ci = 2 pF (preamplifier) + 6 pF (antenna), RL = 25 G and
RH = 10 G . The antenna-atmosphere coupling depends on two parameters, Cs and Rs , but
the measurements only provide one parameter, the atmospheric conductivity σ . Since the
equations describing the instrument input depend on these variables separately, an equivalent radius to a spherical electrode (41.6 cm) is assumed to derive Cs and Rs from the
conductivity.
The signal distortion at each zeroing of the instrument is clearly visible on Fig. 22. Note
that at every restart (switching on and off of the signal processing board), the input of the
preamplifier is connected to ground and then consequently discharged, so the slow response
time of the atmospheric coupling is exaggerated.
In order to correct the signal the response is modelled with the assumptions made previously. The data reconstruction depends on the determination of σ (for the values of Rs
and Cp ) and the initial potential condition i.e. Ue (0). On Mars, σ would have been provided
by the conductivity measurement mode, during calm weather, and the relay operation during
stormy weather (the quick response time allows computation of the ratio before and after
relay operations to easily determine Rs and σ ), but for Earth measurements it is not available
directly, without the post processing. The data reconstruction also relies on the assumption
that, during the 6 s reset time, the electric field is steady, thus allowing a good approximation
of Ue (0) during each 20 min of measurement.
336
R.G. Harrison et al.
Fig. 23 Discontinuity
smoothing with simultaneous
determination of conductivity.
The time axis is local time
Determination of σ has been achieved by iterating its value in the input inversion equations, at every switching cycle of the instrument visible on Fig. 22, so that the derivative
and potential changes before and after the switch are minimized (see Fig. 23). Comparison
of the reconstructed electric field with that obtained from an adjacent commercial field-mill
(Fig. 24) demonstrates the effectiveness of this technique, with a correlation of 0.9 during
calm weather and 0.75 during stormy weather.
As explained, the MicroARES measurement sequence and post-processing also allow
the retrieval of conductivity, which is generally not possible with commercial field-mills
alone running in a standard mode. Figure 25 shows the retrieved conductivity for 9–13th
July 2014; the scale chosen emphasises the variability. The increase of conductivity during
dust events implies that the dust is contributing to the conductivity, as discussed for dense
aerosol concentrations by Dhanorkar and Kamra (1997).
7 Conclusions
Dusty environments on Earth can, as shown here, readily become electrified, and, by analogy
and theoretical consideration, the Martian atmosphere also seems likely to contain electrified dust. Central to developing an improved understanding of this dust electrification is the
understanding of charge exchange between the substances concerned, and possible breakdown effects. Ultimately, bulk measurements made in natural dust storms on Earth and Mars
need to be reconciled with laboratory measurements of electrification at a particle level and
the modelled electric field structures they present under the appropriate electrodynamic conditions.
Measurements of terrestrial dust electrification at the surface and aloft indicate that particle charging in dust storms is a common, and, perhaps, even a universal phenomenon to a
greater or lesser extent. The use of field mill measurements in dust devils shows that electric fields of order 10 kV m−1 are readily observed, which are associated with the pressure
minimum and wind speed maximum (e.g. Fig. 8). Larger values to ∼100 kV m−1 have also
been reported. The magnitude of the field appears to be reduced with increasing humidity
(Fig. 10), which is likely to be related to the properties of the contact electrification process
when water layers at the molecular level are present (Sect. 3.4).
Several modern experimental techniques are proving to be extremely informative with
regard to contact electrification processes, notably, as mentioned, atomic force microscopy
Applications of Electrified Dust and Dust Devil Electrodynamics. . .
337
Fig. 24 Electric field samples made at 1 Hz: comparison during calm (upper panel: 10th July 2014) and
disturbed (lower panel: 12th July 2014) weather. Note that the electric field intensity measured during fair
(calm) weather by MicroARES at 80 cm from the ground is smaller than the value measured by the reference
field mill that was placed at 2 m from the ground. This is in agreement with the expected atmospheric electric
field vertical profile
(AFM) (Matsusaka et al. 2010; Horn et al. 1993, Gady et al. 1998). This approach promises
to provide a more detailed (atomic scale) understanding of electrification and has already
led to the discovery that net surface electrification in fact is made up of a ‘mosaic’ of highly
electrified regions of positive and negative patches (Baytekin et al. 2011). Both processes of
electron and ion transfer can be involved in contact electrification, depending on the surface
conditions or composition. In the terrestrial environment, surface water layers will typically
be expected on particles facilitation ion exchange irrespective of surface composition.
For the electrostatics in dust devils, specific studies aimed at understanding and quantifying the dust or sand electrification interactions are of direct practical use. Several laboratory (and wind tunnel) based studies are currently being performed (Bo et al. 2014;
Lacks and Sankaran 2011, Merrison et al. 2012). To allow modelling of electrification
338
R.G. Harrison et al.
Fig. 25 Conductivity inferred with MicroARES instrument in terrestrial desert conditions from 9–13th July
2014
within a dust devil, a deeper understanding is needed in the size dependence of contact
electrification. This is probably only likely to be achieved by pursuing more advanced
laboratory experiments, applying detailed physical models (possibly involving both electron and ion transfer processes). Detailed and novel field studies could also contribute
here.
Modelling of dust devils has shown that, by linking the flow dynamics produced by
the balance between centripetal and pressure gradient forces (i.e. cyclostrophic balance)
with theoretically derived inter-particle charge exchange processes, leads to electric fields
comparable with those observed (Fig. 12). This includes a rapid initial rise in the electric
field. At extreme electric fields, the dynamics of the charged particles can be influenced by
the field (Kok and Renno 2008a, 2008b), providing a feedback on the electrification or a
means by which the effectiveness of dust uplift is enhanced. Ultimately the growth of the
field is limited by breakdown processes, which are dependent on the chemical composition
of the atmosphere concerned.
Within the Martian atmosphere there are good arguments to expect electrification of dust
devils and dust storms, in the absence of any direct electrical measurements to date. The
DREAMS-MicroARES instrument due to land on Mars in 2016 will provide these first
direct measurements, and with them, an important and substantial change in the knowledge
of its planetary electrical environment.
Acknowledgements KAN acknowledges the support of the UK’s Natural Environment Research Council
through an Independent Research Fellowship (NE/L011514/1). This work was facilitated by a workshop on
dust devils hosted by the International Space Science Institute in Bern.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution,
and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source,
provide a link to the Creative Commons license, and indicate if changes were made.
Applications of Electrified Dust and Dust Devil Electrodynamics. . .
339
References
M.M. Anderson et al., The Allen Telescope Array search for electrostatic discharges on Mars. Astrophys. J.
744, 15 (2012). doi:10.1088/0004-637X/744/1/15
K.L. Aplin, Aspirated capacitor measurements of air conductivity and ion mobility spectra. Rev. Sci. Instrum.
76, 104501 (2005). doi:10.1063/1.2069744
K.L. Aplin, Atmospheric electrification in the solar system. Surv. Geophys. 27, 63–108 (2006).
K.L. Aplin, R.G. Harrison, A self-calibrating programable mobility spectrometer for atmospheric ion measurements. Rev. Sci. Instrum. 72(8), 3467–3469 (2001)
K.L. Aplin, R.G. Harrison, M.J. Rycroft, Investigating Earth’s atmospheric electricity: a role model for planetary studies. Space Sci. Rev. 137, 11–27 (2008). doi:10.1007/s11214-008-9372-x
K.L. Aplin, T. Goodman, K.L. Herpoldt, C.J. Davis, Laboratory analogues of Martian electrostatic discharges.
Planet. Space Sci. 69, 100–104 (2012). doi:10.1016/j.pss.2012.04.002
K.L. Aplin, I.M.P. Houghton, K.A. Nicoll, M. Humphries, A. Tong, Electrical charging of volcanic ash, in
Proc. ESA Annual Meeting on Electrostatics (2014)
K.L. Aplin, A.J. Bennett, R.G. Harrison, I.M.P. Houghton, Electrostatics and in situ sampling of volcanic
plumes, in Volcanic Ash: Methods of Observation and Monitoring, ed. by S. Mackie (Elsevier, Amsterdam, 2016)
M.M. Apodaca, P.J. Wesson, K.J.M. Bishop, M.A. Ratner, B.A. Grzybowski, Contact electrification between
identical materials. Angew. Chem., Int. Ed. Engl. 49, 946–949 (2010)
S.K. Atreya, A.-S. Wong, N.O. Renno, W.M. Farrell, G.T. Delory, D.D. Sentman, S.A. Cummer, J.R. Marshall, S.C.R. Rafkin, D.C. Catling, Oxidant enhancement in Martian dust devils and storms: implications
for life and habitability. Astrobiology 6(3), 439–450 (2006) doi:10.1089/ast.2006.6.439
P.F.H. Baddeley, Whirlwinds and dust-storms of India (Bell and Daldy, London, 1860)
R.A. Bagnold, Auto-suspension of transported sediment; turbidity currents. Proc. R. Soc. Lond. Ser. A
265(1322), 315–319 (1962)
A.G. Bailey, The charging of insulator surfaces. J. Electrost. 51/52, 82–90 (2001)
E.L. Barth, W.M. Farrell, S.C.R. Rafkin, Electric field generation in Martian dust devils. Icarus 268, 253–265
(2015). doi:10.1016/j.icarus.2015.12.020
H.T. Baytekin, A.Z. Patashinski, M. Branicki, B. Baytekin, S. Soh, B.A. Grzybowski, The mosaic of surface
charge in contact electrification. Sci. Express (2011). doi:10.1126/science.1201512
A.J. Bennett, R.G. Harrison, In situ calibration of atmospheric air conductivity measurements. Rev. Sci.
Instrum. 77, 016103 (2006)
A.J. Bennett, R.G. Harrison, Surface measurement system for the atmospheric electrical vertical conduction
current density, with displacement current correction. J. Atmos. Sol.-Terr. Phys. 70, 1373–1381 (2008).
doi:10.1016/j.jastp.2008.04.014
A.J. Bennett, P. Odams, D. Edwards, TH. Arason, Monitoring of lightning from the April–May 2010 Eyjafjallajökull volcanic eruption using a very low frequency lightning location network. Environ. Res. Lett.
5(4), 044013 (2010)
J.-J. Berthelier, R. Grard, H. Laakso, M. Parrot, ARES, Atmospheric relaxation and electric field sensor,
the electric field experiment on NETLANDER. Planet. Space Sci. 48, 1193–1200 (2000). doi:10.1016/
S0032-0633(00)00103-3
C. Bettanini, F. Esposito, S. Debei, C. Molfese, I. Arruego Rodríguez, G. Colombatti, A.-M. Harri,
F. Montmessin, C. Wilson, A. Aboudan, S. Abbaki, V. Apestigue, G. Bellucci, J.-J. Berthelier, J.R. Brucato, S.B. Calcutt, F. Cortecchia, G. Di Achille, F. Ferri, F. Forget, G.P. Guizzo, E. Friso, M. Genzer, P. Gilbert, H. Haukka, J.J. Jiménez, S. Jiménez, J.-L. Josset, O. Karatekin, G. Landis, R. Lorenz,
J. Martinez, L. Marty, V. Mennella, D. Möhlmann, D. Moirin, R. Molinaro, E. Palomba, M. Patel,
J.-P. Pommereau, C.I. Popa, S. Rafkin, P. Rannou, N.O. Renno, P. Schipani, W. Schmidt, S. Silvestro,
F. Simoes, A. Spiga, F. Valero, L. Vázquez, F. Vivat, O. Witasse, R. Mugnuolo, S. Pirrotta E:Marchetti,
The DREAMS experiment on the ExoMars 2016 mission for the study of Martian environment during the dust storm season, in Proceedings of IEEE Metrology for Aerospace (MetroAeroSpace) (2014),
pp. 167–173. doi:10.1109/MetroAeroSpace.2014.6865914
M.A. Bilici, J.R. Toth II., R.H. Sankaran, D. Lacks, Particle size effects in particle-particle triboelectric charging studied with an integrated fluidized bed and electrostatic separator. Rev. Sci. Instrum. 85, 103903
(2014)
T. Bo, H. Zhang, W. Hu, X. Zheng, The analysis of electrification in windblown sand. Aeolian Res. 11, 15–21
(2013a)
T.L. Bo, H. Zhang, W. Zhu, X.J. Zheng, Theoretical prediction of electric fields in wind-blown sand. J. Geophys. Res. 118, 4494–4502 (2013b)
T. Bo, H. Zhang, X. Zheng, Charge-to-mass ratio of saltating particles in wind-blown sand. Sci. Rep. 4, 5590
(2014)
340
R.G. Harrison et al.
E. Chassefière, P. Drossart, O. Korablev, Post-Phobos model for the altitude and size distribution of dust in
the low Martian atmosphere. J. Geophys. Res. 100(E3), 5525–5539 (1995)
C. Cimarelli, M.A. Alatorre-Ibargüengoitia, U. Kueppers, B. Scheu, D.B. Dingwell, Experimental generation
of volcanic lightning. Geology 42(1), 79–82 (2014)
C.F. Clement, R.G. Harrison, The charging of radioactive aerosols. J. Aerosol Sci. 23(5), 481–504 (1992)
W.D. Crozier, The electric field of a New Mexico dust devil. J. Geophys. Res. 69, 24 (1964). doi:10.1029/
JZ069i024p05427
W.D. Crozier, Dust devil properties. J. Geophys. Res. 75, 4583–4585 (1970). doi:10.1029/JC075i024p04583
D.K. Davies, Charge generation on dielectric surfaces. Br. J. Appl. Phys. 2, 1533–1537 (1969)
G.T. Delory, W.M. Farrell, S.K. Atreya, N.O. Renno, A.-S. Wong, S.A. Cummer, D.D. Sentman, J.R. Marshall, S.C.R. Rafkin, D.C. Catling, Oxidant enhancement in Martian dust devils and storms: storm electric fields and electron dissociative attachment. Astrobiology 6(3), 451–462 (2006). doi:10.1089/ast.
2006.6.451
L. Demon, P. DeFelice, H. Gondet, Y. Kast, L. Pontier, Premiers résultats obtenus au cours du printemps
J. Rech. Cent. Natl Rech. Sci. 24, 126–137 (1953)
S.J. Desch, J.N. Cuzzi, The generation of lightning in the solar nebula. Icarus 143, 87–105 (2000)
S. Dhanorkar, A.K. Kamra, Calculation of electrical conductivity from ion aerosol balance equations. J. Geophys. Res. 102(D25), 30147–30159 (1997)
A.F. Diaz, R.M. Felix-Navarro, A semi-quantitative tribo-electric series for polymeric materials: the influence
of chemical structure and properties. J. Electrost. 62, 277 (2004)
N. Duff, D.J. Lacks, Particle dynamics simulations of triboelectric charging in granular insulator systems.
J. Electrost. 66, 51 (2008). doi:10.1016/j.elstat.2007.08.005
C.B. Duke, T.J. Fabish, Contact electrification of polymers: a quantitative model. J. Appl. Phys. 49, 315–321
(1978)
H.F. Eden, B. Vonnegut, Electrical breakdown caused by dust motion in low pressure atmospheres: considerations for Mars. Nature 280, 962–963 (1973)
T. Encrenaz, T.K. Greathouse, F. Lefèvre, S.K. Atreya, Hydrogen peroxide on Mars: observations, interpretation and future plans. Planet. Space Sci. 68(1), 3–17 (2012)
F. Esposito, S. Debei, C. Bettanini, C. Molfese, I. Arruego Rodriguez, G. Colombatti, A.M. Harri,
F. Montmessin, C. Wilson, A. Aboudan, S. Abbaki, V. Apestigue, G. Bellucci, J.J. Berthelier, J.R. Brucato, S.B. Calcutt, F. Cortecchia, F. Cucciarrè, G. Di Achille, F. Ferri, F. Forget, E. Friso, M. Genzer,
H. Haukka, J.J. Jimènez, S. Jimènez, J.L. Josset, O. Karatekin, G. Landis, R. Lorenz, E. Marchetti,
J. Martinez, L. Marty, V. Mennella, D. Möhlmann, D. Moirin, R. Molinaro, E. Palomba, M. Patel,
J.P. Pommereau, C.I. Popa, S. Rafkin, P. Rannau, N.O. Renno, P. Schipani, W. Schmidt, E. Segato,
S. Silvestro, F. Simoes, A. Spiga, F. Valero, L. Vázquez, F. Vivat, O. Witasse, R. Mugnuolo, S. Pirrotta,
in The DREAMS Experiment of the ExoMars 2016 Mission for the Study of Martian Environment During the Dust Storm Season, Eighth International Conference on Mars, LPI Contribution No. 1791, July
14–18, 2014, Pasadena, California (2014), p. 1246
F. Esposito, S. Debei, C. Bettanini, C. Molfese, I. Arruego Rodríguez, G. Colombatti, A.-M. Harri,
F. Montmessin, C. Wilson, A. Aboudan, S. Abbaki, V. Apestigue, G. Bellucci, J.-J. Berthelier, J.R. Brucato, S.B. Calcutt, F. Cortecchia, F. Cucciarrè, G. Di Achille, F. Ferri, F. Forget, E. Friso, M. Genzer, P. Gilbert, H. Haukka, J.J. Jiménez, S. Jiménez, J.-L. Josset, O. Karatekin, G. Landis, R. Lorenz,
M. Mannetta, J. Martinez, L. Marty, V. Mennella, D. Möhlmann, D. Moirin, R. Molinaro, E. Palomba,
M. Patel, J.-P. Pommereau, C.I. Popa, S. Rafkin, P. Rannou, N.O. Renno, P. Schipani, W. Schmidt,
E. Segato, S. Silvestro, F. Simoes, E. Simoncini, A. Spiga, F. Valero, L. Vázquez, F. Vivat, O. Witasse,
R. Mugnuolo, S. Pirrotta, E. Marchetti, in The DREAMS Experiment On-Board the Schiaparelli Lander
of ExoMars Mission. EPSC Abstracts, vol. 10 (2015), EPSC2015-364
F. Esposito, R. Molinaro, C.I. Popa, C. Molfese, F. Cozzolino, L. Marty, K. Taj-Eddine, G. Di Achille, S.
Silvestro, G.G. Ori, The role of atmospheric electric field in the dust lifting process. Geophys. Res. Lett.
(2016, in revision)
M. Faraday, Experimental Researches in Electricity (Dover, New York, 1855)
W.M. Farrell, M.D. Desch, Is there a Martian atmospheric electric circuit. J. Geophys. Res. 106(E4), 7591–
7595 (2001)
W.M. Farrell, M.L. Kaiser, M.D. Desch, J.D. Houser, S.A. Cummer, D.M. Wilt, G.A. Landis, Detecting
electrical activity from Martian dust storms. J. Geophys. Res. 104(2), 3795–3801 (1999)
W.M. Farrell et al., Radio and optical detection of Martian dust storm discharges. Acta Astronaut. 46, 25–36
(2000)
W.M. Farrell, G.T. Delory, S.A. Cummer, J.R. Marshall, A simple electrodynamic model of a dust devil.
Geophys. Res. Lett. 30, 2050 (2003)
W.M. Farrell et al., Electric and magnetic signatures of dust devils from the 2000–2001 MATADOR desert
tests. J. Geophys. Res. 109, E03004 (2004). doi:10.1029/2003JE002088
Applications of Electrified Dust and Dust Devil Electrodynamics. . .
341
W.M. Farrell, N. Renno, G.T. Delory, S.A. Cummer, J.R. Marshall, Integration of electrostatic and fluid
dynamics within a dust devil. J. Geophys. Res. 111, E01006 (2006a)
W.M. Farrell, J.R. Marshall, S.A. Cummer, G.T. Delory, M.D. Desch, A model of the ULF magnetic and
electric field generated from a dust devil. J. Geophys. Res. 111, E11004 (2006b)
W.M. Farrell et al., Martian dust storms as a possible sink to atmospheric methane. Geophys. Res. Lett. 33,
L21203 (2006c)
W.M. Farrell, J.L. McLain, M.R. Collier, J.W. Keller, T.J. Jackson, G.T. Delory, Is the electron avalanche
process in a Martian dust devil self-quenching? Icarus 254, 333–337 (2015)
D.C. Ferguson, J.C. Kolecki, M.W. Siebert, D.M. Wilt, J.R. Matijevic, Evidence for Martian electrostatic
charging and abrasive wheel wear from the Wheel abrasion experiment on the Pathfinder Sojourner
rover. J. Geophys. Res. 104(E4), 8747–8759 (1999)
M. Fillingim, Global electric circuit of Mars. Unpublished report, Florida Institute of Technology (1998)
K.M. Forward, D.J. Lacks, R.M. Sankaran, Particle-size dependent bipolar charging of Martian regolith simulant. Geophys. Res. Lett. 36, L13201 (2009a). doi:10.1029/2009GL038589
K.M. Forward, D.J. Lacks, R.M. Sankaran, Charge segregation depends on particle size in triboelectrically
charged granular materials. Phys. Rev. Lett. 102, 028001 (2009b)
G.D. Freier, The electric field of a large dust devil. J. Geophys. Res. 65(10), 3504 (1960). doi:10.1029/
JZ065i010p03504
M. Fulchignoni, A. Aboudan, F. Angrilli, M. Antonello, S. Bastianello, C. Bettanini, G. Bianchini, G. Colombatti, F. Ferri, E. Flamini, V. Gaborita, N. Ghafoor, B. Hathi, A.-M. Harri, A. Lehto, P.F. Lion Stoppato,
M.R. Patel, J.C. Zarnecki, A stratospheric balloon experiment to test the Huygens atmospheric structure
instrument (HASI). Planet. Space Sci. 52, 867–880 (2004)
B. Gady, R. Reifenberger, D.S. Rimai, Contact electrification studies using atomic force microscope techniques. J. Appl. Phys. 84(1), 219–322 (1998)
R. Greeley, N. Lancaster, S. Lee, P. Thomas, Martian aeolian processes, sediments, and features. Mars 1,
730–766 (1992)
J.S. Gilbert, S.J. Lane, R.S.J. Sparks, T. Koyaguchi, Charge measurements on particle fallout from a volcanic
plume. Nature 349(6310), 598–600 (1991)
R. Grard, Solar photon interaction with the Martian surface and related electrical and chemical phenomena.
Icarus 114, 130–138 (1995)
W. Gringel, Untersuchungen zur elektrischen Luftleitfahigkeit unter Berucksichtihung der Sonnenaktivitat
und der Aerosolteilchenkonzentration bis 35 km Höehe. PhD thesis, Universität Tübingen (1978)
W. Gringel, R. Mühleisen, Sahara dust concentration in the troposphere over the North Atlantic derived from
measurements of air conductivity. Beitr. Phys. Atmos. 51(2), 121–128 (1978)
F.B. Gross, S.B. Grek, I.C. Calle, R.U. Lee, JSC Mars-1 Martian regolith simulant particle charging experiments in low pressure environment. J. Electrost. 53, 257–266 (2001)
B.A. Grzybowski, M. Fialkowski, J.A. Wiles, Kinetics of contact electrification between metals and polymers.
J. Phys. Chem. B 109, 20511–20515 (2005)
D.A. Gurnett, D.D. Morgan, L.J. Granroth, B.A. Cantor, W.M. Farrell, J.R. Espley, Non-detection of impulsive radio signals from lightning in Martian dust storms using the radar receiver on the mars express
spacecraft. Geophys. Res. Lett. 37, L17802 (2010). doi:10.1029/2010GL044368
D. Halleaux, N.O. Renno, Aerosols–climate interactions at the Owens “Dry” Lake, California. Aeolian Res.
15, 91–100 (2014). doi:10.1016/j.aeolia.2014.08.004
W.R. Harper, Contact and Frictional Electrification (Laplacian Press, Morgan Hill, 1998)
D.J. Harris, Electrical effects of the Harmattan dust storms. Nature 214, 585 (1967)
R.G. Harrison, An antenna electrometer system for atmospheric electrical measurements. Rev. Sci. Instrum.
68(3), 1599–1603 (1997)
R.G. Harrison, The global atmospheric electrical circuit and climate. Surv. Geophys. 25(5–6), 441–484
(2004). doi:10.1007/s10712-004-5439-8
R.G. Harrison, The cloud chamber and CTR Wilson’s legacy to atmospheric science. Weather 66(10), 276–
279 (2011). doi:10.1002/wea.830
R.G. Harrison, The Carnegie curve Surv. Geophys. 34(2), 209–232 (2013). doi:10.1007/s10712-012-9210-2
R.G. Harrison, K.L. Aplin, A multimode electrometer for atmospheric ion measurements. Rev. Sci. Instrum.
71(12), 4683–4685 (2000)
R.G. Harrison, H. Tammet, Ions in the terrestrial atmosphere and other solar system atmospheres. Space Sci.
Rev. 137, 107–118 (2008). doi:10.1007/s11214-008-9356-x
R.G. Harrison, K.L. Aplin, F. Leblanc, Y. Yair, Planetary atmospheric electricity. Space Sci. Rev. 137, 5–10
(2008). doi:10.1007/s11214-008-9419-z
R.G. Harrison, K.A. Nicoll, Z. Ulanowski, T.A. Mather, Self-charging of the Eyjafjallajökull volcanic ash
plume. Environ. Res. Lett. 5(2), 024004 (2010). doi:10.1088/1748-9326/5/2/024004
342
R.G. Harrison et al.
H. Hatakeyama, On the disturbance of the atmospheric potential gradient caused by the smoke-cloud of the
Volcano Yake-Yama. J. Geomagn. Geoelectr. 1(2), 48–51 (1949)
C. Helling, R.G. Harrison, F. Honary, D.A. Diver, K.L. Aplin, I. Dobbs-Dixon, U. Ebert, S. Inutsuka, F.J.
Gordillo-Vazquez, S. Littlefair, Atmospheric electrification in dusty, reactive gases in the solar system
and beyond. Surv. Geophys. (2016, in press). doi:10.1007/s10712-016-9361-7
J.R. Herman, J.A. Caruso, R.G. Stone, Radio astronomy explorer (RAE)—1. Observations of terrestrial radio
noise. Planet. Space Sci. 21, 443–461 (1973)
R.G. Horn, D.T. Smith, A. Grabbe, Contact electrification induced by monolayer modification of a surface
and relation to acid-base interactions. Nature 366, 442–443 (1993)
I.M.P. Houghton, K.L. Aplin, K.A. Nicoll, Triboelectric charging of volcanic ash from the 2011 Grímsvötn
eruption. Phys. Rev. Lett. 111, 118501 (2013). doi:10.1103/PhysRevLett.111.118501
J.G. Houser et al., ULF and ELF magnetic activity from a terrestrial dust devil. Geophys. Res. Lett. 30, 1027
(2003)
N. Huang, X. Zheng, A laboratory test of the electrification phenomenon in wind-blown sand flux. Chin. Sci.
Bull. 46, 417–420 (2001)
N. Huang, G. Yue, X. Zheng, Numerical simulations of a dust devil and the electric field in it. J. Geophys.
Res. 113, D20203 (2008)
I.I. Inculet, G.S.P. Castle, G. Aartsen, Generation of bipolar electric fields during industrial handling of powders. Chem. Eng. Sci. 61, 2249–2253 (2006)
T.L. Jackson, W.M. Farrell, Electrostatic fields in dust devils: an analog to Mars. IEEE Trans. Geosci. Remote
Sens. 44(10), 2942–2949 (2006). doi:10.1109/TGRS.2006.875785
T.L. Jackson, W.M. Farrell, G.T. Delory, J. Nithianandam, Effect of dust absorption on the electron
avalanche process occurring within Martian dust storms. Geophys. Res. Lett. 35, L16201 (2008).
doi:10.1029/2008GL034523
T.L. Jackson et al., Martian dust devil electron avalanche process and associated electrochemistry. J. Geophys.
Res. 115, E05006 (2010)
M.R. James, L. Wilson, S.J. Lane, J.S. Gilbert, T.A. Mather, R.G. Harrison, R.S. Martin, Electrical charging
of volcanic plumes, in Planetary Atmospheric Electricity (Springer, New York, 2008), pp. 399–418
A.R. Johnston, H. Kirkham, B. Eng, DC electric field meter with fiber optic readout. Rev. Sci. Instrum. 57,
2746–2753 (1986)
A.K. Kamra, Measurements of the electrical properties of dust storms. J. Geophys. Res. 77(30), 5856 (1972)
J.F. Kok, D.J. Lacks, Electrification of granular systems of identical insulators. Phys. Rev. E 79, 051304
(2009)
J.F. Kok, N.O. Renno, Enhancement of the emission of mineral dust aerosols by electric forces. Geophys.
Res. Lett. 33, L19S10 (2006)
J.F. Kok, N.O. Renno, Electrostatics in wind-blown sand. Phys. Rev. Lett. 100, 014501 (2008a)
J.F. Kok, N.O. Renno, Electrical activity and dust lifting on Earth, Mars, and beyond. Space Sci. Rev. 137,
419–434 (2008b)
J.F. Kok, N.O. Renno, Electrification of wind-blown sand on Mars and its implications for atmospheric chemistry. Geophys. Res. Lett. 36, L05202 (2009). doi:10.1029/2008GL036691
C.E. Krauss, M. Horanyi, S. Robertson, Experimental evidence for electrostatic discharging of dust near the
surface of Mars. New J. Phys. 5, 70.1–70.9 (2003)
C.E. Krauss, M. Horanyi, S. Robertson, Modeling the formation of electrostatic discharges on Mars. J. Geophys. Res. 111, E02001 (2006). doi:10.1029/2004JE002313
W.B. Kunkel, The static electrification of dust particles on dispersion into a cloud. J. Appl. Phys. 21, 820–832
(1950)
M.V. Kurgansky, L. Baez, E.M. Ovalle, A simple model of the magnetic emission from a dust devil. J. Geophys. Res. 112, E11008 (2007). doi:10.1029/2007JE002952
D.J. Lacks, A. Levandovsky, Effect of particle size distribution on the polarity of triboelectric charging in
granular insulator systems. J. Electrost. 65, 107–112 (2007)
D.J. Lacks, R.M. Sankaran, Contact electrification of insulating materials. J. Phys. D. 44, 45 (2011)
J.D. Lacks, N. Duff, S.K. Kumar, Nonequilibrium accumulation of surface species and triboelectric charging
in single component particulate systems. Phys. Rev. Lett. 100, 188305 (2008)
J. Latham, The electrification of snowstorms and sandstorms. Q. J. R. Meteorol. Soc. 90(383), 91–95 (1964)
M.A. Laughton, D.F. Warne (eds.), Electrical Engineer’s Reference Book, 16th edn. (2004), Newnes
K.Y. Law, I.W. Tarnawskyj, D. Salamida, T. Debies, Investigation of the contact charging mechanism between
an organic salt doped polymer surface and polymer-coated metal beads. Chem. Mater. 7(11), 2090–2095
(1995)
F. Lefèvre, F. Forget, Observed variations of methane on Mars unexplained by known atmospheric chemistry
and physics. Nature 460(7256), 720–723 (2009). doi:10.1038/nature08228
Applications of Electrified Dust and Dust Devil Electrodynamics. . .
343
J. Lowell, Constraints on contact charging of insulators: I. Spatial localisation of insulator states. J. Phys. D,
Appl. Phys. 19, 95–104 (1986)
J. Lowell, A.C. Rose-Innes, Contact electrification. Adv. Phys. 29(6), 947–1023 (1980)
R. Markson, The global circuit intensity: its measurement and variation over the last 50 years. Bull. Am.
Meteorol. Soc. 88(2), 223–241 (2007)
G.J. Marlton, R.G. Harrison, K.A. Nicoll, Atmospheric point discharge current measurements using
a temperature-compensated logarithmic current amplifier. Rev. Sci. Instrum. 84, 066103 (2013).
doi:10.1063/1.4810849
T.A. Mather, R.G. Harrison, Electrification of volcanic plumes. Surv. Geophys. 27(4), 387–432 (2006)
K.C. Mathpal, N.C. Varshneya, N. Dass, Precipitation-powered mechanisms of cloud electrification. Rev.
Geophys. 18, 361 (1980)
S. Matsusaka, H. Maruyama, T. Matsuyama, A. Ghadiri, Chem. Eng. Sci. 65, 5781 (2010)
T. Matsusyama, H. Yamamoto, Impact charging of particulate materials. Chem. Eng. Sci. 61, 2230–2238
(2006)
T. Matsuyama, H. Yamamoto, Charge relaxation process dominates contact charging of a particle in atmospheric conditions II. General model. J. Phys. D, Appl. Phys. 30, 2170–2175 (1997)
M.K. Mazumder, R.E. Ware, T. Yokoyama, B.J. Rubin, D. Kamp, Measurement of particle size and electrostatic charge distributions on toners using E-SPART analyzer. IEEE Trans. Ind. Appl. 27, 611–619
(1991)
L.S. McCarty, G.M. Whiteside, Electrostatic charging due to separation of ions at interfaces: contact electrification of ionic electrets. Angew. Chem., Int. Ed. Engl. 47(12), 2188–2207 (2008)
L.S. McCarty, A. Winkleman, G.M. Whiteside, Ionic electrets: electrostatic charging of surfaces by transferring mobile ions upon contact. J. Am. Chem. Soc. 129, 4075–4087 (2007)
O. Melnik, M. Parrot, Electrostatic discharge in Martian dust storms. J. Geophys. Res. 103, 29107–29118
(1998)
J. Merrison, J. Jensen, K. Kinch, R. Mugford, P. Nørnberg, The electrical properties of Mars analogue dust.
Planet. Space Sci. 52, 279–290 (2004)
J.P. Merrison, H.P. Gunnlaugsson, M.R. Hogg, M. Jensen, J.M. Lykke, M. Bo Madsen, M.B. Nielsen, P. Nørnberg, T.A. Ottosen, R.T. Pedersen, S. Pederse, A.V. Sørensen, Factors affecting the electrification of
wind-driven dust studied with laboratory simulations. Planet. Space Sci. 60, 328–335 (2012)
S.M. Metzger, M.R. Balme, M.C. Towner, B.J. Bos, T.J. Ringrose, M.R. Patel, In situ measurements of
particle load and transport in dust devils. Icarus (2011). doi:10.1016/j.icarus.2011.03.013
M. Michael, M. Barani, S.N. Tripathi, Numerical predictions of aerosol charging and electrical conductivity
of the lower atmosphere of Mars. Geophys. Res. Lett. 34, L044201 (2007). doi:10.1029/2006GL028434
M. Michael, S.N. Tripathi, S.K. Mishra, Dust charging and electrical conductivity in the day and nighttime
atmosphere of Mars. J. Geophys. Res. 113, E07010 (2008). doi:10.1029/2007JE003047
S.L. Miller, A production of amino acids under possible primitive Earth conditions. Science 117(3046), 528–
529 (1953). doi:10.1126/science.117.3046.528
A.A. Mills, Dust clouds and frictional generation of glow discharges. Nature 268, 614 (1977). doi:10.1038/
268614a0
H.A. Mizes, E.M. Conwell, D.P. Salamida, Direct observation of ion transfer in contact charging between a
metal and a polymer. Appl. Phys. Lett. 56, 1597–1599 (1990)
G.J. Molina-Cuberos, H. Lichtenegger, K. Schwingenschuh, J.J. Lopez-Moreno, R. Rodrigo, Ion-neutral
chemistry model of the lower ionosphere of Mars. J. Geophys. Res. 107(E5), 5027 (2002). doi:10.1029/
2000JE001447
G.J. Molina-Cuberos, R. Godard, J.J. López-Moreno, M. Hamelin, R. Grard, F. Simões, K. Schwingenschuh, V.J.G. Brown, P. Falkner, F. Ferri, I. Jernej, J.M. Jerónimo, R. Rodrigo, R. Trautner, M.J. Núñez,
N. Ibrahim, C. Groth, M. Fulchignoni, A new approach for estimating Titan’s electron conductivity
based on data from relaxation probe sensors on the Huygens experiment. Planet. Space Sci. 58(14–15),
1945–1952 (2010)
K.A. Nicoll, Measurements of atmospheric electricity aloft. Surv. Geophys. (2012). doi:10.1007/s10712-0129188
K.A. Nicoll, R.G. Harrison, Z. Ulanowski, Observations of Saharan dust electrification. Environ. Res. Lett.
6, 014501 (2011)
S. Nieh, T. Nguyen, Effects of humidity, conveying velocity, and particle size on electrostatic charges of glass
beads in a gaseous suspension flow. J. Electrost. 21, 99–114 (1988)
W.L. Nighan, Electron energy distribution and collision rates in electrically excited N2 , CO, and CO2 . Phys.
Rev. A 2, 1989–2000 (1970)
T. Oguchi, M. Tamatani, Contact electrification phenomena and powder surface treatments. Wear i68, 91–98
(1993)
344
R.G. Harrison et al.
M. Patzold, S. Tellmann, B. Hausler, D. Hinson, R. Schaa, G.I. Tyler, A sporadic third layer in the ionosphere
of Mars. Science 310(5749), 837-9 (2005)
T. Poppe, J. Blum, T. Henning, Experiments on collisional grain charging of micron-sized preplanetary dust.
Astrophys. J. 533, 472–480 (2000)
J.J. Qu, M.H. Yan, G.R. Dong, H.F. Zhang, R.P. Zu, W.Q. Tuo, A.G. Zhao, Z.H. Xiao, F. Li, B. Yang, Wind
tunnel simulation experiment and investigation on the electrification of sandstorms. Sci. China, Ser. D
47, 529–539 (2004)
K.R. Rasmussen, J.F. Kok, J.P. Merrison, Enhancement in wind-driven sand transport by electric fields.
Planet. Space Sci. 57(7), 804–808 (2009)
M. Ravichandran, A.K. Kamra, Spherical field meter to measure the electric field vector—measurements in
fair weather and inside a dust devil. Rev. Sci. Instrum. 70, 2140–2149 (1999). doi:10.1063/1.1149727
N.O. Renno, S.A. Rogacki, US Patent No. 8,536,879. Washington, DC: US Patent and Trademark Office
N.O. Renno, A.S. Wong, S.K. Atreya, I. de Pater, M. Roos-Serote, Electrical Discharges and broadband radio emission by Martian dust devils and dust storms. Geophys. Res. Lett. 30, 2140 (2003).
doi:10.1029/2003GL017879
N.O. Renno, V.J. Abreu, J. Koch, P.H. Smith, o.K. Hartogensis, H.A.R. De Bruin, D. Burose, G.T. Delory, W.M. Farrell, C.J. Watts, J. Garatuza, M. Parker, A. Carswell, MATADOR 2002: a pilot field
experiment on convective plumes and dust devils. J. Geophys. Res. 109, E07001 (2004). doi:10.1029/
2003JE002219
N.O. Renno, J.F. Kok, H. Kirkham, S. Rogacki, A miniature sensor for electrical field measurements in dusty
planetary atmospheres. J. Phys. Conf. Ser. 142, 012075 (2008)
A. Robledo-Martinez, H. Sobral, A. Ruiz-Meza, Electrical discharges as a possible source of methane on
Mars: lab simulation. Geophys. Res. Lett. 39, L17202 (2012)
W.A.D. Rudge, Atmospheric electrification during South African dust storms. Nature 91, 31–32 (1913).
doi:10.1038/091031a0
C. Ruf, N.O. Renno, J.F. Kok, E. Bandelier, M.J. Sander, S. Gross, L. Skjerve, B. Cantor, Emission of
non-thermal microwave radiation by a Martian dust storm. Geophys. Res. Lett. 36, L13202 (2009).
doi:10.1029/2009GL038715
M.J. Rycroft, R.G. Harrison, K.A. Nicoll, E.A. Mareev, An overview of Earth’s global electric circuit and
atmospheric conductivity. Space Sci. Rev. 137, 83–105 (2008). doi:10.1007/s11214-008-9368-6
M.J. Rycroft, K.A. Nicoll, K.L. Aplin, R.G. Harrison, Recent advances in global electric circuit coupling
between the space environment and the troposphere. J. Atmos. Sol.-Terr. Phys. 90–91, 198–211 (2012).
doi:10.1016/j.jastp.2012.03.015
D.S. Schmidt, R.A. Schmidt, J.D. Dent, Electrostatic force on saltating sand. J. Geophys. Res. 103(D8),
8997–9001 (1998)
B.F.J. Schonland, Atmospheric Electricity, (Methuen, New York, 1953)
E. Seran, M. Godefroy, N. Renno, H. Elliott, Variations of electric field and electric resistivity of air caused
by dust motion. J. Geophys. Res. Space Phys. 118, 5358–5368 (2013). doi:10.1002/jgra.50478
Y. Shirakawa, N. Ii, M. Yoshida, R. Takashima, A. Shimosaka, J. Hidaka, Quantum chemical calculation of
electron transfer at metal/polymer interfaces. Adv. Powder Technol. 45, 366–372 (2008)
A.A. Sickafoose, J.E. Colwell, M. Horanyi, S. Robertson, Experimental investigations on photoelectric and
triboelectric charging of dust. J. Geophys. Res. 106, 8343–8356 (2001)
F. Simões, J.J. Berthelier, M. Godefroy, S. Yahi, Observation and modeling of the Earth-ionosphere cavity
electromagnetic transverse resonance and variation of the D-region electron density near sunset. Geophys. Res. Lett. 36(14), L14816 (2009)
A. Sowinski, L. Miller, P. Mehrani, Investigation of electrostatic charge distribution in gas-solid fluidized
beds. Chem. Eng. Sci. 65, 2771–2781 (2010)
C.D. Stow, Dust and sand storm electrification. Weather 24, 134–140 (1969)
V. Sullivan, Variations of the vertical electric field and wind speed on days with airborne dust in Lubbock,
Texas, Master thesis, Texas Tech University, Lubbock, Texas (2013)
K. Tanoue, A. Ema, H. Masuda, Effect of material transfer and work hardening of metal surface on the current
generated by impact of particles. J. Chem. Eng. Jpn. 32, 544–548 (1999)
S. Trigwell, N. Grable, C.U. Yurteri, R. Sharma, M.K. Mazumder, Effects of surface properties on the tribocharging characteristics of polymer powder as applied to industrial processes. IEEE Trans. Ind. Appl.
39, 79 (2003)
S.N. Tripathi, M. Michael, R.G. Harrison, Profiles of ion and aerosol interactions in planetary atmospheres.
Space Sci. Rev. 137, 193–211 (2008). doi:10.1007/s11214-008-9367-7
Z. Ulanowski, J. Bailey, P.W. Lucas, J.H. Hough, E. Hirst, Alignment of atmospheric mineral dust due to
electric field. Atmos. Chem. Phys. 7, 6161–6173 (2007)
C. von Holstein-Rathlou, J.P. Merrison, C.F. Brædstrup, P. Nørnberg, The effects of electric fields on wind
driven particulate detachment. Icarus 220(1), 1–5 (2012)
Applications of Electrified Dust and Dust Devil Electrodynamics. . .
345
R.C. Whitten, I.G. Poppoff, J.S. Sims, The ionosphere of Mars below 80 km altitude. Planet. Space Sci. 17,
243–250 (1971)
J.A. Wiles, B.A. Grzybowski, A. Winkleman, G.M. Whitesides, A tool for studying contact electrification in
systems comprising metals and insulating polymers. Anal. Chem. 75, 4859–4867 (2003)
E. Williams, N. Nathou, E. Hicks, C. Pontikis, B. Russel, M. Miller, M.J. Bartholomew, The electrification
of dust-lofting gust fronts (‘haboobs’) the Sahel. Atmos. Res. 91, 292–298 (2009)
C.T.R. Wilson, Investigations on lightning discharges and on the electric field of thunderstorms. Philos. Trans.
R. Soc. Lond. A 221, 73–115 (1921)
L. Xie, K. Han, Influence of relative humidity on the aeolian electric field. Aeolian Res. 7, 45–50 (2012)
H.F. Zhang, T. Wang, J.J. Qu, M.H. Yan, An experimental and observational study on the electric effect of
sandstorms. Chin. J. Geophys. 47, 53–60 (2004)
X.J. Zheng, N. Huang, Y.H. Zhou, Laboratory measurement of electrification of wind-blown sands and simulation of its effect on sand saltation movement. J. Geophys. Res. 108, 4322 (2003)