はじめに
現在の飛騨山脈の気候環境では，降雪量が融解量を上回ることができない．そのため，飛騨山脈における氷河と多年性雪渓の分布は，吹きだまりやなだれの地形効果がある場所に限定される．樋口（1968）は，多年性雪渓の地形効果による涵養様式を雪渓の分布高度で分類しており，稜線からの標高差と分布高度が小さい雪渓を「吹きだまり型」，稜線からの標高差と分布高度が大きい雪渓を「なだれ型」，稜線からの標高差が小さく，分布高度が大きい雪渓を「混合型」とした．本稿では，「吹きだまり型」を「吹きだまり涵養型」，「なだれ型」と「混合型」を「なだれ涵養型」と呼ぶ.
質量収支の一般的な測定方法は，氷河や雪渓上に雪尺を打ち込み，１年後の雪面の高度変化を測る雪尺法が用いられる．しかしながら，日本の山岳地域は，膨大な涵養量と消耗量のため，雪尺が倒れてしまい実測できない．そこで，日本の多年性雪渓の質量収支観測では三角測量やトラバース測量がおこなわれている．観測実績のある雪渓は三角測量やトラバース測量の実測が可能な「吹きだまり涵養型」の小規模な多年性雪渓に限定されており，氷河を含む「なだれ涵養型」の多年性雪渓の質量収支は明らかでない．
本研究では，セスナ空撮とSfMソフトを使用し，2015～2018年の飛騨山脈北部の氷河と多年性雪渓の質量収支を算出し，「吹きだまり涵養型」と「なだれ涵養型」の違いを考察した．さらに，氷河の可能性が高い唐松沢雪渓において氷厚と流速を測定し，唐松沢雪渓の氷河の可能性について検討した．
研究手法
ⅰ）なだれ涵養型の氷河と雪渓の質量収支
飛騨山脈北部の立山連峰の「御前沢氷河」，「内蔵助氷河」，「三ノ窓氷河」，「小窓氷河」，「はまぐり雪雪渓」，「剱沢雪渓」，後立山連峰の「白馬大雪渓」，「カクネ里氷河」の8つの氷河と雪渓において，2015～2018年の春と秋に小型セスナ機からデジタルカメラで空撮を実施した．空撮画像と2次元の形状から3次元形状を作成するSfMソフトを用いて，多時期の高分解能の数値表層モデル（DSM）を作成した．これらDSMの比較から「なだれ涵養型」の氷河と雪渓の高度変化を算出した．
ⅱ）唐松沢雪渓の氷厚と流動
後立山連峰の唐松岳の北東斜面に位置する「なだれ涵養型」の唐松沢雪渓において，2018年9月に，中心周波数100MHzの地中レーダー（GPR；GSSI社製）を使用し，雪渓の縦断方向に2列と横断方向に６列の側線で測定を実施した．縦断方向と横断方向で反射波のクロスチェックをおこない正確な氷厚を求めた．GPRの解析結果をもとに氷厚の大きい上流部において，アイスドリルを用いて，長さ4.5ｍのステークを5地点に設置し，GEM-1（測位衛星技術社製）で9月末と10月末でGNSS測量を実施し，2時期のステークの位置情報の差から流動量を求めた．9月末に水平に打ち込んだステークは，再測時の10月末でも水平を保っていたことから，積雪のグライドやクリープは，流動に関与していない．
結果
ⅰ）8つの氷河と雪渓の質量収支は，2015/2016年に全域が消耗域となり，2016/2017年に全域が涵養域となり，2017/2018年にパッチ状に涵養域と消耗域が点在する結果であった．2015～2018年の融雪末期の中で最も氷河と雪渓の規模が小さくなったのは，小雪年の2016年である．2016～2018年までの高度変化で，「吹きだまり涵養型」の雪渓と「なだれ涵養型」を比較したところ，「なだれ涵養型」のほうが大きく涵養していた．ⅱ）唐松沢雪渓では，GPRの結果から最深部で30ｍほどの長さ1㎞ある氷体を確認した．流動観測では，ステークは１ヶ月間で斜面方向に20ｃｍほど動いており，得られた氷厚と流速は氷河の流動則とほぼ一致した．
