基于全球板块运动模型分析大西洋扩张变化

基于现代空间测量技术SLR、VLBI和GPS实测资料，解算出大西洋中脊海底扩张速率，其中北大西洋的东西向扩张速率平均为35mm/a，赤道大西洋东西向扩张速度分别为20－25mm/a，南大西洋东西向扩张速率为22－28mm/a，证实全球板块运动的存在及大西洋扩张学说，并基于全球几百万年地质模型NNR－NUVEL1A，北大西洋的东西向扩张速率平均为24．3mm/a，基于最新全球板块运动模型ITRF2000VEL，北大西洋的东西向扩张速率平均为20．8mm/a，总体上大西洋实测东西扩张速度与根据地学资料推出的地球板块运动模型和最新ITRF2000VEL模型的结果基本一致。     

ITRF2008: an improved solution of the international terrestrial reference frame

ITRF2008 is a refined version of the International Terrestrial Reference Frame based on reprocessed solutions of the four space geodetic techniques: VLBI, SLR, GPS and DORIS, spanning 29, 26, 12.5 and 16?years of observations, respectively. The input data used in its elaboration are time series (weekly from satellite techniques and 24-h session-wise from VLBI) of station positions and daily Earth Orientation Parameters (EOPs). The ITRF2008 origin is defined in such a way that it has zero translations and translation rates with respect to the mean Earth center of mass, averaged by the SLR time series. Its scale is defined by nullifying the scale factor and its rate with respect to the mean of VLBI and SLR long-term solutions as obtained by stacking their respective time series. The scale agreement between these two technique solutions is estimated to be 1.05 ± 0.13 ppb at epoch 2005.0 and 0.049 ± 0.010?ppb/yr. The ITRF2008 orientation (at epoch 2005.0) and its rate are aligned to the ITRF2005 using 179 stations of high geodetic quality. An estimate of the origin components from ITRF2008 to ITRF2005 (both origins are defined by SLR) indicates differences at epoch 2005.0, namely: ?0.5, ?0.9 and ?4.7?mm along X, Y and Z-axis, respectively. The translation rate differences between the two frames are zero for Y and Z, while we observe an X-translation rate of 0.3?mm/yr. The estimated formal errors of these parameters are 0.2?mm and 0.2?mm/yr, respectively. The high level of origin agreement between ITRF2008 and ITRF2005 is an indication of an imprecise ITRF2000 origin that exhibits a Z-translation drift of 1.8?mm/yr with respect to ITRF2005. An evaluation of the ITRF2008 origin accuracy based on the level of its agreement with ITRF2005 is believed to be at the level of 1?cm over the time-span of the SLR observations. Considering the level of scale consistency between VLBI and SLR, the ITRF2008 scale accuracy is evaluated to be at the level of 1.2?ppb (8?mm at the equator) over the common time-span of the observations of both techniques. Although the performance of the ITRF2008 is demonstrated to be higher than ITRF2005, future ITRF improvement resides in improving the consistency between local ties in co-location sites and space geodesy estimates.     

DORIS Contribution to ITRF2005

We examine the contribution of the Doppler Orbit determination and Radiopositioning Integrated by Satellite (DORIS) technique to the International Terrestrial Reference Frame (ITRF2005) by evaluating the quality of the submitted solutions as well as that of the frame parameters, especially the origin and the scale. Unlike the previous versions of the ITRF, ITRF2005 is constructed with input data in the form of time-series of station positions (weekly for satellite techniques and daily for VLBI) and daily Earth orientation parameters (EOPs), including full variance–covariance information. Analysis of the DORIS station positions’ time-series indicates an internal precision reaching 15 mm or better, at a weekly sampling. A cumulative solution using 12 years of weekly time-series was obtained and compared to a similar International GNSS Service (IGS) GPS solution (at 37 co-located sites) yielding a weighted root mean scatter (WRMS) of the order of 8 mm in position (at the epoch of minimum variance) and about 2.5 mm/year in velocity. The quality of this cumulative solution resulting from the combination of two individual DORIS solutions is better than any individual solution. A quality assessment of polar motion embedded in the contributed DORIS solutions is performed by comparison with the results of other space-geodetic techniques and in particular GPS. The inferred WRMS of polar motion varies significantly from one DORIS solution to another and is between 0.5 and 2 mas, depending on the strategy used and in particular estimating or not polar motion rate by the analysis centers. This particular aspect certainly needs more investigation by the DORIS Analysis Centers.     

Defining a DORIS core network for Jason-1 precise orbit determination based on ITRF2000: methods and realization

