利用VLBI观测对ERP参数和月地距变化相关性的研究

处理了2000-01～2012-06的VLBI观测数据,提取了ERP地球自转参数信息,通过JPL的DE405星历计算得到了2000-01～2012-06的月地距时间序列;消除日长变化、极移和月地距时间序列趋势项,采用频谱分析的方法求得它们的周期;根据解算的周期推测日长变化和日地距可能存在相关性,并验证了推测.     

利用1985～2010年VLBI数据进行日长变化分析

处理了1985～2010年的VLBI观测数据,提取了日长变化信息,分析了其内符合精度和外符合精度;通过傅立叶分析以及小波分析得出了日长变化的半月周期、月周期、半年周期、周年周期以及更长周期;利用解算出的日长变化信息绘制日长变化积分曲线,分析日长的季节性变化和长期变化,并分析了日长变化因素。     

TOPEX/Poseidon卫星监测博斯腾湖水位变化及其与NINO3 SST的相关性分析

利用TOPEX/Poseidon(T/P)卫星10年（1992~2002）监测博斯腾湖水位变化的GDR数据,经过数据的地球物理改正、编辑、约化和滤波后,分别采用小波分析和傅里叶分析得到了博斯腾湖水位变化的平均周期为1034d,363d,182d和47d。季节性变化是湖泊水位变化的主项,利用最小二乘法得到了周年振幅为0.7196m,半年振幅为0.5322m,因此季节性变化的主周期是周年,次周期是半年。博斯腾湖水位变化变化趋势为0.2012m/a,说明在全球气候变暖趋势下,博斯腾湖水位在升高。由博斯腾湖月水位变化的时间序列与NINO3 SST进行相关性分析,发现在NINO3 SST 发生6个月后,博斯腾湖月水位变化序列与NINO3 SST数据相关性最大(-0.2784),表明厄尔尼诺事件对我国西部气候影响的具有一定的滞后性。     

VLBI和GPS观测联合解算地球自转参数和日长变化

利用GAMIT软件处理了2005年9月、2006年1月和2月每天17个IGS国际GNSS服务(International GNSS Service,IGS)站的GPS观测数据,获得了地球自转参数(Earth rotation parameter,ERP)和日长变化(ΔLOD),并与IGS综合解进行了对比。利用OCCAM 6.2软件处理了相同时间内的甚长基线干涉测量(very long baseline interferometry,VLBI)观测数据,将解算结果与国际VLBI服务(International VLBI Service,IVS)结果进行了对比。分别采用基于内符合精度和IERS 08C04序列的定权方法对VLBI解算结果与GPS解算结果进行了加权平均,获得了VLBI和GPS技术对ERP和ΔLOD的协议结果。研究结果表明,VLBI解算结果与采用的插值方法有较大的关系,基于IERS 08C04序列的加权平均方法达到了利用VLBI解算结果对GPS解算结果进行修正的目的。VLBI和GPS技术的联合弥补了VLBI观测数据密度不够和GPS解算结果不稳定的缺陷,使解算结果的稳定性和可靠性有所提高。     

2012-2016年的VLBI数据解算及分析

针对甚长基线干涉测量技术以其极高的测角精度和稳定性在解算地球定向参数和维持地球参考框架方面发挥着不可替代的作用,但是由于其设备庞大昂贵,因此测站较少且观测不连续的问题,该文首先利用VieVS软件解算了2012—2016年的甚长基线干涉测量数据,获得地球定向参数。提出了采用极坐标插值法对极移参数进行插值,以获取连续的极移参数解。对插值后的极移时间序列进行频谱分析,得到较为明显的周年、半周年等周期项,验证了本文采用插值方法的可行性和可靠性。     

利用2008～2009年VLBI数据进行日长变化的研究

处理了2008-01～2009-03的VLBI观测数据,提取了日长变化信息,并与IVS网站的结果进行了对比;计算了日长变化的月周期长度和半月周期长度;利用近几年VLBI数据中的日长变化信息绘制日长变化曲线,分析了日长的季节性变化及引起日长变化的因素。     

