Hourly Updated Precise Orbit Products of Quad-constellation

Transcription

1 Hourly Updated Precise Orbit Products of Quad-constellation Satellites in IGS Analysis Center at Wuhan University Qile Zhao; Hongyang Ma; Xiaolong Xu; Jing Guo; Min Li Wuhan University Jul Paris

2 CONTENT Current Status of WHU ultra-rapid GNSS Products The Update of PANDA Software Hourly-updated Orbit Outlook

3 Current Status of WHU ultra-rapid GNSS Products : 6 hours update ultra-rapid GPS products. (Provided to IGS) : 1 hour update ultrarapid Quad-constellation products. (Published in IGS Data Center at Wuhan University) : 6 hours update ultra-rapid Quad-constellation products. (Provided to igmas) : 3 hours update ultra-rapid Quad-constellation products. (Published in IGS Data Center at Wuhan University)

4 1. Current Status of WHU ultra-rapid GNSS Products 1 hour update products Products Hour whuwwwwd_hh.sp3 whuwwwwd_hh.clk whuwwwwd_hh.erp GPS week. Hour index, from 00 to 23 GPS day. Publish n n 1 30 min Publish 30 min Publish 30 min Publish 30 min Publish 30 min Publish

5 2. The Update of PANDA Software Least-squares Orbit Integration Post Batch Process Initial State Epoch Initial Time Error Transformation Matrix Observation Orbit IntegrationOrbit Integration Error

6 2. The Update of PANDA Software Post Batch Process Epoch Advantage Disadvantage Observation Determination Prediction High precision Consuming Publish Bad real-time performance. The precision of predicted orbits decreases rapidly over time. Orbit Integration

7 2. The Update of PANDA Software Post least-squares batch Least-squares collocation method Update Fast least-squares batch & Real-time filtering on each epoch Strong constrain between epoch or observation length. Need the precise dynamic models for the GNSS satellite Variable constrain between epoch or observation length. When enough real-time GNSS data are avaliable

8 2. The Update of PANDA Software Fast least-squares batch Observation length Normal Matrix Advantages Strong constrain Consuming Observation Part: Approach or reach the precision of post batch process. Prediction Part: Quick update, less consuming, high precision of user available part.

9 2. The Update of PANDA Software Real-time filtering N11 N12 X1 W1 N 21 N 22 Y W 2 Observation length or epoch Reduce the weight of previous observation. Consuming Constant Parameters Equivalence X1 X 2 0 P Other Parameters N11 N12 X 2 W1 N 21 N 22 Y W 2

10 2. The Update of PANDA Software Initial state transformation Fast least-squares batch Initial State Epoch wise Each batch Transfer Real-time filtering on each epoch Transfer

11 The First Arc e.g. 33 hours The Second Arc e.g. 3 hours Publish User available part Processing 2.5h Delay 30(20:data +10: processing) min Strategy Hour

12 Quick Update Observation Prediction Advantage WHU hourly orbit WHU 3 hours orbit igmas ultra-rapid orbit 30 min Consuming 2h Consuming 2h Consuming 0.5h~1.5h User Available Part 2h~5h User Available Part 2h~8h User Available Part IGU ultra-rapid orbit 3h Consuming 12 3h~9h User Available Part

13 High Precision User available part Advantage The First Arc The Second Arc Observation Consuming Prediction Observation Part: Approach or reach the precision of post batch process. Prediction Part: Quick update, less consuming, high precision of user available part. 13

14 Models Item Model System GPS GLONASS BDS GALILEO Observation divided 33 hours + 3 hours Observation Un-difference Ionosphere-free combination Frequency L1&L2 G1&G2 B1&B2 E1&E5a Sampling rate of data 30s Sampling rate of parameter estimation 300s Elevation GPS/GLS:7 degree BDS/GAL:10 degree Orbit parameter 3 initial position 3 initial velocity 5 BERN solar pressure Integration step 60s Perturbation Non-spherical gravity, three body gravity, solid earth tide, Earth radiation, Solar radiation Satellite phase center GPS/GLS:igs08.atx BDS/GAL:WHU 14

15 Global Station WHU: IGS Data Center igmas Date Center IGS data MGEX data igmas data BETS data CMONOC data 15

16 Test Date: DOY 106~112, 2017 Orbit Precision Reference Orbit: GPS IGS Final Orbit GLONASS IGL Final Orbit BDS GBM Final Orbit GALILEO GBM Final Orbit Arc Comparison: Observation part 0.5h~1.5h in prediction part 2h~5h in prediction part 3h~9h in prediction part 16

17 GPS GPS 1D RMS (mm) DOY 0.5~1.5h 2~5h 3~9h AVE STD AVE STD AVE STD AVERAGE

18 IGU GPS GPS 1D RMS (mm) DOY 0.5~1.5h IGU(3~9h) AVE STD AVE STD AVERAGE

19 GLONASS GLONASS 1D RMS (mm) DOY 0.5~1.5h 2~5h 3~9h AVE STD AVE STD AVE STD AVERAGE

20 BDS GEO BDS GEO 1D RMS (mm) DOY 0.5~1.5h 2~5h 3~9h AVE STD AVE STD AVE STD AVERAGE

21 BDS I/M BDS I/M 1D RMS (mm) DOY 0.5~1.5h 2~5h 3~9h AVE STD AVE STD AVE STD AVERAGE

22 GALILEO GALILEO 1D RMS (mm) DOY 0.5~1.5h 2~5h 3~9h AVE STD AVE STD AVE STD AVERAGE

23 ftp://igs.gnsswhu.cn/pub/whu/mgex/ Brief Summary 23

24 4. Outlook More real-time data, and Data quality control New force models for each constellation, eg. GLONASS, BDS GEO Parallel computing 6 hours 1 hour Orbit renewal frequency 3 hour 15 min Delay 3 hours 24 2 hours 30 min 5 min

25 Thanks for your attention!