In view of the future adoption of the new precise orbit determination (POD) standards for the TOPEX/Poseidon and Jason-1 satellites, we propose a method to evaluate terrestrial reference frames for POD. We applied this method to the ITRF2000 realization of the DORIS network using local geodetic ties, plate motion models, the recent DORIS IGN04D02 cumulative solution and DORIS weekly time-series of coordinates. We propose to adopt a selection of the ITRF2000 realization based on specific criteria that we define here, and to extend it with ground stations for which we propose new coordinates and velocities. Only 13 out of 131 stations were considered to be inappropriate for POD activities. The result is a robust and well-distributed DORIS core network of 118 stations (DPOD2000) suitable for POD during the 1993–2008 period considered here.     

Geocentre motion measured with DORIS and SLR, and predicted by geophysical models

Geocentre motion signals measured by satellite geodesy and those predicted from the observed mass redistribution in the ocean, atmosphere and terrestrial waters over 1993.1–2003.0 are analysed and compared under two viewpoints: the amplitudes and phases of the seasonal components, and the spectral signature of the non-seasonal components. The geodetic signals partly match the geophysical variations in the seasonal band, with possible remaining annual and semi-annual errors in both techniques, at the millimetre level in the equatorial plane for Satellite laser ranging (SLR) and Doppler Orbitography and radiopositioning integrated on Satellite (DORIS), and at the centimetre level in T _z (Z-axis translation) for DORIS. Unlike SLR, the DORIS annual signatures in all three geocentre components have strongly varying amplitudes after 1996. The amplitude of the annual geophysical signal in T _y is slowly growing with time. All three geophysical fluids contribute to this effect. The magnitude of the geophysically derived long-term geocentre motion is of the same magnitude in the T _x , T _y and T _z directions, with a 0.5–1.0 mm Allan standard deviation for the 1-year sampling time, while the geodetic values are 2 mm in the equatorial plane for both SLR and DORIS, 4 mm for SLR and 9 mm for DORIS in the T _z direction. The mismatch of the geodetic signal with the geophysical one in the inter-annual band is suggested to be due partly to excessive geodetic noise and partly to underestimated geophysical signal.     

A corrective model for Jason-1 DORIS Doppler data in relation to the South Atlantic Anomaly

The DORIS Doppler measurements collected by Jason-1 are abnormally perturbed by the influence of the South Atlantic Anomaly (SAA). The DORIS ultra-stable oscillators on-board Jason-1 are not as stable as they should be; their frequency is sensitive both to the irradiation rate and to the total irradiation encountered in orbit. The consequence is that not only are the DORIS measurement residuals higher than they ought to be, but also large systematic positioning errors are introduced for stations located in the vicinity of the SAA. In this paper, we present a method that has been devised to obtain a continuous observation of Jason-1 frequency offsets. This method relies on the precise determination of the station frequency and troposphere parameters via the use of other DORIS satellites. More than 3 years of these observations have then been used to construct a model of response of the oscillators of Jason-1 to the SAA. The sensitivity of the Jason-1 oscillators to the SAA perturbations has evolved over time, multiplied by a factor of four between launch and mid-2004. The corrective performances of the model are discussed in terms of DORIS measurement residuals, precise orbit determination and station positioning. The average DORIS measurement residuals are decreased by more than 7 % using this model. In terms of precise orbit determination, the 3D DORIS-only orbit error decreases from 5 to 4.2 cm, but the DORIS+SLR orbit error is almost unaffected, due to the already good quality of this type of orbit. In terms of station positioning, the model brings down the average 3D mono-satellite monthly network solution discrepancy with the International Terrestrial Reference Frame ITRF2000 from 11.3 to 6.1 cm, and also decreases the scatter about that average from 11.3 to 3.7 cm. The conclusion is that, with this model, it is possible to re-incorporate Jason-1 in the multi-satellite geodetic solutions for the DORIS station network.     

Multi-technique comparison of tropospheric zenith delays derived during the CONT02 campaign

