# Climatology of transient luminous events and lightning observed above Europe and the Mediterranean Sea

Enrico Arnone, József Bór, Olivier Chanrion, Veronika Barta, Stefano Dietrich, Carl-Fredrik Enell, Thomas Farges, Martin Füllekrug, Antti Kero, Roberto Labanti, et al.

► To cite this version:

Enrico Arnone, József Bór, Olivier Chanrion, Veronika Barta, Stefano Dietrich, et al. Climatology of transient luminous events and lightning observed above Europe and the Mediterranean Sea. Surveys in Geophysics, 2020, 41 (2), pp.167-199. 10.1007/s10712-019-09573-5. hal-03239746

HAL Id: hal-03239746

https://hal.science/hal-03239746v1

Submitted on 9 Aug 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Climatology of transient luminous events and lightning observed above Europe and the Mediterranean Sea

This is a post-peer-review, pre-copyedit version of an article published in Surveys in Geophysics. The final authenticated version is available online at: https://doi.org/10.1007/s10712-019-09573-5

Enrico Arnone · József Bór · Olivier Chanrion · Veronika Barta · Stefano Dietrich · Carl-Fredrik Enell · Thomas Farges · Martin Füllekrug · Antti Kero · Roberto Labanti · Antti Mäkelä · Keren Mezuman · Anna Odzimek · Martin Popek · Marco Prevedelli · Marco Ridolfi · Serge Soula · Oscar van der Velde · Yoav Yair · Przemyslaw Zoladek · Torsten Neubert

Received: date / Accepted: date

Enrico Arnone
Dipartimento di Fisica, Università di Torino, Italy
Tel: +39 0116707261
E-mail: enrico.arnone@unito.it

Enrico Arnone · Stefano Dietrich
Istituto di Scienza dell'Atmosfera e del Clima - CNR, Bologna, Italy

József Bór · Veronika Barta
Research Centre for Astronomy and Earth Sciences, GGI, Hungarian Academy of Sciences, Sopron, Hungary

Olivier Chanrion · Torsten Neubert
Technical University of Denmark, DTU Space, Kongens Lyngby, Denmark

Carl-Fredrik Enell
EISCAT Scientific Association, Kiruna, Sweden

Thomas Farges
CEA, DAM, DIF, F-91297 Arpajon, France

Martin Füllekrug
Department of Electronic and Electrical Engineering, University of Bath, U.K.

Antti Kero
Sodankylä Geophysical Observatory, University of Oulu, Sodankylä, Finland

Roberto Labanti
Italian Meteor and TLE Network (IMTN), Bologna, Italy
**Abstract** In 1999, the first sprites were observed above European thunderstorms using sensitive cameras. Since then, EuroODORE campaigns have been conducted to observe sprites and other transient luminous events (TLEs), expanding into a network covering large parts of Europe and coastal areas. In 2009 through 2013, the number of optical observations of TLEs reached a peak of 2000 per year. Because of this unprecedented number of European observations, it was possible to construct a climatology of 8394 TLEs observed above 1018 thunderstorm systems, and study for the first time their distribution and seasonal cycle above Europe and parts of the Mediterranean Sea. The number of TLEs per thunderstorm was found to follow a power law, with less than 10 TLEs for 801 thunderstorms and up to 195 TLEs above the most prolific one. The majority of TLEs were classified as sprites, 641 elves, 280 halos, 70 upward lightning, 2 blue jets and 1 gigantic jet. The climatology shows intense TLE activity during summer over continental areas, and in late autumn over coastal areas and sea. The two seasons peak respectively in August and November, separated by March and April with almost no TLEs, and a relative minimum around September. The observed TLE activity, i.e. mostly sprites, is shown to be largely consistent with lightning activity, with a 1/1000 of observed TLE to lightning ratio in regions with most observations. The overall

Antti Mäkelä

Finnish Meteorological Institute, Helsinki, Finland

Keren Mezuman

Earth and Environmental Sciences, Columbia University, New York, NY, USA

NASA Goddard Institute for Space Studies, New York, NY, USA

Anna Odzimek

Institute of Geophysics, Polish Academy of Sciences, Warsaw, Poland

Martin Popek

Department of Space Physics, station Nýdek, Institute of Atmospheric Physics, CAS, Prague, Czech Republic

Marco Prevedelli · Marco Ridolfi

Department of Physics and Astronomy, University of Bologna, Italy

Marco Ridolfi

Istituto Nazionale di Ottica - CNR, Firenze, Italy

Serge Soula

University of Toulouse/CNRS, Laboratory of Aerology, Toulouse, France

Oscar van der Velde

Universitat Politecnica de Catalunya, Spain

Yoav Yair

School of Sustainability, Interdisciplinary Center (IDC) Herzliya, Israel

Przemysław Zoladek

Comets and Meteors Workshop, Warszawa, Bartycka, Poland
behavior is consistent among individual years, making the observed seasonal cycle a robust general feature of TLE activity above Europe.

**Keywords** Thunderstorms, lightning, transient luminous events, ground-based observations, Europe, climatology.
## Article Highlights

Detailed climatological study of TLEs and lightning above Europe produced for the first time

Seasonal cycle of TLEs and lightning with two active seasons in summer over land and autumn over sea and coastal areas

Power law distribution of number of TLEs per thunderstorm

# 1 introduction

## 1.1 Transient luminous events

Three decades ago, a test low-light camera recorded a sprite Franz et al. (1990), a spectacular discharge extending for tens of km above a thunderstorm. It was the first discovery of a whole family of upper atmosphere electrical processes, collectively known as transient luminous events (TLEs – see reviews by Rodger, 1999; Neubert, 2003; Füllekrug et al., 2006; Neubert et al., 2008; Pasko et al., 2011). TLEs occur in the region between the upper troposphere and the mesopause, between the top of thunderclouds and the lower ionosphere. They are the visible manifestation of the electrical impact of thunderstorms onto the above atmosphere. They can be made of streamers, weakly ionized plasma channels (Petrov and Petrova, 1999; Ebert et al., 2006; Luque and Ebert, 2009), roughly up to 70 km altitude; or be large diffused patches of enhanced ionization at higher altitude, where the dielectric relaxation timescale becomes comparable with that of dissociative attachment (Pasko et al., 1998). Over the years, TLEs have been distinguished in specific classes: Blue jets (Wescott et al., 1995; Boeck et al., 1995; Petrov and Petrova, 1999) are fountain-like streamers directly expanding from their leader core, and injected from the thundercloud top up to 40 km altitude towards the ionosphere (Krehbiel et al., 2008; Pasko, 2008). Sprites (Sentman and Wescott, 1993; Lyons, 1994; Sentman et al., 1995) are luminous discharges that initiate at about 70 km altitude (Stenbaek-Nielsen et al., 2010) as a consequence of the transient quasi-electrostatic field induced by large positive cloud-to-ground (+CG) lightning strokes (Pasko et al., 1997). They extend downwards to 40 km as streamers and upwards to 90 km altitude as diffuse emission, and are tens of km wide (Stenbaek-Nielsen et al., 2000). Sprite halos occur as downward moving diffuse glow at 70–80 km altitude, often accompanying sprite events but in some occasions alone (Barrington Leigh et al., 2001; Wescott et al., 2001). Elves appear as horizontally expanding diffuse emission rings of a few hundreds km diameter at about 90 km altitude, when the electromagnetic pulse of a triggering lightning stroke hits the lower ionosphere (Boeck et al., 1992; Fukunishi et al., 1996; Inan et al., 1997). Much rarer are gigantic jets (GJs), which join together apparent features of blue jets and sprites, and cause a direct connection of the thunderstorm to the lower ionosphere (Pasko et al., 2002; Su et al., 2003; Cummer et al., 2009; van der Velde et al., 2010a). Depending on
the relaxation timescales at the altitude of occurrence, TLEs last from several hundreds of milliseconds (jets) down to a few milliseconds (elves). Sprites are the most commonly observed kind of TLE from ground based video observations, with global occurrence rate of about 1–3ADELMIN-1 (Ignaccolo et al., 2008, and references therein). Satellite observations (Chen et al., 2008) confirmed a global occurrence rate of sprites of about 1 per minute, a similar rate for halos, whereas the much higher efficiency in observing elves led to an estimate of about 30 ADELMIN-1 globally. A high elve/sprite ratio (6:1) was found also by adopting data from ground based photometers, further pointing to a selection bias towards sprites in ground based video observations (Newsome and Inan, 2010).

## 1.2 Eurosphere: a review of European TLE observations

Europe is a unique region for studies of TLEs under extremely varying conditions. The first images of European TLEs were captured by chance in 1999 over the Balkans from an aircraft based camera during the 999 Leonids-MAC airborne campaign (Gardner and Taylor, Second Leonids-MAC workshop, Tel Aviv, 2000). The first dedicated observations were obtained in 2000 with a camera installed at the Observatoire Midi-Pyrénées, located at Pic du Midi in the French Pyrenees (Neubert et al., 2001). The authors recorded 40 sprites over the Alps and Southern France, in connection with cold fronts coming from the Atlantic. In the following years, Eurosphere campaigns were conducted during summer leading to over 700 TLE images being captured in the period from 2000 to 2008 (Neubert et al., 2005; Chanrion et al., 2007; Arnone et al., 2008; Neubert et al., 2008). The network involves also several non-optical correlative measurements widely adopted for further insight into the physics of the observed TLEs: a broad range of correlative measurements including ELF-VLF radio and infrasound recording (Price et al., 2002; Greenberg and Price, 2004; Liszka and Hobara, 2006; Neubert et al., 2008; Farges and Blanc, 2010; Sátori et al., 2013). These campaigns allowed a great number of detailed studies using specific European TLE events and thunderstorms. This includes investigations of radio signatures of TLEs and thunderstorm-induced effects onto the atmosphere (Haldoupis et al., 2004; Mika et al., 2005; Haldoupis et al., 2006; Mika et al., 2006; Farges et al., 2007; Greenberg et al., 2009; Haldoupis et al., 2010; NaitAmor et al., 2010; Füllekrug et al., 2010, 2011); investigations on infrasound signatures of sprites (Farges et al., 2005; Ignaccolo et al., 2008; Farges and Blanc, 2010); meteology of TLEs and of their producing thunderstorm (Ignaccolo et al., 2006; van der Velde et al., 2006; Ganot et al., 2007; Bór et al., 2009; Iwański et al., 2009; Savtchenko et al., 2009; Vadislavsky et al., 2009; Yair et al., 2009; Soula et al., 2010; Mäkelä et al., 2010; van der Velde et al., 2010b; Bór et al., 2018); the morphological aspects of various LoESS types (Bór, 2013). A general collection of research related to Eurosphere was included in Füllekrug et al. (2006) and presented by Neubert et al. (2008).

Eurosphere has since then expanded to become a network that joins the observational activities of tens of observers across Europe and the Mediterranean
sea, exceeding 1000 observations per year. In particular, a coordination effort
over 2009 through 2013 led to the production of a first database of observa-
tions with a broad coverage of regions over Europe and the Mediterranean sea.
The increased number of observations was accompanied by the first detection
of extremely rare phenomena such as the first GJ observed over Europe (van
der Velde et al., 2010a; Kulak and Młynarczyk, 2011; Neubert et al., 2011),
detailed analysis of individual sprite-producing thunderstorms (Soula et al.,
2014, 2015, 2017), of sprite parent lightning (van der Velde et al., 2014), and
elve statistics of a specific region (van der Velde and Montanyà, 2016), high
speed recording of sprites (Montanyà et al., 2010), sprites signatures in radio
waves (Farges and Blanc, 2011; Füllekrug et al., 2013b; Młynarczyk et al.,
2015), long-lasting TLE signatures in the ionosphere induced by rare very
strong lightning (Haldoupis et al., 2012, 2013). More recent studies included
also detailed impact of sprite-producing thunderstorms on the atmosphere
above with consequent relativistic acceleration of electrons (Füllekrug et al.,
2013a) in association with model studies (Chanrion and Neubert, 2010), im-
pact of sprite-producing thunderstorms on the lower ionosphere (especially on
the sporadic E layer) over Central European region (Barta et al., 2017) using
ionosonde data, or may lead to joint studies also with detection of terrestrial
gamma-ray flashes over the Mediterranean Sea (Gjesteland et al., 2015).

1.3 A climatological perspective

Such a large number of observations allows for the first time detailed clima-
tological studies to be performed over Europe and the Mediterranean sea,
therefore studying the geographical distribution and variability of mean TLE
occurrence over a time span of several years rather than looking at individ-
ual events. A Southern Mediterranean example was presented by Yair et al.
(2015) and an introduction to the climatology presented here was published
by Arnone and Dinelli (2016). Despite being inhomogeneously distributed,
the number of ground-based observed sprites largely exceeds that acquired
globally from satellites over equivalent periods of time (Chen et al., 2008).
They give therefore an essential contribution to climatological studies, which
remain fundamental for understanding the global role of TLEs in the atmo-
sphere (see e.g., reviews by Pasko, 2010; Pasko et al., 2011): e.g. answering
questions on the spatial and temporal distribution of their occurrence and
of their impact; opening the way to comparison with other climatological at-
mospheric and climate data for which one-to-one analysis is not meaningful;
allowing calibration of climatologies based on non-optical measurements, such
as through Schumann resonances and detection at extremely low frequency
(ELF, Füllekrug and Reising, 1998; Whitley et al., 2011). The climatology
of TLEs with associated details for individual TLE types and relationship to
lightning opens the way to further interpretation in future studies dealing e.g.
the meteorology and the global electric circuit.

