Inversion of a Thunderstorm Cloud Charging Model Based on a 3D Atmospheric Electric Field

: The measurement of the atmospheric electric ﬁeld is of great signiﬁcance for the study of thunderstorm cloud charge models. Traditional electric ﬁeld meters can only measure the vertical component of the atmospheric electric ﬁeld, and thus it is difﬁcult to invert the structure of the thunderstorm cloud. A three-dimensional atmospheric electric ﬁeld meter was developed to simultaneously measure the horizontal and vertical components of the atmospheric electric ﬁeld in this paper. The effective measurement linearly relates the measured electric ﬁeld to the induced voltage, and the nonlinear equations of the three-dimensional atmospheric electric ﬁeld and the thunderstorm cloud-charging model parameters were derived. The particle swarm optimization algorithm (PSO) and the three-dimensional atmospheric electric ﬁeld were used to invert the thunderstorm clouds. Experimental observations of the three-dimensional electric ﬁeld in a cloud during a thunderstorm were analyzed. Combined with the typical charged structure model, parameters such as the charge and relative distance of the thunderstorm cloud were determined. The results showed that the value of the inversion ﬁtness function reached 0.7288, and the charge structure was even. The measurement of the three-dimensional atmospheric electric ﬁeld provides a new means of observation for the study of atmospheric electricity.


Introduction
The charging model of thunderstorm clouds refers to the use of parameters to describe the characteristics of the charge distribution inside a cloud. Atmospheric electric field is an important parameter of atmospheric physics. The average intensity of the atmospheric electric field near the ground is about 130 V/m on fine days. The atmospheric electric field will increase significantly during the formation of thunderstorm clouds. The inversion of the thundercloud charge model from ground electric field measurement data has important theoretical value for monitoring lightning and early warning research.
The most direct way to study the charge structure of thunderstorm clouds is to use sounding balloons to detect the thunderstorm clouds and obtain the electric field distribution characteristics of the detection equipment on the flight path. Simpson, Marshall, and others proposed the model of the tripole charge structure of a thunderstorm cloud by using sounding balloon observations [1]. Due to the short flying time of the sounding balloon, the measured electric field is the instantaneous charge structure of the thunderstorm cloud, so it is difficult to realize continuous monitoring. In recent years, scholars have used numerical simulation methods combined with the micro-dynamic factors of cloud physics to simulate the charge structure characteristics during the development of thunderstorm clouds [2][3][4][5][6]. However, most of these methods are simulations, and there is no physical quantity that combines actual observations. Atmospheric electric field meters play an important role in the detection of thunderstorm clouds. Studies [7][8][9][10] have looked at the improvement of atmospheric electric field meters and the method of electric field-based lightning warning. However, current techniques can only measure the vertical electric field component, ignoring the horizontal electric field caused by the charge of the thunderstorm cloud. The literature [11][12][13] has proposed the classical thunderstorm cloud charge model and used the atmospheric electric field data to study the structural parameters of the model. The data used are still derived from an atmospheric electric field that can only measure the vertical component [14]. As the horizontal component is neglected, it is necessary to invert the data from the joint observations of atmospheric electric field meters at different locations. However, the atmospheric electric field is greatly affected by the building and surface conditions at the installation site, leading to a difference in each atmospheric electric field measurement. The consistency is not high, and the quality of the data participating in the inversion is poor [15].
This paper attempted to develop a three-dimensional atmospheric electric field meter capable of simultaneously measuring horizontal and vertical components, and more comprehensively establish a functional relationship between the thunderstorm cloud charge structure and the electric field of the observation point. Combined with the particle swarm algorithm, the measurement data inversion using a single three-dimensional electric field meter was used. The structural parameters of thunderstorm clouds provides a new technical means for the study of the charge structure of thunderstorm clouds.
In this paper, the electric field generated by thunderstorm clouds on the ground was observed by using the developed three-dimensional electric field meter, and compared with the radar observation data. The experimental results showed that the trajectory of the thunderstorm cloud could be analyzed by using the three-dimensional electric field data. At the same time, using the particle swarm algorithm (PSO) and the three-dimensional atmospheric electric field data, the thunderstorm cloud-charged structure of the observation point (Nanjing) can exhibit the classic dipole. These conclusions are not achievable with conventional single-dimensional electric field meters.
In Section 2, we analyzed the measurement and calibration principle of the three-dimensional electric field meter we developed. A method for inversion of thunderstorm cloud charge structure using ground three-dimensional electric field data was proposed. In Section 3, the observation experiment of a thunderstorm cloud was carried out by using the developed three-dimensional electric field meter, and the experimental results were analyzed. In Section 4, the advantages of the three-dimensional atmospheric electric field meter in the lightning warning and the inversion of the charged structure parameters were compared with the traditional single-dimensional atmospheric electric field meter. In the last part, the application prospect of 3D electric field meter was prospected, and its limitations were pointed out.

