ELF/VLF signatures of sprite-producing lightning discharges observed during the 2005 EuroSprite campaign

During the summer of 2005, transient luminous events were optically imaged from the French Pyrénées as part of the EuroSprite campaign. Simultaneously, extremely low frequency (ELF: 3–3000 Hz) and broadband very low frequency (VLF: 3–30 kHz) data were recorded continuously at two separate receivers in Israel, located about 3300 km from the area of the parent lightning discharges responsible for the generation of sprites. Additionally, narrowband VLF data were collected in Crete, at about 2300 km away from the region of sprites.The motivation for the present study was to identify the signature of the sprite-producing lightning discharges in the ELF and VLF electromagnetic frequency bands, to qualify and compare their parameters, and to study the influence of the thunderstorm-activated region on its overlaying ionosphere. For the 15 sprites analyzed, their causative positive cloud-to-ground (+CG) discharges had peak current intensities between +8 and +130 kA whereas their charge moment changes (CMC) ranged from 500 to 3500 C km. Furthermore, the peak current reported by the Météorage lightning network are well correlated with the amplitudes of the VLF bursts, while showing poor correlation with the CMCs which were estimated using ELF methods.Additionally, more than one +CG was associated with six of the sprites, implying that lightning discharges that produce sprites can sometimes have multiple ground connections separated in time and space. Finally, for a significant number of events (33%) an ELF transient was not associated with sprite occurrence, suggesting that long continuing current of tens of ms may not always be a necessary condition for sprite production, a finding which influences the estimation of the global sprite rate based on Schumann resonance (SR) measurements.     

Analysis of lightning activity in two thunderstorm systems producing sprites in France

Two sprite-producing mesoscale convective systems above the South-Western part of France are studied. Three sprite events during the first night and seven during the second night were captured. Except for two events, the sprites could be associated with causative positive cloud-to-ground (+CG) discharges in the stratiform region of the storm. The analysis of lightning activity reveals that in both nights sprites occurred when lightning activity decreased rapidly and the ratio of +CGs to the total number of CGs decreased slightly. The average peak current of sprite-producing lightning was lower than 60 kA, in agreement with other observations. The delay times of sprites to their SPCGs varied from 57 to 140 m s and no correlation between events’ delay and shape was established.     

ELF transients associated with sprites and elves in eastern Mediterranean winter thunderstorms

During the northern hemisphere winter of 2005–2006, transient luminous events (TLEs) known as ‘sprites’ and ‘elves’ were imaged over thunderstorm cells in the eastern Mediterranean. Simultaneously, extremely low frequency (ELF) data (ELF: 3–3000 Hz) were recorded at two observation stations in Israel and Hungary in order to qualify and quantify parameters of the parent lightning discharge associated with the transient optical emissions in the upper atmosphere. In this study, we found that for 87% (Israel) and 77% (Hungary) of optically observed TLEs an intense ELF transient event was recorded. These stations are located some 500 and 2100 km, respectively, from the region of the TLEs. All ELF transients that were associated with TLEs were caused by lightning discharges with positive polarity. Calculation of the charge moment change showed values between 600 and 2800 C km with a peak around 1000 C km. Additionally, the time delay between the +CG and ensuing sprite was 76±34 ms and it was displaced up to 50 km from its parent CG.One of our objectives in the present study was to characterize, based on the ELF radiation from lightning, the electromagnetic (EM) waveforms of the lightning discharges which generate TLEs in the time and frequency domains, and to compare them with other lightning discharges occurring in the same thunderstorm cell at approximately the same time, but which did not produce TLEs. The survey for a typical EM waveform showed no unique ELF signature for lightning discharges associated with either sprites or elves.     

Comparison of time delays of sprites induced by winter lightning flashes in the Japan Sea with those in the Pacific Ocean

