A New Approach of 3D Lightning Location Based on Pearson Correlation Combined with Empirical Mode Decomposition

To improve the accuracy of pulse matching and the mapping quality of lightning discharges, the Pearson correlation method combined with empirical mode decomposition (EMD) is introduced for discharge electric field pulse matching. This paper uses the new method to locate the lightning channels of an intra-cloud (IC) lightning flash and a cloud-to-ground (CG) lightning flash and analyzes the location results for the two lightning flashes. The results show that this method has a good performance in lightning location. Compared with the pulse-peak feature matching method, the positioning results of the new method are significantly improved, which is mainly due to the much larger number of positioning points (matched pulses). The number of located radiation sources has increased by nearly a factor of seven, which can significantly improve the continuity of the lightning channel and clearly distinguish the developmental characteristics. In the CG flash, there were three negative recoil streamers in the positive leader channel. After the three negative recoil streamers were finished, taking approximately 1 ms, 12 ms, and 2 ms, respectively, the negative leader channel underwent a K-process. The three negative recoil streamers are not connected to the K-processes in the negative leader channel. We think that the three negative recoil streamers may have triggered the three K-processes, respectively.

Artificial lightning strike tests and coupled electrical–thermal analysis on roofing glazed tiles of ancient buildings

Lightning strike is one of the natural disasters to the roof components of ancient buildings. To investigate the causes and damage effects of lightning strikes on the roofing glazed tiles of ancient buildings, artificial lightning strike tests were carried out on glazed tiles. Based on the experiment results, a coupled electrical–thermal finite element model of mortar-containing glazed tiles was established and the Joule heat effect of lightning current was further investigated. The results show that when the lightning channel is attached to the surface of the enamel and body with a low electrical conductivity, the lightning current is mainly released in the form of surface flashover, and a minor damage is induced along the flashover path; when the lightning channel is attached to the mortar with a high electrical conductivity, the lightning current is injected into the mortar, resulting in significant tile damage. The spatial distributions of the temperature present clear gradient characteristics. The high-temperature area appears in the mortar while the high–thermal–stress area appears in the body connected to the grounding rail. As the peak of the lightning current increases, both the high-temperature and high–thermal–stress areas of the glazed tiles expand. The combination of the experiments and the numerical analysis results demonstrate that the damage mechanism of lightning Joule heat effect to glazed tiles may include two aspects. One is the internal explosive force generated from the sharp vaporization and expansion of the moisture inside the tiles due to rapid temperature increase, and the other is the thermal stress caused by the uneven temperature distribution.

Computation of Lightning Current from Electric Field Based on Laplace Transform and Deconvolution Method

A method for computation of the lightning channel base current from the corresponding vertical component of lightning electric field was presented. The algorithm was developed by applying Laplace transform. The lightning current was estimated from its deconvolution with a special transfer function. The transfer function includes information about geometry and physical properties of entire lightning impulse generation system. The method was verified for a Heidler-type base current and a MTLL model of its propagation within the lightning channel. Research was done for close, middle, and far distance to the lightning strike point. Optimum performance was obtained for the middle distance of several kilometers where the electrostatic, induction, and radiation components of the transfer function were of the same range. An analysis was done for input electric field with and without noise superimposed on its time domain waveform. Relative uncertainties for the electric field and calculated lightning channel base current were similar each other. The presented approach can substantially increase a number of lightning current parameters which can be identified on the basis of its electric field signature. This method can be applied by the lightning location systems using preprocessing which increases the timing efficiency of the transfer function estimation.

Global ground strike point characteristics in negative downward lightning flashes – Part 1: Observations

Information about lightning properties is important in order to advance the current understanding of lightning, whereby the characteristics of ground strike points (GSPs) are in particular helpful to improving the risk estimation for lightning protection. Lightning properties of a total of 1174 negative downward lightning flashes are analyzed. The high-speed video recordings are taken in different regions, including Austria, Brazil, South Africa and the USA, and are analyzed in terms of flash multiplicity, duration, interstroke intervals and ground strike point properties. According to our knowledge this is the first simultaneous analysis of GSP properties in different regions of the world applying a common methodology. Although the results vary among the data sets, the analysis reveals that a third of the flashes are single-stroke events, while the overall mean number of strokes per flash equals 3.67. From the video imagery an average of 1.56 GSPs per flash is derived, with about 60 % of the multiple-stroke flashes striking the ground in more than one place. It follows that a ground contact point is struck 2.35 times on average. Multiple-stroke flashes last on average 371 ms, whereas the geometric mean (GM) interstroke interval value preceding strokes producing a new GSP is about 18 % greater than the GM value preceding subsequent strokes following a pre-existing lightning channel. In addition, a positive correlation between the duration and multiplicity of the flash is presented. The characteristics of the subset of flashes exhibiting multiple GSPs is further examined. It follows that strokes with a stroke order of 2 create a new GSP in 60 % of the cases, while this percentage quickly drops for higher-order strokes. Further, the possibility of forming a new lightning channel to ground in terms of the number of strokes that conditioned the previous lightning channel shows that approximately 88 % developed after the occurrence of only one stroke. Investigating the time intervals in the other 12 % of the cases when two or more strokes re-used the previous lightning channel showed that the average interstroke time interval preceding a new lightning channel is found to be more than twice the time difference between strokes that follow the previous lightning channel.

Modelling of tall-structure lightning return-stroke current using the Electromagnetic Transients Program

