The Energy, Momentum, and Peak Power Radiated by Negative Lightning Return Strokes

: Electromagnetic radiation ﬁelds generated by return strokes transport both energy and momentum from the return stroke to outer space. The momentum transported by the radiation ﬁeld has only a vertical or z component due to azimuthal symmetry (cylindrical symmetry) associated with a vertical return stroke. In this paper, the energy, momentum, and peak power radiated by return strokes as a function of the return stroke current, return stroke speed, and the zero-crossing time of the radiation ﬁelds are studied. The results obtained by numerical simulations for the energy, vertical momentum, and the peak power radiated by lightning return strokes (all parameters normalized by dividing them by the square of the radiation ﬁeld peak at 100 km) are the following: A typical ﬁrst return stroke generating a radiation ﬁeld having a 50 µ s zero-crossing time will dissipate ﬁeld normalized energy of about (1.7–2.5) × 10 3 J / (V / m) 2 and ﬁeld-normalized vertical momentum of approximately (2.3–3.1) × 10 − 6 Kg m / s / (V / m) 2 . A radiation ﬁeld with a zero-crossing time of 70 µ s will dissipate about (2.6–3.4) × 10 3 J / (V / m) 2 in ﬁeld-normalized energy and (3.2–4.3) × 10 − 6 Kg m / s / (V / m) 2 in ﬁeld-normalized vertical momentum. The results show that, for a given peak radiation ﬁeld, the radiated energy and momentum increase with increasing zero-crossing time of the radiation ﬁeld. The normalized peak power generated by a ﬁrst return stroke radiation ﬁeld is about 1.2 × 10 8 W / (V / m) 2 and the peak power is generated within about 5–6 µ s from the beginning of the return stroke. Conversely, a typical subsequent return stroke generating a radiation ﬁeld having a 40 µ s zero-crossing time will dissipate ﬁeld-normalized energy of about (6–9) × 10 2 J / (V / m) 2 and ﬁeld-normalized vertical momentum of approximately (7.5–11) × 10 − 7 Kg m / s / (V / m) 2 . The ﬁeld-normalized peak power generated by a subsequent return stroke radiation ﬁeld is about 1.26 × 10 8 W / (V / m) 2 and the peak power is generated within about 0.7–0.8 µ s from the beginning of the return stroke. In addition to these parameters, the possible upper bounds for the energy and momentum radiated by return strokes are also presented.


Introduction
Electromagnetic fields generated by lightning play an essential role in lightning protection studies since they generate unwanted currents and voltages in electrical and electronic systems, causing them damage and disruptions [1].They can be utilized as a vehicle to study the terrain features by analyzing how these fields are modified when they travel along grounds of different configurations and electrical properties [2][3][4].These fields can also be used in remote sensing of the spatial and temporal distribution of lightning currents by appealing to different types of return stroke models [5].Moreover, upper atmospheric physicists also need information concerning electromagnetic fields generated by lightning because these fields interact with the upper atmosphere generating electrical discharge activity such as sprites and elves [6].Electromagnetic fields from lightning can also modify the temperature of the electrons in the ionosphere and therefore, could change the electromagnetic propagation characteristics of the ionosphere [7].The energy transported by the electromagnetic fields plays a significant role in many of these interactions.Moreover, in physics-based return stroke models, the knowledge concerning the fraction of energy removed from the lightning channel by electromagnetic radiation is an essential parameter in evaluating the energy balance of the lightning return stroke [5].
Taylor and Jean [8] conducted the first study on the electromagnetic energy radiated by lightning flashes.Using typical features of atmospherics (i.e., distant radiation fields) generated by return strokes, they estimated that a return stroke would dissipate about 3 × 10 4 J on average as electromagnetic radiation.Taylor [9] improved the results of the study by using a large sample of atmospherics and came up with 7 × 10 4 J as a typical value for the radiated energy.Recently, armed with more recent data pertinent to electromagnetic fields from return strokes, Krider and Guo [10] estimated that the peak power radiated by typical first and subsequent return strokes to be about (1-2) × 10 10 W and (3-5) × 10 9 W, respectively.In determining the power and energy dissipated by return strokes utilizing the measured electromagnetic fields, both Taylor [9] and Krider and Guo [10] represented the return stroke as an infinitesimal dipole at ground level with a suitable current moment that reproduced the measured electromagnetic field at ground level.Using this current moment in the dipole, they estimated the Poynting vector at any given point in space.The energy and peak power are estimated by integrating the resulting Poynting vector over the surface of a hemisphere with the center at the strike point.Now, replacing the return stroke as a dipole at ground level will simplify the evaluation of the radiated power and energy from the lightning return stroke.Still, it will introduce errors because the return stroke is an extended source through which current pulses are propagating with speeds close to the speed of light.For example, in the dipole approximation, the temporal variation of the electromagnetic radiation fields is independent of the elevation angle of the point of observation with respect to the vertical, while in an extended source, the duration of the electromagnetic field at different elevation angles depends on the propagation speed of the return stroke; this could influence the calculated energy.
Krider [11] improved the energy calculation by utilizing electromagnetic fields pertinent to the transmission line model [12] to derive an expression for the Poynting vector.From this expression, he estimated the angular distribution of the electromagnetic power radiated by lightning.However, the study has not assessed the total energy radiated out by a return stroke.It is important to point out here, that in the transmission line model, the return stroke current is assumed to propagate with constant speed without attenuation.In reality, the current attenuates, and the speed decreases as the return stroke front forge along the return stroke channel.Both these effects could affect the power and energy radiated by return strokes.
The energy dissipated by a return stroke radiation field depends on its initial peak amplitude, which is related to the return stroke peak current, and on its duration.The duration of the radiation field varies from one return stroke to another, and the zero-crossing time of a radiation field is a good measure of this duration.The zero-crossing time of first return stroke radiation fields may vary from 50 µs in the temperate regions to about 90 µs in the tropics [13].In the case of subsequent return strokes, it may vary from about 20 µs to about 60 µs.The goal of this work is to examine how the energy, momentum, and peak power radiated by return stroke radiation fields vary with these parameters.Neither the momentum transported by the lightning electromagnetic fields nor the dependence of radiated energy and momentum on the zero-crossing time of the radiation field has ever been studied previously in the literature.But, this information is valuable in understanding the coupling of the lightning generated electromagnetic fields to the upper atmosphere.This paper is an attempt to fill this gap in the literature.
It is vital to point out here that even for a given initial peak and the zero-crossing time of the radiation field, the radiated energy and momentum may depend on the details of the temporal variation of the return stroke radiation field.The temporal variation of the radiation field may also vary from one return stroke to another.To take this variation into account, in our analysis we will use two types of return stroke models that produce different wave shapes for the radiation fields.The two models provide return stroke radiation fields that differ from each other, and we believe that use of these two models will allow us to obtain more information on how the energy and momentum radiated by these fields vary with the shape of the radiation fields.
In this paper, we will use the term range-normalized electric radiation field to represent the peak value of the electric radiation field corresponding to 100 km distance.In describing the results, we will also use the terms field-normalized energy, field-normalized momentum, and field-normalized power.These quantities are generated by dividing the energy, momentum, and power by the square of the electric radiation field peak at 100 km.
Preliminary results of this paper were presented at the 34th International Conference on Lightning Protection (ICLP) held in Rzeszow, Poland, in 2018 [14].