The measurement of unusual winter sprites in the Hokuriku area (Japan Sea side) was performed as a primary target of the 2006/2007 winter campaign by means of coordinated optical and extremely low frequency (ELF)/very high-frequency (VHF) electromagnetic observations. We have also added the same observations for the sprites in the Pacific Ocean, to be compared with the characteristics of Hokuriku sprites. The following results have emerged from this campaign: (i) the predominance of column sprites in winter has been confirmed not only for the Hokuriku area but also in the Pacific Ocean (with the probability just above 60%), (ii) carrots are much more frequently observed in the Pacific Ocean (with a probability of ～28%) than in the Hokuriku area (～16%), (iii) a very unique property of Hokuriku sprites is the surprisingly long delay (average ～90 ms) of sprites from their parent lightning flashes and the delays for carrots and columns exhibit some significant difference (80 ms for columns and 100 ms for carrots) and (iv) the time delay of Pacific Ocean sprites is much shorter (～43 ms average) than that at Hokuriku, but there is no remarkable difference in delay between carrots and columns. Finally we discussed the importance of time delay studies to understand sprite generations and their parent lightning discharges, because the difference of time delays on the Japan Sea side and in the Pacific Ocean are thought to be causally related to the parameters of parent thunderstorms.     

A study of the morphology of winter sprites in the Hokuriku area of Japan in relation to cloud charge height

Continuous observations of sprites in the Hokuriku area of Japan were performed from two optical sites during the three winter periods. The purpose of this observation is to study the major effect in the appearance of sprites and in determining the morphology of sprites (columns or carrots). Detailed analysis is performed based on the estimation of the height of ?10 °C at the time of sprite occurrence. When the height of ?10 °C is lower than 1800 m, the occurrence of sprites is infrequent, and the dominant shape is column. Then when it is increased (1800–3000 m), a new situation takes place, namely the occurrence of sprites is very enhanced and more spectacular shapes like carrots tend to be frequently observed in addition to column sprites. These sprite characteristics are first compared with those of parent lightning in the Hokuriku area and with our latest computer simulations on sprite initiation.     

Attempts to create ball lightning with triggered lightning

We describe attempts to create ball lightning by directing lightning, triggered from natural thunderclouds using the rocket-and-wire technique, through a variety of materials. Some of the observed phenomena have features in common with natural ball lightning or with laboratory attempts to create it: flame-like luminosity for up to 0.5 s above salt water; constant-luminosity silicon fragments falling for about 1 s under the influence of gravity; a 0.7 m region of stationary luminosity whose bottom was 0.3 m above a stainless steel surface to which arcing had occurred; and a glow for about 0.5 s above pine tree sections.     

Using full-flash narrowband energy for ranging of lightning ground strokes

We demonstrate that narrowband measurements can be used for rudimentary ranging of cloud-to-ground lightning flashes. The system at present responds to both intra-cloud and cloud-to-ground lightning; ranging is demonstrated for a subset of flashes known to be cloud-to-ground lightning. The system uses a ferrite-core antenna with a length of about 4 cm and diameter 4 mm, and operates on a narrow band at about 1 MHz, close to the HF band (3–30 MHz). It downmixes the signal to audio frequencies and operates in a manner which is very similar to an AM radio. The system triggers on all impulses which exceed a given adjustable threshold above the ambient noise level, and records 1 s of data. Such a system was used to collect lightning-caused electromagnetic disturbances during summer 2006 in Finland. The output is compared to two scientifically verified references: a flat-plate broadband antenna measuring the vertical electric field and a commercial lightning location network giving flash location. A key aim of the system is to reduce the information to as few parameters as possible. Peak intensity and full-flash energy were used as simple parameters. It is shown that accurate flash-by-flash ranging is not possible with this method; however, it is shown that the method can be used to track clusters of ground flashes within a range of about 50–100 km with an accuracy of about 10 km.     

Some results of observations of positive lightning discharges and relative phenomena in the east of Siberia

We have considered spatial distributions of positive lightning discharges in the east of Siberia for the summer seasons of 2003–2007 and properties of their electromagnetic signals with the ELF “slow tail”, which, as is known, can be accompanied by sprites. There are two main regions of positive discharges located in the south and west of Yakutsk. Two other “centers” (the northeastern and the eastern) are located in high-mountainous regions. In these regions the positive discharges intensity can exceed the negative discharge intensity.The electromagnetic signals in the ELF range (usually in the form of two half-cycles) were observed after the VLF atmospherics were recorded in the high-latitude regions. The delay of ELF pulses relative to the corresponding atmospherics was 0–7 ms. The long (up to 350 ms) events of quasi-periodic ELF oscillations with the period of about 7 ms (which corresponds to the quasi-period of ELF pulses) were revealed.     

