Next Article in Journal
Data Assimilation of Ambient Concentrations of Multiple Air Pollutants Using an Emission-Concentration Response Modeling Framework
Previous Article in Journal
Subsistence Agriculture Productivity and Climate Extreme Events
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

Department of Electrical Engineering, Uppsala University, 751 21 Uppsala, Sweden
*
Author to whom correspondence should be addressed.
Atmosphere 2020, 11(12), 1288; https://doi.org/10.3390/atmos11121288
Submission received: 21 October 2020 / Revised: 23 November 2020 / Accepted: 23 November 2020 / Published: 29 November 2020
(This article belongs to the Section Meteorology)

Abstract

:
Electromagnetic radiation fields generated by return strokes transport both energy and momentum from the return stroke to outer space. The momentum transported by the radiation field 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 fields 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 field peak at 100 km) are the following: 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) × 103 J/(V/m)2 and field-normalized vertical momentum of approximately (2.3–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) × 103 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 normalized peak power generated by a first return stroke radiation field is about 1.2 × 108 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 field having a 40 μs zero-crossing time will dissipate field-normalized energy of about (6–9) × 102 J/(V/m)2 and field-normalized vertical momentum of approximately (7.5–11) × 10−7 Kg m/s/(V/m)2. The field-normalized peak power generated by a subsequent return stroke radiation field is about 1.26 × 108 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.

1. 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 × 104 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 × 104 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) × 1010 W and (3–5) × 109 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].

2. 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
I ( t , z ) = I b ( t z / v ) A ( z )   for   t z / v
I ( t , z ) = 0   for   t < z / v
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].
i ( t ) = i 01 ( t / τ 11 ) 2 ( t / τ 11 ) 2 + 1 e t / τ 21 + i 02 ( t / τ 12 ) 2 ( t / τ 12 ) 2 + 1 e t / τ 22
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.

3. 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].
e r 1 ( t , r 1 , θ 1 ) = d z 2 π ε 0 cos θ 1 [ 1 c r 1 2 I ( t r 1 / c , z ) + 1 r 1 3 t r 1 / c t I ( τ , z ) d τ ] a r 1
e θ 1 ( t , r 1 , θ 1 ) = d z 4 π ε 0 sin θ 1 [ 1 c 2 r 1 I ( t r 1 / c , z ) t + 1 c r 1 2 I ( t r 1 / c , z ) + 1 r 1 3 t r 1 / c t I ( τ , z ) d τ ] a θ 1
h φ r ( t , r 1 , θ 1 ) = d z 4 π sin θ 1 [ 1 c r 1 I ( t r 1 / c , z ) t + I ( t r 1 / c , z ) r 1 2 ] a φ
The resulting fields due to the image current in the perfectly conducting ground plane is given by
e r 2 ( t , r 2 , θ 2 ) = d z 2 π ε 0 cos θ 2 [ 1 c r 2 2 I ( t r 2 / c , z ) + 1 r 2 3 t r 2 / c t I ( τ , z ) d τ ] a r 2
e θ 2 ( t , r 2 , θ 2 ) = d z 4 π ε 0 sin θ 2 [ 1 c 2 r 2 I ( t r 2 / c , z ) t + 1 c r 2 2 I ( t r 2 / c , z ) + 1 r 2 3 t r 2 / c t I ( τ , z ) d τ ] a θ 2
h φ i ( t , r 2 , θ 2 ) = d z 4 π sin θ 2 [ 1 c r 2 I ( t r 2 / c , z ) t + I ( t r 2 / c , z ) r 2 2 ] a φ
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/r2, and 1/r3 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
e θ ( t , r , θ ) = [ e r 1 · a θ + e θ 1 · a θ + e r 2 · a θ + e θ 2 · a θ ] a θ
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:
E r ( t , r , θ ) = 0 H e r ( t , r , θ ) d z
E θ ( t , r , θ ) = 0 H e θ ( t , r , θ ) d z
H φ ( t , r , θ ) = 0 H { h φ r ( t , r 1 , θ 1 ) + h φ i ( t , r 2 , θ 2 ) } d z
Note that due to azimuthal (or cylindrical) symmetry, the fields are independent of the azimuth angle φ .

4. 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
S ( t , r , θ ) = E θ ( t , r , θ ) × H φ ( t , r , θ )
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 z T as follows [18,19]:
P ( t ) = φ = 0 2 π θ = 0 π / 2 S ( t , r , θ ) r 2 sin θ d θ d φ
U = 0 T P ( t ) d t
M z ( t ) = 1 c φ = 0 2 π θ = 0 π / 2 S ( t , r , θ ) r 2 sin θ cos θ d θ d φ
M z T = 0 T M z ( t ) d t
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.

5. 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 θ .

5.1. 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.
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 × 108 m/s and 2.5 × 108 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 × 108 m/s and H = 9260 m for v = 2.5 × 108 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 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 θ ).

