水工学論文集,第57巻,2013年2月
土木学会論文集B1(水工学) Vol.69,
No.4, I̲1147-I̲1152, 2013.

砂州性蛇行の形成に関する線形安定解析
LINEAR STABILITY ANALYSIS ON THE MEANDER
FORMATION ORIGINATED BY ALTERNATE BARS

島田 龍市1・清水康行2・長谷川和義3・伊賀久晃4
Ryuichi Shimada, Yasuyuki Shimizu, Kazuyoshi Hasegawa and Hisaaki Iga
1学生会員 北海道大学工学院環境フィールド工学部門（〒060-8628 北海道札幌市北区北13条西8丁目）
2正会員 工博 北海道大学工学院教授 （〒060-8628 北海道札幌市北区北13条西8丁目）
3フェロー会員

工博 株式会社北開水工コンサルタント 先端技術開発センター（〒060-0005 北海道札幌市中央区北5条西6丁目2番地）

4正会員 株式会社北開水工コンサルタント 札幌支店（〒062-0052 北海道札幌市豊平区月寒東2条20丁目5番10号）

There are two theories about the origin of river meanders: one is based on the bed-instability theory
and another one is on the bend-instability theory. According to the bed-instability theory, the meanders
will be developed due to alternate bars on the bed. Although this theory was occasionally denied in thebend-instability theory because of reasons that typical wave number of developed alternate bars falls inthe stable range of bend instability and so on, meanders induced by alternate bars are often seen in naturalrivers in Japan. The theoretical verification of this kind of meander generation is required.
In this study, the linear stability analysis is applied to the straight channel with erodible banks in orderto obtain the characteristics of the meanders. The results indicate that meanders induced by alternate bars
will exist and they have longer wavelength than that from previous bar-theory. And, it was found that
Froude number is the one of the parameter that affects the wavelength of these meander types.
Key Words : meanders, alternate bars , bank erosion, linear stability analysis, bend-instability,
bed-instability

理論研究が始められた4)．流れモデルや流砂量の取り扱いに違いがあるが，多くが線形化方程式系における河床不安定を問題として発生条件や卓越波長を導いている．
近年，自然環境の再生が多方面で試みられるようにな
沖積蛇行の成因に関するもう一つの説は，Engelund5)のり，川の蛇行復元もその中の主要課題になっている．ま蛇行流路内平衡底面形状に関する研究に影響を受けて始た，昨年の音更川のように，流路蛇行による大規模河岸
まった池田・日野・吉川6)，Ikeda Parker Sawai7)の平面不浸食の発生が一応の整備を終えた河川の安全を脅かす事
安定説である．この研究は，流路の初期状態に含まれる
例1)も起こっている．これらを背景として，蛇行関連研
微小な湾曲が主流の平面的な偏りを生みだし，河岸を浸究があらためて注目されている．
食することによって湾曲が成長すると考えるもので，そ今日までの沖積蛇行に関する研究は膨大な数にのぼる. の後Blondeaux Seminara8)，Tubino Seminara9)などのイタ
特に成因に関する研究は川がなぜ曲がるのかという問い
リア研究者によってさまざまな発展が試みられた．このに答えるものであり，蛇行流路の処理法を考える上でもうち，Blondeaux Seminara8)の研究では，交互砂州の卓越
重要な視点を提供する．
波長と平面不安定の卓越蛇行波長の比較がなされており，
現在，沖積蛇行の成因に関して二つの有力な説がある． 後者が前者の３倍ほど大きく実測値に近いこと，この波
その１つは，木下2)3)によって提唱された交互砂州（砂礫長に対応する砂州の成長率がほぼゼロ（中立不安定）で堆）成因説であり，直線水路内おいて発生する交互砂州あることを示した．彼らはこの現象を一種の共振現象
が水流の蛇行を喚起し，さらに流路の蛇行を呼び起こす（レゾナンス）とみなしている．これらの成果をもとにと考えるものである．この考えは蛇行研究に非常に大きSeminara10)は，湾曲機構が河川蛇行のさまざまな性質をな影響を与え，1960年代以降世界中で交互砂州に関する説明できる唯一の機構であると結論づけ，交互砂州は流１． はじめに

I_1147

(a) 平面図

図-1 対象水路概要(b) 断面図