A branch of European TLE-related activities that has developed over the
past years and will greatly benefit from these large observational samples is
that studying the impact of TLEs on the atmosphere. This include model and laboratory studies of the discharges themselves (e.g., Luque and Ebert, 2009; Ebert et al., 2010); modeling of their emissions and impact onto the atmosphere and its chemistry at local or global level (Enell et al., 2008; Gordillo-Vázquez, 2008; Arnone et al., 2008; Gordillo-Vázquez et al., 2011; Parra-Rojas et al., 2013; Neubert and Chanrion, 2013; Arnone et al., 2014; Winkler and Notholt, 2014, 2015; Parra-Rojas et al., 2015; Pérez-Invernón et al., 2018a, 2019); observational studies linking chemistry to TLE occurrence (Arnone et al., 2008, 2009; Arnone and Hauchecorne, 2012; Arnone and Dinelli, 2016); a number of TLE parametrized distributions have been adopted by several of these modeling and observational studies based on regional to global lightning observations. Furthermore, investigations on the role of sprites in the global electric circuit were also brought forward (Füllekrug et al., 2006; Rycroft et al., 2007; Rycroft and Odzimek, 2010; Rycroft and Harrison, 2011). More recently, a capability of investigating TLEs with high resolution spectroscopy was also developed (Gordillo-Vázquez et al., 2018) in association with optical observations of sprites, and methods for the spectroscopic diagnostic of elves and halos were developed for space-based instruments (Pérez-Invernón et al., 2018b).

The growth of this field of research in Europe has brought to the development of the Atmosphere-Space Interaction Monitor (ASIM) that was launched on the 2nd of April, 2018, and installed on an external platform of the Columbus Module of the International Space Station. The main objectives of the mission are related to thunderstorm activity by observing associated emissions in the near UV, near-infrared, X-ray and gamma-ray bands. ASIM consists of two main instrument modules, the Modular Multi-spectral Imaging Array (MMIA) with two cameras and three photometers observing in the 180-230 nm, 337 nm and 777.4 nm bands and the Modular X-ray and gamma-ray Sensor (MXGS) observing photons of energy in the range 15 keV to 20 MeV with imaging capability (Neubert, 2009). A further European mission dedicated to the study of thunderstorms is the TARANIS satellite (Tool for the Analysis of RAdiations from lightNINGS and Sprites) which is planned to be launched in 2020 carrying a set of instruments recording thunderstorm emission in the optical, X-ray, gamma-ray and radio bands (Blanc et al., 2007; Lefeuvre et al., 2008). Both space missions will exploit the coordinated ground-based observation systems and the knowledge on TLE distribution and variability discussed in this study.

In this paper, we present the distribution and seasonal cycle of TLE observations from the Euro pórite network and partners during 2009 through 2013, the period with the widest coverage of European observations to date. The paper aims both at giving an overview of the observational capabilities of the network and exploiting its results to obtain for the first time a climatology of TLEs over Europe. The structure of the paper reflects these aims: the coordinated Euro pórite instrumentation and observations are presented in Sect. 2. The data analysis is described in Sect. 3. The climatology is presented in Sect. 4 and discussed in Sect. 5. Conclusions are given in Sect. 6.
## 2 Eurosphere instrumentation

Several partners contribute to the so called Eurosphere network which provided the observations studied in this work. Eurosphere is an umbrella identifying a coordinated observational effort exploiting several optical cameras permanently installed at strategic locations or temporary deployed during field campaigns throughout Europe and coastal areas. The main optical systems involved in the network are shown in Fig. 1 with their estimated coverage (see Sec. 3.1). A number of further optical systems contribute with sporadic observations through their regional networks (see below). As shown, most of Central and Southern Europe is well covered by observations, whereas there is as yet limited coverage of Eastern and Northern Europe. The Mediterranean sea and coastal areas are covered in the Northern parts by systems in Spain, France and Italy, and in the Southern parts by Israel. An overview of the optical systems involved in this study (i.e. active in 2009 to 2013) and their current status is given here below. An overall description of Eurosphere activities and early campaigns can be found in Neubert et al. (2008).

### 2.1 Optical systems

#### 2.1.1 Southern Europe and the Mediterranean Sea

The first operating European sprite-watch system was installed in Southern France at Pic du Midi (42.94° N, 0.14° E, 2877 m) in the French Pyrenees to conduct campaign operations in 2000 and 2003 (Neubert et al., 2001; Neubert et al., 2005). It consisted of remotely controlled optical cameras sensitive to low-light conditions and a photometer. This camera system was associated to other similar systems for summer and fall campaigns, especially in 2005 and 2006 with a camera at Puy de Dome (45.77° N, 2.962° E, 1465 m) and in 2008 on Mount Corona (42.46° N, 8.92° E, 2144 m) in Corsica, and at the Calar-Alto Observatory (37.22° N, 2.58° E) (Neubert et al., 2005; Chanrion et al., 2007; Arnone et al., 2008; Neubert et al., 2008; Soula et al., 2010). See Fig. 1. A new low-light and high-resolution charge-coupled device (CCD) camera (Watec 902H) mounted on a pan-tilt unit remotely controlled through the Internet, was installed in 2009 at Pic du Midi and continuously operated. It is equipped with a 12 mm f/0.8 lens with a 31° field of view (FOV). The typical maximum range of such a high mountain system is 800 km, so that this system commonly observes TLEs over the southeastern France, the Alps, the western Mediterranean Sea and a large part of Spain (Soula et al., 2010, 2014, 2015, 2017; van der Velde et al., 2010a, 2014). Other similar cameras are located in southern France at lower altitude, in Lannemezan (43.13° N, 0.37° E, 592 m) since 2007, in Clermont-Ferrand (45.76° N, 3.11° E; 400 m) since 2010, in Rustrel (43.94° N, 5.48° E, 1020 m) since 2011. These four cameras can be operated throughout the year. In particular, during the period considered for this study, the Pic du Midi system was active in 2010 (mid-July to end of August), 2011 (mid-May to end of November) and 2013 (mid-
August to end of September); the Corsica system in 2009 (August to mid-November), 2010 (mid-July to mid-October) and 2011 (second half of July); the Calar Alto system in 2009 (mid-April to mid-November) and 2011 (end of September to mid-December). In September 2018, a high speed remote camera system was mounted at the top of the Laboratoire Souterrain à Bas Bruit (LSBB) in Rustrel (43.93° N, 5.49° E) and another system at the Observatoire MidiPyrnes is planned for installation in 2019.

Observations taken from Spain were based on a remotely controlled cam-
era installed at Sant Vicenç de Castellet (41.67° N, 1.85° E) and additional
deployable cameras including a high-speed system (see Fig. 1). The remotely
controlled system has an observation range of about 150 to 450 km centered
at the camera location, with minimum distance rarely below 100 km (about
35° high in the sky) and reaching a maximum distance of about 600 km. The
viewing had partial limitation towards Southern France because of hills block-
ing the lowest 12°. The system was moved to Castellgalí (41.67 ° N, 1.83° E)
in November 2013 where it continues operations. The cameras were operated
continuously throughout the year, with observation taken when viewing con-
ditions were clear and manually operated whenever there were storms.

Observations from Italy and Switzerland were recorded by the Italian Me-
teor and TLE Network (IMTN), which groups over 30 cameras with fixed
pointing and a few automated pointing cameras. IMTN fully covers central
and Northern Italy, and adjacent regions, and has some limited sensitivity
down to Southern Italy. Most cameras are close to sea level and close to cities,
so that the observation range of each camera is limited to a few hundred km.
Because of multiple cameras covering the same areas from different location,
it is often possible to triangulate. The most active stations (i.e. with more
cameras and continuity of observations, resulting in more observations taken)
are installed at Ferrara (44.83° N, 11.57° E), Tortoreto (42.66° N, 13.67° E),
Contigliano (42.41° N, 12.76° E), Bologna/Medicina (44.50° N, 11.26° E),
Lugano/Gnosca (46.12° N, 8.84° E), Cuneo (44.39° N, 7.48° E) and Scan-
dicci (43.74° N,11.08° E) The southernmost camera is placed in Casamassima
(41.03° N, 16.82° E) since 2012. In 2010, a research camera with a remotely
controlled automated TLE detection system was mounted at the Italian Cli-
mate Observatory on Mount Cimone (44.19° N, 10.70° E, 2165 m), Modena,
Italy. The system was moved to the Loiano Observatory (44.39° N, 11.19 ° E)
in early 2013 and continues to be operational. The overall coverage from Italy
includes adjacent seas, parts of central Europe and the Balkans (see Fig. 1).
The cameras are operated throughout the year. The network of overlapping
camera viewing demonstrated to be extremely efficient in capturing TLEs,
leading for example to the only GJ captured to date over Europe (van der
Velde et al., 2010a). Despite shortages in the viewing conditions and continu-
ity of the operation of individual cameras, these regions can be considered to
be largely continuously monitored.

Covering the western Mediterranean Sea and adjacent regions observations
were conducted from Israel by ILAN team. The TLE observing system is de-
ployed on the rooftop of the Geophysics department of Tel Aviv University
and is comprised of two panchromatic CCD cameras, mounted on a remotely
controlled pan-and-tilt unit (Ganot et al., 2007) (see Fig. 1). Observations
are conducted during the thunderstorm active period of September-May ev-
ery year, with a line of sight stretching up to Cyprus and Southern Turkey
(Yair et al., 2009). For specific storm events, an alternative site at the Wise
astronomical observatory (30.60° N, 34.76° E) in the Negev desert was used,
extending the coverage to the Nile delta and beyond.

2.1.2 Central Europe

Organized TLE observations in Central Europe started in 2007. An optical
detection site has been set up in Sopron, on the North-Eastern part of Hungary
(47.69° N,16.44° E) with a remotely controlled monochrome analogue video
camera. Due to the vulnerability of the applied pan-tilt unit, the camera used
to be dismounted for the winter and reinstalled in the following year (see
Fig. 1). Observations were taken in May to August, June to August, and
July to October, in 2009, in 2010, and in 2011, respectively. Further details
of TLE-related observations in Hungary can be found in Sátori et al. (2013).
Recently, another TLE observation site has been set up in southern Hungary
at Baja (46.18° N, 19.01 ° E) with similar instrumentation to that in Sopron
on the Noth-Eastern part of the country. The new site covers central Italy, the
Adriatic Sea, and the middle and the North-Western part of the Balkans.

Observations from Czech Republic started in May 2011 with a camera in-
stalled at at Nydek (49.67° N, 18.77° E) and continued thereafter all years
round with only very minor breaks during the period of this study (5/7 to
21/7 in 2011, 5/11 to 26/11 in 2012 and 1/11 to 14/11 in 2013). See e.g. (Mly-
narczyk et al., 2015). Observed TLEs occurred between May and November.
The contribution from this region reached about 700 TLEs, covering Czech
Republic, Slovakia, Germany, Austria, Poland, Hungary, Ukraine and extend-
ing to Italy, Slovenia, Croatia and the Adriatic Sea with distances from the
observer ranging between 100 and 700 km.

First TLE research observations in Poland have been carried out during a
two-week field campaign at Mount Sniezka (50.74° N, 15.74° E), the highest
peak of the Sudetest, in July 2007 organized as part of EuroERMite (Odzimek
et al., 2008; Iwański et al., 2009). Since summer 2011, research observations
have also been made sporadically from Gliwice (50.28° N, 18.65° E) with
azimuth 180-300° and from Swider (52.01° N, 21.39° E) with azimuth 300-50°.
These observations have been carried out using low-light sensitive cameras
pointed manually, covering parts of Germany, Austria, Hungary, the Czech
Republic and Poland (see Fig. 1), as far as Southern Lithuania from Swider.
Large parts of this area overlap with those covered by systems in Hungary
and Czech therefore allowing triangulation. It has also to be noted that tens
of sprites each year have been observed from Poland since 2009 by the Polish
Fireball Network (PFN). PFN operates cameras at fixed direction over the
territory of Poland with stations disseminated over central, North-West and
South-West of Poland, and two stations in the East, with limited and variable
viewing directions for all stations. PFN stations with the longest detection
period or the most of the observed events are Gniewowo (54.58 N, 18.30 E) which made the first PFN detection in 2009, Twardogora (51.37 N, 17.46 E), Nowe Miasto Lubawskie (53.43 N, 19.59 E), Blonie (52.19 N, 20.62 E), and Chelm (51.16 N, 23.50 E).

### 2.1.3 Northern Europe

Optical observations of TLEs in Northern Europe have to deal with the basic problem that the sky is seldom dark during summer, when most Northern European thunderstorms occur. A station has been operated at Esrange, Kiruna to look for sprites over winter thunderstorms over the Atlantic coast of Norway. In Finland there is an extensive network of amateur astronomers and storm chasers related to the Astronomical Association Ursa running automatic detection software for observing bolides and other bright events in the twilight and night sky. In 2009, the northernmost European sprite and elve were recorded (Mäkelä et al., 2010). In 2009-2013, a total of 25 TLEs were observed (24 sprites, 1 blue jet).

## 2.2 Lightning detection

Lightning detections from national networks were used both as guidance for pointing direction during the observations, and in order to geolocate the TLEs with their parent lightning in post-processing of TLE data. In particular, lightning data came from the VLF/LF lightning detection network (LINET, Betz et al., 2009) network in Italy, Spain, Hungary, Poland and Germany. Data from Meteorage and from the European Cooperation for Lightning Detection (EUCLID) were adopted in France and nearby areas. Lightning data from the community lightning detection network BlitzOrtung were used by automated pointing cameras in Italy and Czech Republic. In Poland lightning detection was acquired also from the Institute of Meteorology and Water Management lightning network (PERUN) and Central Lightning Detection Network (CELDN). In Finland, the lightning location data is from the Nordic Lightning Information Systems (NORDLIS, Mäkelä et al., 2010). Additionally, satellite infrared images and local radar images are used both for pointing and data geolocation.

The analysis of lightning over Europe was performed adopting data from the World Wide Lightning Location Network (WWLLN). Detection of lightning in WWLLN is based on the consistent recording of VLF band spheric signals from lightning strokes at 5 stations (as a minimum) of the network. WWLLN processes the signals from its receivers centrally and uses the time-of-arrival method to calculate geographical location and the time of occurrence of the source of the detected signals with 10 km and microsecond average accuracy, respectively (Rodger et al., 2006; Rudlosky and Shea, 2013). Since 2009, WWLLN detectors covering continental Europe and the Mediterranean Sea are placed in Portugal, Hungary, Israel and Northern Finland, whereas a new
detector was activated in 2012 in the U.K. The switch on of the U.K. detector
may lead to an increase in the detection rate since 2012.

