別紙―２
平成26年度度

物理環境物
境要素の計算算結果をを用いた
たPHAB
BSIMにによるシシロザケ
ケ産卵
卵環境のの評価にについいて（独）土木研究所
所 寒地土木木研究所 水環境保全チ
水チーム○矢野 雅昭
渡邉 和好
平井 康幸

本研
研究では、河河床変動計算
算、浸透流計算算により、シシロザケの産産卵環境に関わわる物理環境境値でHABSIMによ る産卵環境のある流流速、水深、平均粒径、浸透流を予測測し、その値値を用いてPHの評価を行っった。数値計
計算結果は現現地の地形、河河床材料、浸浸透流の傾向向を概ね再現すすることがででき、産卵適適地についてても概略ではあるが評価すすることがでできた。
キーワーード：産卵床床、シロザケ、、PHABSIM、数
数値計算

1． はじめめに
北海道におおけるシロザ
ザケの漁獲高はは、全国の約88割を占め、地域ににとって重要なな水産資源となっている。シロザケの漁獲高
高は人工孵化放放流事業と伴伴に増加してき
きたが1)、近年でははシロザケ野生生魚の価値にもも着目され2)、河道改修においてても、産卵環境境に配慮する必要性があるると考えられる。
既往研究ででは、サケ科魚魚類の産卵環境として河床床材料
3）が重要である
ることや シロロザケが産卵環環境として浸浸透流の湧出する箇箇所を利用してていることがが知られている
る4）。
図-1 対象象個所
そのため、河河道改修が産卵卵環境に及ぼぼす影響を評価
価するには、これららの物理環境にに加え、魚類類の遊泳に影響響する流速、水深のの変化も把握し、産卵環境境として評価すする必本研究の対象河
本
河川は、毎年シシロザケの産卵卵が確認されれ要がある。
る北北海道の豊平川川とした。豊平平川は石狩川の一次支川でで、BSIM（Physicaal Habitat Simuulation
著者らは、過去にPHAB流路路延長72.5km、流域面積902kkm2の一級河川川である。解解Model）5)を用用いて、平面的的なシロザケ産
産卵環境の評
評価を
析はは、石狩川とのの合流点をKPP0.0（KP;上流方方向を正、単単m）として、図行ったが、評評価に必要とななる多項目の物理環境値をを、現位km
図-1に示すKPP11.8～13.00の区間を対象と地調査によりり面的に取得しており、調査量が膨大ととなっしたた。なお河床勾
勾配は約1/420で幅は60～70m程
程で、低水路幅た6)。また、河道改修などどにより将来変変化する産卵卵環境度でである。の評価を行う
うには、何らかかの方法で物理環境値の空空間分(2) 解析
布を予測するることが必要となる。これらの問題を解解決す
解る手段としてて、数値計算にによる推定と、それを用いいて産本研究で産卵環
本環境の評価に用用いる物理環境境要素は、流流卵環境を評価価することが考考えられる。
速、水深、河床材材料の平均粒径径、浸透流とした。そのたため、豊平川の対象
本研究ではは、数値計算にによりシロザケの産卵環境
境に影象区間で平成118年に調査されれた河床の状状SIMに態を響する物理環環境値を予測測し、それを用用いてPHABSIを初期条件としし、平成19年44月～平成23年年9月の期間ののより産卵環境境の評価を試みみた。そして、現地調査結結果と河床床変動計算を行行い、計算結果
果の河床地形と水位分布かか比較することとにより、そのの適応性を検討したものでである。 ら浸浸透流計算を行行った。豊平川
川のシロザケには、前・後後期群群があり、産卵卵環境に適するる湧水条件（伏伏流水由来、地下下水由来）が異異なることが指指摘されていいる7)。そのたた2． 方法
め、伏流した河川川水由来や地下下水由来などの浸透路長のの(1) 対象個所
所
いで、傾向が異異なる、水温や温や溶存酸素8)も検討に用いもい違い
Masaaki Yanno、 Kazuyoshi Watanab
be、 Yasuyukki Hirai