5.2. 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.
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) × 103 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) × 103 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) × 102 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) × 102 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) × 104 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) × 104 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) × 103 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 108 m/s for the first and 1.5 × 108 m/s for subsequent strokes. We repeated the calculations by using speeds of 2.0 × 108 m/s for first strokes and 2.5 × 108 m/s for subsequent strokes. The results obtained are shown in Figure 6 and Figure 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.

5.3. 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 × 108 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.
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 Figure 9c,d and Figure 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.
The results presented in Figure 9 and Figure 10 provide the following information: The field-normalized peak power generated by a first return stroke radiation field is about 1.17 × 108 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 × 108 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 × 1010 W and a subsequent return stroke field of 5 V/m will generate a peak power of about 3 × 109 W. These values are in reasonable agreement with the estimations of (1–2) × 1010 W and (3–5) × 109 W for first and subsequent return strokes, respectively, by Krider and Guo [10].

6. 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 Figure 11 and Figure 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 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.
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 s 2 / 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 Figure 9 and Figure 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-xis). For this reason, there is a polarity dependence on the effect of the radiation momentum on the positive and negative return strokes.

7. 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) × 103 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) × 103 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 × 108 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 field having a 40 μs zero-crossing time will dissipate field-normalized energy of about (6–9) × 102 J/(V/m)2 and field-normalized vertical momentum of approximately (7.5–11) × 10−7. The field-normalized peak power generated by a subsequent return stroke radiation field is about 1.26 × 108 W/(V/m)2 and the peak power is generated within about 0.75–0.82 μs from the beginning of the return stroke. From these data provided in this paper, the energy or momentum radiated by a return stroke could be estimated if the initial peak of the radiation field is known or calculated. The paper also provided the upper bounds for the energy and momentum that could be radiated by return strokes.

Author Contributions