We note that the detection efficiency in the WWLLN system is low for
individual lightning strokes when compared to local or regional lightning de-
tection networks (Abarca et al., 2010). This limitation, however, corresponds
to weak discharges. The overall detection efficiency ranges between 15 and
30% depending on geographical location, biased towards strong lightning. It
is known that the appearance of sprites is predominantly triggered by intense
+CG lightning (e.g., Neubert et al., 2008) and that sprite halos and elves are
produced by high peak current lightning strokes (Williams et al., 2012; van
der Velde and Montanyà, 2016). Unfortunately, WWLLN does not report peak
currents nor the type (CG/IC) of the stroke, so its full lightning database, as
it is considered in this study, contains many lightning strokes which cannot
be associated with TLEs. Nevertheless, the good homogeneity of the dataset
in 2009 to 2013 over Europe and the Mediterranean Sea enables it to be used
for comparison with TLE observations as an aid in the interpretation of the
findings.

3 Observations and data analysis

Over 2009 to 2013, thousands of observations were recorded by the network described in the previous section. Since the range of view of TLEs for these low-light cameras can reach several hundreds of kilometers, TLEs were recorded over large parts of Europe and the Mediterranean sea, through the complex coverage discussed in the following section. All observations are made of optical images, consisting of one or two consecutive video frames capturing one or more TLEs, and one or more elements of a sprite or sprite groups, if present. Given the large number of observations, a very vast family of different types of sprites were registered, including carrot, angel and column sprites, with a variety of shapes and at all levels of detail. Because of detection difficulties, a much lower level of detail characterizes the elve and halo images. The details of TLE images were not analyzed in this work favoring the overall climatological aspects of the distribution through the construction of a database (Sect. 3.2) and binning (Sect. 3.3) of the observations.

3.1 Observational coverage

Figure 1 shows the location of the most representative optical systems that con-
tributed with observations during 2009-2013 (white crosses) with an estimate
of their viewing coverage (colored areas). A main issue with ground-based opti-
cal TLE observations is the unevenness of the distribution of the cameras, and
the continuous changes of their viewing due to either atmospheric processes
or experimental changes. Although the main observation hot spots are known
(i.e. the areas covered by cameras capturing most TLEs), the detection effi-
ciency is affected in a way that is not possible to characterize correctly for the
largest part of the available data, for which detailed historical conditions (e.g. viewing conditions, pointing direction, operating period, observational choice of the observer when more thunderstorms are in view) were not recorded. In order to take into consideration these shortages and avoid introducing analysis biases between different cameras, we adopted and brought forward two approaches in parallel.

Firstly, an approximate coverage was evaluated adopting the same discretization as for the data analysis (0.5°, see Sec. 3.3) and based on the viewing of the most representative camera systems (as introduced above). These systems were selected as main systems of each contributing research group or regional network and on the basis of the continuity of their operations, area covered (e.g. if covering areas that no other cameras are reaching), and reported number of TLEs. Each of these cameras was characterized by a circular area of coverage around its location, with a radius of 300 km (e.g., in a city) to 800 km (e.g., on a high mountain) depending on the geographical location and characterization by the observer. As a first approximation, a circular area was assumed for all cameras with no consideration of their actual viewing direction, assigning a score of 1 to each geographical bin within the circular area of a camera. The score was lowered to 0.5 in the closest 100 km range of the camera (where the high zenith angle may prevent observing TLEs) and in the furthest 100 km ring (200 km for high mountain sites) (because of the expected drop in detection efficiency due to great distances). A weight ranging 0 to 1 was then applied to the scores of each camera depending on the continuity of operations: fraction of the 5 years 2009–2013 and fraction of individual year (season). By this characterization, a bin within the viewing range of a camera operated throughout the year would get a score 1 if operated during the whole 5 years, 0.8 if operated during 4 years only. The weight was tuned also depending on known shortages or high efficiency of the camera system (from +0.2 for the most experienced to -0.2 for novel observers, and assigning a default 1 for the main multi-camera systems at Mt Corona, S.V. Castellet, Ferrara and Teramo). The coverage distribution was then calculated integrating the contributions from all cameras, resulting in maximum scores of 5 (i.e., 5 overlapping cameras) indicating the highest probability of detecting TLEs. The resulting distribution reported in Fig. 1 shows the large areas covered and highlights the areas with higher expected observing capability. This coverage was used to interpret the TLE climatology and to simulate the expected TLE distribution based on the lightning distribution. With this approach, the actual detection efficiency of individual cameras remains to be evaluated a posteriori but can be assumed to be partly compensated and similar among all cameras so that the overall geographical distribution can be correctly interpreted.

Secondly, in order to remove the bias produced by the uneven detection efficiency, we studied seasonal changes of the geographical distribution of TLEs in terms of multi-year seasonal averages normalized by the multi-year yearly average (see Sec. 4.2). We also considered monthly averages integrated over the area where most of TLEs were observed (see white contours in Fig. 3),
excluding data from the Eastern Mediterranean Sea in order to obtain a more
homogeneous sample.

## 3.2 The Eurosphere database

The Eurosphere data used in this study include all available 2009–2013 optical images from the broad network of instrumentation described in the previous section. Individual TLE images were collected by the observers into data entries for the Eurosphere database. Each entry consists of the number of TLEs associated to an individual thunderstorm or closely related thunderstorm cells; the geographical area covered by those TLEs; the period of TLE activity; and the TLE type. TLEs recorded in one individual observation (e.g. several sprite elements in one video frame) were typically counted as one, unless they were clearly discernible as separate events. The geographical coverage of the TLEs was estimated by the observers on the basis of the camera field of view, its pointing direction and correlative meteorological data. The geolocation was performed either associating the TLE to its parent lightning stroke obtained through lightning detection networks, whenever possible, or to its parent thunderstorm through the overall lightning and cloud conditions in the region of interest. In the database, the TLE data coverage was reported with a ±1° latitude and longitude uncertainty for most observations, whereas the timing with a ±1 min uncertainty. The ±1° uncertainty reflects both the inaccuracy of the actual spatial extent of small events, and on their geolocation. Observations reported with higher accuracy were rounded to the closest 0.5°, assuming a minimum spatial extent of 1°. An extract of the database with a selection of the most prolific thunderstorms is reported in Table 1. The database with the complete list of observations was made publicly available (Arnone et al., 2019).

## 3.3 Binning of observations

In order to avoid being affected by uncertainties on individual events, we studied the observed TLEs only in terms of climatology, i.e. their mean distribution and variability. Eurosphere data entries were binned over 0.5° latitude and longitude bins, and over 1-month intervals, distributing their contribution among all bins affected by the data entry/thunderstorm area. This implies an oversampling and smoothing was performed of data over a grid finer than the given uncertainty, in order to obtain an increased resolution allowed by higher accuracy observations and overlap of several observations. Overlapping observations of the same events by multiple cameras was either accounted for directly by the observers (this is the case of e.g. the several cameras of the IMTN network), or by removal of the database entry with less counts in case of entries that clearly overlap in space and time. The latter removal was applied only to cases reporting more than 20 TLEs to avoid rejecting cases where
few different TLEs were captured by different cameras over the same region. Seasonal and yearly averages were calculated based on monthly means.

Lightning data from the WWLLN network was also analyzed in a similar fashion, producing a climatology with number of detected strokes within individual geographical bins in a certain season or year. In addition to that, the subset of lightning strokes which occurred during the night was considered separately. Nighttime strokes were selected by considering the time of the local sunset and sunrise. For lightning climatologies, we adopted a discretization of 1° in latitude and longitude in order to be less affected by differences at too fine regional scale and allow an easier interpretation of the TLE distribution main features. In several occasions, we studied lightning data within the limited area of camera coverage or that delimiting the actual TLE observations (see further details below), and performed integration of the overall stroke counts over these regions to obtain monthly means. WWLLN lightning climatologies were used also to evaluate the distribution of TLE to stroke ratio and to simulate the expected TLE distribution based on detected strokes. In these calculations, the size of bins of the TLE climatology was accordingly increased to 1°.

# 4 Results

## 4.1 EuroERMTE observations in 2009–2013: TLE and lightning climatologies

In the years 2009 to 2013, EuroERMTE and partners recorded 8394 TLEs over 1018 thunderstorms. The vast majority of observations were sprites, with 7157 classified events, followed by 641 elves and 280 halos. Also observed were 70 upward lightning processes, 2 blue jets and the first European gigantic jet (the latter recorded in December 2009, west of Corsica – see further details in van der Velde et al., 2010a). The remaining fraction of the events (4%) were reported as unclassified TLEs.

The climatology of TLEs above Europe and Southern Mediterranean Sea for 2009 to 2013 is reported in Fig. 2 as density of observed TLEs (TLEs $\cdot 10^{-3} \text{ km}^{-2} \text{ yr}^{-1}$). All TLE types were merged into the climatology to allow an overall view of the activity of TLE-producing thunderstorms above Europe. The figure reports also individual climatologies for the three types of TLEs giving the largest contributions to the database, i.e. sprites, elves and halos. In a similar way, the 2009–2013 climatologies of WWLLN lightning strokes as total daily counts and nighttime only counts are reported in Fig. 3 (strokes $\text{km}^{-2} \text{ yr}^{-1}$). Since TLEs are observed only during nighttime, the latter should be preferred for comparison, although both maps show relevant information on the geographical distribution of thunderstorm activity over Europe and the Mediterranean sea. In particular, total lightning can be more correctly compared to climate parameters in terms of e.g., temperature, winds and precipitation. To ease comparison, the main area of observed TLEs over Europe is reported in the WWLLN map as reference (see white contours).
The geographical extension of TLE activity closely resembles the estimated observational coverage, therefore supporting the overall viewing range adopted for the cameras. Drops in recorded TLE rates can promptly be associated with decreases in observational coverage. High TLE activity is found in Southern France and around the Balearic Islands, in Italy and adjacent seas, then extending over Austria and the Czech Republic to the North, towards Hungary to the East, and with a separate set of observations around Cyprus. Despite the biases introduced by the location of the observational systems, some main features in the observed TLE climatology can be extracted with consideration of both observational coverage and lightning activity. Within the regions covered by the observations, the TLE main geographical distribution tends to mimic the distribution of thunderstorm activity. This is evident in the northern part of the TLE climatology where large areas with a weak TLE rate are consistent with a drop in lightning activity over vast areas, e.g. in France and Germany, and particularly over Spain contrasting with the adjacent high activity over the Pyrenees. Interestingly, some local drops in TLE rates are consistently seen in lightning, as e.g. over Eastern Italy, and likely Sardinia, a behavior that is related to nighttime lightning and less evidently to total lightning. Further similarities occur in the coastal areas of Northern Spain and Western France, where increased TLE rates are correlated with weakly increased nighttime lightning rates (more visible with a dedicated change of the color scale here not shown) and more evident in total lightning. High TLE rates are recorded at the German-Polish border with partial correlation to lightning activity (see spots of high lightning activity in this region more visible in total lightning), whereas the high activity over Czech Republic is not reproduce by a similar increase in strokes counts. In contrast, the fading of TLE activity in Southern Italy and towards Eastern Europe is due to poor coverage. The poor coverage causes similar low rates over Corsica, North-West Italy and Hungary, where local cameras tend to observe only thunderstorms that occur at a certain distance from the observation spot.

As also evident looking at absolute numbers, sprites dominate the TLE climatology. The dominant role of sprites in our observations is also visible in the maps of individual types of TLEs. The TLE climatology closely resembles that of sprites. On the contrary, elves are predominantly observed over sea and coastal areas where they mostly reinforce features already present in the sprite climatology rather than changing it substantially. The only large change is on the Atlantic Ocean North of Spain where almost no sprites were detected. Halos are observed both over sea and land, again not showing any peculiar features not seen in the sprite climatology. The observation of elves and halos are affected by the much lower detection efficiency found in low-light ground cameras as compared to sprites. Their resulting individual climatologies should be treated with care since they are more sensitive to the different optical system adopted.
### 4.1.1 TLE occurrence rate and TLE to lightning ratio

The peak TLE occurrence rate is close to or exceeds $10^{-3} \text{ km}^{-2} \text{ yr}^{-1}$ in a few hotspots in Southern France, Northern Mediterranean Sea, Italy, the Balkans and the German-Polish border, whereas it is typically around $0.2-0.3 \cdot 10^{-3} \text{ km}^{-2} \text{ yr}^{-1}$ in large adjacent regions. These rates and the main features of the climatology should be largely assigned to sprites, since the other observed TLEs represent a minor fraction of the database. Elves were mostly observed over autumn/winter maritime thunderstorms in the Northern Mediterranean area, reaching rates in the order of $0.02-0.2 \cdot 10^{-3} \text{ km}^{-2} \text{ yr}^{-1}$. Halos peaked at $0.02 \cdot 10^{-3} \text{ km}^{-2} \text{ yr}^{-1}$ in most active regions. Comparison to total lightning rates seen in the lightning climatologies of around $0.2-0.3 \text{ km}^{-2} \text{ yr}^{-1}$ within the same regions (see Fig. 3, left), suggests a factor around 1/1000 in the ratio of observed TLE to detected total lightning on a yearly average (recall the factor $10^{-3}$ in the TLE climatologies). To further inspect this relationship, the TLE to nighttime lightning stroke ratio is shown in Fig. 4 and was evaluated adopting a lower resolution $1^\circ$ climatology of TLEs. Typical values of the ratio are of the order 1 to 10 in regions well covered by observations, and exceed 10 in a few hot spots previously identified, particularly in the area around the German-Polish border where many more TLE per lightning stroke are observed. This ratio greatly changes as a function of season, with average values over the reference region of 0.5/1000 in February, 3/1000 in July and 1.5/1000 in November (the seasonal cycle will be discussed in Sect. 4.2).

Taking into consideration the shortages in the observational coverage, we can evaluate an average TLE to lightning stroke ratio only on those bins which have the highest ratio values. When considering the bins with the top 50% of the ratio values, the mean TLE to nighttime lightning stroke ratio is 2.7/1000; when considering the top 30%, the mean ratio increases to 4.3/1000 (corresponding to a value of 1.3 when total lightning is used in place of nighttime lightning). Note that over the whole climatology, the rate of nighttime to total lightning has a mean value of 0.39 and median value of 0.37, with higher values over coastal areas and sea, and lower ones over continents (not shown).

### 4.1.2 Simulated TLE distribution

