火山周辺でのGPS観測における数値気象モデルを用いた対流圏補正

全文

(1)῍῎ // (,*+*) + +ῌ+, . ῐ ῏. GPS

(2) = > ? @ , 3A B C D 3E F G H ,**3 7 2 8 ,0 9:; ,**3 7 +, 8 +0 9:<. GPS Baseline Solutions with Tropospheric Correction by Using the JMA Numerical Weather Model for Monitoring Volcanoes ῌῌ Akimichi TAKAGIῌ,ῌῌ, Keiichi FUKUIῌ and Yoshinori SHOJIῌ Japan Meteorological Agency (JMA) has installed and is operating GPS networks around major active volcanoes since ,*** in order to monitor volcanic ground deformation. At present, ,, volcanoes are under continuous observation using about +** GPS stations. At most observation points single-frequency receivers are adopted in consideration of power saving and mobility in rugged environments. GPS baseline solutions include errors due to tropospheric delays which are inhomogeneous in the actual atmosphere. Generally, computation of baseline solutions is done by using a simple atmospheric model assuming horizontal homogeneity. But if the adopted model is not consistent with the actual atmosphere, tropospheric delays cannot be accurately estimated, resulting in poor position estimates. Especially with regard to the volcano observation, the errors in the vertical component of baselines become large according to the large troposphere errors. In the case of baseline solutions between receivers with large vertical di#erence, the time variation of the vertical length is superposed by seasonal noise caused by spatial and temporal variations of refractive index of the atmosphere. For accurate monitoring of volcanic activities, more precise positioning in the vertical component is desirable, which should be realized by incorporating more accurate atmospheric model into the analysis procedure. For this purpose, an improved analysis process was developed, based on the JMA’s operational meso-scale numerical weather analysis (MANAL). The MANAL is applied to daily meso-scale numerical weather prediction as initial field. Generally in the di#erential analysis process of positioning, zenith tropospheric delay (ZTD) is estimated by least-squares method together with the positioning. In this case, initial value of ZTD is given from a simple atmosphere model. In our approach, ZTD between both receiver sites is calculated from MANAL, and then the conventional analysis process is done fixing ZTD between sites. In calculation, analysis software package Bernese Ver. /.* was used, while a part of the program was personally modified. This correction strategy using MANAL was applied to the baseline solutions at Asamayama volcano, where ground deformation has been observed associated with the eruption activity from ,**2 to ,**3. Consequently we could approximately eliminate the vertical seasonal noise at a baseline whose vertical di#erence reaches to +./ km. This approach is quite convenient and e#ective for GPS observation at local and steep areas such as volcanoes. Key words : GPS baseline solution, tropospheric correction, JMA numerical weather model, Asamayama volcano. -*/ῌ**/,

(3) +ῌ+ Seismology and Volcanology Research Department, Meteorological Research Institute, +ῌ+ Nagamine, Tsukuba, Ibaraki, -*/ῌ**/, Japan. ! : +**ῌ23/3 "#$%&'()*+ -ῌ,ῌ, ,-./0123456 Present address : Earthquake and Disaster-Reduction Research Division, Research and Development Bureau, Ministry of Education, Culture, Sports, Science and. Technology - Japan. -*/ῌ**/,

(4) +ῌ+ Forecast Research Department, Meteorological Research Institute, +ῌ+ Nagamine, Tsukuba, Ibaraki, -*/ῌ**/, Japan. Corresponding author : Akimichi Takagi e-mail : [email protected].

(5) &'()*+. 2 + . ,-.$/0/ 123 ,**1 456. 645 789/$ :; <=<. >0? @ ABC 0D!"#A. EF$%6 &B &?G'(H. I)*+ J,- K'..$/ L. M NBOP/QR0S1 / GPS ;. TU !3VW