基于VLBI与BDS联合解算ERP的算法研究

利用PANDA软件解算2016年第61~91天的MGEX(Multi-GNSS Experiment)服务站的北斗数据,获得地球自转参数(ERP)。利用VieVS2.2软件处理了同时段的甚长基线干涉测量(VLBI)数据。采用基于IERS 08C04序列的定权方法对BDS和VLBI的解算结果进行加权平均,得到综合的ERP值。结果表明,与IERS比较,极移在X方向差值的RMS为0.249 mas,Y方向差值的RMS为0.296 mas,UT1-UTC差值的RMS为0.053 ms.利用BDS与VLBI数据对ERP参数进行联合解算,弥补了BDS解算结果不稳定和VLBI观测不连续的缺陷,使解算结果的稳定性和可靠性均有所提高。      

基于时变地球主惯性矩的三轴地球的自由Euler运动

在时变地球主惯性矩情形,采用数值方法求解了自由Euler动力学方程和运动学方程.结果表明,地球存在周期约为304.5 d的自由章动;旋转速度的三个分量均出现周期性的增大和减小;地球主惯性矩A、B、C的时变性导致EuIer周期产生复杂的波动,特别是ω3(自转速率)存在周期为22 a、14 a、8 a、4 a、2 a、周年以及更短周期的波动,这表明A、B、C的时变性导致了10 a尺度、周年以及更短周期的日长变化.     

利用重力和VLBI技术检测地球液态核的动力学效应

利用最新的武汉地区高精度重力和全球高精度甚长基线干涉（very long baseline interferometry,VLBI）测量技术观测确定了地球液态地核的自由核章动本征参数（包括本征周期和品质因子等）,同时获得了液核顶部的粘滞系数和液核的真实动力学椭率,使用Monte Carlo方法研究了本征周期的不确定度。结果表明,分别利用重力和VLBI确定的地球自由核章动本征周期值为430 sd（sd指恒星日）左右,两者间的差异为1.47 sd,粘滞系数为1 028 Pa·s。讨论了导致两者差异的主要原因,并将结果与国内外同类研究结果作了对比。     

利用RADIOASTRON模拟估计的EOP对相关参数精度的灵敏度分析

介绍了包括章动参数的空间VLBI数学模型,探讨了该模型用于大地测量解算的可行性;模拟了针对RADIOASTRON卫星的空间VLBI观测量及相关参数;利用所模拟的时延观测量研究分析包含章动参数的数学模型对EOP参数的可估计性,进而给空间VLBI数学模型中除EOP参数以外的参数赋予不同的先验误差,研究EOP参数对这些先验误差的灵敏度。     

Consistent realization of Celestial and Terrestrial Reference Frames

The Celestial Reference System (CRS) is currently realized only by Very Long Baseline Interferometry (VLBI) because it is the space geodetic technique that enables observations in that frame. In contrast, the Terrestrial Reference System (TRS) is realized by means of the combination of four space geodetic techniques: Global Navigation Satellite System (GNSS), VLBI, Satellite Laser Ranging (SLR), and Doppler Orbitography and Radiopositioning Integrated by Satellite. The Earth orientation parameters (EOP) are the link between the two types of systems, CRS and TRS. The EOP series of the International Earth Rotation and Reference Systems Service were combined of specifically selected series from various analysis centers. Other EOP series were generated by a simultaneous estimation together with the TRF while the CRF was fixed. Those computation approaches entail inherent inconsistencies between TRF, EOP, and CRF, also because the input data sets are different. A combined normal equation (NEQ) system, which consists of all the parameters, i.e., TRF, EOP, and CRF, would overcome such an inconsistency. In this paper, we simultaneously estimate TRF, EOP, and CRF from an inter-technique combined NEQ using the latest GNSS, VLBI, and SLR data (2005–2015). The results show that the selection of local ties is most critical to the TRF. The combination of pole coordinates is beneficial for the CRF, whereas the combination of \(\varDelta \hbox {UT1}\) results in clear rotations of the estimated CRF. However, the standard deviations of the EOP and the CRF improve by the inter-technique combination which indicates the benefits of a common estimation of all parameters. It became evident that the common determination of TRF, EOP, and CRF systematically influences future ICRF computations at the level of several \(\upmu \)as. Moreover, the CRF is influenced by up to \(50~\upmu \)as if the station coordinates and EOP are dominated by the satellite techniques.     

