Introduction to iRIC

Yasuyuki Shimizu, Tomoko Kyuka, Takuya Inoue

and Jon Nelson
Sapporo, Hokkaido
 Sapporo

HOKKAIDO

JAPAN
800km

Tokyo

3
 Hokkaido

Area: 83,000km2

SAPPORO Population:
5,700,000

Sapporo
Population:
1,900,000

4
 Where Hokkaido University is

Main campus located in

Sapporo, Hokkaido Hokkaido University
Sapporo campus
Another campus in Hakodate
(Grad. School of Fisheries
Sciences)

Several off-campus facilities

mostly in Hokkaido

Established in 1876
11,000 Under Graduate Students
6,000 Graduate Students
2,000 Teaching Staff
1,800 Non-Teaching Staff
12 Faculties + 10 Institutes
6
7
8
 9
K.Asahi, Y. Yoshida, H.Tsunematsu, Y. Shimizu and
J. Nelson(2012), Development of the iRIC Software
for River analysis
 10
K.Asahi, Y. Yoshida, H.Tsunematsu, Y. Shimizu and
J. Nelson(2012), Development of the iRIC Software
for River analysis
 11
K.Asahi, Y. Yoshida, H.Tsunematsu, Y. Shimizu and
J. Nelson(2012), Development of the iRIC Software
for River analysis
Hokkaido Food?

Welcome Hokkaido !
13
 3
English Engineering Education Program (e )
Hokkaido University, Japan

Japanese Government
Scholarships are available for admission
in October every year. Deadline for application
is 15 December of the preceding year.

http://www.eng.hokudai.ac.jp/e3/
 Sapporo

HOKKAIDO

JAPAN
800km

Tokyo

15
Typhoon passes between 1951 and 2018 of August

16
国立情報学研究所・デジタル台風サイトより
Typhoon passes between 1951 and 2018 of September

国立情報学研究所・デジタル台風サイトより
17
 August 2016 in Hokkaido

Abashiri
網走
留萌
Rumoi
Kitami
Asahi
旭川
kawa 北見

Typhoon #7
台風7号
th AUG
178/17上陸
小樽
Otaru
札幌
Sapporo
帯広
Obihiro 釧路
Kushiro

台風11号#11
Typhoon
8/21上陸
21th AUG
室蘭
Muroran
Typhoon #10 Typhoon #9
台風9号
31th AUG
台風10号 23th AUG
8/23上陸

函館
Hakodate
Precipitation between 1:00 15th AUG and 24:00 24th AUG

First three Typhoons

Precipitation between 1:00 29th AUG and 9:00 31th AUG

Kushinai
29th-31th AUG
515mm
 Sorachi River, AUG 2016
Hokkaido, Japan

動画：北海道開発局
22
Akatani River in Kyushu, 2017
 International
River Interface Corporative

Akatani River in Kyushu, 2017 iRIC

 International

Kinu River Flood, September 2015 River Interface Corporative

iRIC

9月9日 16:00 9月9日 20:00

9月10日 0:00 9月10日 3:00

Strong
26 rain band stayed along the River for 2 days
 International
River Interface Corporative

iRIC
 International
River Interface Corporative

iRIC

28
 International
River Interface Corporative

iRIC
 International

Kinu River Flood of 2015 River Interface Corporative

iRIC
1981 Flood in the Ishikari River, Hokkaido, Japan 31
1981 Flood in Ishikari River 32
Inundation of the Ishikari River, 6 AUG, 1981
33
 34
Omagari point, Ishikari River 1981
Copyright iRIC Project.
35
 家屋被災

2003 Appetsu River, Hidaka, Hokkaido, Japan

36
37
2011, Levee Beach Otofuke, Hokkaido
Seki River, JUN 1995

40
Pekerebetsu River, 2016 Hokkaido