Principles of Measurement and Calibration
The physical and structural details of the three-dimensional electric field sensor we developed are shown in Figure 1. It consisted mainly of a rotating shield; an induction electrode perpendicular to the radial X, radial Y, and axial Z directions; a photoelectric switch; a synchronous blade; and a motor. The sensing electrodes X and Y were composed of two opposite rectangular metal sheets distributed on the side of the circular body. The sensing electrode Z was a metal disk that was equally divided into eight sectors. The three sensing electrodes were perpendicular to each other. When the motor drove the shield to rotate at a constant speed, the electrodes X, Y, Z and the atmospheric electric field were periodically shielded and contacted, and an induced current was generated on the electrodes. The magnitude of the induced current was proportional to the atmospheric electric field, and the small blades were continuously driven by the motor. The shield rotated synchronously, so that the photoelectric switch generated a square wave signal of the same frequency and the same phase, which was convenient for phase sensitive detection of the induced voltage. rotated synchronously, so that the photoelectric switch generated a square wave signal of the same frequency and the same phase, which was convenient for phase sensitive detection of the induced voltage. The traditional atmospheric electric field meter only had an electrode Z and a corresponding rotating shield on the bottom of the cylinder. Because the shield rotated periodically on the bottom surface, the Z electrode was periodically exposed to the vertical electric field, resulting in induced current. This kind of electric field meter could only measure the vertical component of the atmospheric electric field and was suitable for measuring the atmospheric electric field on sunny days, because the horizontal component of the atmospheric electric field was very weak on sunny days. The two pairs of parallel electrodes were added to the side of the cylinder and the shield was extended along the side. When the shield was rotated, the electrodes on both sides and the bottom of the cylinder were periodically exposed to the space electric field, and the induced current was generated at the same time. It was suitable for detecting horizontal and vertical electric fields in thunderstorms.
The electrodes in all three directions were composed of two symmetrical metal pieces, which could generate differential induced currents and suppress common mode noise. The X and Y electrodes were each composed of two rectangular metal sheets distributed at opposite positions on the sides of the cylinder. The width of each metal piece (chord) corresponded to an arc length of 1/8 of the entire circumference. The Z electrode was composed of two metal sheets, A and B (blue fans A and brown fans B in Figure 1). The shielding cover rotation frequency was f ; the atmospheric dielectric constant was 0 ε ; the electric field intensity in the three directions were x E , y E , and z E ; and the induced charges were ( ) q t , and ( ) z q t . The area of a single metal piece of the X and Y electrodes was S . The outer radius of the Z electrode disk was 2 r , and the inner radius was 1 r .
With the rotation of the shield, the three sensing electrodes were gradually exposed from time 0, and the electrodes were completely exposed at time 2 T . During the time 2 T T − , the three electrodes were gradually blocked by the shield. Then, the induced currents ( ) ( ) x y i t and ( ) z i t , generated by one of the metal pieces of the X (Y) electrode and by the A portion of the Z electrode, respectively, were: The traditional atmospheric electric field meter only had an electrode Z and a corresponding rotating shield on the bottom of the cylinder. Because the shield rotated periodically on the bottom surface, the Z electrode was periodically exposed to the vertical electric field, resulting in induced current. This kind of electric field meter could only measure the vertical component of the atmospheric electric field and was suitable for measuring the atmospheric electric field on sunny days, because the horizontal component of the atmospheric electric field was very weak on sunny days. The two pairs of parallel electrodes were added to the side of the cylinder and the shield was extended along the side. When the shield was rotated, the electrodes on both sides and the bottom of the cylinder were periodically exposed to the space electric field, and the induced current was generated at the same time. It was suitable for detecting horizontal and vertical electric fields in thunderstorms.
The electrodes in all three directions were composed of two symmetrical metal pieces, which could generate differential induced currents and suppress common mode noise. The X and Y electrodes were each composed of two rectangular metal sheets distributed at opposite positions on the sides of the cylinder. The width of each metal piece (chord) corresponded to an arc length of 1/8 of the entire circumference. The Z electrode was composed of two metal sheets, A and B (blue fans A and brown fans B in Figure 1). The shielding cover rotation frequency was f ; the atmospheric dielectric constant was ε 0 ; the electric field intensity in the three directions were E x , E y , and E z ; and the induced charges were q x (t), q y (t), and q z (t). The area of a single metal piece of the X and Y electrodes was S. The outer radius of the Z electrode disk was r 2 , and the inner radius was r 1 . With the rotation of the shield, the three sensing electrodes were gradually exposed from time 0, and the electrodes were completely exposed at time T 2 . During the time T 2 − T, the three electrodes were gradually blocked by the shield. Then, the induced currents i x(y) (t) and i z (t), generated by one of the metal pieces of the X (Y) electrode and by the A portion of the Z electrode, respectively, were: It could be seen that when the rotational frequency f of the motor was a constant value, the induced current generated by the X, Y and Z electrodes was linear with the external electric field. Furthermore, the periods of the induced currents in the three directions were identical. The signal was calculated by measuring the induced voltage to determine the electric field intensity.
The structure of the sensor causes distortion of the electric field near the sensing electrode, so the output voltage of the sensing electrode in each direction is not only related to the electric field in the direction perpendicular to the electrode, but is also weakly affected by the electric field in other directions. For the three-dimensional electric fields E x , E y , and E z , the measured induced voltage amplitudes of the three electrodes were U x , U y , and U z , respectively. When the electric field intensity in all three directions was 0, the induced voltages generated by the three electrodes were U x0 , U y0 , and U z0 , respectively. Thus: where the matrix of a i , b i , and c i represented calibration coefficients. To determine the calibration values, a DC voltage was applied to a parallel plate capacitor to generate standard electric fields, which were perpendicular to the X, Y, and Z electrodes, and the actual output voltage of the sensor was measured. Based on the measured voltage and voltage values, the calibration factor was calculated using Equation (3).