路蛇行の主要な成因ではないとしている．しかしながら， る．これらの主張をめぐっては，蛇行の定義も関係する．
本邦では交互砂州を契機とした蛇行や河岸浸食が各所に前者は水流の蛇行も視野に入れた水衝部に重点を置く見見られ，砂州性蛇行の存在は否定できない．砂州から流方であるのに対し，後者は河道全体の湾曲の連なりを問路蛇行への発達を扱った解析研究は非常に少ないが，そ題にしている．両者はいずれも存在すると考えるべきで11)の一つにMosselman があり，側岸浸食を伴う河床不安あろう．
定解析に先鞭をつけた．ただし，粘着性河岸を対象としているほか，定常解析にとどまるなど制約が多い．
３．定式化
本研究は，砂礫河岸からなる直線状水路において河床砂州の発達が河岸の蛇行を引き起こす場合を対象として
線形安定解析を試みたものである．従来の交互砂州解析(1) 基礎式
は，河床起伏の横断方向表現に両岸固定壁を境界とする
本研究では初期状態として平面的に直線かつ側岸に傾固有関数を仮定してきたため，砂州に伴う水流蛇行を河斜を有する開水路を対象に安定解析を行う．図-1に対象
道蛇行に直接結び付けることができなかった．本研究で水路概要および座標系を示す．
は河床起伏を固有関数表示せず，横断方向流速の両岸境
ここで河床変動および側岸の移動は流れの変化と比較
界条件を運動学的条件とすることによって水流蛇行が河し十分に遅いと仮定し，流れの方程式には以下のデカル
道平面形状の変形に反映するようにしている．ただし，
ト座標上の準定常浅水流方程式を用いる．
これらを表現する方程式系が斉次式となるため，直接河
2
2
z *
u *
u *
h*
g*
岸形状を求めることはできないが, 砂州性蛇行解析の第u * *  v* *  g * b*  g * *  2
u *  v* u *  0 (1)
*
*
x
y
x
x C h１歩になるものと考える．なお，線形解析は河床地形等
*
*
*
2
2
z
v
v
h*
g*
に関し初期直線等流状態からの微小な周期変動の発達を
u * *  v* *  g * b*  g * *  2
u *  v* v*  0 (2)xyyy C * h*扱うものであり，その成長率の極大となる波長がその後u *h* v*h*(3)の現象を支配するものと想定している．その妥当性の検

0
x*
y *
証は波長などの実測値との比較による他はなく，本論文においても実験値との比較がなされている.
ここでx*,y*: 流下方向および横断方向座標, u*,v*: x*, y*方向の流速，zb*: 河床高，h*: 水深，g*: 重力加速度, C*:
Chezy係数を表す．また，添え字*は有次元量を表す．
２．砂州性蛇行の概要
また，流砂平衡式については空隙を含む見かけの流砂量を取り扱う以下の式を用いる．
*(4)木下2)3)は，交互砂州の発生・発達、流路蛇行への移行
zb* S x* S y


0
*
*
*
を体系的にまとめ，向きの異なる1組の砂州によって１t
xy
蛇行する「蛇曲」と，1蛇行中3以上の砂州によって特徴
ここでt*: 時間，Sx*, Sy*: 河床における空隙を考慮した流づけられる「迂曲」の存在を示している．概して河床勾下方向および横断方向流砂量を表す．
配がきつく河床材料が粗い本邦河川では，これらの形態
流砂は掃流砂のみを対象とし，流下方向にはMeyerはごく一般的であり，砂州の発生・発達・減衰は河川改Peter・Müller式を，横断方向にはMosselman11)による横断
修にとって重要な視点となっている．本研究ではこの種
方向平均流，流線曲率効果，横断方向河床勾配を考慮しの河道蛇行を「砂州性蛇行」と呼ぶことにする．
た以下の式を用いる．
3sg * d *310)
S x*  8 *   *c  2(5)Seminara は蛇行研究レビューの中で，①交互砂州は，
1 p
蛇行を形成するには移動速度が大きすぎること．②通常
*
*
*
*
v
 *c zb h v
S y*  S x*  *  A * * (6)の交互砂州の無次元波数（πB*/L*， B*=川幅，L*=砂州
u
u

x

k * y * s波長）は0.5ほどであり，湾曲不安定の安定領域に入る
ここでτ*: 無次元掃流力，τ*c: 無次元限界掃流力，s:こと．③交互砂州は弱蛇行流路において共存移動し，砂砂の水中比重，d*: 粒径，λp: 砂の空隙率，A: 2次流影州が河川蛇行の固定砂州に変化するというアイデアと矛
響の重み係数(=7～11)，μs: 静止摩擦係数(=0.6)，μk:
盾することなどを指摘して，交互砂州成因説を退けてい動摩擦係数(=1.0)を表す．

I_1148

(2) 側岸浸食モデル
本研究では側岸浸食による河床への土砂供給について
斜面掃流浸食モデルを用いる．斜面上の掃流砂について
は平野12)や長谷川13)を参考とした次式を用いる．
S

*
x



 8



(4) 基礎式の線形化及び無次元化
平坦床等流状態を基本状態とし，以下のように摂動展開を行う．

u , v , h , z , B , B    u ,0, h , z , B 2 , B 2  u " , v " , h " , z " , B ' , B ' (17)*(7)

ここで S*xφ: 斜面上における見かけの流下方向流砂量，
φ: 側岸の傾斜角，ε: 側岸におけるせん断力の補正係数(=0.6～0.8)，τ*cφ: 斜面上の無次元限界掃流力を表す．
τ*cφは以下のように表すことができる．

 tan   s 

2 *c   *c cos  1  

(8)

掃流浸食モデルでは，側岸浸食による河床への土砂の供給は河床上の横断方向流砂量と類似した挙動をとると仮定し，以下のように表す．
 v*
 *c 
h * v **S bankR
 S x*  *A * *
 tan 
u *


u

x
*
s  k  *R y  BRy  BR(9) v*
 *c 
h * v **S bankL S x*  *A * *
 tan u * k *L 
u

x
*
s
y   BLy   BL


(10)

*
R

*
L

*
0

*





*
SbankRzb*  tan  H *  t *y  BR* 
y*  B*R 

(11)



*
SbankL
zb*


tan




H*
t *

* 
y   BL 
y*  B*
L 

(12)

*

*
0

*

BL
t *

BR
 v* * *y  BRx* y *  B *

(13)

BL
 v* * *y   BLx* y *   B *

(14)

*

y *  B R*

R

*

*

 u*
y *   B L*S*
 S bankR
*

(15)

*
 S bankL(16)*y y *   B*

*
R

*
L

*
b

*

*
0

b

2
z b "u" v" A  2 v"
1  zb "
 G0



0
t
x y  xy f  *0  y 2

(19)

u" 1 z b " h" 2C f u"C f h"  0
x F02
x(20)

v" 1 z b " h"C f u"  0
x F02y(21)

u" v" h"0x y x(22)

BR " G1 tan 
1 zb "
u"t
H f  *0  y BR  tan  t y  BR

(23)

BL " G1 tan1 zb "
u"t
H f  *0  y  BL  tan t y  BL

(24)

3 *0
 *0  cos  *c   tan   
 G0 
, G1   G0
1  
 *0   *c
 *0   *c    s  


u0*
B0*
g*
F0 
,   * , Cf  2
h0
C*g * h0*34

 s  k *0
 *c

(25)
(26)
(27)

(5) 線形安定解析
摂動項が以下に示すフーリエ級数で表されるとする.

u", v", h", zb "  u, v, h, zb exp ikx  ct  c.c. (28)
ここでk: 摂動の無次元波数(=2πB*/L*)，L*: 波長，c: 摂動の複素角周波数, i: 虚数単位を表し，c.c.は共役複素数を表す．またu，v，h，zbはyのみの関数である．これら
を式(19)，(20)，(21)，(22)に代入すると以下が得られる．

y  BR**b

ここで，係数は以下のように表される．

L

S *y

*

2

(3) 境界条件
右岸，左岸において①横断方向流速が水面波の運動学的条件と同様な規定を受ける，②側岸における横断方向
河床流砂量は側岸からの土砂供給量に等しいという合計4つを用いる．これより，
 u*

*

*
0



f  *0  

*

*
0

*
b0

(18)
B * h*
x* , y * , BR* ' , BL* '  B0* x, y, BR ' , BL ', t *  0* 0 t
S x0
ここで*のついていないものは無次元量を表す．式(17)，(18)を基礎式に代入し，摂動項を抽出すると以下の6つの無次元化された線形化基礎式が得られる．

ここで，H*: 天端高を表す．

BR
t *

*
0

u ", v "  u u", v", h ", z ", H   h h", z ", H 

び左岸の座標値を表す．また，添え字R及びLはそれぞ
れ右岸，左岸の値であることを表す．
また，本研究では川幅変化に関する支配方程式として掃流浸食モデルを組み込むParkerら15)の式において，曲率項を0としたものを用いる．

BL*
1

t * tan   zb* y * **b

ここで添え字0は基本状態を表し，プライムは摂動量であり，ダブルプライムはx*, y*, t*の関数，シングルプライ
ムはx*, t*の関数であること表す．また，基本状態においては側岸からの土砂供給はないものと仮定する．さらに，
以下を用いて基礎式の無次元化も行う．

ここでS*bank: 側岸からの土砂供給量，BR*, BL*: 右岸及

BR*
1

t * tan   zb* y * *

*

*

32

2
2
3

  *   tan  
sg *d *    1  *c 3 2
1  p  *c cos     s *L

I_1149d a1 a2'
dy

 a5
0

a10
 0

d

 a13 a14 dy0a  a  dyd 

a7
' d
a11
dy

a8
' d
a11
dy

a13

0

2

'
4

''
4

2

 u  
 v   0
 h  
 zb 
(29)

a1  ikG0 , a2 ' 1 ikA

, a4 '  ikc , a4 "  

1

T1 

f  *0 

a5  ik  2C f , a7 ikik C f , a8  2
F02
F0

a10  ik  C f , a11' 1, a13  ik , a14  1
F02

(30)

整理すると，以下に示すvに関する4階線形微分方程式が
得られる．
d 4v
d 2v

A
 A2v  0
(31)
1
dy 4
dy 2
A1 

1
a11 ' 2 a 4 ' a5  a1 a 7  a1 a8 
a11 ' a 4 ' a 5

 a11 ' a 2 ' a 8 a 5  a 7  a 8   a10 a 8 a 4 " a 7  a 5 A2 

a8 a10a4 ' a7  a5   a1a8 
a11' a4 ' a5

(32)
(33)

これらから，v，zb，uの一般解を以下のように得る．
1 yv  C1e

 C2 e 1 y



 C3e

2 y

 C4 e 2 y

(34)





A 
A 
 y
  y
 y
  y
z b   A3 1  4  C1 e 1  C 2 e 1   A3  2  4  C 3 e 2  C 4 e 2


1 
 2 u   L1 C1e

1 y

 C2 e 1 y L C e
2

2 y

3

 C4 e 2 y (36)

11
1    A1  A12  4 A2  , 2    A1  A12  4 A2 
2



2

a5
a7  a5 a8 a10
A3 , A4 
a8 a11 ' a8  a7  a5 
a8 a11 ' a8  a7  a5 

Lj 

a2 '
a '
A  a"
 j  4  A3  j  4   4 A3  j  A4   j ,a1a1  j  a1
 (35)(37)(38)

 j  1,2 (39)

ここでC1, C2, C3, C4は未定定数である．式(13)，(14)，(15)，(16)の境界条件を適用し，側岸が移動しても川幅の変化はない条件（BR+BL=1）が維持されるとするとC1,
C2, C3, C4に関する以下の方程式が得られる．
 D1

 D5
D
 9
D
 13

D2

D3

D6

D7

D10
D14

D11
D15

D4  C1e 1 2 
D8  C2 e  1 2 0D12  C3e 2 2 
D16  C4 e  2 2 (40)

A  GT L 
1
D1  1  T1 1  c  A3 1  4   1 2 1    

H 
c
1 



A  GT L 
1
D2  1  T1 1  c  A3 1  4   1 2 1    

H 
c
1 


A  GT L 
1
D3  1  T1 1  c  A3  2  4   1 2 2    H c 2 A  GT L 
1
D4  1  T1 1  c  A3  2  4   1 2 2    


H
c


2 


D5  D2 e 1,

D6  D1e

1,

D7  D4 e 2,

D8  D3e2

 1

D15  e

 2

(42)

式(40)が有意な解を持つためには，係数行列のデターミ
ナントが0となる条件が必要となる．これにより摂動の複素角周波数cを得ることができる．式が複雑であり解析的に解くのが困難であるので，数値計算にて解を得る
必要がある．cの実部Re[c]は摂動の移動速度を，無次元波数kとの積の虚部Im[kc]は摂動の増幅率を表す．