Lightning Return Stroke Model Used in the Study
It is necessary to appeal to a return stroke model that specifies the spatial and temporal variation of the return stroke current to calculate the electromagnetic fields from lightning return strokes.However, as we will show later, the energy transported by return stroke radiation fields depends on the shape, amplitude, and duration of the radiation fields.The features of the radiation fields predicted by different models may vary significantly.In the present study, we estimated the energy and momentum transported by return stroke radiation fields as predicted by two return stroke models that belong to the Modified Transmission Line (MTL) type models to identify the effects of these variations on the energy and momentum transported by return strokes.These are the MTL model with exponential current decay (MTLE) and the MTL model with linear current decay (MTLL) [15,16].According to MTLE and MTLL models, the current distribution along the return stroke channel is given by In the MTLE model, A(z) = e −z/λ , and in the MTLL model, A(z) = (1 − z/H), where λ is the current decay height constant and H is the height of the return stroke channel.In the mathematical expressions given above, z is the vertical coordinate or the height along the return stroke channel, I b (t, z) is the return stroke current at the channel base and v is the return stroke speed.Now, as mentioned previously, in addition to the initial peak, the duration of the radiation field is a parameter that influences the amount of energy and momentum radiated by return strokes.The zero-crossing time of the radiation field is a measure of its duration.For this reason, we present the results for the radiated energy and momentum as a function of the zero-crossing time and the initial peak of radiation fields.Such a presentation allows comparing the predictions of different return stroke models on radiated energy and momentum.Moreover, the peak and the zero-crossing time of the radiation fields of the return strokes can be measured in practice, and this permits us to extrapolate the results presented here to the actual return strokes using peak amplitude and zero-crossings as measures.
The current waveform used in the present study to represent the channel base current of first and subsequent return strokes is given by [17].
The parameters of this equation for the subsequent return strokes are i 01 = 10.7 kA, τ 11 = 0.25 µs, τ 21 = 2.5 µs, i 02 = 6.5 kA, τ 12 = 2.0 µs, and τ 22 = 100 µs.The first return stroke is represented using only the first term in Equation 3 with parameters i 01 = 28 kA, τ 11 = 1.5 µs, and τ 21 = 95 µs. Figure 1 depicts these current waveforms.Observe that the individual parameters in the above expression for the current cannot be connected individually to the physical current parameters such as the peak current and current risetime.However, the combination of these parameters gives rise to current waveforms that are in reasonable agreement with observations.Atmosphere 2020, 9, x FOR PEER REVIEW 4 of 18 as the peak current and current risetime.However, the combination of these parameters gives rise to current waveforms that are in reasonable agreement with observations.

Electromagnetic Fields from the Lightning Channel
Figure 2 illustrates the geometrical parameters necessary to describe the electric field and the magnetic field intensity.Let us represent the lightning channel by a straight vertical channel located over a perfectly conducting ground plane.In calculating electric and magnetic fields, an image channel replaces the perfectly conducting ground.Consider an infinitesimal channel element of length d z located at height z along the channel.The electric and magnetic fields generated by this current element at point P is given by [18,19].

Electromagnetic Fields from the Lightning Channel
Figure 2 illustrates the geometrical parameters necessary to describe the electric field and the magnetic field intensity.Let us represent the lightning channel by a straight vertical channel located over a perfectly conducting ground plane.In calculating electric and magnetic fields, an image channel replaces the perfectly conducting ground.Consider an infinitesimal channel element of length dz located at height z along the channel.The electric and magnetic fields generated by this current element at point P is given by [18,19].
The resulting fields due to the image current in the perfectly conducting ground plane is given by  Observe that the field components given above are the time domain electric and magnetic fields of a small current element of length dz that contain 1/r, 1/r 2 , and 1/r 3 terms.The components of the electric field in the direction of unit vectors a r and a θ are then given by e r (t, r, θ) = e r 1 •a r + e θ 1 •a r + e r 2 •a r + e θ 2 •a r ) a r (10) The total electric radiation fields in a r and a θ directions and the total magnetic radiation field intensity in the a ϕ direction can be obtained by summing up the contributions to these field components from all channel elements located along the return stroke channel including the image in the perfectly conducting ground.That is: Note that due to azimuthal (or cylindrical) symmetry, the fields are independent of the azimuth angle ϕ.

Expressions for the Power, Energy, and Momentum of the Electromagnetic Fields
The Poynting vector, S(t, θ, r), which is directed along the vector a r , associated with the fields, is given by This Poynting vector can be utilized to obtain the power P(t), energy U, flux of momentum in the z-direction M z (t) and the total momentum transferred along the z-direction M zT as follows [18,19]: Note that the momentum has a component only in the z-direction due to symmetry.It is worthy to point out here that at sufficiently large distances (exact distance depends on the frequency content of the channel base current and the return stroke height), the electric and magnetic fields become pure radiation.In this case, Equations ( 16)- (19) define the power, energy, and momentum radiated away by the return stroke.At shorter distances, the electromagnetic field terms varying as 1/r 2 and 1/r 3 also contribute to the radial Poynting vector.As we show later, these contributions represent the energy transferred between the source and the nearby space, and vice versa.

Results
One can deduce from the equations given in the previous section, that for a given current wave shape, the energy, momentum, and the peak power radiated by a return stroke is proportional to the square of the peak of the return stroke current.For this reason, the results are obtained only for 12 kA peak current for the subsequent strokes and 30 kA peak current for the first strokes.Using the proportionality as mentioned above, one can easily obtain the amount of energy, momentum, or peak power radiated by return strokes having other values of peak currents.However, these parameters, especially the energy and momentum, depend not only on the return stroke peak current but also on how the temporal features of the radiation fields vary at different elevations on the surface of a hemisphere centered at the point of the strike of the return stroke.For this reason, first, we will illustrate how the return stroke radiation field varies as a function of the elevation angle θ.

Radiation Field at Different Elevation Angles
Consider a hemisphere of radius R with the center at the point of strike of the return stroke.The hemisphere lies in the upper half-space of Figure 2, i.e., above the ground.Let us select the distance R such that the fields at that distance are pure radiation.Consider a cross-section of this hemisphere along the z-x plane.As one moves along this cross-section of the hemisphere, the radial distance to the point of strike remains constant at R while the angle θ varies.The angle θ is equal to π/2 at ground level, and θ = 0 directly above the strike point.We will present how the radiation field varies with an angle θ using the radiation fields of subsequent return strokes of the MTLL model.5), ( 6), ( 8) and ( 9) that the radiation field is proportional to the sin θ (note that at large distances, sin θ 1 = sin θ 2 = sin θ).Thus, the amplitude of the radiation field decreases as θ decreases.In the waveforms presented in Figure 3, this dependency is removed by normalizing the amplitude of all the electric fields to unity.Note how the radiation field becomes narrower as the angle of elevation decreases.Indeed, this effect depends on the speed of propagation of the return stroke.As the speed increases the radiation fields become narrower with decreasing angle of elevation (i.e., angle θ).
Atmosphere 2020, 9, x FOR PEER REVIEW 7 of 18 radiation field is proportional to the sin θ (note that at large distances, ).Thus, the amplitude of the radiation field decreases as θ decreases.In the waveforms presented in Figure 3, this dependency is removed by normalizing the amplitude of all the electric fields to unity.Note how the radiation field becomes narrower as the angle of elevation decreases.Indeed, this effect depends on the speed of propagation of the return stroke.As the speed increases the radiation fields become narrower with decreasing angle of elevation (i.e., angle θ ).