Conceptualization, V.C.; methodology, V.C.; software, V.C.; data acquisition, V.C.; validation, V.C. and A.L.; writing—original draft preparation, V.C.; writing—review and editing, V.C. and A.L.; project administration, V.C.; funding acquisition, V.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by Swedish Research Council grant VR-2015-05026 and partly by the funds from the B. John F. and Svea Andersson donation at Uppsala University.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rachidi, F.; Nucci, C.A.; Ianoz, M. Transient analysis of multiconductor lines above a lossy ground. IEEE Trans. Power Deliv. 1999, 14, 294–302. [Google Scholar] [CrossRef]
  2. Ye, M.; Cooray, V. Propagation effects caused by a rough ocean surface on the electromagnetic fields generated by lightning return strokes. Radio Sci. 1994, 29, 73–85. [Google Scholar]
  3. Cooray, V.; Ming, Y. Propagation effects on the lightning-generated electromagnetic fields for homogeneous and mixed sea-land paths. J. Geophys. Res. 1994, 99, 10641–10652. [Google Scholar] [CrossRef]
  4. Cooray, V. On the various approximations to calculate lightning return stroke generated electric and magnetic fields over finitely conducting ground. In Lightning Electromagnetics; Cooray, V., Ed.; IET: Stevenage, UK, 2012. [Google Scholar]
  5. Cooray, V. Lightning return stroke models for electromagnetic field calculations. In Lightning Interaction with Power Systems, Volume 1; Piantini, A., Ed.; IET publishers: London, UK, 2020. [Google Scholar]
  6. Boccippio, D.J.; Williams, E.R.; Heckman, S.J.; Lyons, W.A.; Baker, I.T.; Boldi, R. Sprites, ELF Transients, and Positive Ground Strokes. Science 1995, 269, 1088–1091. [Google Scholar] [CrossRef]
  7. Inan, U.S.; Piddyachiy, D.; Peter, W.B.; Sauvaud, J.A.; Parrot, M. DEMETER satellite observations of lightning-induced electron precipitation. Geophys. Res. Lett. 2007, 34. [Google Scholar] [CrossRef]
  8. Taylor, W.L.; Jean, A.G. Very-Low-Frequency Radiation Spectra of Lightning Discharges. J. Res. Natl. Bur. Stand.-D Radio Propag. 1959, 63, 199. [Google Scholar] [CrossRef]
  9. Taylor, W.L. Radiation field characteristics of lightning discharge in the band 1 kc/s. J. Res. Natl. Bur. Stand.-D Radio Propag. 1963, 67, 539–549. [Google Scholar] [CrossRef]
  10. Krider, E.P.; Guo, C. The peak electromagnetic power radiated by lightning return strokes. J. Geophys. Res. 1983, 88, 8471–8474. [Google Scholar] [CrossRef]
  11. Krider, E.P. On the electromagnetic fields, Poynting vector, and peak power radiated by lightning return strokes. J. Geophys. Res. 1992, 97, 15,913–15,917. [Google Scholar] [CrossRef]
  12. Uman, M.A.; McLain, D.K. Magnetic field of lightning return stroke. J. Geophys. Res. 1969, 74, 6899–6910. [Google Scholar] [CrossRef]
  13. Cooray, V.; Lundquist, S. Characteristics of the radiation fields from lightning in Sri Lanka in the tropics. J. Geophys. Res. 1985, 90, 6099–6109. [Google Scholar] [CrossRef]
  14. Cooray, V. The power, energy, momentum, and the action of the electromagnetic radiation fields of lightning return strokes. In Proceedings of the 34th International Conference on Lightning protection, Rzeszow, Poland, 2–7 September 2018. [Google Scholar]
  15. Nucci, C.A.; Mazzetti, C.; Rachidi, F.; Ianoz, M. On lightning return stroke models for LEMP calculations. In Proceedings of the 19th International Conference on Lightning Protection, Graz, Austria, 25–29 April 1988. [Google Scholar]
  16. Rakov, V.A.; Dulzon, A.A. A modified transmission line model for lightning return stroke field calculation. In Proceedings of the 9th International Zurich Symposium and Technical Exhibition on Electromagnetic Compatibility, Zurich, Switzerland, 12–14 March 1991; pp. 229–235. [Google Scholar]
  17. Delfino, F.; Procopio, R.; Rossi, M.; Rachidi, F.; Nucci, C.A. An algorithm for the exact evaluation of the underground lightning electromagnetic fields. IEEE Trans. (EMC) 2007, 49, 401–411. [Google Scholar]
  18. Cooray, G.; Cooray, V. Electromagnetic fields of a short electric dipole in free space—Revisited. Prog. Electromagn. Res. Lett. 2012, 131, 357–373. [Google Scholar] [CrossRef] [Green Version]
  19. Jackson, J.D. Classical Electrodynamics; John Wiley & Sons: New York, NY, USA, 1975. [Google Scholar]
  20. Cooray, V.; Orville, R.E. The effects of variation of current amplitude, current rise time and return stroke velocity along the return stroke channel on the electromagnetic fields generated by return strokes. J. Geophys. Res. 1990, 95, 18617–18630. [Google Scholar] [CrossRef]
  21. Visacro, S.; Mesquita, C.R.; De Conti, A.; Silveira, F.H. Updated statistics of lightning currents measured at Morro do Cachimbo station. Atmos. Res. 2012, 117, 55–63. [Google Scholar] [CrossRef]
  22. Lin, Y.T.; Uman, M.A.; Tiller, J.A.; Brantley, R.D.; Beasley, W.H.; Krider, E.P.; Weidman, C.D. Characterization of Lightning Return Stroke Electric and Magnetic Fields from Simultaneous Two-Station Measurements. J. Geophys. Res. 1979, 84, 6307–6314. [Google Scholar] [CrossRef]
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 1. (a) First and (b) subsequent return stroke channel base current waveforms used in deriving the return stroke electromagnetic fields from the MTL models.
Atmosphere 11 01288 g001
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 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 θ .
Atmosphere 11 01288 g002
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 × 108 m/s and (b) 2.5 × 108 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. 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 × 108 m/s and (b) 2.5 × 108 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.
Atmosphere 11 01288 g003
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 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.
Atmosphere 11 01288 g004
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 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.
Atmosphere 11 01288 g005
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. 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.
Atmosphere 11 01288 g006
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 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.
Atmosphere 11 01288 g007
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. (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.
Atmosphere 11 01288 g008
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. 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.
Atmosphere 11 01288 g009
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. 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.
Atmosphere 11 01288 g010aAtmosphere 11 01288 g010b
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 108 m/s; (d) TL model with speed equal to the speed of light.
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 108 m/s; (d) TL model with speed equal to the speed of light.
Atmosphere 11 01288 g011
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 × 108 m/s; (d) TL model with speed equal to the speed of light.
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 × 108 m/s; (d) TL model with speed equal to the speed of light.
Atmosphere 11 01288 g012
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.
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.
Atmosphere 11 01288 g013
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cooray, V.; Lobato, A. The Energy, Momentum, and Peak Power Radiated by Negative Lightning Return Strokes. Atmosphere 2020, 11, 1288. https://doi.org/10.3390/atmos11121288

AMA Style

Cooray V, Lobato A. The Energy, Momentum, and Peak Power Radiated by Negative Lightning Return Strokes. Atmosphere. 2020; 11(12):1288. https://doi.org/10.3390/atmos11121288

Chicago/Turabian Style

Cooray, Vernon, and Andre Lobato. 2020. "The Energy, Momentum, and Peak Power Radiated by Negative Lightning Return Strokes" Atmosphere 11, no. 12: 1288. https://doi.org/10.3390/atmos11121288

APA Style

Cooray, V., & Lobato, A. (2020). The Energy, Momentum, and Peak Power Radiated by Negative Lightning Return Strokes. Atmosphere, 11(12), 1288. https://doi.org/10.3390/atmos11121288

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop