著者
山口 龍彦 萩野 恭子 平 陽介 濱田 洋平 斎藤 仁志 小野寺 丈尚太郎 板木 拓也 星野 辰彦 稲垣 史生
出版者
日本地球惑星科学連合
雑誌
日本地球惑星科学連合2018年大会
巻号頁・発行日
2018-03-14

IODPを始めとした掘削科学のプロジェクトにおいて、珪藻、有孔虫、放散虫、石灰質ナノプランクトンなどの微化石を用いて堆積物の地質年代を正確に決定することは非常に重要である。微化石の抽出と分類群の正確な同定には豊富な経験・知識が必要とされ、同定と堆積物の年代決定には時間を要する。DSDPが開始された1960年代より微化石の専門家は、掘削船に乗船し船上で微化石の分析を行い、掘削科学に貢献してきた。一方で、コンピュータ処理能力の向上により運転など従来人間が担ってきた操作が人工知能(AI)に置き換わられつつある。人工知能の構築に広く用いられているのはディープラーンニングと呼ばれる機械学習技術であり、多層のニューラルネットワークを用いることで、データの特徴を段階的により深く学習させる手法である。近年、国際空港の出入国管理、アミューズメントパークやコンサート会場の入場セキュリティを担う顔認証システムにおいて、ディープラーニングを用いることで、その識別精度を向上させる実用例も報告されている。我々は、この機械学習・画像認識技術の微化石年代測定への応用、さらには自動年代決定AIシステムの開発について研究を一昨年から開始した。昨年の本学会では、比較的単純な2種類の石灰質ナノ化石および珪藻化石について、NECのAIソフト「RAPID機械学習」による自動識別が可能であることを示した。今回の発表では、実用化に向けたさらなる検討として、分類する石灰質ナノ化石の種類を第四紀の堆積物から多産するEmiliania huxleyi, Gephyrocapsa oceanica, Gephyrocapsa ericsoni, Reticulofenestra haqii, Gephyrocapsa caribbeanica, Reticulofenestra productの6タクサ(種)に増やし、自動分類を試みた。ランダムに30枚の教師画像を与えて構築した分類モデルを用いて教師画像とは別の120枚(6種×20枚)の標本の画像の自動識別を行った。この結果、G. oceanicaやR. haqiiの判別的中率は100%、95%であり、自動判別の有効性が示された。一方、その他の4種類の判別的中率は0-35%であった。判別的中率が低い分類群では同定の鍵となるbridgeの認識がうまく行われていないと考えられる。偏向板を回転させたり、画像の解像度を上げた教師画像のみを与えることでより高精度の分類が可能になることが予想される。本発表では、石灰質ナノ化石の他の微化石への応用例も紹介し、自動分類の可能性について議論を深めたい。
著者
長岡 央 Fagan Timothy 鹿山 雅裕 長谷部 信行
出版者
日本地球惑星科学連合
雑誌
日本地球惑星科学連合2016年大会
巻号頁・発行日
2016-03-10