４．解析結果と考察(1) 実験結果・数値計算結果との比較
本安定解析結果の妥当性を検討するために実験結果との比較を行う．比較対象として，Akahoriら14)の複断面直線水路の蛇行実験の結果を用いる．実験は浸食性河岸を有する長さ50m，初期幅40cm，全幅90cmの直線複断面
水路で行われ， 40時間通水した．河床材料粒径は
d*=0.76mmを用いている．表-1に40時間後の水理量の測定値を示す．図-2は時間経過ごとの河床高の測定値を示している．
また，実験では幅の拡幅に制約があるため，初期条件
を実験と同様にした実験的数値計算を行い，実験よりも大きな振幅の蛇行を発生させた．表-2に数値計算の最終結果を，図-3には数値計算における水深の時間変化を示
す．数値計算には河川解析ソフトIRIC2.016)に搭載されている平面二次元の河床変動モデルを用いた．計算格子は
主流方向，横断方向それぞれ491×61分割とした．
実験結果及び数値計算結果のいずれも交互砂州の発達から流路の蛇行へと発展しており，砂州性蛇行の存在を示唆しているものと考えられる．これらの最終結果をもとに本研究で得られた安定解析を試みると図-4，図-5のような結果が得られる．まず実験についての不安定解析
結果では，成長率のピークが無次元波数k=0.45で現れ，
これにより卓越波長を計算するとL*=12.5mとなる．一方，観測された蛇行波長16.0mである．一般に砂州現象は大きなばらつきを示すものであり，その意味では近い値であると判断できる．また，数値計算についての解析結果
は，成長率ピークがk=0.325で現れ，卓越波長を計算す
るとL*=10.4mとなり，数値計算で得られた波長9.3mと近い値となる．したがって，実験，数値計算のいずれの場
合の蛇行波長とも良好に合致する結果となり，安定解析の結果が十分妥当であるといえる．

(41)

G1 L1 tan
G1 L1 tanD9  
  A3 1  A4 , D10   A3 1  A4 
G1 L2 tan
G1 L2 tanD11    A3  2  A4 , D12   A3  2  A4 

D13  e

S*
ik
tan 
, T2 ,   *x 0*
 tan f  *0 
u 0 h0 G1 L1 tan   G L tan   A3 1  A4  , D14  e 1  1 1
  A3 1  A4 








 G1L2 tan  G L tan   A3 2  A4  , D16  e 2  1 2
  A32  A4 








(2) 先行理論との比較
図-6はBlondeux Seminara8)による平面不安定理論および交互砂州理論を用いて前述の実験及び数値計算と同条
件で解析行ったものである．この理論では無次元波数
は半幅で無次元化されている点に注意されたい．
実験(図-6(a)）について平面不安定理論および交互砂州

I_1150表-1 実験終了時(40時間後)の水理量

表-2 数値計算による結果(14時間後)

250[cm]100[cm]250[cm]100[cm]

図-2 Akahori14)らの実験における各時間の河床高

図-3 数値計算による各時間の水深図-4 実験データによる成長率の解析結果

図-5 数値計算データによる成長率の解析結果(a) 実験データによる解析結果(b) 数値計算データによる解析結果
図-6 Blondeux Seminara8)による解析結果
理論で卓越波長を計算すると，それぞれ14.0m，4.0m とい。また，平面不安定理論による卓越波長も18.8mと計なる．平面不安定理論による卓越波長は実測波長に近く， 算波長の2倍程度大きい．一方，本安定解析理論で得ら
交互砂州理論の中立安定域近くに現れている．交互砂州れた卓越波長はより計算値に近い値を示しており，砂州
卓越波長は小さめである．本研究の卓越波長をこれらと性蛇行ゆえの結果を得たと解釈することができる．しか比較すると，波長はBlondeux Seminaraの交互砂州波長のし，多くのデータをもとにさらに解析する必要がある．
3倍ほどであり，平面不安定理論の波長に近い値になっている．Blondeuxらは，側岸固定の交互砂州発生理論に(3) パラメータの変動特性
よる波長が平面不安定理論及び実測の波長の1/3ほどに本研究における不安定解析モデルを用いて，各種水理
しかならないことを主張し，砂州性蛇行成因説を否定し量を変化させその挙動を調べる．図-7は川幅水深比，ている．しかし，本理論によれば側岸浸食を伴う交互砂
Froude数F0，無次元掃流力*0を変化させ，その成長率を
州波長は十分大きくなることができ，蛇行現象が平面不示したものである．
安定の機構に引き継がれ得ることを示している．
図-7，(a)より川幅水深比の変化では蛇行波長，成長率
数値計算データによる結果(図-6(b))では，砂州成長率の大きさともに大きな変化は見られない．Blondeux
が負になっておりシミュレーション結果に対応していな
Seminara8)の砂州成長率曲線と比較しても卓越波長に大I_1151