EOP and scale from continuous VLBI observing: CONT campaigns to future VGOS networks

Continuous (CONT) VLBI campaigns have been carried out about every 3 years since 2002. The basic idea of these campaigns is to acquire state-of-the-art VLBI data over a continuous time period of about 2 weeks to demonstrate the highest accuracy of which the current VLBI system is capable. In addition, these campaigns support scientific studies such as investigations of high-resolution Earth rotation, reference frame stability, and daily to sub-daily site motions. The size of the CONT networks and the observing data rate have increased steadily since 1994. Performance of these networks based on reference frame scale precision and polar motion/LOD comparison with global navigation satellite system (GNSS) earth orientation parameters (EOP) has been substantially better than the weekly operational R1 and R4 series. The precisions of CONT EOP and scale have improved by more than a factor of two since 2002. Polar motion precision based on the WRMS difference between VLBI and GNSS for the most recent CONT campaigns is at the 30 \(\upmu \)as level, which is comparable to that of GNSS. The CONT campaigns are a natural precursor to the planned future VLBI observing networks, which are expected to observe continuously. We compare the performance of the most recent CONT campaigns in 2011 and 2014 with the expected performance of the future VLBI global observing system network using simulations. These simulations indicate that the expected future precision of scale and EOP will be at least 3 times better than the current CONT precision.     

Effect of post-seismic deformation on earth orientation parameter estimates from VLBI observations: a case study at Gilcreek, Alaska

Earth orientation parameters (EOPs) provide a link between the International Celestial Reference Frame (ICRF) and the International Terrestrial Reference Frame (ITRF). Natural geodynamic processes, such as earthquakes, can cause the motion of stations to become discontinuous and/or non-linear, thereby corrupting the EOP estimates if the sites are assumed to move linearly. The VLBI antenna at the Gilcreek Geophysical Observatory has undergone non-linear, post-seismic motion as a result of the M_w=7.9 Denali earthquake in November 2002, yet some VLBI analysts have adopted co-seismic offsets and a linear velocity model to represent the motion of the site after the earthquake. Ignoring the effects of the Denali earthquake leads to error on the order of 300–600 μas for the EOP, while modelling the post-seismic motion of Gilcreek with a linear velocity generates errors of 20–50 μas. Only by modelling the site motion with a non-linear function is the same level of accuracy of EOP estimates maintained. The effect of post-seismic motion on EOP estimates derived from the International VLBI Service IVS-R1 and IVS-R4 networks are not the same, although changes in network geometries and equipment improvements have probably affected the estimates more significantly than the earthquake-induced deformation at Gilcreek.     

International VLBI Service for Geodesy and Astrometry

The International VLBI Service for Geodesy and Astrometry (IVS) regularly produces high-quality Earth orientation parameters from observing sessions employing extensive networks or individual baselines. The master schedule is designed according to the telescope days committed by the stations and by the need for dense sampling of the Earth orientation parameters (EOP). In the pre-2011 era, the network constellations with their number of telescopes participating were limited by the playback and baseline capabilities of the hardware (Mark4) correlators. This limitation was overcome by the advent of software correlators, which can now accommodate many more playback units in a flexible configuration. In this paper, we describe the current operations of the IVS with special emphasis on the quality of the polar motion results since these are the only EOP components which can be validated against independent benchmarks. The polar motion results provided by the IVS have improved continuously over the years, now providing an agreement with IGS results at the level of 20–25 \(\upmu \)as in a WRMS sense. At the end of the paper, an outlook is given for the realization of the VLBI Global Observing System.     

On comparison of the Earth orientation parameters obtained from different VLBI networks and observing programs