2007/5/17 24時
2007/5/22 22時
2007/5/26 11時
2008/5/20 22時
2010/5/25 10時
2010/5/25 23時
2010/12/ 3 22時
2011/9/ 6 3時
2011/9/ 6 16時
2011/9/ 7 5時
2011/9/ 7 18時
2011/9/ 8 7時
2011/9/ 8 20時
2011/9/ 9 9時
2011/9/ 9 22時
2011/9/10 11時
2011/9/10 24時2011/9/11 13時
2011/9/12 2時
2011/9/12 15時
2011/9/13 4時
2011/9/13 17時
2011/9/14 6時
2011/9/14 19時
2011/9/15 8時
2011/9/15 21時
2011/9/16 10時
2011/9/16 23時
2011/9/17 12時
2011/9/18 1時
2011/9/18 14時
2011/9/19 3時
2011/9/19 16時

流量(m3/s）

通過質量百分率（％）

ることが望ましい。しかし本稿では、河床微地形に起因流からの浸透流の流入は平均河床勾配(1/420)と平均透水した浸透路長の短い伏流水の湧出個所を好む前期群のみ係数(57.9m/day)の積である0.14m/dayが流入し、下流端か対象とすることとし7)、水温と溶存酸素は検討に用いなら同一の流速が流出している条件とした。底面の境界条かった。
件は、流出入がない条件とした。これらの浸透流計算の
9)
物理環境の計算結果は、GIS に取り込み、セルサイズ条件を表-2に示す。
2×2mのメッシュデータとした。そして、後述する平成
産卵環境の評価では、河床面での浸透流の湧出・浸透23年に行った現地調査結果と比較し、その精度についての傾向が重要であることから、計算結果の内、河床表面検討した。の鉛直方向成分の浸透流の実流速を用いた。なお、浸透a) 河床変動計算
表-1 河床変動計算条件（iRIC NAYS2D ソルバー10））
流速、水深、平均粒径の予測値は、混合粒径の河床変項目
計算条件
動計算することにより把握した。河床変動計算にはソルバー
Nays2DH10)
iRIC2.3のNays2Dソルバー を用いて、表-1に示す条件でソルバー・タイプ
河床変動計算有効計算を行った。計算格子は、平成18年に国土交通省北海
移流項の差分法
風上差分
河床材料の種類
混合粒径
道開発局札幌開発建設部が、縦断間隔200m間隔で実施
交換層厚0.30m、堆積層厚0.15m
した河川定期横断測量結果から作成した。河床材料の粒
流砂の種類
掃流砂&浮遊砂度分布は、同建設部がKP12~14の区間で1km間隔に左・乱流モデル
ゼロ方程式
右岸・流心の3点で実施した河床材料調査結果の平均値
樹木あり植生の遮断面積0.1を算出し、その値を図-2のとおり近似して用いた。計算
樹木の抵抗係数1.2
境界条件
下流端水位
等流水深条件では、平成20年に国土地理院で撮影された写真11)で格子設定
格子サイズ
約4×4m
樹木が確認される範囲に植生の抵抗を設定した。植生にマニングの粗度係数
河床 0.03012)
河岸護岸 0.025
よる遮断面積と樹木の抵抗係数は、内田ら が豊平川を100実測 平均値（KP12‐14）
対象に河畔林の影響を検討した際の値を用いた（表-1）。
90
計算 入力値80
計算に用いた河川流量は、国土交通省北海道開発局札7060
幌開発建設部の雁来水位流量観測所における観測値13)を5040用いた。また、計算負荷を軽減するため、図-3に示すと3020
おり平成19年4月から平成23年9月までの融雪出水規模程
1030度の流量である200m /s以上を計算に用いた。ただし平成0.111010023年9月の1,000m3/sを超える出水のみは減水期の影響を粒径(mm)
図-2 河床変動計算に用いた粒度分布考慮するため、平水流量の14.4m3/sに減衰するまでの期1200間も含めた。産卵環境の評価に用いた流速、水深の計算
1000
値は、計算により得られた河床地形に14.4m3/sの流量が800600流れている状態の値とした。
400200b) 浸透流計算
0
浸透流の分布を把握するため、浸透流計算を行った。
浸透流計算は3次元計算モデルであるDTRANSU-3D14)を使用した。計算格子は河床変動計算を行った平面格子データを、DTRANSU-3Dに対応した有限要素格子に変換し、 図-3 河床変動計算に用いた流量（雁来水位流量観測所13））
河床面から6m深部まで0.2m間隔で作成した。河床上面
表-2 浸透流計算の条件（DTRANSU３D・EL14））の境界条件として、河床変動計算で得られた河床地形に、
3項目
計算条件
平水流量14.4m /sの条件で計算して得た、平面的な河川
縦断方向：KP11.8～13.0
水位分布を全水頭として設定した。浸透流計算に用いた
解析対象範囲横断方向：低水路左右岸法尻間
透水係数は、後述する現地で実施した水頭落下試験の平鉛直方向：河床面から6m深部まで均値である57.9m/dayを用いた。また、有効間隙率は沖積
約4×4×0.2m
格子（X×Y×Z）-6 -1
一定（57.9m/day）
透水係数
礫層の0.25を用い、比貯留係数は密な砂礫の9.4×10 m0.25
有効間隙率を用いた15)。左右岸、上下流および底部の境界条件は、-6 -1
9.4×10 m
比貯留係数
これらの位置の浸透流調査を行って決定することが望ま
差分法
中央差分
定常・非定常定常計算
しいが、本研究は河床微地形による浸透流の発生に着目
節点流速0.14m/day
上流境界条件しているため、簡単のため次のように仮定した。現地の
節点流速-0.14m/day
下流境界条件
側岸の護岸がモルタル充填されたものであることを踏ま上面境界条件全水頭=河川水位
左右岸・底面境界条件
流出入なし
え、左右岸方向の浸透流の流出入はないものとした。上Masaaki Yano、 Kazuyoshi Watanabe、 Yasuyuki Hirai