Based on previous study of lunar returned samples and meteorites, the main suites of pristine nonmare igneous rocks have been classified into the following four types: (1) ferroan anorthosite (FAN) or ferroan anorthositic-suite (FAS), (2) magnesian suite (Mg-suites), (3) alkali-anorthosite-suite and (4) KREEP basalt and possibly related rocks such as quartz-monzogabbro (QMG) /monzodiorite (QMD), granite and felsite. The latest suite type, the evolved rock samples related to KREEP, may have been derived from residue of the lunar magma ocean (urKREEP), or from low degrees of partial melting or some other process to account for their high incompatible trace element (high-ITE) compositions. Granite and felsite have Th-rich compositions (10 to 60 ppm), and such lunar samples with bulk SiO2 content of >60wt% originated from silicic volcanic or exposed intrusive material. Recent global remote sensing data have presented several candidates of silicic volcanism over the Moon based on indicators such as ITE-rich compositions, dome-like topography, characteristic infrared spectra (Christiansen Feature), and high albedo. Silica-rich, broadly granitic samples have been identified in lunar returned samples and lunar meteorites, but are rare.Lunar meteorite Northwest Africa (NWA) 2727 is a breccia paired with NWA 773 and the other meteorites of the NWA 773 clan. An olivine cumulate gabbro (OC) is common to most of these lunar meteorites within the NWA 773 clan; in fact NWA 2977 and 6950 consist entirely of OC lithology. However, in addition to the OC lithology, several clast types, including in olivine phyric basalt, pyroxene phyric basalt, pyroxene gabbro, ferroan symplectite, and alkali-rich-phase ferroan (ARFe) rocks have been discovered from the NWA 773 clan. The ARFe clasts have K-feldspar and/or felsic glass, a silica phase and minerals rich in incompatible elements such as merrillite. In this work, we characterize a felsic clast in NWA 2727 and compare our results with other lunar samples to discuss silicic volcanism.A polished thin section (PTS) of NWA 2727 was investigated by a combination of petrographic microscopy and electron probe micro-analysis. The NWA 2727 breccia includes a variety of large-scaled lithic clasts (>1mm) including: OC, ferrogabbro, pyroxene-phyric basalt, and the felsic igneous clast. The felsic clast has a modal composition of 37% silica, 34% plagioclase, 14% K-feldspar, 6% high-Ca pyroxene, 5% fayalite, 3% Ca-phosphate, 1% ilmenite, and traces of troilite and chromite. Feldspar compositions of the plagioclase are near An85-90. Two compositional types of pyroxene were identified—one near hedenbergite (Wo46Fs53, Mg#=1 [calculating Mg# as Mg/(Mg+Fe)x100]) and the other with zoning and more magnesian compositions (Wo25-30Fs55-65, Mg#=8~20). The K-feldspar is also zoned with variable concentrations of Ba, clearly detected in elemental X-ray maps (quantitative analyses of Ba are planned). The abundance of silica + feldspars (>80 mode%), the high proportion of K-feldspar to plagioclase, and the very ferroan compositions of mafic minerals attest to the felsic composition of this clast. Subhedral-euhedral olivine crystals up to 0.3 mm in maximum length are preserved, and silica and K- and Ba-feldspar occur in elongate parallel crystals indicating an igneous origin. These observations indicate that this clast was derived from silica-rich magma.Silicic volcanism is also interesting from the viewpoint of landing site candidates for future lunar landing mission. Global gamma-ray observations have presented several high-Th regions in PKT, but the main lithology of the Th-rich regions remains a subject of dispute; possibilities include mafic impact-melt breccia, KREEP basalt, QMD, and felsite/granite. If a lander/rover mission to a high-Th region is equipped for analysis of major elements, in situ analyses on the Moon can be compared with silica-rich samples such as the felsic clast in NWA 2727.
著者
髙濱 聡 溜渕 功史 森脇 健 秋山 加奈 廣田 伸之 山田 尚幸 中村 雅基 橋本 徹夫
出版者
日本地球惑星科学連合
雑誌
日本地球惑星科学連合2016年大会
巻号頁・発行日
2016-03-10

気象庁では,地震調査研究推進本部の施策に基づき,全国の高感度地震計のデータを収集し震源決定等の処理を一元的に行い,その結果を地震カタログとして公表している。 現在の地震カタログは,精査により一定の基準を満たしたものを掲載することとしている。しかし,東北地方太平洋沖地震後の余震域では余震活動は低下してきているものの以前と比べれば活発な状況にあり,処理対象地震の規模の下限を上げた処理を行っていることから,検知されても処理基準未満であるため地震カタログに掲載されない地震がある。 これに対処するため,平成25年度に同本部地震調査委員会の下で検討が行われ,1)これまでの検知能力は維持し,2)検知された地震のすべてを地震カタログへ掲載する,3)精度に段階をつけた品質管理を行う,の3つの方向性を示した報告がまとめられた。 気象庁ではこの報告を踏まえ,自動震源を活用するなど,震源決定処理手順を変更し改善する。具体的には,領域と深さごとに精査を行う地震のMの閾値(以下,Mthと記す)を設定し, Mth以上の地震については,現行通りに精査した震源決定を行い,Mth未満の地震については自動震源を基本とし,検知されても自動震源が求まらない地震については,最大10点程度の観測点を検測する簡易な手順により震源決定を行うことで、処理の効率を高める。精査される震源の目安は、内陸の浅い地震はM2以上とし,海域については陸域(観測網)からの距離に応じてMを上げて最大でM4以上とする。また,処理方法と精度の違いがわかるような登録フラグを新たに設ける。 ここでは、新たな地震カタログを用いて気象庁が作成する震央分布図等の資料について、具体例を紹介する。
著者
高田 亮
出版者
日本地球惑星科学連合
雑誌
JpGU-AGU Joint Meeting 2017
巻号頁・発行日
2017-03-10

カルデラ形成を伴う大規模噴火,またはそれに準ずるVEIで6以上の噴火は,既存のDB(Seibert et al.,2010)によれば,100年に1-2回程度発生している.しかし,平均値とは異なり,インドネシアだけで最近200年間で3回のカルデラ噴火を経験している.本発表では,その中で,比較的最近の情報が多く残るKrakatau1883のカルデラ噴火について,前兆現象を文献(Nishimura, 1980, Yokoyama, 1981, Simkin and Fiske, 1983;Carey et al.,1996; Mandeville et al., 1996)などでレビュウし(Takada, 2010; Takada et al., 2012),また,現地での巡検視察(1996,2006)などをもとに考察する.Krakatau火山は,ジャワ版「列王記」で,西暦416年の噴火記録があるが詳細は定かでない.1596年のスケッチでは,1883年噴火直前の火山島とほぼ同じ形が描かれていた.1680年には爆発的噴火が起こったことが銅版画で残されている,18世紀以後は,1883年の噴火まで噴火記録はなく,樹木に覆われた島は農地や硫黄採掘,温泉などとして利用されていた.クライマックスである1883年8月27日に先立つ3ヶ月半以上前の5月より有感地震が起こった.1883年噴火の前兆現象として,広範囲に噴気孔が広がっていく異常現象がとらえられていた.5月20日には噴煙が約11km上がった.爆発音,振動,少量の降灰が記録されている.5月22日には間歇的な灰噴火が起こり,海面には軽石が浮遊したまた,5月27-28日には,観光ツアーが島を訪れ,北側の噴火口の様子を報告した.その後も,6月に入り役人やオランダの軍人が上陸し観察記録を残している.6月-7月は,2箇所から噴煙が上がっていたが,8月11-12日には,オランダの大尉が上陸すると,約5kmの範囲で3箇所の噴火口と14箇所の噴気孔を確認していた.クライマックスは,8月26日の早朝から始まり,14時に爆裂とともに噴煙中が26km上がったと記録されている.27日には,早朝から4回の噴煙柱が立ち上がり,最高で38kmに達したと言われている.それに伴う津波も発生し,被害を拡大した.火砕流や火砕サージも発生し,対岸の集落を襲った.この噴火で約12km3のマグマが噴出したと報告されている(Carrey et al., 1996).本研究は,原子力規制庁からの受託研究として行われた.
著者
風早 康平 高橋 浩 佐藤 努 高橋 正明 大和田 道子 安原 正也 森川 徳敏
出版者
日本地球惑星科学連合
雑誌
日本地球惑星科学連合2018年大会
巻号頁・発行日
2018-03-14

我々は2003年から阪神地域において,河川,水路等に現れる”赤水”の目視調査および水質調査を行っており,その分布の概要がわかってきたので報告する.”赤水”湧出の特徴として,比較的標高の高い六甲山麓部にも存在し,一部のものは石垣等から染み出している.これらは,ほぼすべて重炭酸型であり,pHが6-7の弱酸性のものが多くCO2が溶解していることを示唆する.”赤水”が赤褐色~オレンジ色を呈する理由は,湧き出し口でリモナイト(FeOOH)が沈殿することによる.この鉄は,地下水中にCO2があった場合に,岩石・鉱物が風化することにより地下水に溶解する.そして,浅所で酸化することにより赤褐色からオレンジ色を生じる.阪神地域では,地下から断層を通じて上昇する有馬型熱水が浅所でCO2を遊離し,気泡となったCO2が上昇し,浅層地下水を炭酸化すると考えられている.したがって,”赤水”の分布域は,有馬型熱水の上昇域に関連していると考えられる.”赤水”の分布が地下に伏在する断層と関連があるかどうか検討をおこなった.その結果,六甲,宝塚,須磨,仮屋沖(延長),芦屋,甲陽,西宮,伊丹の各断層と昆陽池陥没帯に関連づけられる”赤水”が存在していることがわかった.また,断層が存在していても赤水がほぼ見られない地域もあった.講演では,”赤水”調査の概要,産状を紹介し,詳細な分布・特徴について議論したい.
著者
大木 郁也
出版者
日本地球惑星科学連合
雑誌
日本地球惑星科学連合2016年大会
巻号頁・発行日
2016-03-10

