About Resolvability of the Problem of Lightning Coordinates Finding by the Time of Lightning Discharge Arrival ()
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.