The main aim of this PhD work is to advance tall-structure lightning return-stroke current modelling. The Alternative Transients Program (ATP), a version of the Electromagnetic Transients program (EMTP), is used to model the lightning current distribution within a tall structure and the attached lightning channel. The tall structure, namely the CN Tower, is modeled as three or five transmission line sections connected in series. The lightning channel is represented by a transmission line with a continuously expanding length. The presented model takes into account reflections within the tower and within the lightning channel. Locations of reflections, current reflection coefficients and the parameters of the current simulation function are calculated based on the time analysis of the current derivative signal, measured at the tower. The decay parameters of the simulation function are first determined by curve fitting the decaying part of the current obtained from measurement. The other parameters are determined by curve fitting the measured initial current derivative impulse with the derivative of the simulation function, before the arrival of reflections. The simulation results substantially succeeded in reproducing the fine structure of the measured current derivative signal. The model allows for the computation of the lightning current at any point along the current path (the tower and the attached channel), which is required for the calculation of the associated electromagnetic field. Using the three-section model of the tower, the presented return-stroke current model enables the determination of a discrete return-stroke velocity profile, demonstrating that the velocity generally decays with time. Furthermore, based on the five-section model, the proposed approach enables taking into account the existence of upward-connecting leaders, which allowed, for the first time, the determination of upward-connecting leader lengths and return-stroke velocity variation profiles with more details. The return-stroke velocity profile is found to initially increase rapidly with time, reaching a peak, and then decrease less rapidly. The proposed model is also experimentally verified based on the comparison between the computed and measured electromagnetic fields. The simulated electric and magnetic field waveforms are found to reproduce important details of the measured fields, including initial split peaks that appear due to channel-front reflections in the presence of upward-connecting leaders.

Modelling of tall-structure lightning return-stroke current using the Electromagnetic Transients Program

The main aim of this PhD work is to advance tall-structure lightning return-stroke current modelling. The Alternative Transients Program (ATP), a version of the Electromagnetic Transients program (EMTP), is used to model the lightning current distribution within a tall structure and the attached lightning channel. The tall structure, namely the CN Tower, is modeled as three or five transmission line sections connected in series. The lightning channel is represented by a transmission line with a continuously expanding length. The presented model takes into account reflections within the tower and within the lightning channel. Locations of reflections, current reflection coefficients and the parameters of the current simulation function are calculated based on the time analysis of the current derivative signal, measured at the tower. The decay parameters of the simulation function are first determined by curve fitting the decaying part of the current obtained from measurement. The other parameters are determined by curve fitting the measured initial current derivative impulse with the derivative of the simulation function, before the arrival of reflections. The simulation results substantially succeeded in reproducing the fine structure of the measured current derivative signal. The model allows for the computation of the lightning current at any point along the current path (the tower and the attached channel), which is required for the calculation of the associated electromagnetic field. Using the three-section model of the tower, the presented return-stroke current model enables the determination of a discrete return-stroke velocity profile, demonstrating that the velocity generally decays with time. Furthermore, based on the five-section model, the proposed approach enables taking into account the existence of upward-connecting leaders, which allowed, for the first time, the determination of upward-connecting leader lengths and return-stroke velocity variation profiles with more details. The return-stroke velocity profile is found to initially increase rapidly with time, reaching a peak, and then decrease less rapidly. The proposed model is also experimentally verified based on the comparison between the computed and measured electromagnetic fields. The simulated electric and magnetic field waveforms are found to reproduce important details of the measured fields, including initial split peaks that appear due to channel-front reflections in the presence of upward-connecting leaders.

Lightning Detection and Imaging Based on VHF Radar Interferometry

In this study, detection and three-dimensional (3D) imaging of lightning plasma channels are presented using radar interferometry. Experiments were carried out in Leshan, China with a 48.2 MHz VHF radar configured with an interferometric antenna array. The typical characteristics of lightning echoes are studied in the form of amplitude, phase, and doppler spectra derived from the raw in-phase/quadrature (I/Q) data. In addition, the 3D structure of lightning channels is reconstructed using the interferometry technique. The localization results of lightning are verified with the locating results of lightning detection networks operating at VLF ranges, which indicate the feasibility of using VHF radar for lightning mapping. The interpretation of the observational results is complicated by the dendric structure of lightning channel and the overlap between passive electromagnetic radiations and return echoes. Nevertheless, some parts of the characteristics of lightning are still evident. The observational result of return echoes shows good consistency with the overdense assumption of lightning channels. The transition from the overdense channel to the underdense channel in the form of amplitude and phase is clearly observed. This technique is very promising to reveal the typical characteristics of lightning return echoes and structure of lightning propagation processes.

Artificial Lightning Strike Tests and Coupled Electrical–Thermal Analysis on Roofing Glazed Tiles of Ancient Buildings

Lightning strike is one of the natural disasters to the roof components of ancient buildings. To investigate the causes and damage effects of lightning strikes on the roofing glazed tiles of ancient buildings, artificial lightning strike tests were carried out on glazed tiles. Based on the experiment results, a coupled electrical–thermal finite element model of mortar-containing glazed tiles was established and the Joule heat effect of lightning current was further investigated. The results show that when the lightning channel is attached to the surface of the enamel and body with a low electrical conductivity, the lightning current is mainly released in the form of surface flashover, and a minor damage is induced along the flashover path; when the lightning channel is attached to the mortar with a high electrical conductivity, the lightning current is injected into the mortar, resulting in significant tile damage. The spatial distributions of the temperature present clear gradient characteristics. The high-temperature area appears in the mortar while the high–thermal–stress area appears in the body connected to the grounding rail. As the peak of the lightning current increases, both the high-temperature and high–thermal–stress areas of the glazed tiles expand. The combination of the experiments and the numerical analysis results demonstrate that the damage mechanism of lightning Joule heat effect to glazed tiles may include two aspects. One is the internal explosive force generated from the sharp vaporization and expansion of the moisture inside the tiles due to rapid temperature increase, and the other is the thermal stress caused by the uneven temperature distribution.