気候変動による釧路川の治水・環境に関する対策提案を目指した検討
Investigation of an impact of climate change on Kushiro River

中山恵介1

Keisuke Nakayama
1

北見工業大学要旨
気候変動による影響が世界的に現れており，その影響は極に近くなるほど大きくなることが報告されている．そのため，日本の中でも特に北海道地方では，降雨形態の変化，海水面の上昇などにより，これまでの対策では対応できない治水・環境に関する問題が発生することが予想される．北海道道東に位置する釧路川は，沿岸域の開発による人口集中，その背後の豊かな湿原の存在という特徴を持ち，気候変動による治水・環境への影響が大きな流域であると言える．そこで本研究では，気候変動の影響により釧路川が受ける治水・環境への影響を評価することを目的とする．
目的達成のため，本研究では2種類の検討手法を適用した．まず，降雨量変化を評価するための手法として有効である，分布型流出モデルによる流出モデルの開発である．モデルには，表面流，浸透流，河道流が考慮されており，長期的な貯留の効果などを考慮することが出来る点などが大きな特徴である．また，気候変動の影響により海水面が上昇することにより影響を受ける現象として塩水遡上が考えられる．その塩水遡上に対する影響評価を行うための手法として，観測結果を用いて検証を行い，3次元数値モデルの適用を行った．それらの結果，分布型流出モデル及び塩水遡上モデルが良好に現象を再現できることが確認された．
《キーワード：分布型流出モデル，釧路湿原，3次元数値モデル，Fantom3D，塩水遡上》

１．釧路川近年，世界中で気候変動が起こっているのが確認されている．気候変動は強い台風や海面上昇を引き起こしており，洪水や高潮等が原因となる災害が増加していると考えられる．さらに，2007年の研究に基づきIPCC（気候変動に関する政府間パネル）は21世紀の終わりまでに，海水面水位が59cm程上昇する事を予測している．ラムサール条約によって登録されている北海道の釧路湿原には流域に住民が約20万人弱住んでおり，高度に開発されている地域であるといえる．また，釧路湿原は釧路川によって涵養されており，その釧路川の河川長は154km，河川流域面積は2510km2であり勾配は比較的緩やかである（図-1）．そのため洪水や高潮等が原因となる災害が発生する
図-1 釧路川流域図．と，被害が大きくなると予想される．さらに釧路湿原は保存されるべき重要な生態システムを持っており，保護されなければならない．見解の両方の点から気候変動が引き起こす洪水や高潮等による災害により，釧路湿原の生態システムと開発エリアに重大なダメージを引き起こす事が考えられる．

２．対象流域の地形データ
本研究では，気候変動が沿岸域に及ぼす影響を評価することが出来るモデルを開発することを研究の上位目標としている．そこで，流量観測所標茶より下流の流域を解析対象流域とした（図-2）．河川長は40km，河川流域面積は1620km2，湿原面積は183 km2である．対象流域の土壌構成は流域の大部分が，崩積性未熟土壌23%，残積性未熟土壌66%である（図-3，表-1）．また，流域の土地利用区分は畑地25%，森林55%であり流域の大部分を占めている（図-4，表-2）．

図-2 釧路川流域および解析対象流域図．図-3 対象流域土壌分布．

図-4 対象流域土地利用区分．

表-1 釧路川対象流域での土壌分布面積

土壌分布
岩屑性土壌
崩積性未熟土壌
残積性未熟土壌
粗粒残積性未熟土壌
砂丘未熟土壌
粗粒火山拠出物未熟土壌
くろぼく土壌(a)くろぼく土壌(b)面積

割合

9.5km2367.5 km2
1070.75 km2
150 km2
0.5 km2
10 km2
10 km2
1.25 km2

0.6%
22.7%
66.1%
9.3%
0.0%
0.6%
0.6%
0.1%

表-2 釧路川対象流域での土地利用区分面積

土地利用区分
畑地
森林
荒地
市街地
交通用地
その他
河川
ゴルフ場

面積
398.75 km2
896.5 km2
189.5 km2
24.5 km2
3.5 km2
8.75 km2
90 km2
7.5 km2

割合
24.6%
55.4%
11.7%
1.5%
0.2%
0.5%5.6%0.5%

３．河道データ
本本解析で使用する河道網は国土交通省，国土数値情報を用いて作成した．支川数は13，総メッシュ数は520である（図-5）．各支川に与える河道断面情報は国土交通省釧路開発局計測の釧路川横断測量データを使用する（図-6）．
River 0標高 (m)

River 61086420-2-4-6-8-10River 11River 12-100

River 70River 51086420-2-4-6-8-10

River 9River 8

標高 (m)

River 3
River 1River 2River 0-200

図-5 釧路川流量再現計算におけるモデル河道網．

River 30200
400600距離 (m)

River 815131197531-1-3-5標高 (m)標高 (m)

1086420-2-4-6-8-10-5000距離 (m)

500

-400

River 10-2000 (m)200
距離400

River 122523211917151311975

標高 (m)

標高 (m)

14121086420-200300

River 2River 10River 4

100200距離 (m)-1000 (m)100
距離200-200

図-6 釧路川流域および解析対象流域図．

0 距離 (m)
200

400

800

４．流量再現計算とまとめ
再現対象期間は2004年9月5日から9月11日における降雨による流量に対して行った．表-3に流量再現計算の計算条件を示す．流量の再現計算は図-7に示す流量観測所広里，五十石，幌呂，下久著呂において行った．解析に使用した降雨は，気象庁アメダス観測所塘呂観測所における1時間データを使用し，流域に一様に与えている．本研
究では標茶観測所より下流の流域を解析対象としているため，河道の上流端に国土交通省北海道釧路開発局計測の流量観測所標茶における流量データを１時間ごとに与えている．また今回の解析では，浸透層厚や透水係数における詳細なデータが存在しないため，ハイドログラフから推定し，流域に一様にTry and Errorで決定した．

表-3 釧路川流量再現計算条件．

⊿x=⊿y=500（m）⊿t=10(sec)ks=0.2×10-4 (m/sec)
H=5.0(m)図-7 釧路川流量観測所地点図．

図-8に流量観測所広里，五十石，幌呂，下久著呂における流量再現結果示す．図における実線は釧路開発局計側の観測値，破線は計算結果である．

流量観測所 広里

流量観測所 五十石

流量観測所 幌呂

流量観測所 下久著呂
図-8 2004年9月5日から9月11日の釧路川各流量観測所における流量再現結果．
流量再現結果から，各流量観測所においておおまかな再現を示すことが出来ており，特に下流域に位置する流量観測所広里において良好な再現性が得られたことは本モデルの流量の再現性の高さを示している．しかし，その他の流量観測所におけるピーク流量および流量の減衰部分等，今後検討の余地が残る．減衰部分の再現性が低い要因として，地表面の粗度および浸透モデルのパラメータが推測される．今回の解析では土壌分布における浸透のパラメータや土地利用区分における地表面の粗度を反映しておらず，流域に一様に取り扱った．今後，より多くの観測データを集め，河道形状や粗度を再度詳細に検討することで改善が期待できると考えられる．

５．釧路川における塩水遡上に関して
近年の気候変動によって，規模の大きい台風やそれにともなう自然災害が引き起こされており，加えて環境への直接的・間接的な影響も問題となってきている．気候変動は，海面上昇もともなって発生することが報告されているため，沿岸地域が陸地の32％を占め，住民の46％が沿岸地域に住んでいる日本では，その影響が顕著化する恐れがある．その

3.00

塩分計設置：平成21年12月26日（小潮時：平成21年12月25日調査）

2.50（大潮時：平成21年12月7日調査）
2.00

水位 (TP m)

1.50
1.00
0.50
0.00

-0.50-1.00

岩保木(KP11.0)

-1.50

12/1

12/6

12/11

広里(KP7.4)
12/1612/21

釧路港潮位
12/2612/31

図-10 縦断観測を行った大潮時(平成21年12月 7日)および小潮時(平成21年12月25日)における水位と潮位．

広里流量(m3/s)

中で，環境面が大きな影響を受ける領域として，ラム120
豊水
サール条約によって登録され，日本最大の湿地面積を誇
110
流量る釧路湿原を挙げることが出来る3)．釧路湿原は流域人
100
平水
90
口が約230,000人であり，釧路湿原の中心を流れる釧路流量80
川は，河川長154 km，河川流域面積2,510 km2の一級河川70である．
3360
釧路川の平水時流量は30 m /s から50 m /sであり，河50
口からの塩水遡上によって湿原域へ塩水が進入する可能40性がある．釧路湿原における植生はヨシを除き基本的に30
塩分に弱く，海面上昇の影響で塩水がより遡上すること
1999
2001
2003
2005
2007により，その豊かな湿原域が失われる可能性があると懸（年）念されている．しかし，これまでに釧路川を遡上してく
図-9 新釧路川における過去10年間（2009年度から2018
る塩水の形態に関する現地観測や解析結果は存在せず，年度）の豊水流量，平水流量，低水流量，渇水流量．
僅かに塩水の遡上に関して河口域付近での調査が存在す
るのみである．それ故，遡上形態については不明な部分が多く，至急，解明されなくてはならない．
塩水の遡上に関しては，過去に多くの研究例がある．例えば，同じ道東に位置する網走川では遡上により上流側の網走湖の水質環境が大きく影響を受けていることから，遡上形態に関して多くの調査がなされている．それらの論文では，潮汐による影響，風による影響などが指摘されており，釧路川においても同様な影響により塩水遡上が影響を受けている可能性がある．そこで本章以降では，釧路川を遡上してくる塩水に関する現地計測を行い，その結果に基づいて3次元数値計算モデルによる検討を行い，遡上形態について理解することを目的とする．

６．釧路川下流の領域と観測結果
釧路川は，カルデラ湖の屈斜路湖から流れ出し，釧路湿原内を流れる一級河川であり，釧路川流域は，山地等が約68%，牧草地等の農地が約21%，釧路湿原が約8%，宅地等の市街地が約3%という，酪農と湿原という特徴を持
つ．河口には重要港湾の釧路港をもち，上流の屈斜路湖などは阿寒国立公園であり，下流の釧路湿原はラムサー
ル条約登録湿地及び釧路湿原国立公園に指定されているなど経済的にも重要な領域を持ちながら，豊かな自然環境に恵まれている．

2.50
2.00

塩分計点検：
平成22年2月4日

塩分計設置：
平成21年12月26日

塩分計撤去：
平成22年3月3日

1.50

水位 (TP m)

1.00
0.50
0.00

-0.50
-1.00
-1.50

岩保木(KP11.0)

12/20 12/25 12/30 1/4

1/9

広里(KP7.4)

1/14 1/19 1/24 1/29

釧路港潮位

2/3

2/8

2/13 2/18 2/23 2/28

3/5

図-12 2009年12月25日から2010年3月3日までの間で行われた塩分・水温のKP0.0からKP8.0の1km間隔における時系列
観測時の水位と潮位．

水深(m)(1) 釧路川下流の領域の特徴
2009年度における観測結果を示す前に，対象とする塩分濃度(psu)
釧路川下流から河口における外力条件を理解するため，0102030
水文特性，潮位，流量などのデータを調べた．水文
0.0
データについて，過去10年間の雨量データから，年間（小潮時） KP-0.2
総雨量は約800 mmから1250 mmで変動していることが分かった．観測を行った2009年度は他の年に比べて雨0.5量が多く，平年より多雨であったことも確認された．
また，釧路港における最大潮位差は0.51 mから0.93 m1.0であり，高潮時の最高潮位は0.64 mから0.93 mの範囲であることが確認された．顕著な高潮は冬期の低気圧により発生する場合が多く，この時期は河川流量が渇1.5水となる時期と重なることから，冬期間には塩水が遡上し易い環境であることが分かった．
潮位と同じく，塩水遡上の重要な要因である広里に2.0おける河川流量について過去10年間にわたって調査を行った(図-9)．近10ヶ年の流況は，2008年に最小流量である31.02 m3/sが発生していた．また、年総量につ2.5いても2008年に最小値を示しており，渇水年であった
図-11 KP-0.2における平成21年12月25日の小潮時での塩ことが分かる．なお，気象観測結果でも2008年に総降分の鉛直分布．
雨量が小さな値となっており，気象観測結果と同様な結果が得られている．観測対象とした2009年度に関しては，まだ詳細なデータが得られていないが，総降雨量が大きく，流量も平均よりは大きかった点を記しておく．
(2) 観測結果
現地観測で2種類の観測が実施された．一つは，多項目水質計(ALECAAQ1183)を用いた水深・水温・塩分・濁度・Chl-a・DOの観測であり，河口KP-0.2からKP8.0までの間の縦断方向500 m間隔および水深方向10 cm間隔で実施された．もう一方は，KP0.0からKP8.0までの1 km間隔で最深部に設置された水温・塩分計(ALEC CT計)による10

(a)

2.0
1.5

kp-0.2

1.0

kp0.0

0.5

kp0.4

0.0

標高 (m)

-0.5
-1.0
-1.5
-2.0

(b)

-2.5
-3.0
0

50

100

150

200

250

300350距離 (m)

図-13 新釧路川におけるKP-0.2からKP0.4までの横断面
図．

図-14 (a) 数値計算に用いられた現況断面．(b) 河口を急縮した現況断面．

分間隔の塩分遡上観測である．
塩水遡上の状態を観測するための縦断観測は，大潮時(平成21年12月 7日，8：40～14：10)および小潮時(平成21年12月25日，8：50～15：30)を狙って実施された(図-10)．塩分縦断観測は，作業船（FRP製，長さ2.96m，搭乗人員3名，8PS搭載）を用い，最下流地点（KP-0.2）より上流に向かって実施された．KP-0.2地点は、ほぼ海上となるため朝方の凪を利用して観測を開始し，各横断測線上の最深河床部付近1地点で多項目水質計を用いて鉛直方向10cm間隔でデータが取得した．大潮時の観測では，図-3に示されているとおり，雨の影響により増水しており，流量が大きくなっていたため，塩水の遡上は全く計測されなかった．一方，流量が落ち着いてきている小潮時における観測では，KP-0.2地点においてのみ塩水の進入が確認された(図-11)．底層から約40 cmほどであるが，塩水が進入していた．しかし，他の地点においては，大潮時の観測結果と同様に，塩水遡上の影響を確認することはできなかった．
大潮時・小潮時において塩水遡上を確認できなかったため，時間的により詳細に計測を行う必要があると考え，
2009年12月25日から2010年3月3日までの間，水温・塩分計を9つ，KP0からKP8までの1 km間隔で最深部に設置し，10分間隔での計測を行った(図-12)．その結果，KP0地点において2回の塩分濃度上昇を記録していたが，その他の地点では全く塩分上昇が見られなかった．この2回の塩分上昇は，チリ津波による津波の河川への遡上と一致しており，通常はあり得ない変化であったことを記しておく．
観測を行った年は，例年に比較して総雨量が大きく，流量が普段より大きかったため，塩水が遡上しにくかったことが考えられる．また，釧路川河口には砂州が発達しており，河口において急激に川幅が狭くなり，水深が増大している(図-13)．そのため，上流からの河川流量による，急縮部における流速の増大が塩水の遡上を抑制した原因であるとも考えられる．
７．現況における塩水遡上
2009年度における現地観測の結果から，釧路川における明確な塩水遡上が発生しているとは言えないことが示された．しかし，2009年度が多雨の年度であったことから，流量の増加による塩水遡上の抑制が考えられ，他の年度では塩水が遡上している可能性があるかもしれない．そこで本研究では，同じ冬期ではあるが流量が平年の流量に近い状態である2008年1月23日における流量を用いて塩水遡上に関する検討を行うこととした．塩水遡上の計算には，3次元環境流体モデルFantom3D9)を利用することとした．Fantom3Dはオブジェクト指向型環境流体モデ
ルであり，様々な計算スキームの切り替え，計算領域の分割による並列計算が容易であるという特徴を持つモデルである．詳しくは参考文献9)(新谷・中山，2007)をご参照いただきたい．
現状地形における検討を行うために利用する地形データは200 m間隔で横断面測量結果が存在しているのみであり，今回考慮したい河口の地形急縮部のデータが存在しなかった．そこで本研究では，存在する測量データのみ

図-16 現状地形における塩分の最大進入の様子(大潮時)．(a) 現状地形no1である河口急縮なし354時間後．
(b) 現状地形no2である河口急縮あり366時間後．

図-15 現状地形no1，図-14(a)のA-A’断面における塩分の時間変化．(a) 初期状態．(b) 48時間後．(c) 96時間後．

図-17 数値計算に用いられた計画断面．

