Determination of station positions and velocities from laser ranging observations to Ajisai,Starlette and Stella satellites

The positions and velocities of the four Satellite Laser Ranging (SLR) stations: Yarragadee (7090), Greenbelt (7105), Graz (7839) and Herstmonceux (7840) from 5-year (2001–2005) SLR data of low orbiting satellites (LEO): Ajisai, Starlette and Stella were determined. The orbits of these satellites were computed from the data provided by 20 SLR stations. All orbital computations were performed by means of NASA Goddard’s GEODYN-II program. The geocentric coordinates were transformed to the topocentric North–South, East–West and Vertical components in reference to ITRF2005. The influence of the number of normal points per orbital arc and the empirical acceleration coefficients on the quality of station coordinates was studied. To get standard deviation of the coordinates determination lower than 1 cm, the number of the normal points per site had to be greater than 50. The computed positions and velocities were compared to those derived from LAGEOS-1/LAGEOS-2 data. Three parameters were used for this comparison: station coordinates stability, differences from ITRF2005 positions and velocities. The stability of coordinates of LEO satellites is significantly worse (17.8 mm) than those of LAGEOS (7.6 mm), the better results are for Ajisai (15.4 mm) than for Starlette/Stella (20.4 mm). The difference in positions between the computed values and ITRF2005 were little bit worse for Starlette/Stella (6.6 mm) than for LAGEOS (4.6 mm), the results for Ajisai were five times worse (29.7 mm) probably due to center of mass correction of this satellite. The station velocities with some exceptions were on the same level (≈1 mm/year) for all satellites. The results presented in this work show that results from Starlette/Stella are better than those from Ajisai for station coordinates determination. We can applied the data from LEO satellites, especially Starlette and Stella for determination of the SLR station coordinates but with two times lower accuracy than when using LAGEOS data.     

SPOT-2 and TOPEX/Poseidon precise orbit determination from DORIS doppler tracking

The French earth observation satellite SPOT-2 has served as a testbed for precise orbit determination from DORIS doppler tracking in anticipation of the TOPEX/Poseidon mission. Using the most up-to-data gravity field model, JGM-2, a radial orbit accuracy of about 2–9 cm was achieved, with an rms of fit of the tracking data of about 0.64 mm/s. Furthermore, it was found that the coordinates of the ground stations can be determined with an accuracy of the order of 2–5 cm after removal of common rotations, and translations.

Using a slightly different model for atmospheric drag, but the same gravity model, precise orbits of TOPEX/Poseidon from DORIS tracking data were determined with a radial orbit accuracy of the order of 4–5 cm, which is far within the 13 cm mission requirement. This conclusion is based on the analysis of 1-day overlap of successive 11-day orbits, and the comparisons with orbits computed from satellite laser tracking (SLR) and from the combination of SLR and DORIS tracking. Results indicate a consistency between the different orbits of 1–4 cm, 4–20 cm, and 6–13 cm in the radial, cross-track, and along-track directions, respectively. The residual rms is about 4–5 cm for SLR data and 0.56 mm/s for DORIS tracking. These numbers are roughly twice as large as the system noise levels, reflecting the fact that there are still some modeling errors left.     

DPOD2005: An extension of ITRF2005 for Precise Orbit Determination

For Precise Orbit Determination of altimetry missions, we have computed a data set of DORIS station coordinates defined for specific time intervals called DPOD2005. This terrestrial reference set is an extension of ITRF2005. However, it includes all new DORIS stations and is more reliable, as we disregard stations with large velocity formal errors as they could contaminate POD computations in the near future. About 1/4 of the station coordinates need to be defined as they do not appear in the original ITRF2005 realization. These results were verified with available DORIS and GPS results, as the integrity of DPOD2005 is almost as critical as its accuracy. Besides station coordinates and velocities, we also provide additional information such as periods for which DORIS data should be disregarded for specific DORIS stations, and epochs of coordinate and velocity discontinuities (related to either geophysical events, equipment problem or human intervention). The DPOD model was tested for orbit determination for TOPEX/Poseidon (T/P), Jason-1 and Jason-2. Test results show DPOD2005 offers improvement over the original ITRF2005, improvement that rapidly and significantly increases after 2005. Improvement is also significant for the early T/P cycles indicating improved station velocities in the DPOD2005 model and a more complete station set. Following 2005 the radial accuracy and centering of the ITRF2005-original orbits rapidly degrades due to station loss.     