In this paper, a new geometry index of very long baseline interferometry (VLBI) observing networks, the volume of network V, is examined as an indicator of the errors in the Earth orientation parameters (EOP) obtained from VLBI observations. It has been shown that both EOP precision and accuracy can be well described by the power law σ = aV ^c in a wide range of the network size from domestic to global VLBI networks. In other words, as the network volume grows, the EOP errors become smaller following a power law. This should be taken into account for a proper comparison of EOP estimates obtained from different VLBI networks. Thus, performing correct EOP comparison allows us to investigate accurately finer factors affecting the EOP errors. In particular, it was found that the dependence of the EOP precision and accuracy on the recording data rate can also be described by a power law. One important conclusion is that the EOP accuracy depends primarily on the network geometry and to lesser extent on other factors, such as recording mode and data rate and scheduling parameters, whereas these factors have a stronger impact on the EOP precision.     

On the consistency of the current conventional EOP series and the celestial and terrestrial reference frames

Precise transformation between the celestial reference frames (CRF) and terrestrial reference frames (TRF) is needed for many purposes in Earth and space sciences. According to the Global Geodetic Observing System (GGOS) recommendations, the accuracy of positions and stability of reference frames should reach 1 mm and 0.1 mm year\(^{-1}\), and thus, the Earth Orientation Parameters (EOP) should be estimated with similar accuracy. Different realizations of TRFs, based on the combination of solutions from four different space geodetic techniques, and CRFs, based on a single technique only (VLBI, Very Long Baseline Interferometry), might cause a slow degradation of the consistency among EOP, CRFs, and TRFs (e.g., because of differences in geometry, orientation and scale) and a misalignment of the current conventional EOP series, IERS 08 C04. We empirically assess the consistency among the conventional reference frames and EOP by analyzing the record of VLBI sessions since 1990 with varied settings to reflect the impact of changing frames or other processing strategies on the EOP estimates. Our tests show that the EOP estimates are insensitive to CRF changes, but sensitive to TRF variations and unmodeled geophysical signals at the GGOS level. The differences between the conventional IERS 08 C04 and other EOP series computed with distinct TRF settings exhibit biases and even non-negligible trends in the cases where no differential rotations should appear, e.g., a drift of about 20 \(\upmu \)as year\(^{-1 }\)in \(y_{\mathrm{pol }}\) when the VLBI-only frame VTRF2008 is used. Likewise, different strategies on station position modeling originate scatters larger than 150 \(\upmu \)as in the terrestrial pole coordinates.     

Satellite emission radio interferometric earth surveying series—GPS geodetic system

The satellites of the Global Positioning System (GPS) offer an important new geodetic resource making possible a highly accurate portable radio geodetic system. A concept called SERIES (Satellite Emission Radio Interferometric Earth Surveying) makes use of GPS radio transmissions without any satellite modifications. By employing the technique of very long baseline interferometry (VLBI) and its calibration methods, 0.5 to 3 cm three dimensional baseline accuracy can be achieved over distances of 2 to 200 km respectively, with only 2 hours of on-site data acquisition. The use of quasar referenced ARIES Mobile VLBI to establish a sparse fundamental control grid will provide a basis for making SERIES GPS measurements traceable to the time-invariant quasar directions. Using four SERIES stations deployed at previously established ARIES sites, allows the GPS satellite apparent positions to be determined. These apparent positions then serve as calibrations for other SERIES stations at unknown locations to determine their positions in a manner traceable to the quasars. Because this proposed radio interferometric configuration accomplishes its signal detection by cross-correlation, there is no dependence upon knowledge of the GPS transmitted waveforms which might be encrypted. Since GPS radio signal strengths are 10⁵ stronger than quasar signals, a great reduction in telecommunications sophistication is possible which will result in an order of magnitude less cost for a SERIES GPS station compared to a quasar based mobile VLBI system. The virtually all-weather capability of SERIES offers cost-effective geodetic monitoring with applications to crustal dynamics and earthquake research.     

The 2008 DGFI realization of the ITRS: DTRF2008

