第30号

国立防災科学技術センター研究報告

1983年3月

551，578．48

グライ ド現象の数理モデノレ
系内

口恭

明＊

国立防災科学技術センター雪害実験研究所

Mathematica1Modelof

Snow

G1idimg

By
Yas11akiNohguchi肋∫肋肋oグ∫〃o〃α〃北θS肋肋∫，N〃o〃α1地∫鮒肋Cθ〃θγ〃D乞∫α∫伽1〕榊θ〃ガo〃，jVαgαoんα，八吻g励α一冶2〃，940

Abstmct

Thels

not

snow

cover

a1waysuniform

conditions，andafu11depth

it

on

is

aslopebut

known

avalanche

fundamenta1mode1to

istends

tog1ide

decreases

orincreases

g1idingmotionthat

thereleased．The

describe

snow

adown

slow1y．Theaccording

of

snow

purposse

ofglidingmotions

thisglidetois

ve1ocitychanges

acce1erated

paper

is

inc1uding

toof
beforeproposeanonsteadymOtiOn．
Let！andrespectivelyl

contact

R

denoteforce

and

that1ぞisthefunction

equationofmotionSupPoseareaλ，thedrivingaresistant

of

force

on

the

ve1ocitychangeofλbeaground〃andsurface，

the

real

is

∫＝R（λ，o）

Sinceλchanges
dλduringmovement，consider

thedescribedas：＝一力（λ，o）十σ（λ，o）

dオ

where力andαareofλdue

toFor

the

rateofdecrease

formationofnew

rea1contact

example，if力，σ，Rareofλdue

toshearandthatofincrease

area，respective1y，assumed

力一舳，1−／舳ろ■λ）1熾1い一（1・・l1・）λ，
then

the

behaviorofsnow

gliding3−3，and3−4．Figs．6and7show
g1idingisglide

c1assified

into

four

cases：Figs．3−1，3−2，velocity

and

the

disp1acement

motion，respectively．

＊第1研究室一189一

dueto国立防災科学技術センター研究報告

1、第30号

1983年3月

まえがき

斜面積雪の地表面上の静的な滑りはグライドと呼ばれ，斜面上の構造物に及ぼす積雪の力の重要な要素となっている．それと同時に，全層なだれ発生に至るまでのグライドの変化の過程は，全層なだれの先駆現象として見ることにより，なだれ予知のための有効な情報として利用できることが可能であり（渡辺・五十嵐・山田，！978）その実用化が注目されている．

ある斜面について，グライドを1冬期問観察すると，比較的定常な動きをしている期問・加速している期問・減速している期問・静止している期間などが存在するのがわかる．とくに，全層なだれ発生の直前，あるいはクラツク発生の前にグライド速度の増加が見られることが知られている．しかし，グライドの加速が常になだれに結びつくとは限らず，また単純にグライド速度がある大きさになったときになだれが発生するともいい切れない．

一般にグライド速度の大きさは定常なときでcm／dayのオーダーであるが，地形・植生・気象・雪質などにより値は様々である．また，グライドの加速のような過渡現象の時問スケー

ルは，たとえば人工的にかやを張って滑りやすくした斜面上での数分〜数時問（中村・他，1972）というものから，積雪中に笹が埋まっているような笹地斜面での数十日（遠藤・秋田谷，1977）のオーダーのものまである．

このような積雪と地表面との接触の機構に応じたグライド現象の定性を把握するには体系だった数理モデルの枠組みをつくる必要があり，また実験や結果の解析をする場合にもその背後にある数理的な構造に対する洞察が必要となる．

