Tuesday, April 18, 2017

We'll begin today's class by finishing up the short section on the validation of return stroke peak current estimates made by the NLDN.  You'll find that material at the end of the Tuesday, April  11 notes.

In the remainder of the class (and perhaps the first part of Thursday's class) we will examine how time of arrival (TOA) and interferometric systems can be used to locate the sources of VHF radio emissions from lightning.  The wavelength of VHF radiation is much shorter than the VLF signals the signals propagate in a straight line from the source to the receiving antenna.  That means the source of the radiation can be located in three dimensions.  For discharge processes that produce short duration and separate pulses of radiation we are able to observe and follow channel development as it occurs in the cloud.



Most lightning discharge processes are strong emitters of VHF radiation.  This is evident on the figure above (from Hayenga and Warwick (1981)) which shows VHF emissions and a slow E field record of a multi-stroke cloud-to-ground discharge.  When we examine the VHF radiation on a faster time scale however we see that the VHF signals can broadly divided into two groups.  First sequences of isolated pulses of a few microseconds duration and bursts of "quasi-continuous" emission that last for 10s of microseconds or more.
  Examples of both types of signals are shown in the figure below adapted from Proctor et al., (1988).


On Row (a) are examples of the isolated pulses that are approximately 1
μs wide.  Signals of this type are produced by discharge processes that propagate into unionized air and are best located by TOA systems.  Examples of quasi-continuous radiation are shown Row (b) (Proctor refers to these types of signals as "Q-noise").  The two bursts in the figure lasted 30 μs and 45 μs   This type of radiation seems to be produced by return strokes and recoil streamers, i.e. a discharge process that travels back along an existing channel.  Interferometry is best able to locate the sources of this radiation.


Probably the first 3-D TOA location system used for lightning research was D.E. Proctor's South African network.  It is shown below and is described in more detail in Proctor (1971).  This is a long baseline system where the distance between receivers (a few 10s of km in this case) is much greater than the wavelength of the VHF radiation. 

In Proctor's sytem, four 253 MHz receivers (later changed to 355 MHz) were positioned at the ends of two nearly perpendicular baselines.  A fifth receiver was positioned in the center at the intersection of the two baselines.  Signals received at the four outer stations were transmitted to the center station by microwave link.  Because the time of propagation from the outer stations to the center station could be accurately determined, the time of the arrival of signals at the outer stations could be determined very precisely.




The figure below shows, in crossection, how a constant TOA difference between two receiving stations traces out a hyperbola.  A signal originating anywhere on the green line would have the same TOA difference between stations S0 and S2.  The TOA difference between stations S0 and S1 trace out a second hyperbola.  The intersection of the two hyperbolas is the source of the RF radiation.



Though because we are now working in
3 dimensions, each of the hyperbolae should be rotated around the baseline connecting the stations.


So instead of a single point of intersection we actually have an arc, shown in red above.  To locate the RF emissions source in space a second baseline of stations is needed.

An example of RF data from the 5 receivers is shown below (from Proctor (1983)).  The records are approximately 60 microseconds long. 


Waveform digitizers were not fast enough at that time to be able to record data like this.  So the RF signals from the 5 stations were initially recorded on a cathode ray tube and photographed (a specially designed "laser optical recorder" was later used to record the data directly onto film).  In the example above, the data recorded on film have been digitized and displayed on a computer for analysis.

In order to locate an emission source in space you must first identify the correlated RF impulses on the 5 records.  The records can be complex with signals from multiple sources arriving at the different receivers in different order. 
This analysis was originally done by hand and could be quite tedious.  The following quotes are from Proctor (1981) (I added the underlines).

"After about 2 years experience a good data reader can unscramble the pulses emitted by no less than 4 simultaneous flashes, or widely spaced branches of the same flash with a considerable degree of confidence."

"Even those who enjoy reading reading records that can be deciphered easily found that reading the more complicated variety was a mild from of torture, and that it took about one man-month's effort to locate 100 sources correctly."

On average it was possible to locate the source of one pulse in every 70 microseconds of record
.


We'll look at just one example of many interesting results obtained with the South African VHF TOA network.  Proctor (1991) determined where lightning discharges began by finding the centroid of the first 6 to 10 VHF pulse locations emitted in a flash.  The left figure below shows origin heights (above ground level) for 773 cloud-to-ground and intracloud flashes in 13 thunderstorms.  Origin heights for 214 cloud-to-ground flashes are shown in the figure at right.  The network was 1.43 km above sea level.



For the lower altitude peak in the left figure, median heights for individual storms ranged from 4.4 km to 5.7 km above mean sea level (amsl) which corresponded to temperatures ranging from +1o C to -8.5o C.  Median heights for the 2nd higher altitude peak ranged from 7.5 km to 9.7 km amsl where temperatures ranged from -21.1o C to -33o C.


We can refer back to a figure from an earlier lecture (reproduced below) showing the locations of charge centers involved in cloud-to-ground and intracloud discharges (the charge center locations where derived from E field change measurements made at multiple field mill sites at the Kennedy Space Center).


We can see that the higher altitude peak on the origin height histogram corresponds roughly with the high E field region between the main positive charge center and the main negative charge centers.  The lower altitude peak is positioned roughly in another high field region between the main negative charge center and one of the lower positive charge centers.


Next we'll briefly discuss another long baseline VHF TOA network, the Lightning Detection and Ranging (LDAR) system, that began operating at the NASA Kennedy Space Center in the mid 1970s.  The network has undergone several changes and upgrades since it was first installed.  The figure below shows the layout of the system that was being used in the 1990s.  It consisted of a central station and 6 surrounding stations positioned on roughly a 20 km diameter circle. 




The VHF receivers operate at 66 MHz (6 MHz bandwidth) and data from the remote stations were telemetered back to the center station by microwave link where they are digitized.  A signal crossing threshold at the center receiver triggers a roughly 100 μs; long recording at each of the 7 stations.  Correlated pulses on the separate records and differences in the times of arrival are found using pattern recognition and cross correlation techniques.  TOA data from the center station and 3 of the outer stations (stations 0, 1,3, and 5 for example) are used to determine an RF source location.  A second location is determined using the center data and the other 3 remote locations (stations 0, 2, 4 and 6).  If the two locations agree the location is accepted.


Here is an example of RF data from the LDAR network (from Rustan et al., (1980)).  Just over 100 μs of data from the center station and 3 of the outer stations are shown.  Four impulses have been marked on the top trace.  See if you can find the corresponding impulses on the 3 remaining records.  Click here to see what the analysis procedure, now done by computer using pattern recognition and cross correlation techniques, deemed to be the correct associations.

Here is an example of VHF source locations provided by the LDAR network (from the reference in the paragraph above).



The figure shows the beginning of a 3-stroke CG flash that struck a 150 m tall weather tower at the Kennedy Space Center.  This event occurred during the first year of the
Thunderstorm Research International Program and was observed by a variety of types of instruments (see Uman et al., (1978) for more details).

VHF noise began 4.9 ms before the 1st return stroke. 
The first 2.2 ms was considered preliminary breakdown activity and was followed by the stepped leader which lasted 2.7 ms.  The left figure above shows 422 VHF sources located during the first 4.7 ms of the flash.    Clear visualization of the data becomes a concern with the large amounts of data that are acquired with these systems.  The right figure plots locations at 94 μs intervals during the preliminary breakdown (between Pts. A & B) and the stepped leader (below Pt. B).  The channel has been highlighted in red.  Q1 is the location of the charge transported to ground during the first return stroke in this flash (determined from field change records using the field mill network).  Four additional plots like this tracing out activity during the remainder of the discharge can be found in the Rustan et al. (1980) report.


We'll next consider the most recent and most advanced VHF TOA system, the Lightning Mapping Array (LMA) developed by researchers at the New Mexico Institute of Mining and Technology (New Mexico Tech).  The system detects 60-66 MHz radiation (television channel 3).

The system was originally used in Oklahoma in June 1998 (a map of sensor locations can be found here and results from the Oklahoma field experiment can be found here or here) and then in New Mexico in August and September.  The New Mexico network, as described by Rison et. al (1999), consisted of 10 stations that were deployed over an area about 60 km in diameter near Langmuir Laboratory.  Quoting from the Rison et. al article: "The time and magnitude of the peak radiation is recorded during every 100 μs time interval that the RF power exceeds a noise threshold."  The time of the peak signal is recorded with 50 ns resolution.  TOA information from events strong enough to be detected at six or more stations are located in space and time.  This TOA information is determined at each of the remote stations and the information is sent back to a central location (rather than transmitting the raw VHF data back to a base station).

It is a little difficult to keep up with the rapid development and results coming from the LMA.  As best I can tell (2012 source) permanent and operational LMAs can be found in or at the NASA Marshall Space Flight Center in Alabama (see also this link), the National Weather Service in Washington DC, Oklahoma, Texas A&M University near HoustonWest Texas (Texas Tech University),   Displays of real time data are available online at many of these sites.  Here's a link to a useful site with a variety of lightning related resources.

Some photographs of an early version of the system being demonstrated in the Washington DC area are shown below (source of these images)



VHF antenna and open electronics enclosure showing the electronics (a more detailed photograph of the electronics can be found here)



VHF antenna installation at the Sterling VA location
.



Electronics enclosure
.


Here are some animations of lightning activity located with the Oklahoma LMA.  The links and explanation of the format of the data display below have been copied from http://lightning.nmt.edu/nmt_lms/ok98.html.

Here is a link to images of an intracloud and a cloud-to-ground discharge that I will probably show in class in order to get familiar with the display format before showing any of the animations below.


Here are some MPG movies from the analyses of two lightning discharges in Oklahoma. The plots have five windows in them:

The top t vs. z plot shows all the points located in the given time interval. The time is given in seconds, and the altitude z of the radiation source is given in kilometers. The number of points can be limited if desired. For example, two distinct discharges at different locations could occur at the same time, so the operator may choose to look at only one of the dishcarges in detail. Or the operator may narrow in on a particular event in time which he does not want obsured by other data points. The lower four plots show data points from the top plot which have been selected for such reasons. The lower t vs. z plot shows the time development of the altitude of the selected points. The x vs. y plot shows a plan view of the lightning discharge, where x is the distance east or west of the center of the array, and y is the distance north or south of the center. All distances are given in kilometers. The small squares in the x vs. y plot show the locations of the LMA stations. The x vs. z and y vs. z plots show the projections of the points on the xz and yz planes.

Discharge showing distinct charge layers at different elevations. Note the development of the dendritic structure of the discharge.
11 June 1998 06:33:20 UTC (645 KB Animated GIF)

Cloud-to-ground discharge followed by intercloud discharge. Note the well-defined leaders in the CG. The triangles are locations and times of CG strokes as determined by the National Lightning Detection Network. Note that, in the CG, positive streamers propagate in the negative charge region, radially away from the area of initial breakdown. The subsequent IC develops over the top of of the CG.
11 June 1998 06:19:39 UTC (1.5 MB Animated GIF)

Discharge near the end of a storm which appears to be inverted. There are two distinct charge layers, but streamers appear to originate in the upper layer and propogate to the lower layer. This is inverted from most intercloud discharges in which negative streamers originate in the lower negative charge layer and propogate into the upper positive charge layer. Also note how the lower charge region decreases in altitude to the east.
20 June 1998 03:43:45 UTC (1.8 MB Animated GIF)


Here's another link to some Oklahoma research results:
Lightning locations in a tornadic thunderstorm

Finally, be sure to look at this very cool animation of a lightning flash over New Mexico near Langmuir Lab that is on YouTube.

More recently a LMA was used at the ICLRT facility in northern Florida to study triggered lightning.  The figure below at left is a 3-D visualization of the LMA locations from the first 200 ms of the upward propagating leader/initial continuous current  (UPL/ICC) stage at the beginnings of triggered discharges.  The altitudes of the LMA sources are plotted together with the measured channel base current in the figure at right (both figures adapted from Hill et al. (2012)).







The UPL/ICC process in this flash lasted 945 ms and transported 225 C of charge to the ground.  Both are larger than more typical values of 300 ms and about 30 C.  The upward propagating leader that initiates triggered lightning is positively charged.  It has often been observed that positive polarity leaders do not radiate as strongly at VHF as leaders that are negatively charged.   Here VHF radiation was strong enough that impulses with currents as low as 10 A could be located.  Most of the branching occurred above cloud base and is not visible on video records made at the ground.  Speeds of the UPL and branches ranged from 2 to 12 x 104 m/s.  Somewhat surprisingly there is no indication on the ground level current record of when branching occurs above the ground.  The 3 precursor pulses (the left most point is actually two closely spaced points) shown on the figure at right are discharges that begin to move upward from the tip of the triggering wire but aren't able to continuing moving upward on their own.



Finally we will have a brief look at a different way of locating sources of lightning VHF emissions - interferometry.  It's not clear at this point whether we'll have time to cover this material in class on Tuesday.

T
he basic principle involved in interferometry is shown below.



A plane wave of narrow band (single frequency) radiation is approaching from the right in the same direction as a line connecting two antennas on the ground.  The antennas are a distance d apart.  Radiation will arrive at Antenna 2 first.  The radiation arriving at Antenna 1 travels an additional distance, l, before reaching the antenna.  The interferometer measures the phase difference, α, in the signals arriving at the two antennas. 

The elevation angle of the arriving signal can then be determined from the measured phase angle,


In order to have a unique solution for the elevation angle, the phase angle can't vary by more than 2π as elevation angle ranges from 0o to 180o



This puts an upper limit on the distance separating the two antennas


In the case where d = λ /2, the phase angle would be π for a signal approaching from the right at θ = 0o (the signal arrives at Antenna 2 before Antenna 1) and phase angle would be -π for a signal arriving from the left at θ = 180o (the signal arrives at Antenna 1 first).   I'm not sure whether a phase detector can distinguish between a phase angle of π and -π or whether it would just measure π in both cases. 

The figure below is a polar plot of the absolute value of the phase angle versus the elevation angle of the incoming signal (for
d = λ /2). 




The circular rings are phase angle varying from 0 at the center of the plot to π at the outer edge of the plot.  There won't be any phase angle difference for signals arriving from overhead (elevation angle, θ = 90o )  because identical features on the signal would reach both antennas at the same time.  This plot leaves the impression that there are two solutions for a single value of the phase angle.  For example, the phase angle is π/2 at Points A and B.  At Point A the signal is arriving at an elevation angle of 60o from the right.  This is shown below.  Point B is a signal arriving at the same elevation angle but from the left.  If the phase detector is able to distinguish between phase angles of +π and -π it would be able to determine whether the radiation was coming from the right or left. 

To illustrate the problem of multiple elevation angles for a given value of phase angle we draw a figure for a larger antenna
separation, d = 2
λ.




Points A, B, C, and D have phase angles of 4π. 2π. -2π and -4π.  These multiple solutions are often referred to as fringes.  At best, the phase detector wouldn't be able to distinguish between 2π and 4π (assuming it would be able to distinguish between positive and negative phase angles) so a measurement of phase angle will lead to at least 2 ambiguous solutions for the elevation angle.  There would be 4 possible elevation angles for a given phase angle measurement if the phase detector is unable to distinguish positive from negative values.  The four labelled points above are shown on a polar plot below.



Clearly it would seem like two closely spaced antennas would be best.  However, and we won't go into the details, the error in the elevation angle determination is proportional to 1/d.  So, while there won't be any ambiguities in the elevation angle determination when two antennas are closely spaced, the error could be large.  Two more widely spaced antennas would result in less error but there would be multiple elevation angles possible for a given measured phase angle difference.  What is generally done is to add a 3rd antenna to the baseline as shown below.


A distant antenna to reduce elevation angle error and two close antennas to resolve elevation angle ambiguities.

Up to this point we've been assuming that the direction of the incoming radiation is parallel to the baseline connecting the antennas (i.e. zero azimuth angle).  There is no reason for that to be the case.  As the azimuth angle moves from zero, the measured phase difference will begin to decrease.  The phase angle difference will become zero for a signal approaching from a direction that is perpendicular to the baseline connecting the antennas.  So to be able to determine the true direction angle to the emission source we're going to need a 2nd perpendicular baseline.  Antenna 1 above can be part of both baselines.  So we'll end up with something like shown below. 


An antenna array essentially identical to this was used by Rhodes et. al (1994).  Separation distance between the two inner pairs of antennas was λ/2.  4λ separation was used between the outer pairs of antennas. 



Richard et. al (1986) used 6 antennas arranged in a triangular pattern as shown above.  Spacing between the inner antennas was either
λ/2 or λSpacing between antennas in the outer triangle was 10λ (more than 250 ambiguous elevation angle values).


Full citations and links to references mentioned in this lecture.

C.O. Hayenga and J.W. Warwick, "Two-Dimensional Interferometric Positions of VHF Lightning Sources," J. Geophys. Res., 86, 7451-7462, 1981.

Hill, J.D., J. Pilkey, M.A. Uman, D.M. Jordan, W. Rison, and P.R. Krehbiel, "Geometrical and electrical characteristics of the initial stage in Florida triggered lightning," Geophys. Res. Lett., 39, 2012.

Proctor, D.E., "A Hyperbolic System for Obtaining VHF Radio Pictures of Lightning," J. Geophys. Res., 76, 1478-1489, 1971.


Proctor, D.E., "VHF Radio Pictures of Cloud Flashes," J. Geophys. Res., 86, 4041-4071, 1981.

Proctor, D.E., "Lightning and Precipitation in a Small Multicellular Thunderstorm," J. Geophys. Res., 5421-5440, 1983.

Proctor, D.E., R. Uytenbogaardt, B.M. Meredith, "VHF Radio Pictures of Lightning Flashes to Ground," J. Geophys. Res., 93, 12683-12727, 1988.

Proctor, D.E., "Regions Where Lightning Flashes Begin," J. Geophys. Res., 96, 5099-5112, 1991.

Rhodes et al., 1994

Richard et al., 1986.

Rison, W. R.J. Thomas, P.R. Krehbiel, T. Hamlin, and J. Harlin, "A GPS-based Three Dimensional Lightning Mapping System: Initial Observations in Central New Mexico," Geophys. Res. Lett., 26, 3573-3576, 1999.

Rustan, P.L., M.A. Uman, D.G. Childers, and W.H. Beasley, "Lightning Source Locations From VHF Radiation Data for a Flash at Kennedy Space Center," J. Geophys. Res., 85, 4893-4903, 1980.

Uman, M.A., W.H. Beasley, J.A. Tiller, Y-T. Lin, E.P. Krider, C.D. Weidman, P.R. Krehbiel, M. Brook, A.A. Few, J.L. Bohannon, C.L. Lennon, H.A. Poehler, W. Jafferis, J.R. Gulick, J.R. Nicholson, "An Unusual Lightning Flash at Kennedy Space Center," Science, 201, 9-16, 1978.