A VLBI baseline post-adjustment approach for station velocity estimation in Eurasian continent

Baseline lengths and their time-derivatives among 58 geodetic VLBI stations were fitted by using 4439 observing sessions from the International VLBI Service for Geodesy and Astrometry (IVS). First, the velocities of eight stations in Eurasian continent were set as unknown quantities. Then, two standard global solutions from 3523 IVS sessions and 1110 sessions from database code XA, respectively, were applied prior to all-station coordinates and the non-estimated station velocities. Finally, from the relations among the coordinates, velocities, baseline length and its time-derivative, two types of baseline post-adjustment (BPA) were used to estimate the velocities of the eight stations. We discuss the data processing details, including the effect of different prior values for the stations and the optimal solution.     

Studies of the geocenter motion using 16-years DORIS data

An accuracy of geocenter motion estimation is strongly dependent on the geodetic network size and stations distribution over the Earth’s surface. From this point of view DORIS system has an advantage, as its ground network of beacons consists of more than 50 sites, equally distributed over the Earth’s surface. Aiming to study variations of the geocenter movements, the results of DORIS data analysis for the time span 1993.0–2009.0 (inawd06.snx series), performed at the Analysis Centre of the Institute of astronomy of the Russian Academy of Sciences, have been used. DORIS data processing was made with GIPSY/OASIS II software, developed by Jet Propulsion Laboratory and modified for DORIS data processing by Institute Géographique National. Standard deviations of stations coordinates are estimated at the level 0.5–4.0 cm (internal consistency), depending on the number of satellites used in the solution. RMS of estimated components of the DORIS satellites orbits, compared with the solutions of other IDS analysis centres, do not exceed 1–2 cm. Weekly solutions for coordinates have been transformed from free network solutions (inawd06.snx series) to a well defined terrestrial reference frame ITRF2005 with the use of seven parameters of Helmert transformation, which were examined with a view to study variations of the geocenter movements (ina05wd01.geoc time series). In order to estimate linear trend, amplitudes, periods and phases of geocenter variation a method of linear regression was applied. The evaluated amplitudes of annual variations are of the order of 5–7 mm for X and Y components and 27–29 mm for Z component. Semi-annual amplitudes are also noticeable in all components (1–34 mm for X, Y and Z components). Secular trends in the DORIS geocenter coordinates are: −1.2, −0.1 and −0.3 mm/year for X, Y and Z directions respectively.     

基于N-BEATS的单站对流层天顶总延迟预报

高精度对流层延迟先验值有助于加速精密单点定位的快速收敛。基于高精度高分辨率气象数据库,采用深度学习N BEATS算法,进行了单站对流层天顶总延迟的预报试验。试验选取了9个IGS跟踪站,试验弧段从2002年1月至2019年6月共185a。首先基于N BEATS算法,设计了3种预报策略,然后基于前175a针对不同预报策略进行模型训练,并对最后365d的对流层天顶总延迟进行预报。试验结果表明,以该气象数据库为基准,12h以内预报弧段的预报残差均值量级大多可达亚毫米,2h、4h、6h的预报残差的标准差分别约为5mm、9mm、13mm。     

Azimuthal propagation of Pc5 geomagnetic field pulsations in the southern polar cap

