Seasonal and secular positional variations at eight co-located GPS and VLBI stations

Time series of daily position solutions at eight co-located GPS and VLBI stations are used to assess the frequency features in the solutions over various time-scales. This study shows that there are seasonal and inter-annual signals in all three coordinate components of the GPS and VLBI solutions. The power and frequency of the signals vary with time, the station considered and the coordinate components, and between the GPS and VLBI solutions. In general, the magnitudes of the signals in the horizontal coordinate components (latitude and longitude) are weaker than those in the height component. The weighted means of the estimated annual amplitudes from the eight GPS stations are, respectively, 1.0, 0.8 and 3.6 mm for the latitude, longitude and height components, and are, respectively, 1.5, 0.7 and 2.2 mm for the VLBI solutions. The phases of the annual signals estimated from the GPS and VLBI solutions are consistent for most of the co-located stations. The seasonal signals estimated from the VLBI solutions are, in general, more stable than those estimated from the GPS solutions. Fluctuations at inter-annual time-scales are also found in the series. The inter-annual fluctuations are up to ∼5 mm for the latitude and longitude components, and up to ∼10 mm for the height component. The effects of the seasonal and inter-annual variations on the estimated linear rates of movement of the stations are also evaluated.     

利用GPS和GRACE研究澳大利亚地壳垂向季节性变化

为探究重力场恢复与气候实验(gravity recovery and climate experiment,GRACE)卫星与全球定位系统(global positioning system,GPS)两种独立技术获取的因陆地水储量变化引起的地壳垂向季节性位移的一致性,选取澳大利亚27个GPS站点5～10 a的高程时间序...     

Anomalous harmonics in the spectra of GPS position estimates

Prior studies of the power spectra of GPS position time series have found pervasive seasonal signals against a power-law background of flicker noise plus white noise. Dong et al. (2002) estimated that less than half the observed GPS seasonal power can be explained by redistributions of geophysical fluid mass loads. Much of the residual variation is probably caused by unidentified GPS technique errors and analysis artifacts. Among possible mechanisms, Penna and Stewart (2003) have shown how unmodeled analysis errors at tidal frequencies (near 12- and 24-hour periods) can be aliased to longer periods very efficiently. Signals near fortnightly, semiannual, and annual periods are expected to be most seriously affected. We have examined spectra for the 167 sites of the International GNSS (Global Navigation Satellite Systems) Service (IGS) network having more than 200 weekly measurements during 1996.0–2006.0. The non-linear residuals of the weekly IGS solutions that were included in ITRF2005, the latest version of the International Terrestrial Reference Frame (ITRF), have been used. To improve the detection of common-mode signals, the normalized spectra of all sites have been stacked, then boxcar smoothed for each local north (N), east (E), and height (H) component. The stacked, smoothed spectra are very similar for all three components. Peaks are evident at harmonics of about 1 cycle per year (cpy) up to at least 6 cpy, but the peaks are not all at strictly 1.0 cpy intervals. Based on the 6th harmonic of the N spectrum, which is among the sharpest and largest, and assuming a linear overtone model, then a common fundamental of 1.040 ± 0.008 cpy can explain all peaks well, together with the expected annual and semiannual signals. A flicker noise power-law continuum describes the background spectrum down to periods of a few months, after which the residuals become whiter. Similar sub-seasonal tones are not apparent in the residuals of available satellite laser ranging (SLR) and very long baseline interferometry (VLBI) sites, which are both an order of magnitude less numerous and dominated by white noise. There is weak evidence for a few isolated peaks near 1 cpy harmonics in the spectra of geophysical loadings, but these are much noisier than for GPS positions. Alternative explanations related to the GPS technique are suggested by the close coincidence of the period of the 1.040 cpy frequency, about 351.2 days, to the “GPS year”; i.e., the interval required for the constellation to repeat its inertial orientation with respect to the sun. This could indicate that the harmonics are a type of systematic error related to the satellite orbits. Mechanisms could involve orbit modeling defects or aliasing of site-dependent positioning biases modulated by the varying satellite geometry.     

