Precise point positioning for the efficient and robust analysis of GPS data from large networks
TL;DR: This work determines precise GPS satellite positions and clock corrections from a globally distributed network of GPS receivers, and analysis of data from hundreds to thousands of sites every day with 40-Mflop computers yields results comparable in quality to the simultaneous analysis of all data.
Abstract: Networks of dozens to hundreds of permanently operating precision Global Positioning System (GPS) receivers are emerging at spatial scales that range from 10(exp 0) to 10(exp 3) km. To keep the computational burden associated with the analysis of such data economically feasible, one approach is to first determine precise GPS satellite positions and clock corrections from a globally distributed network of GPS receivers. Their, data from the local network are analyzed by estimating receiver- specific parameters with receiver-specific data satellite parameters are held fixed at their values determined in the global solution. This "precise point positioning" allows analysis of data from hundreds to thousands of sites every (lay with 40-Mflop computers, with results comparable in quality to the simultaneous analysis of all data. The reference frames for the global and network solutions can be free of distortion imposed by erroneous fiducial constraints on any sites.
Citations
More filters
TL;DR: MORVEL as discussed by the authors is a new closure-enforced set of angular velocities for the geologically current motions of 25 tectonic plates that collectively occupy 97 per cent of Earth's surface.
Abstract: SUMMARY
We describe best-fitting angular velocities and MORVEL, a new closure-enforced set of angular velocities for the geologically current motions of 25 tectonic plates that collectively occupy 97 per cent of Earth's surface. Seafloor spreading rates and fault azimuths are used to determine the motions of 19 plates bordered by mid-ocean ridges, including all the major plates. Six smaller plates with little or no connection to the mid-ocean ridges are linked to MORVEL with GPS station velocities and azimuthal data. By design, almost no kinematic information is exchanged between the geologically determined and geodetically constrained subsets of the global circuit—MORVEL thus averages motion over geological intervals for all the major plates. Plate geometry changes relative to NUVEL-1A include the incorporation of Nubia, Lwandle and Somalia plates for the former Africa plate, Capricorn, Australia and Macquarie plates for the former Australia plate, and Sur and South America plates for the former South America plate. MORVEL also includes Amur, Philippine Sea, Sundaland and Yangtze plates, making it more useful than NUVEL-1A for studies of deformation in Asia and the western Pacific. Seafloor spreading rates are estimated over the past 0.78 Myr for intermediate and fast spreading centres and since 3.16 Ma for slow and ultraslow spreading centres. Rates are adjusted downward by 0.6–2.6 mm yr−1 to compensate for the several kilometre width of magnetic reversal zones. Nearly all the NUVEL-1A angular velocities differ significantly from the MORVEL angular velocities. The many new data, revised plate geometries, and correction for outward displacement thus significantly modify our knowledge of geologically current plate motions. MORVEL indicates significantly slower 0.78-Myr-average motion across the Nazca–Antarctic and Nazca–Pacific boundaries than does NUVEL-1A, consistent with a progressive slowdown in the eastward component of Nazca plate motion since 3.16 Ma. It also indicates that motions across the Caribbean–North America and Caribbean–South America plate boundaries are twice as fast as given by NUVEL-1A. Summed, least-squares differences between angular velocities estimated from GPS and those for MORVEL, NUVEL-1 and NUVEL-1A are, respectively, 260 per cent larger for NUVEL-1 and 50 per cent larger for NUVEL-1A than for MORVEL, suggesting that MORVEL more accurately describes historically current plate motions. Significant differences between geological and GPS estimates of Nazca plate motion and Arabia–Eurasia and India–Eurasia motion are reduced but not eliminated when using MORVEL instead of NUVEL-1A, possibly indicating that changes have occurred in those plate motions since 3.16 Ma. The MORVEL and GPS estimates of Pacific–North America plate motion in western North America differ by only 2.6 ± 1.7 mm yr−1, ≈25 per cent smaller than for NUVEL-1A. The remaining difference for this plate pair, assuming there are no unrecognized systematic errors and no measurable change in Pacific–North America motion over the past 1–3 Myr, indicates deformation of one or more plates in the global circuit. Tests for closure of six three-plate circuits indicate that two, Pacific–Cocos–Nazca and Sur–Nubia–Antarctic, fail closure, with respective linear velocities of non-closure of 14 ± 5 and 3 ± 1 mm yr−1 (95 per cent confidence limits) at their triple junctions. We conclude that the rigid plate approximation continues to be tremendously useful, but—absent any unrecognized systematic errors—the plates deform measurably, possibly by thermal contraction and wide plate boundaries with deformation rates near or beneath the level of noise in plate kinematic data.
2,089 citations
TL;DR: This paper will describe the approach, summarize the adjustment procedure, and specify the earth- and space-based models that must be implemented to achieve cm-level positioning in static mode and station tropospheric zenth path delays with cm precision.
Abstract: The contribution details a post-processing approach that used undifferentiated dual-frequency pseudorange and carrier phase observations along with IGS procise orbit products, for stand-alone precise geodetic point positioning (static or kinematic) with cm precision. This is possible if one takes advantage of the satellite clock estimates available with the satellite coordinates in the IGS precise orbit products and models systematic effects that cause cm variations in the satelite to user range. This paper will describe the approach, summarize the adjustment procedure, and specify the earth- and space-based models that must be implementetd to achieve cm-level positioning in static mode. Furthermore, station tropospheric zenth path delays with cm precision and GPS receiver clock estimates procise to 0.1 ns are also obtained. © 2001 John Wiley & Sons, Inc.
1,200 citations
TL;DR: Global Positioning System (GPS) measurements in China indicate that crustal shortening accommodates most of India's penetration into Eurasia, but the Tibetan plateau south of the Kunlun and Ganzi-Mani faults is moving eastward relative to both India and Eurasia.
Abstract: Global Positioning System (GPS) measurements in China indicate that crustal shortening accommodates most of India's penetration into Eurasia. Deformation within the Tibetan Plateau and its margins, the Himalaya, the Altyn Tagh, and the Qilian Shan, absorbs more than 90% of the relative motion between the Indian and Eurasian plates. Internal shortening of the Tibetan plateau itself accounts for more than one-third of the total convergence. However, the Tibetan plateau south of the Kunlun and Ganzi-Mani faults is moving eastward relative to both India and Eurasia. This movement is accommodated through rotation of material around the eastern Syntaxis. The North China and South China blocks, east of the Tibetan Plateau, move coherently east-southeastward at rates of 2 to 8 millimeters per year and 6 to 11 millimeters per year, respectively, with respect to the stable Eurasia.
1,019 citations
TL;DR: In this article, the relative velocities of 19 plates and continental blocks were derived from publicly available space geodetic (primarily GPS) data for the period 1993-2000, including an independent and rigorous estimate for GPS velocity uncertainties to assess plate rigidity and propagate these uncertainties to the velocity estimates.
Abstract: [1] We present a new global model for Recent plate velocities, REVEL, describing the relative velocities of 19 plates and continental blocks. The model is derived from publicly available space geodetic (primarily GPS) data for the period 1993–2000. We include an independent and rigorous estimate for GPS velocity uncertainties to assess plate rigidity and propagate these uncertainties to the velocity estimates. The velocity fields for North America, Eurasia, and Antarctica clearly show the effects of glacial isostatic adjustment, and Australia appears to depart from rigid plate behavior in a manner consistent with the mapped intraplate stress field. Two thirds of tested plate pairs agree with the NUVEL-1A geologic (3 Myr average) velocities within uncertainties. Three plate pairs (Caribbean–North America, Caribbean–South America, and North America–Pacific) exhibit significant differences between the geodetic and geologic model that may reflect systematic errors in NUVEL-1A due to the use of seafloor magnetic rate data that do not reflect the full plate rate because of tectonic complexities. Most other differences probably reflect real velocity changes over the last few million years. Several plate pairs (Arabia–Eurasia, Arabia–Nubia, Eurasia–India) move more slowly than the 3 Myr NUVEL-1A average, perhaps reflecting long-term deceleration associated with continental collision. Several other plate pairs, including Nazca–Pacific, Nazca–South America and Nubia–South America, are experiencing slowing that began ∼25 Ma, the beginning of the current phase of Andean crustal shortening.
954 citations
Cites methods from "Precise point positioning for the e..."
...We use GIPSY/OASIS II, Release 5.0 software developed at the Jet Propulsion Laboratory (JPL) and nonfiducial satellite orbit and clock files provided by JPL [Zumberge et al., 1997]....
[...]
...0 software developed at the Jet Propulsion Laboratory (JPL) and nonfiducial satellite orbit and clock files provided by JPL [Zumberge et al., 1997]....
[...]
...We thank R. Stern for suggesting the Philippine plate study; D. Jeffferson, K. Hurst, and M. Heflin at JPL, M. Mader at NGS; and Z. Altamimi, at IERS for technical assistance and advice....
[...]
TL;DR: In this paper, a new model of the last deglaciation event of the Late Quaternary ice age is described and denoted as ICE-6G_C (VM5a), which has been explicitly refined by applying all of the available Global Positioning System (GPS) measurements of vertical motion of the crust that may be brought to bear to constrain the thickness of local ice cover as well as the timing of its removal.
Abstract: A new model of the last deglaciation event of the Late Quaternary ice age is here described and denoted as ICE-6G_C (VM5a). It differs from previously published models in this sequence in that it has been explicitly refined by applying all of the available Global Positioning System (GPS) measurements of vertical motion of the crust that may be brought to bear to constrain the thickness of local ice cover as well as the timing of its removal. Additional space geodetic constraints have also been applied to specify the reference frame within which the GPS data are described. The focus of the paper is upon the three main regions of Last Glacial Maximum ice cover, namely, North America, Northwestern Europe/Eurasia, and Antarctica, although Greenland and the British Isles will also be included, if peripherally, in the discussion. In each of the three major regions, the model predictions of the time rate of change of the gravitational field are also compared to that being measured by the Gravity Recovery and Climate Experiment satellites as an independent means of verifying the improvement of the model achieved by applying the GPS constraints. Several aspects of the global characteristics of this new model are also discussed, including the nature of relative sea level history predictions at far-field locations, in particular the Caribbean island of Barbados, from which especially high-quality records of postglacial sea level change are available but which records were not employed in the development of the model. Although ICE-6G_C (VM5a) is a significant improvement insofar as the most recently available GPS observations are concerned, comparison of model predictions with such far-field relative sea level histories enables us to identify a series of additional improvements that should follow from a further stage of model iteration.
902 citations
References
More filters
TL;DR: In this article, a technique for resolving the ambiguities in the GPS carrier phase data (which are biased by an integer number of cycles) is described which can be applied to geodetic baselines up to 2000 km in length and can be used with dual-frequency P code receivers.
Abstract: A technique for resolving the ambiguities in the GPS carrier phase data (which are biased by an integer number of cycles) is described which can be applied to geodetic baselines up to 2000 km in length and can be used with dual-frequency P code receivers. The results of such application demonstrated that a factor of 3 improvement in baseline accuracy could be obtained, giving centimeter-level agreement with coordinates inferred by very-long-baseline interferometry in the western United States. It was found that a method using pseudorange data is more reliable than one using ionospheric constraints for baselines longer than 200 km. It is recommended that future GPS networks have a wide spectrum of baseline lengths (ranging from baselines shorter than 100 km to those longer than 1000 km) and that GPS receivers be used which can acquire dual-frequency P code data.
578 citations
TL;DR: In this paper, the velocity field in central and southern California using Global Positioning System (GPS) observations from 1986 to 1992 and very long baseline interferometry (VLBI) observations was estimated.
Abstract: We estimate the velocity field in central and southern California using Global Positioning System (GPS) observations from 1986 to 1992 and very long baseline interferometry (VLBI) observations from 1984 to 1991. Our core network includes 12 GPS sites spaced approximately 50 km apart, mostly in the western Transverse Ranges and the coastal Borderlands. The precision and accuracy of the relative horizontal velocities estimated for these core stations are adequately described by a 95% confidence ellipse with a semiminor axis of approximately 2 mm/yr oriented roughly north-south, and a semimajor axis of approximately 3 mm/yr oriented east-west.
252 citations
TL;DR: An overview of applications of the Global Positioning System (GPS) for active measurement of the Earth's atmosphere can be found in this paper, where GPS signals are delayed due to the presence of water vapor.
Abstract: This paper provides an overview of applications of the Global Positioning System (GPS) for active measurement of the Earth's atmosphere. Microwave radio signals transmitted by GPS satellites are delayed (refracted) by the atmosphere as they propagate to Earth-based GPS receivers or GPS receivers carried on low Earth orbit satellites. The delay in GPS signals reaching Earth-based receivers due to the presence of water vapor is nearly proportional to the quantity of water vapor integrated along the signal path. Measurement of atmospheric water vapor by Earth-based GPS receivers was demonstrated during the GPS/STORM field project to be comparable and in some respects superior to measurements by ground-based water vapor radiometers. Increased spatial and temporal resolution of the water vapor distribution provided by the GPS/STORM network proved useful in monitoring the moisture-flux convergence along a dryline and the decrease in integrated water vapor associated with the passage of a midtropospheri...
246 citations
TL;DR: In this article, baseline lengths and geocentric radii have been determined from GPS data without the use of fiducial sites and the geodetic information can be obtained by examining the structure of the polyhedron and its change with time.
Abstract: Baseline lengths and geocentric radii have been determined from GPS data without the use of fiducial sites. Data from the first GPS experiment for the IERS and Geodynamics (GIG '91) have been analyzed with a no-fiducial strategy. A baseline length daily repeatability of 2 mm + 4 parts per billion was obtained for baselines in the Northern Hemisphere. Comparison of baseline lengths from GPS and the global VLBI solution GLB659 (Caprette et al. 1990) show rms agreement of 2.1 parts per billion. The geocentric radius mean daily repeatability for all sites was 15 cm. Comparison of geocentric radii from GPS and SV5 (Murray et al. 1990) show rms agreement of 3.8 cm. Given n globally distributed stations, the n(n - 1)/2 baseline lengths and n geocentric radii uniquely define a rigid closed polyhedron with a well-defined center of mass. Geodetic information can be obtained by examining the structure of the polyhedron and its change with time.
150 citations