Energy and Momentum Radiated by Return Strokes
As mentioned previously, the energy and momentum radiated by return strokes depend also on the duration of the radiation field, and the zero-crossing time can act as a measure of this duration.In general, the zero-crossing time of the radiation field depends on the duration of the return stroke current, how the return stroke current is attenuated along the channel, return stroke speed, and the return stroke channel length [20].In the MTLL and MTLE models, once the channel base current is specified, the zero-crossing time of the radiation field is controlled by the return stroke speed, and how the return stroke current is attenuated along the channel as specified by the parameters H and λ , in the MTLL and MTLE models respectively.In obtaining the results to be presented here, we changed these parameters in the models to have the desired zero-crossing time in the radiation field.As mentioned earlier, one goal of this study is to investigate how the energy and momentum radiated by return strokes vary with the zero-crossing time of the radiation field.For a given return stroke speed, the zero-crossing time of the radiation field of the MTLL model increases with increasing H.
In the MTLE model, the same is true for the parameter λ .However, as the H in the MTLL model and λ in the MTLE model increase there is a slight increase in the radiation field peak due to the reduction of the rate of current attenuation with height in the two models.In order to get the range of zero-crossing times used in this study, the values of H and λ in the two models had to be changed over the range 2 km to about 12 km and 1 km to about 7 km, respectively.
Figure 4i depicts the peak amplitude of the radiation fields of first return strokes at 100 km as a function of the zero-crossing time as obtained using MTLL and MTLE models.Figure 4ii shows the energy radiated by the radiation fields corresponding to MTLL and MTLE models as a function of the zero-crossing time.Note that for a given zero-crossing time, the radiated energy is larger in MTLL

Energy and Momentum Radiated by Return Strokes
As mentioned previously, the energy and momentum radiated by return strokes depend also on the duration of the radiation field, and the zero-crossing time can act as a measure of this duration.In general, the zero-crossing time of the radiation field depends on the duration of the return stroke current, how the return stroke current is attenuated along the channel, return stroke speed, and the return stroke channel length [20].In the MTLL and MTLE models, once the channel base current is specified, the zero-crossing time of the radiation field is controlled by the return stroke speed, and how the return stroke current is attenuated along the channel as specified by the parameters H and λ, in the MTLL and MTLE models respectively.In obtaining the results to be presented here, we changed these parameters in the models to have the desired zero-crossing time in the radiation field.As mentioned earlier, one goal of this study is to investigate how the energy and momentum radiated by return strokes vary with the zero-crossing time of the radiation field.For a given return stroke speed, the zero-crossing time of the radiation field of the MTLL model increases with increasing H.In the MTLE model, the same is true for the parameter λ.However, as the H in the MTLL model and λ in the MTLE model increase there is a slight increase in the radiation field peak due to the reduction of the rate of current attenuation with height in the two models.In order to get the range of zero-crossing times used in this study, the values of H and λ in the two models had to be changed over the range 2 km to about 12 km and 1 km to about 7 km, respectively.
Figure 4a depicts the peak amplitude of the radiation fields of first return strokes at 100 km as a function of the zero-crossing time as obtained using MTLL and MTLE models.Figure 4b shows the energy radiated by the radiation fields corresponding to MTLL and MTLE models as a function of the zero-crossing time.Note that for a given zero-crossing time, the radiated energy is larger in MTLL than in MTLE, since the half-width of the predicted radiation field is larger in the case of MTLL than in the case of MTLE.The vertical momentum radiated by the first return strokes as a function of zero-crossing time obtained using both MTLE and MTLL models are shown in Figure 4c.Observe again, that for a given zero-crossing time, the MTLL predicts higher radiated momentum than the MTLE model.Observe from the equations given in the previous section that the energy and momentum radiated by return strokes is proportional to the square of the electric or magnetic radiation fields.Thus, one can normalize the radiated energy and momentum by dividing them by the square of the radiation field peak at 100 km.The results obtained thus, which we call the field-normalized energy and momentum, are presented in Figure 4d,e.
Atmosphere 2020, 9, x FOR PEER REVIEW 8 of 18 crossing time obtained using both MTLE and MTLL models are shown in Figure 4iii.Observe again, that for a given zero-crossing time, the MTLL predicts higher radiated momentum than the MTLE model.Observe from the equations given in the previous section that the energy and momentum radiated by return strokes is proportional to the square of the electric or magnetic radiation fields.Thus, one can normalize the radiated energy and momentum by dividing them by the square of the radiation field peak at 100 km.The results obtained thus, which we call the field-normalized energy and momentum, are presented in Figure 4iv,v.
Results pertinent to subsequent return strokes are given in Figure 5.Note again that for a given zero-crossing time, the radiation fields of the MTLL model radiate more energy and momentum than the ones pertinent to the MTLE model.The reason is, the half-width of the predicted radiation field is larger with MTLL than with MTLE.The results presented in these diagrams show that a typical first return stroke generating a radiation field having a 50 μs zero-crossing time will radiate field-normalized energy of about (1.7-2.5)× 10 3 J/(V/m) 2 and a field-normalized vertical momentum of about (2.25-3.1)× 10 −6 Kg m/s/(V/m) 2 .A radiation field with a zero-crossing time of 70 μs will radiate about (2.6-3.4) × 10 3 J//(V/m) 2 in fieldnormalized energy and (3.2-4.3)× 10 −6 Kg m/s/(V/m) 2 in field-normalized momentum.On the other hand, a typical subsequent return stroke generating a radiation field having a 40 μs zero-crossing time will radiate field-normalized energy of about (6-9) × 10 2 J//(V/m) 2 and vertical momentum of approximately (0.75-1.1) × 10 −6 Kg m/s/(V/m) 2 .If the zero-crossing time of the subsequent return stroke radiation field is reduced to 25 μs, the field-normalized energy and field-normalized  Results pertinent to subsequent return strokes are given in Figure 5.Note again that for a given zero-crossing time, the radiation fields of the MTLL model radiate more energy and momentum than the ones pertinent to the MTLE model.The reason is, the half-width of the predicted radiation field is larger with MTLL than with MTLE.
radiate energy of about (6-9) × 10 4 J and about (8-11) × 10 -5 Kg m/s in momentum-the energy and momentum radiated by return strokes increases as the square of the radiation field peak.A radiation field with the same peak amplitude but with a zero-crossing time of 70 μs will radiate about (9-12) × 10 4 J in energy and (11-15) × 10 −5 Kg m/s in momentum.On the other hand, a typical subsequent return stroke generating a radiation field having a peak amplitude of 3.5 V/m and a 40 μs zerocrossing time will radiate energy of about (7.3-11) × 10 3 J and vertical momentum of about (9-13) × 10 −6 Kg m/s.In the results presented earlier for energy and momentum, we have set the speed of propagation of the return stroke at 10 8 m/s for the first and 1.5 × 10 8 m/s for subsequent strokes.We repeated the calculations by using speeds of 2.0 × 10 8 m/s for first strokes and 2.5 × 10 8 m/s for subsequent strokes.The results obtained are shown in Figures 6 and 7, respectively, as a function of zero-crossing time.Note that even though there is a large increase in the radiated energy and momentum with increasing speeds, which indeed is caused by the increasing peak amplitude of the radiation field, the increase in the field-normalized energy and momentum is not that significant.The results presented in these diagrams show that a typical first return stroke generating a radiation field having a 50 µs zero-crossing time will radiate field-normalized energy of about (1.7-2.5)× 10 3 J/(V/m) 2 and a field-normalized vertical momentum of about (2.25-3.1)× 10 −6 Kg m/s/(V/m) 2 .A radiation field with a zero-crossing time of 70 µs will radiate about (2.6-3.4) × 10 3 J//(V/m) 2 in field-normalized energy and (3.2-4.3)× 10 −6 Kg m/s/(V/m) 2 in field-normalized momentum.On the other hand, a typical subsequent return stroke generating a radiation field having a 40 µs zero-crossing time will radiate field-normalized energy of about (6-9) × 10 2 J//(V/m) 2 and vertical momentum of approximately (0.75-1.1) × 10 −6 Kg m/s/(V/m) 2 .If the zero-crossing time of the subsequent return stroke radiation field is reduced to 25 µs, the field-normalized energy and field-normalized momentum will be reduced to (4-6.2) × 10 2 J//(V/m) 2 and (0.5-0.8) × 10 −6 Kg m/s/(V/m) 2 , respectively.
In experimental observations, it is possible to collect many electromagnetic fields from lightning return strokes in different geographical regions around the world.Thus, for an extensive collection of radiation fields, it is possible to measure the peak amplitude (normalized to 100 km) and the zero-crossing time.From such data and using the information provided in Figure 4d,e and Figure 5d,e, one can obtain the distributions of radiated energy and momentum by first and subsequent return strokes in different geographical regions.For example, a typical first return stroke generating a radiation field having a range normalized (to 100 km) peak amplitude of 6 V/m and a 50 µs zero-crossing time will radiate energy of about (6-9) × 10 4 J and about (8-11) × 10 −5 Kg m/s in momentum-the energy and momentum radiated by return strokes increase as the square of the radiation field peak.A radiation field with the same peak amplitude but with a zero-crossing time of 70 µs will radiate about (9-12) × 10 4 J in energy and (11-15) × 10 −5 Kg m/s in momentum.On the other hand, a typical subsequent return stroke generating a radiation field having a peak amplitude of 3.5 V/m and a 40 µs zero-crossing time will radiate energy of about (7.3-11) × 10 3 J and vertical momentum of about (9-13) × 10 −6 Kg m/s.In the results presented earlier for energy and momentum, we have set the speed of propagation of the return stroke at 10 8 m/s for the first and 1.5 × 10 8 m/s for subsequent strokes.We repeated the calculations by using speeds of 2.0 × 10 8 m/s for first strokes and 2.5 × 10 8 m/s for subsequent strokes.The results obtained are shown in Figures 6 and 7, respectively, as a function of zero-crossing time.Note that even though there is a large increase in the radiated energy and momentum with increasing speeds, which indeed is caused by the increasing peak amplitude of the radiation field, the increase in the field-normalized energy and momentum is not that significant.1,0x10 5   2,0x10 5   3,0x10 5   4,0x10 5   5,0x10 5   6,0x10 5