Modeling thermal deformation of VLBI antennas with a new temperature model

Temperature variations at very long baseline interferometry (VLBI) sites cause thermal deformations of the VLBI antennas and corresponding displacements of the VLBI reference points. The thermal deformation effects typically contain seasonal and daily signatures. The amplitudes of the annual vertical motion of the antenna reference point can reach several millimeters, depending on the design of the antenna structure, on the material, and on the environmental effects such as global station position, station height and climatology effects. Simple methods to correct this effect use the difference of the environmental temperature with respect to a defined reference temperature, the antenna dimensions, the elevation of the antenna, the material of antenna structure. Applying these simple models for thermal deformation in the VLBI data analysis improves the baseline length repeatability by 3.5%. A comparison of these simple models with local thermal deformation measurements at the antennas in Onsala and Wettzell show that the local measurements and the modeled corrections agree well when the temperature of the antenna structure is used, but agree less good when the surrounding air temperatures are used. To overcome this problem we present a method to model temperature penetration into the antenna structures, that allows to model thermal deformation effects that agree with the observed vertical deformation of the Onsala and Wettzell radio telescopes with a root mean square deviation of 0.07 and 0.13 mm, respectively. Possible implementations in the VLBI analysis are presented, and the definition of an adequate reference temperature is discussed.     

Quality assessment of GPS reprocessed terrestrial reference frame

The International GNSS Service (IGS) contributes to the construction of the International Terrestrial Reference Frame (ITRF) by submitting time series of station positions and Earth Rotation Parameters (ERP). For the first time, its submission to the ITRF2008 construction is based on a combination of entirely reprocessed GPS solutions delivered by 11 Analysis Centers (ACs). We analyze the IGS submission and four of the individual AC contributions in terms of the GNSS frame origin and scale, station position repeatability and time series seasonal variations. We show here that the GPS Terrestrial Reference Frame (TRF) origin is consistent with Satellite laser Ranging (SLR) at the centimeter level with a drift lower than 1 mm/year. Although the scale drift compared to Very Long baseline Interferometry (VLBI) and SLR mean scale is smaller than 0.4 mm/year, we think that it would be premature to use that information in the ITRF scale definition due to its strong dependence on the GPS satellite and ground antenna phase center variations. The new position time series also show a better repeatability compared to past IGS combined products and their annual variations are shown to be more consistent with loading models. The comparison of GPS station positions and velocities to those of VLBI via local ties in co-located sites demonstrates that the IGS reprocessed solution submitted to the ITRF2008 is more reliable and precise than any of the past submissions. However, we show that some of the remaining inconsistencies between GPS and VLBI positioning may be caused by uncalibrated GNSS radomes.     

Impact of seasonal station motions on VLBI UT1 intensives results

UT1 estimates obtained from the very long baseline interferometry (VLBI) Intensives data depend on the station displacement model used during processing. In particular, because of seasonal variations, the instantaneous station position during the specific intensive session differs from the position predicted by the linear model generally used. This can cause systematic errors in UT1 Intensives results. In this paper, we first investigated the seasonal signal in the station displacements for the 5 VLBI antennas participating in UT1 Intensives observing programs, along with the 8 collocated GPS stations. It was found that a significant annual term is present in the time series for most stations, and its amplitude can reach 8 mm in the height component, and 2 mm in horizontal components. However, the annual signals found in the displacements of the collocated VLBI and GPS stations at some sites differ substantially in amplitude and phase. The semiannual harmonics are relatively small and unstable, and for most stations no prevailing signal was found in the corresponding frequency band. Then two UT1 Intensives series were computed with and without including the seasonal term found in the previous step in the station movement model. Comparison of these series has shown that neglecting the seasonal station position variations can cause a systematic error in UT1 estimates, which can exceed 1  $\upmu $ s, depending on the observing program.     