Ten-year study of cloud-to-ground lightning activity in the Iberian Peninsula

The characteristics (annual and diurnal cycle, polarity, multiplicity and first stroke peak current) of ∼4.3×10⁶ cloud-to-ground flashes recorded in the Iberian Peninsula during the first decade of measurements of the lightning detection network installed in Spain are analyzed. The mean monthly variation shows maximum lightning activity between May and September, while minimum values are observed in January and February. The mean diurnal cycle shows maximum values at 1700 LT and minimum values at around 1000 LT. The average maximum flash density (not corrected for detection efficiency) is 2.1 flashes km⁻² year⁻¹. Maximum lightning activity is associated with mountainous areas. The effect of the Mediterranean Sea is also seen. The percentage of positive flashes is 9%, although this changes over the year from 6.5% in June to 22.6% in January. The average multiplicity is found to be 2.0 for the negative flashes and 1.1 for the positive flashes, and the percentages of single-stroke flashes are 53.6% and 89%, respectively. The monthly distribution of multiplicity for negative flashes peaks in the summer and minimum values are found in the winter. The multiplicity of the positive flashes does not seem to be function of the month. The median (mean) first stroke peak currents are found to be 23.5 kA (27.3 kA) for the negative flashes and 35.3 kA (47.1 kA) for the positive flashes. For both polarities, the peak current is higher in the summer than in the winter. The percentage of positive flashes and the mean peak currents for negative flashes are higher over the sea areas than over land.     

Dynamic soil–structure interaction analysis of a telescope at the Javalambre Astrophysical Observatory

This paper presents the dynamic soil–structure analysis of the main telescope T250 of the Observatorio Astrofísico de Javalambre (OAJ, Javalambre Astrophysical Observatory) on the Pico del Buitre. Vibration control has been of prime concern in the design, since astrophysical observations may be hindered by mechanical vibration of optical equipment due to wind loading. The telescope manufacturer therefore has imposed a minimal natural frequency of 10 Hz for the supporting telescope pier. Dynamic soil–structure interaction may significantly influence the lowest natural frequency of a massive construction as a telescope pier. The structure clamped at its base has a resonance frequency of 14.3 Hz. A coupled finite element–boundary element (FE–BE) model of the telescope pier that accounts for the dynamic interaction of the piled foundation and the soil predicts a resonance frequency of 11.2 Hz, demonstrating the significant effect of dynamic soil–structure interaction. It is further investigated to what extent the coupled FE–BE model can be simplified in order to reduce computation time. The assumption of a rigid pile cap allows us to account for dynamic soil–structure interaction in a simplified way. A coupled FE–BE analysis with a rigid pile cap predicts a resonance frequency of 11.7 Hz, demonstrating a minor effect of the pile cap flexibility on the resonance frequency of the telescope pier. The use of an analytical model for the pile group results in an overestimation of the dynamic soil stiffness. This error is due to the large difference between the actual geometry and the square pile cap model for which the parameters have been tuned.     

Characteristics of Japanese winter sprites and their parent lightning as estimated by VHF lightning and ELF transients

Sprites are newly discovered optical emissions in the mesosphere over large thunderstorms. This paper is the observational summary of winter sprites in the Hokuriku area of Japan and their parent lightning in the winter of 2004/2005, by using the coordinated optical and electromagnetic (VHF and ELF) measurements in Japan. As the results of optical observations at two stations, we have found that this campaign has yielded a variety of sprite shapes; V-angle shaped structures have been often observed (25%) in addition to columnar structures familiar for us. All of the sprite events are found to be associated with $+$CG lightning, as seen from the macroscopic information by ELF data at Moshiri. However, examining the microscopic properties of parent lightning as seen from the VHF SAFIR lightning detection network, has suggested very complicated characteristics of parent lightning discharges inducing sprites, as compared with the ELF data. One half of the sprite events are also found to be associated with $+$CG by the SAFIR observation, but another half has yielded rather new results as compared with earlier results. Four events are definitely associated with $-$CG and the remaining three events, inter-cloud flashes. The overall picture for Japanese winter sprites and their parent lightning discharges, is significantly different from that for the summer-time, continental sprites. This is indicative of complexity of winter lightning in the Hokuriku area of Japan and this would provide new information on the sprite generation mechanism.     

