火山 第 57 巻 ( 2012)
第 4 号 219-225 頁

寄 書

一般用デジタルカメラを利用した簡易 SO2カメラの製作
橋 本 武 志＊・寺 田 暁 彦＊＊・江 尻省＊＊＊・中 村 卓 司＊＊＊・阿 保真＊＊＊＊（2012 年 4 月 27 日受付，2012 年 9 月 25 日受理）

A low-cost SO2 Imager with the Use of Digital Cameras of Consumer Use
Takeshi HASHIMOTO＊, Akihiko TERADA＊＊, Mitsumu EJIRI＊＊＊,Takuji NAKAMURA＊＊＊ and Makoto ABO＊＊＊＊We produced an SO2 imager of a low-cost version with the use of digital cameras of consumer use. Generalconfiguration and characteristics of the instrument are presented. Calibration with SO2 cells of known column
concentration confirmed a comparable absorption coefficient to the device which is previously invented by Mori andBurton (2006). Although simultaneous shooting with two cameras is necessary to overcome the temporal change of an
object because of somewhat long exposure time (5 to 10 sec) required, this device operates without external power orcontrol PCs, and thus, is suitable for mobile use. We performed a field test of the instrument at Sakurajima volcano andconfirmed its validity as an SO2 imager. However, further improvement of the optical system or a special care with UV
scattering in front of a plume is necessary for quantitative applications in a field operation.
Key words : SO2 camera, discharge rate, plume, Sakurajima volcano
1．は じ め に

パスフィルタを組み合わせて SO2濃度分布を画像として火山噴煙に含まれる SO2は，白色噴煙や有色噴煙とし得る装置も開発されている（例えば，Mori and Burton,て認識される H2 O や火山灰と異なり，視認することが2006）．この手法の最大の特長は，いうまでもなく SO2難しい．SO2を簡便に可視化できれば，火山活動の監視ガスを画像として視認できることである．COSPEC やや活動評価の指標として有用である．また，その測定にDOAS のような単視線の測定では，特に濃度の低い状況定量性や高時間分解能が付加されれば，噴火準備過程やにおいて背景レベルとの区別が難しいケースがある．こ爆発機構の研究にも大いに貢献できる．

れに対して，画像で濃度分布のスナップショットが得ら

近年，火山ガスの遠隔 SO2測定は，従来型の相関スペ

れる SO2カメラの場合には，パターン認識によって背景クトロメータ (COSPEC) に代わって，DOAS (Differentialとの分離は容易となる．また，測定原理は単純で，特殊Optical Absorption Spectroscopy)と呼ばれる手法を用いな解析技術を必要としない．このような特長を有する

た単視線 の 紫 外 分 光 装 置（例 え ば，Galle et al., 2003;

SO2カメラは，脱ガス機構を含めた基礎研究から，火山McGonigle et al., 2002）による，パニングもしくはトラ活動のモニタリングに至るまで幅広く応用できるであろバース測定が一般的となっている．これとは別に，いわう．しかし，高性能の計測用 CCD を用いた装置は高価

ゆる SO2カメラと呼ばれる，CCD カメラと紫外域バンドであり，一般の火山研究者や行政機関が手軽に入手でき＊〒060-0810 札幌市北区北 10 条西 8 丁目
北海道大学大学院理学研究院付属地震火山研究観測センター
Institute of Seismology and Volcanology, Faculty ofScience, Hokkaido University＊＊〒377-1711 群馬県吾妻郡草津町草津 641-36
東京工業大学火山流体研究センター
Volcanic Fluid Research Center, Tokyo Institute of Technology＊＊＊〒190-8518 東京都立川市緑町 10-3
国立極地研究所
National Institute of Polar Research＊＊＊＊〒191-0065 東京都日野市旭が丘 6-6
首都大学東京システムデザイン研究科
Graduate School of System Design, Tokyo Metropolitan
University
Corresponding author: Takeshi Hashimotoe-mail: hasimoto@mail.sci.hokudai.ac.jp

220

橋本武志・寺田暁彦・江尻省・中村卓司・阿保真