斜面積雪のグライドは積雪に作用する斜面方向の力と積雪底面における地表との問の抵抗力によって規定される．したがってグライドのモデル化という意味は底面低抗力をモデル化することを指す．この底面抵抗力に対するもっとも簡単な表現としては，ふつうの摩擦力のように，斜面の垂直抗力に動摩擦係数（一般に速度に依存する）を乗じたものが考えられるが，このような表現では定常なグライドは表わせないし，積雪と地表面との問の物理的な過程を表現するものでもない．この意味で，笹地斜面におけるグライドの挙動の説明のために底面抵抗力を笹を積雪から引き抜くときの低抗力としてモデル化した遠藤（1980）の研究は，底面低抗カの時問的変化に物理的な解釈を与えるものとして意義がある．

本論文では，いろいろな接触機構に応じた，底面低抗力の時問変化を記述するための基礎となる数理モデルを提案し，いくつかの仮想的な接触機構を与えたときに現われるグライド現象の特徴について考察する．

2．底面抵抗のモデル

一190一

グライド現象の数理モデルー納口

積雪に作用する駆動力を！，底面低抗カをRとするときグライド現象のような静的な運動は力のつりあいの式を満たす．

！＝R
一般に沢は速度oの関数と考えてよいであろう．それと同時に，積雪と地表面との接触をになっているなんらかの物理量の関数である．この量は，たとえば巨視的な意味での積雪と地表面との真の接触面積であったり，あるいは微視的なレベノレでの真の接触面積であったり，また笹地の場合のように積雪内にとり込まれている地表の植物の量（これもある意味では巨視的な意味での真の接触面積と考えてよレ））であったりする．いまこの量をλと表わすことにする．このλは，積雪のグライドにともなって接触が断たれて減少する部分と，その一方で新たな接触が形成されることによって増加する部分を含み，その収支によりλの変化が規定されると考えてよいであろう．単位時問当りの減少量を力，増加量をσとする．これらは一般的にはλとoに依存するであろう．以上をまとめて底面低抗力Rを表現すると次のような形になる．

一阜二㌦）十州／1）

すなわちRの変化は，λの変化を通して表現されることになる．力＞αの場合はλは増加し，力＞αの場合はλは減少する．したがって接触状態が時問の経過につれて変化するという意味で非定常な状態である．力＝σの場合はλは一定で変化せず，したがって接触状態は定常である．

この底面低抗力のモデノレを用いると，駆動力！を与えたときのグライド速度ひを記述する方程式は次のとおりである．

㌃1㌦）十、〜）／・）

力のつり合いの式からoはλの関数として表わされる．したがって力，αはλだけの関数となり，微分方程式はλについて解くことができλの変化を通してoの変化が記述されることになる．

次の式
仁ぐ1む二11）十、（、引伽）／（4）を満す点λ、，o、は微分方程式の平衡点である．すなわちこの値に対してdλ／dチはゼロとな

一191一

国立防災科学技術センター研究報告第30号

1983年3月

（b）

（a）

べ、・、q

へ一・、

A困図1Fig．1

AAeuA（・）安定平衡点と（b）不安定平衡点（α）Stab1e．（ろ）Unstable

り，グライドは定常状態を保つ．ところでλの関数である力，σが図！の（・）のような関係にある場合を考えてみる．ヵとσの交点，すなわちλ＝λ、。ではdλ／dτがゼロでグライドは定常的である．一方λ＜λ、、ではdλ／d1は正となりλは増加し，λ＞λ、、ではdλ／dチは負となり減少する．このような平衡点は安定平衡点であり，なんらかの原因でλが平衡点から外れても，また再び平衡点にもどろうとする．すなわち，この場合定常的なグライドは安定して存在する．
これとは逆に力，σの関係が（b）のような場合にはλ、αは不安定平衡点であり，わずかのずれでλはんαからますます遠ざかるようになる．したがってこのような場合には定常グライドは安定には存在しえないことになる．

つぎに，動摩擦実験のように速度o＝o。を与えた場合を考えると，このモデノレでは低抗力は以下の方程式によって記述されることになる．