図-4 現地地調査地点 ※GG：河床材料調査
査位置、HYP：浸浸透流調査位置
置

流の実流速ととは河床内の連連続した空隙内を流れる流流速でよび浸透流調査
査を行った。
およあり、浸透流流の断面平均流流速を有効空隙率で除したた値でa) 調査方法
調15)ある 。本稿
稿では以降、実実流速を浸透流流速と呼ぶ。
各調査項目別の各の調査地点を図図-4に示す。
KP11.8～13.0のc) PHABSIM
横断測量は、K横の範囲の延長1.2kmにおいてて、PHABSIMにによる評価ををするため、物物理環境の要要素毎縦断断間隔10mで行行った。に、産卵環境境としての適性性値（Suitabiliity Index Valuee：以河床材料調査は
河は、縦断間隔220mで横断測線線を設定し、5)
1横断降、SI）を算算出する適性基基準 を作成しした。適性基
基準を断測線あたりの調査数は、平水時水面下下の主流部をを作成するにはは実際の産卵床床位置における物理環境値値を把概ねね等分した4～5地点と、同測測線上の砂州の水面上の河河握する必要ががある。産卵床床位置の物理理環境値は、著著者ら
岸部部1地点としたた。河床材料はは、産卵床の造造成に影響す
すが過去に行っった方法6)と同同様に、後述すする現地調査
査結果
る河河床表面から河河床下30cm程程度までを採取取した。粒度
度mm未満の成分を内挿補間しして物理環境値値の平面分布を作成し、産産卵床
分析
析は、直径75m分はふるい分け
け分析を行いい、75mmmm以上の成分位置の値を抽抽出することでで把握した。適性基準は環
環境的分は現地で礫径径を計測した。
。ふるい分けけな偏りに対すする補正を行ったもの5)（以以降、SI3）をを式(1)、 分析
析と礫径計測のの結果は、同一一の粒径加積曲線に取りまま(2)にて作成しした。さらにに、SI3を後述すする方法で調調整しとめめ、平均粒径をを算出した。
た適性基準ををSI4とした。そして、数値値計算により予測
浸透流調査の地浸地点設定は、河河床材料調査の水面下の調調16)
axter et al が考した物理環境境値を適性基基準に代入してて、2×2mメッッシュ
査地地点と同一地点点とした。調査査方法は、Bax
考の各セルのSSIを算出した。
案しした方法で行い
い、ピエゾメーータの計測値より動水勾配配を、水頭落下試験
験より透水係数数をそれぞれ求求め、これら
ら(1))の積積から浸透流速速を算出した。
。なお、前述したとおり、
3/(22)
浸透透流計算では、この水頭落下下試験で計測された透水係係数のの平均値を用いいている。
ここで、SI3
3i：第3種適性性基準、SNi：産産卵床の確認
認数、水面下における水る調査地点の中中には調査内容
容により流速速ENi：河道内物
物理環境分布
布のセル数、cmax
大値でや水水深が大きいた
ため、調査がでできないものがあった。そそm ：c の最大ある。なお、添え字のiはは、ある物理環環境の数値帯iiに属のたため、図-4に示示すように数量量が異なる箇所所がある。
PH
HABSIMによる産卵環境のするものでああることを表すす。これらの適性基準によより得の評価結果と実実際の産卵床床られるSIは、値が0～1の指指標で、値が大大きいほど産産卵環
分布
布の比較には、札幌市豊平川
川さけ科学館
館17)が本調査と境に適していいることを表すす。
同年
年度に同区間でで行った産卵床床調査結果を用用いた。なおお、Sにより座標がSI4より算出された各物
物理環境要素ののSIから、式式(3)に産卵卵床位置はGPSが計測されていいる。
より合成適性性値（Composiite Suitability In
ndex Value。以以降、CSI）5)を、各各セルで算出しした。