るものとは言い難い．そこで，著者らは，安価に入手可
能な一般用のデジタルカメラを利用して簡易型の SO2カメラを製作した．後述のように，測定時間分解能や野外測定における定量性にはまだ問題が残されているものの，本装置は誰にでも容易に製作できるうえ，火山噴気に含まれる SO2を可視化するために十分な性能を有している．本稿では，その構成と基本性能について，観測例を示しつつ報告する．
2．装置の基本構成
本研究で製作した装置の基本構成と装置の外観を Fig.
1 に示す．本研究では，受光装置として一般撮影用のデジタル CCD カメラ (Nikon D70) を使用した．撮影は通常の写真撮影と同様に行い，画像は各カメラの内蔵メモリに保存される．このため，外部電源や制御用 PC などの特別な仕組みは不要で，可搬性・軽量性に優れている．
CCD 素子の大きさは 23.7×15.6 mm であり，総画素数は 6.24 メガピクセル，輝度分解能は 12 bit である．一般に，CCD 素子は紫外域から近赤外域までの幅広い分光Fig. 1. (Top) Basic configuration of the SO2 camera developed in the present study. (Bottom) SO2 camera in field
operation. Two CCDs are mounted alongside on a tripod.

感度特性を有する．そのため，可視領域を撮影対象とする一般のデジタルカメラには，モアレ防止等を目的とし両波長帯でほぼ同程度の背景輝度を得ることとした．上た光学式ローパスフィルタが装着されている場合があ記の構成において，自然散乱光を光源とした場合の適正り，紫外線を受光するためにはこれを除去しなければな

な露出時間は 5〜10 秒程度であり，2 つの波長帯を 1 台らない．また，通常の写真撮影に用いられるレンズは紫のカメラで交互に測定しようとすると撮影対象の時間変外線を十分に透過しないため，本研究では，石英ガラス化が無視できなくなる．これを避けるため，2 台のカメで作られた紫外線撮影用のレンズ Nikon UV-105 mm (Fラを用いて 310 nm 帯と 330 nm 帯を同時に，同じ露出時4.5)を使用した．波長 220-900 nm における透過率は公間で撮影する方式を採用した．この方式は既に Kantzas

称値で 70 % である．

et al. (2010) や Kazahaya et al. (2011) でも用いられていCCD が感度を有する波長帯のうち，本研究では SO2のる．得られた 2 枚 1 組の画像（正確には，両波長帯それ吸収帯に相当する 310 nm と，非吸収帯に相当する 330ぞれの背景画像を含めて 4 枚）から SO2濃度を抽出する

nm を使用する．これらは，バンドパスフィルタを用いには，次章で説明する一連の演算を PC 上で後処理する．て抽出することができる．本研究では，朝日分光製の狭なお，本研究で用いたカメラでは，CCD 素子に固有の帯域バンドパスフィルタ（BPF1：半値全幅約 9 nm，透過ダークノイズの減算処理は，撮影時に自動で行う機能を率 65 %，阻止率 10−5）に加えて，このフィルタで阻止し有している．

きれない近赤外域の漏れ込みをカットするために，やや本論文で記述した光学系の構成では，量産品ではない広めの紫外通過バンドパスフィルタ 2 種（BPF2：朝日分UV レンズの価格が全体のコストを決めている．ただ光製，中心波長 320 nm，半値全幅 220 nm，透過率約 98

し，画像の歪みや収差の面で妥協するならば，透過率の%，BPF3：朝日分光製 RR340，中心波長 340 nm，半値全高い石英製の単レンズに光学系を変更することで，極端幅約 80 nm，透過率約 80 %）を併用した．

な低コスト構成とすることも可能である．

310 nm 帯と 330 nm 帯では，CCD の感度や太陽光の強度が異なる．この影響を考慮せずに撮影を行うと，いず3．輝度-SO2カラム濃度特性の実測試験

れかの画像が CCD のダイナミックレンジを超えて露出撮影される画像の輝度情報からは，SO2のカラム濃度，不適正となってしまうが，開口絞りを変えると光学系のすなわち視線方向に積分された濃度が得られる．まず，幾何学的特性に影響を与える可能性があるため不都合でその原理について説明する．SO2 による紫外線吸収がある．そこで，本研究では，330 nm 帯の光学系に ND8Lambert-Beer の法則に従うとすると，入射光と透過光のフィルタを 2 枚追加して光量を約 100 分の 1 に落とし，強度比は以下の式で表すことができる．

一般用デジタルカメラを利用した簡易 SO2カメラの製作
I λI λ=exp −ελcd

(1)

