Electromagnetic scattering in clouds

Techniques used to explain the nature of the optical effects of clouds on the light produced by lightning include a Monte Carlo simulation, an equivalent medium approach, and methods based on Boltzmann transport theory. A cuboidal cloud has been considered using transform methods and a diffusion approximation. Many simplifying assumptions have been used by authors to make this problem tractable. In this report, the cloud will have a spherical shape and its interior will consist of a uniform distribution of identical spherical water droplets. The source will be modeled as a Hertz dipole, electric or magnetic, inside or outside the cloud. An impulsive source is used. Superposition may be employed to obtain a sinusoid within an envelope which describes a lightning event. The problem is investigated by transforming to the frequency domain, obtaining Green's functions, and then using the Cagniard-DeHoop method to symbolically recover the time domain solution

A geometrical result regarding time-of-arrival lightning location

One reason for investigating Lightning Detection And Ranging (LDAR) is to validate data from the Optical Transient Detector (OTD). A Time-Of-Arrival (TOA) procedure may be used with radio wave portions of lighting signatures. An antenna is in place at KSC

Scattering of point source illumination by an arbitrary configuration

The problem of electromagnetic scattering of an incident plane wave by an arbitrary configuration of obstacles was solved by Twersky. In this report, the results are extended to point source incidence corresponding to a Hertz dipole. Knowledge of the response of a fixed configuration of scatterers excited by a point source may provide insight to improve the accuracy of the values of bulk parameters for clouds which have been found using plane wave excitation

The Ground Flash Fraction Retrieval Algorithm Employing Differential Evolution: Simulations and Applications

The ability to estimate the fraction of ground flashes in a set of flashes observed by a satellite lightning imager, such as the future GOES-R Geostationary Lightning Mapper (GLM), would likely improve operational and scientific applications (e.g., severe weather warnings, lightning nitrogen oxides studies, and global electric circuit analyses). A Bayesian inversion method, called the Ground Flash Fraction Retrieval Algorithm (GoFFRA), was recently developed for estimating the ground flash fraction. The method uses a constrained mixed exponential distribution model to describe a particular lightning optical measurement called the Maximum Group Area (MGA). To obtain the optimum model parameters (one of which is the desired ground flash fraction), a scalar function must be minimized. This minimization is difficult because of two problems: (1) Label Switching (LS), and (2) Parameter Identity Theft (PIT). The LS problem is well known in the literature on mixed exponential distributions, and the PIT problem was discovered in this study. Each problem occurs when one allows the numerical minimizer to freely roam through the parameter search space; this allows certain solution parameters to interchange roles which leads to fundamental ambiguities, and solution error. A major accomplishment of this study is that we have employed a state-of-the-art genetic-based global optimization algorithm called Differential Evolution (DE) that constrains the parameter search in such a way as to remove both the LS and PIT problems. To test the performance of the GoFFRA when DE is employed, we applied it to analyze simulated MGA datasets that we generated from known mixed exponential distributions. Moreover, we evaluated the GoFFRA/DE method by applying it to analyze actual MGAs derived from low-Earth orbiting lightning imaging sensor data; the actual MGA data were classified as either ground or cloud flash MGAs using National Lightning Detection Network[TM] (NLDN) data. Solution error plots are provided for both the simulations and actual data analyses

Analytic Perturbation Method for Estimating Ground Flash Fraction from Satellite Lightning Observations

An analytic perturbation method is introduced for estimating the lightning ground flash fraction in a set of N lightning flashes observed by a satellite lightning mapper. The value of N is large, typically in the thousands, and the observations consist of the maximum optical group area produced by each flash. The method is tested using simulated observations that are based on Optical Transient Detector (OTD) and Lightning Imaging Sensor (LIS) data. National Lightning Detection NetworkTM (NLDN) data is used to determine the flashtype (ground or cloud) of the satellite-observed flashes, and provides the ground flash fraction truth for the simulation runs. It is found that the mean ground flash fraction retrieval errors are below 0.04 across the full range 01 under certain simulation conditions. In general, it is demonstrated that the retrieval errors depend on many factors (i.e., the number, N, of satellite observations, the magnitude of random and systematic measurement errors, and the number of samples used to form certain climate distributions employed in the model)

Tests of the Grobner Basis Solution for Lightning Ground Flash Fraction Retrieval