汽ぷ∴）十舳／（5）（3）式のoがλの関数であったのに対して，ここではoは定数o。であり，抵抗力Rがλの関数となる．このような設定で低抗力Rを求めることを以後定速グライド試験と呼ぶことにする．実際にモデルを検討するうえでは，単なる駆動力！を与えるグライド現象ばかりでなく，いろいろな状況設定での現象を記述できる必要があり，またそのことがモデルの有用性につながることにもなる．

この低抗力に定常状態が存在するならば，その時の平衡値R、，λ、は一192一

グライド現象の数理モデルー納口

lfl㍍1（ん，伽）＝O／（6）から求められる．（4）式からわかるようにヵ，σ，Rの関数形が同じである限り，すなわち両方の試験で接触機構が同じである限り沢、とo。の関係は！と伽の関係と同じである．これは直感的には当然のことと思われるかもしれないが，このモデルにおいてもそれが満たされていることになる．しかしそれぞれの場合で平衡値の安定性に関しては別の問題であることに注意しなければならない．

これまで接触をになっている機構が1つの場合だけを表現してきたが，さらに一般化して，いくつかの機構が接触に関与している場合は次のように拡張すればよい．
沢＝R（λ1，ん，ん，……，0）
dん一一カ1（λ、，λ、，λ、・・ …・，0）十σ1（λ1，ん，ん，…・，0）dオ

dん一一カ、（λ1，λ、，λ、，・ …・・，0）十σ。（λ、，ん，ん，…・0）

（7）

dオ

dん一一力、（λ1，λ、，λ、，・・・・，0）十α、（ん，ん，ん，……0）

d左

このモデノレを便ううえで本質的なことは，問題となるような過渡現象の時間スケールを支配するような積雪と地表間の主な接触の機構を設定することである．いいかえると，λを定めて力，α，R，をλとoの関数として表現することである．実際には様々な機構が作用するであろうが，スケールを設定することでそのうちのいくつかは一定として扱ってもよくなるであろう．いずれにせよ現象の本質を的確に表現しようというモデルの目的上，主要な効果と枝葉末節的な効果の選択が重要な問題といえる．

3．簡単なモデルによるグライド現象の考察

簡単ないくつかのモデルを与え，それから導かれるグライド現象の特徴について考察する．

3．1例1
λを真の接触面積とする．低抗カRは一般にはλの増加関数と考えてよいであろう．いまとくにλ＝0では低抗力もゼロとなるようなモデルを考えRを次のように表わす．

R＝（ε。十ε1。）λ（8）ε。，ε1は定数である．力については，静止状態ではゼロと考え，グライドが存在するときはひ一193一国立防災科学技術センター研究報告第30号

1983年3月とλに比例するものとする．一方σは，積雪の粘弾性的な性質を考慮して，λがある一定値λ。。よりも小さいときはその差に比例するように増加し，λ。。よりも大きいときはゼロであるとする．以上をまとめると
カ＝ωλ（9）1−／ろ（∵） （λ≦λ。。）（1o）（λ＞λ。。）となる，α，ろは比例定数である．

いま次のような無次元化をする．
εo1
λ＝A。。〃，0＝一ガ，！＝ε。λ。。戸，オ＝一オ＊ε1

ろ＊印のついた量はそれぞれの無次元量である，このときグライドを記述する基礎方程式は次のようになる．

扶一／■㍗1：1：1・∴：：1一（11）戸＝（1＋〆）λ＊

ただし

α＝oεo
ろε1

〆は〃の関数としてo。

＝

ブーλ＊（12）〃と表わされるから微分方程式は次のように〃だけで閉じた形となる．
V＊qホP＊

f＊
図2Fig．2

1A‡OA‡Of‡グ，ゲ，〆とパの関係
The

re1ationsofグtoλホ，ゲtoλ‡andパtoλ＊．

一194一κグライド現象の数理モデルー納口qホq

ホぴべOべO

鵬：＝沈Vヅ

図3−1

形態2

形態1

図3−2CaseFig．3−2

