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.

https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/96JB03860

Precise point positioning for the efficient and robust analysis of GPS data from large networks

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.

Geologically current plate motions

GEODETIC data, obtained by ground- or space-based techniques, can be used to infer the distribution of slip on a fault that has ruptured in an earthquake. Although most geodetic techniques require a surveyed network to be in place before the earthquake1–3, satellite images, when collected at regular intervals, can capture co-seismic displacements without advance knowledge of the earthquake's location. Synthetic aperture radar (SAR) interferometry, first introduced4 in 1974 for topographic mapping5–8 can also be used to detect changes in the ground surface, by removing the signal from the topography9,10. Here we use SAR interferometry to capture the movements produced by the 1992 earthquake in Landers, California11. We construct an interferogram by combining topographic information with SAR images obtained by the ERS-1 satellite before and after the earthquake. The observed changes in range from the ground surface to the satellite agree well with the slip measured in the field, with the displacements measured by surveying, and with the results of an elastic dislocation model. As a geodetic tool, the SAR interferogram provides a denser spatial sampling (100 m per pixel) than surveying methods1–3 and a better precision (∼3 cm) than previous space imaging techniques12,13.

The displacement field of the Landers earthquake mapped by radar interferometry

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.

Precise Point Positioning Using IGS Orbit and Clock Products

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.

https://science.sciencemag.org/content/sci/294/5542/574.full.pdf

Present-Day Crustal Deformation in China Constrained by Global Positioning System Measurements

The M = 7.1 Hector Mine earthquake ruptured the Lavic Lake fault near Twentynine Palms, CA at 09:46 UTC October 16, 1999. Because it occurred near the eastern edge of the Southern California Integrated GPS Network (SCIGN), a network of permanent, continuously recording GPS receivers for measuring the crustal deformation field around Los Angeles, CA, it was possible to determine the deformation associated with the earthquake with unprecedented speed and reliability. Thirty-four stations recorded displacements over the 3-sigma level. The displacements measured with GPS can be modeled by a fault 46.2±2.6 km long, 8.2±1.0 km wide, striking 330°, dipping 84° east, with 301±36 cm right lateral strike-slip, and 145±36 cm of east-up dip-slip, yielding a potency of 1.3 km³ and geodetic moment of 3.8 × 1026 dyne-cm. The trace and dip of the model fault is consistent with the observed ground rupture and seismic focal mechanisms.

https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2000GL011841

The coseismic geodetic signature of the 1999 Hector Mine earthquake

Use of GPS tracking data from different dual-frequency receiver types (cross-correlating vs. codeless) has revealed satellite-dependent biases in pseudorange observables P1 (Y-code) and C1 (C/A, Clear Acquisition code). These biases can have a direct effect on clock estimates, carrier phase bias fixing, and other parameters estimated in GPS data processing. A set of satellite-specific compensatory pseudorange offsets is calculated, and each is applied to a wee of daily global network analyses in which satlellite, receiver, atmospheric, and Earth rotation parameters are estimated. Results from these analyses are then compared to those from corresponding baseline cases in which no biases were applied. There is also some evidence that suggests that the pseudorange biases differ even among codeless receiver models. Hence, a second set of offsets is computed on a different basis, and compared with the baseline model in a similar manner. A preliminary examination of C1-P1 variations over time is presented. Finally, recommendations are made for the use of the calculated offsets, and consideration is given to a future dissemination of updates to these values as necessary. © 2001 John Wiley & Sons, Inc.

Examining the C1-P1 pseudorange bias

Geodetic time series determined with the Global Positioning System indicate that the geodetic rate of a permanent site in Pasadena, California (JPLM) changed significantly after the 17 January 1994 Northridge California earthquake. Subtracting the pre-quake rate and co-seismic offset leaves 30±4 mm of integrated eastward excess motion observed in the three years following the earthquake. North and vertical components show excess motion of −11 plusmn;3 mm and 25±11 mm respectively. Local surveys to three additional points near JPLM changed by no more than 6 mm E, 3 mm N, and 15 mm V during the two years after the earthquake, ruling out the possibility of a local effect at the JPLM monument. The direction and size of the post-seismic displacements at JPLM are not consistent with additional slip on the fault which ruptured. The most rapid accumulation of excess motion occurs immediately after the earthquake, suggesting a relationship between the two events.

https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/97GL03397

Rate change observed at JPLM after the Northridge Earthquake

Jet Propulsion Laboratory IGS Analysis Center Report

Mars Science Laboratory (MSL) Orbit Determination (OD) met all requirements with considerable margin, MSL OD team developed spin signature removal tool and successfully used the tool during cruise, A novel approach was used for the MSL solar radiation pressure model and resulted in a very accurate model during the approach phase, The change in velocity for Attitude Control System (ACS) turns was successfully calibrated and with appropriate scale factor resulted in improved change in velocity prediction for future turns, All Trajectory Correction Maneuvers were successfully reconstructed and execution errors were well below the assumed pre-fight execution errors, The official OD solutions were statistically consistent throughout cruise and for OD solutions with different arc lengths as well, Only EPU-1 was sent to MSL. All other Entry Parameter Updates were waived, EPU-1 solution was only 200 m separated from final trajectory reconstruction in the B-plane

David C. Jefferson

Papers

The coseismic geodetic signature of the 1999 Hector Mine earthquake

Examining the C1-P1 pseudorange bias

Rate change observed at JPLM after the Northridge Earthquake

Jet Propulsion Laboratory IGS Analysis Center Report

Mars Science Laboratory Orbit Determination