Differential GPS Measurements for Gravity Survey

Introduction

A primary science activity for the Survey Fire was to provide accurate height information for a Gravity Meter Survey to be conducted by the Geosciences Fire.  In the event, both Fires worked together on this project.

The most important consideration was the required accuracy of height measurements. The Gravity meter measures the pull of gravity and the last count on the dial is equivalent to a mere 0.03 metres in height, although accurately reading this is extremely tricky, so we needed to achieve relative heights between successive gravity reading of better than 0.1 metres with 0.03 metres as the goal. 

This level of height accuracy is achievable using the more traditional methods of Spirit Levelling or Trigonometrical Heighting with a Theodolite but is also achievable using the more modern method of GPS (Global Positioning System) using Survey GPS rather than Navigation GPS.  Both Spirit Levelling and Trigonometrical Heighting are very time consuming when compared with GPS measurements and would conflict seriously with the logistics of the expedition, so GPS was chosen.  In order to achieve required height accuracy the technique of DGPS (Differential GPS). 

The difference between Navigation and Survey GPS

Navigation GPS receivers use a Code transmitted by the GPS satellites to calculate the position on the earths surface.  The precision of the position fix using Navigation GPS is to tens of meters as detailed in a subsidiary report.  In order to achieve centimetre accuracy, Survey GPS receivers use the signal carrier wave from the GPS satellites.  This does not provide a position fix, but if two receivers are collecting data simultaneously, at each end of a line (base line), the vector (differences) between the two receivers can be calculated to this high level of precision.  The assumption is that both receivers have the same errors caused by the signals passing through the Ionosphere and the Atmosphere, especially if they are relatively close together (i.e. 10 to 20 Kilometres apart).  Absolute position can be determined using Survey GPS by relating the data back to a base line with one receiver on a “known” point. 

Methodology

DGPS involves the simultaneous collection data by two or more Survey GPS receivers with the “corrections” being related to one of the receives, called the base station.  This can be done in “real time” using a technique known as RTK (Real Time Kinematic) GPS where the base station transmits the “corrections” to the other receiver(s) by radio or GSM phone, or the data can be collected and post processed.  The advantage of working with RTK is that as soon as sufficient data has been collected to resolve the ambiguities (see the end of this report) the survey team can move on, but the disadvantage is that additional equipment (e.g. radios, batteries, etc.) is required in the field adding to logistical considerations, especially with battery charging and radio reception.  The disadvantage of collecting data for post processing is that the GPS receivers have to be recording data simultaneously for enough time for the users to be sure that there is sufficient data to resolve the ambiguities at the processing stage.  This method fitted well with the expedition logistics as the gravity meter readings themselves required a similar time to that needed to be sure that sufficient GPS data had been collected.  The general rule of thumb is two minute for each Kilometre between the receivers and the base line lengths were envisaged to be around 10 Km making occupation time of around 20 minutes, which worked in well with the gravity meter reading times. 

Simultaneous collection of GPS and Gravity Meter data.The GPS and Gravity Meter data is collected simultaneously.

Equipment

The choice of equipment was two Leica SR530 receivers with lightweight aluminium tripods.  These Leica GPS SR530 receivers proved their worth in this harsh environment and worked faultlessly throughout the expedition despite the heat and dust.  The Young Explores were   inexperienced in the use of Survey GPS equipment, having not seen any until they actually arrived in Tanzania, but were able to operate it and achieve excellent results with tuition on-site. 

The data was processed using the LGO (Leica Geo Office) software.

Absolute Position on the Earths Surface

The DGPS data provided accurate relative positions for all the points measured during the expedition.  Although these positions could be related to their location on the earths surface by using the Navigated positions determined as part of the data collection process a more accurate result was desired.  These could be achieved by computing the data collected during the expedition in conjunction with data collected simultaneously at ISG Earth Stations.  These ISG Earth Stations collect GPS data at points on the earth’s surface with locations determined to great precision by a variety of techniques including GPS.  The data can be accessed through the Internet and used to process GPS base lines in the same way as the used to calculate the GPS base lines measured during the expedition.