Fig．3−1！：snow

decelerates

tocoverstop

プCase2：snow COVer

g1idessteadily

ず岬・1・

庭1

肘肘対

O

I1・O

守Vざ

V＊V＊

図3−1
Fig．3−1

形態4

形態！

図3−4CaSe

Fig．3−4 Case 4：snow
acce1erates to！：snow

acce1erates

tocover
ava1anche（o＊＞パ）ordece1e・ates

tostop（パ＜パ）

一195一

cover
ava1anche

国立防災科学技術センター研究報告第30号

1983年3月

43㌧

42

2
3

1001

図4
Fig．4

234

パラメータ空間上の形態の分布
Diagramof

thedistribution・f・・…1−4

以1ヂー／■α（㌃川㍍）・／1∴lllll（13）（13）式の解のふるまいを幾何学的に考えてみる．が，σ＊，がは図2に示すような〃の関数である．〆に関して，λ＊＝Oで〆→o・になるという意味は静的な運動から動的な速度をもつ運動にかわることをさしており，その意味で〃＝0をなだれ発生点を呼ぶことにする．これに対して〃≧！＊ではo＊はゼロであるからA＊＝！＊を静止点と呼ぶことにする．このとき微分方程式の解の形態，いいかえるとグライド現象の形態はがと〆の関係により4種類に分類される．

形態1

はじめどんなグライド速度にあっても減速し静止する．はじめから静止している

場合は静止状態のままである．（図3−1）
形態2

安定な平衡点λ葦をもち，グライド速度がo姜よりも大きいときは減速し，ノ」、さいときは加速しやがて安定な定常グライド速度o姜になる．はじめから静止している場合は静止状態のままである．（図3−2）

形態3不安定な平衡点λ葦をもち，グライド速度がo姜よりも大きいときはますます加速されやがてなだれに至る．o姜よりも小さいときは減速して静止する．はじめから静止している場合は静止状態のままである．（図3−3）

形態4グライドが存在する場合はすべて加速されなだれに至る．はじめから静止している場合は静止状態のままである．（図3−4）

図4はパラメータ空問における形態の分布を示したものである．外的な要因の変化はパラ一196一グライド現象の数理モデルー納口

＊

＊V。

V。

α〉1〜：ム：忌：；1O111（Xf＊0

ト1＋・→・一図5111αf＊

ト1一ト・一1・一駆動カプと平衡速度o、‡の関係

Fig．5

The

relationofequilibrium

glide

ve1ocityぺtodriving

force戸

1515以；O．5以＝O，5○べ山べ廿11｛I1⑭〃り、1010‡㌧

支、㌧、・8㌧

ぐ多1ユ5f＊。1．5戸。O．95戸。1，2デ。095510

図6Fig．6

図7

グライド速度パの時問変化
G1ide

ve1ocityバvs．tinleオ＊

10

ゼゴ

時間ボにおける変位量パ

Fig．7

Displacementn／OtiOn

一197一VS．

戸duetinユeパ

tog1iding国立防災科学技術センター研究報告第30号

1983年3月

メータ空問の点の移動として表現される．この点の移動が形態1と形態2の領域内にあるかぎりグライドの加速が生じたとしてもなだれにはならない．なだれが起りうるのは形態3と形態4の場合だけである．しかも形態3の場合は初期条件として〃＜κが満たされることすなわちグライド速度がo姜を越えることがなだれ発生の条件であり，形態4の場一合は動きが存在するということがなだれ発生の条件となる．

定常グライド速度o姜とそのときのλ乏は次のとおりである．

。姜＝戸■1f14）

1一α戸

κ＝1一α1＊1一α（1l）