Multi-technique comparison of troposphere zenith delays and gradients during CONT08

CONT08 was a 15 days campaign of continuous Very Long Baseline Interferometry (VLBI) sessions during the second half of August 2008 carried out by the International VLBI Service for Geodesy and Astrometry (IVS). In this study, VLBI estimates of troposphere zenith total delays (ZTD) and gradients during CONT08 were compared with those derived from observations with the Global Positioning System (GPS), Doppler Orbitography and Radiopositioning Integrated by Satellite (DORIS), and water vapor radiometers (WVR) co-located with the VLBI radio telescopes. Similar geophysical models were used for the analysis of the space geodetic data, whereas the parameterization for the least-squares adjustment of the space geodetic techniques was optimized for each technique. In addition to space geodetic techniques and WVR, ZTD and gradients from numerical weather models (NWM) were used from the European Centre for Medium-Range Weather Forecasts (ECMWF) (all sites), the Japan Meteorological Agency (JMA) and Cloud Resolving Storm Simulator (CReSS) (Tsukuba), and the High Resolution Limited Area Model (HIRLAM) (European sites). Biases, standard deviations, and correlation coefficients were computed between the troposphere estimates of the various techniques for all eleven CONT08 co-located sites. ZTD from space geodetic techniques generally agree at the sub-centimetre level during CONT08, and??as expected??the best agreement is found for intra-technique comparisons: between the Vienna VLBI Software and the combined IVS solutions as well as between the Center for Orbit Determination (CODE) solution and an IGS PPP time series; both intra-technique comparisons are with standard deviations of about 3?C6?mm. The best inter space geodetic technique agreement of ZTD during CONT08 is found between the combined IVS and the IGS solutions with a mean standard deviation of about 6?mm over all sites, whereas the agreement with numerical weather models is between 6 and 20?mm. The standard deviations are generally larger at low latitude sites because of higher humidity, and the latter is also the reason why the standard deviations are larger at northern hemisphere stations during CONT08 in comparison to CONT02 which was observed in October 2002. The assessment of the troposphere gradients from the different techniques is not as clear because of different time intervals, different estimation properties, or different observables. However, the best inter-technique agreement is found between the IVS combined gradients and the GPS solutions with standard deviations between 0.2 and 0.7?mm.     

Nontectonic signals in reprocessed continuous GPS position time series of CMONOC

The nearly nine-year continuous GPS data collected since 1 March 1999 from the Crustal Motion Observation Network of China (CMONOC) were consistently analyzed. Most of the nonlinear movements in the cumulative position time series produced by CMONOC data center disappeared; and more accurate vertical terms and tectonic signals were extracted. Displacements caused by atmospheric pressure loading, nontidal ocean loading, soil moisture mass loading, and snow cover mass loading using the National Centers for Environmental Prediction (NCEP) Reanalysis I/II models and Estimation of the Circulation and Climate of the Ocean (ECCO) data can explain most of the vertical annual terms at many stations, while only parts can be explained at Lhasa and southern coastal sites, indicating that there are some deformation mechanisms that are still unknown or not modeled accurately. The remarkable differences in vertical position time series for short-baseline sites reveal that GPS stations can be greatly affected by local factors; and attention should be paid when explaining observed GPS velocity vectors.     

Nontectonic signals in reprocessed continuous GPS position time series of CMONOC

The nearly nine-year continuous GPS data collected since 1 March 1999 from the Crustal Motion Observation Network of China(CMONOC) were consistently analyzed.Most of the nonlinear movements in the cumulative position time series pro-duced by CMONOC data center disappeared;and more accurate vertical terms and tectonic signals were extracted.Displacements caused by atmospheric pressure loading,nontidal ocean loading,soil moisture mass loading,and snow cover mass loading using the National Centers for Environmental Prediction(NCEP) Reanalysis I/II models and Estimation of the Circulation and Climate of the Ocean(ECCO) data can explain most of the vertical annual terms at many stations,while only parts can be explained at Lhasa and southern coastal sites,indicating that there are some deformation mechanisms that are still unknown or not modeled accurately.The remarkable differences in vertical position time series for short-baseline sites reveal that GPS stations can be greatly affected by lo-cal factors;and attention should be paid when explaining observed GPS velocity vectors.     