Four primary locations were selected at the camp sites at Longido, Kitumbeine, Lengai Base and Kerimasi Base and these were occupied for up to 12 hours.

Primary Base GPS Station at Kerimasi Base camp One of the Primary Base GPS Station was at Kerimasi Base camp with Ol Doinyo Lengai dominating the horizon.

The data from three ISG Earth Stations was downloaded.  These were based at Malindi on the Kenyan coast, Mbarara in Uganda and Mahe in the Seychelles.  The LGO software has a limit of 400 Km for resolving base lines and Malindi was the only Earth Station in this range so all calculations related back to this ISG Earth Station.  The data from Mbarara and Mahe was used to check the results from the Malindi calculations and used as a guide to determine the accuracy of the calculations.

The data for Malindi and Mahe was good data, but Mbarara had been experiencing difficulties so the results were provisional and only accurate to some 0.3 metres, but this still provided a good check on the results from Malindi.  Had the data from Mbarara been good it would have been a useful Earth Station to use, as its altitude was similar to the points surveyed during the expedition, whereas Malindi and Mahe are at sea level which means that the GPS signals had an additional Kilometre of atmosphere to pass through.

The positional information for the ISG Earth is typically displayed as:
33201M001 MALINDI       GPS MALI  4865366.487  4110737.445  -331121.710
39801M001 MAHE ISLAND   GPS SEY1  3602870.766  5238174.329  -516275.514

The ISG Earth Station at Malindi on the Kenyan coast.The GPS receiver at the ISG Earth Station at Malindi on the Kenyan coast is located on the red concrete pillar to the right of the photograph.

Results

Whilst collecting the GPS data a pattern of “check bases” and “closed traverses” we were used so that it was possible to establish that the precision of the results.  These showed that they were better than 0.1m in both plan and height with many of the points showing much better results than this, so we achieved our goal in an efficient manner. 

During the course of the expedition the Leica  SR530 receiver, aluminium tripod and gravity meter were backpack packed to some very difficult locations.  Not only was it taken to the top of Ol Doinyo Lengai and used to traverse down its precipitous flanks but it was also taken up the western wall of the Rift Valley, by way of Laws Ridge, to the complete the gravity traverse across the valley floor.

The Gravity Meter and GPS points followed the route taken by the expedition across the floor of the Rift Valley starting at Longido in the east across to the region around Ol Doinyo Lengai.  The primary Stations to relate back to the ISG Earth Stations were at Longido, the easternmost point, Kitumbeine, the southernmost point, and Kerimasi and Lengai base camps in the west.

Gravity Meter and GPS points
As Ol Doinyo Lengai was the focus of the expedition, a greater density of points were recorded in this location.  Points were measured down the flank of the volcano, and the two westernmost points  at the top of the western wall of the Rift Valley.

Gravity Meter and GPS points

The data collected from the four primary GPS Stations were first computed to the ISG Earth Station at Malindi to establish their absolute location on the earth’s surface.  Check computations were also made to the ISG Earth Stations at Mbarara and Mahe and these agreed to 0.3 metres in plan and height with the values obtained from the calculations to Malindi.

The base lines observed during the expedition were computed using these four primary stations as the reference stations.  Using check measurements and common points the precision of the survey was determined to be better than 0.1 metres in both plan and height.

Bearing in mind that ambiguities were only resolved to the ISG Earth Station at Malindi and not to those at Mahe or Mbarara, the overall precision of the absolute locations of the GPS points measured during the expedition was achieved to 0.1 metres in plan and height, which was the desired result for providing information for the Gravity Meter survey and for related studies of investigating the accuracy of a Navigation GPS determination and calculating the height of Ol Doinyo Lengai.

Position in Degrees, Minutes and Seconds for the points determined by DGPS