3. 結果と考察結＝(3)(1) 数値計算結果
数と実測の比較較a) 河床高の比較
河（SI4）から得られた、ああるセルのあるる物理環境（（a、b、
図-5aに平成23年
図年に調査したた実測の河床高高を示す。
c…）の適性値値である。
実測
測ではKP12.2付付近の左岸部ととKP12.4付近のの右岸部に、比高高の高い砂州がが確認され、そその対岸側が深深掘れしていい(3) 現地調調査る。図-5bに示すとおり河床変変動計算においいてもこの状状
数値計算にによる物理環境境値の予測精精度を確認するるため、 況はは再現されていいる。この高い
い比高の砂州は、平成23年
年9月の平成23年に対対象区間においいて、横断測量、河床材料料調査の現地調査時時においては樹樹林化しているることが確認
認

ここで、CSISI：あるセルの合成適性値
値、SI：適性性基準

Masaaki Yanno、 Kazuyoshi Watanab
be、 Yasuyukki Hirai

図-5 実測調調査結果と数値値計算結果の比較較している。ままたこの個所はは、河床変動計算においても、平成220年に撮影されれた国土地理院の航空写真
真で、図-5bの網掛掛けの範囲が樹樹林化しているると判断し、樹木を設定した個個所である。そそのため、これらの比高のの高い砂州と対岸部部の深掘れは、、樹木の影響
響で緩流域とななった個所に土砂がが堆積し、そのの影響で対岸部の主流部がが狭窄KP12.6～12.8のして深掘れしした状況と考考えられる。Kの右岸は実測においてては、標高が高高いが、平成23年の調査時
時部は点ででは樹林化してていることを確確認している。
。しかし平成成20年
年の国土地理院
院の航空写真ででは樹木とは判判断できなかかったたため、計算ににおいて植生のの設定を行わなかった。ここの箇箇所のように経経年的に植生がが変化する場合合は、河床変変動ににも影響するた
ため、何等かのの方法で植生変変化も計算にに反映
映すると精度がが向上すると考考えられる。