A comparison of GPS, VLBI and model estimates of ocean tide loading displacements

In recent years, ocean tide loading displacements (OTLD) have been measured using the Global Positioning System (GPS) and Very Long Baseline Interferometry (VLBI). This study assesses the accuracy of GPS measurements of OTLD by comparison with VLBI measurements and estimates derived from numerical ocean tide models. A daily precise point positioning (PPP) analysis was carried out on ∼11 years of GPS data for each of 25 sites that have previous OTLD estimates based on data from co-located VLBI sites. Ambiguities were fixed to integer values where possible. The resulting daily estimates of OTLD, at eight principal diurnal and semi-diurnal tidal frequencies, were combined to give GPS measurements of OTLD at each site. The 3D GPS and VLBI measurements of OTLD were compared with estimates computed (by convolution with Green’s functions) from five modern ocean tide models (CSR4.0, FES2004, GOT00.2, NAO99b and TPXO6.2). The GPS/model agreement is shown to be similar to the VLBI/model agreement. In the important radial direction, the GPS/model misfit is shown to be smaller than the VLBI/model misfit for seven of the eight tidal constituents; the exception being the K2 constituent. Fixing of GPS carrier-phase ambiguities to integer values resulted in a marginal improvement to the GPS/model agreement. Statistically, it is shown there is no significance to the difference between the fit of the GPS and VLBI measurements of OTLD to modelled values. Equally, differences in fit of either the complete set of GPS or VLBI estimates to the five sets of model-derived values cannot be identified with statistical significance. It is thus concluded that, overall, we cannot distinguish between GPS and VLBI measurements of OTLD, and that at the global scale, present ocean tide models are accurate to within the current measurement noise of these techniques.     

Transformation between SLR/VLBI and WGS-84 reference frames

In geodetic and geophysical applications of GPS, it is important to realize the ephemerides of the GPS satellites and the coordinates of station positions in a consistent reference system. At present, more than one reference system is being used by various GPS users depending on their specific applications. The WGS-84 and various reference frames based on satellite laser ranging (SLR), very long baseline interferometry (VLBI), or a combination of SLR and VLBI are the most commonly used in high precision geophysical applications. The WGS-84 is widely used in applications which rely on the GPS broadcast ephemeris. Station coordinates estimated in one system may have to be transformed to another for further use or for evaluation/comparison purposes. This paper presents a seven-parameter transformation between the WGS-84 and SLR/VLBI reference frames. The GPS double-differenced phase measurements for two consecutive weeks from a set of five Defense Mapping Agency (DMA) sites (defined in the WGS-84 frame) and from an augmented set of fifteen CIGNET sites (defined in the SLR/VLBI frame) were processed in a least squares estimation scheme to determine station coordinates, from which the transformation parameters were determined. A scale difference of about 0.2 ppm and an orientation difference in longitude of about 31 milliarcseconds were found to be the only parameters of significance between the adopted SLR/VLBI and the WGS-84 frames. Transformation between WGS-84 and the ITRF90 is also included, in which the scale difference is the same as before but the longitude rotation is about 16 mas.     

基于IGS连续跟踪站的GPS高程时间序列分析

对全球分布的IGS连续跟踪站GPS高程分量时间序列,进行了小波多分辨率分析,发现各台站高程时间序列大都存在着季节性变化,这种以半周年项和周年项为主的季节性变化是由多种地球物理现象引起的,我们从多个方面分析了原因,并着重分析了大气压强的影响。利用分析的结果讨论了地壳垂直变化的规律,得出地球存在以半年和周年为周期的整体性扩张与收缩振荡运动,同时南北半球的变化规律存在差异的结论。     