で地形を再現したものを用意し，現状地形no1と定義した(図-14a)．そして砂州の発達による河口幅の急縮と深掘れの効果を考慮するために現状地形の河口部分を急縮断面へと変化させた地形を現状地形no2と定義した(図-14b)．
この2種類の地形を用意することにより，現状地形において急縮地形がいかにして塩水遡上に影響を及ぼしているかを検討することが出来る．
計算領域は，8000 m×300 m×3 m(縦断方向×横断方向×水深方向)とし，メッシュ数は160×30×22とした．
それほど多くのメッシュではないが，さらに計算速度を向上させるために並列計算を行い，領域を8つに分割し，20×30×22での計算を可能とした．上流端の条件は，2008年1月23日における流量である36.25 m3/sを与え，潮汐条件としては大潮を模擬した振幅0.7 m，小潮を模擬した0.2 m，潮汐振幅なしの3種類を用意した．初期条件では河道を塩水(34 psu)で満たし，上流端と下流端の水深差0.2 mを再現するために，上流の水深を6時間かけて上昇させた．計算時間ステップは1 sとした．計算は実時間で400時間行い，定常状態に達することを確認した．
計算における状態の変化を示し，実行された計算の理解を深めるために，計算開始直後から計算開始96時間後までの横断面中央を取り出した塩分の縦断分布を示す(図-15)．現状地形no1における潮汐条件なしの計算結果を示している．上流端からの淡水の流入により，塩水が下流端から流れ出してゆき，定常状態に達するとほぼ塩水が河道からなくなってしまっていることが分かる．また，およそ100時間で定常状態に達していた．続いて，大潮および小潮における塩水遡上形態について検討を行った(図-16a)．潮汐の大きな大潮時においても塩水遡上距離はそれほど大きくなく，最大で河口から1 km程度の遡上距離となっていた．
砂州の発達による急縮効果を検討するために，現状地形no2における同様な計算を行った．本論文では遡上距離が最も大きかった大潮時における検討結果のみ示す(図-16b)．現地観測の章において考察された通り，急縮部が存在することにより河川流量によるフラッシュの効果が増大し，僅かではあるが塩水が遡上しにくい状態になっ図-18 図-17の計画断面における塩分の最大進入の様子(大潮時)．366時間後．ていることが確認された．これらの計算から，冬期における塩水遡上はそれほど大きくなく，その遡上距離の抑制は主として堆砂によるものであり，砂州の発達による河口幅の急縮効果もその効果を大きくしているのではないかと推測された．
８．計画断面における塩水遡上
前章では，現状における塩水遡上の状態について検討が行われ，塩水遡上距離が現状地形により抑制されていることが示された．そこで本章では，堆砂影響を評価するため，堆砂前の状態における塩水遡上の評価，つまり計画断面における塩水遡上の形態について検討することとした．計画断面における計算地形を図-17に示す．計算は，現状地形と同様に3種類の潮汐条件を与えて行われた．ほぼどのケースも同程度の塩水遡上距離であったが，もっとも張り出していた大潮時における計算結果のみ示す(図-18)．
５．塩水遡上に関するまとめ
釧路川における塩水遡上形態に関する検討を行い，以下のような結論を得た．
(1) 2009年度冬期に多項目水質計を用いた縦断観測および塩分・水温計による時系列観測を行い，明確な塩分遡上が発生していなかったことが確認された．その塩分遡上の抑制には，河口における砂州の発達等の地形の効果が影響を及ぼしていることが分かった．
(2) 通常年における冬期の塩水遡上形態に関する検討を行うため，3次元数値計算モデルを用いた検討を行い，河口閉塞させない現状地形でも塩水が遡上しにくい状態であることが示された．
(3) 河口における砂州の発達の影響を考慮した数値計算を行い，その影響により塩水の遡上が抑制されていることも示された．
(4) 堆砂の影響による塩水遡上の抑制効果を検討するために，堆砂前の計画断面における塩水遡上計算を行い，堆砂が塩水遡上抑制に大きく影響していることが確認された．
今後，冬期のみの影響だけでなく，他の季節における影響評価も検討し，釧路川における塩水遡上形態について解明してゆく必要がある．
参考文献
1) Solomon, S., D. Qin, M. Manning, Z. Chen, M. Marquis, K.B. Averyt, M. Tignor and H.L. Miller (eds.), Contribution ofWorking Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change, Cambridge University
Press, Cambridge, United Kingdom and New York, NY, USA, 2007.
2) Dutta, D., M.S. Babel and A. Das Gupta (2005). An Assessment of the Socio-Economic Impacts of Floods in Large
Coastal Areas, Final Report for APN CAPaBLE Project: 2004-CB01NSY-Dutta, Asian Institute of Technology, ISBN 97493908-0-6, 205 p., 2005.
3) Kushiro Wetland Restoration Committee, Comprehensive Concept, Japan, 2005.
4) 桒島知哉, 高橋克人，宮島滋近,平野道夫，山田正: 塩淡二成層を形成している網走湖の塩水の流出入に関する研究, 水工学論文集, 第37巻, pp.305-312, 1993.
5) 池永均, 向山公人，大島伸介，山田正: 塩淡二成層を形成する汽水湖沼の長期的な界面変動予測手法の開発,
土木学会論文集, 第628巻, pp.77-96, 1999.
6) 池永均、大束淳一、三沢大輔、内島邦秀、樫山和男、山田正：網走川における塩水遡上特性と塩水楔の形態変化に関する研究、水工学論文集第42巻,pp.775-780,1998.
7) 吉川泰弘，渡邊康玄, 結氷時の塩水遡上に関する現地観測 土木学会，河川技術論文集, Vol.12, pp.157-162,
20068) 吉川泰弘，安田浩保，渡邊康玄, 網走川における透過性構造物の塩水遡上への影響 土木学会，河川技術論文集, Vol.13, pp.7-12, 20079) 新谷哲也, 中山恵介, 環境流体解析を目的としたオブジェクト指向型流体モデルの開発と検証, 水工学論文集,第53巻, pp.1267-1272, 2009.