Smart Localization of Thunderstorm Charge for Human 4IR Applications

Collaborative Innovation Center for Meteorological Disaster Prediction and Evaluation, Nanjing University of Information Science and Technology, Nanjing 210044, China Jiangsu Key Laboratory of Meteorological Detection and Information Processing, Nanjing University of Information Science and Technology, Nanjing 210044, China Department of Electrical and Mining Engineering, University of South Africa, Florida 1710, South Africa


Introduction
As one of the changes caused by the global trend of 4IR, climate is an important research field. understorm cloud is a common weather phenomenon, which may cause serious disasters [1][2][3]. For buildings [4,5], the metal devices with poor grounding inside are easy to produce sparks after lightning induction, which is dangerous for warehouses storing flammable products. In the case of lightning, instantaneous high currents will be introduced into the room through transmission cables, communication lines, telephone lines, and metal pipes.
is may cause lightning strikes in computer networks, communication equipment, and industrial control systems due to electromagnetic induction. rough the analysis of lightning data, the incidence of such accidents is high and the consequences are very serious. According to estimates by the relevant departments, about 1,800 lightning strikes per minute around the world are in progress. erefore, it is urgent to study thunderstorm monitoring methods with high early warning rate [6,7]. e atmospheric electric field is often used to analyze thunderstorm activity, which has been widely used in meteorology, communication, and other fields [8][9][10]. To acquire the thunderstorm charge position, a lot of research studies have been done on the field grinding electrometer [11,12]. To determine the degree of accumulation of thunderstorm charge for lightning warning, they use the magnitude of the vertical electric field measured by AEFA. However, it is just a nondirectional detector, and it is impossible to know the specific position of the thunderstorm cloud. For example, Schumann et al. [13] analyzed the internal discharge before the positive cloud flickering to the ground. ey measured the vertical component firstly. Besides, within the range of 3-80 km from the electric field measurement point, they found that the return stroke causes some changes in electric field characteristics. By combining the interplanetary magnetic field, Frank-Kamenetsky et al. [14] studied the relationship between the thunderstorm development and the global electric field. After a long-term research, they found that there were horizontal and vertical components [15,16]. erefore, the exploration of the more real electric field is closely related to the measurement of the 3D electric field component.
From the aspects of material and structure [17][18][19][20][21], the researchers have improved the one-dimensional to threedimensional AEFA. Among them, Wen et al. [19] proposed a three-dimensional measurement method based on coplanar decoupling structure to solve the coupling problem between 3D atmospheric electric field components. Xu et al. [20] used particle swarm optimization (PSO) and 3D atmospheric electric field to invert the thunderstorm clouds and analyzed the electric field variation during a thunderstorm. Yang and Xing [21] presented a data complementary method to implement the network localization of thunderstorm charge.
ese solve the problem of being only able to measure a single electric field component by the former and reduce the electric field component coupling among existing electric field sensors. However, after using the electric field component data to obtain thunderstorm charge position information from the apparatus angle, these studies often ignore the practical significance of obtaining data from the observation angle of the apparatus. is paper focuses on the real-time and accurate calibration of thunderstorm charge smart localization in the actual environment. Based on the classical thunderstorm detection technology, this paper focuses more on implementing omnidirectional calibration relative to the observation point. e primary reasons include the following: (1) Classical detection techniques can only approximate the development of thunderstorm clouds. e root cause is that the classical AEFA can only measure a single vertical electric field component, which cannot meet the thunderstorm charge localization.
(2) ere is a lack of 3D AEFA which can simultaneously measure horizontal and vertical electric field components. In addition, the calibration of the electric field and the neglect of permittivity measurement will affect the calibration. erefore, it is necessary to design a smart localization fusion system.
(3) e data come from the apparatus itself, but the observation point is not necessarily in the position of the apparatus. is means the existence of real-time information dislocation. erefore, we need to implement the angle conversion from apparatus angle to observation angle. e contributions of this paper are as follows: (1) In order to locate and calibrate thunderstorm charge, a 3D electric field measurement model is designed according to the relative position relationship between thunderstorm charge and 3D AEFA, which is different from the classical model. (2) An all-weather calibration algorithm for the thunderstorm charge is proposed. is method focuses on the establishment of the relationship between the position of AEFA and the observation point. In other words, during the thunderstorm monitoring process, the real-time information of charge motion is obtained from the observation point. (3) We design a fusion system for smart localization.
Compared with the traditional system, the system includes 3D electric field measurement unit, electric field calibration unit, and permittivity measurement unit, which provides a guarantee for high-precision localization.

