Engineering topics
Publications and source records attributed to Koshak, William J..
This viewgraph presentation describes the significant progress made in the flash-type discrimination algorithm development. The contents include: 1) Highlights of Progress for GLM-R3 Flash-Type discrimination Algorithm Development; 2) Maximum Group Area (MGA) Data; 3) Retrieval Errors from Simulations; and 4) Preliminary Global-scale Retrieval.
An algorithm has been developed to estimate the altitude distribution of one-meter lightning channel segments. The algorithm is required as part of a broader objective that involves improving the lightning NOx emission inventories of both regional air quality and global chemistry/climate models. The algorithm was tested and applied to VHF signals detected by the North Alabama Lightning Mapping Array (NALMA). The accuracy of the algorithm was characterized by comparing algorithm output to the plots of individual discharges whose lengths were computed by hand; VHF source amplitude thresholding and smoothing were applied to optimize results. Several thousands of lightning flashes within 120 km of the NALMA network centroid were gathered from all four seasons, and were analyzed by the algorithm. The mean, standard deviation, and median statistics were obtained for all the flashes, the ground flashes, and the cloud flashes. One-meter channel segment altitude distributions were also obtained for the different seasons.
A general theory for retrieving the fraction of ground flashes in N lightning observed by a satellite-based lightning imager is provided. An "exponential model" is applied as a physically reasonable constraint to describe the measured optical parameter distributions, and population statistics (i.e., mean, variance) are invoked to add additional constraints to the retrieval process. The retrieval itself is expressed in terms of a Bayesian inference, and the Maximum A Posteriori (MAP) solution is obtained. The approach is tested by performing simulated retrievals, and retrieval error statistics are provided. The ability to retrieve ground flash fraction has important benefits to the atmospheric chemistry community. For example, using the method to partition the existing satellite global lightning climatology into separate ground and cloud flash climatologies will improve estimates of lightning nitrogen oxides (NOx) production; this in turn will improve both regional air quality and global chemistry/climate model predictions.
The Geostationary Lightning Mapper (GLM) is a single channel, near-IR imager/optical transient event detector, used to detect, locate and measure total lightning activity over the full-disk. The next generation NOAA Geostationary Operational Environmental Satellite (GOES-R) series will carry a GLM that will provide continuous day and night observations of lightning. The mission objectives for the GLM are to: (1) Provide continuous, full-disk lightning measurements for storm warning and nowcasting, (2) Provide early warning of tornadic activity, and (2) Accumulate a long-term database to track decadal changes of lightning. The GLM owes its heritage to the NASA Lightning Imaging Sensor (1997- present) and the Optical Transient Detector (1995-2000), which were developed for the Earth Observing System and have produced a combined 13 year data record of global lightning activity. GOES-R Risk Reduction Team and Algorithm Working Group Lightning Applications Team have begun to develop the Level 2 algorithms and applications. The science data will consist of lightning "events", "groups", and "flashes". The algorithm is being designed to be an efficient user of the computational resources. This may include parallelization of the code and the concept of sub-dividing the GLM FOV into regions to be processed in parallel. Proxy total lightning data from the NASA Lightning Imaging Sensor on the Tropical Rainfall Measuring Mission (TRMM) satellite and regional test beds (e.g., Lightning Mapping Arrays in North Alabama, Oklahoma, Central Florida, and the Washington DC Metropolitan area) are being used to develop the prelaunch algorithms and applications, and also improve our knowledge of thunderstorm initiation and evolution.
Lightning optical flash parameters (e.g., radiance, area, duration, number of optical groups, and number of optical events) derived from almost 5 yrs of Optical Transient Detector (OTD) data are compared with peak current and multiplicity observations derived from the US National Lightning Detection Networkm (NLDN). Despite the relatively low lightning geolocation accuracy afforded by OTD, a total of 48,870 NLDN cloud-to-ground (CG) flashes were correlated with OTD flashes, or about 10,000 CGs per year. The median values of the above OTD flash parameters for the 48,870 CGs were, respectively: 0.137 J/square meters/sr/micrometers, 313.7 square kilometers, 0.189 s, 4 optical groups per CG, and 8 optical events per CG. Invoking the multiplicity data, the median number of optical groups per stroke was 2.5, and the median number of optical events per stroke was 5.0. Median values of peak current for negative and positive CGs were -21.6 kA and 17.8 kA, respectively, and as expected, the negative CGs had a larger average multiplicity than the positive CGs. A statistical summary is provided for all CGs, for positive and negative CGs, and for CGs from different seasons. Standard two-distribution hypothesis tests were perfonned to intercompare the population means of the various lightning parameters. In particular, and to greater than the 99% confidence level, it was found that positive CGs are on average more radiant, of greater areal extent, and are longer lasting than negative CGs. Rankings from a complete set of hypothesis tests between CGs of different polarities and from different seasons are also provided. Most notably, wintertime positive CGs tend to be more radiant, of greater areal extent, and longer lasting than any other group of CGs (i.e., negative springtime CGs, positive summertime CGs, etc.).
The Optical Transient Detector (OTD) is a space-based instrument specifically designed to detect and locate lightning discharges as it orbits the Earth. This instrument is a scientific payload on the MicroLab-1 satellite that was launched into a low-earth, 70 deg. inclination orbit in April 1995. Given the orbital trajectory of the satellite, most regions of the earth are observed by the OTD instrument more than 400 times during a one year period, and the average duration of each observation is 2 minutes. The OTD instrument optically detects lightning flashes that occur within its 1300x1300 sq km field-of-view during both day and night conditions. A statistical examination of OTD lightning data reveals that nearly 1.4 billion flashes occur annually over the entire earth. This annual flash count translates to an average of 44 +/- 5 lightning flashes (intracloud and cloud-to-ground combined) occurring around the globe every second, which is well below the traditional estimate of 100 flashes per second that was derived in 1925 from world thunder-day records. The range of uncertainty for the OTD global totals represents primarily the uncertainty (and variability) in the flash detection efficiency of the instrument. The OTD measurements have been used to construct lightning climatology maps that demonstrate the geographical and seasonal distribution of lightning activity for the globe. An analysis of this annual lightning distribution confirms that lightning occurs mainly over land areas, with an average land:ocean ratio of 10:1. A dominant Northern Hemisphere summer peak occurs in the annual cycle, and evidence is found for a tropically-driven semiannual cycle.
The preliminary design of an optical/acoustical instrument is described for making highly accurate real-time determinations of the location of cloud-to-ground (CG) lightning. The instrument, named the Lightning Optical Camera And ThundEr (LOCATE) sensor, will also image the clear and cloud-obscured lightning channel produced from CGs and cloud flashes, and will record the transient optical waveforms produced from these discharges. The LOCATE sensor will consist of a full (360 degrees) field-of-view optical camera for obtaining CG channel image and azimuth, a sensitive thunder microphone for obtaining CG range, and a fast photodiode system for time-resolving the lightning optical waveform. The optical waveform data will be used to discriminate CGs from cloud flashes. Together, the optical azimuth and thunder range is used to locate CGs and it is anticipated that a network of LOCATE sensors would determine CG source location to well within 100 meters. All of this would be accomplished for a relatively inexpensive cost compared to present RF lightning location technologies, but of course the range detection is limited and will be quantified in the future. The LOCATE sensor technology would have practical applications for electric power utility companies, government (e.g. NASA Kennedy Space Center lightning safety and warning), golf resort lightning safety, telecommunications, and other industries.
Laboratory calibration and observed background radiance data are used to determine the effective sensitivities of the Optical Transient Detector and Lightning Imaging Sensor, as functions of local hour and pixel location within the instrument arrays. The effective LIS thresholds, expressed as radiances emitted normal to cloud top, are 4.0 plus or minus 0.7 and 7.6 plus or minus 3.3 micro J/ster/m (sup 2) for night and local noon; the OTD thresholds am 11.7 plus or minus 2.2 and 16.8 plus or minus 4.6 microJ/ster/m (sup 2). LIS and OTD minimum signal to noise ratios occur from 0800 to 1600 local time, and attain values of 10 plus or minus 2 and 20 plus or minus 3, respectively. False alarm rate due to instrument noise yields approximately 5 false triggers per month for LIS, and is negligible for OTD. Flash detection efficiency, based on prior optical pulse sensor measurements, is predicted to be 93 plus or minus 4% and 73 plus or minus 11% for LIS night and noon; 56 plus or minus 7% and 44 plus or minus 9% for OTD night and noon, corresponding to a 12 - 20% diurnal variability and LIS:OTD ratio of 1.7. Use of the weighted daily mean detection efficiency (i.e., not controlling for local hour) corresponds to a sigma = 8 - 9% uncertainty. These are likely overestimates of actual flash detection efficiency due to differences in pixel ground field-of-view across the instrument arrays, which are not accounted for in the validation optical pulse sensor data.
Each year, thousands of lightning electric field disturbances are recorded and archived by the ground-based field mill (FM) network at the NASA Kennedy Space Center (KSC) and USAF Eastern Range (ER). The FM network has a range of several tens of kilometers, and a digital accuracy of 4 V/m. It has provided years of continuous lightning warning surveillance to KSC-ER space vehicle launch operations, and has undergone one major hardware upgrade since its inception in the early 1970s. Additional KSC lightning warning data is derived from a multistation radio time-of-arrival system called Lightning Detection and Ranging (LDAR). This system provides the location and space-time mapping of individual lightning channels (for both cloud and ground flashes). Additional lightning information for the KSC region is available from the National Lightning Detection Network (NLDN) and a 5-station local magnetic direction finder network. In this study, all of the above mentioned data are used to ground-validate data derived from the Lightning Imaging Sensor (LIS) onboard the Tropical Rainfall Measuring Mission (TRMM). The FM network can be used to retrieve the charges deposited in a lightning flash, provided the flash is within a few kilometers of the FM Network. Although it is rare to obtain a TRMM overpass of thunderstorms hat occur this close to the FM network, seven such storms have been found and examined in this study. We compare the times and locations of LIS optical pulses with the spatial-temporal character of the FM, LDAR, and magnetic direction finder data. We also inter-compare LIS optical pulse amplitude data with FM-derived charge magnitudes, number of LDAR radio sources, and peak current values from magnetic direction finder data. Generally speaking, LIS lightning locations and times agree favorably with the KSC ground-based systems for most cases, but little correlation appears to exist between optical pulse amplitude and any of charge, # LDAR sources, peak current), owing possibly to the effects of source complexity and/or cloud multiple scattering.
The design, alignment, calibration, and field deployment of a solid-state lightning detector is described. The primary sensing component of the detector is a potassium dihydrogen phosphate electro-optic crystal that is attached in series to a flat-plate aluminum antenna; the antenna is exposed to the ambient thundercloud electric field. A semiconductor laser diode (lambda = 685 nm), polarizing optics, and the crystal are arranged in a Pockels cell configuration. Lightning-caused electric field changes are related to small changes in the transmission of laser light through the optical cell. Several hundred lightning electric field change excursions were recorded during five thunderstorms that occurred in the summer of 1998 at the NASA Marshall Space Flight Center in northern Alabama.
The Optical Transient Detector (OTD) is a space-based instrument specifically designed to detect and locate lightning discharges (intracloud and cloud-to-ground) as it orbits the Earth. A statistical examination of OTD lightning data reveals that nearly 1.2 billion flashes occurred over the entire earth during the one year period from September 1995 through August 1996. This translates to an average of 37 lightning flashes occurring around the globe every second, which is well below the traditional estimate of 100 flashes per second. An average of 75% of the global lightning activity during the year occurs between 30' S and 30' N. An analysis of the annual lightning distribution reveals that an average of 82% of the lightning flashes occur over the continents and 18% over the oceans, which translates to an average land-ocean flash density ratio of nearly 11.
Since April 1995, lightning activity around the globe has been monitored with the Optical Transient Detector (OTD). The OTD observations acquired during the one year period from September 1995 through August 1996 have been used to statistically determine the number of flashes that occur over the Earth during each hour of the diurnal cycle, expressed both as a function of local time and universal time. The globally averaged local [il,htnina activity displays a peak in late afternoon (1500-1800 local time) and a minimum in the morning hours (0600- 1000 local time) consistent with convection associated with diurnal heating. No diurnal variation is found for oceanic storms. The diurnal lightning distribution (universal time) for the globe displays a variation of about 35% about its mean as compared to the Carnegie curve which has a variation of only 15% above and below the mean.
The Lightning Imaging Sensor (LIS) is a NASA Earth Observing System (EOS) instrument on the Tropical Rainfall Measuring Mission (TRMM) platform designed to acquire and investigate the distribution and variability of total lightning (i.e., cloud-to-ground and intracloud) between q35' in latitude. Since lightning is one of the responses of the atmosphere to thermodynamic and dynamic forcing, the LIS data is being used to detect deep convection without land-ocean bias, estimate the precipitation mass in the mixed phased region of thunderclouds, and differentiate storms with strong updrafts from those with weak vertical motion.
Changes in the surface electric field are frequently used to infer the locations and magnitudes of lightning-caused changes in thundercloud charge distributions. The traditional procedure is to assume that the charges that are effectively deposited by the flash can be modeled either as a single point charge (the Q model) or a point dipole (the P model). The Q model has four unknown parameters and provides a good description of many cloud-to-ground (CG) flashes. The P model has six unknown parameters and describes many intracloud (IC) discharges. In this paper we introduce a new analysis method that assumes that the change in the cloud charge can be described by a truncated multipole expansion, i.e., there are both monopole and dipole terms in the unknown source distribution, and both terms are applied simultaneously. This method can be used to analyze CG flashes that are accompanied by large changes in the cloud dipole moment and complex IC discharges. If there is enough information content in the measurements, the model can also be generalized to include quadrupole and higher order terms. The parameters of the charge moments are determined using a dme-dimensional grid search in combination with a linear inversion, and because of this, local minima in the error function and the associated solution ambiguities are avoided. The multipole method has been tested on computer-simulated sources and on natural lightning at the NASA Kennedy Space Center and U.S. Air Force Eastern Range.
The current National Lightning Detection Network (NLDN) retrieval algorithm is based on a chi sq minimization. This numerical approach finds the lightning ground strike location on an ellipsoidal Earth that optimally agrees with multiple station time of arrival (TOA) measurements and, if available, magnetic bearing data. In addition, an analytic solution for determining lightning ground strike locations on a spherical Earth using only TOA data has been recently introduced. In the current work, a quasi-analytic approach is suggested for determining lightning ground strike locations on an oblate spheroidal Earth by perturbing the spherical model results proposed by Koshak. Latitude, longitude, time, and the associated perturbed quantities are related through terms in a Taylor series. The correction terms may be considered collectively as a vector which may be calculated by an overconstrained inversion. Expressions for the derivatives contained in the Taylor series are given in terms of elliptic integrals.
The chance of having the TRMM satellite pass over east central Florida when there is lightning over the NASA Kennedy Space Center (KSC) and USAF Eastern Range (ER) is small; however, such a condition did occur on September 21, 1998 (Day 264). Starting at about 20:40 GMT, the Lightning Imaging Sensor (LIS) reported 5 flashes during a 90 second interval that the KSC-ER was within the sensor field of view. Ground-based instrumentation, the Lightning Detection and Ranging (LDAR) system and a network of electric field mills (FM), detected 6 flashes in the same interval. In this paper, we will compare the times and locations of the optical pulses that were detected by LIS with the times and locations of RF sources (LDAR) and the charges that were deposited by the flash (FM network). We will show that LIS responded to all flashes that the LDAR and FM network detected; however, two discharges that were separated by less than 1 second in time and by about 10 km in space were grouped as one flash by the LIS data processing algorithm. In spite of the fact that all flashes occurred near the edge of the LIS field of view, the locations of the LIS events were consistent with both the LDAR and FM locations (the latter are usually within 1-2 kilometers of each other and often are co-located). Two of the 5 flashes reported by LIS were shifted north by about 8 km from the corresponding LDAR and FM locations. The LIS flash times tended to be after the first LDAR pulse was detected and before the last, and the integrated light signal (per LIS event) was surprisingly constant over the 5 flashes that were detected by LIS. In the future, we plan to study more correlated events and will try to determine whether and how the LIS light signal is related to the charge transfer in the flash and/or the number and spatial extent of RF sources.
The Optical Transient Detector (OTD) is a space-based instrument specifically designed to detect and locate lightning discharges (intracloud and cloud-to-ground) as it orbits the Earth. A statistical examination of OTD lightning data reveals that nearly 1.2 billion flashes occurred over the entire earth during the one year period from September 1995 through August 1996. This translates to an average of 37 lightning flashes occurring around the globe every second, which is well below the traditional estimate of 100 flashes per second. An average of 75% of the global lightning activity during the year occurs between 30 deg S and 30 deg N. An analysis of the annual lightning distribution reveals that an average of 82% of the lightning flashes occur over the continents and 18% over the oceans, which translates to an average land-ocean flash density ratio of nearly 11.
A linear algebraic solution is provided for the problem of retrieving the location and time of occurrence of lightning ground strikes from an Advanced Lightning Direction Finder (ALDF) network. The ALDF network measures field strength, magnetic bearing and arrival time of lightning radio emissions. Solutions for the plane (i.e., no Earth curvature) are provided that implement all of tile measurements mentioned above. Tests of the retrieval method are provided using computer-simulated data sets. We also introduce a quadratic planar solution that is useful when only three arrival time measurements are available. The algebra of the quadratic root results are examined in detail to clarify what portions of the analysis region lead to fundamental ambiguities in source location. Complex root results are shown to be associated with the presence of measurement errors when the lightning source lies near an outer sensor baseline of the ALDF network. In the absence of measurement errors, quadratic root degeneracy (no source location ambiguity) is shown to exist exactly on the outer sensor baselines for arbitrary non-collinear network geometries. The accuracy of the quadratic planar method is tested with computer generated data sets. The results are generally better than those obtained from the three station linear planar method when bearing errors are about 2 deg. We also note some of the advantages and disadvantages of these methods over the nonlinear method of chi(sup 2) minimization employed by the National Lightning Detection Network (NLDN) and discussed in Cummins et al.(1993, 1995, 1998).