写真提供：(株）パスコ

41
 International
River Interface Corporative

Introduction to iRIC Project iRIC

International River Interface Cooperative

Copyright iRIC Project. 42
 What is iRIC?
International River Interface
Cooperative
iRIC is a public-domain modeling interface and
an associated group of open-source models that
can be used to simulate flow, sediment transport,
channel and bank evolution and habitat in
riverine environments across a wide range of
temporal and spatial scales
Using iRIC does not mean you need to use the
models we present!
Interface developed by USGS for Interface developed by Hokkaido
FASTMECH and STORM Uni. and HRDPRC for NAYS
MD_SWMS RIC-NAYS

iRIC 1.0 (Old Version, 2009)

FASTMECH, NAYS2D, MORPHO2D, STORM

iRIC 2.0 (Released on 25 May, 2012)

NaysCube(3D), NaysFlood, Delft3D, River2D…..
Confluences, Inundation, Multi Inflow, Bank erosion,
Mixture Grain Size, Ground Water, Habitat, Multi Language ・・・・・

iRIC 2.3 (Released on 23 Jun, 2014)

+ Nays2DH, Mflow_02, CERI1D, NaysEddy, ELIMO, SRM・・・・・
44
 International
River Interface Corporative
International River iRIC
Interface Cooperative

iRIC Project
iRIC Japan

・Universities ・USGS(USA), Illinoi Univ.

・Research Institutes ・Dltaras(Netherlands）
・Foundations ・Alverta Univ.(Canada)
・Consulting ・Torento Univ.(Italy)
Companies ・EDF(France, UK)
・IT Companies ・SNU(Korea)
・Tsinghua University (China)
 Outline of iRIC
Import Data Data Exporting

Pre Processing User’s Post Processing

GUI Tools Solvers GUI Tools
Vector, Contour, Shading,
Editing Geographic Data RIC iRIC Stream Line, Particle
Grid Generation and Eidt Libraries Libraries Line Graph…..

2D, 3D Fortran、C ….
Structured, Unstructured 1D, 2D, 3D, Structured,
Flow, Bed and Bank EvolutionUnstructured
Public Domain Solvers
 Solvers Line-Up for iRIC
• FASTMECH (2D Steady)
• STORM (2D Unstructured) Nays2DH
• Nays2D (2D Unsteady)
• Morpho2D (2D Unsteady)
• NaysCube(3D) iRIC Version 1.0
• Nays2DFlood(Inundation)
• Delft3D(Delft Flow), River2D Nay = “River” in Hokkaido’s
• CERI1D(Flow with Ice), Nays1D Original Language
• EvaTrip(Environmental Planning)
• ERIMO(Tsunami) • Morpho2D(Debris Flow)
• Mflow_2(2D FEM) Version 2.3
• SRM(Ruoff Model) (JUN 2014) Version 3.0 (APL 2018)
• NaysEddy(3D LES) Coming Soon!
• RRI(Rain-Runff-Inundation)
• NaysHybrid(1D and 2D Coupling), TELEMAC, Nays2dV
Copyright iRIC Project.
RRI Model (Rainfall-Runoff-Inundation Model)

iRIC-RRI is coming in December 2018

48
 International
River Interface Corporative
http://i-ric.org/en/ iRIC

49
 iRIC version2.0
www.i-ric.org International

iRIC Seminars River Interface Corporative

iRIC

Nepal, MAR 2012

Pittsburgh, USA, APR 2012

Seoul, Korea , JAN 2013

50
 International
iRIC Class in Kairo Egypt, JAN 2011 River Interface Corporative

iRIC
 International
River Interface Corporative

iRIC
Arab Spring International
River Interface Corporative

iRIC
 iRIC version2.0
www.i-ric.org International
River Interface Corporative

iRIC
Tokyo, Japan APR 2013 ２０１３．１１．１５京

Illinois, USA
Delft, the Netherlands
Bangkok Thailand
Taipei, Tiwan
Hanoi Vietnam
…….. 54
 iRIC version2.0
www.i-ric.org International
River Interface Corporative
Jakarta 16-20 SEP 2013, Organized by JICA iRIC

55
 iRIC version2.0
www.i-ric.org International
River Interface Corporative

iRIC

2014.7.3 名古屋大学

2014.7.26 インド工科大学＠ニューデリー
2014.7.30 Muroran Institute of Technology
 iRIC version2.0
www.i-ric.org International
1-3, April 2015 in Bangkok River Interface Corporative

iRIC

57
iRIC version2.0
www.i-ric.org International
River Interface Corporative

iRIC

30-31 AUG, 2016

Colombo Sri Lanka
 iRIC version2.0
www.i-ric.org International
River Interface Corporative
Jakarta 29 SEP – 3 OCT, 2015 iRIC

59
 iRIC version2.0
www.i-ric.org International
River Interface Corporative

2015.11.24-26 Nagoya iRIC

60
 iRIC version2.0
www.i-ric.org International
Katmandu, Nepal, 1 MAR, 2016 River Interface Corporative

iRIC

61
 iRIC version2.0
www.i-ric.org International
River Interface Corporative

iRIC

16-19 April 2016, Tucson Arizona USA

62
 iRIC version2.0
www.i-ric.org International
Beijing, China, NOV 19-23 2014 River Interface Corporative

iRIC

12-14 November 2014 at Kyushuu Univ.

63
 iRIC version2.0
www.i-ric.org International
River Interface Corporative
MAY 11-13 China
iRIC

64
iRIC version2.0
www.i-ric.org International
River Interface Corporative

iRIC

65
 iRIC version2.0
www.i-ric.org International
River Interface Corporative

iRIC
IAHR-APD 2014, Hanoi Vietnam

66
 iRIC version2.0
www.i-ric.org International
River Interface Corporative

iRIC

2014.10.13-17 iRIC ASEN 講習会＠ジャカルタ

67
 iRIC version2.0
www.i-ric.org International
River Interface Corporative

iRIC
iRIC Seminar in Colombia, 17-20 OCT, 2017, Bogota
iRIC version2.0
www.i-ric.org International
River Interface Corporative

iRIC

69
iRIC version2.0
www.i-ric.org International
River Interface Corporative

iRIC

70
 iRIC version2.0
www.i-ric.org International
River Interface Corporative
iRIC Seminars in private companies iRIC

K-Company in Nagoya
71
iRIC version2.0
www.i-ric.org International
River Interface Corporative

iRIC

72
 iRIC version2.0
www.i-ric.org International
River Interface Corporative

iRIC

10th Anniversary Symposium

in Sapporo, 16 NOV 2017
iRIC version2.0
www.i-ric.org International
River Interface Corporative

iRIC
 ➢ More than 40 times iRIC classes in the world in the
past 10 years
➢ 45,000 downloads in the past 3 years, and more
than 20,000 users in the world
JP and EN SITE
45000
40000
35000
30000
25000
20000
15000
10000
5000
0

Download numbers of iRIC GUI and solvers manuals(after 2014)

iRIC user registration numbers

79
 Introduction to
Nays2DH in iRIC

80
 What is Nays2DH?
• Numerical solver for 2D river flow and bed evolution

iRIC Software
Input Data FaSTMECH
Morpho2D
SToRM
Interface
Pre- and Post- Solvers
Output processing
Nays2DH
Data Nays2D Flood
Nays Cube, etc.
 Features of Nays2DH Solver

Original Unsteady + Bed

Some
2D Flow features for
Evolution
inundation
Bank
Nays2D_Flood
Erosion
New
Mixture Grain-
Confluence
size Sediment
Model
Model
 Pre-processing (1)

Import and Setting

Grid Run
Editing Hydraulic etc.
Generation Analysis results
Geographic data Conditions

84
K.Asahi, Y. Yoshida, H.Tsunematsu, Y. Shimizu and
J. Nelson(2012), Development of the iRIC Software
for River analysis
 Pre-process : Data Set
➢ Following data set can be treated in iRIC software

Geographic and Geo-surface data

Geographic Surface

River Channel Flood plain Roughness Vegetation Obstacle etc.

River survey data DEM Photo, Map etc. Polygon

（*.riv） （*.tpo,*anc,*txt） （*.jpg,*.png,*tif）

85
K.Asahi, Y. Yoshida, H.Tsunematsu, Y. Shimizu and
J. Nelson(2012), Development of the iRIC Software * Detail format is opened in http://i-ric.org
for River analysis
 Pre-Processor :
iRIC software can read and edit DEM data

Number of Data

ASTER
http://www.gds.aster.ersdac.or.jp
/gds_www2002/index_e.html

SRTM
(http://srtm.csi.cgiar.org/
SELECTION/inputCoord.asp

X Y Elv.
Copyright (C) 2015 iRIC Project . All Rights Reserved
Pre-Processor :
iRIC software can read and edit DEM data

• Visualization of DEM data of SRTM(Around

Bangkok)

Copyright (C) 2015 iRIC Project . All

Rights Reserved
Pre-Processor :
iRIC software can read and edit
River survey data

Distance Left Bank mark Right bank mark

Mark （X,Y） （X,Y）

Number of points in the

Distance mark index1 index2 index3 index4
cross-section

Distance from the left bank, riverbed elevation

： ：
※index1 – index4
Index1:Number that corresponds to the left bank distance post
Index2:Number that corresponds to the left bank shoreline post
Index3:Number that corresponds to the right bank shoreline post
Index4:Number that corresponds to the right bank distance post
Pre-Processor :
iRIC software can read and edit River survey data

• Visualization with cross-section profile

Copyright (C) 2015 iRIC Project . All

Rights Reserved
Pre-Processor :
iRIC software can read and edit
River survey data

transverse direction
Pre-Processor :
iRIC software can visualize River survey data

• You can edit the river survey data on following

screen

Copyright (C) 2015 iRIC Project . All

Rights Reserved
Pre-Processor :
Creating grids from River survey data

Completion
Pre-Processor :
Creating grids by Polygonal line and width function

DrawCompletion
Center
initial
line
Pre-Processor :
Creating grids for experimental flume etc.

Completion
Pre-Processor :
Creating grids for rectangular area.

Specify Select gridInitial

Completion
parameters generation area
Pre-Processor :
Creating grids for experimental flume etc.

Draw center line and specify

Select
Select
grid
low-water
generation
channel
area
parameters
Initial

Completion
 Grid Generator
Structure grids

Un-structure grids
Straight channel

Sine-generated curve

Kinoshita Meandering curve

Pre-Processor :
Aerial photos and maps can be used as background image
• Aerial photo and map can be used as
background image
 Importing background image from
Geographical Survey Institute of Japan
Import background image from Google
Pre-Processor :
Polygon can be used to specify roughness,
vegetation, obstacles, un-erodible bed etc.
Pre-Processor :
Attributes mapping system
 Pre-process : Mapping
Geographical information

Geographic

Roughness

Vegetation

： above attributes are mapped into each grid nodes and cells
Grid

104
K.Asahi, Y. Yoshida, H.Tsunematsu, Y. Shimizu and
J. Nelson(2012), Development of the iRIC Software
for River analysis
 Unique Feature of Nays2DH

• HOT START
- Users can re-start calculations at arbitrary time-steps.
- This feature helps users to set time steps to be very small
only for specific calculation periods.

1. Output results in the specific Mild Change Rapid Change Mild Change
format
2. Set time steps shorter then use
HOT START function
Discharge

Time
 Cell Attributes in Nays2DH

Vegetation density , Theshio River, Hokkaido Japan

Beige：vege. Density medium
Red ：vege. density coarse

• Various attributes can be added to

every cell.
- Obstacles to block flows
- Vegetation
- Flood plain and low flow channel, etc.
107
 Output
➢ Various
➢ Calculation results can be output various format as text, binary, image, etc.

Calculation Result

Cgn(s) Vtk Image data

Csv format Kml data
（*.csv） format format （*.jpg,*.png,*.bmp,*ppm,
（*.kml）
*.xbm）
（*.cgn） （*.vtk）

Visit Image
Excel Tecplot Viewing
Paraview Google
Text editor MATLAB and Editing
: Earth
: Mathmatica Tools

108
K.Asahi, Y. Yoshida, H.Tsunematsu, Y. Shimizu and
J. Nelson(2012), Development of the iRIC Software
: :
for River analysis
 International
Example of iRIC output (Appetsu River) River Interface Corporative

iRIC
Vector visualization
Nays2D : 2-D flow and morphodynamics

Copyright (C) 2015 iRIC Project . All Rights Reserved

 International
River Interface Corporative

iRIC

Flooding Simulation of Appetsu River in

Japan, SEP 2003
 International
River Interface Corporative

iRIC

Visualizations of calculation results

113
114
115
Alatna River, Allakaket, AK,
United State

116
117
118
119
120
121
 Examples of
Nays2DH and Nays2Flood

122
K.Asahi, Y. Yoshida, H.Tsunematsu, Y. Shimizu and
J. Nelson(2012), Development of the iRIC Software
for River analysis
 Effects of
Vegetation
Ｂ 13KP

14KP
Azuma bridge
Mizuho bridge

Minami ichijo bridge

Ａ Repaired Forest Ｃ

Repaired Forest Ｂ

Repaired Forest A

123
札幌市豊平川 樹木の繁茂度の違いによる洪水時の流況比較
124
2011.9.6 The Otofuke River, Hokkaido, Japan
125
 Channel changes by a flood
Flow
Before the flood (2011)

After the flood (2012)

Chute Cutoff
Destroyed river embankment
Initial Bed Elevation

Chute cutoff

Bank erosion

- Grid size is about 5m

- Based on Laser Profiler data
Initial Bed Surface Status
Vegetated

Unvegetated

- Based on Aerial Photo in 2010

Computational Grid

The destruction point of

the river embankment
- Grid is 512 x 64
- Average size is about 5m
Destruction of Vegetation Area

• Vegetation is removed, when bed is eroded 0.5m

from the initial bed.
Hydrograph
Nasy2DH Computation
Flow, bed and bank evolution by Nays2D

133
Application: Kokuba River, Naha, Japan

RCEM2 134
Grid generation of the confluence
1. Grids of main and branch streams are generated
independently.
2. Export grid files as “.csv” format and combine
them by a spread sheet application.
136
Visualize the result on Google Earth
Mixing Flow an International

Confluence
River Interface Corporative

iRIC

138
Channel Erosion and Bridge Corruption at Chubetsu River
AUG 2011

Damaged Bridge Google Earth, 2010.9.25 download

Site-A

Site-B

Damaged Road

2011年8月24日撮影北海道開発局提供 2011年8月24日撮影北海道開発局提供

139
 600
Discharge Hydrograph
Simulated Result AUG 2011, Chubetsu Riv.

By Satomi Kawamura 400

Chubetsu Riv. Start 8/23 22:00

Discharge(m3/s)
Peak 8/24 4:20
6hr. 20min. After staring
200

Piukenai Riv.
0
0.00 5.00 10.00 15.00 20.00
Time(hr)

Depth
H18 航空写真

140
Observed
Erosion

Brige Damage
Breach - Erosion
- 2m Degradation

From Surveyed Data by Hokkaido Prefectural Governmaent

Simulated
Erosion
Over 0.5m

Bank Erosion
New Channel

Over 2m

H18 航空写真 141

Real Scale Levee Breach Experiments by HDB and
CERI, 2010

Tokachi River, Chiyoda Flood Way

142
 International
River Interface Corporative

iRIC
144
 (Zenno，2011)
145
Grate Flood of Mississippi River 2011

146
Structure of Floodway

Water will come

inside from Inflow
CrevasseHoweve
r usually there are
no gate at there.
So for the
execution, they
have to blow up
the levee.
147
May 5th
(Estimation)
63 feet

62 feet
May 2nd
61 feet (18.6m)

60 feet

59 feet

Execution of Floodway
148
Satellite Image

2011.4.29 2011.5.3

149
 Bomber Cell

On this calculation I define the

cell which call “Bomber
Cell” .First, Bomber Cell behave
as Obstacle Cell. However after
the time has set, water can path
Grid
freely.
150
151
 10m

0m

After 2days
By comparing before and after the execution, it is obvious that water
level has decreased.The red zone becomes smaller which water depth
152
Comparison

10m

153
0m
Comparison

10m

154
0m
The flood simulation in 2011 Chao Phraya Flood by Nays2DFlood
By Toshiki IWASAKI,
Hokkaido univercity

For 3 months

Conditions
➢ Grid mesh:
250m
➢ No. of grid:
25,000
➢ Roughness:
Maninng’s n
= 0.03

155
156
157
158
159
160
Koshi Barrage
Built in 1968
L=1,1150m
52 Gates

161
Chakraborty et al.
2010
162
The channel avulsion
in 2008

163
Channel avulsion, barrage and sedimentation
Sediment supply
from Himalaya

The Koshi
barrage

Channel avulsion

↓
Sedimentation of fine
sediment on floodplain

164
The sedimentation on the floodplain

165
2008
Depth Erosion/Deposition

166
Ｓｉｍｕｌａｔｉｏｎ ｏｆ Ｋｏｓｈｉ Ｒｉｖｅｒ Ｆｌｏｏｄ ｉｎ ２００８

167
 International
River Interface Corporative

iRIC

2015 Myanmar Floods

Image cited from Adventist.org; the Official website of the Seventh-day Adventist world church
https://news.adventist.org/en/all-news/news/go/2015-08-04/myanmar-floods-damage-32-adventist-churches-adventists-assist-flood-victims/
 2015 Myanmar Floods International
River Interface Corporative

• Severe flooding in Myanmar began in July 2015 and continued into iRIC
September, resulting in about 103 deaths and up to 1,000,000 people
affected Movie cited from Reuters
 International
River Interface Corporative

Target basin iRIC

Chindwin River
(largest tributary of
the Irrawaddy)

Irrawaddy

Myittha River
 International
River Interface Corporative

Calculation condition
• <-- topography data created from SRTM data
• Calculation grid size = 90 m
iRIC

• Manning's roughness n = 0.03

• Input hydrograph: Nakayasu's Unit synthetic analysis

85
Elevation [m] 1200
0

Calculation domain
1300

Myittha River

Peak: 2700m3/s
Duration: 8 days
 Results
International
River Interface Corporative

• Comparison of inundated area with calculation results iRIC

Background image: Myanmar: Flood Inundated Area in Kale Townships, Sagaing Region (as of 29 Jul. 2015 12:11)
http://reliefweb.int/map/myanmar/myanmar-flood-inundated-area-kale-townships-sagaing-region-29th-july-2015-1211-am-local
 International
River Interface Corporative

Export to Google Earth iRIC

Kimura and Kang(2018) 175
Result of driftwood deposition with root effect in simulation
Large discharge / Low slope

S=0.0045, Q=0.00100 m3/s

Root part
S=0.0045, Q=0.0010m3/s

Large discharge / High slope

S=0.0070, Q=0.00110 m3/s

S=0.0070, Q=0.0011m3/s
Drift Wood Simulation by I. Kimura(2018)
 International
River Interface Corporative

iRIC

178
 International
River Interface Corporative

iRIC
 International
River Interface Corporative

iRIC

180
 International
River Interface Corporative

iRIC

181
 International
River Interface Corporative

iRIC

Simulation of Free Meandering, with Parker’s chunk block model.

182
 International
River Interface Corporative

iRIC

183
 International
River Interface Corporative

iRIC

184
Great East Japan The Abukuma River
Abukuma River Basin
Earthquake, L=239km
11 MAR 2011 A=5,390km2

Fukushima
20km Nuclear Plant

30km
Cesium 137
Deposit (/m2)

60km

80km
 Sediment size dependency to the density of
radionuclides

Fine sediment has radionuclides of high density

→ The radionuclides will be mainly transported by the transport of
washload and suspended load. 186
187
 Cesium137 Deposition Rate

Cesium 137 (Bq/m2)

Kurikami
River

188
 Cesium 137
Deposit (KBq/m2)
0.05 2.5 5
137Cs(kBq/m2)

Computational result
189
Observation by the helicopter
 Rikuzen Takada City International
River Interface Corporative

iRIC
Before

Pine Tree Foresets

190
資料：東北開発建設部
 International
River Interface Corporative

iRIC
After

2,000 people, 10% of the total population is killed or missing

in this city.
191
資料：東北開発建設部
 International
River Interface Corporative

iRIC

192
 津波遡上の計算 International
River Interface Corporative
釜石市 iRIC
8
(m) 津波水位
6

4
Wave data
2

0
0 1000 2000
Time(s)

地形 粗度 (Manning’s const.)
193
0 10(m) 0.020 0.070
 The simulation of flooding caused by Tsunami runup
International
River Interface Corporative
by S. Kawamura Water Depth 0.0 iRIC
7.0(m

194
iRIC-ELIMO
Easy-performable Long-wave Inundation Model
By Yasunori Watanabe (Hokkaido Univ.)

2011東北地方太平洋沖地震津波 • 手軽に本格的津波計算
• 任意の既往・想定津波を生成
• フリーの海底地形データを使用可能
• さくさく可視化

2003十勝沖津波
 ERIMO by Dr. Yasuharu Watanabe
(Tsunami Model)
• Nonlinear long-wave computation in
global spherical coordinates
• Reliable computing methods with
high-order accuracy
• Polished open boundary conditions:
minimize reflection at the boundary
W z

v
u

Greenwich R
meridian

f y
Download Bathymetry data from
l Japan Oceanographic Data
Center
x Equator

Bathymetry format conversion to

iRIC compatible dataset

spherical coordinate
196
Grid generation on iRIC
 Fault Model Parameters
Input fault parameters for a rectangle
fault model as an initial condition of
tsunami
1. Longitude, latitude of the fault top
2. Rupture length
3. Rupture width
4. Slip length
5. Strike angle
6. Dip angle
7. Slip angle

North

Strike angle
Dip angle
Slip length
Slip angle

197
既往津波 （2011 東北津波）
 想定津波 （北海道500年間隔津波）

北海道500年間隔津波 入力断層

北海道500年間隔津波 北海道庁計算結果
Oiranbuchi in Toyohira River Application of Nays2D Computation
（2015/7/29）
with multi UAV in a bedrock river

20
 調査区間

Toyohira River

Survey Area
Length about800m
Width about100m

Knick Point

20
 Supported by
Hokkaido Consultant Corporation

Precise survey

Ground Area Water Area Supplement

ADCP
GPS
Laser survey

Sonar survey

202
3D data made by Composing all data

精密測量の成果

各測量から得た標高データを合成して作成した鳥
瞰図

20
Simple survey with UAV

We removed this point

because it broke
Setting reference points
at these red points

Phantom3

20
3D data made by UAV Survey
（ moved in
CloudCompare ）

Added about 500 pictures

Into PhotoScan to create 3D data

20
Comparing laser survey and UAV survey

In forest area and upstream where

reference point was removed, the distance
between 2 data is large.

The exposed bedrock area,

the target of this survey, the
distance is about ±15cm.
It’s almost same number.
Yellow and Red…UAV survey data is smaller than Laser data
Blue… UAV survey data is larger than Laser data
Green…almost same number
20
20
 The distance of
knickpoint migration

I made longitudinal sectional view at same

place of knickpoint of Oiranbuchi in iRIC.

２００６ ２０１５

20
209
21
 particles circle in large
depth area.

Water velocity is large at

Knickpoint area

Simulation（iRIC Nays2DH , Depth, particle)

Data：2015Laser＋Sonar1m mesh
Discharge：8m³/s(constant)
21
 When Water quantity is large,
Water flows to right side in downstream area.

Weathered mudstone

Simulation（iRIC Nays2DH , Depth, particle)

Data：2015Laser＋Sonar1m mesh
Discharge：100m³/s(constant)
21
Kinu River Flood of 2015

By. Y.Ishida(2015)
 Akatani River in Kyushu, 2017

By. Y.Ishida(2017)
 Levee Breach and Inundation in Sorachi River, 2018

By. Y.Ishida(2016)
Land Slide and Debris Flow in Hiroshima
2018 撮影：
株式会社パスコ/
国際航業株式会社
2018年7月10日(火)

Simulated
Maxiumum
Deposition

By H.Takebayashi
 iRIC Application Examples in
2016 Hokkaido Typhoon Disaster

画像提供RIC
 August 2016 in Hokkaido

Abashiri
網走
留萌
Rumoi
Kitami
Asahi
旭川
kawa 北見

Typhoon #7
台風7号
th AUG
178/17上陸
小樽
Otaru
札幌
Sapporo
帯広
Obihiro 釧路
Kushiro

台風11号#11
Typhoon
8/21上陸
21th AUG
室蘭
Muroran
Typhoon #10 Typhoon #9
台風9号
31th AUG
台風10号 23th AUG
8/23上陸

函館
Hakodate
Precipitation between 1:00 15th AUG and 24:00 24th AUG

First three Typhoons

 Animation taken from NIED Website
(National Research Institute for Earth Science and Disaster)
Precipitation between 1:00 29th AUG and 9:00 31th AUG

Kushinai
29th-31th AUG
515mm
 Outline of the Ishikari River
Russia - Catchment area: 14,330 km2
China
- Main channel length: 268 km
Urbanized area, Asahikawa City
Main Catchment
Rank River
North
Prefecture area (km2)
Korea

1 Tone Chiba etc. 16,840

South
Korea
Japan 2 Ishikari Hokkaido 14,330
3 Shinano Niigata etc. 11,900
4 Kitakami Miyagi 10,150
5 Kiso Mie 9,100
Hokkaid 6 Tokachi Hokkaido 9,010
o 7 Yodo Osaka 8,240
8 Agano Niigata 7,710
9 Mogami Yamagata 7,040
10 Teshio Hokkaido 5,590

Urbanize area, Sapporo City

Characteristic features Other regions

in Hokkaido
- Population in the Ishikari River 43%
Ishikari River
Basin: More than 3 million Basin
57%
- More than half the population of
Hokkaido live here.
 形 式 アーチ式コンクリートダム
堤 高 102.5m
Taisetsu Dam （1975)
堤 長 305m
堤 体 積 285,000m3

形 式 越流頂式　正面越流型

洪 水 吐
ｹﾞｰﾄ形式･数 鋼製ﾛｰﾗｹﾞｰﾄ～5門
設計洪水流量 1,000m3/s
Mikasa Ponbetsu Dam

常用洪水吐
形 式 管路式　高圧型
(under construction)
ｹﾞｰﾄ形式･数 ﾊｳｴﾙ･ﾊﾞﾝｶﾞｰﾊﾞﾙﾌﾞ～2門
設計洪水流量 140m3/s Chubetsu Dam （2007）
総 貯 水 量 47,100千m3

貯
有効貯 水量 37,100千m3

水
Takisato Dam （2001） 湛 水 面 積 1.5km2
利 用 水 深 37.2m

池
Sunagawa reservoir (1995） 標準断面図

New Katsurazawa Dam

(redevelopment; under
construction)

全景
Kitamura retarding
(under
construction） Present Katsurazawa Dam
Jozankei Dam （1991）
(1957)

Kanayama Dam (1967)

Hoheikyo Dam （1972）
redevelopment
Today, 11 dams and 8 retarding
basen have been built (including
those under construction).
Chitose R., 6 retardings
(under construction）
Izarigawa Dam (1980) Yubari Shuparo Dam (2014)
Chubetsu Dam
25th AUG 2016
Taisetsu Dam
24th AUG 2016
Reservoir operation figure of the Taisetsu Dam

Rain

Water EL

Inflow Discharge

Outflow Discharge
 Kanayama Dam
16:16 31st AUG 2016
 Kanayama Dam Operration Figure, 29th to 31th AUG
5,390 million m3 stored, 1,320 m3/s reduced

0 0 0
時間雨量(mm/h)

累加雨量(mm)
10 10 100
8/28 18:00～9/1 3:00
20 20 200
Total Rainfall 272.9mm,
30 30 300
600 1600 350

400 1400
流入量 Inflow
流入量
200 1200
流入量、放流量（m3/s）

000 1000

貯水位（ｍ）
800 800 340

600 600
貯水位 貯水位
400 400
Water Elevation
of reservor 放流量
Outflow
放流量
200 200

0 0 330
21時

21時

30日 12時
15時

31日 12時
15時

15時

18時

15時

18時
29日 12時
15時

18時

30日 12時
15時

18時

18時

21時

18時

21時

31日 12時

21時

1日 12時

21時

1日 12時
15時

12時
0時

0時

0時

3時

0時

3時

6時

9時

6時

9時
3時

6時

9時

3時

6時

9時

6時

9時

6時

9時

0時

3時

0時

3時
29日

※同時刻の流入量と放流量の
差

※本資料の数値は速報値であるため、今後の調査で変わる場合があります。
 Inundation of the Sorachi River

Levee Breach
Levee Breach
L堤防決壊
= 300m

Breach R38
堤防決壊
L = 150m
富良野 R38
太平橋

空知川

8月31日 6:30 撮影
 Upstream levee
Upstream levee breach point breach

上流側破堤の決壊口上流の写真
撮影：北海道豪雨災害調査団 撮影日：9/1

232
Downstream levee breach
(from Taihei bridge)
Inundation area
Flood flow turned over
the revetment
levee

233
 Discharge Hydrograph at Ikutora Station
1600.00
1400.00
1200.00
Discharge [m3/s]

1000.00
Computed using Distributed
800.00 Hydraulic Model by HDB

600.00
400.00
200.00
0.00
2016/08/29 00:00

2016/08/30 00:00

2016/08/31 00:00

2016/09/01 00:00

2016/09/02 00:00

2016/09/03 00:00

2016/09/04 00:00
236
Numerical Simulation of AUG 30th Flood
Using Nays2DFlood of iRIC

Computation by Y. Yoshida
計算 石田さん HTBイチオシスペシャル(2016年10月1日放送)
Velocity：About 4 [m/s]

8/31 0:40
Downstream levee
breach happened

239
240
 Main Catchment
Rank River
Prefecture area (km2)
1 Tone Chiba etc. 16,840
2 Ishikari Hokkaido 14,330
3 Shinano Niigata etc. 11,900
4 Kitakami Miyagi 10,150
5 Kiso Mie 9,100
6 Tokachi Hokkaido 9,010
7 Yodo Osaka 8,240
8 Agano Niigata 7,710
9 Mogami Yamagata 7,040 Dendritic River
10 Teshio Hokkaido 5,590
Shintoku Town, the Pankeshintoku River
Shimizu Town, Pekerebetsu River
新得町 神社橋
平成28年8月20日からの大雨及び台風10号の大
雨による主な橋梁被災について

【市道】御料十二線橋（名寄市） 【人道橋】江差牛吊橋（上川町） 【町道】菊水12線吊橋（上川町） 【町道】新生橋（遠軽町）

天塩川水系風連別川 石狩川水系石狩川 石狩川水系石狩川 【町道】湯の沢橋
【町道】武利橋（遠軽町）
湧別川水系支湧別川
湧別川水系武利川
【町道】９線橋（美瑛町） 【国道273】高原大橋（上川町）
石狩川水系辺別川 石狩川水系石狩川 【市町村道】厚和橋（北見市） 【道道】光雲橋（鹿追町）
常呂川水系無加川 十勝川水系シイシカリベツ川

【町道】阿満一号橋（美瑛町）
石狩川水系美瑛川 【町道】西上橋（士幌町）
十勝川水系音更川

【町道】神社橋（新得町）
十勝川水系パンケ新得川
【町道】水楽橋（美瑛町）
石狩川水系美瑛川
【国道38】太平橋（南富良野町）
空知川水系空知川
【町道】両泉橋（美瑛町）
石狩川水系美瑛川
【市町村道】北落合橋（南富良野町）
空知川水系シーソラプチ川
【道道】風雲橋（音更町）
十勝川水系然別川
【国道274】新清橋（清水町）
十勝川水系ペケレベツ川
【国道38】元村橋（芽室町）
十勝川水系ピウカ川
【国道274】千呂露橋（日高町）
【国道274】大颱橋 （日高町）
【国道274】日栄橋 （日高町） 【道道】中島橋（帯広市）
【国道274】岩瀬橋 （日高町） 十勝川水系戸蔦別川
【国道274】清瀬橋 （日高町）
【国道274】ニセクシュマナイ橋 （日高町）
【国道274】新岩魚橋（日高町） 【道道】上美生橋（芽室町）
【国道274】鹿鳴橋 （日高町） 十勝川水系美生川
沙流川水系沙流川

【市町村道】戸蔦大橋（中札内村）
【市町村道】更生橋（占冠村） 十勝川水系戸蔦別川
鵡川水系

【市町村道】明星橋（帯広市）
【市町村道】トマム7線橋（占冠村） 【市町村道】千鳥橋（帯広市）
鵡川水系 十勝川水系岩内川

【国道38】小林橋（清水町） 【市町村道】上芽室橋（芽室町） 【町道】石山橋（清水町）

【国道274】日勝大橋（清水町） 十勝川水系小林川 【市町村道】日新橋（芽室町） 【町道】ペケレベツ橋（清水町） 【道道】ヌビナイ橋（大樹町）
十勝川水系 十勝川水系芽室川 【町道】錦橋（清水町） 歴船川水系ヌビナイ川
【町道】新錦橋（清水町）
【町道】地蔵橋（清水町）
【町道】シラウ橋（平取町） 【国道38】清見橋（清水町） 【国R38】芽室橋（芽室町） 【町道】円山橋（清水町） 【国道236】熊の沢橋（広尾町）
沙流川水系シラウ川 十勝川水系ペケレベツ川 十勝川水系芽室川 十勝川水系 豊似川水系

※本資料は速報値であるため、今後の調査で変わる可能性があります。
 北海道開発局提供
国道２７４号 日高町千呂露橋（H28.8.31の状況）
 北海道開発局提供
国道２７４号 日高町千呂露橋（H28.8.31の状況）
 北海道開発局提供
国道２７４号 日高町岩瀬橋（H28.9.2の状況）
 北海道開発局提供
国道２７４号 日高町三国の沢覆道（H28.9.3の状況）
 北海道開発局提供

国道２７４号 清水町日勝峠８合目付近（H28.9.2の状況）
NHK北海道スペシャル(2016年9月30日放送)
The Pekerebetsu River
Photo by RIC
Pekerebetsu River
1㎞

Flow Direction Google

29 JUL 2014

Photo Hokkaido
1 SEP 2016

256
Pasco Co.
1 SEP 2016 257
 400m

Flow Slope Gap

Photo
SEP 1 2016

Slope 1/48 Slope 1/66

258
 Longitudinal Profile
Flow Direction

Slope
1/48

Bridge
Channel Bed Channel Bed
After Flooding Before Flooding

Degradation Damaged
Bridges
About 1㎞

259
 Channel Bed
Before Flooding

After Flooding

Bed Degradation
Channel Bed
Before Flooding

After Flooding
Bedrock
Exposure
 Downstream Part
400m

Flow

263
 Deposition at
Channel Center

Before

After

264
Sediment Budget (Erosion/Deposition) (m3)

3000

Volume（m3/50m）
河道堆積
Channel Deposit
-2000
河道浸食
Channel Erosion

側岸堆積
Bank Deposit

-7000 側岸浸食
Bank Erosion

1/48 1/66
-12000
SP

265
Computation “without” sediment transport
and bed deformation

266
Computation “with” sediment transport and bed deformation
Inundation

Straighten

Erosion

Channel Sedimentation
267
Bank Erosion
Flow Direction

268
Flow Direction

269
Otofuke River
at Kamishihoro Town
 Levee Breach in the Otofuke River
Kamishihoro
N
Google Earthによる最新画像（撮影日時不明)

300m

北開水工コンサルタント伊賀氏による撮影画像（2016年9月13,14日撮影)

Flow Direction Levee

Breach
Otofuke
River
 After 2011 Flood
After 2016 Flood

音更川の上士幌町平和
四十号橋付近より上流を望む
 Time

Straight Slightly Strongly

Meander Meander
275

2011.9.6 The Otofuke River, Hokkaido, Japan

 Destruction of Vegetation Area
• Vegetation is removed, when bed is eroded
0.5m from the initial bed.
Flow, bed and bank evolution by Nays2D
芽室川の河床変動
 The iRIC software
Can be downloaded from following url
http://i-ric.org/

Thank You
280
K.Asahi, Y. Yoshida, H.Tsunematsu, Y. Shimizu and
J. Nelson(2012), Development of the iRIC Software
for River analysis