Based on the above results, we simulated an expected observed TLE distribution scaling the nighttime lightning climatology (see Fig. 3) by the observational coverage (see Fig. 1) and the estimated TLE to nighttime lightning ratio 4/1000. The observational coverage was normalized to a maximum value of 1. The results are shown in Fig. 4, right panel. The simulated TLE distribution now carries information on both lightning distribution and observational coverage, with maximum values scaled by the adopted TLE to lightning ratio. The main characteristics and typical values seen in the observed TLE climatology are found also in the simulated one: e.g., overall extension and shape of the area covered by the observations, main active regions and the overall values, features of the high rate areas and low rate ones. Also, some peculiarities of
the observed climatologies can now be more easily interpreted on the base of
a simultaneous effect of lightning activity and camera sensitivity. This is the
case of the TLE distributions in Northern Spain and Western France (where
the line features could be better reproduced with a change in the color scale
– not shown) or the drop in central Spain; high rates in the Southern France
coast and close to the Balearic Islands; and drop in the rates over Sardinia.
The simulation reproduces correctly the observed drop over coastal Eastern
Italy, Northern France, and Germany, and can also pick out the behavior of
the distribution observed over Hungary. In contrast, the simulation does not
replicate the behavior of the climatology over Czech Republic, Slovakia and
Poland, where the observational coverage can be further refined. This can be
done possibly in dedicated studies where the actual field of view of each camera
is considered. The simulation also predicts higher than observed rates around
Cyprus because of the local very intense lightning activity, although the overall
area is reproduced suggesting lower weights should be estimated via dedicated
studies.

The comparison of the observed and simulated TLE climatologies was evaluated quantitatively by inspecting the probability distribution functions and main statistics of the two datasets, considering a 1.0 ° lower resolution TLE climatology which was found not to differ in its statistics from the original 0.5 ° resolution climatology. The TLE and simulated climatologies have about 530 and 590 points respectively, with TLE rates greater than zero and lead to very similar probability distribution functions (not shown). Mean TLE rates for the observed climatology are 0.07, 0.10 and 0.21 (medians 0.03, 0.07 and 0.19) ·10⁻³km⁻² yr⁻¹ considering only geographical bins with rates above 0, above 0.01 and above 0.1 ·10⁻³km⁻² yr⁻¹ respectively. In comparison, mean TLE rates for the simulated climatology are 0.11, 0.13 and 0.25 (medians 0.05, 0.07 and 0.20) ·10⁻³km⁻² yr⁻¹ respectively. The slightly higher mean values of the simulated climatology can be reduced adopting a TLE to nighttime lightning rate of 3/1000 to 0.08, 0.10 and 0.21 (medians 0.04, 0.06, 0.17) ·10⁻³km⁻² yr⁻¹, although at the expenses of the median which is also slightly reduced, and of a poorer simulation of rate values over some areas (e.g., Spain and France). Overall, the adoption of a 3.5 or 4/1000 ratio seems to be adequate, supporting a posteriori what was previously calculated based on the top 30% of the geographical bins with highest rates.

## 4.2 TLE seasonal cycle above Europe: changes to the mean geographical distribution

Seasonal averages of the geographical distribution of TLEs for individual years 2009 to 2013 are presented in Fig. 5. Dark grey shades show the overall area of observed TLEs during each year. As for the yearly average, the same seasonal averages are reported also for WWLLN nighttime lightning strokes in Fig. 6. Color grading refers to Fig. 2 and 3, respectively. In terms of seasonal change, TLE activity over Europe is concentrated over the sea in winter, then
moves to the coastal areas in spring substantially fading in intensity, it in-
creases and spreads over the continent in summer, and relocates again over
the sea in autumn. This behavior is very consistent among the 5 years of the
sample. The activity in the Southern Mediterranean Sea is consistent with the
onset of maritime thunderstorms in autumn and winter (see Yair et al., 2015,
for further details). Comparison to Fig. 6 shows a remarkable consistency.
Firstly, the general seasonal oscillation of TLE detection between land and
sea is consistent with thunderstorm activity. Secondly, detailed comparison of
the hotspots of observed TLEs within a certain season corresponds in most
cases to regions having the highest lightning rate (see e.g. the coastal and sea
areas close to Southern France, Balearic Sea and Italy), or, vice-versa, lack
of TLEs corresponds to very low lightning rates (e.g. France and Germany as
mentioned above). There are clear deviations from this general winter over-sea
and summer-land behavior, e.g. in the case of autumn thunderstorms observed
at the border of Czech Republic, Germany and Poland both in 2011 and 2012,
or of the large summer activity observed in 2011 around the Balearic islands.
Consistently, in both cases there are signatures of high lightning activity also
in WWLLN data (note that the maps show the yearly rate per km², so that
rates in individual seasons can largely exceed that calculated over the whole
year dataset).

Especially during its first months, 2009 shows a relatively poorer data
coverage as compared to the following years. Despite this difference, the main
features in the distribution of TLEs are repeated over the years, making it
meaningful to average the 5 years together as performed in Fig. 2. In order to
remove the bias induced by the clustering of optical observation systems over
specific regions (see Fig. 1), we normalized the seasonal distributions of the
observed TLE activity by the 2009-2013 yearly average (Fig. 2). Results of the
normalization are reported in Fig. 7 as percent component of the yearly mean.
The magnitude of the seasonal changes in the distribution of TLEs is now more
clearly visible: TLEs are detected largely over the sea in winter followed by a
negligible activity over land in spring. In summer, TLEs are detected largely
over land and the higher alpine regions (see both activity over the Pyrenees
and the Alps) abruptly shifting the activity over the coastal areas and sea in
autumn. The strongest activity over the sea in most areas is seen in autumn,
together with a strong in activity throughout Italy (with maximum activity
spread in autumn and winter). Particularly interesting is the highlight of the
activity over the Southern Mediterranean Sea around Cyprus, shifting from
the southern area in autumn to the northern area in winter.

4.3 TLE seasonal cycle above Europe: a monthly mean perspective

We further analyzed the seasonal evolution of TLEs calculating total average
TLE and lightning stroke counts above Europe per month (see Fig. 8). The
calculation for TLEs was performed without the inclusion of ILAN observa-
tions of the Southern Mediterranean in order to have a more uniform area
and a more consistent seasonal behavior; for consistency, we included only lightning data within the main TLE coverage (see white contours in Fig. 3). Monthly nighttime lightning and TLE counts (Fig. 8) both show a largely consistent seasonal cycle with generally two active seasons (summer and autumn) separated by a deep minimum in late winter-spring (mainly March and April), and a relative minimum in late summer-early autumn (in September for TLEs, and August for lightning but not when including 2013). Maxima are reached in August and November for TLEs, and July and November for lightning. This behavior is fairly consistent among the 5 years, apart from a very active June in 2009 (which was largely affected by the most prolific thunderstorm in the sample, see 4.4). Both in summer and autumn, TLEs appear to peak some time later than lightning, but in both cases November is the most prolific month. In Northern Europe, most of the TLE observations are made in late summer (Aug-Oct) between 19 and 03 UTC. Despite most of the lightning occurs in June to August in this region, the relative number of TLE detections is lower in this period partly because of the bright summer sky but mostly due to the much higher threshold in lightning charge moment change required to initiate daytime sprites (Stanley et al., 2000).

In addition to total TLE counts, bottom panels in Fig. 8 report the individual seasonal evolution of the occurrence of sprites, elves and halos (left to right). The overall behavior of the TLE seasonal cycle is very close to that of sprites, although elves contribute significantly to the November peak, which in the sprite cycle is smaller than the August one. Elves partially contribute also to December and January, making the period November-January the only one with significant activity, consistently with the shift in activity over the sea shown in Fig. 7. Halos tend to mimic the sprite seasonal cycle and lead to no significant effect on the TLE seasonal cycle.

The seasonal evolution of the number of nighttime hours should be taken in consideration when comparing summer and autumn counts since optical observations of TLEs are only available at night. This effect can be very pronounced at high latitude. In order to compensate for the varying length of nighttime conditions, we calculated an average nighttime duration (considering duration for Bari, Budapest, Warsaw and Granada – see dashed grey curve in Fig. 8) and normalized the monthly counts to a constant 12-hour night. The results are reported in Fig. 8 together with monthly mean lightning activity from the whole day (right panels). The overall TLE behavior is confirmed, with two distinct summer and autumn seasons now reaching similar total counts per month. An almost flat distribution between May to November is found for lightning, with an increase during July in some years.

An analysis of the time of occurrence of the observed TLEs is shown in Fig. 9 (top-left panel), reporting for each month both the average start and end of the observations (step lines) and the earliest as well as the latest observation (circles). The average start and end tends to be very close to the reference local midnight (at 11 UTC, dotted black line) in May to August, then increasing the average observation interval with increasing nighttime hours in autumn and winter, with a tendency to last longer during the second part of the night.
First observations of each month generally start 1 or 2 hours after sunset (see bottom of the figure), whereas the last observations typically occur within 1 hour before sunrise, and sunrise (or slightly beyond). This indicates well that nighttime ionospheric conditions change locally only about 20 minutes past the ground-level sunrise time (Satori et al., 2007). Note that the time information for February to April is based on a very limited number of observations.

Figure 9 also reports the seasonal evolution of the average latitude of the observations during the 5 years (top-right panel, compare to Fig. 7). The seasonal changes in the distribution discussed above are now summarized by an average latitude: It cycles between 41–42 ° latitude in autumn and winter, and 46–47 ° latitude in summer, with a fairly consistent behavior among the years, and therefore well characterizing the geographical change above Europe. Since the main seasonal change in TLE activity occurs in latitude, no meaningful information can be extracted from average longitude.

## 4.4 TLE-producing thunderstorms

The seasonal evolution of the number of observed TLE-thunderstorms in 2009 to 2013 is reported in Fig. 9 (bottom-left panel) in a similar fashion as for the number of TLEs in Fig. 8. For consistency, data from Israel were again excluded from the analysis. Here we assume that each entry of the database (i.e. a collection of TLEs observed above the same thunderstorm or closely related thunderstorm cells) can be considered as an individual thunderstorm system, although some database entries may likely extend over several thunderstorm cells which were observed during the observing period. In general, separated sprite-producing periods were reported so the opposite may also be true: Several entries may correspond to the same thunderstorm at different stages. Care should thus be taken to interpret these results correctly. Peaks in number of observed TLE-thunderstorms are reached in July, August and November, consistently with lightning activity, with about 150 TLE-thunderstorms accumulating over the 5 years. Comparison to TLE counts in Fig. 8 (top-left panel) shows that November is characterized by fewer but more prolific TLE-thunderstorms as compared to summer months. Out of 1018 observed thunderstorms, 801 (79%) were reported with less than 10 TLEs each, 921 (90%) with less than 20, 97 (10%) with 20 or more, 56 (6%) with more than 30, 34 (3%) with more than 40. The majority of the observations were related to relatively small thunderstorms and small TLE area. An inspection of the size of the reported TLE area showed in fact about 600 entries reported with a 0-10,000 or 10-20,000 km² area, followed by a second group of about 200 entries with a 30-40,000 km² area and the remaining 200 entries with reported larger areas. Considering that also the possible motion of the thunderstorm during the observation affect the area, this is an indication that the TLE production area was often close to the size of TLEs and of the uncertainty in their geolocation. Clearly, closer analysis of individual thunderstorms will be needed to further asses the actual TLE-production area.
The distribution in number of TLEs observed per thunderstorm is reported in Fig. 9 (bottom-right panel) in a log-log view. The occurrence is defined as the fraction of thunderstorms having a certain number of TLEs; only cases occurring more than once are considered. The distribution of number of thunderstorms as function of number of TLEs per thunderstorm (n) follows a power law proportional to $n^{-k}$ as shown by a best linear fit in the log-log plot with coefficient $k = -1.5$ (dashed black line) and having a correlation of -0.96. The power law fits for individual years 2009 to 2013 are also shown, and have correlation coefficients -0.87, -0.97, -0.91, -0.88 and -0.95 respectively. The linear fits for individual years are similar among them (-0.94 for 2009 and between -1 and -1.2 for the other years), but have a weaker slope as compared to the multiyear average: this is due to the fact that increasing the number of entries, more bins with more than one count enter the curve and these bins are increasing the weight of the fit on the right tail of the distribution. Correlation coefficients and slopes for individual years and for the multiyear distribution show a consistent power law behavior, suggesting the law may be expected to apply for large sample of thunderstorms observed over Europe in future years.

The 30 most prolific thunderstorms are summarized in Tab. 1. In particular, 12 of these thunderstorms were observed during the summer season (or late spring), and 18 during the autumn-winter season. Seven thunderstorms were reported with more than 100 TLEs each, the most prolific one being observed 29 October 2013 with 195 TLEs (193 sprites and 2 halos), followed by a thunderstorm on 10 June 2009 with 147 TLEs (146 sprites and 1 halo, a thunderstorm analyzed in detail by Bór (2013)) and by a thunderstorm on 28 November 2011 with 140 TLEs (131 sprites, 8 elves and 1 halo). The 29 October 2013 case was studied by (Soula et al., 2017), with more than 100 sprite videos reported in terms sprite/stroke pairs. The 10 June 2009 thunderstorm makes up most of the anomalously high June TLE count found in 2009 as shown in Fig. 8 (top-left panel). Most of the prolific thunderstorms were reported with almost only sprites. Seven thunderstorms were reported with a significant (greater than 30) number of elves and were all observed during the autumn-winter season (4 in November, one in December and one in January). The 12 December 2009 thunderstorms included observation of the gigantic jet discussed above, and of an upward lightning. For comparison, data for the 12 November 2011 thunderstorm observed from Israel were included in the table: this is the only case of a very prolific thunderstorm observed to last only 40 minutes, whereas all other prolific thunderstorms observed above Europe and coastal areas lasted several hours each.

# 5 Discussion

