应用iGMAS超快速星历的实时精密单点定位研究

中国主导建设的国际GNSS监测评估系统（iGMAS）相比国际上比较成熟的IGS系统在产品精度等方面存在差别，目前实时精密单点定位应用多采用IGS实时、近实时产品。为改变这一现状，针对iGMAS产品特性以及实时精密单点定位对超快速精密星历的需求，对iGMAS超快速星历的精度和稳定性方面进行评估，设计了iGMAS产品实时/事后下载应用程序，开展了基于iGMAS超快速星历的实时精密单点定位研究，并结合NovAtel OEM617双频接收机进行了GPS实时精密单点定位试验。实验结果表明，在连续观测23min后定位误差即可收敛到分米级，较接收机原始定位精度高一个量级，且稳定性好，最终在E/N/U方向定位误差均方根分别为7.2cm、6.4cm、15.2cm，与应用IGS超快速星历实时PPP试验取得相近的结果。研究实现了iGMAS数据获取、评估和实时PPP应用的一整套方案，验证了iGMAS超快速产品的性能，对推进iGMAS产品的应用提供了借鉴。     

GPS/GLONASS/GALILEO多星座组合导航系统研究

通过对多星座组合导航系统的分析,对多星座组合导航定位算法进行了研究。由于GPS、GLONASS、GALILEO系统分别采用了不同的坐标系,文中利用坐标系变换将不同星座统一到同一个坐标系中,采用增加状态变量的方法实现了组合系统的时间统一,从而实现了多星座间的时空统一;针对多星座组合系统的可见卫星数目大大增加的特点,采用最小二乘的方法实现了用户的导航定位解算,提高了用户的定位精度。通过对GPS/GLONASS/GALILEO多星座组合导航系统的仿真研究,其结果表明:与GPS单星座相比,GPS/GLONASS/GALILEO多星座组合导航系统,同一时间内可见卫星数目大大增加,PDOP值明显减小,有效提高了导航定位精度,具有良好的定位性能和可靠性。     

基于Hadamard方差的导航星座自主时间同步算法研究

文章在总结国内外研究成果基础上,针对GPS铷钟虽然短期稳定性较好,但采用Allan方差描述铷钟频率稳定性时,其钟差状态方程仅为两参数,在较长平滑时间里存在时钟漂移和甚低频噪声的影响,使噪声特性淹没或估值不收敛的缺陷,引入Hadamard方差建立了三参数系统状态误差模型,通过三次采样方差从模型上解决了线形漂移和甚低频噪声的影响问题。在时钟系统状态模型和星间双向测量方程建模基础上,给出了工程实用的标准Kalman基本滤波方程。数值分析仿真表明,采用Hadamard方差描述时钟频率稳定性显著提高星载时钟自主同步精度,从而克服了Allan方差描述产生的频率漂移影响较大和甚低频噪声不收敛的问题。     

基于GPS/北斗组合系统的高轨卫星定位技术研究

针对目前高轨GPS信号可用性差及定位精度低的特点, 对GPS/北斗组合系统的 高轨卫星定位技术进行研究, 对比分析了单GPS系统与GPS/北斗组合系统的卫 星可见性和几何精度因子. 结果表明, GPS/北斗组合系统比单GPS系统的卫星可 见性好, 且定位精度高. 同时通过提出在星载接收机上采用高精度原子钟, 可实现三星定位, 降低对接收机的技术要求.      

差分改正数的更新率对定位精度的影响分析

介绍了GPS(Global Positioning System)定位的误差来源及差分的具体形式,描述了伪距定位及伪距差分定位的数学模型.分析了差分修正后的用户伪距残差和SA(Selective Availability)信号的关系,并从SA的时间相关特性出发,分析了伪距差分中差分改正数的更新率对定位精度的影响.最后,利用RTCM SC 104格式的电文1进行了不同更新率的实验,得出了具有实践意义的结论.     

卫星高度角对GPS单向时间传递性能的影响研究