The Relationship between a Three-Dimensional Atmospheric Electric Field and Thunderstorm Cloud Charge
The three-dimensional atmospheric electric field meter installed in the field is shown in Figure 2a. To solve the relationship between the atmospheric electric field at the observation point and the thunderstorm cloud charge, an observation model (Figure 2b) was established. The X-Y plane was mounted on the ground; A (x, y, z) was an arbitrary observation point; B (x 0 , y 0 , z 0 ) was the equivalent charge point of thunderstorm cloud; and the charge of the thunderstorm cloud was q. The horizontal azimuth of the thunderstorm cloud B relative to the observation point A was α, the elevation angle was γ, and A' and B' were projections of the observation station A and the thunderstorm cloud B on the X-Y plane, respectively. The dielectric constants of region 0 above the surface X-Y and region 1 inside the earth were 0 ε and 1 ε , respectively. The electric field of region 0 was generated by the thunderstorm cloud charge q and the induced charge of the surface. The induced charge was equivalent to the image charge 1 q .
The distance between the observation point A and charge q and its image charge 1 q were 0 R and 1 R , respectively. The electric field of region 1 was generated by charge q and the image charge 0 q of region 0, and the distance of the image charge 0 q from the observation point A was also 0 R . The potentials at any point above the surface ( 0 z > ) and inside the earth ( 0 z < ) are shown in Equations (6) and (7), respectively: According to the geometric relationship of Figure 2b, the horizontal azimuth angle α and the elevation angle γ of the thunderstorm cloud satisfied: The dielectric constants of region 0 above the surface X-Y and region 1 inside the earth were ε 0 and ε 1 , respectively. The electric field of region 0 was generated by the thunderstorm cloud charge q and the induced charge of the surface. The induced charge was equivalent to the image charge q 1 . The distance between the observation point A and charge q and its image charge q 1 were R 0 and R 1 , respectively. The electric field of region 1 was generated by charge q and the image charge q 0 of region 0, and the distance of the image charge q 0 from the observation point A was also R 0 . The potentials at any point above the surface (z > 0) and inside the earth (z < 0) are shown in Equations (6) and (7), respectively: The boundary conditions at z = 0 are shown in Equations (8) and (9): Finally, the image charges can be obtained from Equations (6)- (9): By substituting Equations (10) and (11) into Equation (6), one can obtain the potential ϕ 0 at the observation point A, as shown in Equation (12): The potential scalar function was used to derive the partial derivatives of x, y, and z. The value z = 0 was used to obtain the electric field component produced by the charged q of thunderstorm cloud B (x 0 , y 0 , z 0 ) at the ground observation point; that is, the installation position A (x, y, z) of the three-dimensional electric field meter: where → e x , Based on Equations (4) and (5) for the elevation and azimuth angles of thunderstorm clouds and on Equations (13)-(15) of the three-dimensional electric field components, the horizontal and elevation angles of the charge center B of thunderstorm clouds relative to the ground observation stations could be calculated from the three-dimensional atmospheric electric field component of the observation point (z = 0): Combined with the polarity of the electric field, the new three-dimensional electric field meter determined the horizontal and elevation angles of the thunderstorm cloud charge relative to the observation point by measuring the dielectric constant and electric field component of the installed surface in real time, so as to realize the thunderstorm cloud location in complex surface environments and different weather conditions.