The dataset and climatology presented in the previous section represent the first attempt to have such a coordinated ground-based climatology over Europe and Mediterranean sea. Similar continuous observations were performed over limited regions such as the US High Plains (Lyons, 1996) or in Japan
(Adachi et al., 2005; Suzuki et al., 2011). Considering that only a negligible fraction of satellite observations of TLEs are currently taken over Europe (Chen et al., 2008), this is also the largest European TLE dataset available to date. The overall distribution and seasonal cycle of TLEs present robust features repeated in all the 5 years included in the sample, both as number of observed TLE-producing thunderstorms and as number of observed TLEs. This includes the two-peak seasonal cycle, with maxima in summer and late autumn, separated by two minima in March-April and less pronouncedly in September-October. The two active seasons are due to land driven convection during summer, typically associated with large thunderstorms, and sea driven convection during autumn and winter, found in synoptic-scale weather systems which transport cold continental air masses over the relatively warm Mediterranean Sea water, leading to instability and convection, and generating smaller thunderstorm cells. For example, the thunderstorm producing 54 TLEs on 12 December 2009 including the first European GJ, was composed by several small cells with cloud top height of only 6 km (van der Velde et al., 2010a). Nevertheless, TLEs were observed to be produced by fewer but more prolific thunderstorms in autumn rather than in summer. The shift of TLE activity from the continental areas in summer to sea and coastal areas in autumn and winter is therefore abrupt both in the region of activity and in the leading processes driving the TLE production.

The observed evolution of the TLE distribution is remarkably consistent with lightning activity reported by the WWLLN network within the region of Eurosphere observations, similarly to the migration from land to coastal and maritime regions previously observed over the oceans (Füllekrug et al., 2002). The adoption of an observational sensitivity map based on actual regions covered by each camera and the related simulated TLE climatology allows a better interpretation of the main features of the observed climatologies, assigning them to either observational limitations or lightning activity. Considering the number of shortages that were not included in the calculation, the agreement between the observed and simulated climatologies is remarkable. On the one hand, dedicated studies can be planned to further inspect specific regions and try to resolve the differences in terms of observational capabilities. On the other hand, this comparison also sheds light on possible differences in the behavior of TLE-producing thunderstorms in terms of production rates: can one expect the same TLE to lightning rate for all thunderstorms? Our results point to a very good consistency in the majority of cases but cannot exclude large differences. Discrepancies will therefore need to be further investigated inspecting the +CGs/-CGs ratio or further characteristics of the thunderstorms. Thunderstorm characteristics can be very much different in autumn/winter maritime thunderstorms and in summer continental ones (e.g., McGorman and Rust, 1998). Interestingly, our results for the summer season show a very similar distribution as those found with radio measurements by (Füllekrug and Reising, 1998, their figure 3).

This includes also a larger fraction of high charge moment change and high peak current lightning, respectively needed for sprite and elve generation (e.g.
Pasko et al., 2011). The climatological approach will aid in the comparison to relevant electrical and atmospheric parameters at large scale.

In agreement with previous studies of European TLE-producing thunderstorms (Neubert et al., 2001; Neubert et al., 2005, 2008; Soula et al., 2010; van der Velde et al., 2010b), the majority of TLE-producing thunderstorms in 2009-2013 was relatively small, with a limited TLE area and were reported to produce a small number of TLEs as compared to large thunderstorm observed e.g. in North or South America (Lyons, 1996; Taylor et al., 2008). The power law we found describing the number of thunderstorms producing a certain number of TLEs is novel and not reported before and would certainly need to be proved also in other continents. This power law behavior is not unexpected, since it is often found in several other natural phenomena, such as when describing the occurrence and intensity of e.g., tornadoes, fires, earthquakes (see the review by Pinto et al., 2012). Furthermore, it is also found in describing the distribution of terrestrial gamma ray flashes (TGFs), the second type of exotic emissions produced by thunderstorm activity (Smith et al., 2005). Since the power law was found to be consistent over the 4 most prolific years, it suggests a possible new pathway for modeling the occurrence and distribution of TLEs and TLE-producing thunderstorms above Europe and other regions, so as to further fill the gaps in the observations and allow inclusion in global models (see e.g. the parametrization adopted by Arnone et al., 2014). Given a large sample of TLEs (e.g. from intermittent satellites passages over Europe), the scaling to obtain yearly rates should be performed by clustering the observations over individual thunderstorms as described by the power law. It also allows a new approach for comparing the distribution of TLE-producing thunderstorms over different regions of the globe.

A selection mechanism due to better visibility in winter cannot be excluded to have contributed to observing more elves, halos, upward lightning and blue jets during autumn/winter, besides the expected higher rate of occurrence of elves over the sea. Comparison to Chen et al. (2008) and Newsome and Inan (2010) shows that a strong selection mechanism exists in ground optical observations in favor of sprites: the elve to sprite detection ratio in our climatology is in fact only 1:11 as compared to roughly 6:1 (30:1 accounting for detection efficiency corrections) found by the ISUAL satellite (Chen et al., 2008), and 6:1 found by ground using high time resolution photometer array (Newsome and Inan, 2010). However, the detection of elves almost exclusively over maritime thunderstorms is well in agreement with what found by previous studies (Chen et al., 2008): if limiting to November or January the ratio of elves to sprites ratio drops to 1:3. Besides shortages in the optical sensitivity to elves, the algorithms used for identifying elve images in the video frames are also less effective than for sprites. A similar detection bias occurs for halos, which we observed with a 1:26 halo to sprite ratio, whereas in the ISUAL satellite observations it was reported around 4:5. The origin and consequences of different detection efficiency was discussed also by Williams et al. (2012).

As discussed in Sect. 3.1, a key shortage of the adopted dataset is the inhomogeneity of the observational coverage. Cameras have their own detection
efficiency, they are often moved to different locations or pointing directions, or their views can be affected by the cloud coverage along their lines of sight. Besides, some cameras are used to track the evolution of thunderstorms (i.e. their pointing direction follows peak activity), other cameras constantly cover the same region (i.e. they have a fixed pointing direction). The detection efficiency of our network of cameras therefore vary greatly across the observed regions, and this can be even greater when talking about elves and halos which are at the edge of the camera sensitivity. Given the information recorded, an analysis of the detection efficiency of each system was not possible (see e.g. the approach developed by Ignaccolo et al., 2008). The overlap of different cameras can certainly overcome part of the shortages of the individual systems and the adopted sensitivity function was shown to lead to a broad consistency between observed and simulated TLE distributions. Even though one cannot expect a 100% detection efficiency, the observations of the same TLE-producing thunderstorms recorded by multiple cameras of the network suggest a very high efficiency in regions where the coverage function is maximized. The consistency of the observed TLE seasonal changes with the lightning activity shows that the climatological approach is robust if read in terms of TLE-producing thunderstorms. This further underlines the robustness of a joint TLE climatology, considering together TLE-producing thunderstorms, which are often observed to produce several types of TLEs, several of which have likely been missed. This novel climatology therefore represents an ideal reference for several TLE space missions (e.g. GLIMS, ASIM and TARANIS) operating over the next few years. The joint use of ground and space-based observations will improve our simulation and the possibility of extending the climatology outside its current limitations.

## 6 Conclusions

We presented a first climatology of TLEs over Europe and the Mediterranean Sea based on a coordinated database of optical observations by the Eurosphere ground based network from 2009 to 2013. The main features of the TLE seasonal cycle were found to be robust, repeating over the years and consistent with lightning activity. The main TLE activity shifts from continental areas in summer to coastal and sea areas in late autumn and early winter. The largest number of TLEs per month is recorded in November, aided by the longer nighttime duration. In March and April TLE activity is almost completely halted, which, considering the availability of continuously operated cameras and agreement with lightning activity, cannot be due to a bias in active observational campaigns. Elves are observed almost exclusively over autumn-winter maritime thunderstorms, whereas sprites and halos follow the seasonal changes from land to sea.

Taking into consideration both observational coverage and WWLLN lightning activity in a joint simulation, the consistency of the observed TLEs with lightning activity confirms that nighttime lightning is a good proxy for TLEs
at large scale. The distribution of nighttime WWLLN lightning can be used to
highlight regions where further studies are needed to optimize observational
shortages or investigate possible peculiarities in the TLE production. Our anal-
ysis points to a 1/1000 (4/1000) TLE to (nighttime) lightning ratio. These ra-
tios lead to a TLE global rate of 2.6 TLE/min (largely a sprite global rate in
our case), when they are applied on a global lightning rate of 44 flashes/s. This
value is largely consistent with previous estimates in literature (Ignaccolo et
al., 2008). Further consideration of detection efficiency for individual cameras
and WWLLN lightning data will need to be adopted to further refine such es-
timate since the efficiency of both were neglected in our calculation. The same
considerations apply to a future refinement of the overall climatology. We also
found that given a large sample of TLEs, a power law can well describe how
to distribute the expected number of TLE per thunderstorm. These features
should be taken into account when investigating the global impact of TLEs
onto the atmosphere and parametrizing TLEs into global models. They can
also be a valuable guidance for new ground-based observations and current or
upcoming space missions.

Despite the consistency on these general features, we found some geograph-
ical discrepancies and different rates in some months (e.g. May and September)
comparing TLEs to lightning. This confirms that disclosing the details of TLE
activity requires further investigation in relation to the type of lightning (e.g.
+CG/-CG ratio and intensity) and electrification mechanisms, since of course
different types of TLEs are well known to have different production mecha-
nisms. The proposed climatological approach may lead to further advances also
in this direction, by comparing large TLE samples to key atmospheric param-
eters. The very low rate of elve and halo observations is also evidence of the
strong limitation of ground-based optical cameras for their observation, stim-
ulating once more a needed synergetic approach with space missions such as
ASIM. At the same time, ongoing Eurosphere activities will attempt to remove
the main gaps in the observational coverage of Europe and the Mediterranean
sea, for example including observations from regional all-sky camera networks
which have developed over the last years.
References

Abarca, S. F., Corbosiero, K. L., and Galarneau, T. J.: An evaluation of the Worldwide Lightning Location Network (WWLLN) using the National Lightning Detection Network (NLDN) as ground truth, Journal of Geophysical Research (Atmospheres), 115, D18206, doi:10.1029/2009JD013411, 2010.

Adachi, T., Fukunishi, H., Takahashi, Y., Sato, M., Ohkubo, A., and Ya-mamoto, K.: Characteristics of thunderstorm systems producing winter sprites in Japan, Journal of Geophysical Research (Atmospheres), 110, D11203, doi:10.1029/2004JD005012, 2005.

Arnone, E. and Dinelli, B.: Remote Sensing Advances for Earth System Science, vol. XIII of Springer Earth System Sciences, chap. CHIMTEA - Chemical Impact of Thunderstorms on Earth's Atmosphere, pp. 1-14, Springer International Publishing, 2016.

Arnone, E. and Hauchecorne, A.: Stratosphere NO$_{y}$ Species Measured by MIPAS and GOMOS Onboard ENVISAT During 2002-2010: Influence of Plasma Processes onto the Observed Distribution and Variability, Space Science Reviews, 168, 315-332, doi:10.1007/s11214-011-9861-1, 2012.

Arnone, E., Berg, P., Arnold, N. F., Christiansen, B., and Thejll, P.: An estimate of the impact of transient luminous events on the atmospheric temperature, Advances in Geosciences, 2008.

Arnone, E., Berg, P., Boberg, F., Bór, J., Chanrion, O., Enell, C.-F., Ignaccolo, M., Mika, Á., Odzimek, A., van der Velde, O., Farges, T., Laursen, S., Neubert, T., and Sátori, G.: The Eurosphere 2005 campaign, in: Proceedings of the 33rd Annual European meeting on Atmospheric Studies by Optical Methods (33AM), edited by Arvelius, J., 2008.

Arnone, E., Kero, A., Dinelli, B. M., Enell, C.-F., Arnold, N. F., Papandrea, E., Rodger, C. J., Carlotti, M., Ridolfi, M., and Turunen, E.: Seeking sprite-induced signatures in remotely sensed middle atmosphere NO$_{2}$, Geophys. Res. Lett., 2008.

Arnone, E., Kero, A., Enell, C.-F., Carlotti, M., Rodger, C. J., Papandrea, E., Arnold, N. F., Dinelli, B. M., Ridolfi, M., and Turunen, E.: Seeking sprite-induced signatures in remotely sensed middle atmosphere NO$_{2}$: latitude and time variations, Plasma Sources Science Technology, 18, 034014+, doi:10.1088/0963-0252/18/3/034014, 2009.

Arnone, E., Smith, A. K., Enell, C.-F., Kero, A., and Dinelli, B. M.: WACCM climate chemistry sensitivity to sprite perturbations, Journal of Geophysical Research (Atmospheres), 119, 6958-6970, doi:10.1002/2013JD020825, 2014.

Arnone, E., Bor, J., Chanrion, O., Barta, V., Fullekrug, M., Labanti, R., Mezuman, K., Odzimek, A., Popek, M., Soula, S., Valeri, D., van der Velde, O., Yair, Y., Zanotti, F., Zoladek, P., and Neubert, T.: The Eurosphere 2009-2013 catalogue (Version 1.3) [DataSet], http://doi.org/10.5281/zenodo.3480108, 2019.

Barrington Leigh, C. P., Inan, U. S., and Stanley, M.: Identification of sprites and elves with intensified video and broadband array photometry, J. Geo-
phys. Res., 106, 1741-1750, 2001.

Barta, V., Haldoupis, C., Sátori, G., Buresova, D., Chum, J., Pozoga, M., Berényi, K. A., Bór, J., Popek, M., Kis, Á., and Bencze, P.: Searching for effects caused by thunderstorms in midlatitude sporadic E layers, Journal of Atmospheric and Solar-Terrestrial Physics, 161, 150-159, doi: 10.1016/j.jastp.2017.06.006, 2017.

Betz, H. D., Schmidt, K., Laroche, P., Blanchet, P., Oettinger, W. P., Defer, E., Dziewit, Z., and Konarski, J.: LINET - An international lightning detection network in Europe, Atmospheric Research, 91, 564-573, doi: 10.1016/j.atmosres.2008.06.012, 2009.

Blanc, E., Lefeuvre, F., Roussel-Dupré, R., and Sauvaud, J. A.: TARANIS: A microsatellite project dedicated to the study of impulsive transfers of energy between the Earth atmosphere, the ionosphere, and the magnetosphere, Advances in Space Research, 40, 1268-1275, doi:10.1016/j.asr.2007.06.037, 2007.

Boeck, W. L., Vaughan, Jr., O. H., Blakeslee, R., Vonnegut, B., and Brook, M.: Lightning induced brightening in the airglow layer, Geophysical Research Letters, 19, 99-102, doi:10.1029/91GL03168, 1992.

