Tanzania 2004
 
 

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.03m in height, although accurately reading this is extremely tricky, so we needed to achieve relative heights between successive gravity reading of better than 0.1m with 0.03m 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.

 

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 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 is located on the red concrete pillar to the left 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 western wall of the Gregory Rift Valley from Lengai Base camp, with Laws Ridge to the left.

 

The Gravity Meter Traverse extended east to the first campsite at Longido.  To complete the traverse points were required on the top of the western wall and access was gained via the ridge, which was previously unclimbed.

 

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.

 

 

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.

 

 

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”. 

 

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 o 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 50metrers.

 

 

DGPS

Differential GPS.

 

High relative accuracy.

 

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 50metrers.