Masaaki Yanno、 Kazuyoshi Watanab
be、 Yasuyukki Hirai

110~ 120

120~ 1300.08 ~ 0.09

100~ 1100.07 ~ 0.08

1.4 ~ 1.5

1.3 ~ 1.4

2.0 ~ 2.2

1.8 ~ 2.0

0.06 ~ 0.07

1.6 ~ 1.8

1.1 ~ 1.2

1.2 ~ 1.3

80~ 9090~ 100

0.05 ~ 0.06

1.4 ~ 1.6

1.0 ~ 1.160~ 70

70~ 800.04 ~ 0.05

1.2 ~ 1.4

0.9 ~ 1.0

0.8 ~ 0.9

50~ 601.0 ~ 1.2

0.6 ~ 0.7

0.7 ~ 0.8

0.03 ~ 0.04

40~ 500.02 ~ 0.03

0.5 ~ 0.6

0.8 ~ 1.030~ 400.01 ~ 0.02

0.6 ~ 0.8

0.4 ~ 0.5

20~ 300.00 ~ 0.01

0.4 ~ 0.6

0.3 ~ 0.4

10~ 20‐0.01 ~ 0.00

0.2 ~ 0.4

0.2 ~ 0.3

‐0.02 ~ ‐0.01

0.0 ~ 0.1

0.1 ~ 0.2

0.0 ~ 0.2

産卵床・環境割合，SI値