Boeck, W. L., Vaughan, O. H., Blakeslee, R. J., Vonnegut, B., Brook, M., and McKune, J.: Observations of lightning in the stratosphere, Journal of Geophys. Res., 100, 1465-1476, doi:10.1029/94JD02432, 1995.

Bór, J.: Optically perceptible characteristics of sprites observed in Central Europe in 2007-2009, Journal of Atmospheric and Solar-Terrestrial Physics, 92, 151-177, doi:10.1016/j.jastp.2012.10.008, 2013.

Bór, J., Sátori, G., and Betz, H. D.: Observation of TLEs in Central Europe from Hungary Supported by LINET, in: American Institute of Physics Conference Series, vol. 1118 of American Institute of Physics Conference Series, pp. 73-83, doi:10.1063/1.3137716, 2009.

Bór, J., Zelkó, Z., Hegedüs, T., Jäger, Z., Mlynarczyk, J., Popek, M., and Betz, H. D.: On the Series of +CG Lightning Strokes in Dancing Sprite Events, Journal of Geophysical Research (Atmospheres), 123, 11, doi: 10.1029/2017JD028251, 2018.

Chanrion, O. and Neubert, T.: Production of runaway electrons by negative streamer discharges, Journal of Geophysical Research (Space Physics), 115, A00E32, doi:10.1029/2009JA014774, 2010.

Chanrion, O., Crosby, N. B., Arnone, E., Boberg, F., van der Velde, O., Odzimek, A., Mika, Á., Enell, C.-F., Berg, P., Ignaccolo, M., Steiner, R. J., Laursen, S., and Neubert, T.: The EuroSprite2005 Observational Campaign: an example of training and outreach opportunities for CAL young scientists, Advances in Geosciences, 13, 3-9, 2007.

Chen, A. B., Kuo, C.-L., Lee, Y.-J., Su, H.-T., Hsu, R.-R., Chern, J.-L., Frey, H. U., Mende, S. B., Takahashi, Y., Fukunishi, H., Chang, Y.-S., Liu, T.-Y., and Lee, L.-C.: Global distributions and occurrence rates of transient luminous events, Journal of Geophysical Research (Space Physics), 113, A08306, doi:10.1029/2008JA013101, 2008.
Cummer, S. A., Li, J., Han, F., Lu, G., Jaugey, N., Lyons, W. A., and Nelson, T. E.: Quantification of the troposphere-to-ionosphere charge transfer in a gigantic jet, Nature Geoscience, 2, 617–620, doi:10.1038/ngeo607, 2009.

Ebert, U., Montijn, C., Briels, T. M. P., Hundsdorfer, W., Meulenbroek, B., Rocco, A., and van Veldhuizen, E. M.: The multiscale nature of streamers, Plasma Sources Science Technology, 15, 118–+, doi:10.1088/0963-0252/15/2/S14, 2006.

Ebert, U., Nijdam, S., Li, C., Luque, A., Briels, T., and van Veldhuizen, E.: Review of recent results on streamer discharges and discussion of their relevance for sprites and lightning, Journal of Geophysical Research (Space Physics), 115, A00E43, doi:10.1029/2009JA014867, 2010.

Enell, C.-F., Arnone, E., Chanrion, O., Adachi, T., Verronen, P. T., Seppälä, A., Neubert, T., Ulich, T., Turunen, E., Takahashi, Y., and Hsu, R.-R.: Parameterisation of the chemical effect of sprites in the middle atmosphere, Ann. Geophys., 26, 13–27, 2008.

Farges, T. and Blanc, E.: Characteristics of infrasound from lightning and sprites near thunderstorm areas, Journal of Geophysical Research (Space Physics), 115, A00E31, doi:10.1029/2009JA014700, 2010.

Farges, T. and Blanc, E.: Lightning and TLE electric fields and their impact on the ionosphere, Comptes Rendus Physique, 12, 171–179, doi: 10.1016/j.crhy.2011.01.013, 2011.

Farges, T., Blanc, E., Le Pichon, A., Neubert, T., and Allin, T. H.: Identification of infrasound produced by sprites during the Sprite2003 campaign, Geophys. Res. Lett., 32, 1813–+, 2005.

Farges, T., Blanc, E., and Tanguy, M.: Experimental evidence of D region heating by lightning-induced electromagnetic pulses on MF radio links, Journal of Geophysical Research (Space Physics), 112, A10302, doi: 10.1029/2007JA012285, 2007.

Franz, R. C., Nemzek, R. J., and Winckler, J. R.: Television Image of a Large Upward Electrical Discharge Above a Thunderstorm System, Science, 249, 48–51, 1990.

Fukunishi, H., Takahashi, Y., Kubota, M., and Sakanoi, K.: Elves: Lightning induced transient luminous events in the lower ionosphere, Geophys. Res. Lett., 23, 2157, 1996.

Füllekrug, M. and Reising, S. C.: Excitation of Earth-ionosphere cavity resonances by sprite-associated lightning flashes, Geophyiscal Research Letters, 25, 4145–4148, doi:10.1029/1998GL900133, 1998.

Füllekrug, M. and Reising, S. C.: Excitation of Earth-ionosphere cavity resonances by sprite-associated lightning flashes, Geophys. Res. Lett., 25, 4145–4148, 1998.

Füllekrug, M., Price, C., Yair, Y., and Williams, E. R.: Intense oceanic lightning, Annales Geophysicae, 20, 133–137, doi:10.5194/angeo-20-133-2002, 2002.

Füllekrug, M., Mareev, E. A., and Rycroft, M. J., eds.: Sprites, Elves and Intense Lightning Discharges, vol. 225 of NATO Science Series II. Mathematics, physics and chemistry, Springer Verlag, ISBN 1-4020-4628-6, 2006.
Füllekrug, M., Roussel-Dupré, R., Symbalisty, E. M. D., Chanrion, O., Odzimek, A., van der Velde, O., and Neubert, T.: Relativistic runaway breakdown in low-frequency radio, Journal of Geophysical Research (Space Physics), 115, A00E09, doi:10.1029/2009JA014468, 2010.

Füllekrug, M., Roussel-Dupré, R., Symbalisty, E. M. D., Colman, J. J., Chanrion, O., Soula, S., van der Velde, O., Odzimek, A., Bennett, A. J., Pasko, V. P., and Neubert, T.: Relativistic electron beams above thunderclouds, Atmospheric Chemistry & Physics, 11, 7747–7754, doi:10.5194/acp-11-7747-2011, 2011.

Füllekrug, M., Kolmasova, I., Santolik, O., Farges, T., Bór, J., Bennett, A., Parrot, M., Rison, W., Zanotti, F., Arnone, E., Mezentsev, A., Lan, R., Uhlir, L., Harrison, G., Soula, S., van der Velde, O., Pinçon, J.-L., Helling, C., and Diver, D.: Electron acceleration above thunderclouds, Environmental Research Letters, 8, 035027, doi:10.1088/1748-9326/8/3/035027, 2013a.

Füllekrug, M., Mezentsev, A., Soula, S., Velde, O., and Farges, T.: Sprites in low-frequency radio noise, Geophys. Res. Lett., 40, 2395–2399, doi: 10.1002/grl.50408, 2013b.

Ganot, M., Yair, Y., Price, C., Ziv, B., Sherez, Y., Greenberg, E., Devir, A., and Yaniv, R.: First detection of transient luminous events associated with winter thunderstorms in the eastern Mediterranean, Geophys. Res. Lett., 34, L12801, doi:10.1029/2007GL029258, 2007.

Gjesteland, T., Østgaard, N., Laviola, S., Miglietta, M. M., Arnone, E., Marisaldi, M., Fuschino, F., Collier, A. B., Fabró, F., and Montanya, J.: Observation of intrinsically bright terrestrial gamma ray flashes from the Mediterranean basin, Journal of Geophysical Research (Atmospheres), 120, 12, doi:10.1002/2015JD023704, 2015.

Gordillo-Vázquez, F. J.: Air plasma kinetics under the influence of sprites, Journal of Physics D: Applied Physics, 41, 234016+, doi:10.1088/0022-3727/41/23/234016, 2008.

Gordillo-Vázquez, F. J., Luque, A., and Simek, M.: Spectrum of sprite halos, Journal of Geophysical Research (Space Physics), 116, A09319, doi: 10.1029/2011JA016652, 2011.

Gordillo-Vázquez, F. J., Passas, M., Luque, A., Sánchez, J., van der Velde, O. A., and Montanyá, J.: High Spectral Resolution Spectroscopy of Sprites: A Natural Probe of the Mesosphere, Journal of Geophysical Research (Atmospheres), 123, 2336–2346, doi:10.1002/2017JD028126, 2018.

Greenberg, E. and Price, C.: A global lightning location algorithm based on the electromagnetic signature in the Schumann resonance band, Journal of Geophysical Research (Atmospheres), 109, D21111, doi:10.1029/2004JD004845, 2004.

Greenberg, E., Price, C., Yair, Y., Haldoupis, C., Chanrion, O., and Neubert, T.: ELF/VLF signatures of sprite-producing lightning discharges observed during the 2005 EuroSprite campaign, Journal of Atmospheric and Solar-Terrestrial Physics, 71, 1254–1266, doi:10.1016/j.jastp.2009.05.005, 2009.

Haldoupis, C., Neubert, T., Inan, U. S., Mika, A., Allin, T. H., and Marshall, R. A.: Subionospheric early VLF signal perturbations observed in one-to-one
association with sprites, Journal of Geophysical Research (Space Physics), 109, A10303, doi:10.1029/2004JA010651, 2004.

Haldoupis, C., Steiner, R. J., Mika, Á., Shalimov, S., Marshall, R. A., Inan, U. S., Bösinger, T., and Neubert, T.: “Early/slow” events: A new category of VLF perturbations observed in relation with sprites, Journal of Geophysical Research (Space Physics), 111, A11321, doi:10.1029/2006JA011960, 2006.

Haldoupis, C., Amvrosiadi, N., Cotts, B. R. T., van der Velde, O. A., Chanrion, O., and Neubert, T.: More evidence for a one-to-one correlation between Sprites and Early VLF perturbations, Journal of Geophysical Research (Space Physics), 115, A07304, doi:10.1029/2009JA015165, 2010.

Haldoupis, C., Cohen, M., Cotts, B., Arnone, E., and Inan, U.: Long-lasting D-region ionospheric modifications, caused by intense lightning in association with elve and sprite pairs, Geophys. Res. Lett., 39, L16801, doi: 10.1029/2012GL052765, 2012.

Haldoupis, C., Cohen, M., Arnone, E., Cotts, B., and Dietrich, S.: The VLF fingerprint of elves: Step-like and long-recovery early VLF perturbations caused by powerful CG lightning EM pulses, Journal of Geophysical Research (Space Physics), 118, 5392–5402, 2013.

Ignaccolo, M., Farges, T., Mika, A., Allin, T. H., Chanrion, O., Blanc, E., Fraser-Smith, A. C., and Füllekrug, M.: The Planetary Rate of Sprite Events, Geophys. Res. Lett., 33, L11808, doi:10.1029/2005GL025502, 2006.

Ignaccolo, M., Farges, T., Blanc, E., and Füllekrug, M.: Automated chirp detection with diffusion entropy: Application to infrasound from sprites, Chaos Solitons and Fractals, 38, 1039–1050, doi:10.1016/j.chaos.2007.02.011, 2008.

Inan, U. S., Barrington-Leigh, C., Hansen, S., Glukhov, V. S., Bell, T. F., and Rairden, R.: Rapid lateral expansion of optical luminosity in lightning-induced ionospheric flashes referred to as ‘elves’, Geophysical Research Letters, 24, 583–586, doi:10.1029/97GL00404, 1997.

Iwański, R., Odzimek, A., Clausen, L., Kanawade, V., Cnossen, I., and Edberg, N.: Meteorological study of the first observation of red sprites from Poland, Acta Geophysica, 57, 760–777, doi:10.2478/s11600-009-0008-7, 2009.

Krehbiel, P. R., Riouset, J. A., Pasko, V. P., Thomas, R. J., Rison, W., Stanley, M. A., and Edens, H. E.: Upward electrical discharges from thunderstorms, Nature, 2008.

Kulak, A. and Młynarczyk, J.: A new technique for reconstruction of the current moment waveform related to a gigantic jet from the magnetic field component recorded by an ELF station, Radio Science, 46, RS2016, doi: 10.1029/2010RS004475, 2011.

Lefeuvre, F., Blanc, E., Pinçon, J.-L., Roussel-Dupré, R., Lawrence, D., Sauvaud, J.-A., Rauch, J.-L., de Feraudy, H., and Lagoutte, D.: TARANIS - A Satellite Project Dedicated to the Physics of TLEs and TGFs, Space Sci. Rev., 137, 301–315, doi:10.1007/s11214-008-9414-4, 2008.

Liszka, L. and Hobara, Y.: Sprite-attributed infrasonic chirps - their detection, occurrence and properties between 1994 and 2004, Journal of Atmospheric and Solar-Terrestrial Physics, 68, 1179–1188, 2006.
Luque, A. and Ebert, U.: Emergence of sprite streamers from screening-ionization waves in the lower ionosphere, Nature Geoscience, 2, 757–760, doi:10.1038/ngeo662, 2009.

Lyons, A.: Characteristics of luminous structures in the stratosphere above thunderstorms as imaged by low-light video, Geophys. Res. Lett., 21, 875–878, 1994.

Lyons, W. A.: Sprite observations above the U.S. High Plains in relation to their parent thunderstorm systems, J. Geophys. Res., 101, 1996.

Mäkelä, A., Kantola, T., Yair, Y., and Raita, T.: Observations of TLEs above the Baltic sea on Oct 9 2009, Geophysica, 46, 79–90, 2010.

Mika, Á., Haldoupis, C., Marshall, R. A., Neubert, T., and Inan, U. S.: Subionospheric VLF signatures and their association with sprites observed during EuroSprite-2003, Journal of Atmospheric and Solar-Terrestrial Physics, 67, 1580–1597, doi:10.1016/j.jastp.2005.08.011, 2005.

Mika, Á., Haldoupis, C., Neubert, T., Su, H. T., Hsu, R. R., Steiner, R. J., and Marshall, R. A.: Early VLF perturbations observed in association with elves, Annales Geophysicae, 24, 2179–2189, 2006.