Assessment of second- and third-order ionospheric effects on regional networks: case study in China with longer CMONOC GPS coordinate time series

Higher-order ionospheric (HOI) delays are one of the principal technique-specific error sources in precise global positioning system analysis and have been proposed to become a standard part of precise GPS data processing. In this research, we apply HOI delay corrections to the Crustal Movement Observation Network of China’s (CMONOC) data processing (from January 2000 to December 2013) and furnish quantitative results for the effects of HOI on CMONOC coordinate time series. The results for both a regional reference frame and global reference frame are analyzed and compared to clarify the HOI effects on the CMONOC network. We find that HOI corrections can effectively reduce the semi-annual signals in the northern and vertical components. For sites with lower semi-annual amplitudes, the average decrease in magnitude can reach 30 and 10 % for the northern and vertical components, respectively. The noise amplitudes with HOI corrections and those without HOI corrections are not significantly different. Generally, the HOI effects on CMONOC networks in a global reference frame are less obvious than the results in the regional reference frame, probably because the HOI-induced errors are smaller in comparison to the higher noise levels seen when using a global reference frame. Furthermore, we investigate the combined contributions of environmental loading and HOI effects on the CMONOC stations. The largest loading effects on the vertical displacement are found in the mid- to high-latitude areas. The weighted root mean square differences between the corrected and original weekly GPS height time series of the loading model indicate that the mass loading adequately reduced the scatter on the CMONOC height time series, whereas the results in the global reference frame showed better agreements between the GPS coordinate time series and the environmental loading. When combining the effects of environmental loading and HOI corrections, the results with the HOI corrections reduced the scatter on the observed GPS height coordinates better than the height when estimated without HOI corrections, and the combined solutions in the regional reference frame indicate more preferred improvements. Therefore, regional reference frames are recommended to investigate the HOI effects on regional networks.     

Continental hydrology loading observed by VLBI measurements

Variations in continental water storage lead to loading deformation of the crust with typical peak-to-peak variations at very long baseline interferometry (VLBI) sites of 3–15 mm in the vertical component and 1–2 mm in the horizontal component. The hydrology signal at VLBI sites has annual and semi-annual components and clear interannual variations. We have calculated the hydrology loading series using mass loading distributions derived from the global land data assimilation system (GLDAS) hydrology model and alternatively from a global grid of equal-area gravity recovery and climate experiment (GRACE) mascons. In the analysis of the two weekly VLBI 24-h R1 and R4 network sessions from 2003 to 2010 the baseline length repeatabilities are reduced in 79 % (80 %) of baselines when GLDAS (GRACE) loading corrections are applied. Site vertical coordinate repeatabilities are reduced in about 80 % of the sites when either GLDAS or GRACE loading is used. In the horizontal components, reduction occurs in 70–80 % of the sites. Estimates of the annual site vertical amplitudes were reduced for 16 out of 18 sites if either loading series was applied. We estimated loading admittance factors for each site and found that the average admittances were 1.01 \(\pm \) 0.05 for GRACE and 1.39 \(\pm \) 0.07 for GLDAS. The standard deviations of the GRACE admittances and GLDAS admittances were 0.31 and 0.68, respectively. For sites that have been observed in a set of sufficiently temporally dense daily sessions, the average correlation between VLBI vertical monthly averaged series and GLDAS or GRACE loading series was 0.47 and 0.43, respectively.     

Tropospheric parameters: combination studies based on homogeneous VLBI and GPS data