Peak Power Radiated by Return Strokes
Another parameter which is of interest in evaluating the energy balance of return strokes is the power radiated by return strokes.For a given peak current, the peak power radiated by a return stroke depends on the return stroke speed, which controls the initial peak of the radiation field.However, the zero-crossing time of the field does not significantly influence the peak power.Figure 8 depicts the temporal variation of the radiated power of a first return stroke obtained using both MTLL and MTLE models.In the calculation, the peak current was fixed at 30 kA and the return stroke speed at 1.0 × 10 8 m/s.The model parameters are selected so that the radiation fields of both models have a zero-crossing time of 50 µs.Note that the peak power is released by a first return stroke within the first few microseconds of the return stroke.For this reason, it is not influenced by the zero-crossing time of the radiation field.
8 depicts the temporal variation of the radiated power of a first return stroke obtained using both MTLL and MTLE models.In the calculation, the peak current was fixed at 30 kA and the return stroke speed at 1.0 × 10 8 m/s.The model parameters are selected so that the radiation fields of both models have a zero-crossing time of 50 μs.Note that the peak power is released by a first return stroke within the first few microseconds of the return stroke.For this reason, it is not influenced by the zerocrossing time of the radiation field.
In order to study how the peak power radiated by return strokes are affected by return stroke speed, we evaluated the power radiated by first and subsequent return strokes as a function of return stroke speed using both MTLL and MTLE models.The parameters we evaluated are the peak power radiated and the time at which the peak power is radiated.In this calculation, we set the model parameters equal to the following values: λ = 2 km and H = 7.5 km.
Two reasons cause the variation of peak power with return stroke speed.First, the peak field of the radiation field increases as the speed increases.Since the power is proportional to the square of the peak of the radiation field, the peak power increases with increasing speed.Second, as mentioned previously, the shape of the radiation field at different elevation angles changes with return stroke speed, which could possibly affect the power generation.On the other hand, as mentioned before, the zero-crossing time does not affect the peak power because the peak power is generated within about a few microseconds from the beginning of the return stroke.Figure 9i shows the peak of the first return stroke radiation field at 100 km at ground level as a function of the return stroke speed for MTLL and MTLE model, and Figure 9ii depicts the peak power generated by the return stroke as simulated using MTLE and MTLL models.Figure 9iii shows the peak power divided by the square of the radiation field peak at 100 km, while Figure 9iv depicts the time at which the first return strokes generate the peak power as a function of return stroke speed.The results pertinent to subsequent strokes are shown in Figure 10.In order to study how the peak power radiated by return strokes are affected by return stroke speed, we evaluated the power radiated by first and subsequent return strokes as a function of return stroke speed using both MTLL and MTLE models.The parameters we evaluated are the peak power radiated and the time at which the peak power is radiated.In this calculation, we set the model parameters equal to the following values: λ = 2 km and H = 7.5 km.
Two reasons cause the variation of peak power with return stroke speed.First, the peak field of the radiation field increases as the speed increases.Since the power is proportional to the square of the peak of the radiation field, the peak power increases with increasing speed.Second, as mentioned previously, the shape of the radiation field at different elevation angles changes with return stroke speed, which could possibly affect the power generation.On the other hand, as mentioned before, the zero-crossing time does not affect the peak power because the peak power is generated within about a few microseconds from the beginning of the return stroke.
Figure 9a shows the peak of the first return stroke radiation field at 100 km at ground level as a function of the return stroke speed for MTLL and MTLE model, and Figure 9b depicts the peak power generated by the return stroke as simulated using MTLE and MTLL models.Figure 9c shows the peak power divided by the square of the radiation field peak at 100 km, while Figure 9d depicts the time at which the first return strokes generate the peak power as a function of return stroke speed.The results pertinent to subsequent strokes are shown in Figure 10.Data in Figures 9iii and 10iii illustrate that the field-normalized peak power increases with return stroke speed and the time at which the peak power generation decreases with increasing return stroke speed.The latter change is caused by the narrowing of the radiation field at higher elevation angles with increasing speed.Observe also that the field-normalized peak power of first return strokes is slightly higher for MTLL model than the MTLE model whereas the corresponding values for the subsequent return strokes are almost equal.The reason for this observation is the fact that in the case of subsequent return strokes, the radiation field amplitude of subsequent return strokes at times pertinent to the occurrence of the peak power is not much affected by the differences in the current attenuation curves in the two models.
The results presented in Figures 9 and 10 provide the following information: The fieldnormalized peak power generated by a first return stroke radiation field is about 1.17 × 10 8 W//(V/m) 2 and the peak power is generated within about 5-6 μs from the beginning of the return stroke.The field-normalized peak power generated by a subsequent return stroke radiation field is about 1.26 × 10 8 W//(V/m) 2 and the peak power is generated within about (0.7-0.8) μs from the beginning of the return stroke.
In the measurements conducted by Krider and Guo [5], the peak radiation field, range normalized to 100 km, of the first return stroke is about 10 V/m, and for subsequent strokes, it is about 5 V/m.According to our estimation, a first return stroke of 10 V/m will generate a peak power of 1.2 × 10 10 W and a subsequent return stroke field of 5 V/m will generate a peak power of about 3 × 10 9 W.These values are in reasonable agreement with the estimations of (1-2) × 10 10 W and (3-5) × 10 9 W for first and subsequent return strokes, respectively, by Krider and Guo [5].

