Fixed point theorems of GPS carrier phase ambiguity resolution and their application to massive network processing: Ambizap

[1] Precise point positioning (PPP) has become popular for Global Positioning System (GPS) geodetic network analysis because for n stations, PPP has O(n) processing time, yet solutions closely approximate those of O(n3) full network analysis. Subsequent carrier phase ambiguity resolution (AR) further improves PPP precision and accuracy; however, full-network bootstrapping AR algorithms are O(n4), limiting single network solutions to n < 100. In this contribution, fixed point theorems of AR are derived and then used to develop “Ambizap,” an O(n) algorithm designed to give results that closely approximate full network AR. Ambizap has been tested to n ≈ 2800 and proves to be O(n) in this range, adding only ∼50% to PPP processing time. Tests show that a 98-station network is resolved on a 3-GHz CPU in 7 min, versus 22 h using O(n4) AR methods. Ambizap features a novel network adjustment filter, producing solutions that precisely match O(n4) full network analysis. The resulting coordinates agree to ≪1 mm with current AR methods, much smaller than the ∼3-mm RMS precision of PPP alone. A 2000-station global network can be ambiguity resolved in ∼2.5 h. Together with PPP, Ambizap enables rapid, multiple reanalysis of large networks (e.g., ∼1000-station EarthScope Plate Boundary Observatory) and facilitates the addition of extra stations to an existing network solution without need to reprocess all data. To meet future needs, PPP plus Ambizap is designed to handle ∼10,000 stations per day on a 3-GHz dual-CPU desktop PC.

[1]  James L. Davis,et al.  Assessment of GPS velocity accuracy for the Basin and Range Geodetic Network (BARGEN) , 2003 .

[2]  Michael R. Craymer,et al.  Observation of glacial isostatic adjustment in “stable” North America with GPS , 2007 .

[3]  Robert J. Renka,et al.  Algorithm 772: STRIPACK: Delaunay triangulation and Voronoi diagram on the surface of a sphere , 1997, TOMS.

[4]  Franz Aurenhammer,et al.  Voronoi diagrams—a survey of a fundamental geometric data structure , 1991, CSUR.

[5]  Geoffrey Blewitt,et al.  GPS and Space-Based Geodetic Methods , 2007 .

[6]  P. Teunissen.  The least-squares ambiguity decorrelation adjustment: a method for fast GPS integer ambiguity estimation , 1995 .

[7]  Gerd Gendt,et al.  The International GPS Service: Celebrating the 10th anniversary and looking to the next decade , 2005 .

[8]  G. Blewitt.  Carrier Phase Ambiguity Resolution for the Global Positioning System Applied to Geodetic Baselines up to 2000 km , 1989 .

[9]  Yehuda Bock,et al.  Southern California permanent GPS geodetic array: Spatial filtering of daily positions for estimating coseismic and postseismic displacements induced by the 1992 Landers earthquake , 1997 .

[10]  J. Kruskal.  On the shortest spanning subtree of a graph and the traveling salesman problem , 1956 .

[11]  J. Zumberge,et al.  Precise point positioning for the efficient and robust analysis of GPS data from large networks , 1997 .

[12]  J. Nocquet,et al.  Deformation of the North American plate interior from a decade of continuous GPS measurements , 2006 .

[13]  Zvi Galil,et al.  Data structures and algorithms for disjoint set union problems , 1991, CSUR.

[14]  Y. Bock,et al.  Global Positioning System Network analysis with phase ambiguity resolution applied to crustal deformation studies in California , 1989 .