カルデラ内イグニンブライトは,カルデラ外に比べ,カルデラ形成に伴う噴火の始まりや終わり,その後の情報をより完全に記録している.しかしながら,厚いカルデラフィル堆積物により埋積されるため,カルデラ底部へアクセスできないなどの問題点がある.それゆえ,カルデラ内イグニンブライトの層序及び地質構造,堆積相を連続的に記載した例は少ない. 本研究は,カルデラ内イグニンブライトを連続的に観察できる理想的な例として三途川カルデラをとりあげ,そのカルデラ内イグニンブライト(虎毛山層)の層序や地質構造,堆積相を記載し,年代測定を行った.その結果,三途川カルデラの噴火史及び地質構造発達史について明らかにした.秋田県南部に位置する三途川カルデラは,カルデラ形成時に堆積したとされる軽石流堆積物(虎毛山層)により埋積される.虎毛山層は主に結晶に富むデイサイト質の火山礫凝灰岩,ブレッチャ,凝灰岩からなる.全体の層厚は1500m以上である.更新世(1.2 Ma;K-Ar法)に相当し,女川から西黒沢相当層間の基盤岩を不整合関係に被覆する.本層は5つの岩相からなる.5つの岩相とは,(1)ユータキシティック組織が発達した塊状無層理の火山礫凝灰岩(emLT), (2)塊状無層理の角礫岩(mlBr), (3)斜交層理が発達した火山礫凝灰岩(xsLT), (4)平行層理が発達した凝灰岩(//sT), (5)遍在的に成層構造を持つ火山礫凝灰岩(dsLT)である.emLTは本層の主体をなし繰り返し分布する.mlBr及びxsLTはemLTの下位に発達し,//sT及びdsLTはemLTの上位と中部にそれぞれ発達する.5つの岩相はシャープまたは漸移的に変化する.これらの岩相の特徴及び接触関係はすべてイグニンブライトの岩相を特徴づけるため,本層はイグニンブライトのシーケンスからなると判断できる.イグニンブライトに先行する降下軽石堆積物が欠如することやイグニンブライトの層厚が1500 mを超えること,本層の分布がカルデラ内に限られることから,本噴火は約1.2 Maにマグマ溜まりの天盤崩壊をトリガーとして発生し,プリニー式噴火フェーズのないイグニンブライトを形成するフェーズから始まったと考えられる.さらにイグニンブライト岩相が繰り返し分布することから,本層は7層の火砕流堆積物(PDC-1 to PDC-7)に識別される.7層の火砕流堆積物の存在より,本噴火で発生した大規模火砕流はwaxingとwaningを繰り返しながら少なくとも7回の火砕流パルスを発生させたと解釈できる.大規模火砕流の給源方向は,イグニンブライト中に発達するデューン構造やインブリケーションの方向から推定でき,北東から南西方向であると考えられる.給源方向にある大鳥谷沢では結晶に乏しいことやmlBrが多く狭在することからも支持される.また,虎毛山層は高松岳周辺を中心に環状に分布し,外側へ急傾斜する.この地質構造は高松岳を中心としたドーム状の隆起構造を示唆する.この隆起構造は,後カルデラ期の再生ドームと考えられ,虎毛山層形成後に発達したと考えられる.
著者
大木 郁也 大場 司
出版者
日本地球惑星科学連合
雑誌
日本地球惑星科学連合2014年大会
巻号頁・発行日
2014-04-07

秋田県南部に位置する三途川カルデラは,約1Ma以前に大規模火砕流を伴うカルデラ陥没により形成した.本地域には,カルデラ形成時に堆積したとされる虎毛山層が分布する.虎毛山層は,下位より虎毛山凝灰岩部層,皆瀬川凝灰岩部層からなる.虎毛山凝灰岩部層は,溶結凝灰岩,火山礫凝灰岩,凝灰質砂岩・黒色頁岩・礫岩の互層からなり,層厚は900mに達する.皆瀬川凝灰岩部層は,火山礫凝灰岩,凝灰岩,礫岩からなり,層厚は450mに達する.虎毛山層は,8層の火砕流堆積物(PDC-1~8),土石流堆積物(DF-1),湖成堆積物(LD-1)から構成される.これらの層序は,下位よりPDC-1,DF-1,LD-1,PDC-2~PDC-8からなる.各層の厚さはPDC-1が20m,DF-1が80m,LD-1が140m,PDC-2が50m,PDC-3が250m,PDC-4が200m,PDC-5が340m,PDC-6が160m,PDC-7が90m,PDC-8が30mである.火砕流堆積物は,塊状無層理の火山礫凝灰岩からなり,軽石と異質岩片を含む.しばしば,炭化木片を含み,脱ガスや柱状節理が発達する.PDC-4, 6は火砕流堆積物の基底部はグラウンドサージ堆積物からなる.このグラウンドサージ堆積物には,低角斜交層理が発達し,デューン構造が認められる.このうちPDC-6は,グラウンドサージ堆積物の下位にグラウンドブレッチャーが認められる.このグラウンドブレッチャーは,最大礫径2.5mの異質岩片を含む基質支持礫岩からなる.また,PDC-1, 3, 4, 8は特徴的に溶結相を伴う.溶結相には,ユータキシティック組織やスフェルライトが認められる.土石流堆積物(DF-1)は,層理が発達し,円礫を主体とする礫支持礫岩からなる.礫は平行に配向し,弱く逆級化する.湖成堆積物(LD-1)は,黒色頁岩及び凝灰質砂岩,礫岩の互層からなる.黒色頁岩にはラミナが発達し,凝灰質砂岩には葉理・層理が発達し,礫岩は塊状無層理である.湖成堆積物(LD-1)の上位のPDC-2は,水中環境での堆積を示唆する.8層の火砕流堆積物の存在は,本地域では火砕流が少なくとも8回発生していたことを示唆する.火砕流堆積物(PDC-4)の流向方向は,グラウンドサージ堆積物のデューン構造から,北東から南西方向であると推定でき,給源位置は滝ノ原火口であると推定した.休止期間を示す湖成堆積物(LD-1)が虎毛山層中部に狭在し,カルデラ陥没が少なくとも2回発生したと推定される.地層の走向は石神山周辺を中心とする半同心円構造をなし,その傾斜は半同心円の外側を向く.この構造は,石神山周辺を中心とするドーム状の隆起構造を示唆する.この隆起構造は,再生ドームであると考えられ,カルデラ中心域にあたる小安岳周辺の基盤岩の高まりの原因の一つである.再生ドームの形成と厚い火砕流堆積物の分布と環状割れ目の存在は,三途川カルデラがValles型カルデラである可能性を示唆している.
著者
深尾 良夫 三反畑 修 杉岡 裕子 伊藤 亜紀 塩原 肇 綿田 辰吾 佐竹 健治
出版者
日本地球惑星科学連合
雑誌
日本地球惑星科学連合2018年大会
巻号頁・発行日
2018-03-14

1: 津波地震は地震の規模と比較して異常に津波の規模が大きな地震(Kanamori, PEPI, 1972)で、典型的には海溝すぐ内側に起こる低角逆断層地震が津波地震になりやすい。一方、ほぼ10年に一度繰り返して起こる鳥島近海地震(Mw=5.6-5.7)は、海底カルデラ下で発生する全く別種の津波地震である。最新のイベント(2015年5月2日)であるが、その地震波と津波を震央距離約100km(北北東)の海底に展開した10点の水圧計アレーが記録した。本発表では、津波地震の観測からメカニズム提唱までの概要を報告する。津波解析の詳細は、三反畑らによる別途発表に譲る。2: 10点の絶対圧力計(7000mの水深に相当する圧力を9桁の分解能で測定可能)を最小辺長10km、最大辺長30kmの正三角形を構成するよう配置する。今回は、このアレーを青ヶ島東方沖水深1470-2240mの海底に2014年5月に設置し翌年5月に回収した。このアレーは、周期半日の内部潮汐波(Fukao et al., JpGU, 2017)、周期50-200秒の長周期海洋重力波(Tonegawa et al., JGR, 2018)、周期100-300秒の津波(Sandanbata et al., PAGEOPH, 2018; Fukao et al., Sci. Adv., 2018)など、多様な海洋現象の観測に有効である。3: アレーの観測点配置は、卓越周期200秒の津波がアレーを伝播する間に位相がほぼ1周期ずれ、そのズレを10点で測定することに相当し、高精度な位相解析が可能である。得られた位相速度は周期に依存し理論的な分散曲線と良く一致する。点波源をスミスカルデラ内に仮定し、周期に依存する局所位相速度分布図を用いて破線追跡を行うと、観測された走時と入射方向の周波数依存性をよく説明できる。仮に波源がリムにあるとすると、測定された到達時刻あるいは入射方向との一致は有意に悪くなる。一方、津波の初動は、波形の立ち上がり寸前のゼロ線を切る瞬間として読み取ることができるが、そのアレー通過速度は理論的な長波速度に一致する。この立ち上がりを地震発震時刻まで逆伝播させると、波源域の縁がカルデラリムにあることがわかる。津波波源はカルデラ全体にわたり、それを超えることなかったと推測される。一方、USGS、JMA、GCMTの求めた地震の震央はカルデラサイズを超えて散らばり、津波波源のほうが高精度で求められていることがわかる。4: 津波波源をモデリングするためにt=0の瞬間に海面擾乱を与え、その擾乱の伝播をブジネスク方程式の解として求めた。初期擾乱を軸対称とし、観測波形を最も良く説明する波源域の半径Rと中心隆起の振幅Aをグリッドサーチにより求めた。最適半径A=4kmはスミスカルデラのサイズとよく一致する。最適波高はA=1.5mと求められた。最適モデルに基づいて計算された波形と観測波形との一致は驚くほど良い。アレーで観測された最大振幅は約2cmであるが、八丈島(八重根港)の験潮儀には約60cmの最大振幅が観測されている。八重根港における津波波形を、湾の複雑な地形と津波の非線形効果を考慮して計算すると、計算波形は観測波形と驚くほどの良い一致を示した。5: この地震のメカニズムはT軸がほぼ鉛直のCLVD(Compensated Linear Vector Dipole)であり、Mwは5.7相当である(JMA, GCMT)。しかしこのメカニズムでは、直径8kmの波源域、1.5mの津波波高に相当する海底変位を生ずることはできない。震源を如何に浅くしようとも海底変位は津波の初期波高のたかだか1/20程度しかならない。両者のこの大きな差異は、実際のメカニズムがCLVDではなくHorizontal Tensile Fault(HTF)が鉛直に開口したためであったことを示唆する。このメカニズムが海底下の極浅部に働く場合、(1)遠方長周期地震波の励起効率はゼロに近く、一方で津波の励起効率は最大、(2)励起された地震波は殆ど上側(海側)に放出されるので、断層面上の変位は上側に集中し、水中音波の励起効率も大きく増幅される、(3)震源で体積変化なしと仮定して遠方変位場からメカニズムを求めるとモーメントの不当に小さなVertical-T CLVDが得られる。これら3つの極浅部HTFの特徴は鳥島近海津波地震の特徴と整合する。
著者
姫松 裕志 Freysteinn Sigmundsson 古屋 正人
出版者
日本地球惑星科学連合
雑誌
日本地球惑星科学連合2018年大会
巻号頁・発行日
2018-03-14

陸域におけるプレート発散境界を観測できる地域のひとつとしてアイスランドが挙げられる.プレート境界に沿って活動的な火山がみられ,いくつかの火山は氷河や氷帽などの雪氷に覆われている.これらの火山が噴火した際には火砕流や火山灰の降灰などの他に,山麓のインフラに被害をもたらすラハールなどの融雪型洪水を引き起こす可能性がある.近年のアイスランドでhあ1996年Gjalp火山や2010年Eyjafjallajokull火山の噴火に伴う洪水被害が報告されている.Bardarbunga火山もアイスランドの中で雪氷に覆われている火山のひとつであり,ヨーロッパ最大の氷帽であるVatnajokull氷帽の北西に位置している.2014年8月に始まったBardarbunga火山におけるダイク貫入イベントでは,氷帽から10km離れた地点で震源の移動は終息し,この地点で始まった割れ目噴火は2015年2月まで続いた.このイベントに伴う非雪氷域における地殻変動は衛星SARデータを利用した干渉SAR(InSAR)やピクセルオフセット法によって報告されている (Sigmundsson et al., 2014, Nature; Ruch et al., 2015, Nat. Comm.).これらの結果はリフト帯におけるダイクの貫入時に観測されるグラーベン構造を形成する地殻変動の描像(最大6mの沈降と最大2mの西北西-東南東方向の水平変位)を明らかにした. 航空機高度計による観測から氷帽のダイク経路上に直径500m,深さ20mを超える氷帽の円状の沈降(Ice cauldron)や最大60mのBardarbungaカルデラの沈降が報告されている.ダイク貫入イベント前後の数値標高モデルの差分(DEM difference)は氷帽下における地殻変動の影響を受けた鉛直変位を示した.本研究では衛星SARデータにピクセルオフセット法を適用することにより,氷帽上の変位からダイク見入イベントに伴う氷帽下の地殻変動の抽出を試みた.ダイク貫入イベントの前後に撮像された画像ペアにピクセルオフセット法を適用した結果,氷帽の定常的な表面変動を示すシグナルに加えて,地殻変動に伴う氷帽上の表面変動を示すシグナルを検出した.そこでダイク貫入イベント時のシグナルから任意の重みをつけたイベント前のシグナルを差し引くことにより,氷帽下における地殻変動を推定した.推定された氷帽下における地殻変動データを使用して,従来のダイク開口モデルから改良した新しいモデルについて議論する.
著者
小島 慎也 佐藤 香枝 前田 亮太 呉 宏堯 矢田 拓也 森田 敏明 岩崎 博之
出版者
日本地球惑星科学連合
雑誌
日本地球惑星科学連合2014年大会
巻号頁・発行日
2014-04-07

明星電気株式会社は、小型気象計POTEKA Sta.(ポテカ:Point Tenki Kansoku、以下POTEKA)を開発した。POTEKAは気温・湿度・気圧・感雨・日照を1分間隔で測定でき、従来気象計と比較して安価で、設置が容易なため稠密な設置及びデータ収集が可能である。そのPOTEKAを用いて、伊勢崎市内小中学校及び同市周辺のコンビニ(SAVE ON)に約1.5~4km間隔で計55ヶ所に設置した。本稿では、顕著な観測事例として8月11日に高崎市・前橋市で発生した突風現象の観測結果について紹介する。8月11日18時頃に高崎市から前橋市にかけて突風が発生し、住家の屋根の飛散などの被害がみられた。POTEKAの気温1分値を見ると、最大12分間で-13.9℃の気温低下がみられた。前橋地方気象台発表の突風経路に近いPOTEKAの海面補正した気圧の1分値時系列を下図に示す。気象台の10分値の気圧は徐々に増加していく傾向しか見られないが、POTEKAの1分値では、1~2hPa程度の一時的な上昇がみられた。これはダウンバースト発生時の下降流による一時的な気圧上昇であると示唆される。さらに詳しく見ると、気圧の上昇は2回発生している地点もあり、1回目はガストフロントによるもの、2回目はダウンバーストによる上昇と考えられる(詳細は「地上稠密観測POTEKAによるダウンバーストとガストフロントの識別」を参照のこと)。今回の稠密観測のようなダウンバースト・ガストフロント発生時の地上における気圧変化を、これほど細かい時間的・空間分解能で観測した事例はほとんど見られない。このような稠密観測をすることによって、突風の種類の判別や突風に対する事前の注意喚起が出来る可能性がある。謝辞:本プロジェクト始動にあたり、サンデン(株)殿、(株)セーブオン殿、伊勢崎市教育員会殿にはPOTEKA設置のご協力を頂きました。ここに御礼申し上げます。
著者
村上 亮 古屋 正人 高田 陽一郎 青木 陽介 小澤 拓 島田 政信
出版者
日本地球惑星科学連合
雑誌
日本地球惑星科学連合2018年大会
巻号頁・発行日
2018-03-14

はじめに 合成開口レーダ(SAR)は,全天候性,広域性,非接触性,計算機親和性など,地殻変動観測に適した多くの長所を有している.一般に,植生が広く分布するわが国においては,Lバンド帯が地殻変動の検出に適している.Lバンドで運用する,我が国の衛星であるJERS,ALOS,およびALOS2に搭載されたLバンドSARによる観測は,火山観測に多くの成果を上げており,すでに火山性地殻変動モニタリングの標準的手法となっている.一方,航空機搭載SARは,ある程度自由に照射方向を設定でき,即時対応にも適応性が高いことから,衛星型にない多くの長所を有しているものの,リピートパス干渉法(DInSAR)に関しては,飛行軌跡の制御や位置追跡の精度がボトルネックとなっており,広く用いられる段階には至っていない.我々は,Lバンド航空機SAR干渉技術の高度化を目ざす研究を実施中であるが,その一環として,九州に点在する活火山(桜島火山,霧島火山,雲仙火山)を対象とした観測を宇宙航空技術研究開発機構(JAXA)が運用する航空機SARシステム(Pi-SAR-L2)を用いて実施した.特に,桜島火山と霧島火山に対しては,リピートパス干渉を目的とした観測をJAXAが2013年以降実施中であり,2017年に新たに観測したデータは,これらの既存データとの干渉処理が可能である.本報告では,衛星観測や水準測量の結果から,最近の数年間における膨張性の変動の存在が明らかになっている,霧島火山硫黄山をターゲットとした,航空機SAR干渉解析結果を紹介する.2.Pi-SAR-L2データの干渉処理についてPi-SAR-L2には,高精度なINS-GPSハイブリッド型の航路追跡装置が搭載されている.これにより,高い軌道再現性が実現されているが,予測不可能な気流の変化等の影響によって,完全な同一航路の実現は困難である.その結果,航空機干渉SARの飛翔航路偏差起源の位相差の分布は,衛星のそれに比べてより複雑な形状を呈し,地殻変動情報の有効な抽出には,航跡の複雑性に起因する位相分布の適切な除去が必要である.これまでの予備的な解析から,ペアを構成する主画像(Master)および従画像(Slave)の位置合わせの達成度が干渉性をほぼ支配することが分かっているので,そのプロセスの確実性を高めることに,解析の主眼をおいた.位置合わせは,主従の元画像どうしのピクセルのズレの検出と,同一地点のピクセルどうしが同じ位置に来るように従画像の位置ズレを補正するリサンプルのプロセスで構成される.解析に使用する干渉SAR解析パッケージであるRINCでは,基本的に二次関数によるズレ予測に基づき,成功領域を徐々に拡大するアルゴリズムに基づいており,今回は,その特徴を生かしつつ確実なズレ検出を全画面領域において成立させるため,手動による確認プロセスを挟みながら,繰り返し処理を行い,徐々に成功領域を広げる処理を定式化した.この方法により,解析した6ペア全ての全画面において,ズレの高精度検出に成功した.一方,従画像(Slave)のリサンプリングは,現状では,二次関数による多項式近似を依然として採用しており,より複雑な変化をする飛行(Azimuth)方向において,数ピクセルにおよぶ残差が残留しており,これが,帯状の干渉不良領域を発生させている原因であることも分かった.二次関数では近似できない高周波成分に対しても対応する,より高精度なリサンプリングを実施すれば,干渉度がさらに改善される期待がある.いずれにせよ,位置合わせ手法を改善した結果,斑状の干渉不良領域は残存するものの2014,2016,2017の三時期に実施された三方向から観測から構成した,全6ペアの全てについて,全領域の干渉が達成された.3. 霧島火山の干渉解析結果リサンプリング手法に改良の余地があり,全画像での均一な干渉は実現できなかったが,現時点でも,かなりの領域で良好な干渉が達成されており,衛星SARや水準測量など他観測から地殻変動の存在が確認されている硫黄山において,航空機SARでも火山活動に対応すると考えられるフリンジが確認できた. 全て膨張性の変動が示唆される結果となっており,視線方向の距離変化の大きさは数cm程度であった.航空機干渉SARでは,上空の風速の変化などの影響で,飛行軌跡が小刻みに変動するため,航空機搭載のGPSによる航跡情報では補正しきれない,軌道縞や地形縞が残存する.今回の解析では,空間的長波長のフリンジは,全て軌跡起源と解釈して,一律に取り除く処理をしたため,空間的に長波長の地殻変動も取り除かれる恐れがある.しかし,硫黄山の地殻変動は力源が浅く,局所に偏在する変動と考えられるため,今回の手法でも確実に捉えることができたと考えられる.なお,長波長成分除去後のノイズレベルの見積もりは,概ね+-2cm程度であるなお,航空機SARの特性を生かし,多方向からの観測が行われており,地殻変動の三次元化を実施し,講演時にはその結果についても報告する.謝辞:干渉SAR解析には小澤拓博士が開発したRINC(Ver.0.36)および国土地理院の標高データを使用した.ここに記して感謝する.
著者
加藤 倫平 植田 勇人 吉田 孝紀
出版者
日本地球惑星科学連合
雑誌
日本地球惑星科学連合2016年大会
巻号頁・発行日
2016-03-10

北海道中軸部の空知-エゾ帯の主要な構成メンバーとして南北に渡って広く分布する白亜系堆積物は蝦夷層群と呼ばれ,多くの層序学的・古生物学的研究が重ねられてきた。その層序は従来から下部,中部,上部に区分されている(Matsumoto, 1942など)が,蝦夷層群の下部と中部の境界の層序学的位置は各地で異なり,蝦夷層群の広域的な対比は困難である。一方,蝦夷層群分布域全域に渡って追跡できる珪長質な凝灰岩や火山砕屑性砂岩からなる鍵層が確認され,北海道中央部では丸山層と命名された (Matsumoto, 1942;本山ほか,1991など)。この鍵層の堆積年代は蝦夷層群の広域的層序対比に有効であると考えられる。しかし,空知-エゾ帯北部の天塩中川地区や南部の春別川地区においても,化石の産出が乏しいため,この丸山層に対比されるような地層の特定は困難である。このような化石の産出に乏しい地層の年代決定には,U-Pb法による砕屑性ジルコンの年代値が有効である。よって,本研究では,天塩中川地区と春別川地区の凝灰岩層中の砕屑性ジルコンを用いて,U-Pb法により堆積年代を決定し,蝦夷層群下部~中部に及ぶ層準の広域的層序対比を試みた。ジルコンのU-Pb測定には新潟大学のLA-ICP-MS (Agilent 7500a ICP-MS及びNew Wave UP213 レーザーアブレーションシステム)を用いた。 その結果,天塩中川地区の蝦夷層群中部の白滝層の凝灰岩の堆積年代として,96.3-103.4Maを得た。一方,春別川地区の蝦夷層群中部春別川層の凝灰岩層の堆積年代として,98.5±0.5Maを得た。これらの凝灰岩及び凝灰質砂岩の岩石学的性質は珪長質であることが分かった。すでに報告されている北海道中央部の丸山層の堆積年代は浮遊性有孔虫の検討から,102-105Maである(高嶋ほか, 1997b)。また,北海道中央部の日陰の沢層の凝灰岩層のサニディンから98.98±0.38Maおよび99.16±0.37Maの Ar-Ar年代がすでに報告されている(Obradovich et al., 2002)。これらの凝灰岩層は,96-127MaのK-Ar年代を示す渡島帯の花崗岩類(柴田・山田, 1978など),100.6±3.3MaのAr-Ar年代を示す火山岩類(滝上, 1984など)に由来すると考えられる。
著者
佐藤 勲
出版者
日本地球惑星科学連合
雑誌
日本地球惑星科学連合2015年大会
巻号頁・発行日
2015-05-01

日本語で書かれている外国人の名前には、不適切な表記のものが存在する。科学における代表的な例は、「ハレー彗星」や「コリオリの力」である。これらは、「ハリー彗星」、「コリオリスの力」と表記するのが正しい。このような誤った表記が放置されているのは、これまで学会が外国語の名前の発音調査や規制をきちんと行ってこなかったことが大きな原因である。 食品の虚偽表示問題を受けて、昨年末に景品表示法が改正され、あらゆる商品とサービスを対象として虚偽表示に対して課徴金が課せられるようになった。学会としてもこのような状況に鑑みて、外国語の正確な表記に努める必要があろう。
著者
伊藤 颯矢 菅原 篤弥 操 知希 三野 正太郎
出版者
日本地球惑星科学連合
雑誌
日本地球惑星科学連合2018年大会
巻号頁・発行日
2018-05-11

本研究は,本校卒業生の鈴木湧平氏が,自身の撮影した流星の画像を見て,流星の発光途中で色が変化していることに疑問を持ったところから始まる。そこから発展し,流星の発光に由来する元素の特定,発光の時間変化における各波長の輝線の輝度変化の測定を目的に研究している。研究の手順として,まずデジタル一眼レフカメラと回折格子を用い,流星が縦に写るように角度を調節して連続撮影をする。流星の分光画像が撮影できたら,波長測定の基準を作るために校正作業を行う。そこで得た数値をもとに流星から分光した輝線の波長を求め,文献値と比較して発光に由来する元素を特定する。また,輝度測定の範囲を各輝線上に設定して,輝度変化の様子を複数の輝線で比較する。これまでの結果として,散在流星(2015)とおうし座北流星群(2017)でNaを,ペルセウス座流星群(2016)とふたご座流星群(2017)で複数の元素を特定した。特に,ふたご群ではオーロラで見られる酸素禁制線発光も確認できた。さらに,27例撮影したふたご群において,Naの欠乏と考えられる結果もあった。また,ペルセウス群の流星では,Naの輝度がMgの輝度よりも小さくなる部分があったり,輝度のピークが元素ごとに違っていたりした。
著者
Toshiro Tanimoto Anne Valovcin Vala Hjorleifsdottir Felix Rodríguez Cardozo
出版者
日本地球惑星科学連合
雑誌
日本地球惑星科学連合2018年大会
巻号頁・発行日
2018-03-14

From satellite data, we can learn the location of a hurricane (typhoon) quite accurately. What is hard to get from satellite data is its intensity change and that is why hurricanes hunters (NOAA) still fly into a hurricane to measure pressure and wind near its center. We are exploring an alternative idea of using seismic signals to monitor its intensity change. Our main questions are (i) when a hurricane (typhoon) intensifies, what kind of seismic waves are generated and can be measured on land and (ii) what is the best seismic indicator of an intensifying hurricane?In this study, we take an example of Hurricane Patricia in 2015 which was the strongest hurricane on record in the eastern North Pacific and North Atlantic basins (Kimberlain et al., 2016). We examine the nature of seismic waves that were excited as the hurricane intensified over anomalously warm waters to the south of Mexico and then reached land near Playa Cuixamala (Oct. 23, 23th hour, UTC). This landfall area was surrounded by many seismic stations from the Mexican National Seismological Service (SSN) which has broadband seismic data (STS-2 type velocity sensors).We examined temporal variation of power spectral density (PSD) of seismic data for selected 32 stations using the time-frequency plots. The landfall occurred near the center of 32 stations and depending on the distance from the hurricane, we can categorize the characteristics of time-frequency plots into three types, the northern group, the central group, and the southern group. Stations in the central group recorded the time evolution of hurricane intensity most faithfully and for monitoring purposes, we should focus only on these stations. The southern group showed signals from the hurricane but it also showed strong effects of nearby coastal (trapped) waves which made it problematic for monitoring purposes. The northern group seems too far away in general as seismic signals from the hurricane became much weaker.This hurricane developed from category 1 to category 5 within 24 hours on Oct. 22, 2015, while it was off the coast about 350-400 km; the data in the central group shows a rapid increase of seismic energy around 0.15 Hz (between 0.1 and 0.2 Hz) which matches with intensification of wind speeds of this hurricane. This frequency band is similar to those of secondary microseisms. The integrated power between 0.1 and 0.3 Hz shows a sudden increase of power during the first 6 hours of Oct. 22 and remains at this high level for about a day when the hurricane was the strongest. Seismic amplitudes decay from the coastal stations to the interior, suggesting (obviously) that coastal stations are critical for monitoring purposes.The data in the central group showed a rapid increase of seismic energy for higher frequencies up to 1 Hz too, although the maximum energy was at about 0.15 Hz. We also noted that there were no significant changes in seismic data below 0.1 Hz.These observations indicate that we should examine the frequency band from about 0.1 Hz to 0.4 Hz focusing mainly on the stations within about 500 km from the hurricane center. Since the location of hurricane is critical, seismic monitoring has to be conducted together with satellite data.
著者
津旨 大輔 坪野 考樹 三角 和弘 立田 穣 豊田 康嗣 恩田 裕一 青山 道夫
出版者
日本地球惑星科学連合
雑誌
日本地球惑星科学連合2018年大会
巻号頁・発行日
2018-03-14

A series of accidents at the Fukushima Dai-ichi Nuclear Power Plant following the Great East Japan Earthquake and tsunami of 11 March 2011 resulted in the release of radioactive materials to the ocean by two major pathways: direct release from the accident site and atmospheric deposition. A 6 years, regional-scale simulation of 137Cs activity in the ocean offshore of Fukushima was carried out by the Regional Ocean Model System (ROMS), the sources of radioactivity being direct release, atmospheric deposition, the inflow of 137Cs deposited into the ocean by atmospheric deposition outside the domain of the model, and river discharges.Direct releases of 137Cs were estimated for 6 years after the accident by comparing simulated results and measured activities adjacent to the accident site. In addition, river discharge rates 137Cs were calculated by multiplication between river flow rate and 137Cs activity. River flow rates were simulated by a water circulation analysis model for each catchment. Temporal change of 137Cs activity both of particle and dissolved forms were measured at 8 rivers and normalized by the inventory of 137Cs in each catchment. 137Cs activity in other 4 rivers were estimated by the normalized 137Cs activity and inventories of catchments. After 2013, direct release and river discharge were dominant for input of 137Cs to the ocean. Apparent half-life of direct release and river discharge of were estimated to be about 1 year and 2 years, respectively.Apparent half-life of measured 137Cs activity adjacent to 1F NPP was about 1 year, on the other hand, the ones in the coastal zone away from 1F NPP were about 2 years after 2013. Apparent half-life of simulated results with river discharge was in good agreement with the one in the coastal zone away from 1F NPP. River discharge affected on temporal change of 137Cs activity there. On the other hands, simulated 137Cs activities with river input were one order of magnitudes smaller than observations. This underestimation suggests modifications of river input process, such as estuary mixing process, removal from particle form 137Cs and inputs from small rivers around the 1F NPP.