Nonlinear Inversion Method for Thunderstorm Cloud Charge Parameters
The charge region in the thunderstorm cloud is mainly related to the structural characteristics of the temperature layer in which it is located [16]. The charge region in the cloud is evenly distributed in the horizontal direction, and is only related to the vertical height. It was assumed that the thunderstorm cloud was divided into several equally thick layers, and the charge of each layer was regarded as a point charge, as shown in Figure 3. The ground three-dimensional electric field value was the result of the interaction of multiple charge regions in the thunderstorm cloud. The ground three-dimensional electric field value was determined as in Equations (18)-(20): where (x, y, 0) and (x 0 , y 0 , z 0 ) were the coordinates of the ground three-dimensional electric field meter and the thunderstorm cloud, respectively. h was the thickness of the thunderstorm cloud charge, k was the index that indicated the charge zone layers, Q k was the charge of the k layer charge zone, and n was the total number of charge layers. Based on the functional relationship between the ground three-dimensional atmospheric electric field and the charge of the thunderstorm cloud (Equations (18)- (20)), combined with the classical thunderstorm cloud charge model, the parameters of the charge structure of the thunderstorm cloud-including the charge quantity Q k of each layer and the relative distance between the thunderstorm cloud and observation point |x − x 0 |, |y − y 0 |, z 0 could be inverted by using the ground three-dimensional electric field data.
Generally, the Newton iteration method is a common method for calculating non-linear equations. Each iteration direction is along the direction of the current point function value descending. However, the selection of the initial value of iteration is very important. The least square method is usually used to determine the initial value of iteration. The combination of the two algorithms can transform the non-linear problem into a linear one. Due to the existence of x − x 0 and y − y 0 in Equations (18) and (19) and their high dimensions, it is difficult to convert to the form of linear AX = B, and the initial values of Newton's method cannot be solved by the least square method. Therefore, other non-linear inversion algorithms were considered.
where ( , ,0) x y and 0 0 0 ( , , ) x y z were the coordinates of the ground three-dimensional electric field meter and the thunderstorm cloud, respectively. h was the thickness of the thunderstorm cloud charge, k was the index that indicated the charge zone layers, k Q was the charge of the k layer charge zone, and n was the total number of charge layers.  (18) and (19) and their high dimensions, it is difficult to convert to the form of linear AX = B, and the initial values of Newton's method cannot be solved by the least square method. Therefore, other non-linear inversion algorithms were considered. Genetic algorithm and particle swarm optimization are commonly used in geophysical non-linear inversion. In the process of solving nonlinear inversion, the genetic algorithm does not have a memory function. The solution obtained will change with the charge structure of thunderstorm cloud in the inversion process, and easily fall into local optimum, so it is difficult to obtain truly satisfied structural parameters. In the process of solving parameters of the three-dimensional electric field { , , } x y z E E E , the convergence speed of genetic algorithm is slow.
Particle swarm optimization (PSO) has memory function and can well preserve all the charge structure parameters with the optimal solution probability. It only updates the internal velocity of the particle and has fast convergence speed [17,18]. The problem of slow convergence of the Genetic algorithm and particle swarm optimization are commonly used in geophysical non-linear inversion. In the process of solving nonlinear inversion, the genetic algorithm does not have a memory function. The solution obtained will change with the charge structure of thunderstorm cloud in the inversion process, and easily fall into local optimum, so it is difficult to obtain truly satisfied structural parameters. In the process of solving parameters of the three-dimensional electric field E x , E y , E z , the convergence speed of genetic algorithm is slow. Particle swarm optimization (PSO) has memory function and can well preserve all the charge structure parameters with the optimal solution probability. It only updates the internal velocity of the particle and has fast convergence speed [17,18]. The problem of slow convergence of the three-dimensional inversion process can be well solved. Therefore, the particle swarm optimization (PSO) algorithm is suitable to inverse the charge structure of thunderstorm cloud.
The D-dimensional search space of the function randomly generates solid particles. These particles determine their movement direction by comparing the current best position and the historical best position of themselves or other particles. They approach the target by a constant updating cycle of their own speed and position, and finally find the optimal solution of the objective function as a whole [19]. The steps of using particle swarm optimization to retrieve the structural parameters of the thunderstorm cloud were as follows: (1) Initialization of the charge structure parameters In order to avoid blind searching of structural parameters in the inversion process, the velocity and position of structural parameters need to be initialized. Combined with the classical charged structural parameters of thunderstorm clouds, the inversion parameters are limited to a certain range [−V max , V max ] and [−X max , X max ], the initial velocity of charge parameter Q k should be [−5, 5], the relative distance between the thunderstorm clouds and three-dimensional atmospheric electric field meter should be assumed. The initial speed of the parameter (|x − x 0 |, |y − y 0 |, z 0 ) is [−100, 100]. The initial location of the charge parameter is [−30, 30], and the initial location of the relative distance is [−10,000, 10,000]. (2) Fitness function of the three-dimensional electric field The difference between the measured value E of the ground electric field and the inverted ground electric field value E * is defined as a fitness function. When value Y of the fitness function approaches zero, the optimal solution of the parameters can be obtained: (3) Find the initial extremum of the inversion parameter set Assuming the population size i of the charged structure parameter is 5 and the dimension of search space D is the number of inversion variables, and assuming the thunderstorm cloud is a classical three-layer charged structure, then D = 9. That is, the inversion parameter set is {Q 1 , Q 2 , . . . , Q k , Q 1 , Q 2 . . . , Q k , |x − x 0 |, |y − y 0 |, z 0 }. If it is a two-layer charged structure, then D = 5, that is, the set of inversion parameters is {Q 1 , Q 2 , . . . , Q k , |x − x 0 |, |y − y 0 |, z 0 }. The initial extremum P best of each parameter of the charged structure is the position of the individual when the initial fitness is optimal, and the global initial extremum G best is the position with the best fitness searched in the i parameter sets. Let the initial velocity of the i structural parameters be  (24) and (25): where V k id and V k+1 id are the convergence speeds of the kth and k + 1th structural parameters, respectively, and X k id and X k+1 id are the structural parameter values of the k and k + 1 times, respectively. c 1 and c 2 are acceleration factors that determine the ability of particles to search for global optimal values and individual optimal values. d is the number of inverse evolutions, w is the weight coefficient, and P k id is the ith particle at the kth iteration. P k gd is the global extremum of the ith structural parameter set at the kth iteration.