ここで，I0は入射光の強度，c は SO2の濃度，d はガスAi, j=−ln221Pi, j Pi, j +ofs
Pi, j Pi, j (4)の厚みである．ε は吸光係数で，ガス種に固有の定数でここで，P (i, j) は測定画像上のあるピクセルにおけるあるが，波長依存性をもつため ε (λ) と表記した．両辺の生の輝度値を表し，添字の数字部分は，310 nm または自然対数をとり符号を反転すると，吸光度 A が得られる．

330 nm の画像であることを示す．また，添字の s は SO2A=−ln I λI λ=ελcd(2)ガスを透過してきた測定対象（この実験ではガスセルの窓部分）の画像であることを表し，bg は，ガスセルを除吸光度は，吸収減光がない場合（すなわち I =I0の場去した状態で背景散乱光だけを撮影した画像であること合）にゼロであり，ガス濃度と厚みの積，すなわちカラを表している．なお，これは Mori and Burton (2006) の見ム濃度 Ccolに比例する．これより，カラム濃度は，以下掛け吸光度と同じ定義である．右辺第 1 項の分母と分子の式で与えられる．

は，それぞれ 330 nm と 310 nm 帯における背景輝度を基c≡cd=Aε(3)準とした輝度比を表している．それぞれの背景画像 (bg)で除算しているのは，主にレンズの周辺減光を除去する

実際には，複数の吸収波長を含む有限の幅をもった波ことを目的とした処理であり，原理上は必須ではない．

長帯で測定を行うため，実効的なε の値（以下，実効吸光(2) 式からわかるように，カラム濃度と吸光度の間に係数と呼ぶ）はフィルタ等の光学系に依存する．そこで，は線形関係がある．4 種のガスセルのカラム濃度に対し(3) 式では，ε の代わりに εeff と表記した．機器毎に予めて，セルの窓部分の輝度を平均処理して (4) 式の右辺第この実効吸光係数を実測しておく必要がある．本研究で1 項を求めてプロットし，直線近似したときの傾きが (3)は，カラム濃度既知の複数の SO2ガスセルに晴天時の自式の実効吸光係数 eeffである．ただし，生の輝度画像に然散乱光を透過させ，後述の (4) 式で定義される見掛けは，ガスセルのガラスによる減光が含まれている．ガス吸光度を測定した．カラム濃度に対する見掛け吸光度のセルの窓部分の輝度と bg 画像の輝度を，SO2の吸収がなプロットを Fig. 2 に示す．その傾きからεeffが得られる．

い 330 nm 帯で比較することにより，ガラスによる減光

以下，この校正実験について詳述する．

が 10〜12 % であることが確認できたので，実効吸光係

ガスセルの厚みは 10 mm で，カラム濃度は 90，192，数の計算には，これを補正した輝度値を用いた．なお，

497，および 1011 ppmm の 4 種を用いた．Fig. 1 に示しFig. 3 に示したガスセルの透過特性グラフでは，カラムたように，この校正実験ではレンズの前にガスセルを置濃度によらず長波長側で約 15 % の減光が認められるが，き，ガスセルの窓が画角の中心付近になるようにカメラ

これには，ガスセルのガラスによる減光に加え，セルのを固定する．光源強度の時間変化を避けるため，4 種の窓が分光光度計の光源光路径よりもやや小さかったこと

ガスセルの置き換えは数分間で完了するように撮影をによる影響が含まれている．右辺第 2 項のオフセット行った．参考までに，ガスセルとバンドパスフィルタ，(ofs) 値は，上記近似直線の y 切片に対応するもので，カおよび ND8 フィルタの波長透過特性を Fig. 3 に示す．

ラム濃度が 0 のときに Aa =0 となるように決定した．

一連の校正撮影後に画像を PC に取り込み，セルの窓(4) 式の定義より，オフセット値は理論上ゼロになるべ部分の輝度を読み取る．本研究で用いたカメラの場合，きものであるが，実際の値は-0.02 程度である．なお，実ひとつの画素について RGB の 3 つの素子があり，三原際のフィールド測定では，上記特性試験で得た実効吸光色分解するためにそれぞれにカラーフィルタが施されて係数を用いるが，オフセット値については，特性試験のいる．このカラーフィルタは，可視光を透過するもので値を参考にしながら，噴気から離れた方向のカラム濃度あるが，紫外域にも漏れ特性がある．実験の結果，紫外値がゼロになるよう試行錯誤的に調整している．Fig. 2

