http://www.asyura2.com/11/jisin17/msg/429.html
| Tweet |
10年くらいの巨大地震が集中する期間と、安定期があるが
ポアソン分布では説明できないだろう
プレート内の巨大な断層の崩壊が連動するからか
今後の測定手法の進化とネットワークの高密度化で、地震前後のダイナミクスの詳細な記述が可能になり警報精度も上がる
http://www.sciencemag.org/content/332/6036/1390.full?sid=80d39ae3-33b9-4e27-bea4-7221024ac777
Science 17 June 2011:
Vol. 332 no. 6036 pp. 1390-1391
DOI: 10.1126/science.1206643
Perspective
Geophysics
A Tale of Two Earthquakes
Kosuke Heki
+ Author Affiliations
Department of Natural History Sciences, Hokkaido University, Sapporo 060-0810, Japan.
E-mail: heki@mail.sci.hokudai.ac.jp
Thousands of lives were saved by the “Large Tsunami Warning” announcements of the Japan Meteorological Agency that followed the Tohoku-Oki earthquake of 11 March 2011. A similar alert followed the Central Chile (Maule) earthquake on 27 February 2010. Understanding better how such earthquakes develop might provide even more effective early-warning systems to help mitigate their devastating power. Four papers in this issue―by Sato et al. on page 1395 (1), Simons et al. on page 1421 (2), Ide et al. on page 1426 (3), and Vigny et al. on page 1417 (4)―report on both of these magnitude 9 (M9) earthquakes, illustrating the use of networks of Global Positioning System (GPS) detectors to reveal how the Earth's surface deformed during and after the events.
Four M9 earthquakes hit Kamchatka, Aleutian Islands, Chile, and Alaska in a relatively short period between 1952 and 1964. After remaining silent for more than 40 years, M9 earthquakes have resumed: the 2004 Sumatra-Andaman (M9.2) (5), 2010 Maule (M8.8), and 2011 Tohoku-Oki (M9.0) earthquakes. The last two are the first M9 events to be studied by dense networks of continuous GPS observing stations. The region of the Chilean subduction zone reported on by Vigny et al. is known as the Concepción-Constitución seismic gap (6), because it fits the “seismic gap” hypothesis―that the recurrence history of past earthquakes suggests imminent rupture. Thus, multiple international teams deployed GPS stations in the area, resulting in a dense network of seismic detectors.
Tsunamis are caused by coseismic vertical movements of the sea floor. Coseismic crustal movements in the 2010 Maule earthquake were observed at ∼90 GPS stations and observed continuously at ∼60 sites (4) (see the figure). The earthquake released east-west compressional strain accumulated in South America by the eastward subducting Nazca Plate. Meters of westward displacements occurred along the Chilean coast with smaller displacements extending across the continent to the Atlantic coast. In the 2011 Tohoku-Oki earthquake, eastward displacements extended as far as Kyushu in southwest Japan. In Chile, because the boundary between coseismic uplift and subsidence roughly coincides with the coastline, coastal towns experienced relatively small vertical movements. However, the northeastern Japan shoreline was farther away from the trench, and coastal lowlands already suffer from inundation at high tide caused by the coseismic subsidence of this earthquake. This subsidence will not be compensated by future interseismic uplift (7).
The largest displacement of a GPS station on land was about 5 m in both earthquakes (southwest of Concepción in Chile, and in the Oshika Peninsula in northeast Japan). In Japan, measurements of sea-floor positioning revealed eastward coseismic movement of 24 m near the epicenter, ∼100 km off the coast (1). This, together with the 31-m displacement east-southeast (not included in the figure), measured by the Tohoku University group at a submarine benchmark ∼175 km offshore, would be the world record for measured coseismic displacement.
Large interplate earthquakes are often followed by silent fault slip (afterslip), and this causes postseismic crustal movement in the same direction as that of coseismic jumps (8). Postseismic movement of ∼15 cm was recorded in Concepción in the first 12 days after the Maule earthquake (4). In Japan, postseismic movement in the same period reached 20 cm at Yamada, Iwate. The difference in displacement between co- and postseismic movements suggests that coseismic slips concentrate on discrete asperities (the area where most seismic energy is radiated), whereas afterslip is distributed more uniformly within the ruptured plane. In Chile, such a difference is seen in the station to the north of the fault that moved southwestward during the earthquake but westward after the earthquake.
View larger version:
In this window
In a new window
Download PowerPoint Slide for Teaching
Monitoring movement.
Coseismic horizontal movements of GPS stations of the two M9 earthquakes, i.e., the 2011 Tohoku-Oki earthquake (1–3) (left) (coordinate difference between 10 and 12 March data by Geospatial Information Authority of Japan available at www.gsi.go.jp/ENGLISH; only ∼10% of the whole network is shown for visual clarity) and the 2010 Maule earthquake (right) (4) in the same scale. In addition to terrestrial GPS stations (red arrows), coseismic displacement of the five submarine geodetic benchmarks (black arrows) were added (1). The bottom panel compares the 12 days of postseismic movements (blue arrows) in the two earthquakes. In Japan, the fixed reference is 0747, in southern Ryukyu (outside of the map).
GPS receivers usually record phases of microwave signals from satellites every 30 s. In Chile, high-rate (reading the phase every second) sampling was taken at eight stations. This makes a GPS receiver a seismometer to directly record displacement, which has so far been obtained by integrating acceleration or velocity records of conventional seismometers. High-rate GPS receivers showed the complex history of the movement of the San Javier station, east of Constitución, during the ∼2-min period of the earthquake. The station first moved southwestward, then turned right moving northward, and suddenly turned left ∼1 min after the rupture start, ending up with cumulative movement of ∼3 m toward west-northwest. Such displacement seismograms are useful to infer rupture propagation speed and rupture sequence of multiple asperities. The Japanese GPS network also records 1-Hz sampling routinely at most of the stations. This technique will become a popular tool to study earthquake source mechanics.
A difference between the two earthquakes is that the seismic gap hypothesis was not valid in the Japanese case. Seismologists had anticipated recurrence of the 1978 M7.4 Miyagi-Oki event instead of a M9 event (9). We should not have complete confidence in seismic gaps, but it would be fair to say that its basic concept is not wrong.
Crustal deformation causes several geophysical phenomena. Vertical deformation of layer boundaries with density contrasts [e.g., sea floor and Moho (base of the Earth's crust)] and changes in rock density around the fault produce subtle changes in gravity. The 2004 Sumatra-Andaman earthquake was the first earthquake whose gravity change was detected by satellite gravimetry (10). The Maule earthquake was the second (11), and the Tohoku-Oki earthquake will be the third. The combination of dense detector networks and the suite of new techniques will provide a picture of earthquakes and their dynamics never before available.
Received for publication 25 April 2011.
Accepted for publication 2 June 2011.
References
↵
M. Sato
et al
., Science 332, 1395 (2011); 10.1126/science.1207401.
Abstract/FREE Full Text
↵
M. Simons
et al
., Science 332, 1421 (2011); 10.1126/science.1206731.
Abstract/FREE Full Text
↵
S. Ide,
A. Baltay,
G. C. Beroza
, Science 332, 1426 (2011); 10.1126/science.1207020.
Abstract/FREE Full Text
↵
C. Vigny
et al
., Science 332, 1417 (2011); 10.1126/science.1204132.
Abstract/FREE Full Text
↵
P. Banerjee,
F. F. Pollitz,
R. Bürgmann
, Science 308, 1769 (2005).
Abstract/FREE Full Text
↵
M. Moreno,
M. Rosenau,
O. Oncken
, Nature 467, 198 (2010).
CrossRefMedlineWeb of Science
↵
K. Heki
, Earth Planet. Sci. Lett. 219, 13 (2004).
CrossRef
↵
K. Heki,
S. Miyazaki,
H. Tsuji
, Nature 386, 595 (1997).
CrossRef
↵
R. J. Geller
, Nature 472, 407 (2011).
CrossRefMedlineWeb of Science
↵
S.-C. Han,
C. K. Shum,
M. Bevis,
C. Ji,
C.-Y. Kuo
, Science 313, 658 (2006).
Abstract/FREE Full Text
↵
K. Heki,
K. Matsuo
, Geophys. Res. Lett. 37, L24306 (2010).
CrossRef
The editors suggest the following Related Resources on Science sites
In
The 2011 Magnitude 9.0 Tohoku-Oki Earthquake: Mosaicking the Megathrust from Seconds to Centuries
Mark Simons,
Sarah E. Minson,
Anthony Sladen,
Francisco Ortega,
Junle Jiang,
Susan E. Owen,
Lingsen Meng,
Jean-Paul Ampuero,
Shengji Wei,
Risheng Chu,
Donald V. Helmberger,
Hiroo Kanamori,
Eric Hetland,
Angelyn W. Moore,
and Frank H. Webb
Science 17 June 2011: 1421-1425.Published online 19 May 2011
Abstract
Full Text
Full Text (PDF)
Supporting Online Material
In Science Magazine
The 2010 Mw 8.8 Maule Megathrust Earthquake of Central Chile, Monitored by GPS
C. Vigny,
A. Socquet,
S. Peyrat,
J.-C. Ruegg,
M. Métois,
R. Madariaga,
S. Morvan,
M. Lancieri,
R. Lacassin,
J. Campos,
D. Carrizo,
M. Bejar-Pizarro,
S. Barrientos,
R. Armijo,
C. Aranda,
M.-C. Valderas-Bermejo,
I. Ortega,
F. Bondoux,
S. Baize,
H. Lyon-Caen,
A. Pavez,
J. P. Vilotte,
M. Bevis,
B. Brooks,
R. Smalley,
H. Parra,
J.-C. Baez,
M. Blanco,
S. Cimbaro,
and E. Kendrick
Science 17 June 2011: 1417-1421.Published online 28 April 2011
Abstract
Full Text
Full Text (PDF)
Supporting Online Material
Shallow Dynamic Overshoot and Energetic Deep Rupture in the 2011 Mw 9.0 Tohoku-Oki Earthquake
Satoshi Ide,
Annemarie Baltay,
and Gregory C. Beroza
Science 17 June 2011: 1426-1429.Published online 19 May 2011
Abstract
Full Text
Full Text (PDF)
Supporting Online Material
Displacement Above the Hypocenter of the 2011 Tohoku-Oki Earthquake
Mariko Sato,
Tadashi Ishikawa,
Naoto Ujihara,
Shigeru Yoshida,
Masayuki Fujita,
Masashi Mochizuki,
and Akira Asada
Science 17 June 2011: 1395.Published online 19 May 2011
Abstract
Full Text
Full Text (PDF)
Supporting Online Material
http://www.ep.sci.hokudai.ac.jp/~heki/jindex.htm
測地学会誌 第56巻第4号, pp125-134, 2010 解説・入門講座
GPS-TEC法による地球物理学
日置 幸介1，菅原守1,2，大関優1，岡崎郁也1
1北海道大学理学院自然史科学専攻
2NTTデータ基盤システム事業本部
（2010年6月24日受付，2010年9月14日改訂，2010年9月24日受理，2011年2月25日出版）
Geophysics with GPS-TEC
Kosuke Heki1, Mamoru Sugawara1,2, Masaru Ozeki1 and Ikuya Okazaki1
1Dept. Natural History Sciences, Hokkaido University,
N10 W8, Kita-ku, Sapporo, Hokkaido 060-0810, Japan
2System Engineering Business Unit, NTT Data Corporation
Toyosu Center Bldg. Annex, 3-9, Toyosu 3-chome, Koto-ku, Tokyo 135-9671, Japan
（Received June 24, 2010; Revised Sep. 14, 2010; Accepted Sep. 24, 2010, Published Feb. 25, 2011）
Abstract
Two carrier frequencies, L1 and L2, used in Global Positioning System (GPS) can be used to isolate ionospheric total electron contents (TEC) by taking the phase differences between the two carriers (L4). Here we review technical aspects of geophysical applications of TEC measurements. Topics include conversion of standard format raw GPS data files into L4 and TEC, and numerical techniques to isolate disturbance signals with various time scales. We also introduce typical examples of TEC change signatures of variety of geophysical phenomena, e.g. volcanic eruptions, earthquakes, solar flares, ionospheric hole formation by rockets and missiles, solar eclipse, and elusive precursory TEC changes before earthquakes.
1. はじめに
全地球測位システム(Global Positioning System, GPS)は，衛星が発射するLバンドのマイクロ波を地上局で受信し，その搬送波位相の変化から高精度の測位を行うシステムである．
我が国では地殻変動の計測を目的とした稠密な連続観測網GEONET (GPS Earth Observation Network) が国土地理院によって展開されている．搬送波は1.5 GHz (L1)および1. 2 GHz (L2) の二つの周波数が同時に用いられる．電離圏と呼ばれる超高層大気では，太陽放射によって大気分子・原子の一部が電離しており，電子がマイクロ波の伝搬を遅延（電離圏遅延）させる．GPS受信機は，二周波を同時受信することにより，周波数の二乗に逆比例する電離圏遅延を除去しているのである．その時に用いられるL1とL2の線形結合（Ionospheric-free Combination，しばしばL3と呼ばれる）はそれぞれの搬送波周波数をf1，f2とすると，
L3 = f12/(f12−f22) L1 − f22/(f12−f22) L2 (1)
であらわされる．ここでL1, L2は，観測された位相に波長を掛けて長さ(m)の単位にしている．一方L1とL2の単純な差（ここではL4と呼ぶ）は，電離圏の情報のみを含んでおり，中性大気遅延や衛星位置，局位置などの情報は差をとった時点で除かれている．そのためGeometry-free Combination（幾何学的要素が含まれない線形結合）と呼ばれる．
L4 = L1 – L2 (2)
L4は視線に沿って積分した電子の数に比例するため，主に電離圏研究に用いられてきたが，地震や火山噴火などの固体地球に関連した現象に起因する電離圏擾乱の研究も盛んになりつつある．本稿では，これまで地殻変動の計測のためにGPSを用いていた研究者や大学院生にL4の地球物理学的応用を普及させることを目指し，GPS生データからL4時系列を作成する方法を解説する．また観測対象の時間スケールによって擾乱信号の最適な取り出し方が異なるため，それらの解析手法の詳細について個々に説明する．
2. RINEXファイルからのL4の作成
GPSデータの解析は，衛星や局の位置，地球回転，潮汐，中性大気遅延等の様々な要因をすべて考慮する必要があるため，概して複雑で個人レベルでのプログラム作成は難しい．一方L4では，これらの要因のすべてが(2)式で差を取る時に取り除かれているため，データ解析は本質的に簡単である．そのため，既存のGPS解析ソフトウェアでL4に相当する量を出力させて解析に用いることも可能だが，共通フォーマット（Receiver Independent Exchange format, or RINEX format）のGPS生データのファイルを直接読んでL4の時系列を作るのが手っ取り早い．
第一著者のウェブページでRINEXファイルを入力してL4に変換するフォートランプログラムrdrnx.fを公開している(http://www.ep.sci.hokudai.ac.jp/~heki/software.htm)．そのプログラムでは，最初にRINEXファイルのヘッダー（# / TYPES OF OBSERVの行）を読みこみ，様々
な種類のデータの順番を確認する．例えばオンライン公開されているGEONETのRINEXデータにはL1/L2の位相である“L1，L2”, Pコードを示す“P2”やC/Aコードの“C1”，また昨今はL1/L2の信号対雑音比を表す“S1/S2”なども含まれる．次にデータ本体部分をエポック毎（GEONETでは30秒毎）に読み込み，様々な衛星で得られたL1とL2の位相を取り出す．位相にそれぞれの波長を掛けて単位をラジアンから長さに変換し，(2)式の差をとってL4を求め，衛星毎にソートした時系列を出力する．本プログラムではL4 (単位m) に(3)式のようなファクターを掛けて，視線方向の電子数を積分した全電子数Total Electron Content (TEC，しばしばテックと発音される)に換算してある．
ΔTEC = (1/40.308) f12f22/(f12−f22) ΔL4 (3)
TECの単位はTECU (TEC unit，1 TECUは視線に沿った底面積1 m2の円柱に1016 個の電子が含まれることを意味する)を用いている．日本列島では，太陽静穏時の昼間のTECはおおむね10-20 TECU程度であり，夜間は数TECUに下がる．
位相データには一般的に整数値の不確定性があるため，実際にはL4の絶対値に意味はなく，衛星の観測開始から終了までの時間変化にのみ意味がある．(3)式でΔL4およびΔTECとなっているのはそのためである．上記のページには入力ファイル（RINEX形式）と出力ファイルの見本も置いてある．RINEXファイルに二周波のPコード情報（P1/P2）が含まれている場合は，それらの差を取ることによって整数値不確定性のないL4相当量（P4）を求めることもできるが，電離圏擾乱の大きさに比べてノイズが大きいため，ここで取り上げるような研究目的には実用的ではない．昨今ではサイクルスリップ（整数値不確定性のとび）の発生は稀であるが，このプログラム中ではサイクルスリップの検出も行っている．ただし見つけたスリップを修復はせず，単にサイクルスリップの無い最も長いデータの部分を取り出して出力する．
Figure １にはこうして得られた鹿児島県屋久島のGPS点での2009年7月22日のTEC変化を衛星毎に曲線で示している（衛星毎の時間変化にのみ意味があり，縦軸の絶対値には意味がない）．左半分でU字型の変化が顕著であるが，これは主に昇って沈む衛星の仰角変化に伴って視線ベクトルが貫通する電離圏の厚さが変化するためである（斜めに視線が貫く場合のTECをSlant TECと呼ぶ，詳細は3.2章参照）．U字型の変化の大きい部分が電子密度の高い昼間，小さい部分が電子密度の低い夜間に相当する．この日は皆既日食が起こっており，UTで2時頃を中心にTECが少し下がっているのが，生のTEC変化でも見ることができる（これについては後に詳しく議論する）．
多くの電離圏擾乱は比較的小さいため，衛星の仰角変化やTECの日変化等の長い時間スケールの変化と何らかの方法で区別する必要がある．地震や火山による擾乱のような数分程度の短い時間スケールの変動では，軌道情報を必要としない単純なハイパスフィルターが便利である（3.1章参照）．一方地震の前兆としてのTEC日変化の振幅変化を問題とする
ような場合は，TECの絶対値を最小自乗法で推定する必要が生じる（3.5章参照）．この場合は視線が電離圏を貫く角度の情報を用いるため，地球固定座標系における衛星の位置を知る必要がある．擾乱の信号を取り出すのに適切な手法は時間スケールに応じてかなり異なる．次の章では，対象とする様々な現象について解析手法を解説する．
3.様々な擾乱信号とその抽出
3.1.地震と火山噴火：特定の衛星に生じる短期的振動
火山噴火に伴う電離圏擾乱は，爆発的噴火で発生した音波が熱圏に達してF領域の電離圏に電子密度の濃淡を生じさせたものである．Heki (2006)は2004年9月の浅間山噴火の際に生じた電離圏擾乱から火山爆発のエネルギーを推定している．その後Dautermann et al. (2009)は，西インド諸島Montserrat島火山の2003年の爆発的噴火に関して同様の研究を行っている．2009年十月に発生した桜島火山の爆発的噴火に伴うTEC変化の例をFigure 2 に示す（Figure 1と同様時間変化のみに意味がある）．Peak-to-peakで最大0.2-0.3 TECU程度の擾乱が熱圏の音速で南に伝搬してゆく様子がわかる．擾乱は音波が熱圏に達するのに要する10分程度経過してから現れ，2分弱の周期を持つ．
地震に伴う電離圏の擾乱(CID, Coseismic Ionospheric Disturbance)(e.g. Heki and Ping, 2005)には，表面波（レーリー波）による地表の上下運動が作った音波が熱圏に達してもたらす擾乱と，震源近くの地表の地殻上下変動が作った音波が直接熱圏に達して起こす擾乱の二種類がある．前者は表面波の伝搬速度(4 km/sec程度)で伝わるため，大きな地震では熱圏での音速(1 km/sec以下)で伝わる後者と明瞭に区別できる(Astafyeva et al., 2009)． Figure 3は2008年6月の岩手宮城内陸地震によって生じた電離圏擾乱である．伝搬速度が遅いので，表面波起源のものではなく直達音波による擾乱であることが分かる．地震や津波が励起した大気の内部重力波による電離圏の擾乱は，我が国で発生した地震に関してはまだ検出例はない．
Heki et al. (2006) はスマトラ・アンダマン地震に伴う電離圏擾乱を解析して，断層すべり量や破壊伝搬速度を拘束する試みを行っている．またAstafyeva and Heki (2009) は，直達音波による電離圏擾乱の初動の増加と減少の分布が，地震波のP波の初動の押し引き分布のように地震の発生機構を反映していることを見出した．Figure 3の2008年岩手宮城内陸地震の例は正のTEC変化で始まっているが，これは逆断層地震の特徴である．直達音波による擾乱は距離に伴う幾何減衰が大きく，この擾乱も震源から300 km程度で見えなくなっている．擾乱の周期は約4.5分で火山噴火による擾乱に比べて長い．これは大気の共鳴周波数(e.g. Nishida et al., 2000) の中で最も周期の長いものに一致している．
地震，火山噴火いずれの例でも擾乱の周期は十分短いため，Figure 1でみるような衛星の仰角や太陽天頂角の変化に伴う長周期の変化と混同する危険性は少ない．この場合は，生のTEC変化を多項式近似したものを引き去るという簡便なハイパスフィルターで効率よく
信号を取り出すことができる．Figure 2と3の例では，それぞれ4次および5次の多項式を用いている．多項式の次数はデータを見ながら最適なものを選ぶことになるが，一般に時間窓を長くとるほど高次の多項式による近似が必要となる．なお対象とする擾乱の周期が既知の場合はウェーブレット変換も有効である(Heki and Ping, 2005)．また，地震や火山噴火に伴う電離圏の擾乱が南に選択的に伝搬してゆくのは，超高層での電子の動きが地磁気の方向に拘束されるため，北向きの音波が選択的に減衰されることを反映している(Heki and Ping, 2005)．
3.2. 太陽フレアに伴うTEC変化：全衛星に現れる短期的上昇
太陽面での爆発に伴って増加した太陽放射に対する電離圏の応答として，太陽に照らされた半球でTECが急増することをSITEC (Sudden Increase of TEC) と呼ぶ．Figure 4 に，2006年12月5日と13日に生じた二つのXクラス（静止衛星で監視している太陽からのX線の放射量に基づいて決められるフレアの規模，地震のマグニチュードのようなもの）の太陽フレアについて，フレア発生を含む約１時間のTEC時系列を示す。ここでは高度300 kmの層を貫くときの角度を用いて，鉛直方向に貫いた場合のTEC（斜めに貫いた場合のSlant TECと区別するためにVertical TEC, VTECと呼ばれる）に換算してある．
地震や火山噴火による擾乱では，たまたま視線が震源や火山の上空南側の電離圏を貫くGPS衛星のデータでのみ擾乱信号が見られる．また視線は波面（地表から伝搬してきた音波がつくる面状の電子の濃淡）に浅い角度で交わる必要がある．一方太陽フレアの場合は，昼側全体が同時に影響をうけるため，すべてのGPS衛星でTECが同時刻に増加するという顕著な特徴がある．また太陽フレアによるTEC上昇の時間スケールは数分程度と短いため，地震や火山噴火に伴う擾乱の解析と同じ手法（多項式近似によるハイパスフィルター）が使用可能である．
一般にフレアに伴って硬X線（波長が1オングストロームより短いX線）と極端紫外線(EUV)の放射量が数分の時間スケールで増え(impulsive component)，その後に軟X線を中心とする放射が数十分継続する(slow component)． 12月5日のフレアには12月13日のフレアに見られるような急峻な立ち上がりがみられない．これは，このフレアが太陽の縁辺部で発生したため，impulsiveな成分が太陽大気に吸収されて，軟X線を主成分とするslowな成分に対する応答が選択的に強調されたからである．
Impulsiveな成分による電子密度増加は，通常の電子密度の高度分布の相似性を保ちつつ全高度にわたって生じるが，slow成分による電子密度増加は高度100 km内外のE領域に選択的に生じる(Donnelly, 1976)．日置(2006)は，2003年10月に生じた大規模太陽フレアに関して，昼夜境界線の夜側で境界線から様々な距離にあるGPS局におけるSITECの量を比較することによって，フレアによる電子密度増加の高度分布を求めようとした．境界線から遠い局ほど高い高度の電子密度上昇を反映するという幾何学的な関係を利用したのである．しかし地球大気による吸収のモデル化が難しく精度よく高度分布を求めることができなか
った．今後は，GPS受信機を搭載した低高度衛星によるGPS掩蔽観測(e.g. Kursinski et al., 1997) を併用し，地上のGPS受信機ではわからない電子密度変化の高度分布を求める必要があるだろう．
3.3. ロケット（ミサイル）打上に伴うTEC変化：急激な減少とゆっくりした回復
ロケットや弾道ミサイルの排気ガスに含まれている水分子，水素分子，二酸化炭素等は，電離圏F領域の酸素イオンから正電荷を奪い，電子と解離性結合することによって電子の枯渇をもたらす(Mendillo et al., 1975)．こうしてできる電離圏の「穴」は宇宙開発の初期段階から知られていたが， HII-Aロケットの打ち上げに伴って生じた穴をGEONETが捉えたことから測地学研究者にも知られるようになった(Furuya and Heki, 2008)．また2009年4月5日には，北朝鮮から東に向かって発射された，大陸間弾道ミサイル（北朝鮮政府発表では衛星の打ち上げ）と考えられるテポドン2号の航跡がGEONETによって捉えられ，稠密GPS連続観測網の新たな応用として注目された(Ozeki and Heki, 2010)．
この現象の特徴は，TECの急激な減少（1-2分で減少が完了）と数十分かけたゆるやかな回復である．この回復過程はこれまでに述べた地震や火山噴火，太陽フレアなどによる擾乱に比べて時間スケールが長く，多項式近似によるハイパスフィルターでは回復過程をうまく取り出すことができない．以下に我々がテポドンによるTEC変化の解析に用いた手法を紹介する．まず地球に固定された座標系における衛星の位置を別途求めておく（詳細は4章を参照）．次に時々刻々の視線ベクトルが電離圏を貫く天頂角ζを計算する（電離圏高度は265 kmとしている，詳細はOzeki and Heki (2010)参照）．TEC (Slant TEC)は時刻tとζの関数であるが，それらは鉛直方向のTEC (VTEC)と衛星毎に異なるバイアスdを用いて
Slant TEC(t, ζ) = VTEC(t)/cosζ + d, (4)
と表せる（電離圏の空間構造は無視している）．VTECには仰角変化による変動分がない分，時間変化がSlant TECより単純でモデル化しやすい．実際のVTECは日変化するが，数時間の間の変化であれば時間の二次関数として，
VTEC(t)=at2 + bt + c. (5)
のようにモデル化できる．(5)を(4)に代入したものを観測方程式とし，a, b, c, dを最小二乗法で推定すれば，3.1.や3.2.で用いた単なる時間の多項式より，少ないパラメータで近似度の高いモデルができる．Figure 5で濃いなめらかな曲線で示したモデルはこのようにして求められたもので，TECの約二時間にわたる時間変化を忠実にモデル化している．ただし単一の衛星から推定したa, b, cから(5)式で復元したVTECの値は，パラメータ間の相関が高い場合あまり信頼できない．VTECそのものの絶対値を精度よく求める目的には，次の章で
述べるように複数衛星から共通の値を推定する必要があるが，特定の衛星で観測された擾乱を効率よく取り出すための手法としては本章の方法で十分である．ミサイルによって形成された電子の穴をGPS衛星と受信局を結ぶ視線ベクトルが通過することによって生じるTEC減少は，このようにして推定したモデルと実際のTECの差として取り出すことができる(Figure 5の下部の黒い曲線)．なお，類似の現象は宇宙から地球熱圏に突入する彗星（主成分は水の氷）によっても生じると思われる．小彗星仮説(e.g. Frank and Sigwarth, 1993)によると日本列島上空にも20トン規模の彗星が連日複数個突入していることになるが，まだ明らかな検出例は知られていない．
3.4. TEC絶対値の推定：皆既日食によるTEC減少の例
太陽フレアと対照的に，日食は月の影によって太陽の放射が遮られるため，昼間の半球上で月の影が広がる広範な地域でTECの減少が見られる．皆既日食の起こる地点のTECは，月の半影に入った時からゆっくり減少し，皆既日食時に最小となった後に徐々に回復する．このような数時間かけた現象は，TECの日変化や衛星仰角の変化に伴うゆっくりした変化と分離できない．これらを解析するには，仰角の異なる複数の衛星のTEC変化から，TECの絶対値を推定する必要がある．
前章の(4)式は，ある程度長い時間同一衛星を観測して，視線が電離圏となす角度ζが大きく変化しないと，バイアスdとVTECがパラメータとして分離しにくいことを示唆する．そのため個々の衛星から別々に得られたTECは一般に大きくばらつくが，すべての衛星に共通の値を推定することによって安定な推定が可能になる．また東の空と西の空を飛ぶGPS衛星において，視線が電離圏を通過する点の経度と観測点の経度の差を，時刻差に変換する（同じ空間パターンのTEC分布の下を日本列島が西から東に動いてゆくと仮定している）と，さらなる解の安定化を図ることができる．
VTECの時間変化のモデルとして，前章の例のように対象とする時間幅が数時間と短い場合は(5)式のように時間の多項式が適当である．しかし数時間以上の時間スパンでTECの変化を調べる場合，1時間毎の値を推定して折れ線で時間変化を近似する方法が適している．Figure 6は，皆既日食のあった2009年7月22日とその前後1日の合計3日間のTEC絶対値を一時間毎の折れ線で推定したものである．通常の日変化に加えて，皆既日食を迎えたUT 2時を中心に部分日食が生じた前後数時間にわたって全電子数が減少していることがわかる．
3.5.地震前兆としてのTEC日変化の異常
最後にTEC日変化の異常を地震発生と関連付けた最近のいくつかの研究を紹介する．Liu et al. (2001) は，1999年に台湾で発生した集々地震の4-5日前に，TECの日変化の振幅が異常に小さくなったことを見出し，地震の前兆の可能性を示唆した．地震の数日前にTEC日変化の振幅が小さくなることは，その後2008年四川地震を始め複数の地震で報告されてい
る(Liu et al., 2009)． 磁気赤道において大気潮汐によって生じた東向きの電場と北向きの磁場に、電子が相互作用して鉛直上方に高度千km付近まで移送され(EXBドリフト)，それらが磁場に沿って電離圏F領域に降りてきて赤道の南北両側磁気緯度15度付近に電子密度の濃い地帯を形成するのが赤道異常である．これらの前兆を示した地震はいずれも赤道異常の直下に震源域を持っているが，Liu et al. (2009) は，地震に先立つ岩石破壊が東向きの電場を何らかの原因で弱めることによって，赤道異常の位置と強さが変わりTEC日変化の振幅に異常が生じると考えた．
2008年岩手宮城内陸地震は，日本列島内陸浅部で生じた地震としては比較的大きく，かつ太陽活動の最小期に発生している．我が国において，地震の前兆として台湾のようなTEC変化が見られるかを確認するには絶好の条件がそろっている．ここでは震源近傍の岩手県奥州市胆沢のGPS点において，地震を含む六週間の期間をとり，前の章で述べた手法でTEC絶対値を1時間毎の折れ線で求めた（Figure 7）．一般にRINEXファイルはUTで一日毎に別ファイルになっており，VTECの絶対値の推定も一日毎に行われる．しかしそのようにして求めた時系列を複数日ならべると，零時で不自然な不連続がしばしば見られる．一方日の境界をまたいでも(4)式にある衛星毎のバイアスdは変わらない．したがって，零時のVTEC不連続を防ぐために，複数日のデータを結合し，日をまたぐ衛星のTEC時系列からは一つのバイアス値を推定する必要がある．ここでは，TEC変化を求めたい日の前後の日のデータを結合して三日間のTEC変化を一度に推定し，その真ん中の日のデータだけを取り出した．このようにして求めたデータを並べてゆけば，Figure 7に見られるように日の境界でも自然な連続した変化を見せる．
異常の判定はLiu et al. (2009)に倣い，過去15日間のmedian (中央値)を取り，median とquartileの差の1.5倍を自然な揺らぎと仮定してそれを超えるものを異常値とした．地震の一日後に見える大きな正の異常は，偶然に発生した磁気嵐によるもので地震とは無関係である．Figure 7をみると，Liu et al. (2009)が主張する地震4-5日前のTECの減少は負の異常として確かに見えている．しかしその異常の量はわずかであり，地震と関係ない時にもしばしば起きている程度の異常である．ここでは前兆の有無は議論しないが，条件が良いと思われる2008年岩手宮城内陸地震でもTEC日変化の振幅異常がこの程度の大きさであれば，地上GPSによるTECの監視を我が国で地震予知に役立てるのは難しそうである．
４．視線ベクトルと電離圏の交点
L4はgeometry-freeであり，局位置や衛星位置などの幾何学的なことは「ほぼ」無視できるが，例外的に幾何学的な計算が必要になる場合がある．それは，ある局で受信したあるGPS衛星の電波から求めたTECが，地理的にどのあたりの電離圏を見ているかを求める場合である．実際の電離圏は上下に幅を持っているが，便宜的に最も電子密度の高い高度（日本では300 km程だが，場合によっては電子密度の重心として，より高い400 km程度を仮定することもある）に薄い層を仮定し，衛星と受信機を結ぶ視線ベクトルがその層と交わ
る点(IPP, Ionospheric Penetration Point)，およびその地表への投影点(SIP, Sub-ionospheric Point)を求めて，「緯度何度，経度何度の電離圏を見ている」とすることが多い．
既に紹介したウェブサイト(http://www.ep.sci.hokudai.ac.jp/~heki/software.htm)では，RINEX形式の軌道情報ファイルを読み，そこに含まれているケプラー要素から，諸衛星の地球固定系における位置を三分ごとに計算するプログラムrdeph.fが提供されている．そこから，特定の時刻の衛星位置を内挿して求めれば，あとは局位置と衛星位置を結ぶベクトルが特定の高度の層を貫く点の座標を初等的な数学で計算することができる．Figure 2, 3のSIPの座標はこのようにして求めたものである．なおSIPはその定義から本質的に高精度で求める必要がない．様々な高度に分布した電子を高度300 kmで代表させた点ですでに大幅な近似だからである．従って軌道精度も放送暦以上のものを追求する必要はない．
５．おわりに
測地学の専門知識を持つ教員が所属する日本の大学では，GPSを用いた研究を大学院の学位論文や学部の卒業研究のテーマとすることが多い．それらのテーマは，最近発生した地震に伴う地殻変動に関する研究であったり，地震間の地殻変動の解析から地下における断層固着部分を推定する研究であったりと多彩である．また大気遅延から可降水量を求めて，その時間変化と気象現象の関連を探るGPS気象学的な分野にも様々な興味深い研究テーマがある(小司他, 2009)．しかし，それらの研究を自分で作成した計算機プログラムによって行うことは困難であり，既存のGPSデータ解析ソフトウェアをブラックボックスとして利用することが多い．
L4の解析は本質的に簡単で，大学院生や卒論生が気軽にプログラムを作成してRINEXファイルからそれらの情報を抽出することができ，努力に対する教育効果が高い．また本文中で紹介したプログラムを使用すれば，すぐに解析に取り組むこともできる．データの扱いが簡単であるにも関わらず，本論文で紹介したようにL4の解析による研究対象はバラエティに富み，未開拓の領域もまだ多いと考えられる．今後は様々な大学や研究所でGPS-TEC法による地球物理学研究が普及し，新たな発展を見せることを願っている．
謝辞
匿名の査読者に多くの適切なコメントをいただいた．
参考文献
Astafyeva, E. and K. Heki (2009): Dependence of waveform of neear-field coseismic ionospheric disturbance on focal mechanisms, Earth Planets Space, 61, 939-943.
Astafyeva, E., K. Heki, V. Kiryushkin, E. Afraimovich, and S. Shalimov (2009): Two-mode long-distance propagation of coseismic ionosphere disturbances, J. Geophys. Res., 114, A10307, doi:10.1029/2008JA013853.
Dautermann, T., E. Calais, and G. S. Mattioli (2009): Global Positioning System detection and energy estimation of the ionospheric wave caused by the 13 July 2003 explosion of the Soufrière Hills Volcano, Montserrat, J. Geophys. Res., 114, B02202, doi:10.1029/2008JB005722.
Donnelly, R.F. (1976): Empirical models of solar flare X ray and EUV emission for use in studying their E and F region effects, J. Geophys. Res., 81, 4745-4753.
Frank, L. A. and J. B. Sigwarth (1993): Atmospheric holes and small comets, Rev. Geophys., 31, 1-28.
Furuya, T. and K. Heki (2008): Ionospheric hole behind an ascending rocket observed with a dense GPS array, Earth Planets Space, 60, 235-239.
日置幸介(2006): 太陽フレアに伴う電離層全電子数上昇のGPS 観測, 測地学会誌, 52, 319-328.
Heki, K. (2006): Explosion energy of the 2004 eruption of the Asama Volcano, Central Japan, inferred from ionospheric disturbances, Geophys. Res. Lett., 33, L14303, doi:10.1029/2006GL026249.
Heki, K. and J.-S. Ping (2005): Directivity and apparent velocity of the coseismic ionospheric disturbances observed with a dense GPS array, Earth Planet. Sci. Lett., 236, 845-855.
Heki, K., Y. Otsuka, N. Choosakul, N. Hemmakorn, T. Komolmis, and T. Maruyama (2006): Detection of ruptures of Andaman fault segments in the 2004 Great Sumatra Earthquake with coseismic ionospheric disturbances, J. Geophys. Res., 111, B09313, doi:10.1029/2005JB004202.
Kursinski, E. R., G.A. Hajj, J.T. Schofield, R. P. Linfield, and K. R. Hardy (1997): Observing Earth’s atmosphere with radio occulation measurements using the Global Positioning System, J. Geophys. Res. 102, 23429-23466.
Liu, J. Y., Y. I. Chen, C. H. Chen, C. Y. Liu, C. Y. Chen, M. Nishihashi, J. Z. Li, Y. Q. Xia, K. I. Oyama, K. Hattori, and C. H. Lin (2009): Seismoionospheric GPS total electron content anomalies observed before the 12 May 2008 Mw7.9 Wenchuan earthquake, J. Geophys. Res., 114, A04320, doi:10.1029/2008JA013698.
Liu, J. Y., Y. I. Chen, Y. J. Chuo, and H. F. Tsai (2001): Variations of ionospheric total electron content during the Chi-Chi earthquake, Geophys. Res. Lett.,28, 1383-1386.
Mendillo, M., G. S. Hawkins, and J. A. Klobuchar (1975): A sudden vanishing of the ionospheric F region due to the launch of Skylab, J. Geophys. Res., 80, 2217-2225.
Nishida, K., N. Kobayashi, and Y. Fukao (2000): Resonant oscillations between the solid earth and the atmosphere, Science, 287, 2244-2246.
Ozeki, M. and K. Heki (2010): Ionospheric holes made by ballistic missiles from North Korea detected with a Japanese dense GPS array, J. Geophys. Res., A09314, doi:10.1029/2010JA015531.
小司禎教・岩淵哲也・畑中雄樹・瀬古弘・市川隆一・大谷竜・萬納寺信崇(2009): GPS気象学：GPS水蒸気情報システムの構築と気象学・測地学・水文学への応用に関する研究，
測地学会誌，55, 17-38.
Figures
Figure 1. Time series of raw slant total electron content (TEC) at a station in the Yaku Island (950493), southern Japan, on 22 July, 2009. U-shaped changes reflect variations of GPS satellite elevation as it moves in the sky, i.e. line-of-sight vector penetrates thicker ionosphere when satellite elevation is lower. A small decrease of TEC due to a total solar eclipse is recognized at ~2:00 UT.
Figure 2. Time series of slant TEC changes at nine GPS stations in islands south of Kyushu, southern Japan, with GPS Satellite No.11 on 2 Oct, 2009. The Sakurajima volcano (star in (a)) exploded at ~7.75 UT, and TEC pulsation appeared ~10 minutes after the explosion and propagated southward by ~0.7 km/sec (b). The vertical axis shown at the right-hand side of (b) shows the distances of sub-ionospheric points (SIP) (a) from the volcano (black star). Open and solid squares in (a) shows GPS stations and their SIP at 8:00 UT.
Figure 3. Time series of slant TEC changes at ten GPS stations in Honshu, central Japan, with GPS Satellite No.8, from 8:30 to 9:30 JST, 14 Jun, 2008 (from 23.5 UT 13 Jun. to 0.5 UT 14 Jun.). The 2008 Iwate-Miyagi Nairiku earthquake (M7.2) occurred at ~23.43 UT (a), and TEC oscillation appeared ~10 minutes after the earthquake and propagated southward by ~0.9 km/sec (b). The vertical axis shown at the right-hand side of (b) shows the distances of SIP at ~23.9 UT (a) from the epicenter (black star).
Figure 4. Time series of vertical TEC changes at the Sutherland, South Africa, GPS station (a) and the Chichijima, Japan, GPS station (b), on December 5 and 13, 2006, respectively. Solar flares of X9.0 and X3.4 occurred at ~10.5 UT and ~2.5 UT on these days, respectively. Sudden increases of TEC are seen for all the satellites simultaneously. The December 5 X9.0 flare occurred at the edge of the solar disk, for which we hardly see the impulsive component. The December 13 flare occurred near the center of the solar disk, and impulsive components are well recognized.
Figure 5. Time series of slant TEC changes at three GPS stations (950232, 950198, and 950203) in Honshu, central Japan, on April 5, 2009. Raw TEC time series (light gray curves) were modeled with vertical TEC changing as a quadratic function of time (black smooth curves, see equations 4 and 5). Anomalous TEC decrease (black curves) are defined as the difference between the model and the observed slant TEC. Sudden dips of TEC are caused by ionospheric hole formed by the exhaust plume of the Taepodong-2 ballistic missile from North Korea (Ozeki and Heki, 2010).
Figure 6. Time series of absolute vertical TEC at a GPS station (950493) in Yaku island, southern Japan, over a three-days period July 21-23 (day of the year 202-204), 2009, estimated every hour using all the available GPS satellites together with satellite-specific biases according to the equation (4). In addition to daily TEC changes, we see TEC decrease around ~02:00 UT (11:00 in the morning in local time) on the middle day, when the island experienced a total solar eclipse.
Figure 7. Time series of absolute vertical TEC (open circles connected with black lines) at the Isawa GPS station (970796) in the Iwate Prefecture, northeastern Honshu, over six weeks including the 2008 July 13 Iwate-Miyagi earthquake (day of the year 165 in UT, thick vertical line), estimated every hour with the same method as Figure 6. Thick black curve shows the median of the preceding 15 days with upper and lower bounds of natural variability (taken 1.5 times as far from median as quartiles) shown by thinner curves. Gray and black shades at the bottom show the amount of positive and negative anomalies (amount above/below the upper/lower bounds of natural variability). A large positive anomaly about one day after the earthquake is related to a geomagnetic storm. We can see negative anomalies 4-5 days before the earthquake. However, they are faint, and similar amount of anomalies occur regardless of earthquakes.
この記事を読んだ人はこんな記事も読んでいます（表示まで20秒程度時間がかかります。）
スパムメールの中から見つけ出すためにメールのタイトルには必ず「阿修羅さんへ」と記述してください。すべてのページの引用、転載、リンクを許可します。確認メールは不要です。引用元リンクを表示してください。