Thunderstorm Charge Localization
In this section, we establish a model of 3D electric field measurement and deduce the thunderstorm charge localization formula. In addition, the theory of electric field measurement is further given.

Establishment of 3D Electric Field Measurement Model.
Based on the mirror image theory [7,[22][23][24][25], the 3D electric field measurement model for thunderstorm charge localization as shown in Figure 1 is established. In Figure 1, M(x, y, z) is the negative charge position of the lower part of the thunderstorm cloud. N(0, 0, 0) is the position of AEFA. h represents the sum of the height of the apparatus itself and the elevation of the apparatus. (r, α, β) represents the spherical coordinates of M, and the parameters in (r, α, β) represent the distance from M to N, azimuth angle, and elevation angle in turn. Finally, the electric field intensity of M measured by AEFA is defined as E.
Secondly, the potential distribution φ of charge M at apparatus N is given by In equation (1), charge q ′ is the mirror charge of charge q. In addition, ε 1 and ε 2 represent the permittivity of the air and the ground, respectively.
In Figure 1, the three relationships of E x , E y , and E z are as follows: Mobile Information Systems (2) After performing partial derivative processing on equation (1), we get In general, the following relationships can be obtained [7,22]: It can be seen from Figure 1 that After substituting equation (4) into equation (3), we can get In equation (6), A � − q/4πε 1 and B � (ε 2 − ε 1 )/(ε 2 + ε 1 ). Equation (7) can be obtained by using equation (6): Finally, according to equation (7), M(x, y, z) can be further obtained as follows: x � r cos α cos β, y � r sin α cos β, It can be seen from Figure 1 that when the thunderstorm cloud appears above the electric field apparatus, a large amount of charge will accumulate in the cloud. At this time, the electric field E generated by the thunderstorm charge M will be sensed by the 3D electric field sensor of the apparatus. After that, the coordinates of the thunderstorm charge can be obtained by calculation, which can provide a data source for locating the charge. By contrast, the classical thunderstorm detection model is mainly based on the field mill type atmospheric electric field meter, as shown in Figure 2 [11,12].
In Figure 2, the point S represents the position of the field mill type electric field meter. Different from Figure 1, after the vertical electric field component is obtained by using this model, only the distance between charge M and point S and the degree of charge accumulation in the cloud can be roughly judged. In practical application, the position of the thunderstorm charge cannot be determined by using a onedimensional component, which leads to its poor thunderstorm warning ability. Figure 3 shows the difference between the two. e measurement of the 3D atmospheric electric field is the basis of smart localization. After analyzing the advantages and disadvantages of AEFA [7,8,[19][20][21][22], a 3D atmospheric electric field sensor structure is designed for atmospheric electric field measurement. Details are shown in Figure 4.
In Figure 4, the smart single-axis rotary vane atmospheric electric field sensor is mainly composed of a rotating blade, a shielded rotor, a photoelectric switch, a synchronous blade, a motor, an induction electrode Z, and an induction electrode X, and Y surrounding and parallel to the Z-axis. Among them, there are 8 sensing slices in both the horizontal direction and the vertical direction, and their arrangement positions are shown in Figure 5.
In Figure 5(a), there are 8 sensing slices in the vertical electrode direction. Divide these sensing slices into two groups, each with four slices, and alternately arrange them in a fan shape. At the same time, the electrodes with the same name are electrically connected, that is, Za and Zb are connected separately. In this way, when one set of sensing slices is shielded by the shielded rotor, the other set of sensing slices will be exposed to the electric field to form a differential signal.
In Figure 5(b), there are also 8 sensing slices in the horizontal electrode direction, which are divided into four groups, each with two slices. Among them, the horizontal electrodes are divided into X-direction electrodes and Ydirection electrodes. e electrodes in each direction are Model shown in Figure 1 Based on the three-dimensional atmospheric electric field apparatus Thunderstorm charge localization calibration is studied by using coordinate data Difference between the model shown in Figure 1 and the classical model divided into two groups, a and b. e electrodes with the same name are electrically connected, that is, the two pieces of Xa, Xb, Ya, and Yb are connected, respectively. After that, when one set of sensing slices is exposed to the horizontal electric field and the other set is shielded by the shielded rotor, the two sets of electrodes will form currents of opposite polarity, thus forming a pair of differential signals. In addition, the photoelectric switch is fixed on the back of the Z-direction electrode, and the hollow-shaped small blade without a fixed ring rotates synchronously with the shielded rotor. e outer layer of the 3D atmospheric electric field sensor is the shielding electrode, and its upper surface and side surfaces are hollow cylindrical metal slices. Specifically, there are 4 symmetrical hollowed-out sectors on the upper surface with 45 degrees and 4 symmetrical hollowed-out side cylindrical surfaces with 45 degrees on the circumference. At the same time, the hollow shape corresponds to the distribution of the sensing electrodes. Both the upper surface of the shielding electrode and the lower edge of the side are equipped with fixed rings to avoid vibration or deformation during rotation. When AEFA is working, the grounded shielding electrode on the sensor rotates under the drive of the motor to periodically shield the sensing electrode. At this time, a variable induced charge is generated on the sensing electrode, thereby forming 3 sets of alternating current signals, which are then, respectively, connected to the subsequent circuits.
e electrodes in all three directions are composed of two symmetrical metal pieces, which can generate differential induced currents and suppress common-mode noise. e X and Y electrodes are each composed of two rectangular metal sheets distributed at opposite positions on the sides of the cylinder. e chord width of each metal piece is equal to 1/8 of the arc length of the entire circumference. e Z electrode is composed of two metal sheets, light blue A and orange B. e shielding cover rotation frequency is f, and the induced charges are q x (t), q y (t), and q z (t). e area of a single metal piece of the electrodes X and Y is S. e inner and outer radii of the Z electrode disk are r 1 and r 2 , respectively. With the rotation of the shield, the three sensing electrodes are gradually exposed from time 0, and the electrodes are completely exposed at time T/2. During the time (T/2) − T, the three electrodes are gradually blocked by the shield. en, the induced currents i x(y) (t) and i z (t), generated by one of the metal pieces of electrode X(Y) and by the A portion of the electrode Z, respectively, are In equation (9), when the frequency f is constant, the induced current generated by electrodes X, Y, and Z has a linear relationship with the external electric field. It should be noted that the periods of the induced currents in three directions are the same. e electric field intensity can be determined by measuring the induced voltage.

Analysis of Thunderstorm Charge Localization Performance
In this section, after determining the relationship between the localization parameters and localization error, we further analyze the localization performance. e thunderstorm charge localization performance is related to the measurement error of electric field components E x , E y , and E z , charge amount q, and permittivities ε 1 and ε 2 . When charge amount q and the two permittivities above are constant, the measurement accuracy of the electric field components plays a key role in the localization performance.
Let σ E k be the standard deviation of the components, k � x, y, z, and they are equal to σ E i [26,27]. In addition, let us take the permittivities ε 1 and ε 2 to be equal to 1 and 5, respectively [20,23,24,26,27].
After using equation (7), the measurement errors of r, α, and β caused by error σ E i are [25,28,29] In equation (10), σ rE i , σ αE i , and σ βE i , respectively, are the measurement errors of r, α, and β. In particular, the reduction of error σ E i can effectively improve the localization performance.

Analysis of Ranging
Performance. According to equation (10), the error curve in Figure 6 can be obtained.
In Figure 6, ranging error σ rE i is affected by distance r and error σ E i and is greatly affected by the former. Ranging error σ rE i increases as distance r and error σ E i increase. When error σ E i is between 0 and 1 kV/m, error σ rE i is almost independent of distance r variation. Meanwhile, when σ E i is between 1 and 2 kV/m, error σ rE i increases sharply as distance r increases. In summary, error σ rE i is less than 0.12 km, which is acceptable.

Analysis of Direction-Finding Performance.
e simulation results of the direction-finding error are obtained by using equation (10), as shown in Figures 7 and 8.
In Figure 7, the azimuth angle measurement error σ αE i increases with the increase of distance r and elevation angle β. Among them, when r is less than 1 km, error σ αE i is almost unaffected by the change of elevation angle β, and the maximum error σ αE i is only 0.5 degrees. In contrast, when r is greater than 1 km, as angle β increases, the maximum error σ αE i increases to 1.9 degrees.
In Figure 8, with the increase in distance r and elevation angle β, the elevation angle measurement error σ βE i increases slowly. It is worth noting that when r is between 1 km and 2 km, as elevation angle β increases, error σ βE i slowly climbs to 0.16 degrees.

Thunderstorm Charge Localization Calibration
In this section, in order to obtain the thunderstorm charge position directly from the observation point, we establish a 6 Mobile Information Systems  Mobile Information Systems charge localization calibration model and propose its calibration method. In addition, an elevation calibration method is given preliminarily. In the actual environment, the thunderstorm charge localization cannot be intuitively transmitted to its position information, based on AEFA. To realize the calibration of charge localization, we firstly use the apparatus to measure the thunderstorm charge orientation data. en, the spatial rectangular coordinate system is reestablished based on the observation angle of AEFA. Finally, we introduce a calibration method to achieve the above purpose.

understorm Charge Localization Calibration Model.
Let observation angle p of the electric field apparatus be located in the original coordinate system, and substitute the coordinates of p for (X, Y, Z). Based on the observation angle of AEFA, the localization calibration model is shown in Figure 9. According to Figure 9, a new coordinate system is established with observation point P(X, Y, Z) as the origin, and (X, Y, Z) is converted into (0, 0, 0). Based on the coordinate system of p, thunderstorm charge position M(x, y, z) measured by the apparatus N(0, 0, 0) can be converted into (x − X, y − Y, z − Z), and it is equivalent to M(X′, Y ′ , Z ′ ). Let the thunderstorm charge azimuth angle relative to the observation angle be θ, its elevation angle be ψ, and its distance from charge M to observation point p be R ′ .

understorm Charge Localization Calibration Method.
According to the model shown in Figure 9, the localization calibration is performed, and the corrected thunderstorm charge coordinates are In Figure 9, the azimuth and elevation angles θ and ψ are obtained by using the solid geometry analysis method [30][31][32]: For any spatial position of thunderstorm charge M, if we combine equations (11) and (12), the localization calibration can be realized.
In fact, it is also necessary to study the elevation calibration where the 3D AEFA is located. Generally speaking, the height of the thunderstorm charge from sea level is two orders of magnitude higher than the elevation where the AEFA is located. erefore, when the thunderstorm charge localization is conducted, the AEFA is usually regarded as the coordinate system origin where the charge is located. In fact, the origin should be the intersection of the Z-axis where AEFA is located and the sea level. e elevation where each apparatus is installed is different, so the data without elevation calibration are very likely to directly affect the reliability of the localization result. To this end, the model of elevation calibration is shown in Figure 10.
In addition to having the parameter characteristics of the model shown in Figure 1, in Figure 10, the elevation calibration point S(0, 0, − h) is the intersection of the Z-axis where the apparatus is located and the sea level, and the 3D electric field components after calibration are E x ′ , E y ′ , and E z ′ . e corrected horizontal angle and elevation angle of the charge are α ′ and β ′ , respectively. r ′ is the distance from M to S. e electric field intensity of M at the calibration point S is E ′ . Based on the model shown in Figure 10, the corrected 3D electric field components E x ′ , E y ′ , and E z ′ are From equation (13), we can get Among them, C � (z + h)/ ������ x 2 + y 2 . e thunderstorm charge spherical coordinates M(r ′ , α ′ , β ′ ) after elevation calibration obtained by equation (14) are  Mobile Information Systems From equation (15), it can be preliminarily seen that there is a certain relationship between the elevation h and the error of thunderstorm charge localization. Here, the validation of the elevation calibration method and the comparison with related technologies will be studied in the future work.
Finally, the flowchart of the smart calibration method as shown in Figure 11 is given.

Structure Design of Smart Localization
In order to achieve smart localization, a fusion system is obtained by combining 3D electric field measurement unit, electric field calibration unit, and permittivity measurement unit. e schematic diagram is shown in Figure 12.
In Figure 12, after the electric field calibration is carried out, the two sensors first convert physical signals into electrical signals. en, through signal conditioning and photoelectric conversion, the electrical signals are converted into digital signals. Furthermore, the signals are sent to the main control system for processing. After introducing the smart localization method, the performance of the method is analyzed. Finally, the effectiveness of the method is verified by the actual experiments in different weather conditions. In addition to the electric field measurement unit given in Section 2.2, the following three parts are included: electric field calibration unit, permittivity measurement unit, and system performance parameters.

Electric Field Calibration Unit.
To ensure the accuracy of measured electric field data, we implement the calibration of the electric field apparatus before the experiments. Building a known standard uniform electric field is the key to effective calibration. For this, a uniform electric field is generated by using two parallel plate capacitors as shown in Figure 13.
In Figure 13, the upper and lower parallel electrode plates C 1 and C 2 are both square thin metal plates, and the distance between the two plates is 2 d. e electric field apparatus N is located on the center line L between the two parallel plates, and the distance from its center position to the plate edge is D. After adding DC voltages of V 0 and − V 0 to the two plates, the uniform electric field value generated between the two plates is given by After voltages V 0 and − V 0 are input into the parallel electrode plates, uniform electric fields of different intensity are generated by adjusting the voltages. On the one hand, after the real-time voltage value is measured by the data collector, the theoretical value of the electric field intensity can be calculated between the plates, according to distance 2d between the plates. On the other hand, an electric field sensor is used to convert the electric field value between the plates into a direct current voltage signal. Finally, by comparing the theoretical and measured values, the calibration of the electric field can be achieved.

Permittivity Measurement Unit.
Considering the inconvenience of sampling and real-time measurement, a device which can ensure the measurement accuracy of the permittivity contained in equation (1) is put forward by using coplanar electrode sensor. e equivalent circuit of the device is shown in Figure 14.
In Figure 14, ε r represents the permittivity. e device mainly includes the coplanar capacity sensor with 4 electrodes. Among them, electrodes 1 and 2 are the measurement electrodes and electrodes 3 and 4 are the shielding electrodes. e latter two electrodes are connected to the ground to shield the interference of the environmental electric field. e capacitance C x is the measured capacitance between electrode 1 and electrode 2. On the one hand, the difference ΔC of the capacitance C x and C air is the linear function of measuring material permittivity. On the other hand, the capacitance values C 12 , C 13 , and C 23 are fixed and can be eliminated during measurement. erefore, the permittivity ε r can be calculated as follows:

Performance
Parameter. e smart localization system can locate the thunderstorm charge with a radius of 15 km, and the error of ranging and direction finding is less than 10%. e specific performance parameters are shown in Table 1.

Experimental Measurement Results and Analyses
In Figure 15, we install the smart 3D AEFA on the roof of the School of Electronics and Information Engineering, Nanjing University of Information Science and Technology (NUIST). e apparatus is located about 28 meters above sea level. e 3D electric field components (E x , E y , E z ) (units: kV/m) are recorded by a computer to obtain thunderstorm charge position M(r ′ , α ′ , β ′ ) (units: km, degrees, degrees).
As shown in Figure 15, AEFA and microphone array [21] are used for comparative experiments, and the latter used to locate possible thunder source. Figure 15(c) shows the schematic diagram of AEFA N and the array center N1 being placed at the same point. In the actual environment, AEFA should be placed within 5 meters from the array center. It should be noted that the elevation of the shielded rotor plane on the top of the electric field sensor should be equal to that of the XOY plane in the array. Here, the positive half of the x-axis and the y-axis of the new coordinate system are adjusted to the due south and the due   east directions, respectively. At the same time, the direction of the microphone array is adjusted in the same way. Before the experiments, by using the complementary ensemble empirical mode decomposition with adaptive noise and Savitzky-Golay filtering (CEEMDAN-SG), the atmospheric electric field signal collected by AEFA is denoised [21]. More specifically, after CEEMDAN is used to decompose the electric field signal into a series of intrinsic mode function components, the signal is reconstructed after SG filtering of those noise-dominant components. After the observation angle p is set as (0.03, 0.04, − 0.01) (units: km), the following experiments are conducted, respectively.

Experiments in
understorm Weather. At 10:11 on April 9, 2019, the results of AEFA and microphone array are shown in Table 2. Meanwhile, Figure 16 shows the radar map of Nanjing radar station. e red cross marked in the radar map is the observation point.
In Table 2, the thunderstorm charge coordinates of the observation point are (0.438, − 86.97, 71.09) (unit: km, degrees, degrees). e thunderstorm charge is 86.97 degrees south by west and 0.438 km away from the observation point. It is worth noting that the elevation angle is 71.09 degrees. erefore, it can be inferred that the thunderstorm cloud already exists above the observation point. In combination with Figure 16, the fact that radar echo intensity is over 45 dBZ indicates that a large amount of charge has accumulated in the cloud above the observation point, reflecting the existence of thunderstorm. To sum up, at 10:11 on April 9, 2019, the calibration method is consistent with the radar map data, which indicates that the thunderstorm cloud appears above the observation point. is shows the feasibility of the method.
On the premise that it is feasible to obtain the thunder source position from the microphone array, in order to compare the calibration effect with Ref. [29], the observation point we select for both is (0.03, 0.04, − 0.01) (units: km, km, km). In addition, the array element spacing is set to 0.8 m. At the same time point, the measured results are shown in Table 2.
In Table 2, from the perspective of array, the sound source point coordinates are (0.0072, − 0.1368, 0.4254) (units: km). e sound source is about 86.99 degrees south by west and 0.447 km away from the array center. Besides, its elevation angle is 72.15 degrees, which is almost perpendicular to the z-axis. At the moment, it is considered that the thunder source has occurred above the observation point. is experiment proves the validity of the calibration method well and gives real-time data feedback to the thunderstorm forecast warning.

Experiments in Fair
Weather. On the same day, at 13:26 on April 9, 2019, the results measured by AEFA and microphone array are shown in Table 3. At this time, the radar map is shown in Figure 17.
In Table 3, the charge position obtained from the observation point is (0.860, − 52.32, 21.54) (unit: km, degrees, degrees). At this time, the thunderstorm charge is about 52.32 degrees south by west and 0.860 km away from the observation point. However, its elevation angle is relatively small, only reaching 21.54 degrees. is shows that, at about 13:26, according to the data shown in Table 3, there is no thunderstorm cloud above the observation point. In Figure 17, the radar echo intensity at the observation point is only 15 dBZ, which means that there is no thunderstorm above the point because there is only a small amount of charge accumulation in the cloud. In short, at 13:26 on April 9, 2019, there is no thunderstorm cloud generated at the observation point, and the radar map data are consistent with the calibration method, once again verifying the effect of this method. Besides, combined with the data in Tables 2 and 3, it can be speculated that there is thunderstorm activity nearby and it may move further southward.
At the same time point, we use Ref. [29] to locate the thunder source and compare it with the calibration method again. Similarly, the observation point we select is (0.03, 0.04, − 0.01) (units: km), and the array element spacing is set to 0.8 m. Table 3 shows the measured results.
In Table 3, from the array observation angle, the sound source point coordinates measured by the microphone array are (− 0.0352, 0.0174, − 0.0136) (units: km, km, km). e sound source is about 26.29 degrees north by east, and it is only 0.042 km away from the observation point. In particular, the elevation angle is − 19.09 degrees. is means that the sound source is below the XOY plane. In summary, there is no thunder source generated at the observation point at 13:26 on April 9, 2019. ere is a good consistency between the measurement results of the microphone array, of the radar map, and of the results by the calibration method. Better application effects are once again verified.

Contrast Experiments.
Firstly, it can be seen from Sections 6.1 and 6.2 that when the apparatus is in an open environment, the data obtained by the observation angle can be well matched with the real-time radar map and the microphone array. e results are still in radar range, and the measured data are relatively reliable.
Next, we select the experiments in thunderstorm weather, and the calibration method is compared with references [20,33]. Table 4 shows the contrast results.
As can be seen from Table 4, all of them can measure the vertical component. However, it is obvious that Ref. [33] cannot measure the horizontal component. e most direct

12
Mobile Information Systems reason is that they used a one-dimensional AEFA. In contrast, both this paper and Ref. [20] can simultaneously measure the 3D electric field components. In fact, after measuring the 3D electric field components, the charge position may not be obtained. On the one hand, there is a lack of models that can reflect the spatial relationship between the charge and AEFA. On the other hand, there is no further derivation of the charge localization formula. Compared with references [20,33], combined with Tables 2  and 3, this paper carries out the calibration work on the premise of charge localization, which has good continuity. In terms of precision compensation, unlike Ref. [20], the paper adopts an array-based sound source localization method, which can provide a data source for the distributed multisource thunderstorm monitoring system. Finally, based on the theory of indirect measurement, this paper studies the charge localization performance and further verifies the method. However, references [20,33] perform poorly in validation.
In fact, the calibration method can also be verified by using the data measured by several AEFA which are not in the same place. According to Section 5.1 of Ref. [21], the data shown in Table 5 can be obtained. e data of apparatus N1, N2, and N3 are given in Table 5. We assume that N1 is the main AEFA and N2 and N3 are two different observation points. It can be seen from Ref. [21] that the positions of N2 and N3 are (0, − 1, 0) and (− 1, 1, 0) (units: km), respectively. Based on this, after using the calibration method, the thunderstorm charge positions at observation points N2 and N3 are (0.094, 0.862, 0.533) and (1.094, − 1.138, 0.533) (units: km), respectively. At the same time, according to the electric field data measured by N2 and N3, the thunderstorm charge positions are (0.085, 0.851, 0.486) and (1.102, − 1.127, 0.519) (units: km), respectively. In order to further compare the data obtained by the calibration method with the data measured by AEFA at the observation point, the ranging and direction-finding error rates between them are given in Table 6. It should be noted that the Cartesian coordinate data have been transformed into the spherical coordinate data (unit: km, degrees, degrees) shown in Table 6.
It can be seen from Table 6 that the errors of ranging and direction finding are within 3.40% and 6.24%, respectively, which are acceptable. is not only verifies the validity of the calibration method but also reflects the better performance of AEFA.
e above experiments better reproduce the process from the generation of the thunderstorm cloud to their disappearance. e results of the calibration method are close to those of radar map and microphone array. In general, this method has a good effect.

Conclusions
e meteorological observation data have the characteristics of time series, and the number of these data is huge. e intellectualization of algorithms and equipment for observation is an indispensable part of the 4IR. Based on the 3D AEFA, we propose a calibration method of the thunderstorm charge localization. And then, from the observation point, a new coordinate system is established to conduct the charge localization calibration. Finally, the effectiveness of the method is verified by experiments in thunderstorm and fair weather. At the same time, the relevant references are selected and compared with the method, and the results are in good agreement. e paper provides an effective solution to solve the problem that the one-dimensional AEFA can only measure a single electric field component. In addition to being used for thunderstorm detection, the proposed method can also ensure the importance of obtaining data from the observation angle. It should be noted that the existing research shows that the joint work of multiple devices can produce better results in 4IR, while the proposed method is based on a single AEFA. erefore, whether the monitoring range and accuracy can be improved based on the network observation  Tables 2 and 3 Microphone array Equation (10) Ref. [20] Equations (18) and (19) Equation (20) -Particle swarm optimization - Ref. [33] -0.637 kV/m ---  of thunderstorm with multiple devices will be carried out in the future.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.