(a) 川幅水深比による変動特性(b) Froude数による変動特性(c) 無次元掃流力による変動特性

図-7 各種パラメータ変化に対する成長率の変化(粒径水深比=0.005, 側岸傾斜角=30[deg], 天端高水深比=1.0)
きな変化は見られないことから，側岸移動を考慮した場
合でも卓越波長は固定壁と変わらないと考えられる．
図-7，(b)について本安定解析モデルではFroude数の増大に伴い，波長が増大する傾向が見られた．この結果はBlondeux Seminara8)による交互砂州理論と類似した傾向である．
図-7，(c)について無次元掃流力は成長率の振幅に影響するものの，卓越波長については変化が見られない．無次元掃流力が増大しても洗掘が大きくなるのみで，蛇行
波長の変化は生じないものと考えられる．また，限界掃流力に近くなると，数学的性質上，無限大に発散することから，限界掃流力を用いない浸食モデルを用いる必要性があるといえる．

2)木下良作：河床における砂礫堆の形成について―蛇行実態の一観察―，土木学会論文集，42号，1-21，1957
3)木下良作：石狩川河道変遷調査，科学技術庁資源局資料第36号，1961
4)土木学会水理委員会研究小委員会（代表 芦田和男）：洪水流の三次元流況と流路形態に関する研究，3，4，58-118，
1982
5)Engelund, F. : Flow and bed topography in channel bends, J. of Hyd.
Div., ASCE, Vol.100, No.HY11, 1974
6)池田駿介・日野幹雄・吉川秀夫：河川の自由蛇行に関する理論的研究，土木学会論文報告集，第255号，1976
7)Ikeda, S., G. Parker and K. Sawai : Bend theory of river meanders,
Part1. Linear development, J. of Fluid. Mech., Vol.112, 19818)Blondeaux, P. and G. Seminara : A unified bar-bend theory of rivermeanders, J. of Fluid Mech., Vol.157, 449-470, 1985

５．結論

9)Tubino, M. and G. Seminara : Free-forced interactions in developing
meanders and suppression of free bars, J. of Fluid Mech., Vol.214,

本研究の不安定解析結果から，側岸浸食を伴う交互砂州理論は先行研究による交互砂州理論より波長が伸びる傾向が見られた．この結果は，波長の短い交互砂州から長い波長を持つ平面不安定蛇行への移行が可能であることを示唆する．
また，数値シミュレーションの発達した状態での結果
は，本安定解析ではその特性を解釈できるがBlondeux
Seminaraの平面不安定理論では説明できない. これは,砂州性蛇行の独自の発達を裏付けている.
さらに，本研究における安定解析モデルでは，卓越波
長変化に大きく影響するパラメータはFroude数であるといえる．しかし，この傾向は実現象において必ずしも成り立つとは言えず，この点は検討を要する．

131-159, 1990
10)Seminara, G. : Meanders, J. of Fluid Mech., Vol.554, 271-297, 2006
11)Mosselman, E.: Theoretical investigation on discharge induced
river-bank erosion, Communications on hydraulic and geotechnicalengineering, Delft University of Technology, 1989
12)平野宗夫：拡幅を伴う流路の変動について，土木学会論文
報告集，第210号，1973
13)長谷川和義：非平衡性を考慮した側岸浸食量式に関する研究，土木学会論文報告集，第316号，1981
14) Akahori, R., Yamaguchi, S., Kimura, I., Iwasaki, T., Shimizu, Y.
and Hasegawa, K.: Evolution of channel bed topography in a 50meter laboratory flume under various hydraulic conditions of bar
instability and resonance, River, Coastal and Estuarine
Morphodynamics 2011(CD-ROM), pp. 1998-2009, 2011

参考文献

15) Parker G. Shimizu Y., Willkerson G.V., Eke E.C., Lauer J.W.,

1)小田島大祐・桑村貴志・永田朋紀：音更川における出水時のPaola C., Dietrich W.E. and Voller V.R.: A new framework for

堤 防 一 部 流 出 の 原 因 分 析 に つ い て ，modeling the migration of meandering rivers, EARTH SURFACE

http://www.hkd.mlit.go.jp/topics/gijyutu/giken/h23_pre_intra/pdf_fil
es/aa/aa-16.pdf，2011

AND LANDFORMS, pp. 70-86, 201116) http://i-ric.org/ja/index.html（2012.9.30受付）

I_1152