The effects of thunderstorm-generated atmospheric gravity waves on mid-latitude F-region drifts

Superposed epoch analysis (SEA) was used to examine ionospheric drift velocities measured by a digital ionosonde located at the mid-latitude station Bundoora (145.1°E, 37.7°S geographic), near Melbourne. The control times for the SEA were the times of cloud-to-ground (CG) lightning strokes measured from August 2003 to August 2004 by the World Wide Lightning Location Network (WWLLN). Statistically, regions of concentrated lightning activity migrated from west to east across Bundoora, and the stroke frequency was higher the day prior the activity reached the station, and lower on the day after it passed to the east. For the SEA, CG strokes were separated into four directional quadrants centred on north, south, east and west. No SEA results are shown for the south quadrant due to the relatively low detection frequency of strokes across the Southern Ocean (6% of all events). The strongest downward vertical perturbations in F-region drifts, ?4.5 m s^?1, were found for lightning located towards the west during ?30 to ?16 h (i.e., the afternoon prior the activity passed near the station at t=0 h). The downward perturbation decreased in amplitude to ?1.5 m s^?1 for lightning located towards the north during ?6–+6 h, and was weakest (?0.7 m s^?1) for lightning located towards the east during +16–+28 h (i.e., the next afternoon). There were directionally consistent perturbations in the drift azimuths associated with the lightning located in their respective quadrants; lightning located to the west of the station caused eastward azimuth enhancements, northward lightning caused southward enhancements, and eastward lightning caused westward enhancements. Velocity magnitudes and fluctuations tended to increase during the passage of lightning. The observed responses were stronger when the SEA was performed with data selected using time windows of <2 min on either side of each lightning stroke. However, they persisted at longer time scales and were strong when thunderstorm onsets (instead of lightning times) were used as controls. Our results can be explained by thunderstorm-generated atmospheric gravity waves (AGWs) which subsequently gave rise to medium-scale travelling ionospheric disturbances (MSTIDs), with the lightning strokes acting merely as a proxy for this coupling. The prevailing thermospheric winds were flowing from east to west across the study region, and may have acted as a directional ‘filter’ for the MSTIDs, allowing waves generated in the west quadrant to reach the station and preventing those generated in other quadrants. Displacement of the MSTIDs in the direction anti-parallel to mean neutral wind flow has been observed by (Waldock, J.A., Jones, T.B., 1986. HF Doppler observations of medium-scale travelling ionospheric disturbances at mid-latitudes. Journal of atmospheric and terrestrial physics 48(3), 245–260).     

On the electric field at the tip of dart leaders in lightning flashes

The results obtained in this study show that as the dart leader tip passes a given point on the defunct return stroke channel the electric field increases within a fraction of a microsecond to values larger than the critical electric field necessary for the initiation of cold electron runaway in low-density air comprising the channel. These results are in support of the hypothesis that cold runaway electron breakdown may play a role in the emission of X-ray bursts by dart leaders. The calculations also show that the peak power dissipated by a typical dart leader is about 300–500 MW/m and the energy dissipated within the first 10 μs or so is about 500–600 J/m. Furthermore, the minimum resistance and the maximum radius of the core of a typical dart leader are estimated to be about 3 Ω/m and 0.003 m, respectively.     

Thunderstorm-associated responses in the vertical motion of the mid-latitude F-region ionosphere