The combination of tropospheric parameters derived from different space-geodetic techniques has not been of large interest in geodesy so far. However, due to the high correlation between station coordinates and tropospheric parameters, the latter should not be neglected in combinations. This paper deals with the comparison and combination of tropospheric parameters derived from global positioning system (GPS) and very long baseline interferometry (VLBI) observations stemming from a 15-day campaign of continuous VLBI observations in 2002 (CONT02). The observation data of both techniques were processed homogeneously to avoid systematic differences between the solutions. We compared the tropospheric estimates of GPS and VLBI at eight co-location sites and found a very good agreement in the temporal behavior of the tropospheric zenith path delays (ZPD), reflected by correlation factors up to 0.98. Following this, a combination of the tropospheric parameters was performed. We demonstrate that the combination of tropospheric parameters leads to a stabilization of combined station networks. This becomes visible in the improvement of the repeatabilities of the station height components. Furthermore, the potential use of independent data from water vapor radiometers (WVRs) to validate space-technique-derived tropospheric parameters was investigated. Correlation coefficients of 0.95 or better were estimated between the tropospheric parameters of WVR and GPS or VLBI. Additionally, the utility of the tropospheric parameters for validation of local tie vectors was investigated. Both tropospheric zenith delays and tropospheric gradients were found to be very suitable to validate the height component and the horizontal components of the local tie, respectively.     

Empirical modelling of site-specific errors in continuous GPS data

Continuous global positioning system (GPS) stations propagate biases and spurious signals into the derived parameter time series when the measurements are subject to site-specific effects, such as multipath. This is a particular problem in the investigation of geophysical and atmospheric phenomena where signals may be small in magnitude. A methodology to remove these erroneous signals from long-term time series will significantly increase the usefulness of the derived time series. This work provides the theoretical basis for use of an empirical site model (ESM) derived from post-fit phase residuals to mitigate unmodelled site-specific errors. Additionally, we also investigate the effectiveness of applying an ESM to a regional GPS network and a short baseline solution. Under most observing scenarios, we show that the ESM approach is predicted to improve the precision and accuracy of the site coordinates. However, it is important to note that we found some scenarios where the ESM can introduce a bias. For instance, when the antenna is mounted close to the ground. In this scenario, for a short baseline, we observed the introduction of a 4-mm bias in height. Use of an ESM for the same short baseline with an uncalibrated radome substantially improves the results by removing a large bias of over 10 mm in height. Similarly, application of an ESM derived from historic data yields similar improvements. This demonstrates that the ESM can be a powerful tool when applied to appropriate site-specific configurations and could potentially be implemented in routine GPS analysis for a broad range of applications.     

Strategies to mitigate aliasing of loading signals while estimating GPS frame parameters

Although GNSS techniques are theoretically sensitive to the Earth center of mass, it is often preferable to remove intrinsic origin and scale information from the estimated station positions since they are known to be affected by systematic errors. This is usually done by estimating the parameters of a linearized similarity transformation which relates the quasi-instantaneous frames to a long-term frame such as the International Terrestrial Reference Frame (ITRF). It is well known that non-linear station motions can partially alias into these parameters. We discuss in this paper some procedures that may allow reducing these aliasing effects in the case of the GPS techniques. The options include the use of well-distributed sub-networks for the frame transformation estimation, the use of site loading corrections, a modification of the stochastic model by downweighting heights, or the joint estimation of the low degrees of the deformation field. We confirm that the standard approach consisting of estimating the transformation over the whole network is particularly harmful for the loading signals if the network is not well distributed. Downweighting the height component, using a uniform sub-network, or estimating the deformation field perform similarly in drastically reducing the amplitude of the aliasing effect. The application of these methods to reprocessed GPS terrestrial frames permits an assessment of the level of agreement between GPS and our loading model, which is found to be about 1.5 mm WRMS in height and 0.8 mm WRMS in the horizontal at the annual frequency. Aliased loading signals are not the main source of discrepancies between loading displacement models and GPS position time series.     

Vertical deformations from homogeneously processed GRACE and global GPS long-term series