0~ 10環境割合
前期群 産卵床割合
KP12.0の周辺とその下流およびKP12.8の周辺とその上前期群 第3種適性基準（SI3)
前期群 第4種適性基準（SI4)
流は、実測に対する計算の精度は低いが、これは計算区1間の上下流端であるためと考えられる。
a 0.80.6b) 河床材料粒径の比較
0.4
図-5cに平成23年に実測した河床材料調査結果と河床0.20変動計算結果による平均粒径を示す。KP12.2下流の深掘れ部とKP12.2付近の砂州頂部の細粒状況や、KP12.3～平均粒径（mm)12.4付近の砂州頂部とKP12.4～12.8付近の粗粒状況につい
1b 0.8て、現地の実測と計算結果の傾向が一致している。ただ0.6し、平成23年の現地調査では平均粒径100mmを超える粗0.40.2粒部が確認されるのに対し、計算では確認されない。平
0成23年の現地調査では、産卵環境の評価を行うため、粗粒化している可能性がある河床表面の粒度を調査したのに対し、河床変動計算に用いた粒径は、平成18年度に通浸透流速（cm/s)
常の河床材料調査に基づき、表層を剥いで調査された結1c0.8果である。そのため、数値計算では、初期の粗粒化した0.60.4河床面の状態が反映されていない。このことが、100mm0.20を超える大きな粒径の再現性に影響した可能性がある。
c) 浸透流の比較
図-5dに平成23年に実測した浸透流調査結果と浸透流流速（m/s)
計算による予測結果を示す。浸透流速の計算結果は現地
d 10.8
調査結果と比較すると全体的に小さい。また、現地調査
0.60.4においてはKP12.3～12.4付近の砂州頂部やKP12.8のやや0.20上流の砂州において、現地調査での粒径が大きい個所では浸透流速が大きいが、計算ではその傾向はみられない。
透水係数は、例えばTerzaghiの式によると有効粒径の2乗
水深（m)に比例して大きくなる15)。現地の粗粒化した河床材料は
図-6 シロザケ産卵床の適性基準
透水係数にも影響していると考えられるが、計算では、河床材料（透水係数）の分布を考慮せず、一定の透水係指摘があり4)、これが大きいほど産卵環境に寄与するこ数を用いている。そのため、現地調査で確認された粗粒とが考えられる。そこで、SI4ではSI3の凸部を結び凹部個所の大きな浸透流速は、計算で表れなかったと考えらを埋め、SI3の最大値より大きい数値帯のSIは1.0とした。
れる。しかし、KP12.2～12.3付近の砂州頂部における河流速のSI3は、図-6cに示すとおり、0.0～0.2m/s、1.4～川水の伏流と前縁部における湧出は確認され、河床微地1.6m/sで凸となった。SI4では、凸部を結び凹部を埋めた形に依存する浸透流速は把握できると考えられる。
SI分布とした。ただし、既往研究では流速0.8m/sを超えると雌は産卵行動をほとんど行わなくなるとの指摘があ(2) PHABSIMによる産卵環境の評価
る18)。そのため、1.4～1.6m/sの凸を無視することし、0.0a)適性基準SI3の補正～0.2m/sの数値帯のSIが1となるように全体を補正した
SI3を図化したものを図-6a～dに示す。SI3は横軸に対(図-6c紫線)。そして、1.4～1.6m/sの凸を無視し、上下のして滑らかではなく、凹凸があり、そのまま使用するこ数値帯の適性基準を結んだ。とができないと判断される。このため、次に述べる方法
水深のSI4は、SI3の凸部を結び、凹部を埋めたSI分布
でSI4を作成した。とした。ただし、0.0～0.1mの数値帯は、水深が浅く、シ平均粒径のSI4は、図-6aに示すようにSI3の凸部を結び， ロザケの遊泳が困難であると考えSIを0.0とした。
凹部を埋めたSI分布とした。ただし、サケ科魚類が使用b) 数値計算結果によるPHABSIMと産卵床の対応できる50%粒径は体長の1/10との知見があるように3)、大PHABSIMにより算出したCSIと産卵床の分布を図-7にきすぎる粒径は産卵床造成の支障となることから、100示す。産卵床が比較的多いKP12.2～12.4付近など、砂州～110mmの適性値の凸は無視し、その上下の数値帯の適の前縁部付近にCSIが高い個所が確認される。KP12.0付性基準値を結んだ。
近やKP12.8付近においても産卵床が多く、CSIが高い個
浸透流速のSI3は、図-6bに示すように0.01～0.02cm/sお所がみられるが、前述したとおり、これらの個所では計よび0.06～0.07cm/sで凸状に値が大きい。シロザケは浸透算区間の上下流端付近であり、河床変動計算の再現精度
流が湧出傾向の個所を産卵環境として利用しているとのが低いため、産卵床とCSI分布がずれている。
Masaaki Yano、 Kazuyoshi Watanabe、 Yasuyuki Hirai

図-7 PHABSIM によよる合成適性値
値と産卵床との対対応