GPS单向时间传递可用于完成本地时间频率参考的校准及向协调世界时的高精度溯源。然而,卫星高度角的变化对时间频率传递的性能会产生一定的影响。概述了GPS卫星单向时间传递的基本原理,分析了不同卫星高度角对时间频率传递性能产生影响的原因,给出了时间频率传递性能的评估方法。最后,通过实验测试GPS卫星高度角变化对时间频率传递性能的影响,给出了实验结果并进行了分析。实验结果表明,随着卫星截止高度角的增加,GPS单向时间频率传递的性能将会提高,但是可视卫星的数量会减少。工程应用中,如果对时间频率传递长期稳定性要求较高,可将截至高度角设为25°;如只对短期稳定性要求较高,可将截止高度角设为5°。     

GPS软件接收机原理样机设计与实现

为了满足GPS/INS(Global Positioning System/Inertial Navigation System)超紧组合导航系统研究的需要,克服硬件接收机参数固定,适应性差的弱点,设计了一种参数可调、灵活控制的GPS软件接收机.采用GPS L1频率的中频采样信号,运用FFT(Fast Fourier Transform)频域捕获算法和锁相环与锁频环相互辅助的载波环路,实现了信号捕获、码环和载波环路跟踪、导航电文提取与解码、伪距及导航定位解算,并与NovAtel公司的FlexPax型硬件接收机进行了比较.跑车测试结果表明,该GPS软件接收机捕获迅速、跟踪准确,导航定位精度小于10m,动态抗干扰能力明显优于一般GPS硬件接收机,适合于GPS/INS超紧组合导航系统的应用.     

基于相位偏差的实时PPP模糊度固定

采用CNES发布的实时相位偏差数据，实现包含模糊度固定的实时精密单点定位.对全球10个IGS测站10天观测数据进行RTPPP解算，分别统计模糊度首次固定时间和定位精度，结果显示利用实时相位偏差数据能在平均30min内实现模糊度首次固定，模糊度固定时水平位置误差由6cm迅速降低至2cm左右，三维位置误差由10cm迅速降低至5cm左右，同时RTPPP模糊度固定在3h观测内可保持水平3cm、三维5cm左右的定位精度.通过分析得出，基于相位偏差的RTPPP模糊度固定技术具有较高的定位精度和定位稳定性，能够快速实现cm级定位.      

基于GEO／HEO混合星座的区域卫星定位系统性能分析

星座方案选择对卫星定位系统的性能具有很大影响。本文结合GNSS－2星座方案中的阿基米德计划，设计了针对我国区域的HEO／GEO混合星座，首先从可见性角度对GEO卫星定点参数和HEO卫星轨道参数进行了优选，之后，以GDOP和PDOP因子为衡量标准，在选定覆盖区域内，对混合星座的定位性能进行了分析，仿真结果表明，在覆盖区域内，该星座能够达到与GPS相当的定位精度。     

电离层不规则结构对GPS性能的影响

电离层不规则结构的存在可引起无线电信号的幅度和相位发生随机起伏, 这 种电离层闪烁现象会影响全球定位系统(Global Positioning System, GPS)的 性能, 降低定位精度, 严重时导致信号失锁. 电离层不规则结构对GPS性能的 影响涉及电离层物理、接收机设计和表征卫星几何分布的精度衰减因子(Dilution of Positioning, DOP)等多方面因素. 本文通过对表征电离层不规则结构参数 的分析, 根据GPS接收机跟踪环路和闪烁信号模型, 综合研究了电离层闪烁对 GPS接收机载波跟踪环和码跟踪环跟踪误差的影响; 结合实际观测, 评述 了电离层不规则结构对单频和双频GPS接收机定位性能的影响, 在此基础上 提出了有待深入研究的问题及具体建议.      

Performance assessment of real-time multi-GNSS integrated PPP with uncombined and ionospheric-free combined observables

