1. Introduction
The article (Lozbin, Shpadi, & Inchin, 2016) proposes an algorithm for calculating the location and time of a lightning discharge using the Time-of-Arrival (ToA) method for three stations. However, this algorithm is not sufficiently mathematically substantiated. This shortcoming is addressed in this article. In addition, this article proposes a new algorithm for solving the problem, which is simpler for numerical calculations than the previous one (Lozbin, Shpadi, & Inchin, 2016).
The mathematical model of the problem under consideration is implemented on the unit sphere. In this regard, we note some provisions of the geometry on the sphere, which we will use (Krantz, 2007; Volynsky, 1977).
A sphere, its radius is one which is called a unit sphere.
A great circle of a sphere is a circle cut by a sphere from a plane passing through its center.
A spherical segment is an arc of a circle of a great circle. The length of the spherical segment of the unit sphere is measured by a central angle, the sides of which pass through the ends of the spherical segment.
The distance between two points of a sphere is the length of a smaller spherical segment connecting these points.
A spherical triangle is a part of a sphere bounded by three spherical segments with common vertices.
A spherical triangle, the length of each side of which is strictly less than the length of the semicircle of the great circle, is called Euler’s triangle. For an Euler triangle, the following statement is true: the sum of the lengths of its two sides is greater than the length of the third side.
A spherical triangle whose area is zero is called a degenerate spherical triangle. All of its vertices are located on the same circumference of a great circle and can coincide. In this article, a degenerate spherical triangle appears as a special case in problems with parameters that determine the values of its sides and angles.
2. Formulation of the Problem
Consider a unit sphere with a center coinciding with the origin of the Cartesian system Oxyz. We introduce a coordinate system on the sphere, where and are longitude and latitude of the points of the sphere respectively, , . In accordance with the concepts of a geographic coordinate system, points with latitude and we will be called the North and South poles, and the arcs connecting them as meridians. In addition, let the equator be points of a circle on a sphere with latitude .
Consider Euler’s triangle, the vertices of which are at the points , . Let be an unknown point on the sphere, the position of which is to be found. Let , , be spherical segments of a given length going out from points in the direction of an unknown point . The length of each segment is also denoted by . Such a combination of concepts will do not cause a semantic contradiction.
We denote as the increment of the lengths for all segments . This increment , which is the same for all segments , can be a positive or negative number, as well as zero.
The considered problem is a determination on the sphere the point and the value of increment of the segments in accordance with the requirement: the spherical distance from the vertices to the point must be equal to the sum , .
3. Task Analysis
To solve the problem, we will use vector algebra (Beklemishev, 2008). The one-to-one correspondence between Cartesian and spherical coordinates of points of a unit sphere is determined by the following equations:
(1)
We introduce unit radius vectors of the vertices of triangle and unit radius vector of the target point. The criterion for the non-degeneracy of the triangle is the non-zero value of the triple product of three vectors . Therefore the condition
(2)
must be met, where the symbol “ “ and “ “ are scalar and vector products of vectors.
Since the triangle is by definition Eulerian, the following inequalities hold:
, , , (3)
where determines the angle between unit vectors and , and distance between points and equal to this angle.
Thus, the task is reduced to solving the system of equations:
(4)
with three unknown , and .
Next, we assume that the numbering of the vertices of the triangle is selected taking into account inequalities:
. (5)
Theorem 1. In order for the system of Equations (4) to have some solution in the conditions of the above problem, it is necessary to fulfill the next conditions:
(6)
Proof. First, we note that condition (5) implies inequalities
(7)
Suppose that as a result of the solution of the system (4), the position of the point is determined. In this case, for each possibly degenerate, triangle the inequalities must be satisfied:
(8)
Substituting in the right part of inequalities (8) the values of arc cosines from (4), we get a system of inequalities (6). ■
Theorem 2. For solvability of the system of Equations (4), the lengths of the segments must satisfy the conditions
(9)
Proof. The non-negativity of the difference in (9) follows from inequalities (5). In case of , the corresponding inequality in (6) will contradict conditions (3). ■
4. Solution Algorithm
We consider the algorithm of solution of system (4), assuming that the segments are satisfied the inequalities (5) and (9), and the inequalities (6) are satisfied in a strict sense, that is:
(10)
Without losing generality, we will assume that the point is at the North pole, that is . In this case, the radius vector of the point is
. (11)
If the point does not coincide with the North pole, then we should move the triangle around the sphere by orthogonal transformation:
, (12)
which must be applied to radius vectors . The matrices and are defined by the longitude and latitude of the point and are written as:
, . (13)
The matrix rotates the triangle points about the axis by an angle , and the matrix rotates about the axis by an angle . As the result of these rotations, the triangle moves along the surface of the sphere, with preserving the lengths of the sides and the values of the angles, as well as the direction of the border so that the point is at the North pole.
After solving the system (4) and finding the point , the inverse transformation
, (14)
applied to the radius vector will restore the location of the point relative to the original position of the spherical triangle . The symbol “T” in (14) means the transposition of matrices.
Subtracting the first equation of system (4) from the second and third, we obtain a system of two equations
. (15)
in which only two spherical coordinates and of the points are unknown.
Taking to account (11), we obtain:
, . (16)
The system (15) taking into account (16) be as
(17)
Calculating the cosines of the left and right parts of the equations in (17), we obtain:
(18)
Let us expand the scalar products of vectors in (18)
,
and let us convert the right parts in (18) by formula of cosine of two angles sum
.
We write the equations of system (18) in the form of equalities to zero of two linear combinations of functions and
(19)
Since the functions and do not simultaneously vanish for any values of due to the identity , then for the existence of a solution to system (19) its determinant must be equal to zero, i.e.,
. (20)
Calculating the determinant in (20) and selecting in it a linear combination of functions and , we obtain the equation for :
(21)
Let us introduce the following notation
, (22)
, (23)
, (24)
and write an Equation (21) as
. (25)
We show, that
. (26)
Suppose that the opposite is true, that is, that the equality
. (27)
From the conditions (10) and equations
, (28)
we get that
. (29)
If we now consider Equations (27) taking into account (22)-(23) as a linear homogeneous system of equations in which expressions , , are unknown, then in view of (29) its determinant must be zero, that is,
.
This equality means that the projection of the vector product on the axis is zero. It follows that the plane of the vectors and passes through both poles of the sphere. As a result, we get that all three vectors , and belong to the same plane, and this contradicts the condition (2) of non-degeneracy of the triangle . Thus, assumption (27) is not true and inequality (26) under conditions (10) be valid for any triangle .
By defining a value of parameter from a system of equations
(30)
and assume, that
. (31)
reduce Equation (25) to form
. (32)
In case of , the general solution of the Equation (32) will be as
, (33)
from which we need to select values which are in the range . In case of , there will be two such solutions and, in the case of only one solution. If , then the problem has no solutions.
After finding longitude , the corresponding latitude value of the target point , belonging to the interval is determined from any equation of the system (19), that is,
, where or . (34)
From (29) it follows that , . In this case, Equation (34) with unknown has only one solution on the interval .
For calculating directly, it is convenient to use equivalent formulas
, , (35)
since in computing environments for computers the arc cotangent function is usually absent.
The value is determined from any equation of the system (4), that is,
, . (36)
Theorem 3. The inequality
(37)
is always satisfied, when , and are defined by expressions (22)-(24) and conditions (5), (6) and (9) are met.
Proof. The inequality (37) is proved by directly calculating its left part using the expressions (22)-(24) and transforming them identical. We have
(38)
Applying identities
,
,
,
and taking into account (38), we have
Then, after a series of identical transformations taking into account (24), we have the final result
(39)
Inequalities (6) under conditions (9) are equivalent to inequalities
from which the inequality (37) is coming. ■
Corollary. System of Equations (15) always has two solutions if, under the conditions of Theorem 3, inequalities (10) are satisfied in the strict sense. Indeed, in case of satisfaction (10), the right part of (39) is positive and, hence, there is a strict inequality, from which, taking into account (31), the is follows.
5. Exceptional Cases of Solving the Problem
Those cases of solution of Equation (25) when the equalities
, (40)
were called exceptional. Equality (40) occurs if at least one of the equalities
(41)
is fulfilled.
Theorem 4. If, with some i and j, , the equality (41) is fulfilled, then the vector lies in a plane of vectors and .
Proof. First, we emphasize that due to the non-degeneracy of the triangle, the vectors and , , , are not collinear. Therefore,
. (42)
In this case, each pair of vectors and uniquely defines a certain plane. Note also that, a strict inequality , , follows from (5), (41) and (42).
Calculating the cosines of the left and right sides in (41), we obtain . Subtracting the i-th equation of system (4) from the j-th, , and replacing the difference with , we obtain the equality
. (43)
Consider a trihedral angle with edges , and . Based on the theorem from stereometry (Kiselev, 2004) on the properties of the sum of two plane angles of a trihedral angle, from equality (43) we obtain that the trihedral angle with edges , and is degenerate and, therefore, the vector lies in the plane of the vectors and .
Consider the proof of Theorem 4 based on vector algebra and find a formula for calculating the vector . We introduce the notations , and , from which we find
(44)
From (43) and (44) we obtain
, . (45)
We expand the vector in terms of the system of non-coplanar unit vectors
, (46)
where
, .
Let us calculate the coefficients in (46). Multiplying both sides of equality (46) scalarly by vectors , and , also taking into account their unity, we obtain the system of equations
with unknowns . Solving this system, we find,
, ,
.
Simplifying the right-hand side of the last equality, we get .
Thus,
, (47)
that is, is a linear combination of vectors and and, therefore, lies with them in the same plane. ■
Let us consider the case when in (39) only one factor of the right-hand side vanishes. The solution of the problem is unique and can be calculated using the Formulas (22)-(24), (30), (31), (33) and (35). In this case, in Formula (35), you need to choose a denominator that is not equal to zero.
If two factors vanish on the right-hand side of (39), then the solution to the problem is the common vector of these factors. In this case, the desired point coincides with the vertex of the triangle, which is the end of this vector.
The third case, when all three factors in the right part of (39) become zero, cannot be realized, since otherwise the triangle will be degenerate.
6. The Problem of Lightning Discharge Coordinates and Time Finding
One of the current actual problems is monitoring of cloud-to-ground lightning discharges (Hu, Zhao, & Zhu, 2010; Koshak & Solakievicz, 2001; Kuterin, Bulatov, & Shlyugaev, 2014; Ramachandran, Kumar, & Kishore, 2003; Shpadi, Inchin, Lozbin, & Ayazbayev, 2019). To solve it, the technology of deploying a network of stations registering a low-frequency radio wave coming from lightning discharges, is often used.
If a station located at a point , , , on the surface of the Earth recorded at the moment a radio signal from a lightning discharge that occurred at an instant at a point , then this event leads to the equations
, (48)
in which and are the unit radius vectors of the i-th station and the points of the lightning discharge, is the average radius of the Earth, is the speed of light. Introducing the notation
, , (49)
we obtain one of the equations of system (4).
If the number N of stations in the network is more than three, then they are grouped by three stations and the calculation is performed for each group. The calculation results are summarized in order to determine the optimal location and moment of the lightning discharge, taking into account the error in measuring the time and coordinates of the stations.
The method based on the idea of recording the moment of receipt of a radio signal is called Time-of-Arrival (ToA or TOA) in the scientific literature.
Let’s consider a test case illustrating calculation of parameters of a lightning discharge for a network of three stations. Table 1 shows the location of stations in Kazakhstan and the arrival time of lightning radio signals.
We assume that the lightning strike occurred at a point (near the city of Aktobe, Kazakhstan) at moment . The latitude of the stations in Table 1 and the latitude of the test point are indicated in the range , that is, in accordance with geographic norms. The latitude of these points in the calculation is converted to the range . The calculation of the fixation time , , of lightning radio signals was carried out using equality (48). Mathematically exact values of this calculation are rounded to the practically achievable value of 10−8 s. Also, Table 1 shows the values of the parameters , which are calculated using Formulas (49).
Table 1. The initial data of the problem.
Using the data in Table 1 and applying Formulas (1), we obtain:
, (50)
, (51)
Based on expressions (50) and (51), we conclude that conditions for solvability of problem are satisfied and the problem has two solutions.
The transformation of stations coordinates is performed by the matrix (12)
. (52)
The results of the transformation of coordinates are shown in Table 2. The station near the city of Taraz has been moved to the North Pole.
Further, using Formulas (22)-(24), (30), (31), (33) and (35), we successively obtain:
, (53)
, , (54)
, (55)
. (56)
Finally, after performing the inverse transformation of coordinates of desired points and transforming their latitude to the range , we find solution to the problem. The calculation results are shown in Table 3.
Table 2. Stations location and parameters.
Table 3. Calculation results of the problem.
The first lightning in Table 3 corresponds to the test , and the second point is optional. The moments of discharges were calculated using right-hand Formula (49).
All calculations are performed in the MATLAB environment.
7. Conclusion
1) For the geometric problem considered in the article, the necessary and sufficient conditions for it are established. Subject to their observance, this task usually has two solutions, in exceptional cases, only one.
2) The above problem can be used as a mathematical model in solving the problem of time and place of cloud-to-ground type of lightning discharge finding.
3) The values of the parameters can be arbitrary, taking into account that the solution of the problem depends on their differences. It follows that the time reference point in the Time-of-Arrival method can be chosen arbitrarily, but, naturally, it should be the same for the entire network of stations.
Acknowledgements
This article was created in the frame of the scientific and technical program BR05336383 of the Ministry of Digital Development, Innovations and Aerospace Industry of the Republic of Kazakhstan.
[1] Beklemishev, D. V. (2008). Course of Analytical Geometry and Linear Algebra. Moscow: Fizmatlit. (In Russian)
[2] Hu, Z. X., Zhao, W. G., & Zhu, H. P. (2010). A Source Location Algorithm of Lightning Detection Networks in China. Annales Geophysicae, 28, 1981-1991.
https://doi.org/10.5194/angeo-28-1981-2010
[3] Kiselev, A. P. (2004). Geometry. Moscow: Fizmatlit. (In Russian)
[4] Koshak, W. J., & Solakiewicz, R. J. (2001). TOA Lightning Location Retrieval on Spherical and Oblate Spheroidal Earth Geometries. Journal of Atmospheric and Oceanic Technology, 18, 187-199.
https://doi.org/10.1175/1520-0426(2001)018<0187:TLLROS>2.0.CO;2
[5] Krantz, P. (2007). Spherical Trigonometry. Moscow: LKI Publishing House. (In Russian, Translation from German)
[6] Kuterin, F. A., Bulatov, A. A., & Shlyugaev, Yu. V. (2014). The Development of the Lightning Detection Network Based on BoltekStormTracker Hardware. XV International Conference on Atmospheric Electricity, Norman, 15-20 June 2014, 71-74.
[7] Lozbin, A., Shpadi, Yu., & Inchin, A. (2016). Triangles Technique for Time and Location Finding of the Lightning Discharge in Spherical Model of the Earth, Journal of Geoscience and Environment Protection, 4, 125-135.
https://doi.org/10.4236/gep.2016.44016
[8] Ramachandran, V., Kumar, S., & Kishore, A. (2003). Cloud-to-Ground Lightning Location Using TOGA of Sferics. The South Pacific Journal of Natural and Applied Sciences, 21, 63-65.
https://doi.org/10.1071/SP03013
[9] Shpadi, Yu. R., Inchin, A. S., Lozbin, A. Yu., & Ayazbayev, G. M. (2019). The Method of Differentiable Functional in a Problem of Lightning Discharge Localization and Time Moment Determination. Journal of Physics: Conference Series, 1352, 012049.
https://doi.org/10.1088/1742-6596/1352/1/012049
[10] Volynsky, B. A. (1977). Spherical Trigonometry. Moscow: Nauka. (In Russian)