Climate change impact on flood hazard over Italy
Adriano Fantini
3rd year PhD Course in Earth Science, Fluid Dynamics, and Mathematics
Supervisor: Erika Coppola
ADRIANO.FANTINI@phd.units.it
Aims
- Flood hazard mapping over Italy
- Scientific, reliable approach
- Future projections
Models
- ICTP RegCM Regional Climate Model
- CHyM hydrological model
- CA2D hydraulic model
Project overview
Participants
- Erika Coppola
- Rita Nogherotto
- Filippo Giorgi
- Adriano Fantini
- Francesca Raffaele
Methodology
Precipitation:
- Observations
- RCM output
Gridded netCDF:
- River network
- Discharges
hydrological model
For each RP, cell:
- Gumbel distr.
- Hydrographs
- Extreme Q
Statistical analysis
For each RP, cell:
- Flood extent
- Flood depth
(multiple simulations)
- RCM output
- Discharges
- Floods
Validation and change for
CA2D hydraulic model
Based on Maione et al., 2003
(over nine domains)
1.0 - Flood hazard: overview
RISK
=
HAZARD
×
EXPOSURE
×
COPING FACTOR
Event probability, Return Period
Goods, people and services exposed
Emergency plans, adaptation strategies
1.1 - Flood hazard: overview
Methods used to obtain hazard maps:
- Historical discharge and flooding records
- Documentary evidence from past events
- Local surveys
- Hydrological + hydraulical modelling
Advantages:
- Possible on any region/basin, also large scale
- Extreme value analysis extends to any Return Period
- Climate projections can drive the models (future change)
- Uncertainty can be assessed via ensemble modelling
PRECIPITATION
DISCHARGE
FLOOD
HYDROLOGICAL MODEL
HYDRAULICAL MODEL
1.2 - Flood hazard: overview
Flood hazard over Italy (PRESENT DAY)
ISPRA (2018), Dissesto idrogeologico in Italia: pericolosità e indicatori di rischio
Return Period:
20-50 yrs
Return Period:
100-200 yrs
1.3 - Flood hazard: overview
Flood hazard over Italy (FUTURE CHANGE)
- Several studies available
- No study specific to Italy
- Relatively low resolution (>5km)
- No flood extents, only flood proxies (extreme discharge)
Results over Italy:
- General increase in flood proxies by the end of the century
- Especially for Northern Italy
- Low resolution does not resolve smaller basins
Rojas et al. (2012); Thober et al. (2018); Donnelly et al. (2017); Alfieri et al. (2015)
2.0 - The GRIPHO dataset
Two main goals:
- Driving the hydrological model with high resolution data
- Validating the regional climate simulations
Data provided by Marco Verdecchia (CETEMPS):
- 2001 - present
- Hourly
- 3712 stations
- No quality check
GRIPHO
(GRidded Italian Precipitation Hourly Observations)
2.1 - The GRIPHO dataset
Issues:
- Varying station availability and density across the domain
- Inconsistencies, outliers and data errors
- Only ~15 years or data
Strong points:
- High station density (~1 station/100km²)
- High temporal resolution (1 hour)
- Only high resolution station-based dataset covering all Italy
FIRST-STAGE FILTERING
FLAGGING
FLAG CHECKING
MANUAL DATA CLEANING
CLEANED DATASET!
2.2 - The GRIPHO dataset
Cleaning procedure
2.3 - The GRIPHO dataset
Gridding
- 12km Lambert Conformal Conic grid
-
Gridding method based on Delaunay polygons using SciPy's interpolate.griddata
- Simple, fast method which minimizes smoothing
- Similar to Norway's KLIMAGRID dataset (Mohr, 2008, 2009)
- NetCDF CF-compliant output format
Mohr M., 2008: New Routines for Gridding of Temperature and Precipitation Observations for seNorge.no
Mohr M., 2009: Comparison of versions 1.1 and 1.0 of gridded temperature and precipitation data for Norway
Velasquez N. et al., 2011: Rainfall distribution based on a Delaunay triangulation method
2.4 - The GRIPHO dataset
Validation (mean PR, detail)
2.5 - The GRIPHO dataset
Validation (R95ptot, detail)
2.6 - The GRIPHO dataset
Conclusions
- Performance in the North similar to other high resolution datasets (EURO4M-APGD, ARCIS)
- Much finer details compared to HMR and E-OBS, especially for extremes
- Only station-based dataset over Italy providing hourly precipitation
Fantini A., Coppola E., Verdecchia M. and Giuliani G.:
‘GRIPHO: a gridded high-resolution hourly precipitation dataset over Italy’, in preparation
Issues:
- Further data cleaning must be performed (esp. in Sicily)
- Some areas are left missing in some time periods, could be filled in?
Some examples...
Outliers
2.x - The GRIPHO dataset
2.x - The GRIPHO dataset
Validation
2.x - The GRIPHO dataset
Validation
2.x - The GRIPHO dataset
Validation (PDF, detail)
2.x - The GRIPHO dataset
2.x - The GRIPHO dataset
Automated filtering
- Removal of extreme values > 200 mm/h
- Removal of isolated > 100 mm/h reports
- Removal of complete months with > 1800 mm (2.5 mm/h)
- Removal of continuous identical values
Rauthe et al. 2013; A Central European precipitation climatology–Part I: Generation and validation of a high-resolution gridded daily data set (HYRAS)
Isotta et al. 2013; The climate of daily precipitation in the Alps: development and analysis of a high‐resolution grid dataset from pan‐Alpine rain‐gauge data
Perry et al., 2009; The generation of daily gridded datasets of temperature and rainfall for the UK
Hiebl et al., 2017; Daily precipitation grids for Austria since 1961—development and evaluation of a spatial dataset for hydroclimatic monitoring and modelling
2.x - The GRIPHO dataset
Flagging
- Total % of valid values
- Total % of valid values ≠ 0
- % of events > mean + {10,15,25} SD
- % of events > median + {10,15,25} IQR
- % of events with 50 < pr < 100mm/h
- % of events with pr > 100mm/h
- % of the top {1,5,10} most common values
Manual filtering on a monthly basis
2.x - The GRIPHO dataset
Before | After | Diff | |
---|---|---|---|
Total valid values | 250M | 243M | -2.6% |
Total flags | 324468 | 58008 | -82.1% |
pr > Mean + 20SD | 3240 | 2538 | -21.7% |
pr > Median + 20IQR | 49753 | 22519 | -54.7% |
pr > 100 mm/h | 221822 | 711 | -99.7% |
Other softer flags | 49646 | 32240 | -35.1% |
2.x - The GRIPHO dataset
Results on warning flags
FILTERED
ORIGINAL
3.0 - Regional Climate Model
Two RegCM 4.6.1 12km EURO-CORDEX simulations run on ICTP's Argo and CINECA's Marconi clusters:
- ERA-Interim driven 1979-2016 historical simulation
- HadGEM driven 1971-2099 scenario simulation (RCP8.5)
- 135 3-year tuning experiments
- 6000 runtime hours
- 3 million core-hours
- >100TB disk usage
Thanks to James Ciarlo` for running part of the HadGEM driven simulation!
3.1 - Regional Climate Model
MEAN PR VALIDATION
3.2 - Regional Climate Model
MEAN PR
3.3 - Regional Climate Model
R95ptot
3.4 - Regional Climate Model
R99ptot
3.5 - Regional Climate Model
Conclusions
- Model precipitation and temperature generally in line with observations
- Increased projected average precipitation by the end of the century in winter in the north; decrease in the south and isles in summer
- Precipitation extremes projected to strongly increase
3.x - Regional Climate Model
3.x - Regional Climate Model
PDF change
3.x - Regional Climate Model
PDF change
3.x - Regional Climate Model
R95ptot validation
3.x - Regional Climate Model
R99ptot validation
3.x - Regional Climate Model
TAS
PR
3.x - Regional Climate Model
PDF validation
3.x - Regional Climate Model
PDF validation
3.x - Regional Climate Model
Temperature validation
3.x - Regional Climate Model
Temperature change
4.0 - Cetemps Hydrological Model
CHyM Is a distributed (gridded) hydrological model. Peculiarities:
- Can build DEM from various sources, smoothing by cellula automata algorithms
- Can use several kind of inputs, such as station observations, gridded model data, etc.
- Designed to work on any domain
- Used daily at CETEMPS for operational forecasts
- Hourly NetCDF output
4.1 - Cetemps Hydrological Model
- 9 simulated domains
- Tested several Digital Elevation Models, chose HydroSHEDS
- Specific tuning for each region
- 300-900m resolution
- Argo and Marconi clusters
- 3000 runtime hours
- 100k core-hours
- ~35TB
Three drivers:
- GRIPHO (MM5) 2001-2016
- RegCM-ERA 1980-2016
- RegCM-HAD 1972-2099
4.2 - Cetemps Hydrological Model
4.3 - Cetemps Hydrological Model
MEAN DISCHARGE
4.4 - Cetemps Hydrological Model
MEAN ANNUAL MAXIMUM DISCHARGE
4.5 - Cetemps Hydrological Model
100-YEAR ESTIMATED DISCHARGE
GRIPHO
HadGEM
present
HadGEM
2020-2049
HadGEM
2070-2099
100-year discharges estimated following:
Maione et al., 2003: Regional estimation of synthetic design hydrographs
Thanks to Francesca Raffaele for the research work!
4.6 - Cetemps Hydrological Model
Conclusions
- Model discharges in line with observations, for the few available stations
- Average discharge projected to slightly increase in the north; mixed changes in the rest of the domain/seasons
- Extreme discharges projected to strongly increase by the end of the century (often >1.5x)
- Signal very dependent on the region and season
4.x - Cetemps Hydrological Model
96 river network reconstruction tests for each region
(manual and automatic, with distance measures)
4.x - Cetemps Hydrological Model
CHyM-OP reproduced domains:
4.x - Cetemps Hydrological Model
CHyM-OP reproduced domains:
4.x - Cetemps Hydrological Model
4.x - Cetemps Hydrological Model
So far:
- Performed ~2000 model simulations to find the best configuration for the river network reconstruction
- Identify and compare the reconstruction of the Po river with different metrics: mean distance, basin area, distance Q95...
4.x - Cetemps Hydrological Model
5.0 - CA2D hydraulic model
2D flood inundation model from Dottori and Todini, 2010, 2011, modified by Rita Nogherotto to run in parallel
- DEM
- River network
- SDH
5.1 - CA2D hydraulic model
- 90m resolution
- HydroSHEDS DEM and river channels
- "Virtual stations" cover all domains every 5km
- Normalization with 1.5-year return period flood level
- Still ongoing, currently completed only for CHyM (GRIPHO)
Satellite images from COSMO-SkyMed, November 2016 event
5.2 - CA2D hydraulic model
5.3 - CA2D hydraulic model
Nogherotto R., Fantini A., Raffaele F., Coppola E. and Giorgi F.:
´An integrated hydrological and hydraulic modelling approach for the flood risk assessment over Po river basin: a case study for the ALLIANZ Insurance Company´ (in preparation)
Conclusions
- Encouraging initial results
- Good results in the case study
- General agreement with ISPRA maps
Issues:
- Ignoring coastal flooding, dams and water management
- Lack of data available for validation
- Future flood hazard not yet computed
5.x - CA2D hydraulic model
5.x - CA2D hydraulic model
"Virtual stations"
6.0 - Summary
New data:
- New state of the art GRIPHO precipitation dataset
- Two new RCM EURO-CORDEX simulations
- Three new high resolution CHyM hydrological simulations
Ongoing work and future improvements:
- Complete and analyse all the CA2D simulations, projections
- Perform ensemble analysis to assess uncertainty
Flood hazard:
- We can produce flood hazard maps via a model chain
- Much higher resolution than previous studies
- Scientific, reproducible result
- The methodology can be applied anywhere
- Continent-scale studies are possible
- Strong projected increase in extreme pr and flood proxies
6.x - Visualization
An R/Leaflet tool for flood, river, DEM, basin and station visualization
6.x - Visualization
Click_edit: an R/Shiny tool for WYSISYG editing of NetCDF files
Thanks for your attention!
adriano.fantini@phd.units.it
Precipitation:
- Observations
- RCM output
Gridded netCDF:
- River network
- Discharges
hydrological model
For each RP, cell:
- Gumbel distr.
- Hydrographs
- Extreme Q
Statistical analysis
For each RP, cell:
- Flood extent
- Flood depth
(multiple simulations)
- RCM output
- Discharges
- Floods
Validation and change for
CA2D hydraulic model
Based on Maione et al., 2003
(over nine domains)
ESFM 2017-2018 report (III year)
By odineidolon
ESFM 2017-2018 report (III year)
Presentation for III year of PhD, 20 min
- 539