A new realization of the International Terrestrial System was computed at the ITRS Combination Centre at DGFI as a contribution to ITRF2008. The solution is labelled DTRF2008. In the same way as in the DGFI computation for ITRF2005 it is based on either normal equation systems or estimated parameters derived from VLBI, SLR, GPS and DORIS observations by weekly or session-wise processing. The parameter space of the ITRS realization comprises station positions and velocities and daily resolved Earth Orientation Parameters (EOP), whereby for the first time also nutation parameters are included. The advantage of starting from time series of input data is that the temporal behaviour of geophysical parameters can be investigated to decide whether the parameters can contribute to the datum realization of the ITRF. In the same way, a standardized analysis of station position time series can be performed to detect and remove discontinuities. The advantage of including EOP in the ITRS realization is twofold: (1) the combination of the coordinates of the terrestrial pole—estimated from all contributing techniques—links the technique networks in two components of the orientation, leading to an improvement of consistency of the Terrestrial Reference Frame (TRF) and (2) in their capacity as parameters common to all techniques, the terrestrial pole coordinates enhance the selection of local ties as they provide a measure for the consistency of the combined frame. The computation strategy of DGFI is based on the combination of normal equation systems while at the ITRS Combination Centre at IGN solutions are combined. The two independent ITRS realizations provide the possibility to assess the accuracy of ITRF by comparison of the two frames. The accuracy evaluation was done separately for the datum parameters (origin, orientation and scale) and the network geometry. The accuracy of the datum parameters, assessed from the comparison of DTRF2008 and ITRF2008, is between 2–5?mm and 0.1–0.8?mm/year depending on the technique. The network geometry (station positions and velocities) agrees within 3.2?mm and 1.0?mm/year. A comparison of DTRF2008 and ITRF2005 provides similar results for the datum parameters, but there are larger differences for the network geometry. The internal accuracy of DTRF2008—that means the level of conservation of datum information and network geometry within the combination—was derived from comparisons with the technique-only multi-year solutions. From this an internal accuracy of 0.32?mm for the VLBI up to 3.3?mm for the DORIS part of the network is found. The internal accuracy of velocities ranges from 0.05?mm/year for VLBI to 0.83?mm/year for DORIS. The internal consistency of DTRF2008 for orientation can be derived from the analysis of the terrestrial pole coordinates. It is estimated at 1.5–2.5?mm for the GPS, VLBI and SLR parts of the network. The consistency of these three and the DORIS network part is within 6.5?mm.     

VLBI在探月卫星定位中的应用分析

中国实施的"嫦娥"探月工程中,探月卫星的定轨测控系统由我国现有的S频段航天测控网(USB)和甚长基线干涉测量(VLBI)系统组成。系统中,VLBI技术主要为绕月卫星定轨提供卫星的角位置。本文分析了在探月项目中,VLBI单点定位的必要性。探讨了VLBI技术用于探月卫星单点定位的基本原理及其实现方法。通过算例对模拟数据进行处理,检验了方法的正确性。对结果进行分析,得出一些结论。     

IGS contribution to the ITRF

We examine the contribution of the International GNSS Service (IGS) to the International Terrestrial Reference Frame (ITRF) by evaluating the quality of the incorporated solutions as well as their major role in the ITRF formation. Starting with the ITRF2005, the ITRF 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. Analysis of time series of station positions is a fundamental first step in the ITRF elaboration, allowing to assess not only the stations behavior, but also the frame parameters and in particular the physical ones, namely the origin and the scale. As it will be seen, given the poor number and distribution of SLR and VLBI co-location sites, the IGS GPS network plays a major role by connecting these two techniques together, given their relevance for the definition of the origin and the scale of the ITRF. Time series analysis of the IGS weekly combined and other individual Analysis Center solutions indicates an internal precision (or repeatability) <2 mm in the horizontal component and <5 mm in the vertical component. Analysis of three AC weekly solutions shows generally poor agreement in origin and scale, with some indication of better agreement when the IGS started to use the absolute model of antenna phase center variations after the GPS week 1400 (November 2006).