Temporal variations in the geographic distribution of surface mass cause surface displacements. Surface displacements derived from GRACE gravity field coefficient time series also should be observed in GPS coordinate time series, if both time series are sufficiently free of systematic errors. A successful validation can be an important contribution to climate change research, as the biggest contributors to mass variability in the system Earth include the movement of oceanic, atmospheric, and continental water and ice. In our analysis, we find that if the signals are larger than their precision, both geodetic sensor systems see common signals for almost all the 115 stations surveyed. Almost 80% of the stations have their signal WRMS decreased, when we subtract monthly GRACE surface displacements from those observed by GPS data. Almost all other stations are on ocean islands or small peninsulas, where the physically expected loading signals are very small. For a fair comparison, the data (79 months from September 2002 to April 2009) had to be treated appropriately: the GPS data were completely reprocessed with state-of-the-art models. We used an objective cluster analysis to identify and eliminate stations, where local effects or technical artifacts dominated the signals. In addition, it was necessary for both sets of results to be expressed in equivalent reference frames, meaning that net translations between the GPS and GRACE data sets had to be treated adequately. These data sets are then compared and statistically analyzed: we determine the stability (precision) of GRACE-derived, monthly vertical deformation data to be ~1.2 mm, using the data from three GRACE processing centers. We statistically analyze the mean annual signals, computed from the GPS and GRACE series. There is a detailed discussion of the results for five overall representative stations, in order to help the reader to link the displayed criteria of similarity to real data. A series of tests were performed with the goal of explaining the remaining GPS–GRACE residuals.     

Multi-technique comparisons of 10?years of wet delay estimates on the west coast of Sweden

We present comparisons of 10-year-long time series of the atmospheric zenith wet delay (ZWD), estimated using the global positioning system (GPS), geodetic very long baseline interferometry (VLBI), a water vapour radiometer (WVR), radiosonde (RS) observations, and the reanalysis product of the European Centre for Medium-Range Weather Forecasts (ECMWF). To compare the data sets with each other, a Gaussian filter is applied. The results from 10 GPS–RS comparisons using sites in Sweden and Finland show that the full width at half maximum at which the standard deviation (SD) is a minimum increases with the distance between each pair. Comparisons between three co-located techniques (GPS, VLBI, and WVR) result in mean values of the ZWD differences at a level of a few millimetres and SD of less than 7?mm. The best agreement is seen in the GPS–VLBI comparison with a mean difference of ?3.4?mm and an SD of 5.1?mm over the 10-year period. With respect to the ZWD derived from other techniques, a positive bias of up to ~7?mm is obtained for the ECMWF reanalysis product. Performing the comparisons on a monthly basis, we find that the SD including RS or ECMWF varies with the season, between 3 and 15?mm. The monthly SD between GPS and WVR does not have a seasonal signature and varies from 3 to 7?mm.     

温度变化对中国大陆三维周年位移的影响

除了地表质量重新分布外,地表温度变化是影响地表周年变化的另一重要成因。本文利用全球温度变化数据,基于三维全空间热弹性形变模型,计算温度变化在中国大陆引起的地表热弹性形变,并讨论它对中国大陆三维周年位移的影响。结果表明,温度变化引起的地表周年变化振幅在毫米量级。中国大陆构造环境监测网络(简称:陆态网络)GPS台站受地表温度变化影响最大的台站是HLAR(海拉尔),东向、北向以及垂向的周年振幅矢量和为~2.293mm;影响最小的台站是HIYS(永暑礁),东向、北向以及垂向的周年振幅矢量和为~0.177mm。为了说明温度变化对地表周年形变的影响,本文联合GRACE以及物质负荷模型(MODEL)研究中国大陆地表三维周年位移。考虑温度变化后的MODEL、GRACE获取的地表周年形变在东向、北向、垂向的周年信号分别改进了6%、6%、2%;16%、5%、15%。结果表明,温度变化是物质负荷以外引起大陆地表形变的重要因素。