Mid-latitude Digisonde Doppler velocities, auroral electrojet (AE) indices and cloud-to-ground (CG) lightning strokes during August 2003–2004 were used to study the perturbations in the F-region vertical drift associated with terrestrial thunderstorms. A superposed epoch analysis (SEA) showed that the F-region vertical drifts V_z had a net descent of ～0.6 m s^?1 peaking ～3 h after lightning. Stronger downward perturbations of up to ～0.9 m s^?1 were observed in the afternoon on the day prior to lightning days. The perturbations were less significant on the day after and insignificant during the remaining intervals up to 144 h on either side of the lightning. The stronger responses on the day before are consistent with causality because the lightning times were merely proxies for the physical mechanisms involved. The actual causes are unclear, but we discuss the possible roles of lightning-induced ionisation enhancements, intense electric fields penetrating upward from electrified clouds, and atmospheric gravity waves (AGWs) radiated from thunderstorms or from the accompanying tropospheric fronts. There is no doubt that the behaviour of the mid-latitude F-region is controlled by the thermospheric winds and the solar wind-magnetosphere electrical generators, but our results suggest that electrified clouds also account for a significant, albeit relatively small component of the ionospheric variability.     

Narrow positive bipolar radiation from lightning observed in Sri Lanka

Narrow positive bipolar pulses (NPBPs), whose origin largely remains unknown as yet, have occasionally been noticed in Sri Lanka. These discharges are found to be opposite in polarity to that of negative return strokes, and are found to occur at the beginning and active stage of thunderstorm activities in Sri Lanka. They are emitted from the thunderstorms that produce other activities also and are relatively narrower, bipolar and isolated in nature. They are neither preceded by leader-type pulses nor succeeded by the subsequent activities and hence, could not be associated with any other known activity. Similar radiations have previously been identified and termed as NPBPs. Their features are indeed unique and different from the other known thunderstorm electrical processes. These events have been found to be accompanied by the HF radiations at 5 and 10 MHz, at the trailing part of the pulse. However, such narrow, bipolar and isolated events have not been observed in Sweden, during a lightning measurement campaign at Uppsala, in 2006. So, it is speculated that the meteorological conditions to be responsible for the thunderstorm activity to produce such pulses. Features of such pulses recorded on different thunderstorm days in 2005 and 2006 have been analyzed and presented in this study.The average rise time (10–90%) (Tr), of the pulses was found to be 2.6 μs, the average zero crossing time (Tz) was found to be 5.85 μs, the average duration of slow front (Ts) was found to be 1.86 μs, and the average ratio of amplitude of overshoot to the corresponding peak amplitude (Os/Pa) of these pulses was found to be 0.39.     

Estimation of the Free Core Nutation parameters from SG data: Sensitivity study and comparative analysis using linearized least-squares and Bayesian methods

The Free Core Nutation (FCN) is investigated with the help of its resonance effect on the tidal amplitudes in Superconducting Gravimeter (SG) records of the GGP network. The FCN resonance parameters are combined in a resonance equation involving the Earth's interior parameters. The sensitivity of the FCN parameters to the diurnal tidal waves demonstrates that the quality factor of the FCN is strongly dependent on the accuracy of the imaginary part estimates of the gravimetric factors close to the resonance. The weak amplitude of Ψ₁ tidal wave on the Earth, which is the closest in frequency to the FCN, in addition to errors in ocean loading correction, explains the poor determination of the quality factor Q from surface gravimetric data. The inversion of tidal gravimetric factors leads to estimates of the period, Q and resonance strength of the FCN. We show that, by inverting log(Q) instead of Q, the results using the least-squares method optimized using the Levenberg–Marquardt algorithm are in agreement with the Bayesian probabilistic results and agree with the results obtained from VLBI nutation data. Finally, a combined inversion of 7 GGP European SG data is performed giving T = 428 ± 3 days and 7762 < Q < 31,989 (90% C.I.). An experimental estimate of the internal pressure Love number is also proposed.     

Sensitivity of ground-motion scenarios to earthquake source parameters in the Tehran metropolitan area,Iran