A statistical analysis of low frequency geomagnetic fluctuations at the two Antarctic stations Mario Zucchelli Station (geographic coordinates: 74.7°S, 164.1°E; corrected geomagnetic coordinates: 80.0°S, 306.8°E) and Dumont D’Urville (geographic coordinates: 66.7°S, 140.0°E; corrected geomagnetic coordinates: 80.4°S, 236.0°E) is shown. The analysis focuses on power spectra, coherence and phase difference between the stations, which are both located in the polar cap, with a 5-h magnetic local time displacement along a geomagnetic parallel; in this situation, the phase difference between geomagnetic fluctuations indicates the direction of their azimuthal propagation. Coherent fluctuations have been found to occur preferably when both stations are on the same side (dawnward or duskward) with respect to the polar cusp; moreover, around local magnetic midnight, they occur essentially during open magnetospheric conditions. The phase difference for coherent fluctuations indicates a propagation direction away from local geomagnetic noon and midnight. Also the analysis of three individual pulsation events, occurring at different times during the day, is shown; they are characterized at the two stations by simultaneous, coherent fluctuations, whose phase difference finds correspondence with the statistical behaviour.     

Impact of loading effects on determination of the International Terrestrial Reference Frame

The International Terrestrial Reference Frame (ITRF), as a realization of the International Terrestrial Reference System (ITRS), is represented by a set of station positions and linear velocities. They are intended to be used as regularized coordinates to which some corrections should be added to access instantaneous coordinates. The latest ITRS realization is the ITRF2005, which has integrated time series of station positions to form long-term solutions for the four space geodetic techniques. Currently, a purely linear model is used to parameterize station displacements in the estimation process, plus occasional discontinuities in case of earthquakes or equipment changes. However the input data have been derived without applying surface loading models and so surface loading effects are supposed to be embedded in the coordinates as measured quantities. We evaluate the effect of applying a posteriori loading corrections, which include the effect of atmospheric, non-tidal ocean, and continental water loading, to time series of positions estimated from Satellite Laser Ranging (SLR), Very Long Baseline Interferometry (VLBI), and Global Positioning System (GPS) data. We notice that they reduce about 50% or more of the annual signals in the translation and scale parameter time series of the SLR and VLBI techniques, except in SLR Z translation. In general, the estimated secular frame definition is negligibly affected and estimated positions and velocities are not significantly modified for stations that have accumulated a large number of observations. A multi-technique combination of such derived frames allows concluding that, for some cases, loading model corrections might degrade co-located station coordinates almost as much as they benefit them. However, most significant improvement of the estimated secular coordinates is observed for stations with less than 100 estimated positions as demonstrated with a multi-technique combination.     

一种应用于实时深空干涉测量的电离层时延修正方法

针对"嫦娥4号"中继星任务S频段信标信号的高精度实时干涉测量需求,结合深空测控干涉测量系统采用的稀疏标校工作模式,研究验证了一种面向测控模式实时干涉测量的电离层时延修正方法。首先分析了电磁波经电离层传播的延迟机理及特性;基于深空站历史观测数据,通过自相关函数分析验证了天顶向TEC的周日特性;在此基础上,结合深空干涉测量中心数据处理设备软件系统,讨论了电离层时延修正方法;通过任务期间的实测数据处理分析,验证了所提方法可以将实时测量精度提升1～3 ns,对低仰角跟踪弧段,该技术方法优势更为明显。该方法为后续深入推进深空测控干涉测量系统在任务中的实时应用提供了技术储备。     

利用"火星快车"三程多普勒跟踪数据限定局部洛伦兹不变性

目前航天器的三程多普勒跟踪技术已经在深空探测的控制与导航领域起到了重要作用。利用包含了对局部洛伦兹不变性(LLI)以及局部位置不变性(LPI)原理有破坏的三程多普勒跟踪理论,研究分析了"火星快车"(MEX)三程多普勒跟踪数据的残差。这些多普勒观测于2009年8月7日和8日进行,利用了欧洲航天局(ESA)在澳大利亚新诺舍(New Norcia)的上行站和三个分别在中国上海、昆明以及乌鲁木齐的下行站。我们发现,这些观测结果给出的LLI上限在10-2的量级。但由于各观测站本身对频率测量的精度有限,这些数据并不适合于检验LPI。     

Comparison of GPS S/C orbits determined from GPS and SLR tracking data

