skip to main content


Title: Application of Ensemble Empirical Mode Decomposition in Low-Frequency Lightning Electric Field Signal Analysis and Lightning Location
The application of empirical mode decomposition (EMD) in the analysis and processing of lightning electric field waveforms acquired by the low-frequency e-field detection array (LFEDA) in China has significantly improved the capabilities of the low-frequency/very-low-frequency (LF/VLF) time-of-arrival technique for studying the lightning discharge processes. However, the inherent mode mixing and the endpoint effect of EMD lead to certain problems, such as an inadequate noise reduction capability, the incorrect matching of multistation waveforms, and the inaccurate extraction of pulse information, which limit the further development of the LFEDA's positioning ability. To solve these problems, the advanced ensemble EMD (EEMD) technique is introduced into the analysis of LF/VLF lightning measurements, and a double-sided bidirectional mirror (DBM) extension method is proposed to overcome the endpoint effect of EMD. EEMD can effectively suppress mode mixing, and the DBM extension method proposed in this article can effectively suppress the endpoint effect, thus greatly improving the accuracy of a simulated signal after a 25-500-kHz bandpass filter. The resulting DBM_EEMD algorithm can be used in the LFEDA system to process and analyze the detected electric field signals to improve the system's lightning location capabilities, especially in terms of accurate extraction and location of weak signals from lightning discharges. In this article, a 3-D image of artificially triggered lightning obtained from an LF/VLF location system is reported for the first time, and methods for further improving the location capabilities of the LF/VLF lightning detection systems are discussed.  more » « less
Award ID(s):
1720600
NSF-PAR ID:
10197472
Author(s) / Creator(s):
; ; ; ; ; ; ; ;
Date Published:
Journal Name:
IEEE Transactions on Geoscience and Remote Sensing
ISSN:
0196-2892
Page Range / eLocation ID:
1 to 15
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Multi‐resolution analysis methods can reveal the underlying physical dynamics of nonstationary signals, such as those from lightning. In this paper we demonstrate the application of two multi‐resolution analysis methods: Ensemble Empirical Mode Decomposition (EEMD) and Variational Mode Decomposition (VMD) in a comparative way in the analysis of electric field change waveforms from lightning. EEMD and VMD decompose signals into a set of Intrinsic Mode Functions (IMFs). The IMFs can be combined using distance and divergence metrics to obtain noise reduction or to obtain new waveforms that isolate the physical processes of interest while removing irrelevant components of the original signal. We apply the EEMD and VMD methods to the observations of three close Narrow Bipolar Events (NBEs) that were reported by Rison et al. (2016,https://doi.org/10.1038/ncomms10721). The ΔE observations reveal the occurrence of complex oscillatory processes after the main NBE sferic. We show that both EEMD and VMD are able to isolate the oscillations from the main NBE, with VMD being more effective of the two methods since it requires the least user supervision. The oscillations are found to begin at the end of the NBEs' downward fast positive breakdown, and appear to be produced by a half‐wavelength standing wave within a weakly‐conducting resonant ionization cavity left behind in the wake of the streamer‐based NBE event. Additional analysis shows that one of the NBEs was likely initiated by an energetic cosmic ray shower, and also corrects a misinterpretation in the literature that fast breakdown is an artifact of NBE‐like events in interferometer observations.

     
    more » « less
  2. Abstract

    The lightning data products generated by the low‐frequency (LF) radio lightning locating system (LLS) deployed during the Remote sensing of Electrification, Lightning, and Mesoscale/Microscale Processes with Adaptive Ground Observation (RELAMPAGO) field campaign in Argentina provide a valuable data set to research the lightning evolution and characteristics of convective storms that produce high‐impact weather. LF LLS data sets offer a practical range for mesoscale studies, allowing for the observation of lightning characteristics of storms such as mesoscale convective systems or large convective lines that travel longer distances which are not necessarily staying in range of regional VHF‐based lightning detection systems throughout their lifetime. LF LLSs also provide different information than optical space‐borne lightning detectors. Lightning measurements exclusive to LF systems include discharge peak current, lightning polarity, and lightning type classification based on the lightning‐emitted radio waveform. Furthermore, these measurements can provide additional information on flash rates (e.g., positive cloud‐to‐ground flash rate) or narrow bipolar events which may often be associated with dynamically intense convection. In this article, the geolocation and data processing of the LF data set collected during RELAMPAGO is fully described and its performance characterized, with location accuracy better than 10 km. The detection efficiency (DE) of the data set is compared to that of the Geostationary Lightning Mapper, and spatiotemporal DE losses in the LF data set are discussed. Storm case studies on November 10, 2018, highlight the strengths of the data set, which include robust flash clustering and insightful flash rate and peak current measures, while illustrating how its limitations, including DE losses, can be managed.

     
    more » « less
  3. Abstract

    A particular strength of lightning remote sensing is the variety of lightning types observed, each with a unique occurrence context and characteristically different emission. Distinct energetic intracloud (EIC) lightning discharges—compact intracloud lightning discharges (CIDs) and energetic intracloud pulses (EIPs)—produce intense RF radiation, suggesting large currents inside the cloud, and they also have different production mechanisms and occurrence contexts. A Low‐Frequency (LF) lightning remote sensing instrument array was deployed during the RELAMPAGO field campaign in west central Argentina, designed to investigate convective storms that produce high‐impact weather. LF data from the campaign can provide a valuable data set for researching the lightning context of EICs in a variety of subtropical convective storms. This paper describes the production of an LF‐CID data set in RELAMPAGO and includes a preliminary analysis of CID prevalence. Geolocated lightning events and their corresponding observed waveforms from the RELAMPAGO LF data set are used in the classification of EICs. Height estimates based on skywave reflections are computed, where prefit residual data editing is used to improve robustness against outliers. Even if EIPs occurred within the network, given the low number of very high‐peak current events and receiver saturation, automatic classification of EIPs may not be feasible using this data. The classification of CIDs, on the other hand, is straightforward and their properties, for both positive and negative polarity, are investigated. A few RELAMPAGO case studies are also presented, where high variability of CID prevalence in ordinary storms and high‐altitude positive CIDs, possibly in overshooting tops, are observed.

     
    more » « less
  4. One of the most basic pieces of information gained from dynamic electromyography is accurately defining muscle action and phase timing within the gait cycle. The human gait relies on selective timing and the intensity of appropriate muscle activations for stability, loading, and progression over the supporting foot during stance, and further to advance the limb in the swing phase. A common clinical practice is utilizing a low-pass filter to denoise integrated electromyogram (EMG) signals and to determine onset and cessation events using a predefined threshold. However, the accuracy of the defining period of significant muscle activations via EMG varies with the temporal shift involved in filtering the signals; thus, the low-pass filtering method with a fixed order and cut-off frequency will introduce a time delay depending on the frequency of the signal. In order to precisely identify muscle activation and to determine the onset and cessation times of the muscles, we have explored here onset and cessation epochs with denoised EMG signals using different filter banks: the wavelet method, empirical mode decomposition (EMD) method, and ensemble empirical mode decomposition (EEMD) method. In this study, gastrocnemius muscle onset and cessation were determined in sixteen participants within two different age groups and under two different walking conditions. Low-pass filtering of integrated EMG (iEMG) signals resulted in premature onset (28% stance duration) in younger and delayed onset (38% stance duration) in older participants, showing the time-delay problem involved in this filtering method. Comparatively, the wavelet denoising approach detected onset for normal walking events most precisely, whereas the EEMD method showed the smallest onset deviation. In addition, EEMD denoised signals could further detect pre-activation onsets during a fast walking condition. A comprehensive comparison is discussed on denoising EMG signals using EMD, EEMD, and wavelet denoising in order to accurately define an onset of muscle under different walking conditions. 
    more » « less
  5. Abstract

    The very‐low frequency (VLF) and low frequency (LF) waves from ground transmitters propagate in the ionospheric waveguide, and a portion of their power leaks to the Earth's inner radiation belt and slot region where it can cause electron precipitation loss. Using Van Allen Probes observations, we perform a survey of the VLF and LF transmitter waves at frequencies from 14 to 200 kHz. The statistical electric and magnetic wave amplitudes and frequency spectra are obtained at 1 < L < 3. Based on a recent study on the propagation of VLF transmitter waves, we divide the total wave power into ducted and unducted portions, and model the wave normal angle of unducted waves with dependences onLshell, magnetic latitude, and wave frequency. At lower frequencies, the unducted waves are launched along the vertical direction and the wave normal angle increases during the propagation until reaching the Gendrin angle; at higher frequencies, the normal angle of unducted waves follows the variation of Gendrin angle. We calculate the bounce‐averaged pitch angle and momentum diffusion coefficients of electrons due to ducted and unducted VLF and LF waves. Unducted and ducted waves cause efficient pitch angle scattering atL = 1.5 and 2.5, respectively. Although the wave power from ground transmitters at frequencies higher than 30 kHz is low, these waves can cause the pitch angle scattering of lower energy (2–200 keV atL = 1.5) electrons, which cannot resonate with the VLF transmitter waves at frequencies below 30 kHz, lightning generated whistlers, or plasmaspheric hiss.

     
    more » « less