Atmospheric Electrification in Dusty, Reactive Gases in the Solar System and Beyond

Surveys in Geophysics, Apr 2016

Detailed observations of the solar system planets reveal a wide variety of local atmospheric conditions. Astronomical observations have revealed a variety of extrasolar planets none of which resembles any of the solar system planets in full. Instead, the most massive amongst the extrasolar planets, the gas giants, appear very similar to the class of (young) brown dwarfs which are amongst the oldest objects in the Universe. Despite this diversity, solar system planets, extrasolar planets and brown dwarfs have broadly similar global temperatures between 300 and 2500 K. In consequence, clouds of different chemical species form in their atmospheres. While the details of these clouds differ, the fundamental physical processes are the same. Further to this, all these objects were observed to produce radio and X-ray emissions. While both kinds of radiation are well studied on Earth and to a lesser extent on the solar system planets, the occurrence of emissions that potentially originate from accelerated electrons on brown dwarfs, extrasolar planets and protoplanetary disks is not well understood yet. This paper offers an interdisciplinary view on electrification processes and their feedback on their hosting environment in meteorology, volcanology, planetology and research on extrasolar planets and planet formation.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

Atmospheric Electrification in Dusty, Reactive Gases in the Solar System and Beyond

Surv Geophys Atmospheric Electrification in Dusty, Reactive Gases in the Solar System and Beyond Christiane Helling 0 1 2 3 4 5 6 7 8 9 R. Giles Harrison 0 1 2 3 4 5 6 7 8 9 Farideh Honary 0 1 2 3 4 5 6 7 8 9 Declan A. Diver 0 1 2 3 4 5 6 7 8 9 Karen Aplin 0 1 2 3 4 5 6 7 8 9 Ian Dobbs-Dixon 0 1 2 3 4 5 6 7 8 9 Ute Ebert 0 1 2 3 4 5 6 7 8 9 Shu-ichiro Inutsuka 0 1 2 3 4 5 6 7 8 9 Francisco J. Gordillo-Vazquez 0 1 2 3 4 5 6 7 8 9 Stuart Littlefair 0 1 2 3 4 5 6 7 8 9 0 Department of Physics, Lancaster University , Lancaster LA1 4YB , UK 1 Department of Meteorology, The University of Reading , Reading RG6 6BB , UK 2 SUPA, School of Physics and Astronomy, University of St Andrews , North Haugh KY16 9SS , UK 3 Department of Physics and Astronomy, University of Sheffield , Sheffield S3 7RH , UK 4 Instituto de Astrof ́ısica de Andaluc ́ıa , P.O. Box 3004, 18080 Granada , Spain 5 Department of Physics, Nagoya University , Nagoya, Aichi 464-8602 , Japan 6 Centre for Mathematics and Computer Science , PO Box 94079, NL-1090 GB Amsterdam , The Netherlands 7 NYU Abu Dhabi , P.O. Box 129188, Abu Dhabi, UAE 8 Department of Physics, University of Oxford , Denys Wilkinson Building, Keble Road, Oxford OX1 3RH , UK 9 SUPA, School of Physics and Astronomy, University of Glasgow , Glasgow G12 8QQ , UK Detailed observations of the solar system planets reveal a wide variety of local atmospheric conditions. Astronomical observations have revealed a variety of extrasolar planets none of which resembles any of the solar system planets in full. Instead, the most massive amongst the extrasolar planets, the gas giants, appear very similar to the class of (young) brown dwarfs which are amongst the oldest objects in the Universe. Despite this diversity, solar system planets, extrasolar planets and brown dwarfs have broadly similar global temperatures between 300 and 2500 K. In consequence, clouds of different chemical species form in their atmospheres. While the details of these clouds differ, the fundamental physical processes are the same. Further to this, all these objects were observed to produce radio and X-ray emissions. While both kinds of radiation are well studied on Earth and to a lesser extent on the solar system planets, the occurrence of emissions that potentially originate from accelerated electrons on brown dwarfs, extrasolar planets and protoplanetary disks is not well understood yet. This paper offers an - interdisciplinary view on electrification processes and their feedback on their hosting environment in meteorology, volcanology, planetology and research on extrasolar planets and planet formation. Dust charging Discharging Solar system Extrasolar planets Electrification processes Electrical phenomena Moon 1 Introduction The Earth and the solar system planets were the only planetary objects known until the discovery of the first brown dwarf GD165B (Becklin and Zuckerman 1988) and the first extrasolar planet in 1992 (orbiting the pulsar PSR1257?12, Wolszczan and Frail (1992) ). Earth, Jupiter and Saturn are cloudy solar system planets for which atmospheric discharges in the form of lightning are confirmed observationally in radio and in optical wavelengths. Space exploration and ground-based observations have shown that lightning is a process universal in the solar system, but also that charge and discharge processes occur in a large diversity on solar system planets. Charging and discharging processes are essential for our understanding of the origin of our planet and maybe even for the origin of life: It is believed that charged dust is required to form planets and that lightning opens chemical paths to the formation of biomolecules. The purpose of this paper is to point out overlapping interests in electrifying media that contain liquid and solid particles in meteorology, volcanology, solar system objects, extrasolar planets, brown dwarfs and protoplanetary disks. We therefore provide a selective overview of atmospheric electrification processes and related electrical phenomena based on knowledge from solar system and Earth observations, and on laboratory-based research in combination with relevant findings and development in research on extrasolar planets, brown dwarfs and protoplanetary disks. We hope to stimulate a closer interaction between these communities. The last few decades have taken us from a Universe with only a single planetary system known, to one with thousands, and maybe millions, of such systems. We are now entering the time when we explore theories and results derived for the solar system and for Earth in application to unknown worlds. Figure 1 places Jupiter, one of the solar system giant gas planets, into the astrophysical context: Jupiter (right) is compared to the coolest stellar objects (M-dwarfs and brown dwarfs). Brown dwarfs bridge the stellar (represented by the Sun in Fig. 1) and the planetary regime as their atmospheres can be as cold as those of planets, but they form like stars. The Sun (left) is surrounded by hot plasma (corona), while planets are enveloped in a cold cloud-forming atmosphere some of which exhibit electrical phenomena as part of a global electric circuit. The Sun is intensively studied by satellites like SOHO1 and HINODE2 leading to efforts like SWIFF for space weather forecasting (Lapenta et al. 2013) . Comparable high-resolution monitoring is neither feasible for solar system planets, moons or comets nor for extrasolar objects. Instead, experimental work on Earth, Earth observation, modelling and comparative studies for the solar system and as an extrasolar objects need to be combined; examples for Earth studied as extrasolar planet are, e.g., given in Kitzmann et al. (2010), Be´tre´mieux and Kaltenegger (2013) and Hodosa´n et al. (2016). 1 2 Figure 2 compares images, spectra (disk-integrated radiation flux), atmospheric (Tgas, pgas)-structures and the local degrees of gas ionisation for Earth, Saturn and two types of brown dwarfs (L-type (pink)–hotter, and T-type (purple)–cooler). All data for Earth are from observations, the Saturn data are derived from Cassini3 spacecraft observation, the brown dwarf spectra are observed with SpeX on IRTF4 (Cushing et al. 2005) , and the (Tgas, pgas)—and the fe-structure are results from atmosphere simulations. fe refers to the local degree of ionisation and is defined as fe ¼ pe=pgas with pe and pgas the local electron and the local gas pressure, respectively. The Earth image is a photograph taken from the International Space Station. The Saturn image is a visible light image taken by the Cassini spacecraft, and the brown dwarf image is an artist’s impression based on atmosphere simulations. No direct image exists for any brown dwarf because the nearest brown dwarfs (the binary system Luhman 16) are 6.59 light years away from Earth. All three classes of objects have chemically and dynamically active atmospheres that form clouds and that may be undergoing local charge and discharge events. Their local atmospheric conditions differ, including the chemical composition, as result of their formation history and the irradiation received from a host star. Interdisciplinary research combining plasma physics, meteorology, volcanology, solar system exploration and astrophysics as suggested in F u¨llekrug et al. (2013) is required to study weather phenomena on Earth, solar system planets and on 3 4*spex/. Earth (NASA) Saturn (ESA) Brown Dwarfs radiation fluxes F (λ) |H--2-0--| CH4 |----| |-|H--C--2--H0----4-|| |C--i-A-||---CH---C---O-2--H|--|-4--| CH4 |---| L4 T4.5 1 1.5 2 Wave2le.5ngth (m3icrons)3.5 4 4.5 log(g)=3.0 solar metal icity log(g)=3.0 solar metal icity 104 TTeeff == 11840000 KK ((LT ddwwaarrff)) 500 1000 1500 2000 2500 3000 3500 Tgas[K] Tef = 1400 K (T dwarf) Tef = 1800 K (L dwarf) 10-18 10-16 10-14 10-12 10-10 10-8 10-6 fe atmospheric (Tgas, pgas)-structures: pe local degree of ionization fe = pgas 1.6e-14 1.4e-14 1.2e-14 )da (xFlbam 1e-14 lfu 8e-15 ve ilta 6e-15 eR 4e-15 2e-15 10-10 10-8 10-6 ]a 10-4 p [hsa 10-2 pg 100 102 10-10 10-8 10-6 ]a 10-4 p [hsa 10-2 pg 100 102 104 extrasolar planetary objects also in view of upcoming space missions like CHEOPS5, PLATO6 and JWST7. Plasma and discharge experiments are essential in providing a controlled environment in contrast to observation of atmospheric phenomena. Such experiments can involve the three different mass components constituting an atmospheric gas: electrons, ions and dust particles with their masses me \mion\md. The mass differences result in different spatial effects like ion acoustic waves and plasma crystals. An atmospheric environment that is only partially ionised may show plasma character on only local scales compared to the global scale of the comet, moon, planet, brown dwarfs or protoplanetary disk. One potentially far-reaching example for the origin of life on Earth is volcanoes (Johnson et al. 2008) which can produce significant electrostatic charging and subsequent lightning during eruption (Sect. 3.4), maybe also on Jupiter’s moon Io, for example. In volcanoes but also in terrestrial clouds, particles of similar mass govern the charge and discharge processes and plasmas form during violent discharges only. Understanding dust charging processes is important for space exploration because the local ionisation changes as a result of the variability of the solar wind hitting the moon’s or an asteroid’s surface. A spacecraft landing, like Philae, the Rosetta lander, has a very similar effect (Sect. 4). In situ measurements from the chemically active Earth’s atmosphere offer insight into charge and discharge processes, their local properties and their global changes (Sect. 3.1). While plasma experiments are conducted in a controlled laboratory environment, measurements inside the uncontrollable Earth’s natural atmospheric environment lead to an understanding of the vertical and horizontal ionisation where the relative importance of electrons, ions and dust, hence their total mass relation, changes with atmospheric height. For example, the fair weather current is carried by ions only due to the lack of free electrons between 0 and 60 km. Understanding the Wilson Global circuit (Sect. 3.3) helps the understanding of the Earth’s weather and climate. Such observations allow an understanding of atmospheric processes on Earth that can only be gained for solar system and extrasolar bodies by intensive modelling efforts guided by observations and experiments. Section 2 provides a short background summary on charge processes of discrete solid or liquid surfaces in atmospheric gases, the link to laboratory works and an example of related plasma technology development. Section 2 further sets the stage for this interdisciplinary paper by defining terms used in later sections. Section 3 summarises charging and discharging processes in the terrestrial atmosphere, including processes in the atmospheres of other solar system planets. Section 4 reviews charging processes on the Moon and asteroids in the presence of solar wind and space plasmas, but without substantial neutral atmospheres. Section 5 provides insight into astronomical observations that suggest that mineral-cloud forming atmospheres of brown dwarfs and extrasolar planets are also electrically active, that different ionisation processes will electrically activate different parts of such atmospheres, and that similar processes are expected to act in protoplanetary disks. Section 6 concludes this paper. Each section ends with a list of future works/ open questions where suitable. 5 6 7 2 Setting the Stage for Interdisciplinary Exchange This section outlines the key concept of this interdisciplinary paper, and it provides definitions of terms used in Sects. 3–5. This section links to laboratory experiments which have driven the understanding of ionised atmosphere gases that contain or form dust particles or liquid droplets. One example of plasma technology development is included to demonstrate the impact of this paper’s theme also beyond academic research. This section deals with the smallest scales where charge processes act, and later sections will address topics related to successively larger-scale charge processes in the terrestrial atmosphere, on the Moon and asteroids, and also outside the solar system in extrasolar planets, brown dwarfs and protoplanetary disks. 2.1 Fundamental Charging Processes The key concepts in this paper depend on the accumulation and dissipation of electrical charge on discrete solid or liquid surfaces suspended in atmospheric gases. The free charge on the surfaces can arise from two primary mechanisms (in the planetary atmosphere context): processes involving (1) friction (triboelectric charging); and (2) the transport of free charge (plasma processes). More details on processes specific to various environments like Earth’s atmosphere, volcanoes or extrasolar planets are provided in the respective subsections (e.g., Sects. 3.1, 3.2 and 3.4). 2.1.1 Classical Frictional Charging Transiently contacting surfaces can lead to charge accumulation, by producing either a surplus or a deficit of electrons compared to the neutral case. Indeed, there is evidence that fragments of polymer chains can be exchanged by colliding particles (Saunders 2008), leaving net charges on the surfaces. This process is termed triboelectric charging and has a very long history of practical application (Galembeck et al. 2014), even if the underlying processes are still not entirely resolved. Originally, contact electrification was used to refer to electrostatic charge transfer resulting from contact, including contact modes such as detachment, sliding, rolling, impact. The specific charge processes related to rubbing were only later termed as triboelectrification. Such charging is an inevitable consequence of the frictional interaction between hard surfaces: electrons transfer (by some process) from one surface to the other, leading to charged surfaces. For example, dust entrained in strong, collisional flows (such as volcanic eruptions or mineral clouds in extrasolar planets, Sects. 3.4 and 5) will acquire charges of different polarity (negative and positive) directly from the intergrain collisions themselves. Such macroscopic particles can include ice crystals in atmospheric clouds, where the diversity of growth rates (and consequent dynamics) of crystals influences the polarity of charge transfer and leads to such clouds becoming charge separated by the relative drift of the charged particles (Saunders 2008). Charge accumulation and separation can lead to energetic relaxation, in the form of lightning. 2.1.2 Plasma Charging There is an additional mechanism for forcing charge onto a surface, in possibly much larger quantities than can be acquired by triboelectric or contact processes: plasma charging. A plasma is a gas in which a fraction of the molecules are ionised, leading to an abundance of free charge existing as an additional ‘‘gas’’ component. Though neutral overall, there is a natural scale-length over which the plasma can create large potential differences caused by charge population fluctuations: this is because free electrons are light and mobile compared to the heavier positive ions, and therefore the electrons can temporarily escape their charged counterparts, leading to charge densities appearing for short intervals, and over restricted distances (this is explained in detail in subsequent sections below). Should an isolated solid (dust or crystal) or liquid (aerosol) surface be introduced into this plasma, these natural fluctuations in the charge distribution will cause such surfaces to acquire surplus free charge, forced onto it by the action of the plasma itself. Isolated surfaces exposed to plasma will quickly (typically on a microsecond timescale or less) charge up to reach the plasma or floating potential (Khrapak et al. 2012; Khrapak and Morfill 2008; Hutchinson and Patacchini 2007) , by the action of a continuous electron current to the surface from the ambient plasma, which rapidly establishes a negative charge before the compensating positive ion current can respond. Ultimately, there is a balance reached, but one that reflects the relative electron mobility over the ions. Since there is so much more free charge available in a plasma compared to triboelectric processes, there is an enhanced capacity for dust exposed to plasma discharges to store considerable surface charge in comparison with purely collisional interactions between grains: since the plasma surface charge reflects the plasma conditions, and not just the grain chemistry and collisionality, then the plasma is an independent and effective agent for creating charged particles. 2.1.3 Defining General Terms After a summary of the principal mechanisms for charging surfaces in gases in Sects. 2.1.1 and 2.1.2, the most important vocabulary used throughout the paper is defined below to allow a better understanding of the links between the interdisciplinary topics in Sects. 3–5. The Appendix provides an glossary. Dust particles, aerosols, droplets. An important feature of many charging processes is the presence of macroscopic particles such as dust, aerosols or droplets. These are macroscopic particles large enough to move under the influence of gravity. The particle sizes can vary by orders of magnitude. They can be liquid or solid. They can be composed of a mix of different materials that changes with temperature. Aerosols are suspended particles of either phase. Dust is predominant on the Moon and asteroids, in volcanic lightning and mineral clouds of extrasolar planets and brown dwarfs, and as building blocks for planets in protoplanetary disks. Also hydrometeors (droplets, graupel and ice particles, snowflakes . . .) could fall into this category, but are considered aerosols in geoscience. Macroscopic particles such as dust and aerosols can be electrically charged which de-mobilises the charge that previously resided in the gas in the form of electrons or ions. Dust, for example, will acquire a negative total charge in the absence of external influence like stellar UV radiation. Ionisation is the process of dissociating neutrals into charged species, due to a variety of mechanisms: electron impact ionisation, Penning ionisation (ionisation through chemical reactions), direct dissociation by strong electric fields, UV-photoionisation. The total electric charge is conserved during ionisation, but once the charges are free they can move independently. In air (the atmospheric gas on Earth with its electronegative oxygen component) free electrons are very short lived in the absence of strong electric fields. Ionised air in the Earth’s troposphere and stratosphere consists of positive and negative ions. The fair weather currents on Earth are ion currents (see Sect. 3.3). Plasma is a gas consisting of charged particles. It is often restricted to charged particle gases where collective phenomena, like plasma oscillations, are more important than collisional phenomena. A plasma is created if there is sufficient ionisation of neutrals that the charged particle density becomes significant. A plasma is characterised by the capacity to produce a collective self-field that is significant when compared to any imposed field (such as that produced by external electrodes, or induced by collapsing magnetic fields, or by impinging electromagnetic radiation). An electrically neutral medium is created that can respond to an external electromagnetic field, but there is no spontaneous charge separation in equilibrium on scale-lengths greater than the Debye length.8 There is a significant distinction between plasmas which are collisionless, and those which are collisional9: (1) Collisionless plasmas consist mainly of positively charged ions and of electrons or negatively charged ions, depending on the electronegativity of the ionised gas. They interact through electromagnetic fields rather than through mechanical collisions. Examples are the magnetosphere and the interplanetary plasma (Sect. 4) where the assumption of ideal magnetohydrodynamics (MHD) holds. (2) In a collision dominated plasma, the motion of charged particles is dominated by collisions with neutral atoms and molecules, rather than by the direct electromagnetic interaction with other charged particles. The transiently existing plasmas in the terrestrial tropo-, strato- and mesosphere up to the E layer of the ionosphere are mostly collision dominated plasmas, except for the highly ionised and hot lightning return stroke channel. Charging or Charge Separation will be used for the process where macroscopic particles like dust or aerosols are charged. This can occur in particle collisions (in thundercloud electrification, dust devils in deserts, volcanic lightning) in non-ionised atmospheres or in vacuum or by attaining charge from a plasma (e.g., in dusty plasmas) spontaneously due to the different mobility of the charged species, in ambipolar diffusion, for example. If mechanical forces (gravity, convection) that act on the charged dust particles are stronger than the electric forces, charges can be separated over a certain distance. An electric potential builds up that can discharge by lightning and the related transient luminous events. Electrification is understood as the processes leading to charging of dust or other macroscopic particles obeying both polarity and charge conservation. As a result, a macroscopic electric field can build up. Electrification is sometimes used synonymously with Charging or Charge separation. Discharging is the process where the electric potential is released by electric currents. This can happen continuously or through a rapid transition like the rapid growth of discharge channels in lightning discharges. Emission of high-energy radiation can be associated with the rapid channel growth. 8 The Debye length is the length beyond which the Coulomb force of a charge cannot affect other charges. Strictly, the Debye length is the e-folding distance within which charge neutrality is not guaranteed, because thermal fluctuations can displace electrons relative to positive ions, leaving a small net charge. 9 These terms refer to approximations made in the plasma kinetic gas theory where the Boltzmann equation describes the evolution of the particle distribution function f ðx; v; tÞ. Neglecting the collisional source term of the Boltzmann equation leads to the collisionless Boltzmann equation (Vlasov equation) from which then the MHD equations are derived, and the electric and magnetic field strengths are derived as macroscopic quantities. In a collisional plasma, the full Boltzman equation is to be solved. 2.2 Charged Dust in Experimental Work Dust in plasmas has a long history—one which is even more relevant in contemporary planetary exploration. This section explores the phenomena associated with dust interacting with ionisation in the ambient atmosphere to ensure non-equilibrium processes (both physics and chemistry) have a significant and enduring influence on the evolution of the atmosphere in general, including the dust itself. The discussion here ranges over the impact of charged dust imposing a long-range order in confined plasmas, through to microdischarges arising from binary encounters between freely floating charged aerosols, releasing low-energy free electrons into the ambient atmosphere, with all the possibilities that this entails for molecular activation by dissociative attachment and radical formation. The common theme throughout is the capacity—literally— for dust to retain the electrostatic memory of ambient discharges via free-charge acquisition and for that discharge legacy to be reshaped and realised in potent form by harnessing hydrodynamical forces on fluid timescales, rather than plasma ones. In this way, transient plasma effects can be stored, reconfigured and released on meaningful scales in such a way as to have a tangible influence on large-scale evolution of planetary atmospheres. The following sections discuss dust–plasma interactions in Sect. 2.2.1. laboratory plasma dust, where floating particulates can be a help or a hazard in plasma applications, including plasma crystals, and in Sect. 2.2.2. the dynamic evolution of charged aerosols, where fluid deformation and evaporation can moderate the evolution of encapsulated targets. 2.2.1 The Plasma Laboratory: Dusty Plasmas and Plasma Crystals Dusty plasmas have been studied in laboratory experiments for several decades. Langmuir et al. (1924) reported the observation of minute solid particles and aggregates in a laboratory streamer discharge and suggested the dust could play a role in ball lightning (see also Rakov and Uman 2003 for a review) . ‘‘Dusty plasmas’’ are sometimes referred to as ‘‘complex plasmas’’ although the latter description is more wide-ranging and can include other types of constituents and features such as sheaths (Phelps and Allen 1976) , quantum effects and dust. Dusty plasma is referred to in cases when collective behaviour of dust becomes important resulting in new types of waves and instabilities. This occurs when the Debye length and interparticle distance are of the same order and the effects of neighbouring particles cannot be neglected, as opposed to the case when the Debye length is much less than the typical interparticle distance (isolated charged dust). The experimental research on dusty plasmas in laboratories has (1) been aimed at increasing fundamental understanding and (2) also been strongly motivated by the need to control the behaviour of dust in plasmas that are used in industrial applications. Dust deposited from within the plasmas that are involved in the semiconductor component fabrication and materials processing industries can damage the components and significantly affect the productivity of these industries. In contrast to the need to mitigate the potentially harmful effects of dust in industrial plasma etching and deposition, the capability to form and control dust in plasmas is being exploited in the production of nanoparticles for the expanding nanoscience industry. Fundamental research programmes have explored phenomena such as dust crystallisation and wave propagation within dusty laboratory plasmas where a stationary and fully ionised gas is considered. In laboratory experiments, the Earth’s gravitational field influences the dusty plasma behaviour and, while the vast majority of experiments have been carried out in laboratories on the surface of the Earth, there have been some experiments on dusty plasmas carried out in the near-weightless conditions within the International Space Station. Whereas at sea level 2D dust crystals can be produced, the low-gravity conditions are usually needed to produce 3D dust crystals. Several types of waves, including longitudinal electron plasma waves and ion acoustic waves (Allen and Phelps 1977) , can propagate in dust-free plasmas formed from ionised gas and containing electrons and ions as well as some neutral atoms and molecules. Additional wave propagation modes appear if a magnetic field is applied to the plasma. While all of these waves are usually damped as they propagate, it is also possible for them to become growing waves, or instabilities (Allen and Phelps 1977; Kuhn et al. 1981) , when appropriately excited. For example, ion acoustic waves (Allen and Phelps 1977) can be driven unstable by passing a current through the plasma, i.e. they are triggered by a drift motion of the electrons relative to the ions. In a dusty plasma, the charged, massive dust particles can produce new types of wave motion: The dust-ion acoustic wave (DIAW) is a modified ion acoustic wave, where the ions continue to provide the inertia and the presence of the quasi-stationary charged dust particles modifies the normal ion acoustic wave dispersion. In contrast to the DIAW, in the dust acoustic wave (DAW) the dust particles move and provide the inertia rather than the ions. Both the DIAW and the DAW can be observed because their frequencies are low enough for camera systems to resolve the images of the wave propagation. Measurement of dusty plasmas in the laboratory and comparison with simulations using particle in cell (PIC) codes allows these codes to be benchmarked against the laboratory experimental observations. PIC code simulation of laboratory plasma experiments and comparison with space measurements has proven successful in the case of auroral kilometric radiation (Speirs et al. 2008; McConville et al. 2008) because of their capability to simulate the onset and dynamics of microinstabilities in dusty plasmas. The use of PIC codes to simulate the behaviour of dusty plasmas in space should prove equally fruitful in obtaining detailed explanations of the formation, properties and consequences in astrophysics (Shukla and Mamun 2002; Fortov and Morfill 2010). 2.2.2 Delivering Charges to Microscopic Particles The evolutionary processes governing the dynamics and stability of charged macroscopic water droplets in a discharge plasma are part of an innovative collaborative project on bacteria detection (Rutherford et al. 2014; Maguire et al. 2015) . The technique of using droplet evaporation as a moderator for charge deposition provides a method to precisely deliver a known amount of charge to microscopic particles such as bacteria cells or (cloud) condensation seeds. For that, aerosolised bacteria samples will be passed through a discharge plasma to acquire significant electrical charge which can be measured in the laboratory. If the charge-carrying aerosol evaporates, its surface area decreases, but the aerosol retains the charge. Ultimately, if the Coulomb force overcomes the surface tension, then the droplet expels charge to bring the retained charge back into the stability limit (the Rayleigh limit QrðtÞ), which is a function of its radius. Hence the droplet continues to track the Rayleigh limit10 as it evaporates. Once all the fluid has gone and the interior seed 10 The Rayleigh limit, QrðtÞ, gives the limiting size of the surface electric field that balances the surface tension: the latter provides the restoring force to return the droplet to its equilibrium spherical shape and so causes the perturbed droplet to oscillate. If the distorted outer surface of the droplet carries sufficient electric charge, then the local surface field may oppose the effect of surface tension and thus prolong the restoration particle (bacterium or grain) is revealed, the charge placed on it is known. This is the charge consistent with the Rayleigh limit at the radius of the grain. The charging mechanism can be described as follows (Maguire et al. 2015). Water droplets entering a plasma will form a sheath between the droplet surface and the plasma, as a simple consequence of the disparity in mobility between electrons and ions. Electrons will collide more frequently with the drop surface and remain there, causing it to acquire a negative surface charge. The charged droplet will then attract positive ions from the plasma until the electron and ion currents to the surface of the droplet reach equilibrium; at this point, the droplet is at the plasma potential. Suppose an initially stable water droplet has acquired charge by passing through a plasma (or indeed by an alternative charging mechanism; green vertical line in Fig. 3) and is now floating freely in air, having left the plasma behind. If the initial droplet charge is less than the initial Rayleigh limit, Qr0 , of the droplet, then the droplet is stable. As evaporation proceeds outside the plasma, the droplet charge stays roughly constant, while the Rayleigh limit, QrðtÞ, evolves according to QrðtÞ ¼ bðtÞQr0 ; with Qðt ¼ 0Þ ¼ a Qr0 ; a\1 being the initial charge on the droplet, and bðtÞ\1 for all t [ 0. The initial values for the results in Fig. 3 are: aðt ¼ 0Þ ¼ 0:0025, r0 ¼ rðt ¼ 0Þ ¼ 10 lm, Qðt ¼ 0Þ ¼ 104 electronic charges (e) because Rayleigh limit is 4 106 e. b ¼ 1 at t ¼ 0; b is not shown in Fig 3. If QrðtÞ decreases far enough that QrðtÞ QðtÞ, then the droplet will become unstable and emit sufficient charge to restore the stability condition of QrðtÞ [ QðtÞ. Evaporation continues until once again the stability condition is broken and more charge is emitted back into the ambient gas. This feedback loop continues until the entire droplet has evaporated. As the droplet evaporates, both the droplet radius r(t) and the Rayleigh limit for the charges on the droplet, QrðtÞ, decrease. If the droplet encapsulates a bacteria or dust grain, the evaporation cannot proceed beyond a minimum radius rm. The final charge on the droplet of size rm at a final time, tf , is then Qðtf Þ bðtqfÞ ffiffiffiffiffiffiffiffi¼ffiffiffi Qrðtf Þ Qr0 8p ce0rm3: ð1Þ ð2Þ The upper limit of final droplet charge depends only on the minimum radius of the particle, c surface tension of the droplet, rm, left behind once the droplet has evaporated, irrespective of the starting charge. This is assuming that the Rayleigh limit is encountered at some intermediate point in the evaporative evolution of the water mantle that forms the drop encapsulating a bacteria or dust grain. This is a valuable process, since grains processed in this way carry the electrostatic legacy of the plasma environment encountered earlier in their history. Such charged particles can either act as a source of low-energy free charge injected into the atmosphere to produce non-equilibrium electron-moderated chemical evolution of the latter (for example, Footnote 10 continued to equilibrium profile, i.e. reduce the oscillation frequency. If there is sufficient surface charge, then the deformation persists, and the oscillation frequency is formally zero which defines the Rayleigh limit. Exceeding the Rayleigh limit means that the droplet is unstable to perturbation and is forced to eject charge and mass. dissociative attachment producing radicals) or indeed a constraining electrostatic environment stable over fluid length and time scales. 3 Electrification and Discharging in Terrestrial and Planetary Atmospheres When we aim to understand electrification and electric phenomena in weakly ionised atmospheres of extrasolar planets, a characterisation of the phenomena on Earth and in the atmospheres of solar system planets can provide guidelines and inspiration. This section therefore starts with an overview of the main electrical processes in the terrestrial atmosphere up to the ionosphere, the fair weather currents and the thunderstorms with transient luminous events and terrestrial gamma-ray flashes. Then we continue with lightning phenomena in volcanic ash plumes and review a few processes in the atmospheres of other solar system planets. For more details see Rakov and Uman (2003) , Leblanc et al. (2008), Dwyer and Uman (2014) , Betz et al. (2009), F u¨llekrug et al. (2006), Ebert and Sentman (2008) . Ionisation and electric currents in the terrestrial atmosphere are driven by two main mechanisms: (a) The atmosphere is very weakly ionised by external sources like Cosmic Rays and radioactivity (Sect. 3.1). The resulting conductivity supports the fair weather currents that relax electric potentials in atmospheric regions far from thunderstorms. (b) Thunderclouds play a particular role in separating electric charges and in building up large electric potentials (Sect. 3.2). Cloud particles first exchange charge during collisions and are then separated due to mechanical forces (such as gravity and convection) larger than the attractive electric forces between particles of opposite polarity. For this reason, meteorologists use lightning flashes as indicators for strong turbulent convection in the atmosphere. When these electric potentials suddenly discharge, a variety of ionised and conducting channels is formed through localised ionisation processes (collisional, thermally driven or photon impact). In the first stage of a discharge, these ionisation reactions are driven by strong electric fields and local field enhancements and are dominated by the impact of fast electrons on neutral atoms or molecules, while at later stages Ohmic heating and thermal equilibration create temperature driven ionisation reactions. 3.1 Ionisation of the Terrestrial Atmosphere Outside Thunderstorm Regions In common with other solar system atmospheres (Harrison et al. 2008), the Earth’s lower atmosphere outside thunderstorm regions is made electrically conductive by the ionising action of high-energy charged particles generated within the heliosphere (e.g., solar energetic particles, SEPs) and beyond (e.g., galactic cosmic rays, GCRs). A consequence of the terrestrial atmosphere’s small but finite conductivity ( 10 14 S m 1 in surface air, see also Fig. 6) is that current flows can occur through the atmosphere, between disturbed weather and fair weather regions. Similar circumstances occur in other atmospheres, depending on the existence of charge separation processes and the atmospheric conductivity. Ion production in the Earth’s lower atmosphere (i.e. the troposphere and stratosphere) results from a combination of terrestrial and extraterrestrial sources. Near the planet’s continental surfaces, the effects of natural radioactivity contained within the soil and rocks, or released in the form of radioactive gases such as radon, provide the dominant source of ion production. At heights from 3 to 5 km above, the continents (i.e. above the boundary layer where eddy diffusion of radon isotopes occurs which depend on orography), or over the oceans, extraterrestrial sources, principally GCRs dominate the ion production, while SEPs and UV irradiation dominate the ionisation in the ionosphere, but typically do not have sufficient energy to reach the troposphere. Balloon-borne Measurements Vertical soundings of the ion production rate in the troposphere and stratosphere (i.e. to about 35 km) can be made using balloon-carried instruments11. Historically, this was the original airborne platform through which the existence of the cosmic source of ionisation was confirmed, in a manned balloon flight made by Victor Hess on 7 August 1912 (Hess 1912). This flight carried ionisation chambers and fibre electrometers, in which the rate of decay of the charged fibre was recorded visually and the ion production rate inferred (Pfotzer 1972) . Hess found that the ion production rate initially diminished with height, but then began to increase (Fig. 4A). 11 The atmosphere above this altitude is sometimes called ignorosphere, because above balloon and below satellite altitudes it is very difficult to explore. In particular, the density of free electrons in the lower ionosphere can now be measured only indirectly through the pattern of electromagnetic radiation that is emitted by lightning strokes and reflected by the ionosphere (Lay et al. 2010; Shao et al. 2013) . This subsequent increase indicated that ionisation was originating from above. Figure 4B shows a profile of the ion production rate per unit volume at standard temperature and pressure, qSTP, made using a modern balloon-carried Geiger counter (or Geigersonde) launched from a midlatitude site (details are given in Harrison et al. 2014). This shows the same increase in ionisation observed by Hess at the lower altitudes, but the modern balloons extend the measurements to greater altitudes. A characteristic feature is the maximum in ionisation at about 20km, first observed by Regener and Pfotzer (1935) . The presence of the Regener–Pfotzer maximum results from a balance between the energy of the incoming particles, and the density of the atmosphere. A long series of regular Geigersonde measurements has been made by the Lebedev Institute in Moscow, using a variety of sites including Moscow, Murmansk and Mirny (Antarctica). The value of this stable long-term measurement series is considerable, as, by taking advantage of the different geomagnetic latitudes of the sites concerned, it allows features of the cosmic ray ionisation to be established. Cosmic rays follow the geomagnetic field lines, and the lower energy particles are able to enter at higher latitudes (which is expressed as a lower geomagnetic rigidity). The high-energy CR particles survive for longer in the Earth’s atmosphere, while the low-energy CR particles are completely absorbed soon after they enter the atmosphere. Figure 5 shows a long times series of Geigersonde measurements made at the Regener–Pfotzer maximum, from sites with different rigidity (Stozhkov et al. 2013) . The 11-year (Schwabe) cycle in solar activity is clearly present through the inverse response in GCRs, and, at the high-latitude sites, the exceptional nature of the cosmic ray maximum in 2010/11 associated with the deep solar minimum is particularly apparent. Atmospheric Conductivity Cosmic ray ionisation in the terrestrial atmosphere sustains a steady source of cluster ions, which provide the finite conductivity of air. The total conductivity, rt, is given by rt ¼ e ðlþnþ þ l n Þ ð3Þ where l represents the mean mobility of positive or negative ions present, n the associated bipolar ion number concentrations and e is the elementary charge. Ions are removed by attachment to aerosol particles and water droplets, reducing the conductivity in these regions. Both the mobility and concentration vary with atmospheric properties and composition. The mobility of ions depends on the environmental temperature and pressure, and the ion concentration is strongly affected by attachment to aerosol particles and water droplets, reducing the conductivity accordingly where the aerosols are abundant. This means that, in the Earth’s environment, where aerosols are generated both naturally and through human activities, the local air conductivity can show an anthropogenic influence (Harrison 2006; Silva et al. 2014), allowing early indirect conductivity measurements to provide an insight into historical air pollution (Harrison 2006; Aplin 2012) . Together with variations in the source rate, qSTP, these lead to a variation in the conductivity with height (e.g., Harrison and Carslaw 2003). At the heights of the lower ionosphere, where photoionisation also contributes appreciably, the conductivity becomes substantially larger than in the lower atmosphere. Figure 6 shows a vertical profile of the air’s conductivity, and a calculation of the relaxation timescale, defined by 0=rt. This is the e-folding timescale for the discharge of an isolated particle in a conductive medium. This provides an indication of how active (in terms of the rate of charge separation) a charging process needs to be at different heights in the atmosphere. In comparison with lower troposphere air with a typical conductivity of 10 14 S m 1 as reviewed by Rycroft et al. (2008), the Fig. 4 Vertical profile of the ionisation rate in the terrestrial atmosphere, as (A) originally obtained by Hess (7th August 1912), with ionisation at each height shown relative to the measured surface ionisation and, (B) from a series of balloon flights (colours used to identify individual flights) made from Reading, UK, during 2013. qSTP ½cm 3 s 1 the ion production rate per unit volume, for air at standard temperature and pressure (STP) A planetary surface has a greater electrical conductivity, of at least 10 8 S m 1. This means the air represents a low-conductivity region sandwiched between upper and lower boundaries having much greater conductivity. 3.2 Thundercloud Electrification, Lightning and Transient Luminous Events Ionic Conductivity and Ionic Plasmas in the Terrestrial Atmosphere Most electric phenomena in the terrestrial atmosphere are carried by ions and aerosols; only in the strong transient electric fields of an evolving discharge or in the ionosphere are more electrons free and not attached to electronegative atoms, molecules or larger compounds consisting, e.g., of water molecules clustering around ions, other aerosols, up to droplets from microto millimetre size. Cosmic rays and radioactivity are external sources of ionisation (Sect. 3.1); they first create electron ion pairs, and then the electrons rapidly attach to electronegative molecules (mostly to oxygen) leaving the positive and negative ions in the atmosphere behind which carry the fair weather currents (Sect. 3.3). The Electric Field in Thunderclouds builds up in two stages. In the first stage, macroscopic particles are electrically charged, and in the second stage particles of different polarity are separated by gravitation or other (mechanical) forces; in order to separate particles with different polarities, these forces need to be stronger than the electric attraction between charges of different polarity, since otherwise the electric forces would counteract the growth of the electric field. The possible charging mechanisms at work within normal terrestrial thunderclouds are reviewed, e.g., by Jayaratne et al. (1983) and Saunders (2008). An important conclusion of these reviews is that charge is efficiently separated between particles only in direct collisions. Liquid droplets cannot experience collisions or fracture as a charging process as they would typically merge on contact, and hence they do not charge easily. However, frozen particles can collide and exchange charge. Therefore, terrestrial water clouds get electrified mostly in regions below the freezing temperature (Mason 1953), more precisely at temperatures between 0 and 40 C. The dominant charging mechanism is thought to occur when graupel and ice particles collide. Saunders (2008) reviews the evidence from Krehbiel’s (1986) cloud measurements in 1986 ‘‘that ice crystals rebounding from riming graupel12 in the presence of super-cooled water is a requirement of the charge transfer process’’. This observation is consistent with laboratory measurements of Saunders (2008) that during collision essentially ‘‘fast growing ice surfaces charge positively, and conversely, sublimating (graupel) surfaces charge negatively’’. However, further dependencies on growth velocities, etc. need to be taken into account. The particle collisions are mediated by gravity acting on large particles and by turbulent convection within the cloud. Gravity will also move the heavy positively charged graupel particles downward, while the light positive ice crystals move upward with the convective flow of the cloud air, creating charge centres and electric fields within the cloud. This particular charging mechanism is based on the intrinsic polarisation of water molecules. Macroscopic particles of different material can charge quite efficiently, too, and create electric fields and discharges. Both volcanic ash plumes, socalled dust devils in terrestrial deserts and various granular media in the laboratory, support discharges, as is discussed further in Sect. 3.4. The understanding of charging processes in volcanic ash plumes might inspire further progress on the long-standing question of charging normal thunderclouds (Yair 2008) . Such normal water clouds mixed with dust have recently been observed to exhibit particularly strong and exceptional discharges (Fu¨ llekrug et al. 2013) . Due to the attachment of ions to water droplets, electric charges in clouds are particularly immobile. The conductivity in the remaining gas phase is therefore low before lightning activity starts. This low conductivity (hence low degree of ionisation, see also Fig. 2) supports a high electric field up to the moment of discharging. The Stages of Lightning Lightning is the sudden release of the electric potential energy through the fast growth of a disperse network of ionised channels. On average, 44 5 lightning flashes (intracloud and cloud-to-ground combined) occur around the globe every second (Christian et al. 2003) . Moreover, according to OTD (Optical Transient Detector) measurements, lightning occurs mainly over land areas with an average land/ocean ratio of approximately 10:1 (Christian et al. 2003) . The visible growing channels are called lightning leaders; their path is prepared by streamer coronae. While streamers are spacecharge-driven ionisation fronts, leaders maintain their internal conductivity by increased temperature, molecular excitations and ionszation reactions in the discharge channel. If a conducting channel connects cloud and ground, the return stroke carries the largest current and is visible and audible as the lightning stroke; but intra- and intercloud lightning are much more likely. The stages of lightning have been described in many articles, with varying emphasis on phenomena or physical mechanisms. A few recent ones are by Bazelyan et al. (2009), Cooray (2003) , Rakov and Uman (2003) , Betz et al. (2009), Dwyer and Uman (2014) and Cooray (2015) . A long-standing question is how lightning can be initiated because the observed electric fields are below the classical breakdown field (where electron impact ionisation overcomes electron attachment to oxygen in the Earth’s atmosphere; e.g., Treumann et al. 2008; Helling et al. 2013) , and free electrons are not available anywhere in the atmosphere. Gurevich et al. (1992) suggested that cosmic particle showers could supply free electrons and that relativistic run-away electron avalanches could develop in an electric field below the classical breakdown value. Gurevich and Karashtin (2013) recently suggested that the 12 Riming graupel is a graupel particle coated with water droplets that froze immediately when they collided with the ice surface of the graupel. The surface structure of graupel deviates from a perfect crystalline structure (e.g., Blohn et al. 2009). interplay of a cloud particle with Cosmic Rays could start the discharge. A quantitative analysis confirming this scenario is presented by Dubinova et al. (2015) . Lightning occurs not only between cloud and ground, but also within and between clouds. Also the ‘‘bolt from the blue’’ is a phenomenon where a lightning strike seems to appear out of a blue sky next to a thundercloud. These strikes are an indication that lightning leaders can leave the cloud also at its upper edge or in the sideward direction and then turn downwards. Transient Luminous Events The full-scale discharge activity associated with terrestrial water clouds became known in the scientific literature only after 1989 when the first Transient Luminous Events were described (for article collections, see F u¨llekrug et al. (2006), Ebert and Sentman (2008) ). Basically, electric potential stored in a cloud can also discharge in the upward direction as a jet up into the stratosphere or as a gigantic jet that extends into the mesosphere. The primary lightning can drive secondary discharges, namely elves, halos and sprites in the E layer of the ionosphere and in the night-time mesosphere (where the D layer of the ionosphere is located during day time)13. Elves and halos are responses of the lower edge of the ionospheric E layer to the electromagnetic pulse and the quasi-static potential of the parent lightning stroke, while sprites propagate downward from the ionosphere into the mesosphere (so-called column sprites) and sometimes back up again (carrot sprites; Stenbaek-Nielsen and McHarg 2008, Luque and Ebert 2009) . Due to similarity relations between discharges at different atmospheric densities (Pasko 2007; Ebert et al. 2010) , tens of kilometres-long sprite discharge channels in the tenuous upper atmosphere are physically similar to cm-size streamer discharges at normal temperature and pressure up to corrections due to different electron attachment and detachment reactions that can explain long-delayed sprites (Luque and Gordillo-Va´zquez 2012). Sprites are pure streamer discharges (Liu and Pasko 2004b, a) and therefore are less complex than lightning strokes with their streamer, leader and return stroke stages, evolving on very different scales of space, time and energy. Due to the efforts of many authors in the past 20 years, the models for streamer discharges are now becoming more quantitative, so that we now approach the quantitative understanding of sprite discharges through detailed modelling and experimental efforts (Nijdam et al. 2014). Gamma-Ray Flashes and Other High-Energy Emissions from Thunderstorms In 1994, the BATSE14 satellite detected gamma radiation from Earth, and it was recognised that this radiation came from thunderstorms (Fishman et al. 1994; Fishman and Meegan 1995) . Later also beams of electrons (Dwyer et al. 2008) and even positrons (Briggs et al. 2011) were discovered by satellites. The Fermi Gamma-Ray Space Telescope detected a clear positron annihilation signal over Egypt from a thunderstorm over Zambia where the two events were connected in space and time through a geomagnetic field line (that electrons and positrons follow sufficiently high in the ionosphere where collisions with air molecules are negligible; Briggs et al. 2011). High-energy X-rays were also detected from lightning leaders approaching ground and from long sparks in the laboratory, see, e.g., Kochkin et al. (2012). We refer to the review by Dwyer and Uman (2014) . It is clear that electrons are accelerated into the run-away regime within the electric fields inside and above the thunderstorm, where they continuously gain more energy from the field than they can lose in collisions with neutral air molecules. These collisions with molecules result in X- or gamma-ray emission (Bremsstrahlung). The gamma rays are ionising radiation and 13 The electron density at these altitudes is an important parameter for discharge modelling. Only recently a method was developed to determine it partially and indirectly (Lay et al. 2010; Shao et al. 2013) . 14 generate electron positron pairs or liberate neutrons or protons in photonuclear reactions (Babich et al. 2014). There are two basic mechanisms discussed in the literature for the primary electron acceleration: either galactic cosmic rays with sufficient energy to penetrate deep into the atmosphere and to generate relativistic run-away electrons avalanches (RREAs) in the electric fields inside the thundercloud, or the acceleration of low-energy free electrons into the high-energy run-away regime at the tip of a lightning leader where electric fields are very high. The review by Dwyer and Uman (2014) favours the RREA mechanism, in agreement with the previous model development by the first author. The alternative is the runaway of thermal electrons at the leader tip suggested by Xu et al. (2012) . Such detailed models depend on the model parameters for the background cloud field and its geometry, on the altitude of the lightning leader, but also on the collision cross-sections at the required energies that are not reliably available. F u¨llekrug et al. (2013) reported on the observation of two consecutive positive lightning discharges where the first positive lightning discharge initiates sprite streamers which discharge the lightning electromagnetic field above the thundercloud. This was seen as a pulsed discharge event followed by a high-energy electron beam. A small number of stratospheric, charged aerosols were probably present as result of a Sahara dust storm and forest fires in Spain, providing a collimating electric field geometry that accelerated the electrons. This is the first simultaneous detection of radio signatures from electrons accelerated to thermal and relativistic energies above thunderclouds. 3.3 The Wilson Global Circuit The vertical structure of conductivity in the atmosphere, with the upper and lower conducting regions each able to sustain a local potential, allows a vertical potential difference to exist between the two regions. Investigations using balloon measurements from the late 1800s showed a variation in potential with atmospheric height (Nicoll 2012), with the upper conducting region being about 250kV positive with respect to the lower conducting region. The finite conductivity of the intermediate atmosphere between these charged regions allows a vertical current to flow. This current was observed directly by CTR Wilson (Wilson 1906) in fair weather conditions with no local charge separation. CTR Wilson concluded that the current flow was likely to be sustained by charge separation in distant disturbed weather regions. Evidence supporting this is that the diurnal variation in Universal Time (UT) near-surface electric field, measured under fair weather conditions, is independent of where it is measured globally and shows strong similarities with the diurnal variation in active global thunderstorm area (Whipple and Scrase 1936) . This diurnal variation in surface atmospheric electric field is known as the Carnegie curve, after the sailing vessel on which the original defining measurements were made (Harrison 2013). The conceptual model that described the electrical transport across the planet between disturbed weather and fair weather zones—the global atmospheric electric circuit (Wilson 1921, 1929) —has provided a fruitful description for investigation of terrestrial atmospheric electrification, which may offer useful insights for other atmospheres (Aplin et al. 2008). Although the original reasoning used to identify the global circuit was based on current flow considerations, the wide range of timescales of the contributing processes leads to a distinction being made conventionally between the AC and DC global circuit (Rycroft and Harrison 2012). The AC Global Circuit The upper and lower conducting regions of the terrestrial atmosphere form a simple waveguide, in which electromagnetic waves can propagate, as originally predicted by Schumann (Schumann 1952). Lightning provides a source of such electromagnetic radiation to excite waves in this cavity oscillator, and natural resonances with a fundamental mode at about 8 Hz as predicted were first observed at the Earth’s surface in the 1960s (Balser and Wagner 1960; Rycroft 1965) . These natural resonances in the Earth-ionosphere cavity (Q resonator) constitute the AC global electric circuit. Somewhat surprisingly, resonances at 8, 14, 20 Hz are also observed on satellites at altitudes of several hundred km, above the ionosphere (Simo˜ es et al. 2011; Dudkin et al. 2014) . Although the electric field measured is much smaller at a satellite platform compared with ground-based measurements (three orders of magnitude smaller for the first Schumann peak), the fact that it is detectable at all offers the possibility for fly-by measurements at other planetary bodies. The DC Global Circuit Figure 7 summarises the DC current flow in the Wilson global circuit. Charge separation in disturbed weather regions leads to current flow within the ionosphere, fair weather regions and the planetary surface. The vertical conduction current density, Jc, in fair weather regions is 2pA m 2, where the resistance of a unit area column of atmosphere, Rc, is about 100 to 300 PXm2 (Rycroft et al. 2000) . If horizontal layers of cloud or particles are present, the electrical conductivity is reduced because of the removal of the ions providing the conductivity by the particles. Hence, for a passive particle layer, this means that the layer also defines a region of reduced conductivity. If a current passes vertically through the passive particle layer (PPL), charging will result at the step change in conductivity at the upper and lower layer boundaries. The charging can be simulations as well as more realistic representations of the spacecraft geometry, so that shadowing effects can be studied more carefully. 5 Charge Processes in Extrasolar Atmospheric Environments Charge processes and their effects occur in many astrophysical environments. This section focuses on extrasolar objects where charge and discharge processes introduce feedback cycles similar to those discussed previously. This section summarises a field of research on cool extrasolar objects which starts to emerge as the result of recent progress in X-ray and radio observations of brown dwarfs. Brown dwarfs are objects with mass intermediate between stars and planets (Fig. 1; for a review see Helling and Casewell 2014) . Since they are not sufficiently massive for hydrogen burning in their core, they cool during their entire life time. Brown dwarf atmospheres therefore evolve from the state of a warm stellar atmosphere into an atmosphere as cool as the atmospheres of solar system planets. The oldest brown dwarfs are amongst the oldest objects in our Universe. Very-low mass stars and brown dwarfs are collectively known as ultracool dwarfs. Charge processes are important also in star and during planet formation in protoplanetary disks. Ionisation processes are suggested to help the first steps of planet formation as demonstrated by microgravity coagulation experiments at the International Space Station (Konopka et al. 2005; see also Sect. 3.4) and to allow the star to continue to accrete mass through the propolanetary disk. Planets and stars do have magnetic fields. The magnetic field strength and geometry differ for different stars (Donati and Landstreet 2009) which has implications, e.g., the size of a planetary magnetopause or the high-energy radiation impact into a planetary atmospheres (Vidotto et al. 2014; See et al. 2014) . Stellar magnetism changes with mass and rotation (as an indicator of the age of some stars) of the objects (Donati et al. 2008; Morin et al. 2008, 2010; Vidotto et al. 2014) introducing an additional complexity in the astrophysical context of atmospheric electrification. Brown dwarfs can have magnetic field strengths of 1000G ð¼ 0:1 TÞ. This section first summarises recent multi-wavelength observation of brown dwarfs as the best detectable ultracool and planet-like objects (Sect. 5.1). Section 5.2 addresses ionisation mechanisms in ultracool atmospheres, and Sect. 5.3 summarises recent ideas for ionising protoplanetary disks through which stars grow and planets form. 5.1 Multi-Wavelength Observations of Activity on Ultracool Dwarfs Below the mid-M spectral type stars (Teff \3200 K; for definitions see the Appendix), a strongly declining Ha emission indicates a weakening of chromospheric activity16 (Gizis et al. 2000; Kirkpatrick et al. 2000; Liebert et al. 2003; Reiners and Basri 2008; Williams et al. 2014) . X-ray observations support this finding: While X-ray detections are common for late-M spectral types (Teff [ 3400 K), the X-ray luminosity declines steeply for L-type brown dwarfs (Teff \2000 K; Williams et al. 2014). By contrast, brown dwarfs are very bright radio emitters, but no correlation between X-ray and radio emission exists as is 16 Chromospheric activity in the form of Ha, X-ray or Ca II K&K line emission results from the interaction of the stellar radiation field with a hot plasma above the atmosphere of a stellar object. The hot plasma that forms the chromosphere is the result of magnetic wave dissipation into thin gases and/or the deposition of excess radiation energy. established for stars and solar events (G u¨del-Benz relation, Benz and Gu¨ del 1994). Since the discovery of the first radio emitting brown dwarf by Berger et al. (2001), numerous surveys have detected radio emission at GHz frequencies in nearly 200 objects (Berger 2002, 2006; Berger et al. 2010; Hallinan et al. 2006, 2007, 2008; McLean et al. 2012), including emission in the coolest brown dwarfs with spectral types as late as T6.5 (Route and Wolszczan 2012) . In twelve of these objects, the radio emission is highly polarised, coherent and pulses on the rotation period of the dwarf. These properties suggest that the source of the radio emission is the electron cyclotron maser instability (CMI; Wu and Lee 1979) . The electron cyclotron maser mechanism has been shown to be responsible for the auroral kilometric radiation on Earth (see, Trakhtengerts and Rycroft 2008; Speirs et al. 2008; Vorgul et al. 2011) . Figure 13 shows a light curve of a M-dwarf (TVLM 513-46546) of spectral type M9.5 which shows clear and periodically repeating emission peaks at 8.44 GHz. In the solar system, planets have been extensively shown to be closely associated with auroral emission, caused when electrons moving along the magnetic field lines impact the atmosphere. Nichols et al. (2012) show that a model designed to explain Jupiter’s aurora (Cowley and Bunce 2001) is able to explain the observed radio fluxes in the ultracool dwarf ‘‘pulsars’’17 of order MW Hz 1. It could therefore be possible that the radio emissions of some ultracool dwarfs are powered by auroral currents. If this is true, there are profound implications for the importance of ionisation processes on ultracool dwarfs. The current system described by Nichols et al. (2012) requires, however, both a seed ionisation in the atmosphere and a plasma in the magnetosphere to operate. In turn, the impact of auroral electrons on the atmospheres is likely to be dramatic; while Jupiter’s aurora increase the atmospheric conductivity by a factor of 1000 (e.g., Strobel and Atreya 1983; Millward et al. 2002), the radio power of ultracool pulsars is 10,000 times that of Jupiter. Jupiter’s seed plasma is largely driven by the solar wind and the volcanically active Jupiter moon Io. Brown dwarfs will not have such external plasma sources, unless they are part of a binary system with mass transfer. It is likely that these aurorae are linked to optical variability seen in the ultracool dwarf pulsars. This association is suggested by the fact that five of the six ultracool dwarf pulsars that have been observed at optical wavelengths show periodic variability on the same period as the radio emission. For comparison, 5 % of randomly chosen ultracool dwarfs show periodic variability. While the exact mechanism producing this optical variability is not yet clear, multi-colour observations of one ultracool dwarf pulsar have ruled out starspots as the cause (Littlefair et al. 2008). Near-IR Signature of Chromospheric Activity in Brown Dwarfs Recent near-IR observations with the AKARI satellite can only be explained under the assumption that a chromosphere comparable to the solar chromosphere is present. Theoretical studies of brown dwarf atmospheres predict that such low-temperature atmospheres are dominated by molecules and dust and that they can be well modelled by simple radiative equilibrium assuming thermodynamic equilibrium. However, AKARI observations in the near-infrared wavelength range suggest that also chromospheric activity plays an important role for the atmospheric structure, in particular for early-type brown dwarfs (Sorahana et al. 2014). Deviations between theoretical model spectra and observed spectra around 3.0 and 4:5 lm, sensitive to the upper atmospheric structure of brown dwarfs, suggest an additional heating 17 Classically, a pulsar (a pulsating radio star) is a neutron star that is highly magnetised and rapidly rotating. The emitted radiation can only be observed when the beam is pointing towards Earth. The term ‘‘ultracool dwarf pulsar’’ borrows this idea of beamed, lighthouse-like radiation. source in the upper atmosphere. The comparison of the model spectrum with the observed spectrum for a L4.5-type brown dwarf with moderate Ha emission, 2MASS J2224-0158, is shown in Fig. 14 as an example. Sorahana et al. (2014) construct a simple model that includes heating due to chromospheric activity which results in a dramatic change of the chemical structure of the atmosphere. The resulting model spectra of early-type brown dwarfs with chromospheric heating considerably improve the match with the observed spectra. This result suggests that chromospheric activity is essential to understanding the near-infrared spectra of brown dwarfs and that MHD processes can heat the upper atmosphere. A similar conclusion was reached by Schmidt et al. (2015) who photometrically examined a sample of 11820 M7-L8 dwarfs. Rodriguez-Barrera et al. (2015) have used a grid of model atmosphere simulations to demonstrate that it is reasonable to expect the formation of an ionosphere and, therefore, also a chromosphere in ultracool atmospheres such as on brown dwarfs. Irradiated Brown Dwarf Atmospheres Only a handful of systems are known where a brown dwarf is heated by a hot companion. These brown dwarfs have close orbits of a few days, and they transit their host star, giving a measure of the brown dwarf’s radius, which is inflated by the energy input from its star. White dwarf—brown dwarf binaries provide a case where the brown dwarfs are not outshone by their companions, and therefore an opportunity to study irradiated brown dwarfs. In five of these systems, WD0137-349B (Maxted et al. 2006), NLTT5306 (Steele et al. 2013) , SDSS141126þ200911 (Beuermann et al. 2013) , WD0837þ185 (Casewell et al. 2012) and GD1400B (Farihi and Christopher 2004) , the brown dwarf is known to have survived a phase of common envelope 2MASS J2224-0158(L4.5) [Tcr / log g / Teff / Tconst]=[1800 K / 5.0 / 1700 K / 1445 K (Teff×0.85)] ] y J m [ ν F Wavelength [µm] Fig. 14 Comparison of the model spectrum (red and green smooth lines) with the observed spectrum for the L4.5 brown dwarf 2MASS J22240158 (thin black spiky line), which is well explained by the heating model atmosphere (red line) taking into account the heating in the upper atmospheres evolution.18 WD0137-349 is the best studied system. It is photometrically variable at all wavelengths from the V band though to 8:0 lm. These variations are at the orbital period of the system and peak at 4:5 lm. Converting the dayside and nightside flux to brightness temperature (peak blackbody temperature with that flux at that wavelength) shows a temperature difference between the two hemispheres of 500 K, and a possible temperature inversion in the atmosphere. Casewell et al. (2015) compare the observed photometry fluxes (i.e. radiative fluxes measured in a certain wavelength interval) to models of irradiated brown dwarfs and show that the data are best fit by models that incorporate fullenergy circulation around the brown dwarf, but do not contain a temperature inversion. However, at 2 lm (K band) and at 4:5 lm, the flux of the brown dwarf is still much brighter than the model predicts. Casewell et al. (2015) suggest that UV irradiation can cause photochemical reactions in the upper brown dwarf’s atmosphere that produce large hydrocarbon molecules causing the brown dwarf to be brighter at 2 and 4:5 lm as were demonstrated for cosmic ray impact by Rimmer et al. (2014) . 5.2 Ionisation Processes in Ultracool Atmospheres Brown dwarf and planetary atmospheres are spectroscopically characterised by a rich ensemble of molecules (e.g., SiO, TiO, VO, CO, H2O, FeH) which led to the conclusion that such atmospheres are too cool for thermal ionisation to significantly influence the local chemistry or energy content. But brown dwarfs and extrasolar planets exist in a larger diversity and, hence, are exposed to very different environments: The cosmic ray flux will be different in an interstellar environment than in the solar system (compare Sect. 3). The chemical composition of the gas in atmospheres outside the solar system causes the formation of mineral clouds where the cloud particles are composed of a mix of silicates, iron and metal oxides (Helling 2003, 2009; Helling and Rietmeijer 2009), very similar to volcanic ash. Extrasolar mineral clouds are much larger than terrestrial clouds due to the larger extension of the atmospheres of brown dwarfs and giant gas planets. Also such 18 The phase of common envelope evolution in the binary star evolution involves the brown dwarf being engulfed by, and immersed in, the expanding atmosphere of the white dwarf progenitor as it evolves away from the main sequence (Fig. 15). These systems are very close and tidally locked resulting in one side of the brown dwarf continually being heated by its companion. Wavelength (microns) clouds are susceptible to charge and discharge processes through cosmic rays (Rimmer and Helling 2013) and turbulence-enhanced dust-dust collisions (Helling et al. 2011) . The study of the breakdown condition in non-solar system atmospheres suggests that different intracloud discharge processes dominate at different heights inside mineral clouds: local coronal (point discharges) and small-scale sparks at the bottom region of the cloud where the gas density is high, and flow discharges and large-scale sparks near, and maybe above, the cloud top (Helling et al. 2013) . Bailey et al. (2014) apply scaling laws to demonstrate that discharges will propagate farther in brown dwarf atmospheres than in atmospheres of giant gas planets. Brown dwarfs can be irradiated by a binary companion (Maxted et al. 2006; Casewell et al. 2013, Sect. 5.1) resulting in similar global circulation patterns as demonstrated for irradiated giant gas planets (e.g., Knutson et al. 2007; Dobbs-Dixon and Agol 2013) . If such a wind of sufficient high speed hits a sufficiently pre-ionised gas, Alfve´n ionisation can produce bubbles of gas with a degree of ionisation of 1 (Stark et al. 2013) . The surrounding atmosphere remains in its low degree of ionisation leading to an atmosphere with a time-dependent and spatially intermittent degree of ionisation. A sufficient degree of ionisation is the precondition to understand the magnetic coupled atmospheric gas responsible for radio emission in brown dwarfs (Sect. 5.1), and from magnetically driven mass loss of extrasolar planetary atmospheres (Tanaka et al. 2014). Other mechanisms for planets to lose mass are related to their host star’s radiation field (Murray-Clay et al. 2009). An important input for understanding ionisation processes is the global atmospheric structure, and results from radiative-hydrodynamic simulations are therefore discussed in more detail below. Large-Scale Modelling of Globally Circulating Extrasolar Atmospheres A prominent sub-class of extrasolar planets are the short-period gaseous planets. They are (and due to observational constraints will remain) the best characterised of all extrasolar planets. The short periods ( 3 days), gaseous nature and largely circular orbits suggest that the rotation rate of these planets is tidally locked to their orbital period. The result is a stationary day–night heating pattern across the surface. The proximity to their host stars means that the hot daysides will be highly irradiated, reaching temperatures of several thousand degrees. The large longitudinal temperature gradient between the day and night hemispheres drives atmospheric dynamics that transports heat from day to night sides. The efficiency of this advective transport is a subject of extensive multi-dimensional radiative hydrodynamical modelling. Phase-curve observations, consisting of infrared measurements of the planetary flux throughout the entire orbit, do suggest that the night side of the planet is somewhat cooler than the day side (e.g., Knutson et al. 2007). However, the night-side temperatures are still on the order of a thousand degrees, much larger then one would expect simply from the internal cooling of the planet. This suggests that the atmosphere is actually fairly good at transporting energy across the entire planetary surface. The winds driven by the extreme temperatures on short-period planets are unlike any seen in any solar systems. The coupling between the (slow) planetary rotation and the temperature differential results in the development of a broad, supersonic, super-rotating equatorial jets (Tsai et al. 2014) . Gas velocities at pressure levels of 0.1 bars in the wellknown HD189733b can reach 5 km/s (e.g., Dobbs-Dixon and Agol 2013) . Thus, longitudinal transport and mixing, in contrast to vertical convection/turbulence as is important in Jupiter, plays a much larger role. The complexity of the atmospheric dynamics requires coupling together a dynamical model (solving the fluid equations) to a radiative model (involving the detailed opacities). Currently, models utilise molecular opacities, primarily due to species such as CO, H2O, and CH4. However, transit observations (Pont et al. 2013; Sing et al. 2014) taken for different wavelength bands suggest a cloud coverage throughout the atmosphere. While in hindsight based on observations of brown dwarfs or our gas giant planets this is not surprising, it makes the coupled problem significantly more complex. As the efficiency of energy transport by the gas flow depends on both the fluid velocity and the cooling timescale, the growth of cloud layers and the associated change in opacity will modify the radiative timescale and may have important consequences for energy redistribution. Conversely, the changes in the dynamics will alter the growth efficiency of clouds. Unfortunately, as with clouds on Earth, it is not at all clear if this will result in a net cooling or heating of various regions. The precision with which the local atmospheric properties, like gas temperature, cloud properties, gas composition, are modelled is crucial, and efforts are ongoing in the community by, for example, improving the treatment of the equation of state that covers a temperature range of 250K . . . 6000K and a gas pressure range of 10lbar to 10Mbar. The equation of state provides the abundances of opacity species (ions, atoms, molecules, cloud particles) that influence the local temperature through radiative transfer effects. Theoretical calculations now suggest that clouds should be very prevalent throughout these atmospheres (Lee et al. 2015). MHD simulations allow first insights into magnetic coupling effects despite containing much less information about atmospheric processes than radiative-hydrodynamic circulation models. Such MHD simulations have also been inspired by studies of protoplanetary disks (site of planet formation) and solar physics (e.g., Rogers and Showman 2014) . Tanaka et al. (2014), for example, have used an open magnetic flux-tube model in their ideal MHD simulation to demonstrate that under this conditions the planet could form a wind driven by Alfve´n waves. Murray-Clay et al. (2009) present an extensive study of UV and X-ray driven mass loss from irradiated extrasolar giant gas planets. 5.3 Discharges in Protoplanetary Disks Discharge in Magneto-Hydrodynamically Turbulent Gases Magneto-hydrodynamical turbulence is suggested as a mechanism to sustain the ionisation in a cool atmospheric gas. The energy dissipated from MHD turbulence is to ionise the gas. If this energy is large enough, a positive feedback loop develops where this local ionisation serves as driving mechanisms for magneto-hydrodynamical turbulence. Magneto-hydrodynamical turbulent motion in weakly ionised media continuously creates local electric fields even in the comoving frame of the media. If the electric field is larger than the critical electric field for electron avalanche, discharge occurs in such a weakly ionised media. This idea has been proposed for the turbulence driven by magneto-rotational instability in protoplanetary disks by Si and Sano (2005) who found that energetic electrons in a Druyvesteyn distribution may produce impact ionisation in particular conditions of the dusty gaseous disks. The atmospheres of gaseous planets and brown dwarfs are non-uniform in chemical composition for various reasons (formation of dust grains and their sedimentation, temperature stratification, occasional accretion of planetesimals, etc., see Sect. 5.2). The convection of those objects can be ‘‘double-diffusive convection’’ where the structure is destabilised by the diffusion of an elemental abundance gradient. Recent numerical simulations show that the double-diffusive convection evolves into a multi-layer structure where double-diffusive convection is confined into thin layers and usual convection occupies most of the volumes, which results in very small energy flux in the radial direction (e.g., Rosenblum et al. 2011a, b) . The effect of magnetic field can be important in the thin layers, since the magnetic tension force is inversely proportional to the length scale of the eddy. If the ionisation degree is kept high enough in the thin layers, the magnetic field possibly lowers the convective energy flux even further and hence slows down the gravitational contraction of those objects. Therefore, any processes that may increase the ionisation degree is important in the theory of long-term evolution of very cool brown dwarfs and gaseous planets. A simple energetics argument based on the order of magnitude calculations shows that the energy required for keeping the ionisation degree sufficiently high (magnetic Reynolds number Fig. 16 HCOþ (3-2) line profiles [Jy] at 267.56 GHz (1.1205 mm) without (left) and with (right) lightning. Lightning is here understood in analogy to Earth lightning as a large-scale discharge process. Different line styles (N, T, DP, R) indicate different breakdown models. The disk inclination of 7deg at a distance of 56 pc is similar to TW Hya. A minimum-mass solar nebula was applied, and a lightning region 25AU . . .50 AU was considered. The line flux considerably increases if lightning occurs. The ALMA sensitivity limit is 0.01 Jy for this line -150 -100 -50 0 05 001 150 -150 -100 -05 0 05 001 150 -150 -100 -50 0 05 001 150 -150 -100 -50 0 05 001 150 -150 -100 -50 0 05 001 150 -150 -100 -05 0 05 001 150 -150 -100 -50 0 05 001 150 -150 -100 -50 0 05 001 150 -150 -100 -50 0 05 001 150 -150 -100 -50 0 05 001 150 Fig. 17 Simulated integrated emissions maps [Jy beam 1 km s 1] for the HCOþ line at 267.56 GHz (1.1205 mm) without (top) and with (middle and bottom row, classic Townsend breakdown model) lightning. The lightning is considered to occur in a region of 25AU. . .50AU for the middle row, and of 50AU. . .100AU for the bottom row. These ALMA channel maps were simulated for 10 mJy per beam of 0’’.65x0’’.44. HCOþ appears in a larger fraction of a protoplanetary disk through the effect of lightning [1) is substantially small (\10-4) compared to the available energy as a turbulent convective motion. The critical electric field for impact ionisation in the astrophysical dusty plasma is calculated in detail as a function of gas density in Okuzumi and Inutsuka (2014) who also show that the resultant Ohm’s law is highly nonlinear and requires a new method to handle the magneto-hydrodynamics in particular regime. To determine the viability of the proposed process theoretically, magneto-hydrodynamics numerical simulations incorporating the microphysics of electron impact ionisation are required (e.g., Muranushi et al. 2012, 2013) . Observation of Lightning in Protoplanetary Disks by Ion Lines Lightning, a largescale discharge process in analogy to Earth-lightning, in protoplanetary disks has been studied as a candidate mechanism for chondrule formation, and it provides a unique window to probe the electromagnetic state of the protoplanetary disks. Evidence for strong (500–1000 G = 0.05-0.1T; Wasilewski and Dickinson 2000) , transient magnetic fields is found in meteorites. As a consequence, multiple lightning models have been proposed for protoplanetary disks (Gibbard et al. 1997; Weidenschilling 1997; Desch and Cuzzi 2000; Muranushi 2010) . Muranushi (2010) calculate the charge distribution of dust in a protoplanetary disk where a magneto-resonance instability produces an electromagnetic field. If the electric field potential is large enough for an ensemble of insulated but charged dust particles, a field breakdown will occur similar to the field breakdown in dust clouds of brown dwarfs and extrasolar planets (Helling et al. 2013) . Different breakdown models (Townsend, Druyversteyn-Penning, Runaway) can be tested which lead to different values for the breakdown field influencing the shape of the line profile (Fig. 16, N—no field). The electric field accelerates the free electrons that ionise the surrounding gas but also the positively charged ion species to the energy comparable to the electrons. Because the ionisation energy is a universal constant for each individual species, each ion will move with a characteristic, constant velocity in the lightning zone that is larger than the thermal gas velocity. This will be unique observational feature to detect and distinguish breakdown models in protoplanetary disks. In this model, it is assumed that the fractional abundances of HCOþ relative to H2 is 9 10 9. The value is taken from r = 100 au, z = 2H (z is the height about the midplane, H is the pressure scale height), from the X-ray and ultraviolet disk chemistry model of Walsh et al. (2012) . The underlying assumption is that HCOþ has been gradually produced and, although the HCOþ molecules may experience a sudden accelerated by the lightning electric fields, this does not contribute to the change of the number density of HCOþ in the present model. The velocity distribution of the ion species (e.g., HCOþ, DCOþ and N2Hþ) can now be derived and the line profiles simulated (Fig. 16 for HCOþ). The two-dimensional positionvelocity images with lightning assumed to occur in a certain disk region (middle, bottom) is shown as simulated ALMA19 channel maps in Fig. 17. The change in the line profiles depending on the presence of a large-scale discharge is demonstrated in Fig. 16. We found that lightning features of 10–100 mJy appear in the line profile. Using ALMA, full-disk lightning will produce 100r signals at 56 pc (TW Hya, Fig. 17) and 20r signals at 140 pc (Orion nebula; see Muranushi et al. 2015). 5.4 Future Studies Combining the expertise available on solar system and terrestrial atmosphere research and electrical phenomena therein will help in answering some of the following questions, but might also be inspired by these questions: How far can terrestrial and solar system lightning statistics guide our expectation for extrasolar environments like on extrasolar planets, brown dwarfs and in protoplanetary disks? How are brown dwarf (and extrasolar planet) atmospheres affected by the irradiation of their companion (or host star)? Going beyond the solar system, clouds are present in brown dwarf and exoplanetary atmospheres. What kind or which combination of atmospheric electricity phenomena could explain the required levels of ionised atmosphere to provide an explanation for the continuous radio emissions, the 656 nm and X-ray emissions? What would be the possible atmospheric optical and chemical signatures in the case that lightning activity exists in exoplanetary and brown dwarfs atmospheres? 19 Atacama Large Millimetre Array, 6 Conclusions Electrification processes and electrical phenomena are ubiquitous: dust charging and discharging are linked to the electrical breakdown of gas in planetary atmospheres inside and beyond the solar system where it is involved in global circuits and the occurrence of plasma processes. Charge processes play a major role in modifying the ambient chemical composition and the transport properties of neutral gas also in protoplanetary disks where planets form. Charged relativistic grains are suggested as potential primary particles for ultrahigh-energy Cosmic Rays (Hoang et al. 2015). The following set of challenges has emerged as common for the themes of this paper which have been guided by the Workshop ‘‘Electrification in dusty atmospheres inside and outside the solar system’’ held in September 201420: An increased population of ions, free electrons and radicals, lowers the chemical potential for gas-species reactions, leading to potentially observable spectroscopic fingerprints, and also increases the thermal and electrical conductivity of the gas to a certain threshold, enabling more energetic phenomena such as lightning to take place or accretion to proceed during star and planet formation. The presence of free charge may be transient (as in lightning), but the electrostatic influence can endure: Charging of dust and aerosols can influence the electrostatic character of the ambient atmosphere on longer length and time scales to produce small non-thermal populations of energetic particles. Finite enhanced electrical conductivity can allow magnetic relaxation, and access to stored magnetic energy as a general source of excitation which in unavailable to neutral gases. Non-thermal electrons may facilitate chemical reactions in ways that are classically (i.e. gas-thermodynamically) unlikely: for example, dissociative electron attachment can produce oxygen radicals at little energetic cost, leading to oxidative reactions proceeding at a rate inconsistent with ambient temperatures, or the formation of complex carbohydrates. Charged dust may evolve differently from neutral dust: Long-range organisation produced by electrostatic effects could produce coherent dust dynamics that would not be possible if only fluid mechanics dominates—surface charging, leading to elongated growth (nonzero eccentricity: polarisation of light is observable), or destruction of part of the grain population by Coulomb explosion. These themes lead to the need for: Further research on dust charging mechanisms in the context of volcanic lightning, the atmospheres of brown dwarfs and exoplanets, and protoplanetary disks. New instruments for future space missions to test new findings about electrical activity on solar system planets. The results might provide models that could mimic electrical activity on brown dwarfs and exoplanets. Further research about the role of dust in electrical discharges in the upper atmosphere of the Earth. New key insights into charging mechanisms of the Moon and asteroid fine dust grains by interacting with the solar wind and UV flux. This contributes to the fundamental knowledge of the Moon’s electric environment and will be very useful for further human missions to the Moon and unmanned missions to the asteroids. 3D simulations of extrasolar atmospheres including chemical and electrical feedbacks of clouds in a magnetised gas. Acknowledgments Christian Helling highlights the financial support of the European Community under the FP7 by an ERC starting Grant 25473. Declan Diver gratefully acknowledges support from EPSRC via Grant Numbers EP/K006142/1 and EP/K006088/1. Francisco Gordillo-Vasquez thanks the Spanish Ministry of Economy and Competitiveness (MINECO) under projects FIS2014-61774-EXP and ESP2013-48032-C55-R and the EU through the FEDER programme. We thank all the participants at the Workshop Electrification in dusty atmospheres inside and outside the solar system held in September 2014 in the Scottish Highlands for their input and inspiration. We thank the Royal Astronomical Society, the European Research Council and the Institute of Physics Electrostatics Group for financial support. We thanks Sarah Casewell and Isabel Rodrigues-Barrera for their help in preparing Table 1. Gabriella Hodosa´n is thanked for helping with the literature collection. Keri Nicoll is thanked for her inspiring feedback on Sect. 3.4. Most of the literature search has been performed using ADS. Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, 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. Appendix: Glossary AC: alternating current. Asteroid: small rocky bodies of the inner solar system, ranging in size from 10 to 900 m in diameter. Aurora: large diffuse light-emitting structures in the lower ionosphere ([90 km) generated when energetic charged particles precipitating from the magnetosphere collide with ground state neutral species and excite them. The excited species (oxygen atoms and nitrogen molecules) emit light when returning to their ground state. Carbonaceous Compound: material rich in carbon; in an astrophysical context, such compounds usually are associated with primitive solar system remnants. Conduit, Volcano Conduit: the pipe that carries magma from the magma chamber, up through the crust and through the volcano itself until it reaches the surface. Chondrule: molten or partially molten droplets that appear as spherical, solid inclusions of different chemical composition from the the matrix of their parent asteroid. They represent one of the oldest solid materials within the solar system. Cosmic Rays: ionised nuclei and electrons that are distinguished by their high energies. The ionised nuclei have energies ranging from 106 eV to greater than 1020 eV and comprise 99 % of cosmic rays. They originate from either the Sun or outside the solar system, probably from Super Novae or Gamma Ray Bursts. Cyclotron Maser Instability: the mechanism whereby a population of relativistic electrons move in ambient magnetic fields, producing coherent radiation that reflects depends upon the magnetic field strength. DC: direct current. Debye Length: the scale-length associated with the violation of charge neutrality in a plasma, due to thermal fluctuations causing charge separation. Double-Diffusive Convection: a form of convection (i.e. hydrodynamic bulk motion) that is driven by two distinct gradients in fluid composition arising from two different species abundances. Near-IR: electromagnetic radiation in the near infrared, with wavelengths between 800 nm and 5 microns. Ohm’s Law: relates the electrical current flowing between two points to the potential difference between those same points. Plasma Void: a finite region in a dusty plasma which is dust free. Protoplanetary Disk: a region of dust and rocks orbiting a young star from which planets could be formed. Regolith: layer of unconsolidated dust and fragmented rock that covers a terrestrial planet. Shock Tube: a device designed to create shocks (i.e. sharp density and pressure discontinuities) in gases, usually in order to produce ionisation fronts. Sounding: a method to measure local temperature, humidity, wind, etc., in the Earth’s atmosphere by means of radio sondes, laser beams (optical) or sound waves. Sprite: a large (50 km high and 10–20 km wide) electrical discharge that occurs above thunderclouds at altitudes around 50–85 km with a diffuse region (above 70–75 km) and a filamentary (streamer-like) region (below 70–75 km) [values are given for Earth]. Stratosphere: major layer of the Earth’s atmosphere, lying above the troposphere (0–11 km) and below the mesosphere (50–90 km). Thunderstorm: storm characterised by the presence of lightning. Triboelectrification: the process whereby two surfaces can acquire or lose charges by mutual collision. Troposphere: lowest layer of the Earth’s atmosphere, lying between 0 and 11km, in which most of the weather phenomena occur. Turbulence: chaotic flow in which the pressure and gas velocity change rapidly in space and time. Volcanic Ash: fragments of rock created during a volcanic eruption, usually 2 mm or less in diameter. Volcanic Conduit: passage or tube created by the flow of magma in a volcano. Volcanic Plume: the gas and ash cloud ejected into the atmosphere by a volcanic eruption. Harrison RG, Carslaw KS (2003) Ion-aerosol-cloud processes in the lower atmosphere. Rev Geophys 41:1012. doi:10.1029/2002RG000114 Harrison RG, Mather TA (2006) The electrification of volcanic plumes and volcanic lightning. In: AGU fall meeting abstracts, p A286 Harrison RG, Nicoll KA (2014) Note: active optical detection of cloud from a balloon platform. Rev Sci Instrum 85(6):066104. doi:10.1063/1.4882318 Harrison RG, Aplin KL, Leblanc F, Yair Y (2008) Planetary atmospheric electricity. Space Sci Rev 137:5–10. doi:10.1007/s11214-008-9419-z Harrison RG, Nicoll KA, Ulanowski Z, Mather TA (2010) Self-charging of the Eyjafjallajo¨kull volcanic ash plume. Environ Res Lett 5(2):024004. doi:10.1088/1748-9326/5/2/024004 Harrison RG, Nicoll KA, Aplin KL (2014) Vertical profile measurements of lower troposphere ionisation. J Atmos Sol Terr Phys 119:203–210. doi:10.1016/j.jastp.2014.08.006 Helling Ch (2003) Circuit of dust in substellar objects (with 10 figures). In: Schielicke RE (ed) Reviews in modern astronomy, vol 16, p 115 Helling Ch (2009) Cloud formation in substellar atmospheres. In: Stempels E (ed) 15th Cambridge workshop on cool stars, stellar systems, and the sun, American Institute of Physics Conference Series, vol 1094, pp 162–171. doi:10.1063/1.3099087 Helling Ch, Casewell S (2014) Atmospheres of brown dwarfs. Astron Astrophys Rev 22:80. doi:10.1007/ s00159-014-0080-0 Helling C, Rietmeijer FJM (2009) Glittery clouds in exoplanetary atmospheres? Int J Astrobiol 8:3–8. doi:10.1017/S1473550408004382 Helling Ch, Jardine M, Mokler F (2011) Ionization in atmospheres of brown dwarfs and extrasolar planets. II. Dust-induced collisional ionization. Astrophys J 737:38. doi:10.1088/0004-637X/737/1/38 Helling Ch, Jardine M, Stark C, Diver D (2013) Ionization in atmospheres of brown dwarfs and extrasolar planets. III. Breakdown conditions for mineral clouds. Astrophys J 767:136. doi:10.1088/0004-637X/ 767/2/136 Hess VF (1912) U¨ ber Beobachtungen der durchdringenden Strahlung bei sieben Freiballonfahrten. Phys Z 13:1084 Hess SLG, Sarrailh P, Mate´o-Ve´lez JC, Forest J, Jeanty-Ruard B, Cipriani F (2015) Simulation of the electrostatic charging of Philae on 67P/Churyumov-Gerasimenko and of its interaction with the dusts. In: Martins F, Boissier S, Buat V, Cambre´sy L, Petit P (eds) SF2A-2015: Proceedings of the Annual meeting of the French Society of Astronomy and Astrophysics, pp 57–61 Hoang T, Lazarian A, Schlickeiser R (2015) On origin and destruction of relativistic dust and its implication for ultrahigh energy cosmic rays. Astrophys J 806:255. doi:10.1088/0004-637X/806/2/255 Hodosa´n G, Helling Ch, Asensio-Torres R, Vorgul I, Rimmer PB (2016) Lightning climatology of exoplanets and brown dwarfs guided by solar system data. Mon Not R Astron Soc (submitted) Hora´nyi M, Hartquist TW, Havnes O, Mendis DA, Morfill GE (2004) Dusty plasma effects in Saturn’s magnetosphere. Rev Geophys 42:RG4002. doi:10.1029/2004RG000151 Houghton IMP, Aplin KL, Nicoll KA (2013) Triboelectric charging of volcanic ash from the 2011 Gr´ımsvo¨tn eruption. Phys Rev Lett 111(11):118501. doi:10.1103/PhysRevLett.111.118501 Hutchinson IH, Patacchini L (2007) Computation of the effect of neutral collisions on ion current to a floating sphere in a stationary plasma. Phys Plasmas. doi:10.1063/1.2431584 Inutsuka S, Sano T (2005) Self-sustained ionization and vanishing dead zones in protoplanetary disks. Astrophys J Lett 628:L155–L158. doi:10.1086/432796 Irwin PGJ, Teanby NA, de Kok R, Fletcher LN, Howett CJA, Tsang CCC, Wilson CF, Calcutt SB, Nixon CA, Parrish PD (2008) The NEMESIS planetary atmosphere radiative transfer and retrieval tool. J Quant Spectrosc Radiat Transf 109:1136–1150. doi:10.1016/j.jqsrt.2007.11.006 James MR, Lane SJ, Gilbert JS (1998) Volcanic plume monitoring using atmospheric electric potential gradients. Geol Soc Lond J 155:587–590. doi:10.1144/gsjgs.155.4.0587 James MR, Lane SJ, Gilbert JS (2000) Volcanic plume electrification: experimental investigation of a fracture-charging mechanism. J Geophys Res 105(16):641. doi:10.1029/2000JB900068 Jayaratne ER, Saunders CPR, Hallett J (1983) Laboratory studies of the charging of soft-hail during ice crystal interactions. Q J R Meteorol Soc 109:609–630. doi:10.1002/qj.49710946111 Johnson AP, Cleaves HJ, Dworkin JP, Glavin DP, Lazcano A, Bada JL (2008) The Miller volcanic spark discharge experiment. Science 322:404 Khrapak SA, Morfill GE (2008) An interpolation formula for the ion flux to a small particle in collisional plasmas. Phys Plasmas 15(11):114503. doi:10.1063/1.3035913 Khrapak SA, Tolias P, Ratynskaia S, Chaudhuri M, Zobnin A, Usachev A, Rau C, Thoma MH, Petrov OF, Fortov VE, Morfill GE (2012) Grain charging in an intermediately collisional plasma. Europhys Lett. doi:10.1209/0295-5075/97/35001 McConville SL, Speirs DC, Ronald K, Phelps ADR, Cross AW, Bingham R, Robertson CW, Whyte CG, He W, Gillespie KM, Vorgul I, Cairns RA, Kellett BJ (2008) Demonstration of auroral radio emission mechanisms by laboratory experiment. Plasma Phys Controll Fusion 50(7):074010. doi:10.1088/07413335/50/7/074010 McLean M, Berger E, Reiners A (2012) The radio activity-rotation relation of ultracool dwarfs. Astrophys J 746:23. doi:10.1088/0004-637X/746/1/23 McNutt SR, Williams ER (2010) Volcanic lightning: global observations and constraints on source mechanisms. Bull Volcanol 72:1153–1167. doi:10.1007/s00445-010-0393-4 Michel P, Barucci MA, Cheng AF, Bo¨hnhardt H, Brucato JR, Dotto E, Ehrenfreund P, Franchi IA, Green SF, Lara LM, Marty B, Koschny D, Agnolon D (2014) MarcoPolo-R: near-Earth asteroid sample return mission selected for the assessment study phase of the ESA program cosmic vision. Acta Astronaut 93:530–538. doi:10.1016/j.actaastro.2012.05.030 Millward G, Miller S, Stallard T, Aylward AD, Achilleos N (2002) On the dynamics of the Jovian ionosphere and thermosphere: III. The modelling of auroral conductivity. Icarus 160:95–107. doi:10.1006/ icar.2002.6951 Miura T, Koyaguchi T, Tanaka Y (2002) Measurements of electric charge distribution in volcanic plumes at Sakurajima Volcano, Japan. Bull Volcanol 64:75–93. doi:10.1007/s00445-001-0182-1 Moore LE, Mendillo M, M u¨ller-Wodarg ICF, Murr DL (2004) Modeling of global variations and ring shadowing in Saturn’s ionosphere. Icarus 172:503–520. doi:10.1016/j.icarus.2004.07.007 Morin J, Donati JF, Petit P, Delfosse X, Forveille T, Albert L, Aurie`re M, Cabanac R, Dintrans B, Fares R, Gastine T, Jardine MM, Lignie`res F, Paletou F, Ramirez Velez JC, The´ado S (2008) Large-scale magnetic topologies of mid M dwarfs. Mon Not R Astron Soc 390:567–581. doi:10.1111/j.1365-2966. 2008.13809.x Morin J, Donati JF, Petit P, Delfosse X, Forveille T, Jardine MM (2010) Large-scale magnetic topologies of late M dwarfs. Mon Not R Astron Soc 407:2269–2286. doi:10.1111/j.1365-2966.2010.17101.x Moses JI, Be´zard B, Lellouch E, Gladstone GR, Feuchtgruber H, Allen M (2000) Photochemistry of Saturn’s atmosphere. I. Hydrocarbon chemistry and comparisons with ISO observations. Icarus 143:244–298. doi:10.1006/icar.1999.6270 Muranushi T (2010) Dust–dust collisional charging and lightning in protoplanetary discs. Mon Not R Astron Soc 401(4):2641–2664 Muranushi T, Okuzumi S, Inutsuka S (2012) Interdependence of electric discharge and magnetorotational instability in protoplanetary disks. Astrophys J 760:56. doi:10.1088/0004-637X/760/1/56 Muranushi T, Okuzumi S, Inutsuka S (2013) Erratum: ‘‘Interdependence of electric discharge and the magnetorotational instability in protoplanetary disks’’ /abs/2012ApJ...760...56M (2012, ApJ, 760, 56). Astrophys J 771:138. doi:10.1088/0004-637X/771/2/138 Muranushi T, Akiyama E, Inutsuka SI, Nomura H, Okuzumi S (2015) Development of a method for the observation of lightning in protoplanetary disks using ion lines. Astrophys J 815:84. doi:10.1088/0004637X/815/2/84 Murray-Clay RA, Chiang EI, Murray N (2009) Atmospheric escape from hot Jupiters. Astrophys J 693:23–42. doi:10.1088/0004-637X/693/1/23 Nichols JD, Burleigh MR, Casewell SL, Cowley SWH, Wynn GA, Clarke JT, West AA (2012) Origin of electron cyclotron maser induced radio emissions at ultracool dwarfs: magnetosphere–ionosphere coupling currents. Astrophys J 760:59. doi:10.1088/0004-637X/760/1/59 Nicoll KA (2012) Measurements of atmospheric electricity aloft. Surv Geophys 33:991–1057. doi:10.1007/ s10712-012-9188-9 Nicoll KA (2013) Note: A self-calibrating electrometer for atmospheric charge measurements from a balloon platform. Rev Sci Instrum 84(9):096107. doi:10.1063/1.4821500 Nicoll KA (2014) Space weather influences on atmospheric electricity. Weather 69:238–241. doi:10.1002/ wea.2323 Nicoll KA, Harrison RG (2010) Experimental determination of layer cloud edge charging from cosmic ray ionisation. Geophys Res Lett 37:L13802. doi:10.1029/2010GL043605 Nijdam S, Takahashi E, Markosyan AH, Ebert U (2014) Investigation of positive streamers by double-pulse experiments, effects of repetition rate and gas mixture. Plasma Sources Sci Technol 23(2):025008. doi:10.1088/0963-0252/23/2/025008 Okuzumi S, Inutsuka S (2014) The nonlinear Ohm’s Law: plasma heating by strong electric fields and its effects on the ionization balance in protoplanetary disks. Astrophys J 800(1):47 Oppenheimer C (2003) Climatic, environmental and human consequences of the largest known historic eruption: Tambora volcano (Indonesia) 1815. Prog Phys Geogr 27:230–259. doi:10.1191/ 0309133303pp379r Rycroft MJ, Odzimek A, Arnold NF, Fu¨llekrug M, Kułak A, Neubert T (2007) New model simulations of the global atmospheric electric circuit driven by thunderstorms and electrified shower clouds: the roles of lightning and sprites. J Atmos Sol Terr Phys 69:2485–2509. doi:10.1016/j.jastp.2007.09.004 Rycroft MJ, Harrison RG, Nicoll KA, Mareev EA (2008) An overview of Earth’s global electric circuit and atmospheric conductivity. Space Sci Rev 137:83–105. doi:10.1007/s11214-008-9368-6 Saunders C (2008) Charge separation mechanisms in clouds. Springer, New York, pp 335–353. doi:10.1007/ 978-0-387-87664-1_22 Schmidt SJ, Hawley SL, West AA, Bochanski JJ, Davenport JRA, Ge J, Schneider DP (2015) BOSS ultracool dwarfs. I. Colors and magnetic activity of M and L dwarfs. Astron J 149:158. doi:10.1088/ 0004-6256/149/5/158 Schumann WO (1952) U¨ ber die strahlungslosen Eigenschwingungen einer leitenden Kugel die von einer Luftschicht und einer Ionospha¨renhu¨lle umgeben ist. Z Naturforsch Teil A 7:149 See V, Jardine M, Vidotto AA, Petit P, Marsden SC, Jeffers SV, do Nascimento JD (2014) The effects of stellar winds on the magnetospheres and potential habitability of exoplanets. Astron Astrophys 570:A99. doi:10.1051/0004-6361/201424323 Sen S, Fukuyama A, Honary F (2010) Rayleigh Taylor instability in a dusty plasma. J Atmos Sol Terr Phys 72:938–942. doi:10.1016/j.jastp.2010.05.001 Shalygin EV, Markiewicz WJ, Basilevsky AT, Titov DV, Ignatiev NI, Head JW (2015) Active volcanism on Venus in the Ganiki Chasma rift zone. Geophys Res Lett 42:4762–4769. doi:10.1002/2015GL064088 Shao XM, Lay EH, Jacobson AR (2013) Reduction of electron density in the night-time lower ionosphere in response to a thunderstorm. Nat Geosci 6(1):29–33. doi:10.1038/NGEO1668 Shukla PK, Mamun AA (2002) Introduction to dusty plasma physics. Institute of Physics Publishing, Philadelphia Shukla PK, Silin VP (1992) Dust ion-acoustic wave. Phys Scr 45:508. doi:10.1088/0031-8949/45/5/015 Silva HG, Conceica˜o R, Melgo˜a M, Nicoll K, Tlemcani M, Reis AH, Harrison RG (2014) Atmospheric electric field measurements in urban environment and the pollutant aerosol weekly dependence. Environ Res Lett 9(114):025 Simo˜es F, Pfaff R, Freudenreich H (2011) Satellite observations of Schumann resonances in the Earth’s ionosphere. Geophys Res Lett 38:L22101. doi:10.1029/2011GL049668 Sing DK, Wakeford HR, Showman AP, Nikolov N, Fortney JJ, Burrows AS, Ballester GE, Deming De (2014) HST hot-Jupiter transmission spectral survey: detection of potassium in WASP-31b along with a cloud deck and Rayleigh scattering. ArXiv e-prints Sorahana S, Suzuki TK, Yamamura I (2014) A signature of chromospheric activity in brown dwarfs revealed by 2.5–5.0 lm AKARI spectra. Mon Not R Astron Soc 440:3675–3684. doi:10.1093/mnras/ stu479 Speirs DC, McConville SL, Gillespie KM, Ronald K, Phelps ADR, Cross AW, Bingham R, Robertson CW, Whyte CG, Vorgul I, Cairns RA, Kellett BJ (2008) Numerical simulation of auroral cyclotron maser processes. Plasma Phys Control Fusion 50(7):074011. doi:10.1088/0741-3335/50/7/074011 Stark CR, Helling Ch, Diver DA, Rimmer PB (2013) Ionization in atmospheres of brown dwarfs and extrasolar planets. V. Alfve´n Ionization. Astrophys J 776:11. doi:10.1088/0004-637X/776/1/11 Steele PR, Saglia RP, Burleigh MR, Marsh TR, Ga¨nsicke BT, Lawrie K, Cappetta M, Girven J, Napiwotzki R (2013) NLTT 5306: the shortest period detached white dwarf?brown dwarf binary. Mon Not R Astron Soc 429:3492–3500. doi:10.1093/mnras/sts620 Stenbaek-Nielsen HC, McHarg MG (2008) High time-resolution sprite imaging: observations and implications. J Phys D Appl Phys. doi:10.1088/0022-3727/41/23/234009 Stozhkov YV, Okhlopkov V, Makhmutov V, Logachev A (2013) Solar activity, cosmic rays, and global climate changes. In: Proceedings of 33rd international cosmic ray conference (ICRC) Strobel DF, Atreya SK (1983) Ionosphere. In: Dessler A (ed) Physics of the Jovian Magnetosphere, vol 1. Cambridge University Press, Cambridge, pp 51–67 Taddeucci J, Scarlato P, Montanaro C, Cimarelli C, Del Bello E, Freda C, Andronico D, Gudmundsson MT, Dingwell DB (2011) Aggregation-dominated ash settling from the Eyjafjallajokull volcanic cloud illuminated by field and laboratory high-speed imaging. Geology 39:891–894. doi:10.1130/G32016.1 Tanaka YA, Suzuki TK, Inutsuka S (2014) Atmospheric escape by magnetically driven wind from gaseous planets. Astrophys J 792:18. doi:10.1088/0004-637X/792/1/18 Thomas RJ, Krehbiel PR, Rison W, Edens HE, Aulich GD, Winn WP, McNutt SR, Tytgat G, Clark E (2007) Electrical activity during the 2006 Mount St. Augustine volcanic eruptions. Science 315:1097. doi:10. 1126/science.1136091 Todd JF, Barber SJ, Wright IP, Morgan GH, Morse AD, Sheridan S, Leese MR, Maynard J, Evans ST, Pillinger CT, Drummond DL, Heys SC, Huq SE, Kent BJ, Sawyer EC, Whalley MS, Waltham NR Airey MW , Mather TA , Pyle DM , Glaze LS , Ghail RC , Wilson CF ( 2015 ) Explosive volcanic activity on Venus: the roles of volatile contribution, degassing, and external environment . Planet Space Sci 113 : 33 - 48 . doi: 10 .1016/j.pss. 2015 . 01 .009 Aizawa K , Yokoo A , Kanda W , Ogawa Y , Iguchi M ( 2010 ) Magnetotelluric pulses generated by volcanic lightning at Sakurajima volcano , Japan. Geophys Res Lett 37 : L17301 . doi: 10 .1029/2010GL044208 Allen JE , Phelps ADR ( 1977 ) Waves and microinstabilities in plasmas-linear effects . Rep Prog Phys 40 : 1305 - 1368 . doi: 10 .1088/ 0034 -4885/40/11/002 Anderson R , Bjornsson S , Blanchard DC , Gathman S , Hughes J , Jonasson S , Moore CB , Survilas HJ , Vonnegut B ( 1965 ) Electricity in volcanic clouds . Science 148 : 1179 - 1189 . doi: 10 .1126/science.148. 3674.1179 Anuar AK , Honary F , Hapgood M , Roussel JF ( 2013 ) Three-dimensional simulation of dust charging and dusty plasma using SPIS . J Geophys Res Space Phys 118 : 6723 - 6735 . doi: 10 .1002/jgra.50599 Aplin KL ( 2006 ) Atmospheric electrification in the solar system . Surv Geophys 27 : 63 - 108 . doi: 10 .1007/ s10712-005-0642-9 Aplin KL ( 2012 ) Smoke emissions from industrial western Scotland in 1859 inferred from Lord Kelvin's atmospheric electricity measurements . Atmos Environ 50 : 373 - 376 . doi: 10 .1016/j.atmosenv. 2011 . 12 . 053 Christian HJ , Blakeslee RJ , Boccippio DJ , Boeck WL , Buechler DE , Driscoll KT , Goodman SJ , Hall JM , Koshak WJ , Mach DM , Stewart MF ( 2003 ) Global frequency and distribution of lightning as observed from space by the optical transient detector . J Geophys Res Atmos 108 : 4005 . doi: 10 .1029/ 2002JD002347 Cimarelli C , Alatorre-Ibargu¨engoitia MA , Aizawa K , D´ıaz Marina AI , Yokoo A , Kueppers U , Mueller S , Scheu B , Dingwell DB ( 2015 ) Volcanic lightning: in nature and in the lab , vol 17 . In: EGU General Assembly conference abstracts, p 7194 Cooray V ( 2003 ) The lightning flash . IEE , London Cooray V ( 2015 ) An introduction to lightning . Springer, Dordrecht Cowley SWH , Bunce EJ ( 2001 ) Origin of the main auroral oval in Jupiter's coupled magnetosphereionosphere system . Planet Space Sci 49 : 1067 - 1088 . doi: 10 .1016/S0032- 0633 ( 00 ) 00167 - 7 Cushing MC , Rayner JT , Vacca WD ( 2005 ) An infrared spectroscopic sequence of M, L, and T dwarfs . Astrophys J 623 : 1115 - 1140 . doi: 10 .1086/428040 D'Angelo N ( 1993 ) The Rayleigh-Taylor instability in dusty plasmas . Planet Space Sci 41 : 469 - 474 . doi: 10 . 1016/ 0032 - 0633 ( 93 ) 90007 - O Demory BO , Gillon M , Deming D , Valencia D , Seager S , Benneke B , Lovis C , Cubillos P , Harrington J , Stevenson KB , Mayor M , Pepe F , Queloz D , Se´gransan D , Udry S ( 2011 ) Detection of a transit of the super-Earth 55 Cancri e with warm Spitzer . Astron Astrophys 533 : A114 . doi: 10 .1051/ 0004 -6361/ 201117178 Desch S , Cuzzi J ( 2000 ) The generation of lightning in the solar nebula . Icarus 143 ( 1 ): 87 - 105 Dobbs-Dixon I , Agol E ( 2013 ) Three-dimensional radiative-hydrodynamical simulations of the highly irradiated short-period exoplanet HD 189733b . Mon Not R Astron Soc 435 : 3159 - 3168 . doi: 10 .1093/ mnras/stt1509 Donati JF , Landstreet JD ( 2009 ) Magnetic fields of nondegenerate stars . Annu Rev Astron Astrophys 47 : 333 - 370 . doi: 10 .1146/annurev-astro- 082708 -101833 Donati JF , Morin J , Petit P , Delfosse X , Forveille T , Aurie`re M, Cabanac R , Dintrans B , Fares R , Gastine T , Jardine MM , Lignie`res F , Paletou F , Ramirez Velez JC , The´ado S ( 2008 ) Large-scale magnetic topologies of early M dwarfs . Mon Not R Astron Soc 390 : 545 - 560 . doi: 10 .1111/j.1365- 2966 . 2008 . 13799 .x Dubinova A , Rutjes C , Ebert U , Buitink S , Scholten O , Trinh GTN ( 2015 ) Prediction of lightning inception by large ice particles and extensive air showers . Phys Rev Lett 115 ( 1 ): 015002 . doi: 10 .1103/ PhysRevLett.115.015002 Dubrovin D , Luque A , Gordillo-Vazquez FJ , Yair Y , Parra-Rojas FC , Ebert U , Price C ( 2014 ) Impact of lightning on the lower ionosphere of Saturn and possible generation of halos and sprites . Icarus 241 : 313 - 328 . doi: 10 .1016/j.icarus. 2014 . 06 .025 Dudkin D , Pilipenko V , Korepanov V , Klimov S , Holzworth R ( 2014 ) Electric field signatures of the IAR and Schumann resonance in the upper ionosphere detected by Chibis-M microsatellite . J Atmos Solar Terr Phys 117 : 81 - 87 . doi: 10 .1016/j.jastp. 2014 . 05 .013 Dwyer JR , Uman MA ( 2014 ) The physics of lightning . Phys Rep 534 ( 4 ): 147 - 241 . doi: 10 .1016/j.physrep. 2013 . 09 .004. Dwyer JR , Grefenstette BW , Smith DM ( 2008 ) High-energy electron beams launched into space by thunderstorms . Geophys Res Lett 35 : L02815 . doi: 10 .1029/2007GL032430 Ebert U , Sentman DD ( 2008 ) Editorial Review: Streamers, sprites, leaders, lightning: from micro- to macroscales . J Phys D Appl Phys 41 ( 23 ): 230301 . doi: 10 .1088/ 0022 -3727/41/23/230301 Ebert U , Nijdam S , Li C , Luque A , Briels T , van Veldhuizen E ( 2010 ) Review of recent results on streamer discharges and discussion of their relevance for sprites and lightning . J Geophys Res Space Phys. doi:10.1029/2009JA014867 Farihi J , Christopher M ( 2004 ) A possible brown dwarf companion to the white dwarf GD 1400 . Astron J 128: 1868 - 1871 . doi: 10 .1086/423919 Farrell WM , Desch MD ( 2001 ) Is there a Martian atmospheric electric circuit ? J Geophys Res Planets 106 : 7591 - 7596 . doi: 10 .1029/2000JE001271 Fillingim M ( 1986 ) Global electric circuit of Mars . In: Florida institute of technology (unpublished report) Fishman GJ , Meegan CA ( 1995 ) Gamma-ray bursts . Annu Rev Astron Astrophys 33 : 415 - 458 . doi: 10 .1146/ annurev.aa. 33 .090195.002215 Fishman GJ , Bhat PN , Mallozzi R , Horack JM , Koshut T , Kouveliotou C , Pendleton GN , Meegan CA , Wilson RB , Paciesas WS , Goodman SJ , Christian HJ ( 1994 ) Discovery of intense gamma-ray flashes of atmospheric origin . Science 264 : 1313 - 1316 . doi: 10 .1126/science.264.5163.1313 Fletcher LN , Hesman BE , Achterberg RK , Irwin PGJ , Bjoraker G , Gorius N , Hurley J , Sinclair J , Orton GS , Legarreta J , Garc´ ıa-Melendo E , Sa ´nchez- Lavega A , Read PL , Simon-Miller AA , Flasar FM ( 2012 ) The Owens MJ , Scott CJ , Lockwood M , Barnard L , Harrison RG , Nicoll K , Watt C , Bennett AJ ( 2014 ) Modulation of UK lightning by heliospheric magnetic field polarity . Environ Res Lett 9 , 115009 . doi: 10 .1088/ 1748 -9326/9/11/115009 Parra-Rojas FC , Luque A , Gordillo-Va ´zquez FJ ( 2013 ) Chemical and electrical impact of lightning on the Earth's mesosphere: the case of sprite halos . J Geophys Res Space Phys 118 : 5190 - 5214 . doi: 10 .1002/ jgra.50449 Parra-Rojas FC , Luque A , Gordillo-Va ´zquez FJ ( 2015 ) Chemical and thermal impacts of sprite streamers in the Earth's mesosphere . J Geophys Res Space Phys 120 : 8899 - 8933 . doi: 10 .1002/2014JA020933 Pasko VP ( 2007 ) TOPICAL REVIEW: Red sprite discharges in the atmosphere at high altitude: the molecular physics and the similarity with laboratory discharges . Plasma Sources Sci Technol 16 : 13 . doi: 10 .1088/ 0963 -0252/16/1/S02 Pfotzer G ( 1972 ) History of the use of balloons in scientific experiments . Space Sci Rev 13 : 199 - 242 . doi: 10 . 1007/BF00175313 Phelps ADR , Allen JE ( 1976 ) A floating electrostatic sheath in a thermally produced plasma . R Soc Lond Proc Ser A 348 : 221 - 233 . doi: 10 .1098/rspa. 1976 .0034 Pont F , Sing DK , Gibson NP , Aigrain S , Henry G , Husnoo N ( 2013 ) The prevalence of dust on the exoplanet HD 189733b from Hubble and Spitzer observations . Mon Not R Astron Soc 432 : 2917 - 2944 . doi: 10 . 1093/mnras/stt651 Rakov VA , Uman MA ( 2003 ) Lightning: physics and effects . Cambridge University Press, Cambridge Rao NN ( 1993 ) Dust-magnetoacoustic waves in magnetized dusty plasmas . Phys Scr 48 : 363 - 366 . doi: 10 . 1088/ 0031 -8949/48/3/015 Rao NN ( 1995 ) Magnetoacoustic modes in a magnetized dusty plasma . J Plasma Phys 53 : 317 . doi: 10 .1017/ S0022377800018237 Rao NN , Shukla PK ( 1990 ) Dust-acoustic waves in dusty plasmas . Planet Space Sci 38 : 543 - 546 Regener E , Pfotzer G ( 1935 ) Vertical intensity of cosmic rays by threefold coincidences in the stratosphere . Nature 136 : 718 - 719 . doi: 10 .1038/136718a0 Reiners A , Basri G ( 2008 ) Chromospheric activity, rotation, and rotational braking in M and L dwarfs . Astrophys J 684 : 1390 - 1403 . doi: 10 .1086/590073 Rimmer PB , Helling Ch ( 2013 ) Ionization in atmospheres of brown dwarfs and extrasolar planets . IV. The effect of cosmic rays . Astrophys J 774 :108. doi: 10 .1088/ 0004 -637X/774/2/108 Rimmer PB , Helling Ch ( 2015 ) A chemical kinetics network for lightning and life in planetary atmospheres . ArXiv e-prints Rimmer PB , Walsh C , Helling Ch ( 2014 ) Cosmic rays, UV photons, and haze formation in the upper atmospheres of hot Jupiters . In: Booth M , Matthews BC , Graham JR (eds) IAU symposium , vol 299 , pp 303 - 304 . doi: 10 .1017/S1743921313008703 Rodriguez-Barrera MI , Helling C , Stark CR , Rice AM ( 2015 ) Reference study to characterise plasma and magnetic properties of ultra-cool atmospheres . Mon Not R Astron Soc 454 : 3977 - 3995 . doi: 10 .1093/ mnras/stv2090 Rogers TM , Showman AP ( 2014 ) Magnetohydrodynamic simulations of the atmosphere of HD 209458b . Astrophys J Lett 782 : L4 . doi: 10 .1088/2041-8205/782/1/L4 Rosenberg M ( 1993 ) Ion- and dust-acoustic instabilities in dusty plasmas . Planet Space Sci 41 : 229 - 233 . doi: 10 .1016/ 0032 - 0633 ( 93 ) 90062 - 7 Rosenblum E , Garaud P , Traxler A , Stellmach S (2011a) Erratum: ''Turbulent mixing and layer formation in double-diffusive convection: three-dimensional numerical simulations and theory'' /abs/ 2011ApJ... 731 ... 66R ( 2011 , ApJ, 731 , 66). Astrophys J 742 :132. doi: 10 .1088/ 0004 -637X/742/2/132 Rosenblum E , Garaud P , Traxler A , Stellmach S (2011b) Turbulent mixing and layer formation in doublediffusive convection: three-dimensional numerical simulations and theory . Astrophys J 731 :66. doi: 10 . 1088/ 0004 -637X/731/1/66 Route M , Wolszczan A ( 2012 ) The arecibo detection of the coolest radio-flaring brown dwarf . Astrophys J Lett 747 : L22 . doi: 10 .1088/2041-8205/747/2/L22 Rutherford D , McDowell D , Mariotti D , Mahony C , Diver D , Potts H , Bennet E , Maguire P ( 2014 ) Impact of plasma induced liquid chemistry and charge on bacteria loaded aerosol droplets . In: APS meeting abstracts Rycroft MJ ( 1965 ) Resonances of the Earth-ionosphere cavity observed at Cambridge, England. Radio Sci J Res Natl Bur Stand D 69 : 1071 - 1081 Rycroft MJ , Harrison RG ( 2012 ) Electromagnetic atmosphere-plasma coupling: the global atmospheric electric circuit . Space Sci Rev 168 : 363 - 384 . doi: 10 .1007/s11214-011-9830-8 Rycroft MJ , Israelsson S , Price C ( 2000 ) The global atmospheric electric circuit, solar activity and climate change . J Atmos Solar Terr Phys 62 : 1563 - 1576 . doi: 10 .1016/S1364- 6826 ( 00 ) 00112 - 7 ( 2007 ) Ion trap mass spectrometry on a comet nucleus: the Ptolemy instrument and the Rosetta space mission . J Mass Spectrom 42 : 1 - 10 Trakhtengerts VY , Rycroft MJ ( 2008 ) Whistler and Alfve´n mode cyclotron masers in space . Cambridge University Press, Cambridge Treumann RA , Zbigniew K , Parrot M ( 2008 ) Physics of elelctric discharges in atmospheric gases: an informal introduction . In: Leblanc F, Aplin KL , Yair Y , Harrison RG , Lebreton JP , Blanc M (eds) Planetary atmospheric electricity , Space Science Series of ISSI , vol 137 , p 133 Tsai SM , Dobbs-Dixon I , Gu PG ( 2014 ) Three-dimensional structures of equatorial waves and the resulting super-rotation in the atmosphere of a tidally locked hot Jupiter . Astrophys J 793 :141. doi: 10 .1088/ 0004 -637X/793/2/141 Van Eaton AR , Muirhead JD , Wilson CJN , Cimarelli C ( 2012 ) Growth of volcanic ash aggregates in the presence of liquid water and ice: an experimental approach . Bull Volcanol 74 : 1963 - 1984 . doi: 10 .1007/ s00445-012-0634-9 Ververka J et al ( 2013 ) Return to Comet Tempel 1: overview of Stardust-NExT results . Icarus 222 : 424 - 435 Vidotto AA , Gregory SG , Jardine M , Donati JF , Petit P , Morin J , Folsom CP , Bouvier J , Cameron AC , Hussain G , Marsden S , Waite IA , Fares R , Jeffers S , do Nascimento JD (2014a) Stellar magnetism: empirical trends with age and rotation . Mon Not R Astron Soc 441 : 2361 - 2374 . doi: 10 .1093/mnras/stu728 Vidotto AA , Jardine M , Morin J , Donati JF , Opher M , Gombosi TI ( 2014b ) M-dwarf stellar winds: the effects of realistic magnetic geometry on rotational evolution and planets . Mon Not R Astron Soc 438 : 1162 - 1175 . doi: 10 .1093/mnras/stt2265 von Blohn N , Diehl K , Mitra SK , Borrmann S ( 2009 ) Riming of graupel: wind tunnel investigations of collection kernels and growth regimes . J Atmos Sci 66 : 2359 . doi: 10 .1175/2009JAS2969.1 Vorgul I , Kellett BJ , Cairns RA , Bingham R , Ronald K , Speirs DC , McConville SL , Gillespie KM , Phelps ADR ( 2011 ) Cyclotron maser emission: stars, planets, and laboratory . Phys Plasmas 18 ( 5 ): 056501 . doi: 10 .1063/1.3567420 Walsh C , Nomura H , Millar TJ , Aikawa Y ( 2012 ) Chemical processes in protoplanetary disks. II. On the importance of photochemistry and X-ray ionization . Astrophys J 747 :114. doi: 10 .1088/ 0004 -637X/ 747/2/114 Wasilewski P , Dickinson T ( 2000 ) Aspects of the validation of magnetic remanence in meteorites . Meteorit Planet Sci 35 : 537 - 544 . doi: 10 .1111/j.1945- 5100 . 2000 .tb01434.x Weidenschilling S ( 1997 ) Production of chondrules by lightning in the solar nebula? Not so easy! In: Lunar and Planetary Science Conference , vol 28 , p 1515 Whipple EC ( 1981 ) Potentials of surfaces in space . Rep Prog Phys 44 : 1197 - 1250 . doi: 10 .1088/ 0034 -4885/ 44/11/002 Whipple FJW , Scrase FJ ( 1936 ) Point discharge in the electric field of the Earth . Geophys Mem Met Off Lond 38 : 1 Williams ER , McNutt SR ( 2005 ) Total water contents in volcanic eruption clouds and implications for electrification and lightning . In: Pontikis C ( ed) Recent progress in lightning physics . Research Signpost, Kerala, pp 81 - 94 Williams PKG , Cook BA , Berger E ( 2014 ) Trends in ultracool dwarf magnetism. I. X-ray suppression and radio enhancement . Astrophys J 785:9 . doi: 10 .1088/ 0004 -637X/785/1/9 Wilson CTR ( 1906 ) On the measurement of the Earth-air current and on the origin of atmospheric electricity . Proc Camb Philos Soc 13 : 363 - 382 Wilson CTR ( 1921 ) Investigations on lightning discharges and on the electric field of thunderstorms . R Soc Lond Philos Trans Ser A 221 : 73 - 115 . doi: 10 .1098/rsta. 1921 .0003 Wilson CTR ( 1929 ) Some thundercloud problems . J Frankl Inst 208 : 1 - 12 Witte S , Helling Ch , Barman T , Heidrich N , Hauschildt PH ( 2011 ) Dust in brown dwarfs and extra-solar planets . III. Testing synthetic spectra on observations . Astron Astrophys 529 : A44 . doi: 10 .1051/ 0004 - 6361/ 201014105 Wolszczan A , Frail DA ( 1992 ) A planetary system around the millisecond pulsar PSR 1257 ? 12 . Nature 355 : 145 - 147 . doi: 10 .1038/355145a0 Wu CS , Lee LC ( 1979 ) A theory of the terrestrial kilometric radiation . Astrophys J 230 : 621 - 626 . doi: 10 . 1086/157120 Xu W , Celestin S , Pasko VP ( 2012 ) Source altitudes of terrestrial gamma-ray flashes produced by lightning leaders . Geophys Res Lett. doi:10.1029/2012GL051351 Yair Y ( 2008 ) Charge generation and separation processes . In: Leblanc F, Aplin KL , Yair Y , Harrison RG , Lebreton JP , Blanc M (eds) Planetary atmospheric electricity . Space Science Series of ISSI , vol 137 , p 119 Yair Y ( 2012 ) New results on planetary lightning . Adv Space Res 50 : 293 - 310 . doi: 10 .1016/j.asr. 2012 . 04 . 013

This is a preview of a remote PDF:

Christiane Helling, R. Giles Harrison, Farideh Honary, Declan A. Diver, Karen Aplin, Ian Dobbs-Dixon, Ute Ebert, Shu-ichiro Inutsuka, Francisco J. Gordillo-Vazquez, Stuart Littlefair. Atmospheric Electrification in Dusty, Reactive Gases in the Solar System and Beyond, Surveys in Geophysics, 2016, 705-756, DOI: 10.1007/s10712-016-9361-7