光に感度のあるのは，RGB のうち R チャンネルのみで，では，この ofs 値を加味した吸光度 Aaを示しており，実G チャンネルおよび B チャンネルには殆ど有効な感度

効吸光係数は約 1.5×10−4 ppmm−1であった．

がないことが確認された．このため，以下に示す輝度とは，すべて R チャンネルについてのものである．
画像上のピクセル毎の見掛け吸光度を，(2) 式を拡張した次式により定義する．

この結果から，輝度 1 ビットの変化に対応するカラム濃度値の分解能を計算することができる．カラム濃度分解能は，背景輝度をどのレベルに設定するかに依存し，かつ，どの輝度範囲を対象としているかによっても変化する．また，カラム濃度が低い部分ほど分解能が高い．階調 8 ビットの JPEG フォーマットを用い背景輝度を 200

222

橋本武志・寺田暁彦・江尻省・中村卓司・阿保真

4 倍である 6.1×10−4 ppmm−1 が得られる．実測値と理論値がこのように大きく食い違っているのは，CCD の波長感度特性がバンドパスフィルタの帯域幅の中で平坦ではないことによる可能性が高い．一般に CCD の感度特性は短波長側で急激に低下する一方，SO2による吸収は短波長側で強くなる．このため，CCD の感度特性を含んだ実際の総合特性では実効吸光係数が理論値よりも小さくなっているものと考えられる．
4．野外動作試験
Fig. 2. Apparent absorbance with respect to the column
concentration of calibration cells. Error range indicates 1σ of 100×100 pixels. Asterisks indicate thevalues from Mori and Burton (2006).

次に著者らは，本研究で製作した SO2カメラの野外性能を確認するために，実際の火山噴気を対象とした動作試験を行った．野外試験の対象として，ほぼ安定して多量の SO2ガスが放出されている桜島火山の昭和火口を選定した．2011 年 3 月 31 日の午前 10 時から 12 時にかけての約 2 時間，昭和火口から南南西約 3 km の位置にある古里運動公園において分光撮影を行った．この位置からは，昭和火口を直視することはできないが，火口から100 m 程度上昇したところから，稜線を越えて噴気が視認できる位置関係にある．当日は快晴で，北西の弱い風が吹いていた（鹿児島地方気象台における風向風速は北東 2〜3 m/s，湿度 46 %，気温 14〜16℃）．適正な条件で分光撮影を行ったのは 11 時 30 分〜12 時 00 分の約 30分間である．レンズ絞りは開放，露出時間は 8 秒間とし，Fig. 3. Spectral transmission characteristics of the twoband-pass filters (BPF1 for 310 nm and 330 nm) and of
the calibration cells with known column concentrations(90, 192, 497, 1011 ppmm). These are shown as solidlines with the scale at left-hand side. Transmittance of
the ND8 filter (two-ply) is also shown as a broken linewith the scale at right-hand side.

撮影は 1 分間隔で行った．一連の噴気撮影の最初と最後に背景画像 P(i,j)310bg，P(i,j)330bgを取得した．背景画像としては，火口方向から仰角を変えずにカメラを水平回転させ，噴気のない方向の青空を撮影した．
なお，撮影時間帯において昭和火口で爆発は発生しなかった．火口から立ち上る噴気は透明度の高い白色であり，目視で判断する限り，火山灰の混入しない状態でほぼ連続的に放出されていた．測定時間中，噴気の流れるとした場合では，理論上，輝度 200 付近で 34 ppmm/bit，方向や速度には殆ど変化が見られなかったが，噴気水蒸輝度 100 付近で 67 ppmm/bit である．CCD のノイズがあ気の濃さには若干の揺らぎが繰り返されていた．

るために，実際上の分解能はこれよりも悪化するが，画11 時 35 分のデータを Fig. 4 に例示する．まず，2 枚像の平滑化処理を行うと，実用上問題のない程度まで改の紫外画像を比較すると，310 nm 帯にのみ，SO2の吸収善される（Fig. 2 の誤差範囲は，100×100 ピクセルあたによる有意な減光があることが視認できる．310 nm とりの標準偏差±1 σで表現している）．なお，周辺減光を330 nm の 2 枚の画像（それぞれの背景画像を含めると計問題にしない場合には，(4)式で背景画像による除算を行4 枚）から，(4) 式によってピクセル毎の見掛け吸光度 Aaわない方が，ノイズはより押さえられる．Fig. 2 には，を求め，校正実験(Fig. 2)から得られた実効吸光係数 εeff

