高时间分解能GPSGNSS.ppt

高時間分解能GPS/GNSS 衛星時計推定 PPP応用 Estimation of high-rate GPS/GNSS Satellite clock and application to the kinematic-PPP 技術 技術 高須 知二高須 知二 Tomoji TAKASUTomoji TAKASU 地球惑星科学連合2006年大会 http:/gpspp.sakura.ne.jp 背景 PPP(精密単独測位)精度 (24H) :水平4mm,垂直10mm :水平1-2cm,垂直3-4cm 精密暦高品質化 IGS Final : 実効5mm程度精度 補正高度化 GPS/GNSS衛星時計補間誤差問題 PPP制約 http:/gpspp.sakura.ne.jp PPP (300s) MIZU(IGS) : 2004/12/25 0:00:00-2:55:00GPST IGSIGR CODEJPL 10cm E-W N-S U-D http:/gpspp.sakura.ne.jp PPP (30s) MIZU(IGS) : 2004/12/25 0:00:00-2:59:30GPST IGS (Interpolated)IGR (Interpolated) CODEJPL 10cm E-W N-S U-D http:/gpspp.sakura.ne.jp PPP (1Hz) MIZU(IGS) : 2004/12/25 0:00:00-2:59:59GPST IGS (Interpolated)IGR (Interpolated) CODE (Interpolated)JPL (Interpolated) 10cm E-W N-S U-D http:/gpspp.sakura.ne.jp HR-衛星時計推定(1) IGS Sat Clock (300s) Clock Estimation IGS Sat Orbit (15min) QC (HR) (300s) HR Sat Clock HR Rcv Clock HR Sat Clock Interpolation Static-PPP(300s) ZTDPOSBIAS Rcv Clock Step1 Step2 OBS http:/gpspp.sakura.ne.jp HR-衛星時計推定(2) IGS Sat Clock(300s) IGS Sat Orbit(15min) Rcv Clock ZTDBIAS OBS (300s) HR Rcv Clock HR Sat Clock Step1 Step2 OBS (HR) POS Interpolation http:/gpspp.sakura.ne.jp HR-衛星時計推定(3) EKF/S (PPP:3pass, CLK:1pass) 電離層遅延:電離層線形結合 対流圏遅延: ZTD推定+NMF/GMF : IGS01.PCV/IGS05.ATX 局位置変動: 固体地球潮汐,海洋荷重,極運動潮汐 Phase Windup補正 http:/gpspp.sakura.ne.jp GT0.6.2 http:/gpspp.sakura.ne.jp HR-IGS局 (2004/12, 58 Stas) http:/gpspp.sakura.ne.jp HR-衛星時計(1) PRN01 : 2004/12/25 0:00:00-23:59:59 GPST wrt IGS Final (300s) http:/gpspp.sakura.ne.jp HR-衛星時計(2) wrt IGS Final Interpolated (1Hz) PRN01 : 2004/12/25 0:00:00-23:59:59 GPST http:/gpspp.sakura.ne.jp 衛星時計安定度 2004/12/25 0:00:00-23:59:59 GPST Block II/IIA Cs Block II/IIA Rb Block IIR Rb Ref: AMC2 (HM) 150s15s Clock Interpolation Error 300s : 0.05-0.3ns 30s : 0.01-0.03ns http:/gpspp.sakura.ne.jp HR-PPP HR-衛星時計使高時間分解能 PPP 高精度化/高安定化 位相 対流圏遅延推定(ZTD/Gradient) Mapping Function 食衛星取扱(Block II/IIA) Sidereal Filter http:/gpspp.sakura.ne.jp 1Hz HR-PPP測位解例 MIZU(IGS) : 2004/12/25 0:00:00-2:59:59GPST CODE (Interpolated) 1Hz HR-PPP E-W N-S U-D 10cm http:/gpspp.sakura.ne.jp Sidereal Filter MIZU(IGS) : 2004/12/25 0:00:00-2:59:59GPST without Sidereal Filter with Sidereal Filter E-W N-S U-D 10cm http:/gpspp.sakura.ne.jp 相対測位比較 http:/gpspp.sakura.ne.jp 解析実行時間例 Pentium4 3.2GHz, Win XP, GT0.6.2, Matlab6.5.1 http:/gpspp.sakura.ne.jp HR-PPP応用例 新潟県中越地震 地殻変動計測 (2004/10/23, M6.8) 地震 地震波観測 (2004/12/26, M9.0) http:/gpspp.sakura.ne.jp 新潟県中越地震 地殻変動(1) 新潟県中越地震 M6.8, 2004/10/23 8:56 UTC http:/gpspp.sakura.ne.jp 新潟県中越地震 地殻変動(2) 電子基準点 位置変動 E-W : 2004/10/23 GPST 10cm http:/gpspp.sakura.ne.jp 新潟県中越地震 地殻変動(3) 電子基準点 位置変動 N-S : 2004/10/23 GPST 10cm http:/gpspp.sakura.ne.jp 新潟県中越地震 地殻変動(4) 電子基準点 位置変動 U-D : 2004/10/23 GPST 10cm http:/gpspp.sakura.ne.jp 新潟県中越地震 地殻変動(3) 電子基準点 小千谷(950240) : 2004/10/23 GPST 8:00 9:0010:0011:00 -50 -40 -30 -20 -10 0 10 20 30 M6.8 M6.3 M6.0 M6.5M5.7 E-W N-S U-D Displacements (cm) 10cm http:/gpspp.sakura.ne.jp 新潟県中越地震 地殻変動(3) -10 0 10-10 0 10 -10 0 10 (cm) E N S W 電子基準点 小千谷(950240) 2004/10/23 GPST 8:00:008:56:00 8:56:30 8:57:00 9:00:00 9:12:00 9:12:0 9:34:00 :34:30 :35:00 11:00 M6.8 8:56 M6.3 9:03 M6.0 9:11 M6.5 9:34 M5.7 10:45 http:/gpspp.sakura.ne.jp 8:00 9:0010:0011:00 -50 -40 -30 -20 -10 0 10 20 30 Displacements (cm) 新潟県中越地震 地殻変動(4) 電子基準点 小千谷(950240) : 2004/10/23 GPST E-W N-S U-D: GSI F2解 10/22 : GSI F2解 10/2410cm http:/gpspp.sakura.ne.jp 地震 地震波 (1) Sumatra - Andaman Islands Earthquake M9.0, 2004/12/26 0:58:53UTC http:/gpspp.sakura.ne.jp 地震 地震波観測 (2) (2004/12/26 1:00-2:00 GPST) JOGJ 2024km IISC 2283km DGAR 2861km WUH2 3586km POL2 4841km USUD 5672km ARTU 6754km NRIL 7351km CHUR 13076km 10cm E-W N-S U-D http:/gpspp.sakura.ne.jp 地震 地震波観測 (3) POL2(IGS) : 2004/12/26 GPST R:3.4km/s=4841km/1440s 0:59 24min R L http:/gpspp.sakura.ne.jp 地震 地震波観測 (4) POL2(IGS) - AAK(GSN-IRIS) (d=17.1km) http:/gpspp.sakura.ne.jp 地震 地震波観測 (5) ARTU(IGS) - ARU(GSN-IRIS) (d=0.3km) http:/gpspp.sakura.ne.jp 地震 地震波観測 (6) DGAR(IGS) - DGAR(GSN-IRIS) (d=18.2km) http:/gpspp.sakura.ne.jp 地震 地震波観測 (7) DGAR(IGS) - DGAR(GSN-IRIS) (d=18.2km) (0.01Hz- ) http:/gpspp.sakura.ne.jp 地震 地震波観測 (8) DGAR(IGS) - DGAR(GSN-IRIS) (d=18.2km) (0Hz-0.01Hz) http:/gpspp.sakura.ne.jp HR-GPS/GNSS衛星時計推定 HR-PPP手法 HR-PPP応用例 地殻変動解析、地震波観測 今後課題: 、長周期低減、高速化 http:/gpspp.sakura.ne.jp http:/gpspp.sakura.ne.jp IaLdPgSkVnYq$t*w-A1D5G8JbNeQhTlWo#r%u(y+B2E6H9LcOfRjUmXp!s&w)z0C4F7IaMdPhSkVnZq$t*x-A2D5G8KbNeQiTlWo#r%v(y+B3E6H9LcOgRjUmYp!s&w)z1C4F7JaMdPhSkWnZq$u*x-A2D5H8KbNfQiTlXo#s%v(y0B3E6I9LdOgRjVmYp!t&w)z1C4G7JaMePhSkWnZr$u*x+A2D5H8KcNfQiUlXo#s%v)y0B3F6I9LdOgSjVmYq!t&w-z1D4G7JbMePhTkWoZr$u(x+A2E5H9KcNfRiUlXp#s%v)y0C3F6IaLdOgSjVnYq!t*w-z1D4G8JbMeQhTkWoZr%u(x+B2E5H9KcOfRiUmXp#s&v)z0C3F7IaLdPgSkVnYq$t*w- A1D4G8JbNeQhTlWoZr%u(y+B2E6H9KcOfRjUmXp!s&v)z0C4F7IaMdPgSkVnZq$t*x-A1D5G8KbNeQiTlWo#r%v(y+B3E6H9LcOfRjUmYp!s&w)z0C4F7JaMdPhSkVnZq$u*x-A2D5G8KbNfQiTlXo#r%v(y0B3E6I9LcOgRjVmYp!t&w)z1C4G7JaMePhSkWnZr$u*x+A2D5H8KbNfQiUlXo#s%v(y0B3F6I9LdOgRjVmYq!t&w-z1C4G7JbMePhTkWnZr$u(x+A2E5H8KcNfRiUlXp#s%v)y0C3F6IaLdOgSjVmYq!t*w-z1D4G7JbMeQhTkWoZr$u(x+B2E5H9KcNfRiUmXp#s&v)y0C3F7IaLdPgSjVnYq$t*w- A1D4G8JbNeQhTlWoZr%u(x+B2E6H9KcOfRiUmXp!s&v)z0C3F7IaMdPgSkVnYq$t*x-A1D5G8JbNeQiTlWo#r%u(y+B3E6H9LcOfRjUmYp!s&w)z0C4F7JaMdPhSkVnZq$t*x-A2D5G8KbNeQiTlXo#r%v(y+B3E6I9LcOgRjUmYp!t&w)z1C4F7JaMePhSkWnZq$u*x+A2D5H8KbNfQiUlXo#s%v(y0B3E6I9LdOgRjVmYp!t&w-z1C4G7JaMePhTkWnZr$u*x+A2E5H8KcNfQiUlXp#s%v)y0B3F6IaLdOgSjVmYq!t*w-z1D4G7JbMeQhTkWoZr$u(x+A2E5H9KcNfRiUlXp#s&v)y0C3F6IaLdPgSjVnYq!t*w-A1D4G8JbMeQhTlWoZr%u0B3F6IaLdOgSjVmYq!t&w- z1D4G7JbMePhTkWoZr$u(x+A2E5H9KcNfRiUlXp#s&v)y0C3F6IaLdPgSjVnYq!t*w-A1D4G8JbMeQhTlWoZr%u(x+B2E6H9KcOfRiUmXp#s&v)z0C3F7IaLdPgSkVnYq$t*w-A1D5G8JbNeQhTlWo#r%u(y+B2E6H9LcOfRjUmXp!s&w)z0C4F7IaMdPhSkVnZq$t*x-A1D5G8KbNeQiTlWo#r%v(y+B3E6H9LcOgRjUmYp!s&w)z1C4F7JaMdPhSkWnZq$u*x-A2D5H8KbNfQiTlXo#s%v(y0B3E6I9LcOgRjVmYp!t&w)z1C4G7JaMePhSkWnZr$u*x+A2D5H8KcNfQiUlXo#s%v)y0B3F6I9LdOgSjVmYq!t&w- z1D4G7JbMePhTkWoZr$u(x+A2E5H8KcNfRiUlXp#s%v)y0C3F6IaLdOgSjVnYq!t*w-z1D4G8JbMeQhTkWoZr%u(x+B2E5H9KcOfRiUmXp#s&v)z0C3F7IaLdPgSjVnYq$t*w-A1D4G8JbNeQhTlWoZr%u(y+B2E6H9KcOfRjUmXp!s&v)z0C4F7IaMdPgSkVnZq$t*x-A1D5G8KbNeQiTlWo#r%u(y+B3E6H9LcOVnYq$t*w-A1D4G8JbNeQhTlWoZr%u(y+B2E6H9KcOfRjUmXp!s&v)z0C4F7IaMdPgSkVnZq$t*x-A1D5G8JbNeQiTlWo#r%u(y+B3E6H9LcOfRjUmYp!s&w)z0C4F7JaMdPhSkVnZq$u*x- A2D5G8KbNfQiTlXo#r%v(y0B3E6I9LcOgRjUmYp!t&w)z1C4F7JaMePhSkWnZq$u*x+A2D5H8KbNfQiUlXo#s%v(y0B3F6I9LdOgRjVmYq!t&w-z1C4G7JbMePhTkWnZr$u(x+A2E5H8KcNfQiUlXp#s%v)y0B3F6IaLdOgSjVmYq!t*w-z1D4G7JbMeQhTkWoZr$u(x+B2E5H9KcNfRiUmXp#s&v)y0C3F7IaLdPgSjVnYq!t*w-A1D4G8JbMeQhTlWoZr%u(x+B2E6H9KcOfRiUmXp!s&v)z0C3F7IaMdPgSkVnYq$t*x-A1D5G8JbNeQiTlWo#r%u(y+B2E6H9LcOfRjUmXp!s&w)z0C4F7IaMdPhSkVnZq$t*x- A2D5G8KbNeQiTlXo#r%v(y+B3E6I9LcOgRjUmYp!t&w)z1C4F7JaMePhSkWnZq$u*x-A2D5H8KbNfQiTlXo#s%v(y0B3E6I9LdOgRjVmYp!t&w-z1C4G7JaMePhTkWnZr$u*x+A2E5H8KcNfQiUlXp#s%v)y0B3F6I9LdOgSjVmYu*x-A2D5H8KbNfQiTlXo#s%v(y0B3E6I9LdOgRjVmYp!t&w-z1C4G7JaMePhTkWnZr$u*x+A2E5H8KcNfQiUlXo#s%v)y0B3F6I9LdOgSjVmYq!t&w-z1D4G7JbMePhTkWoZr$u(x+A2E5H9KcNfRiUlXp#s&v)y0C3F6IaLdPgSjVnYq!t*w-z1D4G8JbMeQhTkWoZr%u(x+B2E5H9KcOfRiUmXp#s&v)z0C3F7IaLdPgSkVnYq$t*w- A1D5G8JbNeQhTlWo#r%u(y+B2E6H9LcOfRjUmXp!s&v)z0C4F7IaMdPgSkVnZq$t*x-A1D5G8KbNeQiTlWo#r%v(y+B3E6H9LcOgRjUmYp!s&w)z1C4F7JaMdPhSkWnZq$u*x-A2D5G8KbNfQiTlXo#r%v(y0B3E6I9LcOgRjVmYp!t&w)z1C4G7JaMePhSkWnZr$u*x+A2D5H8KcNfQiUlXo#s%v)y0B3F6I9LdOgRjVmYq!t&w-z1C4G7JbMePhTkWnZr$u(x+A6I9LcOgRjVmYp!t&w)z1C4G7JaMePhSkWnZr$u*x+A2D5H8KcNfQiUlXo#s%v(y0B3F6I9LdOgRjVmYq!t&w-z1C4G7JbMePhTkWnZr$u(x+A2E5H8KcNfRiUlXp#s%v)y0C3F6IaLdOgSjVnYq!t*w- z1D4G7JbMeQhTkWoZr$u(x+B2E5H9KcNfRiUmXp#s&v)y0C3F7IaLdPgSjVnYq$t*w-A1D4G8JbNeQhTlWoZr%u(y+B2E6H9KcOfRjUmXp!s&v)z0C3F7IaMdPgSkVnYq$t*x-A1D5G8JbNeQiTlWo#r%u(y+B3E6H9LcOfRjUmYp!s&w)z0C4F7JaMdPhSkVnZq$u*x-A2D5G8KbNeQiTlXo#r%v(y+B3E6I9LcOgRjUmYp!t&w)z1C4F7JaMePhSkWnZq$u*x+A2D5H8KbNfQiUlXo#s%v(y0B3F6I9LdOgRjVmYp!t&w-z1C4G7JaMePhTkWnZr$u*x+A2E5H8KcNfQYp!t&w)z1C4F7JaMePhSkWnZq$u*x+A2D5H8KbNfQiUlXo#s%v(y0B3E6I9LdOgRjVmYp!t&w- z1C4G7JaMePhTkWnZr$u*x+A2E5H8KcNfQiUlXp#s%v)y0B3F6IaLdOgSjVmYq!t*w-z1D4G7JbMePhTkWoZr$u(x+A2E5H9KcNfRiUlXp#s&v)y0C3F6IaLdPgSjVnYq!t*w-A1D4G8JbMeQhTlWoZr%u(x+B2E6H9KcOfRiUmXp!s&v)z0C3F7IaLdPgSkVnYq$t*w-A1D5G8JbNeQhTlWo#r%u(y+B2E6H9LcOfRjUmXp!s&w)z0C4F7IaMdPhSkVnZq$t*x-A2D5G8KbNeQiTlWo#r%v(y+F7IaLdPgSkVnYq$t*w-A1D5G8JbNeQhTlWo#r%u(y+B2E6H9LcOfRjUmXp!s&w)z0C4F7IaMdPhSkVnZq$t*x- A1D5G8KbNeQiTlWo#r%v(y+B3E6H9LcOgRjUmYp!s&w)z1C4F7JaMdPhSkWnZq$u*x-A2D5H8KbNfQiTlXo#s%v(y0B3E6I9LcOgRjVm