For precise position services, the real-time precise point positioning (PPP) is a promising technology. The real-time PPP performance is expected to be improved by multi-system combination. The performance of real-time multi-system PPP needs to be periodically investigated, with the increasing number of available satellites and the continuously improved quality of real-time precise products of satellite clocks and orbits. In this study, a comprehensive performance assessment is conducted for the four-system integrated real-time PPP (FSIRT-PPP) with GPS, BDS, Galileo and GLONASS in both static and kinematic modes. The datasets from 118 stations spanning approximately a month are used for analysis, and the real-time stream CLK93 is employed. The superior performance of FSIRT-PPP is validated by comparing with the results of GPS/BDS, GPS/Galileo, GPS/GLONASS, GPS-only, BDS-only, Galileo-only and GLONASS-only cases. The FSIRT-PPP using ionospheric-free (IF) combined observables can achieve a convergence time of 10.9, 4.8 and 11.8 min and a positioning accuracy of 0.4, 0.5 and 0.7 cm in the static mode in the east, north and up directions, respectively, while the derived statistic is 15.4, 7.0 and 16.4 min, and 1.6, 1.2 and 3.4 cm in the kinematic mode in the three directions, respectively. Moreover, we also compare the position solutions of real-time PPP adopting IF combined and uncombined (UC) observables, and prove the mathematical equivalence between the two PPP models in the converged stage, provided that there are no external ionospheric corrections or constraints given to the estimated ionospheric delays in the UC model. The difference between the fully converged positioning accuracy of IF-based and UC-based real-time PPP is marginal, but the UC-based real-time PPP has longer convergence time due to the influence of the significant unmodeled time-varying errors in the real-time precise products as well as the different parameterization between them. For completeness, the real-time kinematic PPP results in harsh environments and the post-processed PPP results are also presented.     

The realization and convergence analysis of combined PPP based on raw observation

In order to speed up Precise Point Positioning (PPP)’s convergence, a combined PPP method with GPS and GLONASS which is based on using raw observations is proposed, and the positioning results and convergence time have been compared with that of single system. The ionospheric delays and receiver’s Differential Code Bias (DCB) corrections are estimated as unknown parameters in this method. The numerical results show that the combined PPP has not caused significant impacts on the final solutions, but it greatly improved Position Dilution of Precision (PDOP) and convergence speed and enhanced the reliability of the solution. Meanwhile, the convergence speed is greatly influenced by the receiver’s DCB, positioning results in horizontal which are better than 10 cm can be realized within 10 min. In addition, the ionosphere and DCB products can be provided with high precision.     

Comparison of convergence time and positioning accuracy among BDS,GPS and BDS/GPS precise point positioning with ambiguity resolution

Precise point positioning (PPP) usually takes about 30?min to obtain centimetre-level accuracy, which greatly limits its application. To address the drawbacks of convergence speed and positioning accuracy, we develop a PPP model with integrated GPS and BDS observations. Based on the method, stations with global coverage are selected to estimate the fractional cycle bias (FCB) of GPS and BDS. The short-term and long-term time series of wide-lane (WL) FCB, and the single day change of narrow-lane (NL) FCB are analysed. It is found that the range of GPS and BDS non-GEO (IGSO and MEO) WL FCB is stable at up to a 30-day-time frame. At times frame of up to 60?days, the stability is reduced a lot. Whether for short-term or long-term, the changes in the BDS GEO WL FCB are large. Moreover, BDS FCB sometimes undergoes a sudden jump. Besides, 17 and 10 stations were used respectively to investigate the convergence speed and positioning errors with six strategies: BDS ambiguity-float PPP (Bfloat), GPS ambiguity-float PPP (Gfloat), BDS/GPS ambiguity-float PPP (BGfloat), BDS ambiguity-fixed PPP (Bfix), GPS ambiguity-fixed (Gfix), and BDS/GPS ambiguity-fixed (BGfix). The average convergence speed of the ambiguity-fixed solution is greatly improved compared with the ambiguity-float solution. In terms of the average convergence time, the Bfloat is the longest and the BGfix is the shortest among these six strategies. Whether for ambiguity-float PPP or ambiguity-fixed PPP, the convergence reduction time in three directions for the combined system is the largest compared with the single BDS. The average RMS value of the Bfix in three directions (easting (E), northing (N), and up (U)) are 2.0?cm, 1.5?cm, and 5.9?cm respectively, while those of the Gfix are 0.8?cm, 0.5?cm, and 1.7?cm. Compared with single system, the BDS/GPS combined ambiguity-fixed system (BGfix) has the fastest convergence speed and the highest accuracy, with average RMS as 0.7?cm, 0.5?cm, and 1.9?cm for the E, N, U components, respectively.     

Performance assessment of RTPPP positioning with SSR corrections and PPP-AR positioning with FCB for multi-GNSS from MADOCA products