In October 2002, 15 continuous days of Very Long Baseline Interferometry (VLBI) data were observed in the Continuous VLBI 2002 (CONT02) campaign. All eight radio telescopes involved in CONT02 were co-located with at least one other space-geodetic technique, and three of them also with a Water Vapor Radiometer (WVR). The goal of this paper is to compare the tropospheric zenith delays observed during CONT02 by VLBI, Global Positioning System (GPS), Doppler Orbitography Radiopositioning Integrated by Satellite (DORIS) and WVR and to compare them also with operational pressure level data from the European Centre for Medium-Range Weather Forecasts (ECMWF). We show that the tropospheric zenith delays from VLBI and GPS are in good agreement at the 3–7 mm level. However, while only small biases can be found for most of the stations, at Kokee Park (Hawaii, USA) and Westford (Massachusetts, USA) the zenith delays derived by GPS are larger by more than 5 mm than those from VLBI. At three of the four DORIS stations, there is also a fairly good agreement with GPS and VLBI (about 10 mm), but at Kokee Park the agreement is only at about 30 mm standard deviation, probably due to the much older installation and type of DORIS equipment. This comparison also allows testing of different DORIS analysis strategies with respect to their real impact on the precision of the derived tropospheric parameters. Ground truth information about the zenith delays can also be obtained from the ECMWF numerical weather model and at three sites using WVR measurements, allowing for comparisons with results from the space-geodetic techniques. While there is a good agreement (with some problems mentioned above about DORIS) among the space-geodetic techniques, the comparison with WVR and ECMWF is at a lower accuracy level. The complete CONT02 data set is sufficient to derive a good estimate of the actual precision and accuracy of each geodetic technique for applications in meteorology.     

相似变换模型在参考框架基准实现中的应用

详细探讨了相似变换在地球参考框架的基准定义实施中的相关模型及方法。利用四种空间大地测量技术（GPS、SLR、DORIS、VLBI）在2000．0历元的站坐标组文件,对ITRF2005的内符精度进行了实验评价。     

Systematic errors in the Z-geocenter derived using satellite tracking data: a case study from SPOT-4 DORIS data in 1998

Within the scope of the Global Geodetic Observing System, Doppler Orbit Determination and Radiopositioning Integrated by Satellite – as a geodetic technique – can provide precise and continuous monitoring of the geocenter motion related to mass redistribution in the Earth, ocean and atmosphere system. We have reanalyzed 1998 DORIS/SPOT-4 (Satellite pour l’ Observation de la Terre) data that were previously generating inconsistent geocenter positions (?65 cm offset). We show here that this error is due to an incorrect phase center correction provided with the DORIS preprocessed data resulting from a +12 cm offset in the cross-track direction that has been confirmed since. We also conclude that a 1 mm error in the cross-track offset of non-yawing sun-synchronous SPOT satellites will generate a ?6.5 mm error in the derived Z-geocenter. Other non-yawing satellites would also be affected by a similar effect whose amplitude could be easily estimated from the orbit inclination     

A model of present-day tectonic plate motions from 12 years of DORIS measurements

In the frame of the International DORIS Service (IDS), the Laboratoire d’Etudes en Géophysique et Océanographie Spatiales (LEGOS)/Collecte Localisation Satellites (CLS) Analysis Center (LCA) processes DORIS measurements from the SPOT, TOPEX/Poseidon and Envisat satellites and provides weekly station coordinates of the whole network to the IDS. Based on DORIS measurements, the horizontal and vertical velocities of 57 DORIS sites are computed. The 3D positions and velocities of the stations with linear motion are estimated simultaneously from the 12-year (1993–2004) combined normal equation matrix. We include 35 DORIS sites assumed to be located in the stable zones of 9 tectonic plates. For the motion of these plates, we propose a model (LCAVEL-1) of angular velocities in the ITRF2000 reference frame. Based on external comparison with the most recent global plate models (PB2002, REVEL, GSRM-1 and APKIM2000) and on internal analysis, we estimate an average velocity error of the DORIS solution of less than 3 mm/year. The LCAVEL-1 model presents new insights of the Somalia/Nubia pair of plates, as the DORIS technique has the advantage of having a few stations located on those two plates. We also computed (and provide in this article) the horizontal motion of the sites located close to plate boundaries or in the deformation zones defined in contemporary models. These computations could be used in further analysis for these particular regions of the Earth not moving as rigid plates.     

DORIS and the Determination of the Earth’s Polar Motion