Three-Dimensional Atmospheric Electric Field and Radar Observation Experiments
More than one year of observation experiments were carried out using our three-dimensional atmospheric electric field meter [20]. Figure 4 is a three-dimensional electric field curve for the period from 14:35 p.m. to 17:15 p.m. on 9 October 2018. It can be seen from the figure that the change of the ground electric field caused by the thunderstorm cloud can be divided into four stages. In the first stage (14:35 p.m.-15:25 p.m.), the atmospheric electric field E z on the ground appeared to fluctuate continuously and the amplitude began to increase. The charge carried by the thunderstorm cloud entered the detection range of the three-dimensional atmospheric electric field meter, but was also far away. In the second stage (15:30 p.m.-15:50 p.m.), due to the discharge phenomenon, the atmospheric electric field E z oscillated greatly, the polarity was reversed multiple times, and the charging and discharging process was obvious. The thunderstorm cloud is located directly above the observation point. In the third stage (15:50 p.m.-16:25 p.m.), the atmospheric electric field E z gradually weakens, and the horizontal electric field and electric field amplitude reach a peak multiple times. It was judged that the thunderstorm cloud was located in the east direction of the observation point and the distance was relatively close. In the fourth stage, after 16:25 p.m., E x , E y , and E z gradually decreased and tended to be stable, and the thunderstorm cloud gradually moved away from the detection range. gradually weakens, and the horizontal electric field and electric field amplitude reach a peak multiple times. It was judged that the thunderstorm cloud was located in the east direction of the observation point and the distance was relatively close. In the fourth stage, after 16:25 p.m., x E , y E , and z E gradually decreased and tended to be stable, and the thunderstorm cloud gradually moved away from the detection range. At the same place, the radar map of the same location from 15:06 p.m. to 16:12 p.m. are shown in Figure 5. The central point of the circle was the observation point of the three-dimensional atmospheric electric field meter A. It can be clearly seen that at 15:06 p.m., the thunderstorm cloud was far away from the observation point, and the amplitudes of E x , E y and E z were smaller. At 15:30 p.m., the thunderstorm cloud approaches the observation point A gradually, and the E z component increases. At 15:48 p.m., the thunderstorm cloud is located directly above the three-dimensional atmospheric electric field meter, and the E z component measured is strong. At this time, thunderstorm occurrence is predicted. At 16:12 p.m., the thunderstorm cloud over the region moves southeast, still within the detection range, and E x and E y increase gradually. The movement of thunderstorm clouds in this period is consistent with that measured by three-dimensional atmospheric electric field meter.