Mlynarczyk, J., Bór, J., Kulak, A., Popek, M., and Kubisz, J.: An unusual sequence of sprites followed by a secondary TLE: An analysis of ELF radio measurements and optical observations, Journal of Geophysical Research (Space Physics), 120, 2241–2254, doi:10.1002/2014JA020780, 2015.

Montanyà, J., van der Velde, O., Romero, D., March, V., Solà, G., Pineda, N., Arrayas, M., Trueba, J. L., Reglero, V., and Soula, S.: High-speed intensified video recordings of sprites and elves over the western Mediterranean Sea during winter thunderstorms, Journal of Geophysical Research (Atmospheres), 115, A00E18, doi:10.1029/2009JA014508, 2010.

NaitAmor, S., AlAbdoadaim, M. A., Cohen, M. B., Cotts, B. R. T., Soula, S., Chanrion, O., Neubert, T., and Abdelatif, T.: VLF observations of ionospheric disturbances in association with TLEs from the EuroSprite-2007 campaign, Journal of Geophysical Research (Space Physics), 115, A00E47, doi:10.1029/2009JA015026, 2010.

Neubert, T.: On Sprites and Their Exotic Kin, Science, 300, 747–749, 2003.

Neubert, T.: ASIM - an Instrument Suite for the International Space Station, in: American Institute of Physics Conference Series, vol. 1118 of *American Institute of Physics Conference Series*, pp. 8–12, doi:10.1063/1.3137718, 2009.

Neubert, T. and Chanrion, O.: On the electric breakdown field of the mesosphere and the influence of electron detachment, Geophys. Res. Lett., 40, 2373–2377, doi:10.1002/grl.50433, 2013.

Neubert, T., Allin, T. H., Stenbaek-Nielsen, H., and Blanc, E.: Sprites Over Europe, Geophys. Res. Lett., 28, 3585–3588, 2001.

Neubert, T., Allin, T. H., Blanc, E., Farges, T., Haldoupis, C., Mika, A., Soula, S., Knutsson, L., van der Velde, O., Marshall, R. A., Inan, U., Sátori, G., Bór, J., Hughes, A., Collier, A., Laursen, S., and Rasmussen, I.: Co-ordinated observations of transient luminous events during the EuroSprite2003 campaign, J. Atmos. and Solar-Terr. Phys., 67, 807–820, doi:
10.1016/j.jastp.2005.02.004, 2005.

Neubert, T., Rycroft, M., Farges, T., Blanc, E., Chanrion, O., Arnone, E., Odzimek, A., Arnold, N., Enell, C.-F., Turunen, E., Bösinger, T., Mika, Á., Haldoupis, C., Steiner, R. J., van der Velde, O., Soula, S., Berg, P., Boberg, F., Thejll, P., Christiansen, B., Ignaccolo, M., Füllekrug, M., Verronen, P. T., Montanya, J., and Crosby, N.: Recent Results from Studies of Electric Discharges in the Mesosphere, Surveys in Geophysics, 29, 71, doi:10.1007/s10712-008-9043-1, 2008.

Neubert, T., Chanrion, O., Arnone, E., Zanotti, F., Cummer, S., Li, J., Füllekrug, M., Soula, S., and van der Velde, O.: The properties of a gigantic jet reflected in a simultaneous sprite: Observations interpreted by a model, Journal of Geophysical Research (Space Physics), 116, A12329, doi:10.1029/2011JA016928, 2011.

Newsome, R. T. and Inan, U. S.: Free-running ground-based photometric array imaging of transient luminous events, Journal of Geophysical Research (Space Physics), 115, A00E41, doi:10.1029/2009JA014834, 2010.

Odzimek, A., Clausen, L. B. N., Kanawade, V., Cnossen, I., Edberg, N. J. T., Faedi, F., Del Moro, A., Ural, U., Byckling, K., Krzaczkowski, P., Iwański, R., Struzik, P., Pajek, M., and Gajda, W.: SPARTAN Sprite-Watch 2007 Campaign, in: Young Scientists 15th Proceedings, edited by Choliy, V. Y. and Ivashchenko, G., pp. 64–67, 2008.

Parra-Rojas, F. C., Luque, A., and Gordillo-Vázquez, F. J.: Chemical and electrical impact of lightning on the Earth mesosphere: The case of sprite halos, Journal of Geophysical Research (Space Physics), 118, 5190–5214, doi:10.1002/jgra.50449, 2013.

Parra-Rojas, F. C., Luque, A., and Gordillo-Vázquez, F. J.: Chemical and thermal impacts of sprite streamers in the Earth’s mesosphere, Journal of Geophysical Research (Space Physics), 120, 8899–8933, 2015.

Pasko, V. P.: Blue jets and gigantic jets: transient luminous events between thunderstorm tops and the lower ionosphere, Plasma Physics and Controlled Fusion, 50, 124050, doi:10.1088/0741-3335/50/12/124050, 2008.

Pasko, V. P.: Recent advances in theory of transient luminous events, Journal of Geophysical Research (Space Physics), 115, A00E35, 2010.

Pasko, V. P., Inan, U. S., and Bell, T. F.: Sprite as evidence of vertical gravity wave structures above mesoscale thunderstorms, Geophys. Res. Lett., 24, 1735–1738, 1997.

Pasko, V. P., Inan, U. S., and Bell, T. F.: Spatial structure of sprites, Geophys. Res. Lett., 25, 2123–2126, 1998.

Pasko, V. P., Stanley, M. A., Mathews, J. D., Inan, U. S., and Wood, T. G.: Electrical discharge from a thundercloud top to the lower ionosphere, Nature, 416, 152–154, 2002.

Pasko, V. P., Yair, Y., and Kuo, C.-L.: Lightning Related Transient Luminous Events at High Altitude in the Earth’s Atmosphere: Phenomenology, Mechanisms and Effects, Space Science Reviews, p. 287, doi:10.1007/s11214-011-9813-9, 2011.
Pérez-Invernón, F. J., Luque, A., and Gordillo-Vázquez, F. J.: Modeling the Chemical Impact and the Optical Emissions Produced by Lightning-Induced Electromagnetic Fields in the Upper Atmosphere: The case of Halos and Elves Triggered by Different Lightning Discharges, Journal of Geophysical Research (Atmospheres), 123, 7615–7641, doi:10.1029/2017JD028235, 2018a.

Pérez-Invernón, F. J., Luque, A., Gordillo-Vázquez, F. J., Sato, M., Ushio, T., Adachi, T., and Chen, A. B.: Spectroscopic Diagnostic of Halos and Elves Detected From Space-Based Photometers, Journal of Geophysical Research (Atmospheres), 123, 12,917–12,941, doi:10.1029/2018JD029053, 2018b.

Pérez-Invernón, F. J., Gordillo-Vázquez, F. J., Smith, A. K., Arnone, E., and Winkler, H.: Global Occurrence and Chemical Impact of Stratospheric Blue Jets Modeled With WACCM4, Journal of Geophysical Research (Atmospheres), 124, 2841–2864, doi:10.1029/2018JD029593, 2019.

Petrov, N. I. and Petrova, G. N.: Physical mechanisms for the development of lightning discharges between a thundercloud and the ionosphere, Journal of Technical Physics, 44, 472–475, doi:10.1134/1.1259327, 1999.

Pinto, C. M. A., Mendes Lopes, A., and Machado, J. A. T.: A review of power laws in real life phenomena, Communications in Nonlinear Science and Numerical Simulations, 17, 3558–3578, doi:10.1016/j.cnsns.2012.01.013, 2012.

Price, C., Asfur, M., Lyons, W., and Nelson, T.: An improved ELF/VLF method for globally geolocating sprite-producing lightning, Geophys. Res. Let., 29, 1031, doi:10.1029/2001GL013519, 2002.

Rodger, C. J.: Red sprites, upward lightning, and VLF perturbations, Rev. Geophys., 37, 317–336, 1999.

Rodger, C. J., Werner, S., Brundell, J., Thomson, N.R. and Lay, E., Holzworth, R., and Dowden, R.: Detection efficiency of the VLF World-Wide Lightning Location Network (WWLLN): Initial case study, Annales. Geophys., 24, 3197–3214, 2006.

Rudlosky, S. D. and Shea, D. T.: Evaluating WWLLN performance relative to TRMM/LIS, Geophysical Research Letters., 40, 2344–2348, doi: 10.1002/grl.50428, 2013.

Rycroft, M. J. and Harrison, R. G.: Electromagnetic Atmosphere-Plasma Coupling: The Global Atmospheric Electric Circuit, Space Sci. Rev., p. 305, doi:10.1007/s11214-011-9830-8, 2011.

Rycroft, M. J. and Odzimek, A.: Effects of lightning and sprites on the ionospheric potential, and threshold effects on sprite initiation, obtained using an analog model of the global atmospheric electric circuit, Journal of Geophysical Research (Space Physics), 115, A00E37, doi:10.1029/2009JA014758, 2010.

Rycroft, M. J., Odzimek, A., Arnold, N. F., Füllekrug, M., Kułak, A., and Neubert, T.: New model simulations of the global atmospheric electric circuit driven by thunderstorms and electrified shower clouds: The roles of lightning and sprites, Journal of Atmospheric and Solar-Terrestrial Physics, 69, 2485–2509, doi:10.1016/j.jastp.2007.09.004, 2007.
Sátori, G., Rycroft, M., Bencze, P., Märcz, F., Bór, J., Barta, V., Nagy, T., and Kovács, K.: An Overview of Thunderstorm-Related Research on the Atmospheric Electric Field, Schumann Resonances, Sprites, and the Ionosphere at Sopron, Hungary, Surveys in Geophysics, 34, 255–292, doi:10.1007/s10712-013-9222-6, 2013.

Savtchenko, A., Mitzeva, R., Tsenova, B., and Kolev, S.: Analysis of lightning activity in two thunderstorm systems producing sprites in France, Journal of Atmospheric and Solar-Terrestrial Physics, 71, 1277–1286, doi: 10.1016/j.jastp.2009.04.010, 2009.

Sentman, D. D. and Wescott, E. M.: Observations of Upper Atmospheric Optical Flashes Recorded from an Aircraft, Geophys. Res. Lett., 20, 2857–2860, 1993.

Sentman, D. D., Wescott, E. M., Osborne, D. L., Hampton, D. L., and Heavner, M. J.: Preliminary results from the Sprites94 aircraft campaign: 1. Red sprites, Geophys. Res. Lett., 22, 1205–1208, 1995.

Smith, D. M., Lopez, L. I., Lin, R. P., and Barrington-Leigh, C. P.: Terrestrial Gamma-Ray Flashes Observed up to 20 MeV, Science, 307, 1085–1088, doi: 10.1126/science.1107466, 2005.

Soula, S., van der Velde, O., Palmiéri, J., Chanrion, O., Neubert, T., Montanyà, J., Gangneron, F., Meyerfeld, Y., Lefeuvre, F., and Lointier, G.: Characteristics and conditions of production of transient luminous events observed over a maritime storm, Journal of Geophysical Research (Atmospheres), 115, D16118, doi:10.1029/2009JD012066, 2010.

Soula, S., Iacovella, F., van der Velde, O., Montanyà, J., Füllekrug, M., Farges, T., Bór, J., Georgis, J.-F., NaitAmor, S., and Martin, J.-M.: Multi-instrumental analysis of large sprite events and their producing storm in southern France, Atmospheric Research, 135, 415–431, 2014.

Soula, S., Defer, E., Füllekrug, M., van der Velde, O., Montanya, J., Bousquet, O., Mlynarczyk, J., Coquillat, S., Pinty, J.-P., Rison, W., Krehbiel, P. R., Thomas, R., and Pedeboy, S.: Time and space correlation between sprites and their parent lightning flashes for a thunderstorm observed during the HyMeX campaign, Journal of Geophysical Research (Atmospheres), 120, 11, doi:10.1002/2015JD023894, 2015.

Soula, S., Mlynarczyk, J., Füllekrug, M., Pineda, N., Georgis, J.-F., van der Velde, O., Montanyà, J., and Fabró, F.: Dancing sprites: Detailed analysis of two case studies, Journal of Geophysical Research (Atmospheres), 122, 3173–3192, doi:10.1002/2016JD025548, 2017.

Stenbaek-Nielsen, H. C., Moudry, D. R., Wescott, E. M., Sentman, D. D., and Sabbas, F. T. S.: Sprites and possible mesospheric effects, Geophysical Research Letters, 27, 3829–32, 2000.

Stenbaek-Nielsen, H. C., Haaland, R., McHarg, M. G., Hensley, B. A., and Kanmae, T.: Sprite initiation altitude measured by triangulation, Journal of Geophysical Research (Space Physics), 115, A00E12, 2010.

Su, H. T., Hsu, R. R., Chen, A. B., Wang, Y. C., Hsiao, W. S., Lai, W. C., Lee, L. C., Sato, M., and Fukunishi, H.: Gigantic jets between a thundercloud and the ionosphere, Nature, 423, 974–976, 2003.
Suzuki, T., Hayakawa, M., and Michimoto, K.: Small Winter Thunderstorm with Sprites and Strong Positive Discharge, IEEJ Transactions on Fundamentals and Materials, 131, 723–728, doi:10.1541/ieejfms.131.723, 2011.

Taylor, M. J., Bailey, M. A., Pautet, P. D., Cummer, S. A., Jaugey, N., Thomas, J. N., Solorzano, N. N., Sao Sabbas, F., Holzworth, R. H., Pinto, O., and Schuch, N. J.: Rare measurements of a sprite with halo event driven by a negative lightning discharge over Argentina, Geophys. Res. Lett., 35, L14812, doi:10.1029/2008GL033984, 2008.

Vadislavsky, E., Yair, Y., Erlick, C., Price, C., Greenberg, E., Yaniv, R., Ziv, B., Reicher, N., and Devir, A.: Indication for circular organization of column sprite elements associated with Eastern Mediterranean winter thunderstorms, Journal of Atmospheric and Solar-Terrestrial Physics, 71, 1835–1839, doi:10.1016/j.jastp.2009.07.001, 2009.

van der Velde, O. A. and Montanyà, J.: Statistics and variability of the altitude of elves, Geophys. Res. Lett., 43, 5467–5474, doi:10.1002/2016GL068719, 2016.