DORIS (Détermination d’Orbite et Radiopositionnement Intégrés par Satellite) is a system used for precise orbit determination (POD) and ground-station positioning. It has been implemented on-board various satellites: the SPOT (Système pour l’Observation de la Terre) remote sensing satellites SPOT-2, SPOT-3, SPOT-4, SPOT-5, TOPEX/Poseidon and more recently on its successors Jason-1 and ENVISAT. DORIS is also a terrestrial positioning system that has found many applications in geophysics and geodesy; in particular, it contributes to the realization of the International Terrestrial Reference Frame, ITRF2000 and the forthcoming ITRF2005. Although not its primary objective, DORIS can bring information on Earth orientation monitoring, mainly polar motion and length of day (LOD) variations that complement other astrogeodetic techniques. In this paper, we have analyzed various recent polar motion solutions derived from independent analysis centers using different software packages and applying various analysis strategies. Comparisons of these solutions to the International Earth Rotation and Reference Systems Service (IERS) C04 solution are performed. Depending on the solutions, the accuracy of DORIS polar components are in the range of 0.5–1 mas corresponding to a few centimeters on the Earth’s surface. This is approximately ten times larger than results derived from GPS, which are typically 0.06 mas in both components. This does not allow DORIS results to be taken into account in the IERS–EOP combinations. A gain in the precision could come from technical improvements to the DORIS system, in addition to improvement of the orbit, tropospheric, ionospheric and Earth gravity field modeling.     

Plate kinematics of Nubia–Somalia using a combined DORIS and GPS solution

We have used up to 12 years of data to assess DORIS performance for geodynamics applications. We first examine the noise characteristics of the DORIS time-series of weekly station coordinates to derive realistic estimates of velocity uncertainties. We find that a combination of white and flicker noise best explains the DORIS time-series noise characteristics. Second, weekly solutions produced by the Institut Géographique National/Jet Propulsion Laboratory (IGN/JPL) DORIS Analysis Centre are combined to derive a global velocity field. This solution is combined with two independent GPS solutions, including 11 sites on Nubia and 5 on the Somalia plate. The combination indicates that DORIS horizontal velocities have an average accuracy of 3 mm/year, with best-determined sites having velocity accuracy better than 1 mm/year (one-sigma levels). Using our combined velocity field, we derive an updated plate kinematics model with a focus on the Nubia–Somalia area. Including DORIS data improves the precision of the angular velocity vector for Nubia by 15%. Our proposed model provides robust bounds on the maximum opening rates along the East African Rift (4.7–6.7 mm/year). It indicates opening rates 15 and 7% slower than values predicted by NUVEL-1A for the southern Atlantic Ocean and Indian Ocean, respectively. These differences are likely to arise from the fact that NUVEL-1A considered Africa as a single non-deforming plate, while here we use a more refined approach.     

Comparison of a GPS-defined global reference frame with ITRF2000

The Global Positioning System is a constellation of 24–28 satellites, which can be used to define a global terrestrial reference frame. Daily offsets between a GPS defined frame and ITRF2000 have been estimated using more than a decade of GPS observations from 1990–2001. A linear fit to the full span of data shows agreement between the two frames at the level of –1 ppb and –0.1 ppb/year for scale, 5 mm and 0 mm/year for the X component of center of mass, –2 mm and –3 mm/year for the Y component, and 4 mm and 6 mm/year for the Z component. GPS is a viable tool for defining the global reference frame either alone, or in combination with other geodetic techniques. Electronic Publication     

Evaluation of the ITRF2008 GPS vertical velocities using satellite antenna z-offsets

We develop a method to evaluate the terrestrial reference frame (TRF) scale rate error using Global Positioning System (GPS) satellite antenna phase center offset (APCO) parameters and apply it to ITRF2008. We search for the TRF in which z-APCO parameters have the smallest drift. In order to provide realistic error bars for the z-APCO drifts, we pay attention to model periodic variations and auto-correlated noise processes in the z-APCO time series. We will show that the GPS scale rate with respect to a frame is, as a first approximation, proportional to the estimated mean z-APCO trend if that frame is used to constrain station positions. Thus, an ITRF2008 scale rate error between ?0.27 and ?0.06 mm/yr depending on the GPS analysis center can be estimated, which demonstrates the high quality of the newly constructed ITRF2008. We will also demonstrate that the traditional estimates of the GPS scale rate from 7-parameter similarity transformations are consistent with our newly derived GPS scale rates with respect to ITRF2008 within two sigmas. We find using International GNSS Service (IGS) products that the traditional approach is relevant for scale rate determination even if some of the z-APCO values supplied by the IGS were not simultaneously calibrated. As the scale rate is related to the accuracy of vertical velocities, our estimates supply a conservative evaluation that can be used for error budget computation.     

Analysis of elevation changes in relation to surface characteristics for six glaciers in Northern Labrador,Canada using advanced space-borne thermal emission and reflection radiometer imagery