形態2の場合（！＊＜1／αかつ戸＞1）（14）式は安定に存在する定常グライド速度を表わし，形態3の場合（戸＞1／αかつ戸＜1）は加速してなだれに至るために瞬問的にもたなければならない最小の速度を表わす．戸とo姜の関係を定性的に図5に示す．α＜1のときに存在するo套は安定であり，駆動力戸の増大につれてそれに対応するグライド速度はどこまでも大きくなり，戸が1／αのところで形態4と隣接する．一方，α＞1において存在するo姜は不安定であり，戸の増大につれて値は小さくなる．すなわちなだれは起りやすくなり，戸が1を越えると運動しているものはすべて加速しなだれに至る．

図6，7は同じ初期条件を与えた場合についてグライド速度とグライドによる変位量パのろε1時間変化を計算したものである．ただし〆は実際のグライド量五の無次元量（〆＝一∫）εoである．∫＊二〇．95の曲線は形態1に属し，減速して静止する．戸＝1．2，1．5，1．8の曲線は2以＝0．5右控

α＝11以＝2O051001t‡

図8Fig．8

図9静止状態におけるパの時問変化

2V8Fig．9

速度0oホと平衡抵抗力R、‡の関係
Equi1ibriumgiVen

λ＊vs．timeガfor〆二〇一198一

resistant

Ve1OCityOo‡force

R、申vs．

グライド現象の数理モデルー納口

形態2に属し，はじめは減速あるいは加速するが，やがてそれぞれに応じた定常グライド速度になる．戸＝2．5の曲線は形態4に属し，加速してなだれに至る．

どの形態においても〃＞ブでは静止のままである．しかしその場合でもλについてはλ。。よりも小さいかぎり増加することになる．図8はこの〃の増加を示したもので，これはちょうど接触時問の長さに応じて積雪の粘弾性的な性質のために接触面積が増していくような現象に対応している．この場合時定数は1／ろである．

つぎに，このモデルで定速グライド試験をした場合の低抗力Rのふるまいについて考える．この場合，基礎方程式は次のようになる．

芝ぷ・｛1∴二二：1二／（16）○きは一定なので微分方程式はすでに〃だけで閉じた形になっている．この方程式は平衝点
1

κ＝1＋ωo非（17）をもち，しかも常に安定である．これに対応する抵抗力飛は1＋oさ

飛＝1＋αoo ホ（18）となる．平衡値に達するまでの過渡現象の時定数は，A＊≦1の場合（ろ十〃。）■1であり，＞1の場合1／ω・である．したがって定速グライド試験では過渡現象の時問スケールは眺が大きくなるほど小さくなる．

R；とo葛の関係を図9に示す．前にも指摘したとおり，これは戸とo姜の関係（図5）と同じである．αが1よりも大きい場合は抵抗力の平衡値硝は速度の増大につれて減少し，αが1より小さい場合は鴉は増大する．そして速度がますます大きくなるとR妻は1／αに近づく．

3．2例2
例1のモデルでは，定常グライド速度はいくらでも大きな値をとることができた．そこでつぎに，定常グライド速度の大きさがある一定値以上にはならないようなモデルを示す．

例1ではヵがλと〃に比例するとして考えた．そこでつぎに力がλにはよらずひだけに比例する場合を考える．すなわち（19）

カ＝ωである．ただしσとRについては例1と同じとする．このとき例1と同様の無次元化をしてグライドの基礎方程式を示すと次のようになる．

一199一

︵＃曽︑sV︷s︶き弓〜①揚ω①i＝b0さ︵㌔ぎ︿も︶丑o畠扁畠〇一ω8．〜﹄①司oo〜﹄①＞◎oきo冒ω一N①ω〜O
N糧津㈹−；．喜﹄N−oH国．︵廿曽︑§V§︶q◎揚◎一ω①一d﹄①司o①oF㌔︷．

−H撞酷．〜﹄2①oo〜﹄㊤＞◎o参◎冒ω一H①ω〜O

きぐ嵩ぎ︿も一①︷畠扁奏8ω8を◎uω8①竃O辿① O⊆市−d＞O
HOつm−〇一国