Point Type Date Latitude S Longitude E Height
LONGIDOA Reference 07/19/2004 09:00    2° 43′ 46.65576″ S   36° 43′ 27.08642″ E 1342.88
LONGIDOB Averaged 07/19/2004 06:12    2° 43′ 46.62238″ S   36° 43′ 27.05667″ E 1342.98
STN002 Reference 07/21/2004 09:45    2° 43′ 56.24684″ S   36° 41′ 44.03650″ E 1326.90
STN003 Reference 07/21/2004 10:55    2° 43′ 15.75038″ S   36° 39′ 30.96916″ E 1317.45
STN004 Reference 07/21/2004 12:19    2° 42′ 21.62493″ S   36° 32′ 39.33202″ E 1427.09
STN005 Reference 07/21/2004 13:41    2° 45′ 48.85434″ S   36° 25′ 44.33227″ E 1268.15
KITUMBEINE Reference 07/22/2004 07:43    2° 50′ 00.50408″ S   36° 18′ 42.19313″ E 1200.51
STN007 Measured 07/22/2004 08:36    2° 46′ 02.01005″ S   36° 16′ 07.52006″ E 1082.77
STN008 Measured 07/22/2004 10:15    2° 42′ 55.32037″ S   36° 11′ 12.96449″ E 1031.68
LEN005 Measured 07/30/2004 05:24    2° 45′ 53.56600″ S   35° 55′ 01.68902″ E 2804.40
LEN006 Measured 07/30/2004 06:42    2° 45′ 35.30156″ S   35° 54′ 51.34849″ E 2818.59
LEN007 Measured 07/30/2004 07:14    2° 45′ 38.27198″ S   35° 54′ 46.89467″ E 2802.67
LEN008 Measured 07/30/2004 08:29    2° 45′ 30.70083″ S   35° 54′ 39.46206″ E 2709.24
LEN009 Measured 07/30/2004 09:38    2° 45′ 28.03403″ S   35° 54′ 34.88932″ E 2594.46
LEN010 Measured 07/30/2004 10:36    2° 45′ 23.84674″ S   35° 54′ 26.26206″ E 2409.58
LEN011 Measured 07/30/2004 11:15    2° 45′ 22.59607″ S   35° 54′ 19.01902″ E 2261.07
LEN012 Measured 07/30/2004 12:13    2° 45′ 17.40242″ S   35° 54′ 02.42677″ E 1944.18
LENTRAN Measured 07/31/2004 07:13    2° 43′ 48.97792″ S   35° 52′ 54.20871″ E 1107.61
GPSBL001 Measured 07/31/2004 12:26    2° 43′ 49.83288″ S   35° 52′ 57.44591″ E 1107.67
GPSBL002 Measured 07/31/2004 12:31    2° 43′ 49.83287″ S   35° 52′ 57.44599″ E 1107.68
GPSBL003 Measured 07/31/2004 19:22    2° 45′ 14.73348″ S   35° 53′ 55.05897″ E 1841.07
GPSBL004 Measured 07/31/2004 20:07    2° 45′ 08.57603″ S   35° 53′ 40.30068″ E 1649.55
GPSBL005 Measured 07/31/2004 20:59    2° 44′ 37.40312″ S   35° 53′ 22.52740″ E 1348.63
GPSBL006 Measured 07/31/2004 21:41    2° 44′ 06.36231″ S   35° 53′ 07.90498″ E 1174.21
GPSBL007 Measured 08/01/2004 06:03    2° 43′ 48.20376″ S   35° 52′ 57.35539″ E 1103.56
LENGAI_RO Measured 08/01/2004 06:24    2° 43′ 49.23083″ S   35° 52′ 57.52020″ E 1106.74
GPSBL009 Measured 08/02/2004 06:44    2° 43′ 31.73961″ S   35° 53′ 36.96484″ E 1063.00
GPSBL010 Measured 08/02/2004 07:45    2° 43′ 27.27959″ S   35° 53′ 55.97718″ E 1058.87
GPSBL011 Measured 08/02/2004 08:28    2° 43′ 24.92299″ S   35° 54′ 05.99968″ E 1053.18
GPSBL013 Measured 08/03/2004 05:47    2° 44′ 24.69355″ S   35° 51′ 39.53965″ E 1708.63
GPSBL014 Measured 08/03/2004 07:00    2° 44′ 37.11617″ S   35° 51′ 39.72796″ E 1731.14
KERIMASI Reference 08/05/2004 07:27    2° 48′ 04.96857″ S   35° 58′ 25.64845″ E 1029.98
GPSBL016 Measured 08/05/2004 09:55    2° 45′ 20.32756″ S   36° 05′ 18.27383″ E 1018.57
GPSBL017 Measured 08/05/2004 11:16    2° 47′ 05.26248″ S   36° 03′ 14.04968″ E 915.09
GPSBL018 Measured 08/05/2004 12:45    2° 47′ 13.52203″ S   36° 00′ 29.33248″ E 946.62
LENGAI Reference 09/06/2004 11:12    2° 43′ 49.24388″ S   35° 52′ 57.48465″ E 1106.74
STN009 Measured 09/06/2004 11:27    2° 44′ 22.39308″ S   36° 06′ 45.41933″ E 1103.43
GM1 Measured 07/31/2004    2° 43′ 45.78896″ S   35° 52′ 53.26600″ E 1106.58
GM2 Measured 07/31/2004    2° 43′ 40.58746″ S   35° 53′ 10.03887″ E 1084.45
GM3 Measured 07/31/2004    2° 43′ 37.77218″ S   35° 53′ 22.07219″ E 1084.59
GM4 Measured 07/31/2004    2° 43′ 48.47234″ S   35° 52′ 41.30072″ E 1090.99

 