The last two Global Positioning System (GPS) spacecraft launched on August 30, 1993 and March 10, 1994 carried identical laser retroreflector arrays that allow laser ranging from ground stations. Tracking of these targets is a low priority for most of the ground stations due to the present mandatory support of missions with no alternative tracking. The available data is sparse and concentrated primarily on the first of the two spacecraft, GPS-35 (PRN 5). Despite this, analysis of this data set indicates that there is great potential for engineering and scientific experiments in a synergistic application of Satellite Laser Ranging (SLR) and radiometric GPS techniques.     

基于球场重建的球员运动数据分析

足球比赛中球员运动数据分析对增加观众的观看体验和辅助教练进行球员评估有着重要意义。球员运动数据分析的难点在于如何定位球员在球场上的坐标,即如何确定足球视频中单帧画面出现的缺损球场与标准二维球场之间的映射关系。针对如何在足球比赛中克服相机的高速移动和视角剧烈变化,设计并提出了利用球场重建与球员跟踪来进行球员运动数据分析的方法。球场重建方面,将足球视频中的球场分组为左中右3部分,每组通过球场分割、球场直线检测、球场直线分组、球场中圈点集合识别和球场关键点匹配来实现缺损球场到标准球场的映射;球员跟踪采用核相关滤波(KCF)跟踪算法,得到了球员运动数据统计的可视化结果。结合球场重建和球员跟踪算法定位球员的标准坐标,统计球员的一系列运动数据并进行可视化分析。提出的球员运动数据分析方法能够准确而快速地统计出球员的运动数据,包括球员坐标、运动轨迹、奔跑速度、活动范围和球员间距。球场重建方面采用图像交并进行评估,交并比达到87%,相比于传统的基于字典查询的方法(交并比为83.3%)准确度提升了3.7%。实验结果表明：所提出的球场重建方法能够更准确地表示球场映射关系,为球员运动数据分析统计提供更好的支持...     

Re-analysis of the cosmic ray ground level enhancement of 4 May 1960.

The relativistic solar particle event of 4 May 1960, resulting in a cosmic ray ground level enhancement, occurred well before modern analysis techniques were available. We have located surviving data from 23 neutron monitors and have used these to estimate the spectrum, mean arrival direction and particle pitch angle distribution as the event progressed. We find that the apparent particle arrival direction was at equatorial latitudes, over northern South America, in contrast to contemporary analyses that proposed it to be over North America. Our modified power law spectra are broadly consistent with earlier results. Data from stations above sea level need to be corrected for altitude using a two-attenuation length technique. The standard method involves comparison of data from two relatively close stations at significantly different altitude. We have shown that this method may be unreliable in cases, such as this, of quite sharp anisotropy.     

星载GPS低轨卫星几何法精密定轨研究

本文讨论了星载ＧＰＳ接收机为单频情形的代轨卫星几何法定轨,包括载波一相对定轨法和动态网定轨法,并利用Ｔｏｐｅｘ／Ｐｏｓｅｉｄｏｎ卫星星载ＧＰＳ实测数据中Ｌ１载波相位观测值进行验证,结果表明,载波相位相对定轨精度与地面基准站的观测质量有关,其三轨道位置精度为分米级;载波相位动态风定轨精度介于各基准站皮相位相对定轨之间,它相当于在各基准站相对定轨之间加权均衡,而且提高了定轨的可靠性。     

Spectral and temporal characteristics of aerosol optical depth over a wet tropical location in North East India

A network of multi wavelength solar radiometer (MWR) stations has been in operation since the 1980s in India for measurement of aerosol optical depth (AOD). This network was augmented recently with the addition of a large number of stations located across the length and breath of India covering a variety of climate regimes. The spectral and temporal variations of aerosol optical depths observed over Dibrugarh located in the North East of India (27.3°N, 94.5°E) are investigated by analyzing the data obtained from a MWR during October 2001–September 2003 using the Langley technique. AOD varies with time of the day, month of the year and season. From January to April and October to December, aerosol optical depth decreases with wavelength whereas during May–September aerosol optical depth has been found to be nearly independent of wavelengths. AOD is higher during pre-monsoon season (March–May) and lower in the monsoon (June–September) season at about all wavelengths. The temporal variation of AOD over Dibrugarh have also been compared with those reported from selected locations in India.     