Discussion
In the data presented in the previous section, we have utilized both MTLL and MTLE models in evaluating the energy and momentum radiated by first and subsequent return strokes.The energy and momentum radiated by return strokes depend to some extent on the return stroke models.The reason for this dependency is the fact that different models assume different current attenuation curves in calculating the electromagnetic radiation fields.In general, the larger the rate of attenuation of current with height the smaller the energy and momentum generated by return strokes.For this reason, the possible upper limit of the energy and momentum that can be radiated by return strokes can be obtained by appealing to the transmission line model where the return stroke current   9c,d and 10c,d illustrate that the field-normalized peak power increases with return stroke speed and the time at which the peak power generation decreases with increasing return stroke speed.The latter change is caused by the narrowing of the radiation field at higher elevation angles with increasing speed.Observe also that the field-normalized peak power of first return strokes is slightly higher for MTLL model than the MTLE model whereas the corresponding values for the subsequent return strokes are almost equal.The reason for this observation is the fact that in the case of subsequent return strokes, the radiation field amplitude of subsequent return strokes at times pertinent to the occurrence of the peak power is not much affected by the differences in the current attenuation curves in the two models.

Data in Figures
The results presented in Figures 9 and 10 provide the following information: The field-normalized peak power generated by a first return stroke radiation field is about 1.17 × 10 8 W//(V/m) 2 and the peak power is generated within about 5-6 µs from the beginning of the return stroke.The field-normalized peak power generated by a subsequent return stroke radiation field is about 1.26 × 10 8 W//(V/m) 2 and the peak power is generated within about 0.7-0.8µs from the beginning of the return stroke.
In the measurements conducted by Krider and Guo [10], the peak radiation field, range normalized to 100 km, of the first return stroke is about 10 V/m, and for subsequent strokes, it is about 5 V/m.According to our estimation, a first return stroke of 10 V/m will generate a peak power of 1.2 × 10 10 W and a subsequent return stroke field of 5 V/m will generate a peak power of about 3 × 10 9 W.These values are in reasonable agreement with the estimations of (1-2) × 10 10 W and (3-5) × 10 9 W for first and subsequent return strokes, respectively, by Krider and Guo [10].

Discussion
In the data presented in the previous section, we have utilized both MTLL and MTLE models in evaluating the energy and momentum radiated by first and subsequent return strokes.The energy and momentum radiated by return strokes depend to some extent on the return stroke models.The reason for this dependency is the fact that different models assume different current attenuation curves in calculating the electromagnetic radiation fields.In general, the larger the rate of attenuation of current with height the smaller the energy and momentum generated by return strokes.For this reason, the possible upper limit of the energy and momentum that can be radiated by return strokes can be obtained by appealing to the transmission line model where the return stroke current propagates upwards without attenuation.For this reason, we have estimated the possible upper limits of the radiated energy and momentum as a function of zero-crossing time using the TL model.The zero-crossing time of the radiation field as predicted by the TL model depends on the return stroke speed and the return stroke channel height.In the results to be presented, we have changed the return stroke channel height to obtain the desired zero-crossing time.Moreover, the radiated energy and momentum increase with increasing return stroke speed, and the ultimate maximum values are obtained when the return stroke speed is equal to the speed of light.For this reason, in addition to calculating the radiated energy and momentum by TL model using return stroke speeds typical to those of first and subsequent return strokes, we have calculated these parameters also for return stroke speeds equal to the speed of light.The results are presented in Figures 11 and 12 for first and subsequent return strokes respectively.For comparison purposes, the results obtained previously for first and subsequent return strokes using MTLL and MTLE models are also shown in these diagrams.Note that the results obtained using the TL model significantly higher than the ones obtained using the MTLE or MTLL models.The results obtained using the TL model represent the possible upper limit for the energy and momentum that could be radiated by return strokes with typical speeds.On the other hand, the data obtained using return stroke speed equal to the speed of light represent the ultimate upper bound for these parameters as a function of the zero-crossing time of the radiation fields.
Atmosphere 2020, 9, x FOR PEER REVIEW 14 of 18 propagates upwards without attenuation.For this reason, we have estimated the possible upper limits of the radiated energy and momentum as a function of zero-crossing time using the TL model.The zero-crossing time of the radiation field as predicted by the TL model depends on the return stroke speed and the return stroke channel height.In the results to be presented, we have changed the return stroke channel height to obtain the desired zero-crossing time.Moreover, the radiated energy and momentum increase with increasing return stroke speed, and the ultimate maximum values are obtained when the return stroke speed is equal to the speed of light.For this reason, in addition to calculating the radiated energy and momentum by TL model using return stroke speeds typical to those of first and subsequent return strokes, we have calculated these parameters also for return stroke speeds equal to the speed of light.The results are presented in Figures 11 and 12 for first and subsequent return strokes respectively.For comparison purposes, the results obtained previously for first and subsequent return strokes using MTLL and MTLE models are also shown in these diagrams.Note that the results obtained using the TL model significantly higher than the ones obtained using the MTLE or MTLL models.The results obtained using the TL model represent the possible upper limit for the energy and momentum that could be radiated by return strokes with typical speeds.On the other hand, the data obtained using return stroke speed equal to the speed of light represent the ultimate upper bound for these parameters as a function of the zero-crossing time of the radiation fields.The amount of energy radiated by return strokes presented in the paper is actually the energy radiated into infinity by return strokes.This radiated energy is manifested through the action of propagates upwards without attenuation.For this reason, we have estimated the possible upper limits of the radiated energy and momentum as a function of zero-crossing time using the TL model.The zero-crossing time of the radiation field as predicted by the TL model depends on the return stroke speed and the return stroke channel height.In the results to be presented, we have changed the return stroke channel height to obtain the desired zero-crossing time.Moreover, the radiated energy and momentum increase with increasing return stroke speed, and the ultimate maximum values are obtained when the return stroke speed is equal to the speed of light.For this reason, in addition to calculating the radiated energy and momentum by TL model using return stroke speeds typical to those of first and subsequent return strokes, we have calculated these parameters also for return stroke speeds equal to the speed of light.The results are presented in Figures 11 and 12 for first and subsequent return strokes respectively.For comparison purposes, the results obtained previously for first and subsequent return strokes using MTLL and MTLE models are also shown in these diagrams.Note that the results obtained using the TL model significantly higher than the ones obtained using the MTLE or MTLL models.The results obtained using the TL model represent the possible upper limit for the energy and momentum that could be radiated by return strokes with typical speeds.On the other hand, the data obtained using return stroke speed equal to the speed of light represent the ultimate upper bound for these parameters as a function of the zero-crossing time of the radiation fields.The amount of energy radiated by return strokes presented in the paper is actually the energy radiated into infinity by return strokes.This radiated energy is manifested through the action of The amount of energy radiated by return strokes presented in the paper is actually the energy radiated into infinity by return strokes.This radiated energy is manifested through the action of distant radiation field terms that decreases with distance as 1/r.However, at close distances, there are field terms that vary with distance as 1/r 2 and 1/r 3 .Some of these field terms are directed along the a θ direction.The magnetic field has field terms that vary as 1/r 2 and directed along the vector a φ .These terms also contribute to the Poynting vector (E θ × H φ ) at close distances, and hence to the energy transfer from the source to space.However, they do not transport energy to infinity.
In order to illustrate the effect of these components, consider the transmission line model of the return stroke where the current pulse is assumed to move upwards without attenuation.Let H be the height of the return stroke channel.For brevity, let us denote the electric field components in the a θ direction as E 1 + E 2 .Here, E 1 contains the field terms varying as 1/r, and E 2 contains field terms varying as 1/r 2 and 1/r 3 .Similarly, the azimuthal component of the magnetic field can be written as H 1 + H 2 with H 1 containing the field terms varying as 1/r, and H 2 containing field terms varying as 1/r 2 .With this notation, the component E 1 H 1 transport energy to infinity, while the term E 1 H 2 + E 2 H 1 transport energy into the vicinity of the channel.
For illustrating the energy transported by these fields, we calculated the total energy transported through hemispheres of different radii (but all centered at the point of strike) by the Poynting vector of first and subsequent return strokes simulated using the transmission line model.Figure 13 shows the results.Note that more energy is transported across hemispheres of lower radii than the larger one.This means that some energy is stored in the space between these hemispheres during and after the return stroke.
written as