Ellipsoid and Orthometric Heights

The surface of the earth is not a truly mathematical shape so is approximated to an ellipsoid which is a mathematical figure formed by rotating an ellipse about one of its axis.  In the case of the earth, the shortest (north south) axis is used.

Ellipsoidal heights are the heights above this theoretical figure and are the values determined by GPS observations as these are related to an ellipsoid.  To complicate the matter further, different ellipsoids are used for different parts of the earths surface to get the “best fit”. 

The difference between Ellipsoidal and Orthometric heights.Orthometric heights are heights above the Geoid.  The Geoid is the true shape of the earth and is described as the equipotential surface best fitted to mean sea level.  This is a surface where the pull of gravity is the same for all points on the earths surface.  As there are an infinite number of such surfaces, the one closest to mean sea level is chosen.

The difference between Ellipsoidal and Orthometric heights.
h  =  the height of P above the Ellipsoid, shown in green.
H  =  the height of P above the Geoid, shown in orange.

As can be seen from the data the difference between the Geoid and the Ellipsoid in the region surveyed is only around 10 metres.  This varies around the globe in the United Kingdom the Geoidal Separation, as it is called, is some 50 metrers.

DGPS

Differential GPS has a high relative accuracy.

Differential GPS.GPS data is collected simultaneously at both ends of a line providing the difference in co-ordinates between the two points rather than the actual location of a single point.

Ambiguity Resolution

Survey accuracy GPS uses the phase difference in the carrier wave from the GPS satellites.  The phase difference at each receiver can be measured, but the “unknown” is the number of whole cycles of the carrier waves.  As the satellites move across the receivers of field of view the change in number of whole cycles as well as the phase differences can be recorded.  The “ambiguity” that needs to be resolved is the number of whole cycles.  

In theory, position is three “unknowns”, x, y and z, so in theory the data from three satellites is required to resolve these ambiguities.  However, despite the use of high precision time keeping devices, actual time is also an unknown so a minimum of four satellites is required.  The data observed simultaneously by each receiver from a minimum of four satellites is used to resolve the ambiguities by the method of least squares.

Seen from the data the difference between the Geoid and the Ellipsoid in the region surveyed is only around 10 metres.  This varies around the globe in the United Kingdom the Geoidal Separation, as it is called, is some 50 metrers.

Back to the Top of the Page