The state-space representation (SSR) product of satellite orbit and clock is one of the most essential corrections for real-time precise point positioning (RTPPP). When it comes to PPP ambiguity resolution (PPP-AR), the fractional cycle bias (FCB) matters. The Japan Aerospace Exploration Agency (JAXA) has developed a multi-GNSS (i.e., global navigation satellite system) advanced demonstration tool for orbit and clock analysis (MADOCA), providing free and precise orbit and clock products. Because of the shortage of relevant studies on performance evaluation, this paper focuses on the performance assessment of RTPPP and PPP-AR by real-time and offline MADOCA products. To begin with, the real-time MADOCA products are evaluated by comparing orbit and clock with JAXA final products, which gives an objective impression of the correction. Second, PPP tests in static and simulated kinematic mode are conducted to further verify the quality of real-time MADOCA products. Finally, the offline MADOCA products are assessed by PPP and PPP-AR comparisons. The results are as follows: (1) Orbit comparisons produced an average error of about 0.04–0.13 m for the global positioning system (GPS), 0.14–0.16 m for the global navigation satellite system (GLONASS), and 0.07–0.08 m for the quasi-zenith satellite system (QZSS). The G15 satellite had the most accurate orbit, with a difference of 0.04 m between the JAXA orbit products and MADOCA’s counterpart, while the R07 satellite had the least accurate orbit with a difference of 0.16 m. Clock products had an accuracy of 0.4–1.3 ns for GPS, 1.4–1.6 ns for GLONASS, and 0.7–0.8 ns for QZSS in general. The G15 satellite had the most accurate clock with a difference of only 0.40 ns between the JAXA clock products and MADOCA products, and the R07 satellite had the least accurate clock with a difference of 1.55 ns. The orbit and clock products for GLONASS performed worse than those of GPS and QZSS. (2) After convergence, the positioning accuracy was 3.0–8.1 cm for static PPP and 8.1–13.7 cm for kinematic PPP when using multi-GNSS observations and precise orbit and clock products. The PFRR station performed the good performance both in static and kinematic mode with an accuracy of 2.99 cm and 8.08 cm, respectively, whereas the CPNM station produced the worst static performance with an error of 8.09 cm, and the ANMG station produced the worst kinematic performance with a counterpart of 13.69 cm. (3) The PPP-AR solution was superior to the PPP solution, given that, with respect to PPP, post-processing PPP-AR improved the positioning accuracy and convergence time by 13–32 % (3–89 %) in GPS-only mode by 2–15 % (5–60 %) in GPS/QZSS mode. Thus, we conclude that the current MADOCA products can provide SSR corrections and FCB products with positioning accuracy at the decimeter or even centimeter level, which could meet the demands of the RTPPP and PPP-AR solutions.     

GLONASS-based precise point positioning and performance analysis

Current precise point positioning (PPP) techniques are mainly based on GPS which has been extensively investigated. With the increase of available GLONASS satellites during its revitalization, GLONASS observations were increasingly integrated into GPS-based PPP. Now that GLONASS has reached its full constellation, there will be a wide interest in PPP systems based on only GLONASS since it provides a PPP implementation independent of GPS. An investigation of GLONASS-based PPP will also help the development of GPS and GLONASS combined PPP techniques for improved precision and reliability. This paper presents an observation model for GLONASS-based PPP in which the GLONASS hardware delay biases are addressed. In view of frequently changed frequency channel number (FCN) for GLONASS satellites, an algorithm has been developed to compute the FCN for GLONASS satellites using code and phase observations, which avoids the need to provide the GLONASS frequency channel information during data processing. The observation residuals from GLONASS-based PPP are analyzed and compared to those from GPS-based PPP. The performance of GLONASS-based PPP is assessed using data from 15 globally distributed stations.     

Contribution analysis of QZSS to single-frequency PPP of GPS/BDS/GLONASS/Galileo