E H E H +
transport energy into the vicinity of the channel.
For illustrating the energy transported by these fields, we calculated the total energy transported through hemispheres of different radii (but all centered at the point of strike) by the Poynting vector of first and subsequent return strokes simulated using the transmission line model.Figure 13 shows the results.Note that more energy is transported across hemispheres of lower radii than the larger one.This means that some energy is stored in the space between these hemispheres during and after the return stroke.
The physical explanation of this energy transfer is the following: In the TL model, as the return stroke progresses, the positive charge starts accumulating at the cloud end of the return stroke channel, i.e., at height H .At the end of the return stroke, an amount of charge equal to the time integral of the channel base current, say Q, will be accumulated at the cloud end of the channel.This charge gives rise to a static electric field, and it contains static energy per unit volume given by , where s E is the magnitude of the static field at the point of observation.Thus, at the end of the return stroke, the static energy fills the empty space around the return stroke channel.This energy is transported by the field terms changing as 2 1 / r and 3 1 / r .This is what we observe when we use the TL model where the net effect of a return stroke is to transport a positive charge from ground to the upper end of the channel.However, in reality, at the beginning of the return stroke, a negative charge equivalent to −Q will be located at the top end of the channel (in reality, this is distributed in space).As the return stroke progresses, this negative charge starts decreasing with time, and at the end of the return stroke, the charge at the top end of the channel goes to zero.In this scenario, energy located in the space at the beginning of the return stroke will flow back from space into the source as the return stroke progresses, making the net energy in space at the end of the return stroke equal to zero.Actually, with dipoles, one can easily observe this back and forth energy flow from space to the source, and vice versa, caused by The physical explanation of this energy transfer is the following: In the TL model, as the return stroke progresses, the positive charge starts accumulating at the cloud end of the return stroke channel, i.e., at height H.At the end of the return stroke, an amount of charge equal to the time integral of the channel base current, say Q, will be accumulated at the cloud end of the channel.This charge gives rise to a static electric field, and it contains static energy per unit volume given by ε 0 E 2 s /2, where E s is the magnitude of the static field at the point of observation.Thus, at the end of the return stroke, the static energy fills the empty space around the return stroke channel.This energy is transported by the field terms changing as 1/r 2 and 1/r 3 .This is what we observe when we use the TL model where the net effect of a return stroke is to transport a positive charge from ground to the upper end of the channel.However, in reality, at the beginning of the return stroke, a negative charge equivalent to −Q will be located at the top end of the channel (in reality, this is distributed in space).As the return stroke progresses, this negative charge starts decreasing with time, and at the end of the return stroke, the charge at the top end of the channel goes to zero.In this scenario, energy located in the space at the beginning of the return stroke will flow back from space into the source as the return stroke progresses, making the net energy in space at the end of the return stroke equal to zero.Actually, with dipoles, one can easily observe this back and forth energy flow from space to the source, and vice versa, caused by 1/r 2 field terms, as the charge at the end of the dipole oscillates in time [19].
Observe that, as shown in Figures 9 and 10, for a given return stroke speed, the peak power generated by first return strokes is larger than the peak power radiated by subsequent return strokes.This is the case because the peak radiation field generated by a first return stroke is larger than that of a subsequent return stroke located at the same distance, due to the larger currents in the first return strokes.The peak power scales upwards as the square of the peak of the radiation field and this gives rise to a larger power in the first return strokes.When the power is normalized by dividing it by the square of the range normalized radiation field peak, this dependency of power on the radiation field peak is removed.It is interesting to see that the field-normalized peak power is larger in the subsequent return strokes in comparison to the first return strokes.We believe that the reason for this is the shorter risetime of the radiation field in the subsequent return strokes which makes the power transported at different elevation angles to add more efficiently than in the case of first return strokes where the risetime of the radiation field is longer.
It is important to point out that the momentum calculated in this paper is the net momentum associated with the radiation.Due to azimuthal (or cylindrical) symmetry, it contains only a z component.However, it should be noted that the electromagnetic fields moving in any radial direction also carry momentum, which is directed along the radial direction.This momentum could be imparted onto electric charges with which the electric field interacts.For example, if the field interacts with an azimuthally isotropic ionosphere (an ideal case), after interacting with the electromagnetic field, the charged particles in the ionosphere will experience a radial momentum.However, due to azimuthal symmetry, the net momentum of all the particles, after the interaction with the electromagnetic field, is directed along the z-direction to satisfy the conservation of momentum.It is also important to mention here that an equal but opposite momentum identical to the ones carried out by the radiation field is experienced by the accelerating electrons that give rise to the electromagnetic radiation.Observe that the momentum generated by a return stroke is always directed along the positive z-axis irrespective of whether the return stroke is negative or positive.On the other hand, the direction of movement of the electrons in the return stroke channel is opposite in the case of positive (directed towards positive z-axis) and negative polarity (directed along negative z-axis).For this reason, there is a polarity dependence on the effect of the radiation momentum on the positive and negative return strokes.