Advanced space-borne thermal emission and reflection radiometer imagery and Digital Elevation Models were used to analyse surface elevation changes of six glaciers in Northern Labrador. Results indicate an average surface thinning of0.94 ± 0.49 m y^?1 (water equivalent) between 2000 and 2009. Three glaciers had an average elevation change of ?1.16 ± 0.55 m y^?1 (water equivalent) whichis three times the thinning rate found in a study from 1981 to 1983 ?0.36 ± 0.10 m y^?1 water equivalent). Analysis of surface characteristics in relation to elevation changes shows expected results of rapid thinning in bare ice areas and near zero change in accumulation areas. Debris covered areas of three glaciers show expected results of moderate thinning, but three other glaciers indicate high rates of thinning. Variability in thinning rates suggests possible influences in the type ofdebris and/or variations in climate such as increased rainfall.     

First results of DORIS data analysis at Geodetic Observatory Pecný

In a cooperation between the Astronomical Institute, University of Bern (AIUB), the Geodetic Observatory Pecný (GOPE), and the Institut Géographique National (IGN), DORIS data analysis capabilities were implemented into a development version of the Bernese GPS software. The DORIS Doppler observables are reformulated such that they are similar to global navigation satellite system (GNSS) carrier-phase observations, allowing the use of the same observation models and algorithms as for GNSS carrier-phase data analysis with only minor software modifications. As such, the same algorithms may be used to process DORIS carrier-phase observations. First results from the analysis of 3 weeks of DORIS data (September 2004, five DORIS-equipped satellites) at GOPE are promising and are presented here. They include the comparison of station coordinates with coordinate estimates derived by the Laboratoire d’Etudes en Géophysique et Océanographie Spatiale/Collecte Localisation Satellites analysis centre (LCA) and the Institut Géographique National/Jet Propulsion Laboratory (IGN/JPL), and the comparison of Earth orientation parameters (EOPs) with the International Earth Rotation and Reference Frames Service (IERS) C04 model. The modified Bernese results are of a slightly lower, but comparable, quality than corresponding solutions routinely computed within the IDS (International DORIS Service). The weekly coordinate repeatability RMS is of the order of 2–3 cm for each 3D station coordinate. Comparison with corresponding estimates of station coordinates from current IDS analysis centers demonstrates similar precision. Daily pole component estimates show a mean difference from IERS-C04 of 0.6  mas in X _p and  ? 0.5  mas in Y _p and a RMS of 0.8  mas in X _p and 0.9  mas in Y _p (mean removed). An automatic analysis procedure is under development at GOPE, and routine DORIS data processing will be implemented in the near future.     

Present-day crustal deformation in China relative to ITRF97 kinematic plate model

 A global plate motion model is established based on the ITRF97 velocity fields and geological model NUVEL1. Sub-plate models are estimated by using the velocity fields derived from 45 global positioning system (GPS) sites under the ITRF97 reference frame in China. Comparisons between space geodesy and geological models are given. It is found that the Euler vector of the AFRC–EURA pair has an obvious discrepancy between space geodetic and geological models. The motion patterns of tectonic blocks in China predicted by GPS are consistent with those of geological data on the whole. Received: 9 November 2000 / Accepted: 17 September 2001     

Frontal recession of Parkachik Glacier between 1971-2015, Zanskar Himalaya using remote sensing and field data

Parkachik Glacier is located in the Suru sub-basin of the Upper Indus River, Zanskar Himalaya. The Glacier has been analysed using Corona KH-4B (1971), Landsat-TM (1999), field survey (2015), Google Earth^TM (2015) and ASTER GDEM (2015) for frontal recession and area changes. Overall, from 1971 to 2015, the Glacier has retreated by 127 ± 0.09 m i.e. (0.75 ± 0.07%) at a rate of 2.9 ± 0.004 ma^?1 with a simultaneous decrease in area from 49.5 to 48.8 km² i.e. 740 ± 0.7 m² (1.5 ± 0.09%) at a rate of 74 ± 0.7 m²a^?1. However, during recent decade (1999–2015), the rate of glacier recession of 3.9 ± 0.004 ma^?1 with a corresponding area loss of 500 ± 0.74m² (1 ± 0.1%) was higher than the retreat rate of 2.3 ± 0.001 ma^?1 and an area loss of 240 ± 0.02m² (0.48 ± 0.08%) during 1971–1999. In the field, the evidences of glacier recession are present in the form of separated dead ice blocks from the main Glacier, recessional dumps/moraines, active ice calving activity and a small proglacial pond/lake at the terminus/snout of the Glacier. However, the recession over the studied period has been very slow and is controlled by its topographic configuration, particularly the large altitudinal range (6030–3620 m), almost northerly aspect and steep slope (average ~ 30°).