The Quasi-Zenith Satellite System (QZSS) established by the Japan Aerospace Exploration Agency mainly serves the Asia-Pacific region and its surrounding areas. Currently, four in-orbit satellites provide services. Most users of GNSS in the mass market use single-frequency (SF) receivers owing to the low cost. Therefore, it is meaningful to analyze and evaluate the contribution of the QZSS to SF precise point positioning (PPP) of GPS/BDS/GLONASS/Galileo systems with the emergence of GNSS and QZSS. This study compares the performances of three SF PPP models, namely the GRoup and PHase Ionospheric Correction (GRAPHIC) model, GRAPHIC with code observation model, and an ionosphere-constrained model, and evaluated the contribution of the QZSS to the SF PPP of GPS/BDS/GLONASS/Galileo systems. Moreover, the influence of code bias on the SF PPP of the BDS system is also analyzed. A two-week dataset (DOY 013–026, 2019) from 10 stations of the MGEX network is selected for validation, and the results show that: (1) For cut-off elevation angles of 15, 20, and 25°, the convergence times for the static SF PPP of GLONASS + QZSS are reduced by 4.3, 30.8, and 12.7%, respectively, and the positioning accuracy is similar compared with that of the GLONASS system. Compared with the BDS single system, the convergence times for the static SF PPP of BDS + QZSS under 15 and 25° are reduced by 37.6 and 39.2%, the horizontal positioning accuracies are improved by 18.6 and 14.1%, and the vertical components are improved by 13.9 and 21.4%, respectively. At cut-off elevation angles of 15, 20, and 25°, the positioning accuracy and precision of GPS/BDS/GLONASS/Galileo + QZSS is similar to that of GPS/BDS/GLONASS/Galileo. And the convergence times are reduced by 7.4 and 4.3% at cut-off elevation angles of 20 and 25°, respectively. In imitating dynamic PPP, the QZSS significantly improves the positioning accuracy of BDS and GLONASS. However, QZSS has little effect on the GPS-only, Galileo-only and GPS/BDS/GLONASS/Galileo. (2) The code bias of BDS IGSO and MEO cannot be ignored in SF PPP. In static SF PPP, taking the frequency band of B1I whose multipath combination is the largest among the frequency bands as an example, the vertical component has a systematic bias of approximately 0.4–1.0 m. After correcting the code bias, the positioning error in the vertical component is lower than 0.2 m, and the positioning accuracy in the horizontal component are improved accordingly. (3) The SF PPP model with ionosphere constraints has a better convergence speed, while the positioning accuracy of the three models is nearly equal. Therefore the GRAPHIC model can be used to get good positioning accuracy in the absence of external ionosphere products, but its convergence speed is slower.     

BDS-2/BDS-3 uncalibrated phase delay estimation considering the intra-system bias

Intra-system biases (ISBs) between BDS-2 and BDS-3 are of critical importance when combining observations from the BDS-2 and BDS-3 systems, which is meaningful to fully take advantage of the BDS positioning capability. Meanwhile, ISBs should also be considered in the estimation of BDS uncalibrated phase delays (UPDs). In this research, we present a BDS-2/BDS-3 joint-processing scheme, as well as a method for estimating BDS UPDs. The characteristics of ISBs and the quality of BDS UPDs are analyzed based on 30-day data from 130 multi-GNSS experimental (MGEX) stations. Our results indicate that the ISBs are related to the type and version of the receiver. The ISBs can be regarded as constant across the course of a given day, and the mean standard deviation (STD) values of ISBs over one month for different types of receivers are generally within 0.2 m. Moreover, to assess the quality of UPD products, the residuals of the estimated UPDs and the utilization rates of the observation data are computed. The results show that the quality of BDS UPDs can be improved by correcting the satellite-induced pseudo-range variations, and by estimating the wide-lane (WL) UPD difference between BDS-2 and BDS-3. The average RMS values of the estimated residuals of WL UPD and narrow-lane (NL) UPD are 0.07 and 0.09 cycles, respectively; moreover, the utilization rate of the observation data of WL UPD and NL UPD can reach above 90 %. The performance of BDS precise point positioning (PPP) and PPP ambiguity resolution (PPP-AR) is analyzed in terms of positioning accuracy and convergence performance in both the static and kinematic modes. Compared with PPP ambiguity-float solutions, the positioning accuracy of PPP-AR is significantly improved, especially in the east direction. The impact of ISBs on PPP and PPP-AR is also analyzed, and the results indicate that ISBs can improve the convergence speed of float PPP, but can be disregarded in PPP-AR.     

Comparative analysis of four different single-frequency PPP models on positioning performance and atmosphere delay retrieval