Conclusions
The results based on numerical simulations presented in this paper provide estimations of energy, momentum, and peak power generated by the radiation fields of first and subsequent return strokes.For a given peak value of the radiation field, the energy and momentum transported by the radiation fields increase with increasing zero-crossing time of the radiation fields.The energy transported by typical radiation fields in the tropics can be larger than the ones in temperate regions, not only because the tropical radiation fields can have larger peak amplitudes, as judged by the larger peak current amplitudes in the tropical regions [21], but also because they can also maintain larger zero-crossing times compared to those in other geophysical regions [13,22].For example, A typical first return stroke generating a radiation field having a 50 µs zero-crossing time will dissipate field-normalized energy of about (1.7-2.5)× 10 3 J/(V/m) 2 and field-normalized vertical momentum of approximately (2.25-3.1)× 10 −6 Kg m/s/(V/m) 2 .A radiation field with a zero-crossing time of 70 µs will dissipate about (2.6-3.4) × 10 3 J/(V/m) 2 in field-normalized energy and (3.2-4.3)× 10 −6 Kg m/s/(V/m) 2 in field-normalized vertical momentum.The results show that, for a given peak radiation field, the radiated energy and momentum increase with increasing zero-crossing time of the radiation field.The field-normalized peak power generated by a first return stroke radiation field is about 1.17 × 10 8 W/(V/m) 2 and the peak power is generated within about 5-6 µs from the beginning of the

Figure 1 .
Figure 1.(a) First and (b) subsequent return stroke channel base current waveforms used in deriving the return stroke electromagnetic fields from the MTL models.

Figure 2 .
Figure 2. Geometry pertinent to the evaluation of electromagnetic fields from the lightning return stroke channel.The diagram also shows the vectors necessary for writing down the electromagnetic fields from a small channel element.The vector ϕ a is in the direction of r θ × a a .

Figure 1 .
Figure 1.(a) First and (b) subsequent return stroke channel base current waveforms used in deriving the return stroke electromagnetic fields from the MTL models.

Figure 2 .
Figure 2. Geometry pertinent to the evaluation of electromagnetic fields from the lightning return stroke channel.The diagram also shows the vectors necessary for writing down the electromagnetic fields from a small channel element.The vector ϕ a is in the direction of r θ × a a .

Figure 2 .
Figure 2. Geometry pertinent to the evaluation of electromagnetic fields from the lightning return stroke channel.The diagram also shows the vectors necessary for writing down the electromagnetic fields from a small channel element.The vector a ϕ is in the direction of a r × a θ .

Figure 3
Figure 3 depicts the radiation fields of subsequent return stroke at two different elevation angles (90 • and 20 • ) and for two different return stroke speeds (1.5 × 10 8 m/s and 2.5 × 10 8 m/s) at 100 km distance from the strike point.The results are based on the MTLL model and we have selected the parameter H in such a way that the radiation fields at ground level (i.e., θ = 90 • ) for the two return stroke speeds have the same temporal variation (i.e., H = 3730 m for v = 1.5 × 10 8 m/s and H = 9260 m for v = 2.5 × 10 8 m/s respectively).Observe from Equations (5), (6), (8) and (9) that the radiation field is proportional to the sin θ (note that at large distances, sin θ 1 = sin θ 2 = sin θ).Thus, the amplitude of the radiation field decreases as θ decreases.In the waveforms presented in Figure3, this dependency is removed by normalizing the amplitude of all the electric fields to unity.Note how the radiation field becomes narrower as the angle of elevation decreases.Indeed, this effect depends on the speed of propagation of the return stroke.As the speed increases the radiation fields become narrower with decreasing angle of elevation (i.e., angle θ).

Figure 3 .
Figure 3. Electric radiation field of subsequent return stroke, based on the MTLL model, at two different elevation angles (a) θ = 90 0 and (b) θ = 20 0 .The results are shown for two different return stroke speeds with (a) 1.5 × 10 8 m/s and (b) 2.5 × 10 8 m/s.The parameter H in the MTLL model is changed with changing return stroke speed so that the temporal variation of the waveforms for the two speeds are identical at 90°.Observe that all field amplitudes are normalized to unity.The angle θ is defined in Figure 2.

Figure 3 .
Figure 3. Electric radiation field of subsequent return stroke, based on the MTLL model, at two different elevation angles (a) θ = 90 • and (b) θ = 20 • .The results are shown for two different return stroke speeds with (a) 1.5 × 10 8 m/s and (b) 2.5 × 10 8 m/s.The parameter H in the MTLL model is changed with changing return stroke speed so that the temporal variation of the waveforms for the two speeds are identical at 90 • .Observe that all field amplitudes are normalized to unity.The angle θ is defined in Figure 2.

Figure 4 .
Figure 4. Results pertinent to first return strokes simulated by MTLL and MTLE models with a peak return stroke current of 30 kA and a return stroke speed equal to 8 1 0 m/s.(a) The range-normalized (to 100 km) peak amplitude of the radiation field as a function of the zero-crossing time; (b) The energy radiated as a function of the zero-crossing time; (c) The momentum radiated as a function of the zero-crossing time (d) The field-normalized energy (as defined by dividing the energy by the square of the radiation field peak at 100 km); (e) The field-normalized momentum (as defined by dividing the momentum by the square of the radiation field peak at 100 km).Curve (a) corresponds to the MTLL model, and curve (b) corresponds to the MTLE model.

Figure 4 .
Figure 4. Results pertinent to first return strokes simulated by MTLL and MTLE models with a peak return stroke current of 30 kA and a return stroke speed equal to 10 8 m/s.(a) The range-normalized (to 100 km) peak amplitude of the radiation field as a function of the zero-crossing time; (b) The energy radiated as a function of the zero-crossing time; (c) The momentum radiated as a function of the zero-crossing time (d) The field-normalized energy (as defined by dividing the energy by the square of the radiation field peak at 100 km); (e) The field-normalized momentum (as defined by dividing the momentum by the square of the radiation field peak at 100 km).Curve (a) corresponds to the MTLL model, and curve (b) corresponds to the MTLE model.

Figure 5 .
Figure 5. Results pertinent to subsequent return strokes simulated by MTLL and MTLE models with a peak return stroke current of 12 kA and a return stroke speed equal to 8 1.5 10 × m/s.(a) The rangenormalized (to 100 km) peak amplitude of the radiation field as a function of the zero-crossing time; (b) The radiated energy as a function of the zero-crossing time; (c) The radiated momentum as a function of the zero-crossing time; (d) The field-normalized energy (as defined by dividing the energy by the square of the radiation field peak at 100 km); (e) The field-normalized momentum (as defined by dividing the momentum by the square of the radiation field peak at 100 km).Curve (a) corresponds to the MTLL model, and curve (b) corresponds to the MTLE model.

Figure 5 .
Figure 5. Results pertinent to subsequent return strokes simulated by MTLL and MTLE models with a peak return stroke current of 12 kA and a return stroke speed equal to 1.5 × 10 8 m/s.(a) The range-normalized (to 100 km) peak amplitude of the radiation field as a function of the zero-crossing time; (b) The radiated energy as a function of the zero-crossing time; (c) The radiated momentum as a function of the zero-crossing time; (d) The field-normalized energy (as defined by dividing the energy by the square of the radiation field peak at 100 km); (e) The field-normalized momentum (as defined by dividing the momentum by the square of the radiation field peak at 100 km).Curve (a) corresponds to the MTLL model, and curve (b) corresponds to the MTLE model.

Figure 6 .
Figure 6.Results pertinent to first return strokes simulated by MTLL and MTLE models with a peak return stroke current of 30 kA and a return stroke speed equal to 8 2 10 × m/s.

