9th international conference urban drainage...
TRANSCRIPT
![Page 1: 9th International Conference URBAN DRAINAGE MODELLINGhikom.grf.bg.ac.rs/stari-sajt/9UDM/Presentations/238_PPT.pdf · Evaluating the impact of climate change on urban scale extreme](https://reader036.vdocument.in/reader036/viewer/2022081406/5f0daf8b7e708231d43b945d/html5/thumbnails/1.jpg)
Evaluating the impact of climate change on urban scale extreme rainfall events:
Coupling of multiple global circulation
models with a stochastic rainfall generator
Gonzalo Andrés PEÑA CASTELLANOS Assela PATHIRANA
September 5th2012
9th International Conference
URBAN DRAINAGE MODELLING
![Page 2: 9th International Conference URBAN DRAINAGE MODELLINGhikom.grf.bg.ac.rs/stari-sajt/9UDM/Presentations/238_PPT.pdf · Evaluating the impact of climate change on urban scale extreme](https://reader036.vdocument.in/reader036/viewer/2022081406/5f0daf8b7e708231d43b945d/html5/thumbnails/2.jpg)
What’s new
• Attempt to get extreme rainfalls (urban level) right.
• Number of GCMs (12) + Scenarios (3) + Periods (2)
Content
• Introduction
• Methodology
• Results
• Conclusion
![Page 3: 9th International Conference URBAN DRAINAGE MODELLINGhikom.grf.bg.ac.rs/stari-sajt/9UDM/Presentations/238_PPT.pdf · Evaluating the impact of climate change on urban scale extreme](https://reader036.vdocument.in/reader036/viewer/2022081406/5f0daf8b7e708231d43b945d/html5/thumbnails/3.jpg)
1. Introduction
![Page 4: 9th International Conference URBAN DRAINAGE MODELLINGhikom.grf.bg.ac.rs/stari-sajt/9UDM/Presentations/238_PPT.pdf · Evaluating the impact of climate change on urban scale extreme](https://reader036.vdocument.in/reader036/viewer/2022081406/5f0daf8b7e708231d43b945d/html5/thumbnails/4.jpg)
Earth Movement
Solar Radiation
State of Climatic Variables
Over an extended period of time
Climate and Climate Change In
trod
uctio
n
1 / 11
![Page 5: 9th International Conference URBAN DRAINAGE MODELLINGhikom.grf.bg.ac.rs/stari-sajt/9UDM/Presentations/238_PPT.pdf · Evaluating the impact of climate change on urban scale extreme](https://reader036.vdocument.in/reader036/viewer/2022081406/5f0daf8b7e708231d43b945d/html5/thumbnails/5.jpg)
Earth Movement
Solar Radiation
State of Climatic Variables
Over an extended period of time
Climate and Climate Change In
trod
uctio
n
Composition of Atmosphere
(CO2, H2O, CH4) Gas emissions
1 / 11
![Page 6: 9th International Conference URBAN DRAINAGE MODELLINGhikom.grf.bg.ac.rs/stari-sajt/9UDM/Presentations/238_PPT.pdf · Evaluating the impact of climate change on urban scale extreme](https://reader036.vdocument.in/reader036/viewer/2022081406/5f0daf8b7e708231d43b945d/html5/thumbnails/6.jpg)
Modeling Climate Change In
trod
uctio
n
1 / 11
Gas Emissions Scenarios
(SRES)
![Page 7: 9th International Conference URBAN DRAINAGE MODELLINGhikom.grf.bg.ac.rs/stari-sajt/9UDM/Presentations/238_PPT.pdf · Evaluating the impact of climate change on urban scale extreme](https://reader036.vdocument.in/reader036/viewer/2022081406/5f0daf8b7e708231d43b945d/html5/thumbnails/7.jpg)
Intr
oduc
tion
Modeling Climate Change
2 / 11
3.75˚ 417 km
2.5˚ COARSE 278 km SCALE
Global Circulation Models & Regional Models
- Planet - Region
![Page 8: 9th International Conference URBAN DRAINAGE MODELLINGhikom.grf.bg.ac.rs/stari-sajt/9UDM/Presentations/238_PPT.pdf · Evaluating the impact of climate change on urban scale extreme](https://reader036.vdocument.in/reader036/viewer/2022081406/5f0daf8b7e708231d43b945d/html5/thumbnails/8.jpg)
Intr
oduc
tion
Modeling Climate Change
2 / 11
3.75˚ 417 km
2.5˚ COARSE 278 km SCALE
Global Circulation Models & Regional Models
- Planet - Regions - Cities??
![Page 9: 9th International Conference URBAN DRAINAGE MODELLINGhikom.grf.bg.ac.rs/stari-sajt/9UDM/Presentations/238_PPT.pdf · Evaluating the impact of climate change on urban scale extreme](https://reader036.vdocument.in/reader036/viewer/2022081406/5f0daf8b7e708231d43b945d/html5/thumbnails/9.jpg)
The Objective: Urban Scale In
trod
uctio
n
3 / 11
Total Precipitation Flux [Daily Output]
(Extreme) Precipitation [Subdaily]
SPATIAL RESOLUTION Urban Scale
[Point Scale]
Global Scale [200-400 km]
TEMPORAL RESOLUTION
![Page 10: 9th International Conference URBAN DRAINAGE MODELLINGhikom.grf.bg.ac.rs/stari-sajt/9UDM/Presentations/238_PPT.pdf · Evaluating the impact of climate change on urban scale extreme](https://reader036.vdocument.in/reader036/viewer/2022081406/5f0daf8b7e708231d43b945d/html5/thumbnails/10.jpg)
2. Methodology
time spent
task size
wins
loses runs script
does it manually
does it manually
gets annoyed
Makes fun of geek´s Complicated method
writes script to automate
![Page 11: 9th International Conference URBAN DRAINAGE MODELLINGhikom.grf.bg.ac.rs/stari-sajt/9UDM/Presentations/238_PPT.pdf · Evaluating the impact of climate change on urban scale extreme](https://reader036.vdocument.in/reader036/viewer/2022081406/5f0daf8b7e708231d43b945d/html5/thumbnails/11.jpg)
Poisson cluster process: Neyman Scott Rectangular Pulse (NSRP)
time
intensity [mm/hr]
Generation of Synthetic Rainfall Using a Stochastic Process
1.) Generate a random number of storm origins
Met
hodo
logy
4 / 11
![Page 12: 9th International Conference URBAN DRAINAGE MODELLINGhikom.grf.bg.ac.rs/stari-sajt/9UDM/Presentations/238_PPT.pdf · Evaluating the impact of climate change on urban scale extreme](https://reader036.vdocument.in/reader036/viewer/2022081406/5f0daf8b7e708231d43b945d/html5/thumbnails/12.jpg)
Poisson cluster process: Neyman Scott Rectangular Pulse (NSRP)
time
intensity [mm/hr]
Generation of Synthetic Rainfall Using a Stochastic Process
2.) Each storm generates a random number of cells and random cell origins
Met
hodo
logy
4 / 11
![Page 13: 9th International Conference URBAN DRAINAGE MODELLINGhikom.grf.bg.ac.rs/stari-sajt/9UDM/Presentations/238_PPT.pdf · Evaluating the impact of climate change on urban scale extreme](https://reader036.vdocument.in/reader036/viewer/2022081406/5f0daf8b7e708231d43b945d/html5/thumbnails/13.jpg)
Poisson cluster process: Neyman Scott Rectangular Pulse (NSRP)
time
intensity [mm/hr]
Generation of Synthetic Rainfall Using a Stochastic Process
3.) A random duration of each cell is generated
Met
hodo
logy
4 / 11
![Page 14: 9th International Conference URBAN DRAINAGE MODELLINGhikom.grf.bg.ac.rs/stari-sajt/9UDM/Presentations/238_PPT.pdf · Evaluating the impact of climate change on urban scale extreme](https://reader036.vdocument.in/reader036/viewer/2022081406/5f0daf8b7e708231d43b945d/html5/thumbnails/14.jpg)
Poisson cluster process: Neyman Scott Rectangular Pulse (NSRP)
time
intensity [mm/hr]
Generation of Synthetic Rainfall Using a Stochastic Process
4.) A random intensity is generated
Met
hodo
logy
4 / 11
![Page 15: 9th International Conference URBAN DRAINAGE MODELLINGhikom.grf.bg.ac.rs/stari-sajt/9UDM/Presentations/238_PPT.pdf · Evaluating the impact of climate change on urban scale extreme](https://reader036.vdocument.in/reader036/viewer/2022081406/5f0daf8b7e708231d43b945d/html5/thumbnails/15.jpg)
Poisson cluster process: Neyman Scott Rectangular Pulse (NSRP)
time
intensity [mm/hr]
Generation of Synthetic Rainfall Using a Stochastic Process
5.) Total intensity is the sum of the active cells
Met
hodo
logy
4 / 11
![Page 16: 9th International Conference URBAN DRAINAGE MODELLINGhikom.grf.bg.ac.rs/stari-sajt/9UDM/Presentations/238_PPT.pdf · Evaluating the impact of climate change on urban scale extreme](https://reader036.vdocument.in/reader036/viewer/2022081406/5f0daf8b7e708231d43b945d/html5/thumbnails/16.jpg)
Poisson cluster process: Neyman Scott Rectangular Pulse (NSRP)
time
intensity [mm/hr]
Synthetic rainfall
Generation of Synthetic Rainfall Using a Stochastic Process
Met
hodo
logy
4 / 11
Calibrated from: Mean variance coefficient of variation dry spell duration log autocorrelation
![Page 17: 9th International Conference URBAN DRAINAGE MODELLINGhikom.grf.bg.ac.rs/stari-sajt/9UDM/Presentations/238_PPT.pdf · Evaluating the impact of climate change on urban scale extreme](https://reader036.vdocument.in/reader036/viewer/2022081406/5f0daf8b7e708231d43b945d/html5/thumbnails/17.jpg)
Combining available precipitation data Historical Simulation
Future Simulation
?
Observed Historical Data
Met
hodo
logy
Models: 12 GCM Period: 1974-1999 Scenario: Past
Models: 12 GCM Period: 2046-2065 2081-2100 Scenario: SRES A1b SRES A2 SRES B1
Case study: Japan, Kochi Period: 1974-1999
Factor of Change Of Statistics
5 / 11
![Page 18: 9th International Conference URBAN DRAINAGE MODELLINGhikom.grf.bg.ac.rs/stari-sajt/9UDM/Presentations/238_PPT.pdf · Evaluating the impact of climate change on urban scale extreme](https://reader036.vdocument.in/reader036/viewer/2022081406/5f0daf8b7e708231d43b945d/html5/thumbnails/18.jpg)
Combining available precipitation data M
etho
dolo
gy
5 / 11
?
Factor of Change Of Statistics
Bayes Theorem Solved numerically
Markov Chain Monte Carlo MCMC
Probability of change of the the statistic analyzed
![Page 19: 9th International Conference URBAN DRAINAGE MODELLINGhikom.grf.bg.ac.rs/stari-sajt/9UDM/Presentations/238_PPT.pdf · Evaluating the impact of climate change on urban scale extreme](https://reader036.vdocument.in/reader036/viewer/2022081406/5f0daf8b7e708231d43b945d/html5/thumbnails/19.jpg)
Historical State Future State 2000 1970 2030 2060
Now
Statistic Value
Time
Combining available precipitation data
Change Factor
NSRP can be recalibrated to take into account the effect of climate change by including the factor of change that results from the bayesian ensemble
Met
hodo
logy
5 / 11
![Page 20: 9th International Conference URBAN DRAINAGE MODELLINGhikom.grf.bg.ac.rs/stari-sajt/9UDM/Presentations/238_PPT.pdf · Evaluating the impact of climate change on urban scale extreme](https://reader036.vdocument.in/reader036/viewer/2022081406/5f0daf8b7e708231d43b945d/html5/thumbnails/20.jpg)
Methodology description
5. Modified Parameters Future: 2046-2065 Future: 2081-2100
2. Evaluation of Daily Statistics From GCM ensemble
(Daily Statistics >= 24 hours)
3. Bayesian Ensemble: Factor of Change (Extract mean factor of change)
4. Extend factor of change finer scale (Subdaily scale: < 24 hours)
1. Stochastic Rainfall Generator Parameters (Past: 1974-2000)
Met
hodo
logy
6 / 11
6. Run an ensemble of simulations and perform extreme values analysis
(200 simulations per location, scenario and future period)
![Page 21: 9th International Conference URBAN DRAINAGE MODELLINGhikom.grf.bg.ac.rs/stari-sajt/9UDM/Presentations/238_PPT.pdf · Evaluating the impact of climate change on urban scale extreme](https://reader036.vdocument.in/reader036/viewer/2022081406/5f0daf8b7e708231d43b945d/html5/thumbnails/21.jpg)
4. Results
![Page 22: 9th International Conference URBAN DRAINAGE MODELLINGhikom.grf.bg.ac.rs/stari-sajt/9UDM/Presentations/238_PPT.pdf · Evaluating the impact of climate change on urban scale extreme](https://reader036.vdocument.in/reader036/viewer/2022081406/5f0daf8b7e708231d43b945d/html5/thumbnails/22.jpg)
Resu
lts
GCM models Used
Scenarios Used
Periods Used
12 3 2 NCAR CCSM3.0, MRI CFCM2 3.2A, MPI ECHAME 5, MIUB ECHO G, MIROCS 2 MEDRES, IPSL CM4, INMCM3.0, GISS MODEL ER, GFDL CM2.1, CSIRO MK3.5, CNRM CM3, CCCMA CGCM3.1
SRES A1 SRES A2 SRES B1
2046-2065
2081-2100
Final ensemble
Methodology applicable to any location in the world where hourly precipitation series are
available 7 / 11
![Page 23: 9th International Conference URBAN DRAINAGE MODELLINGhikom.grf.bg.ac.rs/stari-sajt/9UDM/Presentations/238_PPT.pdf · Evaluating the impact of climate change on urban scale extreme](https://reader036.vdocument.in/reader036/viewer/2022081406/5f0daf8b7e708231d43b945d/html5/thumbnails/23.jpg)
Extreme value representation by the NSRP
Resu
lts
Good fit for return periods <= 10 years
8 / 11 Under estimation higher for higher aggregation
![Page 24: 9th International Conference URBAN DRAINAGE MODELLINGhikom.grf.bg.ac.rs/stari-sajt/9UDM/Presentations/238_PPT.pdf · Evaluating the impact of climate change on urban scale extreme](https://reader036.vdocument.in/reader036/viewer/2022081406/5f0daf8b7e708231d43b945d/html5/thumbnails/24.jpg)
Resu
lts
9 / 11
Factors of change Bayesian ensemble results
Only the mean factor of change
was used
Kochi 2045-2065. Hourly mean
Computational Intensive process
Parallel processing
was used
![Page 25: 9th International Conference URBAN DRAINAGE MODELLINGhikom.grf.bg.ac.rs/stari-sajt/9UDM/Presentations/238_PPT.pdf · Evaluating the impact of climate change on urban scale extreme](https://reader036.vdocument.in/reader036/viewer/2022081406/5f0daf8b7e708231d43b945d/html5/thumbnails/25.jpg)
Idf Curves Re
sults
10 / 11
![Page 26: 9th International Conference URBAN DRAINAGE MODELLINGhikom.grf.bg.ac.rs/stari-sajt/9UDM/Presentations/238_PPT.pdf · Evaluating the impact of climate change on urban scale extreme](https://reader036.vdocument.in/reader036/viewer/2022081406/5f0daf8b7e708231d43b945d/html5/thumbnails/26.jpg)
Conclusions
![Page 27: 9th International Conference URBAN DRAINAGE MODELLINGhikom.grf.bg.ac.rs/stari-sajt/9UDM/Presentations/238_PPT.pdf · Evaluating the impact of climate change on urban scale extreme](https://reader036.vdocument.in/reader036/viewer/2022081406/5f0daf8b7e708231d43b945d/html5/thumbnails/27.jpg)
Conc
lusio
ns
11 / 11
Conclusions A methodology based on the use of a weather generator in climate impact studies was extended to include several GCMs. Several scenarios and future periods were used to evaluate change in extreme rainfall events at urban scale. Method itself is OK < 10 year events. However, GCM results are all over the place!
The uncertainty in the use of different GCMs output could be assessed by implementing a Montecarlo type simulation (Even more computationally intensive!!!!)
![Page 28: 9th International Conference URBAN DRAINAGE MODELLINGhikom.grf.bg.ac.rs/stari-sajt/9UDM/Presentations/238_PPT.pdf · Evaluating the impact of climate change on urban scale extreme](https://reader036.vdocument.in/reader036/viewer/2022081406/5f0daf8b7e708231d43b945d/html5/thumbnails/28.jpg)
Conc
lusio
ns
11 / 11
Conclusions Large uncertainties still exist inherent to the Bayesian Ensemble approach (assumption of independence between GCMs and the mismatch between the grid cell size)
The worked methodology was applied to a series of GCMs but could be equally used to Regional circulation models (RCMs)
Freely available and open source tools were successfully explored with the additional
benefit of a simple parallelization scheme.
![Page 29: 9th International Conference URBAN DRAINAGE MODELLINGhikom.grf.bg.ac.rs/stari-sajt/9UDM/Presentations/238_PPT.pdf · Evaluating the impact of climate change on urban scale extreme](https://reader036.vdocument.in/reader036/viewer/2022081406/5f0daf8b7e708231d43b945d/html5/thumbnails/29.jpg)
Questions
Authors: Gonzalo Andrés PEÑA CASTELLANOS
Assela PATHIRANA [email protected]
9th International Conference URBAN DRAINAGE MODELLING
? Thank you!
Acknowledgements Dr. Simone FATICHI Dr. Chris FONNESBECK Dr. Abraham FLAXMAN Dr. Damir BRDJANOVIC COLFUTURO Foundation DUPC Publication fund
![Page 30: 9th International Conference URBAN DRAINAGE MODELLINGhikom.grf.bg.ac.rs/stari-sajt/9UDM/Presentations/238_PPT.pdf · Evaluating the impact of climate change on urban scale extreme](https://reader036.vdocument.in/reader036/viewer/2022081406/5f0daf8b7e708231d43b945d/html5/thumbnails/30.jpg)
Extra info Downscaling methods
Working infrastructure For scientific computing
Spatial and temporal scales
GCM models and Scenarios used
Variability of total precipitation in GCM
Stochastic fit
Gas emissions scenarios
![Page 31: 9th International Conference URBAN DRAINAGE MODELLINGhikom.grf.bg.ac.rs/stari-sajt/9UDM/Presentations/238_PPT.pdf · Evaluating the impact of climate change on urban scale extreme](https://reader036.vdocument.in/reader036/viewer/2022081406/5f0daf8b7e708231d43b945d/html5/thumbnails/31.jpg)
Downscaling Method
Statistical Advantages: Comparatively cheap and computationally efficient Can provide local scale climatic variables from GCM-scale output
Disadvantages: Dependent on GCM boundary forcing; affected by biases in underlying GCM
(Fowler, 2007)
Transfer Functions
Weather Typing
Stochastic
Dynamic
Weather Generators (Advanced Weather
Generator AWE-GEN) (Fatichi, 2011)
Method used: Rainfall Generator
Extr
a in
form
atio
n Home
![Page 32: 9th International Conference URBAN DRAINAGE MODELLINGhikom.grf.bg.ac.rs/stari-sajt/9UDM/Presentations/238_PPT.pdf · Evaluating the impact of climate change on urban scale extreme](https://reader036.vdocument.in/reader036/viewer/2022081406/5f0daf8b7e708231d43b945d/html5/thumbnails/32.jpg)
Working Infrastructure Python for Scientific Computing
Spyder
Python
NumPy
Scipy Matplotlib netCDF4 PyMC
Parallel-Python
Packages and subpackages Development Environment
pySide
Graphical User Interface (under development)
Extr
a in
form
atio
n Home
![Page 33: 9th International Conference URBAN DRAINAGE MODELLINGhikom.grf.bg.ac.rs/stari-sajt/9UDM/Presentations/238_PPT.pdf · Evaluating the impact of climate change on urban scale extreme](https://reader036.vdocument.in/reader036/viewer/2022081406/5f0daf8b7e708231d43b945d/html5/thumbnails/33.jpg)
Flooding Urban drainage
Needed resolution
Water supply
Irrigation
GCM/RCM Output
Climate change
River discharge
Erosion
Rainwater harvest
Century Decade
Year Month
Week Day
Hour Minute Second
Spatial scale
Tim
e
Nor
th A
tlant
ic
Osc
ilatio
n
Forn
tal
Syst
ems
Cell
clus
ter
Rain
cel
l
Rain
dro
p
Hydr
olog
y
Data needs
Measurements
Extr
a in
form
atio
n Spatial and temporal Scales Home
![Page 34: 9th International Conference URBAN DRAINAGE MODELLINGhikom.grf.bg.ac.rs/stari-sajt/9UDM/Presentations/238_PPT.pdf · Evaluating the impact of climate change on urban scale extreme](https://reader036.vdocument.in/reader036/viewer/2022081406/5f0daf8b7e708231d43b945d/html5/thumbnails/34.jpg)
GCM Models and scenarios used Home
![Page 35: 9th International Conference URBAN DRAINAGE MODELLINGhikom.grf.bg.ac.rs/stari-sajt/9UDM/Presentations/238_PPT.pdf · Evaluating the impact of climate change on urban scale extreme](https://reader036.vdocument.in/reader036/viewer/2022081406/5f0daf8b7e708231d43b945d/html5/thumbnails/35.jpg)
Fitting of Stochastic Process Rainfall statistics at different aggregation intervals:
i [mm/hr] 12 hours 1 hour n hours
Observed statistics Bound Constrained
Optimization (l-bfgs)
Limited memory Broyden, Fletcher, Goldfarb, Shanno
Calculated statistics
Set of model parameters NSRP
![Page 36: 9th International Conference URBAN DRAINAGE MODELLINGhikom.grf.bg.ac.rs/stari-sajt/9UDM/Presentations/238_PPT.pdf · Evaluating the impact of climate change on urban scale extreme](https://reader036.vdocument.in/reader036/viewer/2022081406/5f0daf8b7e708231d43b945d/html5/thumbnails/36.jpg)
Fitting of Stochastic Process
Calculated statistics
Bound Constrained Optimization (l-bfgs)
x P(Change Factor)
Set of model parameters NSRP
Home
Rainfall statistics at different aggregation intervals: i [mm/hr] 12 hours 1 hour n hours
Observed statistics
![Page 37: 9th International Conference URBAN DRAINAGE MODELLINGhikom.grf.bg.ac.rs/stari-sajt/9UDM/Presentations/238_PPT.pdf · Evaluating the impact of climate change on urban scale extreme](https://reader036.vdocument.in/reader036/viewer/2022081406/5f0daf8b7e708231d43b945d/html5/thumbnails/37.jpg)
Variability of precipitation results in the same GCM grid box.
GCM grid values for Kochi location Period: 2081-2100
GCM grid values for Kochi location Period: 2046-2065
24 h
ourly
rain
fall
mea
n [m
m]
Home
![Page 38: 9th International Conference URBAN DRAINAGE MODELLINGhikom.grf.bg.ac.rs/stari-sajt/9UDM/Presentations/238_PPT.pdf · Evaluating the impact of climate change on urban scale extreme](https://reader036.vdocument.in/reader036/viewer/2022081406/5f0daf8b7e708231d43b945d/html5/thumbnails/38.jpg)
Gas emissions scenarios Home