一逆H−oH国一−昌．靱竃

2o060−胃．碁﹄−

ちぢO一−︒つ撞亀kらh㌧

蟹o響≧O⊆ωFσ恒ちoF﹄①＞◎O

二一U

も■を

1983年3月
第30号
国立防災科学技術センター研究報告由司

5○勺oグライド現象の数理モデルー納口

2n＝O．5f‡。O．50＝O．5

戸＝1．05

㌧

㌧

0151
べ

図11形態1の場合のグライド速

図12

形態2の場合のグライド速度がの時問変化度〆の時問変化
Fig．11Glide

ve1ocity

パincase1o‡vs．time

Fig．12Glide
ボinvelocityパvs．time
case2

43
も2321OO

1234n

図13
Fig．1310t‡パラメータ空問上の形態の分布
Diagramof

thedistribution

一201一ofcases1−3国立防災科学技術センター研究報告第30号

1983年3月

紫一／一プ：ll：1：1・／∵：1：lll（20）。＝ 戸一〃〃〃αεo1

ただしβ＝・一ろε1 λ。。
この場合はグライドの形態は次の3種類に分類される．

形態1

不安定平衡点κが存在し，〆＞o姜、では加速しなだれに至る．o姜砒＞ガでは減速して静止する．はじめから静止している場合は静止状態のままである．（図10−1）

形態2

安定平衡点κ、と不安定平衡点λ妻αが同時に存在し〆＞o芸、では加速しなだれれに至る．o姜＞が＞o姜。では減速し定常グライド速度〃姜。になる．同様にo姜。＞〆では加速し定常グライド速度o姜。になる．はじめから静止している場合は静止状態のままである．（図10−2）

PイO座；

f斗府

ホVc

一一一一一一一一一一・

図14

構造不安定な形態

Fig．14

Transient casebetween case2andcase3．
V＊一202一グライド現象の数理モデルー納口V二O〈1T1川コ
六Φ＞

⊥、…

讐一」

f‡O

トー一1一一1・一ト・一
図15

β＜1のときの駆動カプと平衡速度ぺの関係．

Fig．15The

re1ationof

equi1ibrium

g1ide

veloci−ty・、＊t・d・i・i㎎f・…戸i・・・…fβ〈1．

2n＝113

ぜ1∩＝112

n：1n＝2OO

12Vま

図16Fig．16

速度o。＊と平衡抵抗力R、幸の関係．

Equilibrium

resistantgiVen乙b

VelOCityt

203一force

凡＃vs．

国立防災科学技術センター研究報告第30号

1983年3月

形態3動きがある場合はすべて加速しなだれに至る，はじめから静止している場合は静止状態のままである．（図10−3）

図1ユと図ユ2はそれぞれ形態1と形態2について初期条件をいろいろかえた場合のグライド速度の時問変化を言十算したものである．安定なグライド速度と不安定なグライド速度があるのがわかる．

パラメータ空問戸一β上での形態の分布は図13のとおりである．このモデルで特徴的なことは，力とσが2点で交わる場合があるということである．すなわち形態2のようにふたつの平衡点を持つ場合があるということである．これらの平衡点を式で示すと（素1）一β十1±平（、1）であり，このときのグライド速度は（簑1）一1一β手乎

（2、）である．図13における形態2と形態3の境界線は上式の根号内がゼロになるところである．
このときの形態は図14のとおりであり構造不安定な状態にある．いまβが1よりもノJ、さいときの！＊とo姜の関係を図15に示す．安定な定常グライド速度は戸が！以下では存在しな

22㌻

らO

1O20t‡
図17

グライド速度〆の時問変化
G1ide2ゼ図18（力＝ωλ，σ＝0）

Fig．171
グライド速度パの時間変化（カ＝ω，σ＝0）

velocity〆vs．time1‡（力＝〃A，σ＝O）

Fig．18Gユidevelocityパvs．timeヂ（力＝ω，σ二〇）一204一