Figure 7 .
Figure 7. Results pertinent to subsequent return strokes simulated by MTLL and MTLE models with a peak return stroke current of 12 kA and a return stroke speed equal to8 2.5 10 × m/s.(a) The rangenormalized (to 100 km) peak amplitude of the radiation field as a function of the zero-crossing time; (b) The energy radiated as a function of the zero-crossing time; (c) The momentum radiated as a function of the zero-crossing time; (d) The field-normalized energy (as defined by dividing the energy by the square of the radiation field peak at 100 km); (e) The field-normalized momentum (as defined by dividing the momentum by the square of the radiation field peak at 100 km).Curve (a) corresponds to the MTLL model, and curve (b) corresponds to the MTLE model.

Figure 6 .
Figure 6.Results pertinent to first return strokes simulated by MTLL and MTLE models with a peak return stroke current of 30 kA and a return stroke speed equal to 2 × 10 8 m/s.(a) The range-normalized (to 100 km) peak amplitude of the radiation field as a function of the zero-crossing time; (b) The energy radiated as a function of the zero-crossing time; (c) The momentum radiated as a function of the zero-crossing time; (d) The field-normalized energy (as defined by dividing the energy by the square of the radiation field peak at 100 km); (e) The field-normalized momentum (as defined by dividing the momentum by the square of the radiation field peak at 100 km).Curve (a) corresponds to the MTLL model, and curve (b) corresponds to the MTLE model.

Figure 6 .
Figure 6.Results pertinent to first return strokes simulated by MTLL and MTLE models with a peak return stroke current of 30 kA and a return stroke speed equal to 8 2 10 × m/s.

Figure 7 .
Figure 7. Results pertinent to subsequent return strokes simulated by MTLL and MTLE models with a peak return stroke current of 12 kA and a return stroke speed equal to8 2.5 10 × m/s.(a) The rangenormalized (to 100 km) peak amplitude of the radiation field as a function of the zero-crossing time; (b) The energy radiated as a function of the zero-crossing time; (c) The momentum radiated as a function of the zero-crossing time; (d) The field-normalized energy (as defined by dividing the energy by the square of the radiation field peak at 100 km); (e) The field-normalized momentum (as defined by dividing the momentum by the square of the radiation field peak at 100 km).Curve (a) corresponds to the MTLL model, and curve (b) corresponds to the MTLE model.

2 Figure 7 .
Figure 7. Same as in Figure 6 except the results are pertinent to subsequent return strokes simulated by MTLL and MTLE models with a peak return stroke current of 12 kA and a return stroke speed equal to 2.5 × 10 8 m/s.

Figure 8 .
Figure 8.(a) Temporal variation of the radiated power as a function of time.Curve (a) MTLL model with a radiation field having 50 μs zero-crossing time.Curve (b) MTLE model with a radiation field having 50 μs zero-crossing time.(b) The electric field at 100 km at ground level corresponding to the power curves.

Figure 8 .
Figure 8.(a) Temporal variation of the radiated power as a function of time.Curve (a) MTLL model with a radiation field having 50 µs zero-crossing time.Curve (b) MTLE model with a radiation field having 50 µs zero-crossing time.(b) The electric field at 100 km at ground level corresponding to the power curves.

Figure 9 .
Figure 9. Data pertinent to the peak power radiated by a first return stroke supporting a channel base peak current of 30 kA.(a) The range-normalized (to 100 km) peak radiation field at ground level as a function of the return stroke speed; (b) The peak power generated by the return stroke as a function of the return stroke speed; (c) The field-normalized peak power obtained by dividing the peak power by the square of the radiation field peak at 100 km; (d) The time at which the power peaks as a function of the return stroke.(a) MTLL model; (b) MTLE model.

Figure 9 .
Figure 9. Data pertinent to the peak power radiated by a first return stroke supporting a channel base peak current of 30 kA.(a) The range-normalized (to 100 km) peak radiation field at ground level as a function of the return stroke speed; (b) The peak power generated by the return stroke as a function of the return stroke speed; (c) The field-normalized peak power obtained by dividing the peak power by the square of the radiation field peak at 100 km; (d) The time at which the power peaks as a function of the return stroke.(a) MTLL model; (b) MTLE model.

Figure 9 . 8 (Figure 10 .
Figure 9. Data pertinent to the peak power radiated by a first return stroke supporting a channel base peak current of 30 kA.(a) The range-normalized (to 100 km) peak radiation field at ground level as a function of the return stroke speed; (b) The peak power generated by the return stroke as a function of the return stroke speed; (c) The field-normalized peak power obtained by dividing the peak power by the square of the radiation field peak at 100 km; (d) The time at which the power peaks as a function of the return stroke.(a) MTLL model; (b) MTLE model.

Figure 10 .
Figure 10.Data pertinent to the peak power radiated by a subsequent return stroke supporting a channel base peak current of 12 kA.(a) The range-normalized (to 100 km) peak radiation field at ground level as a function of the return stroke speed; (b) The peak power generated by the return stroke as a function of the return stroke speed; (c) The field-normalized peak power obtained by dividing the peak power by the square of the radiation field peak at 100 km; (d) The time at which the power peaks as a function of the return stroke.(a) MTLL model; (b) MTLE model.

Figure 10 .
Figure 10.Data pertinent to the peak power radiated by a subsequent return stroke supporting a channel base peak current of 12 kA.(a) The range-normalized (to 100 km) peak radiation field at ground level as a function of the return stroke speed; (b) The peak power generated by the return stroke as a function of the return stroke speed; (c) The field-normalized peak power obtained by dividing the peak power by the square of the radiation field peak at 100 km; (d) The time at which the power peaks as a function of the return stroke.(a) MTLL model; (b) MTLE model.

Figure 11 .Figure 12 .
Figure 11.The field-normalized energy and momentum radiated by a first return stroke simulated by the transmission line model.The peak channel base current is 30 kA.(a) Field-normalized energy.(b) Field-normalized momentum.(a) MTLL model; (b) MTLE model; (c) TL model with speed equal to 10 8 m/s; (d) TL model with speed equal to the speed of light.

Figure 11 .
Figure 11.The field-normalized energy and momentum radiated by a first return stroke simulated by the transmission line model.The peak channel base current is 30 kA.(a) Field-normalized energy.(b) Field-normalized momentum.(a) MTLL model; (b) MTLE model; (c) TL model with speed equal to 10 8 m/s; (d) TL model with speed equal to the speed of light.

Figure 11 .Figure 12 .
Figure 11.The field-normalized energy and momentum radiated by a first return stroke simulated by the transmission line model.The peak channel base current is 30 kA.(a) Field-normalized energy.(b) Field-normalized momentum.(a) MTLL model; (b) MTLE model; (c) TL model with speed equal to 10 8 m/s; (d) TL model with speed equal to the speed of light.

Figure 12 .
Figure 12.The field-normalized energy and momentum radiated by a subsequent return stroke simulated by the transmission line model.The peak channel base current is 12 kA.(a) Field-normalized energy.(b) Field-normalized momentum.(a) MTLL model; (b) MTLE model; (c) TL model with speed equal to 1.5 × 10 8 m/s; (d) TL model with speed equal to the speed of light.

2 1 /Figure 13 .
Figure 13.Energy transported through hemispheres of different radii all having their center at the strike point.(a) First return strokes; (b) subsequent return strokes.