Inversion of the Charge Model
According to the observation data of radar and atmospheric electric field meters, the classical three-layer charging model and the structural inversion of the two-layer charging model were respectively performed on the thunderstorm cloud at a certain time point of a second stage thunderstorm. The sampling interval of the three-dimensional atmospheric electric field meter was 1 s, and the measured data of the three-dimensional electric field at two adjacent sampling moments were 7.299 kV/m, 2.219 kV/m, −29.560 kV/m and 7.222 kV/m, 2.183 kV/m, −29.207 kV/m. Due to the extremely short time, the observation point was considered. The relative distance (|x − x 0 |, |y − y 0 |, z 0 ) from the thunderstorm cloud was approximately constant. A set of nonlinear equations between the three-dimensional electric field value and each charge structure parameter was established as per Equations (26)-(31): According to the inversion method in the previous section, the parameters of the three-layer charged model of thunderstorm clouds were inverted by using the three-dimensional atmospheric electric field measurements. The convergence curve is shown in Figure 6. The inversion results are shown in Table 1, where Q 1 , Q 2 , Q 3 , Q 1 , Q 2 , and Q 3 are the charges of the thunderstorm clouds, |x − x 0 |, |y − y 0 |, and z 0 are the relative distances between the observation points of the ground three-dimensional electric field meter and the thunderstorm clouds. It can be seen from the inversion results that the charged structure of the thunderstorm cloud had a positive charge on the uppermost layer, a negative charge on the middle layer and the bottom layer, and (Q 1 , Q 2 , Q 3 ) were approximately equal to and (Q 1 , Q 2 , Q 3 ), respectively, which also indicated a very short time. The charge of the internal thunderstorm cloud did not change substantially, and the charge structure was an even-polar structure with a negative polarity, which was consistent with the classical dipole structure.    It can be seen from the inversion results that the charged structure of the thunderstorm cloud had a positive charge on the uppermost layer, a negative charge on the middle layer and the bottom layer, and ( 1 ' Q , 2 ' Q , 3 ' Q ) were approximately equal to and ( 1 Q , 2 Q , 3 Q ), respectively, which also indicated a very short time. The charge of the internal thunderstorm cloud did not change substantially, and the charge structure was an even-polar structure with a negative polarity, which was consistent with the classical dipole structure.

Comparative Analysis of Electric Field Data
According to the lightning activity in Nanjing and surrounding areas, Nanjing Meteorological Bureau has established a monitoring network covering the whole city by using the traditional atmospheric electric field meter, as shown in Figure 7. Point A in the figure is the three-dimensional electric field meter (AEFM), which is located at Nanjing University of Information Science and Technology (118°43′30.00″ E, 32°11′56.40″ N). Pukou, Liuhe, Nanjing meteorological Bureau, and Kunlun are the observation stations of several conventional atmospheric electric field meters that are closest to point A. The positive direction of X is the east direction, and the positive direction of Y is the north direction.

Comparative Analysis of Electric Field Data
According to the lightning activity in Nanjing and surrounding areas, Nanjing Meteorological Bureau has established a monitoring network covering the whole city by using the traditional atmospheric electric field meter, as shown in Figure 7. Point A in the figure is the three-dimensional electric field meter (AEFM), which is located at Nanjing University of Information Science and Technology (118 • 43 30.00 E, 32 • 11 56.40 N). Pukou, Liuhe, Nanjing meteorological Bureau, and Kunlun are the observation stations of several conventional atmospheric electric field meters that are closest to point A. The positive direction of X is the east direction, and the positive direction of Y is the north direction. The distance between each station and the three-dimensional electric field meter is known in Figure 7. So the vertical electric field components of each station can be calculated based on the structural parameters of thunderstorm cloud retrieved in Table 1. The vertical electric field values obtained by the three-dimensional electric field meter were compared with the official measurements as shown in Table 2.
It can be seen that the vertical electric field values at the other locations obtained by the three-dimensional electric field meter were consistent with the official measurements. The distance between each station and the three-dimensional electric field meter is known in Figure 7. So the vertical electric field components of each station can be calculated based on the structural parameters of thunderstorm cloud retrieved in Table 1. The vertical electric field values obtained by the three-dimensional electric field meter were compared with the official measurements as shown in Table 2. It can be seen that the vertical electric field values at the other locations obtained by the three-dimensional electric field meter were consistent with the official measurements.

Discussion and Conclusions
The three-dimensional electric field meter we developed overcomes the shortcomings of a traditional electric field meter that can only measure the vertical electric field components at the ground. Observational experiments have confirmed that thunderstorm clouds produce horizontal electric field components that cannot be ignored near the ground. This provides more comprehensive electric field information for the tracking, warning, and charge structure of thunderstorm clouds. By establishing an observation model of the three-dimensional electric field, it was proven that the orientation information of the thunderstorm cloud can be determined by the ratio of the three-dimensional electric field values of the observation points. Therefore, the three-dimensional electric field data can be used to continuously monitor the movement trajectory of the thunderstorm cloud, which can be used as a useful supplement to radar data.
Lightning warning is an important application of atmospheric electric field meters. The traditional atmospheric electric field meter has a certain leakage alarm because it only measures the vertical atmospheric electric field and ignores the horizontal component. For example, when the three-dimensional electric field meter in the experiment measured the thunderstorm cloud in the second stage (15:50 p.m.-16:20 p.m.), the horizontal component electric fields E x and E y were obviously stronger than the vertical electric field E z . The amplitude gave a lightning warning. However, if it was measured by a conventional electric field meter, it will be impossible to give an early warning because the electric field E z in the vertical direction is weak.
The amount of electric field information measured by traditional electric field meters is too small, and it is difficult to invert the structural parameters of the thunderstorm cloud. Multiple networking electric field meters must be used to calculate the positional relationship of each electric field meter and the measured vertical electric field. However, the installation environment of each electric field meter in the grid including the terrain, surrounding buildings, etc., will cause distortion of the atmospheric electric field, resulting in poor consistency of the measurement data, and it is difficult to correct. Based on incorrect electric field data, the inversion of atmospheric electric fields is difficult to carry out. The three-dimensional electric field and the charged structure parameters are nonlinear, and the particle swarm optimization algorithm is suitable for nonlinear inversion problems. However, in the process of finding the optimal solution in the inversion iteration, the position of the particle must be limited according to the typical thunderstorm cloud charge structure parameter. The inversion method assumes that the thunderstorm cloud is a multi-layer structure, and the charge center of the thunderstorm cloud is treated as a point charge. The thunderstorm cloud-charged structure of the experimental inversion observation site (Nanjing) was evenly polar, which was consistent with the results of the comprehensive analysis of the electric field and radar data in the literature [21].
The measurement of three-dimensional atmospheric electric field provides more comprehensive data for the study of lightning structure and lightning warning of thunderstorm clouds. In addition to ground detection, the developed three-dimensional electric field sensor can also be used in the field of high-altitude detection to measure the horizontal electric field component that cannot be ignored, which is impossible for current vertical electric field meters. It should be noted that the proposed inversion method was based on the thunderstorm cloud charge as a point charge model and was based on a classical uniform multilayer structure. Some studies believe that thunderstorm cloud-charge structures are geographically related and further work is recommended to validate this process.