グライド現象の数理モデルー納口

いが1よりも大きいときには戸の増大につれて大きくなりo姜、とo芝が一致した点すなわ

ち形態2と形態3の境界で消失している．これは，はじめ定常グライド運動していたものが戸のわずかの変化に対して突然加速を開始しなだれに至るような現象を表わしている．このときの臨界定常グライド速度o姜は。

1一βo。＝2β（23）である．

図16は，このモデルの場合の定速グライド試験の結果得られる低抗力の平衝値馬と速度○ざの関係を示したものである．

3．3

その他の例

例1，例2では

σ＝ろ（λ。。一λ）と表わして積雪の粘弾性的性質によるλの増加を表現してきたが，たとえば遠藤（1980）の笹の引き抜き低抗にみられるように，積雪中の笹が抵抗力の主なにない手である場合には，グライドにともなって笹は抜けて減少していく一方である．このような場合にはσはゼロとして扱うことになる．

いま例1においてσをゼロとした場合を考える．このときグライドの方程式はつぎのとおりである．

（24）

ただしλ一！λ・，。一旦が，チーε1ヂε0ε1αεO

図17はこの場合のグライド速度の時問変化を計算したものである．同様に図18は，例2の場合でσをゼロとしたのを計算したものである．ただしこの場合の時間スケーノレはチニ（ε1ルε蕎）・ヂである．これらの場合は当然のことながら定常グライド速度はなく，静止状態か加速かのいずれかである．

4．

あとがき

ここで提案したグライドのモデルの考え方は，ちょうど物質の粘弾性的性質を表現するための力学的モデルの考え方に対応している．ダッシュポットやスプリングのいろいろな結合によってさまざまな粘弾性を表現するように，積雪と地表面の問の接触の機構のちがいによ一205一

国立防災科学技術センター研究報告第30号！983年3月

るグライド現象の性質や底面抵抗力の挙動を力，σ，沢の選択によって表現する方法を示したものである．したがってグライド現象を定量的にモデル化しようとする場合，まず現象を定性的に表現できる力，σ，灰を選び，それからグライドや定速グライド試験のような実際現象と対応させて定量化することになる．

これまで，なだれの発生と関連づけて斜面積雪の安定不安定を判断する場合，理論的には力のバランスによる表現のものがほとんどであった．しかし，ここで述べたグライドモデルのようにグライド現象を微分方程式の平衡点の安定性あるいは微分方程式自体の構造安定性と結びつけることによってグライドによる斜面積雪の安定性の表現方法が得られたことになる．なおここで用いたような真の接触面積による底面抵抗の（8）式の表現は，摩擦凝着説による摩擦低抗の数理モデル的な表現と考えることもできる．

ここでは斜面積雪の移動を地表面と積雪の関係だけでみてきたが，これは積雪の動きが均質の場合だけであって動きが不均質でクリープの効果などがきいてくる場合についての考察や実際の現象への応用については別報でおこなう．

本論文の校閲をしていただいた山田穣研究室長に感謝いたします．

参考1）

文献

秋田谷英次（1974）：斜面積雪の挙動の研究III一接地面での移動（グライド）．低温科学，物理篇，32，98−104．

2）

秋田谷英次（1975）：斜面積雪の挙動の研究V一斜面積雪のグライドとクラックの発生，1974〜1975冬．低温科学，物理篇，33，103−108．
3） 遠藤八十一（1980）：笹地斜面における積雪のグライド機構II．低温科学，物理篇，89，81−89．
4） 中村勉他（1972）：かや張り実験斜面で発生した小なだれ．昭和47年度日本雪永学会予稿集．
5）

渡辺興亜・五十嵐高志・山田穣（1978）：1976−1977年冬期の新潟県を中心とする地方の広域積雪現象について．国立防災科学技術センター研究速報，第29号．
（1982年1！月29日原稿受理）

一206一