Mori and Burton (2006) の装置特性も併せて示したが，本の値を用いて，(3) 式によりカラム濃度に変換する．

研究で製作した装置は，実効吸光係数では彼らのものとほぼ同程度であることがわかる．

2 台のカメラの画角には僅かにずれがあるため，前処理として片方の画像に回転と平行移動を施す．また，カ

ところで，実測された校正用セルの吸収特性とバンドラム濃度値はピクセル毎のばらつきが大きいため，ここパスフィルタの透過特性 (Fig. 3) にもとづいて，310 nm

では，5×5 ピクセルの平滑化処理を行う．得られたカラ帯で期待される実効吸光係数を計算すると，実測値の約ム濃度画像（Fig. 4 右下）を見ると，噴気には 300〜600

一般用デジタルカメラを利用した簡易 SO2カメラの製作

223

Fig. 4. An example of spectral images of a plume from Showa crater at Sakurajima volcano viewed from Furusato
athletic park (c.a. 3 km SSW from the crater). Upper-left : 330 nm band; Upper-right : 310 nm band. Lowerleft panel shows the visible image at the same time (11h35 m JST). Yellow inset corresponds to the frame of
the upper panels. Lower-right shows the column concentration image produced from the upper images with acolumn concentration profile along the line AB.

Fig. 5. Column concentration snapshots of the plume from Showa crater on March 31, 2011 from 11 : 32 to 12 : 00 JST.
Each picture frame is the same as the right-bottom panel of Fig. 4 but reduced in 90 % widthwise for saving the space.
ppmm 程度の SO2ガスが含まれているのがわかる．濃度究の装置が火山噴気 SO2の可視化カメラとして機能しての高い部分が，昭和火口から画面の右上方向に向かっていることがわかる．

伸びていることが明瞭に認識できる．これにより，本研

次に SO2の放出率を以下の手順で見積もり，定量性に224

橋本武志・寺田暁彦・江尻省・中村卓司・阿保真

ついて検討する．まず，Fig. 4 右下図において，噴気にほぼ直交する直線 AB に沿ってカラム濃度を積分する（5式）．ここでは，AB の距離は，噴気の幅や中心軸が時間とともに変化することを考慮して，視認できる噴気の幅よりも広く設定している．直線 AB は，火口の直上を通りかつ視線方向に対して直交する位置に置き，その実距
離は，レンズの画角と火口までの距離に基づいて算出する．Δ l は，直線 AB に沿った画像上の 1 ピクセルあたりの実距離である．
C≡∑ c i, j Δl(5)さらに，カラム濃度画像から SO2 放出率への変換は，Fig. 6. Discharge rate of SO2 gas obtained from thesnapshots shown in Fig. 5. Ascent velocity is assumed to be 10-15 m/s. Error range corresponds tothat of the ascent velocity.

次の計算式を用いる．
Φ=CuMf(6)

M は SO2のモル質量，f は気体のモル体積で，ここで3に理解できる．本研究の装置で測定器としての定量性を確保するには，散乱の影響が小さい波長帯に光学系を変更するか，Kazahaya et al. (2011) が行っているようにはそれぞれ 0.064 [kg/mol]，0.0224 [m /mol]とする．u はDOAS との同時観測を併用して距離減衰の補正を施すな噴気の上昇速度である．上昇速度は，連続写真から噴気どの工夫が必要であろう．の特徴的な部分を追跡することにより 10〜15 m/s と見積もった．
上記の方法によって得た 2 分毎のカラム濃度画像を5．まとめ

著者らは，一般用のデジタルカメラを利用して比較的Fig. 5 に示す．濃度の高い部分が画面の中央部に集中し低コストの SO2カメラを製作した．本装置は誰にでも容て柱状に見える時間帯と，周辺部にも拡散して見える時易に製作できるうえ，火山噴気に含まれる SO2を可視化間帯があることがわかる．おそらく，後者の場合は噴気するために十分な性能を有している．また，光学系の選の上昇力が弱い状態で，ガスが滞留しがちになるために，択によっては，さらに低コストの構成も実現可能であり，濃度分布としてはやや高くなり，噴気の幅も広がったよ基礎研究以外への応用可能性も秘めている．