Satellite lightning imagers such as the NASA Tropical Rainfall Measuring Mission Lightning Imaging Sensor (TRMM/LIS) and the future GOES-R Geostationary Lightning Mapper (GLM) are designed to detect total lightning (ground flashes + cloud flashes). However, there is a desire to discriminate ground flashes from cloud flashes from the vantage point of space since this would enhance the overall information content of the satellite lightning data and likely improve its operational and scientific applications (e.g., in severe weather warning, lightning nitrogen oxides studies, and global electric circuit analyses). A Bayesian inversion method was previously introduced for retrieving the fraction of ground flashes in a set of flashes observed from a satellite lightning imager. The method employed a constrained mixed exponential distribution model to describe the lightning optical measurements. To obtain the optimum model parameters (one of which is the ground flash fraction), a scalar function was minimized by a numerical method. In order to improve this optimization, a Grobner basis solution was introduced to obtain analytic representations of the model parameters that serve as a refined initialization scheme to the numerical optimization. In this study, we test the efficacy of the Grobner basis initialization using actual lightning imager measurements and ground flash truth derived from the national lightning network

On the Retrieval of Lightning Radio Sources from Time-of-Arrival Data

We examine the problem of retrieving three-dimensional lightning locations from radio frequency Time-Of-Arrival (TOA) measurements. Arbitrary antenna locations are considered. By judiciously differencing measurements that are related to the location of the antennas and their excitation times, the problem is converted from the initial spherical nonlinear form to a system of linear equations. In the linear formalism, the source location and time-of-occurrence is viewed geometrically as an intersection of hyperplanes in the four-dimensional Minkowski space (x,y,z,t). The linear equations are solved to obtain explicit analytic expressions for the location and time variables. Retrieval errors are not interpreted with conventional Geometrical Dilution of Precision (GDOP) arguments as discussed by Holmes and Reedy (1951), but with more recent inversion analyses considered by Twomey (1977). Measurement errors are propagated analytically so that the specific effect of these errors on the solution is clarified. The sensitivity of the solution on the number of antennas used, antenna network geometry, source position, and measurement differencing schemes are discussed in terms of the eigenvalues of the linear system

A Grobner Basis Solution for Lightning Ground Flash Fraction Retrieval

A Bayesian inversion method was previously introduced for retrieving the fraction of ground flashes in a set of flashes observed from a (low earth orbiting or geostationary) satellite lightning imager. The method employed a constrained mixed exponential distribution model to describe the lightning optical measurements. To obtain the optimum model parameters, a scalar function was minimized by a numerical method. In order to improve this optimization, we introduce a Grobner basis solution to obtain analytic representations of the model parameters that serve as a refined initialization scheme to the numerical optimization. Using the Grobner basis, we show that there are exactly 2 solutions involving the first 3 moments of the (exponentially distributed) data. When the mean of the ground flash optical characteristic (e.g., such as the Maximum Group Area, MGA) is larger than that for cloud flashes, then a unique solution can be obtained

SCATTERING BY AN OBSTACLE IN A HALF-SPACE BOUNDED BY A PENETRABLE INTERFACE.

SCATTERING BY AN OBSTACLE IN A HALF-SPACE BOUNDED BY A PENETRABLE INTERFACE

Modeling the Diffuse Cloud-Top Optical Emissions from Ground and Cloud Flashes

A number of studies have indicated that the diffuse cloud-top optical emissions from intra-cloud (IC) lightning are brighter than that from normal negative cloud-to-ground (CG) lightning, and hence would be easier to detect from a space-based sensor. The primary reason provided to substantiate this claim has been that the IC is at a higher altitude within the cloud and therefore is less obscured by the cloud multiple scattering medium. CGs at lower altitudes embedded deep within the cloud are more obscured, so CG detection is thought to be more difficult. However, other authors claim that because the CG source current (and hence luminosity) is typically substantially larger than IC currents, the greater CG source luminosity is large enough to overcome the effects of multiple scattering. These investigators suggest that the diffuse cloud top emissions from CGs are brighter than from ICs, and hence are easier to detect from space. Still other investigators claim that the detection efficiency of CGs and ICs is about the same because modern detector sensitivity is good enough to "see" either flash type no matter which produces a brighter cloud top emission. To better assess which of these opinions should be accepted, we introduce an extension of a Boltzmann lightning radiative transfer model previously developed. It considers characteristics of the cloud (geometry, dimensions, scattering properties) and specific lightning channel properties (length, geometry, location, current, optical wave front propagation speed/direction). As such, it represents the most detailed modeling effort to date. At least in the few cases studied thus far, it was found that IC flashes appear brighter at cloud top than the lower altitude negative ground flashes, but additional model runs are to be examined before finalizing our general conclusions