Tuesday Feb. 22, 2011
click here to download today's notes in a more printer friendly format

We came back to a "problem" mentioned at the end of class last Thursday.  Most thunderstorms have the dipolar (or tripolar) charge distribution as shown below at left.  The majority of cloud to ground (CG) discharges from these types of storms carry negative charge to ground.  Just a few percent of the CG discharges are positively charged (we'll see though that these positive CG discharges sometimes have very large peak currents).

A much larger percentage (50% or more) of positive CG discharges have been observed coming from some Central Plains storms in recent field experiments.  This would suggest the clouds might have had an inverted charge distribution like that shown below at right.  And one might wonder whether a different cloud electrification mechanism is at work (different from the Reynolds, Brook, Gourley process discussed last Thursday).



Other than the high incidence of positive CG discharge, E field soundings in thunderstorms have provided evidence of the inverted charge structure.  This is discussed in a fairly recent paper that you can find here.  The two figures below (on a class handout) illustrate E field soundings in "normal" clouds and in clouds with inverted charge distributions.



Follow up studies seem to indicate that the graupel-ice crystal collision, non-inductive electrication mechanism that we have been studying is cable of explaining clouds with inverted charge distributions.  This is shown on the figure below


In the center of the figure is the graph that we examined in some detail last Thursday.  It shows the polarity of charge acquired by a graupel particle colliding with ice crystals in a laboratory simulation of a cloud environment containing supercooled water droplets.  The liquid water content on the vertical axis is a measure of the supercooled droplet concentration.

Liquid water content in a typical cloud would be found near the level of the lower dotted line on the graph.  The tripolar cloud charge distribution in a typical cloud is shown at left.  High in the cloud where the temperature is cold, the graupel particle acquires negative charge and the ice crystals positive charge (region B on the graph).  The ice crystals are carried upward and form the main positive charge center.  The heavier graupel particle descends to form the main negative charge center.  Somewhat lower in the cloud where temperatures are warmer than TR (a reversal temperature) but still below freezing the polarity of the charging changes (region C on the graph).  The graupel ends up with positive charge and the ice crystals with negative charge.  The positively charged graupel form the lower positive charge centers.

It seems that the unusual Central Plains storms have very high liquid water contents.  If we look at the level of the upper dotted line on the charging graph in the center of the figure we see that the graupel particle always ends up with positive charge in this high LWC environment.  There is no charge reversal temperature.  The Reynolds, Brook, Gourley mechanism can explain clouds with inverted charge distributions.

The figure below was included on a class handout just to indicate that some very complex cloud charge distributions are found in the stratiform regions of Mesoscale Convective System (MCS) clouds (we didn't and won't spend anytime discussing the structure and formation of MCS clouds)





The main topic of the day was looking a how we can use E field measurements at the ground to learn about the location and amounts of charge involved in lightning discharges.  Most of the figures below were on handouts distributed in class.

The data and results that we will be discussing come largely from experiments conducted at the Kennedy Space Center (KSC).  A large network of electric field mills has been installed and is operating at KSC to identify and warn of thunderstorm and lightning hazards.  The pre-1995 configuration of the field mill network is shown in the figure below (the network was upgraded in 1995).

The network that was operating up to 1995 has 0.1 second time resolution.  That is fast enough to resolve a lightning flash, but not the individual return strokes that make up cloud to ground flashes.  The dynamic range was -15 kV/m to +15 kV/m and E field signals were digitized with 30 V/m accuracy.  The overall accuracy of an individual field mill was about 10%.


Examples of fields recorded during 6 small storms.  Storms in this category lasted from 35 up to 85 minutes.  Individual storms produced 16 to 82 discharges.  The maximum flashing rate during a small storm is 1 flash/minute.

A single example of a large storm E field record is shown below

Storms like this are often broken into initial, very active, and end-of-storm-oscillation portions.  Note the very slow but large amplitude oscillations in the final EOSO portion. 

Large storms in the Livingston and Krider study had durations ranging from 75 to 265 minutes, produced 515 to 1212 discharges, and maintained flashing rates of 5 to 10 flashes per minute for 50 to 90 minutes.  Number of flashes per 5 minute interval in a reprsentative large storm are shown in the histogram below.





Fields at the ground are usually lower than fields measured just a few hundred meters above the ground.  Space charge from corona discharge at the ground limits the amplitudes of fields at the ground and affects the shape of the field recovery between flashes.  You can get an accurate measurement of the field change at the ground (it occurs so quickly that space charge can't be created or move quickly enough to affect the field change value)



Measured field change values as a function of distance.  Note that small field change values are sometimes observed very close to a discharge.  This observation will come up later in today's notes.

Now we will begin to look at what you can do with simultaneous measurements of field change values at multiple field mill sites.


We'll consider the simplest model of a cloud-to-ground discharge and assume that the charge neutralized by the flash comes from a uniform sphere of charge in the cloud.  You can treat this as a single point charge (Point 1).  We will try to use the delta E measurements to determine the location (x, y, and z coordinates) and magntitude of the neutralized charge (delta Q). 

We show the location of one of the field mill sites in the figure.  At that location we will have a measurement of the field change.  We can also calculate the field change that the charge neutralized at Point 1 would produce; the expression for this calculated field change is shown at Point 2.

We have pairs of measured and calculated for all or most of the field mill sites.  The idea then is to adjust x, y, z, and delta Q until the chi-squared function at Point 4 is minimized.



Some of the results obtained using the single charge model.  Fig. 1 shows the magnitude of the charge neturalized.  Fig. 2 shows that this charge was largely found between 6 and 9 or 10 km altitude, in the -10 to -35 C temperature range.  Fig. 3 shows that 75% of cloud-to-ground discharges strike the ground at a distance, D, of 5 km (3 miles) or less from the center of the charge neutralized by the flash.  95% of the discharges strike the ground within 8 km (5 miles) of the charge center.  The distance D is shown in the figure below.



Instead of getting larger as you get closer to the storms center, field charges at the ground are often either very small or of opposite polarity.  This suggested that a lower volume of positive charge might be involved in cloud-to-ground discharges.  This led to the development of a 2 charge model, illustrated below.



This is the most general form of the two charge model.  The two charges can have different magnitudes and can have completely different locations in space ( the point dipole model assumes the two charges are of equal amplitude but opposite polarity, another model assumes the two charges are aligned vertically.  The calculated field change expression now contains 8 unknowns.

The next two figures shows some of the results obtained with this arbitrary 2 charge model.


The circles indicate charge neutralized during lightning discharges.  Cross hatched circles are positive, open circle negative charge.  The radius of the circle gives you an idea of the amount of charge neutralized.

Both intracloud and cloud to ground discharges are plotted.  An example of each is shown.  Intracloud discharges involve positive charge in the main positive charge center located in the upper portion of the thundercloud and negative charge in the main negative charge center in the middle of the cloud.  Cloud-to-ground discharges seem to almost always involve negative charge in the middle of the cloud and positive charge in one of the lower positive charge centers.

The next figure is from a little more active thunderstorm cell.


Note the tendency for the altitudes of the positive charges neutralized in intracloud discharges to be found at higher altitude toward the middle and perhaps the end of the storm.  Also the amounts of charge neutralized in discharges appears to increase with time (larger diameter circles toward the end of the storm).


The chi-squared procedure for determining charge center locations and charge magnitudes using multi-station field change measurements has also be used with Slow E field antenna systems.  Slow E records have faster time resolution and can faithfully record the field changes produced by the separate return strokes in a cloud to ground flash.

A description of an 8-station slow E antenna and an example of a CG flash slow E field record  (containing 5 or 6 return strokes) is shown in the figure below (on a class handout).


The network was installed near Socorro New Mexico and operated by the New Mexico Institute of Mining and Technology (I was an undergraduate student there).  The figure below shows locations of the charges neutralized by the separate return stroke in 4 cloud-to-ground flashes.  The simple point charge model was assumed.


Notice that each of the return strokes "taps" charge from a slightly different location in the cloud.