The Effect of the Earth's Oblate Spheroid Shape on the Accuracy of a Time-Of-Arrival
Total Page:16
File Type:pdf, Size:1020Kb
.-. THE EFFECT OF THE EARTH'S OBLATE SPHEROID SHAPE ON THE ACCURACY OF A TIME-OF-ARRIVAL. LIGHTNING GROUND STRIKE?LOCATING SYSTEM Paul W.Casper, PE and Rodney B. Bent, Ph.D. Atmospheric Research Systems, Inc. 2350 Commerce Park Drive, NE,Suite 3 Palm Bay, Florida 32905 U.S.A. ABSTRACT The algorithm used in previous technology time-d-arrival lightning mapping systems was based on the assumption that the earth is a perfect spheroid. 'I'hese systems yield highly-accurate lightning locations, which is their major strength. However, extensive analysis of tower strike data has revealed occasionally significant (one to two kilometers) systematic offset errors which are not explained by the usual error sources. It has been determined that these systematic errors reduce dramatically (in some cases) when the oblate shape of the earth is accounted for. The oblate spheroid correction algorithm and a case example is presented in this paper. INTRODUCTION it becomes essential to base all mathematics on an extremely accurate earth model. Lightning ground strike tracking systems based on the time-of-arrival (TOA) technique, in Tracking systems using terrestrial timing syn- combination with wideband waveform detection chronization signals (e.g., LORAN) are doubly have been in operation for almost a decade. The affected by the earth's oblate shape. The follow- accuracy of these systems has steadily improved ing section identifies the affected parameters and as more has been learned about the fine char- discusses the methods for improvement. Subse- acte ristics of timing synchronization signals, quent sections contain illustrations of the magni- propagation effects on lightning waveforms and tude of each effect. software processing methods for accuracy en- hancement. One of the more recent improve- MATHEMATICAL PROCEDURE ments has been accomplished by refinement of the mathematics to more accurately accommo- The oblate shape of the earth spheroid directly date the oblate shape of the earth spheroid. affects TOA system accuracy in two ways: Earlier versions of TOA lightning tracking sys- a. Calculation of time clock offsets due to timing tems used mathematics based on a spherical reference signal differential propagation delays approximation of the earth. This is usually not seriously in error, especially if the accurate b. Calculation of lightning stroke coordinates earth's radius for the region of interest is used in given a set of accurate receiver time differences the approximation. However, when accuracies otherwise are approaching a few hundred meters, Different mathematical processes are involved in each of these steps. Derivation of the equations 81-1 is beyond the scope of this paper (consult refer- figure provided by the U.S. Coast Guard for the ences [l,2, 3]), but the methodology for each is LORAN navigational system. This 100 KHz described in the following two subsections. velocity figure is determined as follows: CALCULATION OF TIME CLOCK OFFSETS V' = Vln (1) Receiver sites and a terrestrial timing reference where: V' = 100 KHz ground wave velocity signal transmitter are at fixed locations. There- fore, the clock offsets due to timing signal dif- n = index of refraction at the earth's ferential propagation delays need only be calcu- surface for 100 KHz (1.000338) lated one time and stored in the system software. In order to determine these differential offsets, V = free space velocity (299,792,458 very accurate geodesic distances between the me terslsecond) timing transmitter and each LPATS receiver must be first computed. A geodesic is defined as Accurate time clock offsets are then easily calcu- the curve of minimum length between two points lated: on the surface of a spheroid [l]. T12 = V' (D1 - D2) (2) Geodesic distances can be estimated fairly accu- rately by using a perfect sphere model of the where: D1 = geodesic distance from receiver 1 earth with a radius equivalent to the earth radius to the timing transmitter at the mean latitude of the network. However, data to be presented in later sections shows that D2 = geodesic distance from receiver 2 to achieve systematic location errors of less than to the timing transmitter several hundred meters, the actual shape of the earth must be properly accounted for. Equation (2) produces the time offset between clocks in receivers 1 and 2. This figure for each The non-iterative solution by Sodano [l] pro- clock pair in the system is stored by the central vides a highly-accurate means of calculating software and used to correct the time differences geodesic distances. Sodano developed a system actually reported by the receivers. of equations that are easily programmed and solved in double precision using a personal com- CALCULATION OF STROKE puter equipped with a mathematic co-processor. COORDINATES FROM TIME DIFFER- The input to this set of equations is only the lati- ENCES tude and longitude (accurate to the fourth deci- mal place) of the two points of interest. Naviga- The second part of the problem is computation tion receivers using the satellite-based Global of the stroke coordinates given accurate input Positioning System are most convenient for time differences. Razin [2] describes the math- determining the coordinates of receivers to the ematics for accomplishing this computation for required accuracy. both the spherical earth approximation and for the more exact oblate spheroid case. Once accurate geodesic distances are available, Razin does not present the equations for the they must be converted to an equivalent propaga- oblate case. Fell [3] fills in this void. The full set tion time. For this, we need an accurate ground of equations is quite complex and extensive, and wave propagation velocity figure, applicable for the reader is referred to both Razin and Fell for the high-energy frequencies radiated during the the details. The computational procedure will be first ten microseconds of a lightning ground described here to provide general understanding stroke. A figure for frequencies in the 50 to 300 of the process. KHz range is appropriate, such as the 100 KHz 81-2 The first step in the process is to map the re- actual tower strikes. Although the cluster pattern ceiver coordinates from the spheroid (earth) is identical to the tower pattern, there is an onto an osculating sphere, internal to the apparent systematic offset to the southwest. Af- spheroid, which is tangent to the spheroid at a ter investigation of all the usual sources of point Po. This tangent point is selected near the systematic error (receiver coordinates, computa- center of the receiver network. Figure 1 illus- tional error, etc.), it was determined that the er- trates this arrangement in two dimensions for ror was substantially due to time clock propaga- clarity. Once all the receiver coordinates (PR) tion offsets and the spherical approximation are mapped onto the osculating sphere (PB), the stroke coordinate solution mathematics. Here solution process proceeds as if the earth is a per- we will examine the error due to time clock off- fect sphere, using the Pk set of receiver coordi- sets and address mathematical error in the next nates. When the stroke coordinates on the section. osculating sphere are computed (PIs), they are mapped back to the sphtvoid (Ps) using a very I:igure 4 is a blow-up of the area around towers simple equation pair. #1 and #2. A 200m x 200m grid has been su- perimposed to aid in scaling distance. The cen- Note that the osculating sphere receiver coordi- troid of the cluster near tower #1 is seen to be nates need only be computed one time. They can approximately 590 meters southwest of the then be stored as fixed constants in the LPATS tower. The same is true of the #2 tower and central computer software and used with the re- stroke cluster. The circles around the tower indi- ported time differences to calculate each stroke's cate the approximate attractive radius of each coordinates. The burden on the real time central tower due to its height. It was determined from computer software is substantially unchanged the data base that receiver triad 1-4-5 was used to from the all-spherical case. locate virtually all the strokes in the Bithlo tower clusters, which is in fact the triad the system ERROR DUE TO TIME CLOCK OFFSETS should have used since it is the optimum combination (least affected by timing errors for Offsets in the receiver time clocks, because of the Bithlo area). We will concentrate on this re- spherical approximation error in the offset cor- ceiver triad to examine the spherical approxima- rection constants, can be significant. Here we tion time clock offset error. will examine an actual case to illustrate the point. The centroid of the actual tower #1 cluster is The five Florida LPATS network receiver loca- shown in Figure 5 (point A). Point B is a theo- tions are shown in Figure 2, plus the Jupiter, retical point, computed as follows: Florida, LORAN transmitter that is used for synchronization. This network has been in op- a. A stroke location at the exact coordinates of eration for several years, and a substantial the tower was assumed. archive of data has accumulated. In an effort to objectively assess the accuraj of this network, b. The exact propagation times to receivers 1,4 the data base was searched over a complete light- and 5 were computed using the Sodano ning season in the vicinity of three very tall, method. attractive objects ("Bithlo towers") to determine if an expected pattern of strikes existed. Each of c. nerelative stroke time differences seen by these towers are over 1,000 feet in height. The each receiver pair were then computed using location of the three Bithlo towers relative to the the propagation times.