うに見えているものと考えられる．現有のデータでは，この装置について，濃度既知のガスセルを用いた校正画像毎に上昇速度を精度良く求めることが困難であった実験を行ったところ，Mori and Burton (2006) の装置とほので，本研究では上昇速度は上記 10〜15 m/s の中央値でぼ同程度の実効吸光係数が得られた．著者らは，この装固定し，誤差範囲を±2.5 m/s として画像毎の放出率に換置を用いて 2011 年 3 月 31 日に桜島で野外動作試験を行算した (Fig. 6)．

い，実際の火山噴気を対象とした SO2可視化装置としてその結果，この時間帯における放出率は，3.4±0.7 kg/sの有効性を確認した．一方で，定量性を確保するためにと求められた．これを 1 日あたりの放出率に換算する

は，噴気手前側の紫外線散乱の影響が小さい波長帯に光と，300±59 t/d となる．気象庁の測定によれば，2011 年

学系を変更するか，DOAS 等との同時観測を併用し，カ3 月の SO2 放出率として 700〜1,700 t/d の値が報告されラム濃度値の距離依存性を適宜補正するなどの工夫が必ている（気象庁，2011）．気象庁の報告値は，桜島南岳火要であることが確認された．

口と昭和火口からの放出率を分離して測定したものではなく，両者の和に相当するのに対して，本研究の測定は謝昭和火口からの放出率のみを測定している．また，本研校正用のガスセルは，京都大学地球熱学研究施設火山究の測定は非爆発時に行われたために SO2 放出率が低研究センターよりお借りしました．対応いただいた吉川かった可能性があることに加えて，測定距離が約 3 km慎氏に感謝いたします．ガスセルおよびバンドパスフィと長いため，Mori et al. (2006) が指摘したように，噴気ルタの透過特性測定には，国立極地研究所の光学実験室より手前の大気による紫外線散乱の効果が 310 nm 帯にを利用しました．北海道大学電子科学研究所の竹内繁樹

強く表れ，カラム濃度推定値が過小評価されている可能氏には，ND フィルタの透過特性測定にご協力いただい性が高い．これらのことを考慮すると，本研究の測定値たほか，研究全般にわたって重要な助言をいただきましが気象庁の報告値の約 20〜40 % と小さいことは定性的た．八神辞樹氏には，桜島でのフィールド測定の補助を

一般用デジタルカメラを利用した簡易 SO2カメラの製作していただきました．森俊哉氏と風早康平氏には，注意深い査読と建設的な助言をいただきました．両者の査読により本稿は大幅に改善されました．なお，本研究は，科学研究費基盤研究 (B) 21340122 による支援を受けました．
引用文献
Galle, B., Oppenheimer, C., Geyer, A., McGonigle, A.J.S.,
Edmonds, M. and Horrocks, L. (2003) A miniaturised
ultraviolet spectrometer for remote sensing of SO2 fluxes: anew tool for volcano surveillance. J. Volcanol. Geotherm.
Res., 119, 241-254.
Kantzas, E. P., McGonigle, A.J.S., Tamburello, G., Aiuppa, A.
and Bryant, R. G. (2010) Protocols for UV camera volcanic
SO2 measurements. J. Volcanol. Geotherm. Res., 194, 5560.
Kazahaya, R., Mori, T., Takeo, M., Ohminato, T., Urabe, T.
and Maeda, Y. (2011) Relation between single very-long-225

period pulses and volcanic gas emissions at Mt. Asama,Japan. Geophys. Res. Lett., 38, L11307, doi : 10. 1029/
2011GL047555.
気象庁 (2011) 日本の主な火山活動. 地震・火山月報（防災編）平成 23 年 3 月，41-54.
McGonigle, A.J.S., Oppennheimer, C., Galle, B., Mather, T.
and Pyle, D. (2002) Walking traverse and scanning DOAS
measurements of volcanic gas emission rates. Geophys. Res.
Lett., 29, 1985, doi : 10.1029/2002GL015827.
Mori, T. and Burton, M. (2006) The SO2 camera : A simple,
fast and cheap method for ground-based imaging of SO2 involcanic plumes. Geophys. Res. Lett., 33, L24804, doi : 10.
1029/2006GL027916.
Mori, T., Mori, T., Kazahaya, K., Ohwada, M., Hirabayashi, J.
and Yoshikawa, S. (2006) Effect of UV scattering on SO2emission rate measurements. Geophys. Res. Lett., 33,L17315, doi : 10.1029/2006GL026285.
（編集担当青木陽介）