(6) 45"# 6X 7Y. 52. - %Z0J 8N[9MB \ . @]^_*$/VWD! à:. +/ ;bc 452 : ?G@. <Bde I)*+ AB&/,-JGC. f=0">N0g !3. 0 GPS D!"# h@ 5iG. ,*** ;AB ?06. 0459/$ GPS AB;C D%EF 5. 09 "#. EU Galileo, jkl GLONASS. ,, 45 H +** :; ng$. GPS

(7) GNSS . / .- ; 52 I=opq/J. .- GPS NKrY$L]GNNG. G ,**/ ]m. MN/@ OI=M + PQs RSt u ,

(8) . $/ 456 "vwx= AL Ty U Y@ GPS Sz VWQ JX XY. ,ῌ+ GPS

(9) . TU ABZ&{ [/45N$ \5. |]J G J} GPS J,- %^. 9/$_~; ng \5. Z& ,/02 m 6. *+ABI)*+`B I)*+ J{ %^. 45MB !3 L 5a - : !TUM 5. *+ J{AB

(10) b[/ L cd GPS R0. e - : GPS AB;C fK$/ L.-. @SGIQSz RSt 'gr-rRS0. Y$ ,**0 ;5hi Z& ,.-- m : GPS. /N/R0fK jP-.0/kH IQlV. mAB;: M . n (Fig. +) + PQ. o 0pq N0[! rs p. -* txwx=]uv ; , - u. qwmGNM-. IQ lVo [. ng GPS RSt xI!y MG,+** . ! rs B!z{|. GPS D!. }K. ~YGrs pq G. MG ,+/- xk. ¡ q/;op @ ). IQ lVo AB J|{ ¢. XM. ,**0 ; . ,0 r- ,**2 ; /. ! Y£ !¤F¥- 0¦w§. +/ 1/* . k¨x %^gKABI)gK[/ @. X. ++3 Mg 09 ,**2 ; 2 T4©. I)*+ J{[ª0 452 A;: Z&{AB[/« 9P GPS XYJ nJ} , ;: ! z{ ¬ N$/0P.£ I)* + J,-. AB)`N0 GPS «. D!, I)*+ ®¯0;P $/J} °rP . [!z{. ¡ ¬ g$D!$/0/N ± $/ ) C0[!"# Y£ Z 0[!j²³# ´|$lVo ¬ B0-Jµ¶ n © D!"# " # conventional analysis N·G ]_[!´| "#N¸0J} >0lVo -.¹ /J, 0-0/ Hobiger et al. (,**2).

(11) !"# $ - %. ZQAB GPS IQ lVo µ¶ , J "#. I)*+ °0-. ¹ %^*+YGJ,- º(HM r - %ZQ"# |]JD!» ¼ $ "N½MB +0[!o H. Fig. +. Distribution of the GPS receiver sites around Asamayama volcano. Solid circles and squares are GPS stations installed by MRI and JMA, respectively. A temporary station M. near the summit was newly installed for this study..

(12) GPS !"

(13) #$ % & Table +.. 3. Description of GPS stations around Asamayama volcano.. ' ()*+

(14) ,-./0 123456 789: ;

(15) ;<=9 !" GPS + KWS >?@A B +**0 m # CD ;

(16) BE5 +.,1 m ;

(17) 5 1*3, m Table + F , #GHI GPS J. KL#MN KWS 5 , OP ;

(18) QR5 + 12S#$ : + 12=5 L + (+/1/.., MHz) T7 + 12=5 L+ F=#UN F==N ,ῌ, ῑῒ῏ῐ῎ῌ῍. VWXY Z[\N !" 5 !" $Y ]!#^ _9 F$Y

(19) 5 "`#$ #%a=N#`

(20) (GSM) = GSM b c#&de=9'#fLghNi

(21) (MSM) FF5 j()kl 8mj nopqrst:u#GY

(22) = vw #xy$Y5z{|} *+ z{. Fig. ,. Horizontal distribution of MANAL’s grid points around Asamayama volcano. Open diamonds show +* km-grid MANAL. Solid squares show doublefrequency GPS stations which were used for the ionospheric correction. A square box over Mt. Asamayama indicates the area of GPS observation in Fig. +.. |}5~,!-#-D .$ : :0: F,!-# Y^. &$ 125 !" MSM z{. #QR= ,!/n#j9. # 0$Y 1^_9iQR. z{|}=:QR#^_N QR. (MANAL) 12 (GPV) 2I. 5 3". B4 pqs 56. +* km 72#89$p

(23) #. ¡ ¢£¤ !"¥¦:uf 12. +1§=9$ 129¨©§t ª«¬ ,**,. !"®ps1*1¯1¡ " . GPS °' : 12. f : ,**3 ; . ±O'52I. ¬²< ,**3 #

(24) #³´z{. / km =7

(25) #+1§=9$9¨©§t. |}g0no9$>µ¶8·¸S . ¹x9 '?@? ,**3 O7 MANAL .

(26) A9$QR# Y=y(B 129¨ . GPV # MANAL =º»N. 0 C !"D +331. MANAL ¼E4I5µ85:h "Fh½.

(27) 9:;<=7>?. 4. Fig. -. Schematic figures showing zenith tropospheric delays (ZTD) of two GPS-receiver sites. (a) In the case of small vertical di#erence between two sites, atmosphere equally influences both ZTDs. (b) In the case of large vertical di#erence between two sites, the part of atmospheric between Receiver+ and Receiver, influences only ZTD of the lower sites (ZTD+).. @AB@#CDE 0&8 6. Ver/.* (Dach et al., ,**1) F 0 Bernese GHI. CD!J0&KLM#E JN0&8. OPQR Fortran SN- TUVUGHIOP. "=CDWXY &CD9# Z#. [U\ ]^0&_$ `$ab,=% cd0&.

(28) KWS-.

(29) `$& + '\efghijkl. M. (9) +./ km CD. ++ *m. +. n$ +,- 8 `$F 0opGqpI r -* .20

(30) sr , s. Fig. , 2t_u MANAL v/. 12 z+ 6 345 {. 20 w v/x0

(31) y#M . 8 , ' 6. 3' 7#|}8= n ~w. J Wz+ KWS 9!:3 6. 8 - t

(32) ;Aw<=. n. *+ w 9"vw >81_u -. ,ῌ- MANAL B. 2"!" ? !@A+6

(33) . 0l# C+22W. GPS '()!6_uwDUe B. 40& , m AWA+ GPS MLDL-. 0. 3' 0E$.MLn#. Xl#2F+ % 6 M. . % & GPS '()*+ ,- ./. w 0 wa Onset Computer Corpo-. 0 - 1

(34) (ZTD) 234 (Fig. -a). ration G HOBO Weather Station UG MR. jklab 6 Receiver+ ( H20. 0*** F 0 + sopGqpIw¡ I¢#. 6 Receiver, (22W62 6. 0 ,-LwDUeWJ , £ GPS DUe22. 1

(35) WKH¤20& ,-L "s +. W£¥0 ,**0 J 3 ¦ ,+ L ^§0 ,**2 J. Receiver+ 2 Receiver, 4 M"F9 78X. / ¦ +/ L6 0*, L ,*, L¨L n #. (wGH©ª«NL-w. _ + L +, s O&¨ nL ¨L20. CDabPB¬I-. `$F 0#8,220. -

(36) 56-78 (Fig.. QR 6wDUe ®;S¯YTX!. -a) 0M0 6F9 !8X wCD. A DUe + sU F 0 !A DUe.

(37) !"#$. ¨#8. t°±l²VN1 W+ (Fig. -b) =Xjklab ³ 0wCDML. -

(38) . Y Z'()!6_u1

(39) ab0&_. GPS `$´pµ`$[©¶·¸¹ Bernese. ! t- º»U20&¼½Z#

(40) (.

(41) !" GPS # $%

(42) &'()*+ ,-.. 5. /0&1 &()*+2 34 Saastamoinen (+31-) 5&+26 ˘ ´ ıma , 789 - : (m) + 3, Bauers (+32-) ;+<.6 ZTDSAASpTe . *.**,,11 ῐ ῌ+,// ῎ ῒ *.*/ eBtan,z dR p ῏ ΐ cos z ῑ ῍ T (+). 6, p, T, e =>=?@A&B (hPa) &C (K) DEFG&B (hPa) +H6 z +IJ. K (radian) B dR L%M B NO% dR KNO% M P, p T e ?@A0 )QR+. Fig. .. Schematic figure showing how to estimate ZTD. ZTD+-, (di#erence between ZTD+i and ZTD, i) is calculated from MANAL. ZTD and position at receiver sites are estimated fixing ZTD+-,.. S Berg (+3.2) 7TU+ VN6,O% 3 [:\=. 6 W!]+WY^Z [_R`Qab. (+) 7 ,cde

(43) +fgh=. WX WYZ. i"j ZTD+i ZTD,i 7 (-) kf\=,l. mnT ()*.oT=4# #. $:pq%rs+3, 7 (/).

(44) MANAL GPV +tuvwxyz ]{+|z, Fig. . [}!& Receiver+ ! Receiver, . M

(45) Receiver, :de

(46) ZTD, + *. MANAL GPV +~)6,&'* p,, T,, e, +7 (+) ./6,

(47) + . ZTD+iZTD,ifixZTD+,. i WYZ. (/). $ i '(_R`Qab. i+H T$ GPS IJ89+, +-$ Vj:+ 686:[ , ",%0. ZTD,*ZTDSAASp,T,e,. (,). +# .1 M (Hofmann-Wellenhof et al., ,**.) 686 + ",. ;2 Receiver+ :} Receiver, . GPS )QR

(48) +-:. T . ZTD+, MANAL GPV. ,c + " , # 1 3 4 % , ", GPS # . +~)6,&'*+ Vkf1M l$. )QR+5 6V. MG,+,* 6. M Hopfield (+303) $[E Essen and Froome (+3/+). +& +5 6, = " , ", GPS #. 89. )QR89IJ89+-:+c +- N+. ῌ ῍. ZTD+,+*0 ῌ hj 11.0. j+. MAP _* 7 (*) 7 (l) (t) . p+j e+j +,.30 T+j T+j. e+j ῎ -.1+2 · +* T+j,῏ /. +-:Q* + , 6 =+ 1h 6 + ", GPS # $IJ8 (-). 9+-:+c .+0jM ;8 GPS +5 +-:O9+. \= N Receiver+ Receiver, . %+¡¢ £¤Q¥*T+-¦§. ¨()*:% p+j T+j e+j Receiver+ ;©. ª «Q¤§¬<®¯iRQ °±= 2. ²³. j '(()*:$&B &C FG&BM. # ´+hµV6V . hj Receiver+ ;© j-+ '( j '(()*:. !>+¶·¤Q¸*T?¹ Tv +-¦§ª. º WuW© =>= Receiver + DE Re-. +jM. ceiver, ()*:º\T .ῌ ῏ ῒ ῑ ῎. Receiver+ :de

(49) ZTD+* + ZTD+*ZTD,*ZTD+,. (.). .ῌ+. ZTD ῍ΐῐ. .v + , , »¼]{[ ½.

(50) )*+,-./01234. 6. Fig. /. Temporal changes of ZTD between KWS and M., evaluated by the conventional analysis (a), by using MANAL (b), by using measured weather data (c), and those refer to temporal change by the conventional analysis.. 567"8$ , 9 ZTD !:; Fig. /. <

(51) ='($>?@56A !7" Fig. 0a-b. $ ,**0 ,**2 B KWS A M . 9 ZTD + CD. EFGHI Fig. 0d-e JKHI7" 0LM. NO7" Fig. /a EF?@PAQR$ST

(52). U9VW JKXYBZH[\]^ _. %&'($> ZTD, Fig. /b MANAL GPV . $àbcde %&'( fghAW . cm ij.

(53) ZTD Fig. /c , . BZklVm

(54) (Fig. 0d) V ='($>?@5.

(55) ZTD !"

(56) #$ %&'(. 6 nopqrcs (Fig. 0e) EFGHI tr. ST

(57) 56Au Fig. /d $7

(58) MANAL . gvs gs (Fig. 0a-b). GPV

(59) ZTD wxy$ , cm ijtrhg. AVz {

(60) , . .ῌ-

(61) . |$ =}~QR$

(62) , . b

(63) ZTD wxy$ + . cm 1bhge BZH. , 9 (-). [V csAVz A. OT (/) cEF?@

(64) (Fig. 0f) ?@. MANAL GPV

(65) P56 %&'(. $ MU9VWBZH[V . <

(66) ?@$>56A g5. ,**1 BAXY?@56$ ,- cm. 6Vb. klVAVz %&'($> . cm ij. .ῌ,

(67) . klBZ\]^ (Fig. 0d) >ebcs. K ='(<sc T ¡R + C ,. R9. W MANAL ¢<

(68) '( (Fig. 0e) $ £g. <

(69) ?@567" GPS + CY. m

(70) {W

(71) . {A#c?@$<

(72) V ZTD ST , R9¤$. Saastamoinen (+31-) ¥ῌ$¦§ cP. :mcs (,) ZTD,* (-) ZTD+, + C. :;>eW ¨©ygª« ¨©t. DNO<

(73) . ¥ῌ$¦§ cP :;¬V56®¯V. Fig. 0 $ KWS E©A

(74) M. ~%&'(. >sAV°±bcs ²³´ ,**, A. $>EF?@56A µO MANAL ¢. OV¶·¸¹º»t$x"¼½¯.

(75) GPS

(76) . 7. Fig. 0. Time series of baseline length (a-c) and height (d-f) of M. refer to KWS. (a) and (d) are calculated with the conventional analysis. (b) and (e) are calculated with this study method using MANAL. (c) and (f) are calculated using measured weather data. In the case of using MANAL, vertical seasonal noise could be approximately eliminated (e).. !"#$%&' GPS ( ). *+,-. /0123% , 4567 + . 7 5689:( ;<=$. 456%>?"@A 7# BC4=DE. % FGHIJ K$%LM'". N 74 OP-QRST=UGK. .ῌ.. , . G VW + X 48YZ , X 4['". \7 + 4] GPS 56 ^_à'"H. ST=FGbcGde# $ ) f + 4. IJ K$%L' g , 4] GPS 56. 56 h ij8k,-. lm7 , 45. '"Zn' wdx (q GPS rs%'" y z. 6 op%>?"(qe $ rst'"7 / uvc (q {| lm7

(77) - cm }| lm7

(78) , cm. ~K , 4] r. \7 , 456Zn KG k,. r 30*0+. wd 3/-. -. /0 HI7 + 456 n@A 7. m

(79)

(80) TAN ,. #!. wd ,+3* m %' (Fig. 1c) \[!" 30*0+. . .ῌ/

(81) . = TAN DE Nij! , rs7. $p 7 ,. % 56'"ij!. 202* m wdx7 +,-1 m K ¡7 ,**1 k . ¢. £IL' + ¤ ¥¦%'"ij'£I K. § ,**3 k - ¢ , k K. G MANAL 7 - % ij aiZX. , 456=ij£I Fig. 1 a , 4. ,. ij¨©ª«¨8YZ¬ ij. 56 ® L+ ¯ + 456 °['ij£. $%8Y ¬¡ ±#(²³´µ q. I ¶·° Fig. 1b L= h \ £I. ) NbcGde#¸GK%¹º o$. Z op»p 9¼L' \8° k. \"½

(82) +, ¨©ª.

(83) BCDEFGHI. 8. Fig. 1. Comparison of seasonal noise’s elimination from vertical positioning between doublefrequency analysis (a) and single-frequency analysis (b). (c) Distribution of GPS receivers around Mt. Ontakesan. Positioning of TAN refer to 30*0+. was calculated.. JK LMN# ,. >?9:OPQ6RST' +, >?9:UVKW XYZQ[\&)]^# MANAL >? 9_ à)b<c.d GPS . eX fg6chi jklm nogpq p*6chi * =&'d -* GPS !"# 6&'V>?9: rs t2 EU MANAL - >?uvw* 4x eX5? ZTD $%&'()*+ + ,-). / #yz + cm {|U-a) 786ca'# , cm. }. ~*d-. MANAL TwpZ rs6 '

(84) vw U- /

(85). ?U#eX ta) M.

(86) 6 Fig. 2. Comparison of seasonal noise’s elimination from vertical positioning between ,.-hours analysis (a) and +,-hours analysis (b) by using MANAL. M. refer to KWS.. 6 - 5 6 - 5 GPS eX 01&' * 23 KWS. 4/56&)78.

(87) 0 5 eX5 # --.+.,1 m ,**0 J / vw ,**3 J 1 Ù . 9:6c*w @ >. Fig. 3 6=. jk 0 >?U 9: () Fig. 2 6. M.

(88) G ++02 m. M. 9:;< = +, >?9:#@A&) 0 >?. TAK 31+ m U RS#q6U-. 9:# GPS c.gpqp)9`wp. * W??# o ¡\& ,**2 J 2 . 78¢va) U **U#£&) +, >?9:. ,**3 J , F¤¡¥&) *6¦a. (Fig. 2b) # ,. >?9: Fig. 2a, Fig. 0e §¨ 6©ª. '«odeXT) 6c9:;<vw. 'VKW LMN ¬f*WT) ;<. KWS. #´ µU-.VKW ¶ ·¸¹º»&) E. -a) ,**2 J 2 ¼vw ,**3 J , ¼Ù#½. 4/ o#®¯W6#°±²³U.

(89) ! GPS "#$%&'()*+ ,

(90) ῍-.ῐ/01. 9. Fig. 3. Time series of height of all stations refer to KWS around Asamayama volcano from ,**0 to ,**3. (a) Conventional analysis. (b) By the method of this study using MANAL. Those of all stations could be approximately eliminated the vertical seasonal noise, especially G and TAK where are high.. 234 55 6῎78!9:;. -<=>=?@ABCDEFG5 HI . JKLM!N-OPQRSTUV5<%W . X ,**/ YZ[\]!^$῕_àb! c. de Fig. 3a ! G, TAK BffPQΐgPQh. i$4f Z[!=>=jYW4. k$l-<!mΐn_op

(91) qKrs4. tdfB<-CD4. !;J

(92) $ Z[ῑu"#AB\v!. !k[!w\PQ

(93) 6῎78#xyz5;J. ,**2 O 0 {f PQ|}5 ,**3 O + {Y! -. (Fig. 3a) fK w\PQ!mΐn_opKrsK4. ~Pxf I+ , (Mogi, +3/2)

(94) ῍5<. 6῎78#xyz5;J

(95) $ῌ!7 . ῔am

(96) 9:$K

(97) tῒ5<+ ,. ῔a !m

(98) (Fig. +*b) (K. m

(99) d MaGCAP-V X ,**/

(100) ῍- . 0178

(101) ῍5C x ,.1 km. [ +** m @P}ΐ![ *.++*0 m- . ¡3$ cAK¢£] 1.- km ¤ - ¥!cA$_. ¦8$W §!;J Z[¨! 0km !. W Y@©ªΐK¤ / ¥! +0..*+*0 m- _. ¢£] ,.3 km «¬ -.*+* m !@©ª4Be". !$῔ac®@fB<5Y῏¯K. #

(102) (Fig. +*. °±k[ ,**2 O 0 {,**3 O + { $ ²³<-. a) !CDK ,**. OZ[´$µ«¶PQ$%. _-·]¡3!OPQRST!¸!¹º»¼½k. 0. -.

(103) !"#$. 10. Fig. +*. Estimation of the inflation pressure source of Asamayama volcano by GPS observation from June, ,**2 to January, ,**3. (a) Result estimated by --D component of displacement with the method of this study using MANAL. (b) Result estimated by --D displacement with the conventional analysis. (c) Result estimated by only horizontal component of displacement with the conventional analysis.. %&'()*

(104) +, -./0(12345678&.

(105) Fig. 3a TAK G OPQ. BC2( GPS [\(]^=-._VWX>?. 9:; <=>?@AB(: BC2(DEF GPS GH5IJK7LMN )R) BSTU>?@AVWXY@UZ +ῌ -- ` _a56bZ +ῌ-/ `cEd ) >. U GHJKefO&gFB)h ijJKk*. ?@AVWXY@_a56bl(Umnopq. rsABtuF Fig. +*c OijJKk*7vM. wxyO zB {(VWXY@7|}~(. VWX>? VWXo . Co)h? 4 ,.3 km (?

(106) . UGHJK7STijI +., km (. }~(C - JK5IK)hVWX. Y@U4 *.2 km :;Z + ῌ - cEY@( . o7>? z VWX_a56bU. _abU +ῌ- +.*+*0 m- : BSTU. -.1+*0 m- :; |}~(C)h>?@A. 34JK7 - JK5IK)h>? U. -.*+*0 m- Z , &g=athA .

(107) ) GPS *+,-. / 01. 11. 2+ 3 ,..+* m 45$. 67 !" !%#"#) ,849:;<. , =$>?%%@6. ;&'( 503 )*>6,AB. 0. -. +*+,,-ABC.D E5:F . />. +ῌ-- 0 +%1G). +ῌ-/ 0H. .D39 0 +923I456+4E67. 0 7J&0 AB8K 9:LM5084. AB;NO.1PAB.<Q0RST=#. 9:;< 86@;@ O.U>= GPS

(108). >6 V31PAB.<)WX45. 8K LN<)%?')Y,84 @. , Z6[M O.\]^ A_O.1P. . `?%%,84 AB)* LN<)BC D)Y, ; 86) GPS P*+,aEF. ΐ. ῒ. 5G ;@ HU9U0IRS b$ C. ,c *0& 3 J&L:6. Q0IRSO.\]^K ;< L<. MN O PQd <Ref STUV. & GPS O.1 9:?%%,gW5&'(.

(109) . hX4 8KYi?')Y,. Zj[\F]k&0 ^lKm,c . ;8'( hXE,0b;. $ n o_ 0. 9Y$ p`9aD)Y, q4b

(110) bcr. stcduvefgwhistx. yz{y|UjG)Y, }k -* lryz{. st~m*0&n5o 0 pqst-. y|. !") ,. bc

(111) >=9U0849:;<. V - !"r$0b0&. 9 Zs +, bcDbc y

(112). :tuovw0 xy-z{|# }~. )F5P>=9F MANAL %b - bc. y 2) o. o4

(113) &bcjU0 GPS

(114) 9. 0 { ywI !" >. }E,j69Y$ Z5". &0 - GPS.

(115) 9) ,Uryz{y|) GPS. !" %. w dQf #* 5< . %Nq4b - !"rJb6). I*0&o_0 )Y,SN . xuv^

(116) . &pÈ,aD9Y, Y, ¡

(117) )Y, MANAL . !"b0&$N. !"¢9. 4Vo#<& & ¡. ) ,)*#Z - bc;;, #<& {/". £>6 88¤&¥5>¦T¢§. ]¨hX&'(#, GPS P

(118) E,. £E.

(119) b©9 - bcP Y,0 0 bcP MANAL. ,c. -*+,- ^lKm{¤s. ¥4E,pq (MSM) E,aD9Y$. tª¦st «§{.!¢y¬y|#,{¤. Z6 0

(120) &®¯p`9aD)Y,.

(121) ,**0-B-*0 . 0ῌ ῎ ῍ ῏. -±²!/¡

(122) - !")Y, MANAL 0& O. U>=P´¯E,'. ( µ° 8'( +2w±HU²9K 03*0&'(9Y$ ,c *+, ¶ )¥;+®1 \]^9R·³¸) , GPS. !"#$% + !"&#$'(. F)Y,849:;< ¦b*- !".

(123) 4´¹&% MANAL

(124) 39>=U0º(9J@6 MANAL bcj. E, GPS

(125) bcj U, p`9®¯aD)Y, &'( 0

(126) %4 ,**2 ®;@,c {³1 AB )*48» ,c. µ¶ 0 km ½¼D. - km {I) - +*0 m- +%½9Y. ῐ ῕ ῔ ῑ ¨©ª° « $ n uQ ¬ RV (,**/) ,**.ῌ,**/ ®,c ;N{³1 /* /1/ῌ/2.. Bauer˘ sı ´ma, I. (+32-) NAVSTAR/Global Positioning System (GPS), II., Mitteilungen der Satelliten-Beobachtungsstation Zimmerwald. +*, Astronomical Institute, University of Berne. Berg, H. (+3.2) Allgemeine Meteorologie, D¨ ummler’s Verlag, Bonn. Dach R., U. Hugentobler, P. Fridez and M. Meindl ed. (,**1) Bernese GPS Software Version /.*. Astronomical Institute of University of Berne, 0+,p. Essen, L. and K.D. Froome (+3/+) The refractive indices and dielectric constants of air and its principal constituents at ,.*** Mc/s. Proceedings of Physical Society, 0. (B), 20,ῌ21/. · V &¸¹ ºU©» ¼L.

(127) 12. ῏ !"#$%. &' (,**/) ()*+,-./01234567 DEF (,) +GῐHIJK῎JLMNL OP V*//P-*--. Hobiger, T., Ichikawa, R., Takasu, T., Koyama, Y. and Kondo, T. (,**2) Ray-traced troposphere slant delays for precise point positioning. Earth Planets Space, 0*, e+ῌ e.. UVW5XYZ9[!\]^_àbcd efgh ijklm ")n; opq῏ (,**1) rstuCvwxNyz{|} GEONET ~s /0* +JL OP +*2 P-2. Hofmann-Wellenhof, B., Lichtenegger, H. and Collins, J. (,**.) Global Positioning System : theory and practice, /th, rev. ed., Springer Wien New York, 33ῌ+*/. ¢£¤¥ῑ£ ¦ (,**3) §¨ῌJ>3 . ?@ABC ª «¬ ,+ ®

(128) ¯°Y5 .

(129) 0/ῌ1+. Hopfield, H.S. (+303) Two-quadratic tropospheric refractivity profile for correcting satellite data. J. Geophys. Res., 1., ..21ῌ..33.. 89:;"< = (,**,) >3 . ?@ABC

(130) Q

(131) RST .2 -1ῌ/3.

(132) (+331)

(133) Q

(134) RST .- +00pp. Mogi, K. (+3/2) Relations between the eruptions of various volcanoes and the deformations of the ground surfaces around them. Bull. Earthq. Res. Inst., -0, 33ῌ+-.. Saastamoinen. J. (+31-) Contribution to the theory of atmospheric refraction. Bulletin Ge óde ´sique, +*1, +-ῌ-.. "#$%dq\]^_ ;Z9[! N ῒ (,**3) GPS J : GPS

(135) Y J +JJ *vKx +JL // +1ῌ-2 a (,**,) GPS {~ ¡¢ {©L ,03p. )± « (,**/) o²()³ ῍´ µ¶³ ·]¸v () /* S1ῌS+2. ¹Pº» ¼½.

(136)