Single-frequency precise point positioning (SF-PPP) has attracted increasing attention due to its high precision and cost effectiveness. With various strategies to handle the dominant error, i.e., ionosphere delay, the ionosphere-float (IF), ionosphere-free-half (IFH), ionosphere-corrected (IC), and ionosphere-weighted (IW) SF-PPP models are certain to possess different characteristics and performance levels. This study is dedicated to assessing and comparing the four models from model characteristics, positioning performance, and atmosphere delay retrieval. The model comparison shows that IC and IW models are full-rank while IF and IFH models have a rank deficiency of size one that will result in biased estimations, which means the better solvability of IC and IW models. The experiments are carried out based on the 7-day Global Positioning System (GPS) observations collected at 57 global Multi-GNSS Experiment (MGEX) stations and Global Ionosphere Map (GIM) products. The results indicate that the IW model can accelerate SF-PPP convergence and achieve higher positioning accuracy compared to the other three SF-PPP models, especially in kinematic mode. With convergence criteria of 0.25 m in horizontal and 0.5 m in vertical, the east/north/up convergence times of IW model are 0.5/15.0/25.0 min and 0.5/16.0/36.5 min for static and kinematic modes, respectively. The IW model is able to achieve an instantaneous positioning accuracy of 0.28/0.35/0.75 m. In addition, a real kinematic test also demonstrates the best positioning solutions of IW model. Regarding troposphere delay retrieval, the IF, IFH, and IW models obtain a comparable daily accuracy of 3.0 cm on average, while the IC model achieves the worst accuracy of 8.0 cm. For precise ionosphere delay estimation, IW model only needs an average initialization time of 34.3 min, but a longer initialization time of 51.6 min is required for IF model. The daily precision of ionosphere delay estimation for IW model can reach up to 10.8 cm. At the present accuracy of GIM products, it is suggested that the IW model should be adopted for SF-PPP first due to its superior performance in positioning and atmosphere delay retrieval.     

An improved approach to model regional ionosphere and accelerate convergence for precise point positioning

Given the severe effects of the ionosphere on global navigation satellite system (GNSS) signals, single-frequency (SF) precise point positioning (PPP) users can only achieve decimeter-level positioning results. Ionosphere-free combinations can eliminate the majority of ionospheric delay, but increase observation noise and slow down dual-frequency (DF) PPP convergence. In this paper, we develop a regional ionosphere modeling and rapid convergence approach to improve SF PPP (SFPPP) accuracy and accelerate DF PPP (DFPPP) convergence speed. Instead of area model, ionospheric delay is modeled for each satellite to be used as a priori correction. With the ionospheric, wide-lane uncalibrated phase delay (UPD) and residuals satellite DCBs product, the wide-lane observations for DF users change to be high-precision pseudorange observations. The validation of a continuously operating reference station (CORS) network was analyzed. The experimental results confirm that the approach considerably improves the accuracy of SFPPP. For DF users, convergence time is substantially reduced.     

Towards PPP-RTK: Ambiguity resolution in real-time precise point positioning

Integer ambiguity resolution at a single station can be achieved by introducing predetermined uncalibrated phase delays (UPDs) into the float ambiguity estimates of precise point positioning (PPP). This integer resolution technique has the potential of leading to a PPP-RTK (real-time kinematic) model where PPP provides rapid convergence to a reliable centimeter-level positioning accuracy based on an RTK reference network. Nonetheless, implementing this model is technically subject to how rapidly we can fix wide-lane ambiguities, stabilize narrow-lane UPD estimates, and achieve the first ambiguity-fixed solution. To investigate these issues, we used 7 days of 1-Hz sampling GPS data at 91 stations across Europe. We find that at least 10 min of observations are required for most receiver types to reliably fix about 90% of wide-lane ambiguities corresponding to high elevations, and over 20 min to fix about 90% of those corresponding to low elevations. Moreover, several tens of minutes are usually required for a regional network before a narrow-lane UPD estimate stabilizes to an accuracy of far better than 0.1 cycles. Finally, for hourly data, ambiguity resolution can significantly improve the accuracy of epoch-wise position estimates from 13.7, 7.1 and 11.4 cm to 0.8, 0.9 and 2.5 cm for the East, North and Up components, respectively, but a few tens of minutes is required to achieve the first ambiguity-fixed solution. Therefore, from the timeliness aspect, our PPP-RTK model currently cannot satisfy the critical requirement of instantaneous precise positioning where ambiguity-fixed solutions have to be achieved within at most a few seconds. However, this model can still be potentially applied to some near-real-time remote sensing applications, such as the GPS meteorology.