4．まとめ
本研究ではは、数値計算ににより取得した物理環境値値からPHABSIMにより産卵環境境を評価し、現現地調査結果果と比ことが明らかとなった。
較し、次のこ①計算区間のの上下流端付近
近や、計算対対象期間に現地地で植生が変化した
た個所以外は、、河床変動計
計算により河床床高と粒径分布の傾傾向を概ね再現現できた。
②浸透流計算算の結果、現地地の浸透流速速分布と比べ小小さな値となっていいる個所があるる。これは、計算では一定定の透水係数を用いいたため、現地地で確認された大きい河床床材料に起因する高高い透水係数のの影響を加味味していないたためと考えられる。しかし、砂州
州地形における頂部での河河川水び前縁部での湧湧出といった、河床微地形形に依の浸透および存する浸透流流速は再現でききている。
③PHABSIMににより産卵環環境を評価した
た結果、地形形形状の再現性が低低い計算区間の上下流端以以外は、CSIと産卵2～12.3付近の床の分布が概概ね一致したた。特にKP12.2の産卵床が多く確認認された砂州前前縁部付近の合成適性値がが高くなり、産卵床床の分布と一致致した。
ABSIMによる
今後、PHAる産卵環境の予予測・評価精精度を上げるためにには、河床変動
動計算において、植生変化化や深度方向の河床床材料分布（特特に河床表面の粗粒層）をを反映することや、浸透流計算ににおいて、河床河床内部の透水水係数慮することが考考えられる。の分布を考慮
定期横断測量デデータおよびび河床材料調査
査デー
謝辞：河川定タを提供してて戴いた札幌札札幌開発建設部に、ここ記記して謝意を表す。
参考文献
1) 水産総合研研究センター北北海道区水産研究究所：さけますす来遊速 報 （ 平 成25年 度 ） .
URL
http://salmon.frfra.affrc.go.jp/zousyyoku/H25salmon/hh25salmon.htm
2) 永沢亨：ササケ科魚類のプ
プロファイル-10 サケ. SALMOON 情報 6: 22-25,20123) Kondolf
KM.G. & Wolman G.M. : The sizes of Saalmonid Spawningg
Grravels. Water Resouurce Research vol.2.29 (7), pp2275-22285, 1993,1.
4) Geeist D.R., Hanrahan
n T.P., Arntzen E.V
V., McMichael G. A., Murray C.J. &Chhien Y. : Physicoch
hemical characteriststics of the hyporheeic zone affect redddsitee selection by chuum salmon and faall chinook salmon
n in the Columbiaa
Rivver. North Americcan Journal of Fishheries Managemen
nt, Vol.22, pp1077-10085, 2002.
5) 玉井信行ら：河玉
河川生態環境評
評価法，東京大大学出版，pp95-1002,2000,3.
6) 矢野雅昭ら：豊平矢
平川中流部におおける PHABSIM
M によるシロザ
ザケ産卵環境の評
ケ評価，河川技術術論文集，Vol.19，pp507-512，20013,6.
7) 鈴木俊哉：自然鈴
然再生産を利用し保全への取り組
組したザケ資源保
み.
み SALMON 情報 2: pp3-5,2008.
8) Malcolm
MI.A., Soulssby C., Youngson A.F. & Hannah D.M.
D ：Spatial anddtemmporal variability ofo groundwater-sururface water interacctions in an uplandd
sallmon-spawning strream: implications for egg survival. Hydrology:
HSciencee& Practice for the 21tth Century, Vol.2, ppp130-138, 2004.
9) QGIS，http://wwwQw.qgis.org/
10) iRICiProject，httpp://i-ric.org/ja/11) 国 土 地 理 院 ， 地 図 ・ 空 中 写 真 閲 覧 サ ー ビ ス ，httptp://mapps.gsi.go.jpp/12) 内田有子ら：豊豊平川の洪水時時における河畔林林における影響
響に関する研究，土に土木学会北海道道支部，第 67 号，B-29,2010.
号13) 国土交通省，水国
水文水質データベ
ベース，http://wwww1.river.go.jp/14) 西垣
西誠，三菱菱マテリアル株
株式会社，株式会
会社ダイヤコンン
サルタント：オイサ
イラリアン・ララグランジアン・3 次元飽和・不飽和浸透流－移不移流分散解析プ
プログラムデーータ入力マニュュアル，2001,10.
ア15)土
土木学会 水理公公式集. 社団法人
人 土木学会,1999.
16)B
Baxter,C. Hauer, R.F.
R & Woessner W.W.： Measurring Groundwater-Strream Water Exchaange: new Techniiques for Installingg Minipiezometerss
andd Estimating Hyd
draulic Conductibbity, Transactions of the Americann
Fissheries Society 132
2: 493-502, 2003.
17)札札幌市豊平川さけ科学館，http:/://www.sapporo-paark.or.jp/sake/18)TTiffan F.K., HaskelllC.A & Kock T.J.：Quantifying the behavioral
ressponse of spawning
g chum salmon to eelevated dischargees from BonnevilleDaam, Columbia riverr, USA. River Resesearch and Applicat
ations 26: 8710
01.2010.

Masaaki Yanno、 Kazuyoshi Watanab
be、 Yasuyukki Hirai