The northern Tehran fault (NTF) is a principal active fault of the Alborz mountain belt in the northern Iran. The fault is located north of the highly populated Metropolitan Area of Tehran. Historical records and paleoseismological studies have shown that the NTF poses a high seismic risk for the Tehran region and the surrounding cities (e.g. Karaj). A series of ground-motion simulations are carried out using a hybrid kinematic-stochastic model to calculate broadband (0.1–20 Hz) ground-motion time histories for deterministic earthquake scenarios (M7.2) on the NTF. We will describe the source characteristics of the target event to develop a list of scenario earthquakes that are probably similar to a large earthquake on the NTF. The effect of varying different rupture parameters such as rupture velocity and rise time on the resulting broadband strong motions has been investigated to evaluate the range of uncertainty in seismic scenarios. The most significant parameters in terms of ground-shaking level are the rise time and the value of the rupture velocity. For the worst-case scenario, the maximum expected horizontal acceleration, and velocity at rock sites in Tehran range between 128 and 1315 cm/s/s and 11–191 cm/s, respectively. For the lowest scenario, the corresponding values range between 102 and 776 cm/s/s and 12 to 81 cm/s. Nonlinear soil effects may change these results but are not accounted for in this study. The largest variability of ground motion is observed in neighborhood of asperity and also in the direction of rupture propagation. The calculated standard deviation of all ground-motion scenarios is less than 30% of the mean. The capability of the simulation method to synthesize expected ground motions and the appropriateness of the key parameters used in the simulations are confirmed by comparing the synthetic peak ground motions (PGA, PGV and response spectra) with empirical ground-motion prediction equations.     

Site response analysis of Vartholomio W-Greece from singular spectrum analysis of microtremor and weak motion data

Twenty six sites were instrumented in the city of Vartholomio following the December 2, 2002 Ms 6.0 earthquake. Thirty weak events from the aftershock sequence as well as microtremors were used to identify amplifications due to geological site effects. Horizontal-to-vertical spectral ratios (HVSR—Nakamura estimates) and weak events ratios were calculated and the singular spectrum analysis (SSA) method was used. The results showed that the effects of SSA on the stability of the frequency peak and amplitude distribution of HVSR for both weak motion and microtremors. The data analysis confirms the role of near surface geology in causing locally significant variations of the predominant frequencies and amplitudes of ground shaking as already inferred from the distribution of damages. The site response spectra exhibited significant peaks within the range of 1.5–2.6 Hz and the amplification factor did not exceed 6.5. Finally the parts of the HVSR ratios from ～0.2 up to 10 Hz were used, in order to create an automatic optimal zonation of the study area using a genetic algorithm. This procedure resulted in the division of the city into 2 main zones.     

1-g Experimental investigation of bi-layer soil response and kinematic pile bending

The effect of soil inhomogeneity and material nonlinearity on kinematic soil–pile interaction and ensuing bending under the passage of vertically propagating seismic shear waves in layered soil, is investigated by means of 1-g shaking table tests and nonlinear numerical simulations. To this end, a suite of scale model tests on a group of five piles embedded in two-layers of sand in a laminar container at the shaking table facility in BLADE Laboratory at University of Bristol, are reported. Results from white noise and sine dwell tests were obtained and interpreted by means of one-dimensional lumped parameter models, suitable for inhomogeneous soil, encompassing material nonlinearity. A frequency range from 0.1 Hz to 100 Hz and 5 Hz to 35 Hz for white noise and sine dwell tests, respectively, and an input acceleration range from 0.015 g to 0.1 g, were employed. The paper elucidates that soil nonlinearity and inhomogeneity strongly affect both site response and kinematic pile bending, so that accurate nonlinear analyses are often necessary to predict the dynamic response of pile foundations.     

Observations of narrow bipolar events in East China

Narrow bipolar events (NBEs) are a distinct class of intra-cloud lightning discharge. In this paper we present observations of 10 negative and 67 positive such events in East China. Positive NBEs occurred at 7–12 km altitude above mean sea level (MSL) with a mean altitude of 9.5 km, and negative NBEs occurred at 14–16 km altitude. Electrical/channel characteristics of these events were derived from NBE pulse waveforms based on the transmission-line model. On average, the peak current moment and the charge moment change of a NBE event is 15 kA km, and 0.12 C km, respectively. The mean time for the propagation of current front along the channel is 2.2 μs. The upper limit on channel length for NBEs in this study is 510–1060 m, the lower limit on discharge current amplitude is 12.5–43.2 kA, and the minimum charge transfer is 0.1–0.3 C.