van der Velde, O. A., Mika, Á., Soula, S., Haldoupis, C., Neubert, T., and Inan, U. S.: Observations of the relationship between sprite morphology and in-cloud lightning processes, Journal of Geophysical Research (Atmospheres), 111, D15203, doi:10.1029/2005JD006879, 2006.

van der Velde, O. A., Bór, J., Li, J., Cummer, S. A., Arnone, E., Zanotti, F., Füllekrug, M., Haldoupis, C., NaitAmor, S., and Farges, T.: Multi-instrumental observations of a positive gigantic jet produced by a winter thunderstorm in Europe, Journal of Geophysical Research (Atmospheres), 115, D24301, doi:10.1029/2010JD014442, 2010a.

van der Velde, O. A., Montanyà, J., Soula, S., Pineda, N., and Bech, J.: Spatial and temporal evolution of horizontally extensive lightning discharges associated with sprite-producing positive cloud-to-ground flashes in northeastern Spain, Journal of Geophysical Research (Space Physics), 115, A00E56, doi: 10.1029/2009JA014773, 2010b.

van der Velde, O. A., Montanyà, J., Soula, S., Pineda, N., and Mlynar-czyk, J.: Bidirectional leader development in sprite-producing positive cloud-to-ground flashes: Origins and characteristics of positive and negative leaders, Journal of Geophysical Research (Atmospheres), 119, 12, doi: 10.1002/2013JD021291, 2014.

Wescott, E. M., Sentman, D., Osborne, D., Hampton, D., and Heavner, M.: Preliminary results from the Sprites94 aircraft campaign. 2. Blue jets, Geophys. Res. Lett., 22, 1209–12, 1995.

Wescott, E. M., Stenbaek Nielsen, H. C., Sentman, D. D., Heavner, M. J., Moudry, D., and São Sabbas, F. T.: Triangulation of sprites, associated halos and their possible relation to causative lightning and micrometeors, Journal of Geophysical Research, 106, 10467–77, 2001.

Whitley, T., Füllekrug, M., Rycroft, M., Bennett, A., Wyatt, F., Elliott, D., Heinson, G., Hitchman, A., Lewis, A., Sefako, R., Fourie, P., Dyers, J., Thomson, A., and Flower, S.: Worldwide extremely low frequency magnetic field sensor network for sprite studies, Radio Science, 46, RS4007,
doi:10.1029/2010RS004523, 2011.

Williams, E., Kuo, C.-L., Bór, J., Sátori, G., Newsome, R., Adachi, T., Boldi, R., Chen, A., Downes, E., Hsu, R. R., Lyons, W., Saba, M. M. F., Taylor, M., and Su, H. T.: Resolution of the sprite polarity paradox: The role of halos, Radio Science, 47, RS2002, doi:10.1029/2011RS004794, 2012.

Winkler, H. and Notholt, J.: The chemistry of daytime sprite streamers - a model study, Atmospheric Chemistry & Physics, 14, 3545-3556, doi: 10.5194/acp-14-3545-2014, 2014.

Winkler, H. and Notholt, J.: A model study of the plasma chemistry of stratospheric Blue Jets, Journal of Atmospheric and Solar-Terrestrial Physics, 122, 75-85, doi:10.1016/j.jastp.2014.10.015, 2015.

Yair, Y., Price, C., Ganot, M., Greenberg, E., Yaniv, R., Ziv, B., Sherez, Y., Devir, A., Bór, J., and Sátori, G.: Optical observations of transient luminous events associated with winter thunderstorms near the coast of Israel, Atmospheric Research, 91, 529-537, 2009.

Yair, Y., Price, C., Katzenelson, D., Rosenthal, N., Rubanenko, L., Ben-Ami, Y., and Arnone, E.: Sprite climatology in the Eastern Mediterranean Region, Atmospheric Research, 157, 108-118, doi:10.1016/j.atmosres.2014.12.018, 2015.
<figure><img src="image_2.png" /><figcaption>**Fig. 1** Location (white crosses) and estimated observational coverage (color scale) of the main optical camera systems involved in the Eurosphere network. The coverage is reported in terms of number of cameras observing a specific location weighted by the continuity or shortages of the contribution of each camera. See text for details.</figcaption></figure>
<figure><img src="image_3.png" /><figcaption>**Fig. 2** Top: climatology of observed TLEs (TLEs $\cdot 10^{-3} \text{km}^{-2} \text{yr}^{-1}$) for 2009 to 2013. Bottom: climatology of observed sprites, elves and halos (left to right). All panels adopt the same color scale.</figcaption></figure>
<figure><img src="image_4.png" /><figcaption>**Fig. 3** Climatology of lightning strokes detected by WWLLN (strokes km$^{-2}$yr$^{-1}$) for 2009 to 2013 considering the whole day (left) and nighttime only (right) counts. The white thick contours delimit the approximate area of TLE observations used in monthly mean calculations of WWLLN data.</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>**Fig. 4** Left – Ratio of observed TLEs to $10^3$ WWLLN nighttime strokes for 2009 to 2013. Right – expected TLEs based on WWLLN nighttime strokes (see Fig. 3), observational coverage (see Fig. 1) and an average 4/1000 TLE to nighttime stroke ratio (see left panel).</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>**Fig. 5** Climatology of observed TLEs (TLEs $\cdot 10^{-3} \text{km}^{-2} \text{yr}^{-1}$) for individual years 2009 to 2013 (top to bottom) and season winter (DJF: December, January, February), spring (MAM), summer (JJA) and autumn (SON) (left to right). The area of data coverage of each individual year is shaded in dark grey. Color grading refers to color bar in Fig. 2.</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>**Fig. 6** Climatology of nighttime lightning strokes detected by WWLLN (strokes km$^{-2}$yr$^{-1}$) for individual years 2009 to 2013 (top to bottom) and season winter (DJF: December, January, February), spring (MAM), summer (JJA) and autumn (SON) (left to right). White thick contours of the approximate TLE coverage (Fig. 3) are shown to ease comparison to TLE maps. Color grading refers to color bar in Fig. 3.</figcaption></figure>
<figure><img src="image_8.png" /><figcaption>**Fig. 7** Average seasonal variation of TLEs (see labels). Data are reported as seasonal fractional component of the 2009-2013 average of Fig. 2. The area of data coverage of the complete data set is shaded in dark grey.</figcaption></figure>
<figure><img src="image_9.png" /><figcaption>**Fig. 8** Cumulative number of observed WWLLN strokes (top), TLEs (center) and individual TLE types (bottom) per month including data from the start of the sample (January 2009) to the end of 2009, 2010, 2011, 2012 and 2013 (respectively in grey, blue, green, yellow and red). Lightning is shown in terms of nighttime strokes only (left) and whole day strokes (right). Total TLEs are shown as observed (left) and normalized to a constant 12-hour nighttime (right). The adopted average seasonal cycle of nighttime hours is shown in dashed grey. WWLLN data are scaled by a factor $10^4$ and averages calculated over the rectangular shapes shown in Fig. 3. Note the change in scale for individual TLEs in bottom panels.</figcaption></figure>
<figure><img src="image_10.png" /><figcaption>**Fig. 9** Top panels show the seasonal evolution of the average time (left) and latitude (right) for 2009, 2010, 2011, 2012 and 2013 (respectively in grey, blue, green, yellow and red). Time is shown as average start and average end of observations (step lines) and first and last observation (circles) of each month. Dashed lines show the seasonal evolution of the average sunrise and sunset (dark grey), and of an indicative 1 hour twilight (light grey). Bottom panels show the cumulative number of TLE-producing thunderstorms per month (left, color grading corresponds to Fig. 8) and the distribution of the number of TLEs per thunderstorms (right). A power law fit with correlation coefficient -0.96 is shown (dashed black line), together with power law fits for individual years 2009 to 2013 (colors as in top panels).</figcaption></figure>
<table><caption>**Table 1** Summary of the 30 most prolific TLE-thunderstorms observed in 2009 to 2013.</caption><thead><tr><th>Date</th><th>Time</th><th>Latitude</th><th>Longitude</th><th>TLE</th><th>Sprite</th><th>Elve</th><th>Halo</th></tr></thead><tbody><tr><td>29.10.2013</td><td>17:48 / 5:32</td><td>39.0 / 40.5</td><td>3.0 / 7.5</td><td>195</td><td>193</td><td>0</td><td>2</td></tr><tr><td>10.06-2009</td><td>21:00 / 2:00</td><td>45.0 / 47.0</td><td>19.0 / 23.0</td><td>147</td><td>146</td><td>0</td><td>1</td></tr><tr><td>28.11-2011</td><td>18:26 / 2:48</td><td>37.0 / 40.0</td><td>3.0 / 7.0</td><td>140</td><td>131</td><td>8</td><td>1</td></tr><tr><td>30.08-2012</td><td>21:23 / 4:05</td><td>37.0 / 40.0</td><td>0.0 / 2.5</td><td>129</td><td>123</td><td>0</td><td>6</td></tr><tr><td>26.07-2013</td><td>20:28 / 2:23</td><td>44.0 / 44.5</td><td>-1.0 / 1.5</td><td>111</td><td>111</td><td>0</td><td>7</td></tr><tr><td>12.10-2012</td><td>20:15 / 4:40</td><td>38.0 / 41.0</td><td>1.0 / 4.5</td><td>107</td><td>98</td><td>4</td><td>5</td></tr><tr><td>06.08-2013</td><td>19:37 / 0:49</td><td>48.0 / 50.0</td><td>12.5 / 14.5</td><td>101</td><td>84</td><td>0</td><td>16</td></tr><tr><td>28.11-2012</td><td>16:39 / 4:02</td><td>40.0 / 44.0</td><td>15.0 / 17.0</td><td>79</td><td>79</td><td>0</td><td>0</td></tr><tr><td>30.11.2009</td><td>23:14 / 6:07</td><td>38.0 / 40.0</td><td>4.0 / 6.0</td><td>79</td><td>14</td><td>65</td><td>0</td></tr><tr><td>29.01.2012</td><td>21:18 / 6:8</td><td>36.5 / 39.0</td><td>0.0 / 5.0</td><td>78</td><td>18</td><td>59</td><td>4</td></tr><tr><td>08.10.2009</td><td>18:17 / 0:47</td><td>43.0 / 44.0</td><td>2.0 / 5.0</td><td>77</td><td>77</td><td>0</td><td>2</td></tr><tr><td>12.11.2011</td><td>22:19 / 23:01</td><td>33.0 / 35.0</td><td>31.0 / 33.0</td><td>75</td><td>75</td><td>0</td><td>0</td></tr><tr><td>27.05.2009</td><td>21:10 / 2:18</td><td>42.0 / 46.0</td><td>12.0 / 16.0</td><td>69</td><td>69</td><td>0</td><td>0</td></tr><tr><td>20.08.2012</td><td>20:26 / 2:37</td><td>46.5 / 49.5</td><td>9.0 / 14.5</td><td>62</td><td>61</td><td>0</td><td>1</td></tr><tr><td>22.11.2013</td><td>1:31 / 5:54</td><td>41.0 / 42.5</td><td>4.0 / 8.0</td><td>62</td><td>10</td><td>50</td><td>2</td></tr><tr><td>25.12.2013</td><td>21:34 / 6:21</td><td>43.5 / 47.0</td><td>-7.0 / -1.0</td><td>61</td><td>14</td><td>47</td><td>5</td></tr><tr><td>09.11.2010</td><td>19:59 / 5:30</td><td>43.5 / 44.5</td><td>-5.5 / -0.5</td><td>58</td><td>1</td><td>57</td><td>0</td></tr><tr><td>02.07.2012</td><td>20:15 / 1:32</td><td>49.5 / 52.0</td><td>12.0 / 16.5</td><td>54</td><td>51</td><td>0</td><td>4</td></tr><tr><td>12.12.2009</td><td>22:26 / 3:05</td><td>41.0 / 42.0</td><td>6.0 / 7.0</td><td>54</td><td>47</td><td>2</td><td>3</td></tr><tr><td>20.06.2013</td><td>21:26 / 1:14</td><td>50.5 / 54.0</td><td>11.0 / 16.5</td><td>53</td><td>53</td><td>0</td><td>0</td></tr><tr><td>02.08.2009</td><td>21:15 / 1:00</td><td>49.0 / 51.0</td><td>15.0 / 19.0</td><td>52</td><td>52</td><td>0</td><td>0</td></tr><tr><td>09.01.2012</td><td>0:26 / 5:30</td><td>41.5 / 43.5</td><td>14.0 / 16.0</td><td>50</td><td>50</td><td>0</td><td>0</td></tr><tr><td>24.09.2012</td><td>17:47 / 22:24</td><td>45.0 / 47.5</td><td>12.5 / 18.0</td><td>50</td><td>50</td><td>0</td><td>0</td></tr><tr><td>02.07.2012</td><td>20:09 / 1:32</td><td>50.5 / 51.5</td><td>14.0 / 15.5</td><td>50</td><td>49</td><td>0</td><td>0</td></tr><tr><td>24.11.2013</td><td>18:55 / 5:31</td><td>39.0 / 41.0</td><td>6.5 / 7.5</td><td>49</td><td>8</td><td>39</td><td>4</td></tr><tr><td>28.11.2010</td><td>19:31 / 3:45</td><td>43.5 / 45.0</td><td>-6.0 / -1.0</td><td>49</td><td>0</td><td>49</td><td>0</td></tr><tr><td>06.11.2011</td><td>18:01 / 23:25</td><td>42.0 / 44.0</td><td>13.0 / 16.0</td><td>46</td><td>2</td><td>0</td><td>0</td></tr><tr><td>25.01.2010</td><td>22:00 / 22:00</td><td>35.0 / 37.0</td><td>-3.0 / 0.0</td><td>44</td><td>40</td><td>4</td><td>0</td></tr><tr><td>20.06.2013</td><td>20:38 / 0:43</td><td>51.5 / 53.0</td><td>14.0 / 15.0</td><td>43</td><td>38</td><td>0</td><td>5</td></tr><tr><td>20.08.2012</td><td>19:18 / 3:04</td><td>47.0 / 49.0</td><td>9.5 / 13.5</td><td>42</td><td>42</td><td>0</td><td>0</td></tr></tbody></table>