双站跟踪模式下“嫦娥4号”中继星定轨仿真分析

位于地月平动点的探测器因为较差的观测几何,需要地基USB/UXB与天文VLBI长时间的联合跟踪数据获取稳定精确的轨道。提出了利用中国深空网双站共视跟踪平动点探测器,获取双程、三程测距及VLBI测量数据,解算探测器精确轨道的模式。以"鹊桥"卫星为分析对象,首先评估中国深空网对"鹊桥"的跟踪能力。然后分析不同观测组合模式下的定轨计算精度。结果表明:双站共视约束下,深空站每天对"鹊桥"跟踪弧长大于5 h;使用长于6 h的双站跟踪数据进行定轨,系统差的解算更有利于轨道精度提升;跟踪时长超过2天时,必须在轨道解算的同时估计光压系数,并有望实现优于百米的轨道精度。     

Polar motion and earth rotation series from LAGEOS

As an Associate Analysis Centre for the MERIT campaign ESOC has processed the complete data set of LAGEOS laser ranging normal points for the 14 months period of the campaign (September 1983 to October 1984), using independently developed software. Geocentric coordinates have been obtained for a network of 20 laser stations, and series for the pole, UT1R and length of day have been derived with 5 day resolution. Discrepancies in the node and inclination estimates between successive monthly arcs are attributed to mismodelled tidal effects. Detailed comparisons have been made with solutions obtained at other processing centres based on VLBI and SLR techniques.     

一种非线性航迹自适应跟踪控制方法

针对小型无人飞行器航迹跟踪精度和飞行品质问题,提出了一种基于非线性航迹的自适应跟踪控制方法.应用五阶B样条拟合航点,构建非线性期望航迹;建立基于非线性期望航迹Serret-Frenet坐标系下的位置和运动航向误差方程;根据误差方程设计渐近稳定收敛的自适应运动航向控制律.并应用此方法进行了外场飞行实验,实验结果表明自适应航迹跟踪控制方法有效且能保证航迹跟踪精度.     

一种新的四路激光跟踪坐标测量系统的基点标定方法

在激光跟踪三维坐标测量系统中,系统首先要经过标定,才可以进行实际测量。因此系统标定的准确度会直接影响系统最终测量准确度。为了提高标定准确度,提出一种新的标定方法,即利用直线法原理,在优化布局下,采用两路激光干涉仪测量出基点的三维坐标。仿真结果表明目标点相对于基点呈对称分布时,测得的基点坐标准确度比非对称分布时高6到7倍,并通过实验得出了在该优化布局下的基点坐标。     

Near real-time PPP-based monitoring of the ionosphere using dual-frequency GPS/BDS/Galileo data

Ionosphere delay is very important to GNSS observations, since it is one of the main error sources which have to be mitigated even eliminated in order to determine reliable and precise positions. The ionosphere is a dispersive medium to radio signal, so the value of the group delay or phase advance of GNSS radio signal depends on the signal frequency. Ground-based GNSS stations have been used for ionosphere monitoring and modeling for a long time. In this paper we will introduce a novel approach suitable for single-receiver operation based on the precise point positioning (PPP) technique. One of the main characteristic is that only carrier-phase observations are used to avoid particular effects of pseudorange observations. The technique consists of introducing ionosphere ambiguity parameters obtained from PPP filter into the geometry-free combination of observations to estimate ionospheric delays. Observational data from stations that are capable of tracking the GPS/BDS/GALILEO from the International GNSS Service (IGS) Multi-GNSS Experiments (MGEX) network are processed. For the purpose of performance validation, ionospheric delays series derived from the novel approach are compared with the global ionospheric map (GIM) from Ionospheric Associate Analysis Centers (IAACs). The results are encouraging and offer potential solutions to the near real-time ionosphere monitoring.