SlideShare a Scribd company logo
Addressing climate change uncertainty with a monte
carlo version of TIMES
Socrates Kypreos (PSI - Switzerland), Antti Lehtilä (VTT - Finland) & Evangelos Panos (PSI - Switzerland)
James Glynn (MaREI-UCC – Ireland)
RESEARCH FELLOW PhD, M.Sc. M.A.Econ, B.Eng
@james_glynn | james.glynn@ucc.ie
4th June 2018 | 38th International Energy Workshop | Paris FRANCE
Outline & Research Question
• Output from an collaborative IEA-ETSAP project lead by Paul Scherrer Institute to incorporate monte carlo
analysis (MCA) functionality in the TIMES source code.
• Methodological proof of concept
• Technical report will be available this week after the IEA-ETSAP ExCo meeting on Thursday & MCA in next source code update.
• https://iea-etsap.org/index.php/documentation
• Methodological paper & application with ETSAP-TIAM.v2.ucc.19.5 for; BASE, 2°C & 1°C OS scenarios in progress
• Today we show an example with ETSAP-TIAM for 1000 states of world for BASE, 2C and 1p5C targets.
• Typically energy system model (ESM) & IAM scenario ensembles are presented as statistically distributed
results; however IAM ensemble sets do not typically aim to fully explore and statistically represent the solution
space within their models reach.
• We have a lack of data availability for probability distribution of uncertain input data parameters
• How does the lack of presentation of uncertainty impact policy makers interpretations of single scenario results?
• How does this provide robust and resilient insights for policy design.
• Are we exploring the dominant uncertain variables in our typical scenario runs?
• Discount Rates? Climate Sensitivity?
• We develop a monte carlo analysis framework into the TIMES source-code & workflow: applying latin
hypercube sampling to uncertain probability distributions of inputs and develop a probabilistic method to
compute in parallel and post-process 1000’s of TIMES cases per scenario.
Outline & Research Question
• Output from an collaborative IEA-ETSAP project lead by Paul Scherrer Institute to incorporate monte carlo
analysis (MCA) functionality in the TIMES source code.
• Methodological proof of concept
• Technical report will be available this week after the IEA-ETSAP ExCo meeting on Thursday & MCA in next source code update.
• https://iea-etsap.org/index.php/documentation
• Methodological paper & application with ETSAP-TIAM.v2.ucc.19.5 for; BASE, 2°C & 1°C OS scenarios in progress
• Today we show an example with ETSAP-TIAM for 1000 states of world for BASE, 2C and 1p5C targets.
• Typically energy system model (ESM) & IAM scenario ensembles are presented as statistically distributed
results; however IAM ensemble sets do not typically aim to fully explore and statistically represent the solution
space within their models reach.
• We have a lack of data availability for probability distribution of uncertain input data parameters
• How does the lack of presentation of uncertainty impact policy makers interpretations of single scenario results?
• How does this provide robust and resilient insights for policy design.
• Are we exploring the dominant uncertain variables in our typical scenario runs?
• Discount Rates? Climate Sensitivity?
• We develop a monte carlo analysis framework into the TIMES source-code & workflow: applying latin
hypercube sampling to uncertain probability distributions of inputs and develop a probabilistic method to
compute in parallel and post-process 1000’s of TIMES cases per scenario.
Nice that Resources for
the Future (RFF) will
generate statistical
distribtuions of
uncertain model divers
Context: Uncertainty across social, economic & technical inputs
& the interdependant impact on the Paris goals solution space
Method: ETSAP-TIAM overview
• 15 Region linear programming bottom-up energy system
model of IEA-ETSAP
• Integrated model of the entire energy system
• Prospective analysis on medium to long term horizon (2100)
• Demand driven by exogenous energy service demands
• SSP2 from OECD Env-LINKS CGE model
• Regional Structural detail of the economy
• Partial and dynamic equilibrium
• Price-elastic demands
• General Equilibrium with MACRO
• Minimizes the total system cost
• Or Maximises Consumption/Utility
• Hybrid General Equilibrium MSA
• Optimal energy technology selection
• Environmental constraints
• GHG,
• Local Air Pollution & Health Damages
• Integrated Simple Climate Model (CMIP5 or FAIR calib)
• Nordhaus type climate damage function (NEW)
• Myopic and Stochastic, Robust Optimisation run options
• + Monte Carlo Assessment (new)
Climate
Module
Atm. Conc.
ΔForcing
ΔTemp
Used for
reporting &
setting
targets
Biomass
Potential
Renewable
Potential
Nuclear
Fossil Fuel
Reserves
(oil, coal, gas)
Extraction
Upstream
Fuels
Trade
Secondary
Transformation
OPEC/
NON-OPEC
regrouping
Electricity
Fuels
Electricity
Cogeneration
Heat
Hydrogen production
and distribution
End Use
Fuels
Industrial
Service
Composition
Auto Production
Cogeneration
Carbon
capture
CH4 options
Carbon
sequestration
Terrestrial
sequestration
Landfills Manure
Bio burning, rice,
enteric ferm
Wastewater
CH4 options
N2O options
CH4 options
OI****
GA****
CO****
Trade
ELC***
WIN SOL
GEO TDL
BIO***
NUC
HYD
BIO***
HETHET
ELC
ELC
SYNH2
BIO***
CO2
ELC
GAS***
COA***
Industrial
Tech.
Commercial
Tech.
Transport
Tech.
Residential
Tech.
Agriculture
Tech.
I***
I** (6)
T** (16)R** (11)C** (8)A** (1)
INDELC
INDELC
IS**
Demands
IND*** COM***AGR*** TRA***RES***
Non-energy
sectors (CH4)
OIL***
Moving forward from typical modelling 3-20 scenarios to
incorporating uncertainty into scenario modelling workflow
• Typically most TIMES scenario analysis only vary 1 or a few constraints and a
few parameter variants.
• The probability of an input our outcome is not known
• Interpreting deterministic scenarios makes decision making difficult
• How do we know if our scenario results and policy insights are robust & resilient?
• The cost of uncertainty is generally unknown and unquantified.
• We apply 3 scenarios in 1000 States of World (SOW) * 18 LHS perturbation
distributions for dominant variables in ETSAP-TIAM
• 3*1000*18 >> 54,000 uncertain variables
0.00
1.00
2.00
3.00
4.00
5.00
6.00
7.00
8.00
1950 2000 2050 2100 2150 2200 2250 2300 2350 2400 2450
Temperature(°C)
Typical 3 Scenario w Climate model temperature constraints
BASE_SSP2 2C_SSP2 1p5C_OS_SSP2
Moving forward from typical modelling 3-20 scenarios to
incorporating uncertainty into scenario modelling workflow
• Typically most TIMES scenario analysis only vary 1 or a few constraints and a
few parameter variants.
• The probability of an input our outcome is not known
• Interpreting deterministic scenarios makes decision making difficult
• How do we know if our scenario results and policy insights are robust & resilient?
• The cost of uncertainty is generally unknown and unquantified.
• We apply 3 scenarios in 1000 States of World (SOW) * 18 LHS perturbation
distributions for dominant variables in ETSAP-TIAM
• 3*1000*18 >> 54,000 uncertain variables
Correlated Input distributions drive statistical results outputs
Correlated Input distributions drive statistical results outputs
INPUTS OUTPUTS
Correlated Input distributions drive statistical results outputs
Note that these results show
the full distribution including
the tails before the 5%
interval and beyond 95%
INPUTS OUTPUTS
Scenarios + Latin Hyper Cube Sample of Uncertain Variables
• Typical BASE_SSP2 – “middle of the road” trends & rates of change and without climate or emission constraints
• 2C upper bound on DELTA_ATM from 2020-2400 with the coupled climate model (not carbon budget)
• 1.5C_OS_SSP2 upper bound from 2100-2400; (again not a cumulative carbon budget constraint)
Scenarios + Latin Hyper Cube Sample of Uncertain Variables
• Typical BASE_SSP2 – “middle of the road” trends & rates of change and without climate or emission constraints
• 2C upper bound on DELTA_ATM from 2020-2400 with the coupled climate model (not carbon budget)
• 1.5C_OS_SSP2 upper bound from 2100-2400; (again not a cumulative carbon budget constraint)
Results: Root summary scenario results (t°C, PES, CO2)
0
1
2
3
4
5
6
7
8
2005
2040
2080
2140
2220
2300
2380
DELTA-ATM(t°C)
BASE_SSP2
2C_SSP2
1p5C_OS_SSP2
0
200
400
600
800
1,000
1,200
1,400
2005
2030
2060
2090
2010
2040
2070
2100
2020
2050
2080
BASE_SSP2 2C_SSP2 1p5C_OS_SSP2
PrimaryEnergySupply(ExaJoules[EJ])
ALLBIOENERGY
Renewable except hydro and biomass
Hydro
Nuclear
Gas
Oil
-20
-10
0
10
20
30
40
50
60
70
2005
2030
2060
2090
2010
2040
2070
2100
2020
2050
2080
BASE_SSP2 2C_SSP2 1p5C_OS_SSP2
Fossil&IndustrialEmissions(GtCO2)
GBL
WEU
USA
SKO
ODA
MEX
MEA
JPN
IND
FSU
EEU
CSA
CHI
CAN
AUS
AFR
Results: Resilient Mitigation Technology options
Results: CO2 emissions as a f(temp (tC), CS, Price,DiscountRate)
Key take away messages
KM1
•The de-carbonization levels under SBA scenarios for the 2DS policy case have statistical mean carbon costs starting at $2005 196/tCO2
in 2030 to increase around $2005 380/tCO2 by the end of the century. This cost are high already by 2030, as the model prefers to start
with immediate policy actions to restrain temperature change. The shadow price of carbon control increases with the stringency of
temperature change and lowers with rapid technological change.
KM2
• Another interesting conclusion is the fact that countries with high economic development and stabilized energy consumption are selected
in the optimization process to reduce carbon emissions going to negative net emissions and give space to the former LDCs of high growth
to continue fuelling their economic development using fossil fuels in the end-use markets.…
KM3
• Most interesting is also the fact that the cost benefit index (CBI) is below one at almost zero discount rates which means the benefits of
climate change control that takes place in the last decades, compensates the energy systems cost increase of the 2DS scenario relative
to BASE that starts very early.… IMMEDIATE RAPID CLIMATE ACTION HAS A NET POSITIVE BENEFIT at ZERO DISCOUNT
RATES
KM4
•The “optimal” tax trajectory needed to avoid a high risk of overshooting critical temperature limits strongly depends on the applied PDF
describing the ECS parameter, the portfolio of technological options to be made available in the future, and the discount rate applied in
the utility objective function
KM5
•Feasibility of the policy target of 2 °C is achieved when immediate policy actions, successful penetration of advanced technologies to
decarbonize the energy and the transportation sectors, and net negative global emissions in the next decades are achieved too.
17
“Unlocking the potential of our
marine and renewable energy
resources through the power of
research and innovation”
18

More Related Content

PPTX
Icept 2017 abanihi
PPTX
NCPC IE Conference 2017 P Mukoma
PDF
Coal-Fired Boiler Fault Prediction using Artificial Neural Networks
PDF
65 sutterlueti using_advanced_pv_and_bo_s_modelling_and_algorithms_to_optimiz...
PDF
Proof energy@work midih oc2-demo_day
PPTX
NASPI_BPA_EPG_LSE_SDVCA_031816-F
PDF
62 friesen field_data_requirements_for_the_validation_of_pv_module_performanc...
PDF
Techniques for Running Large Numbers of Scenarios in TIMES
Icept 2017 abanihi
NCPC IE Conference 2017 P Mukoma
Coal-Fired Boiler Fault Prediction using Artificial Neural Networks
65 sutterlueti using_advanced_pv_and_bo_s_modelling_and_algorithms_to_optimiz...
Proof energy@work midih oc2-demo_day
NASPI_BPA_EPG_LSE_SDVCA_031816-F
62 friesen field_data_requirements_for_the_validation_of_pv_module_performanc...
Techniques for Running Large Numbers of Scenarios in TIMES

What's hot (20)

PDF
53 aron p_dobos_recent_and_planned_improvements_to_the_system_advisor_model_sam
PDF
Introducing external costs for Local Atmospheric Pollution in TIAM-MACRO to s...
PDF
Wind power forecasting: A Case Study in Terrain using Artificial Intelligence
PDF
63 matthiss comparison_of_pv_system_and_irradiation_models
PDF
IRJET-Load Forecasting using Fuzzy Logic
PDF
Uniform Load Distribution on Three Phase Alternater and Energy Audit of AMGOI
PDF
Wind power forecasting accuracy and uncertainty in Finland
PPTX
AQ modelling with real world traffic emissions
PDF
Sustainability for Multi-site
PPT
Thesis Presentation
PDF
How costs affect deployment of low carbon technologies - analysis with JRC-EU...
PPTX
Synthesis and Refinement of Artificial HVAC Sensor Data Intended for Supervis...
PDF
66 ueda system_performance_and_degradation_analysis_of_different_pv_technologies
PDF
IEA DSM Task 17: Integration of DSM, DG, RES and ES – Outcome of Phase 1 and ...
PDF
11 schroedter homscheidt_satellite_and_camera
PDF
Cuckoo search algorithm based for tunning both PI and FOPID controllers for ...
PDF
power Reliability pollard
53 aron p_dobos_recent_and_planned_improvements_to_the_system_advisor_model_sam
Introducing external costs for Local Atmospheric Pollution in TIAM-MACRO to s...
Wind power forecasting: A Case Study in Terrain using Artificial Intelligence
63 matthiss comparison_of_pv_system_and_irradiation_models
IRJET-Load Forecasting using Fuzzy Logic
Uniform Load Distribution on Three Phase Alternater and Energy Audit of AMGOI
Wind power forecasting accuracy and uncertainty in Finland
AQ modelling with real world traffic emissions
Sustainability for Multi-site
Thesis Presentation
How costs affect deployment of low carbon technologies - analysis with JRC-EU...
Synthesis and Refinement of Artificial HVAC Sensor Data Intended for Supervis...
66 ueda system_performance_and_degradation_analysis_of_different_pv_technologies
IEA DSM Task 17: Integration of DSM, DG, RES and ES – Outcome of Phase 1 and ...
11 schroedter homscheidt_satellite_and_camera
Cuckoo search algorithm based for tunning both PI and FOPID controllers for ...
power Reliability pollard
Ad

Similar to Addressing climate change uncertainty with a monte carloversion of TIMES (20)

PDF
Welcome to ETSAP Workshop 2024 in Singapore
PDF
Impact of integrating societal factors on the accuracy of optimization-based ...
PDF
Impacts of scenario definitions on CO2 mitigation cost in energy system models
PDF
Welcome to ETSAP Summer Workshop 2025, Prof. Brian Ó Gallachóir
PDF
Modelling global macroeconomic impacts of a carbon constrained energy system ...
PDF
Does myopic foresight modeling better capture the real-world energy transitio...
PDF
Analysis of the relative roles of supply-side and demand-side measures in tac...
PDF
IEA-GHG project and findings to date
PDF
Addressing Uncertainty in TIMES Using Monte Carlo Analysis
PDF
"Taking on TIAM" a new user´s experience and lessons learned
PDF
Using the TIMES-Ireland Model (TIM) to understand Ireland's Carbon Budget imp...
PDF
Global decarbonisation pathways: Contribution of different options in CO2 red...
PDF
Analysis of the relative roles of demand-side measures in tackling the global...
PDF
Abortion pills in Kuwait ௵ +918133066128**) Cytotec & mifepristone use, medic...
PDF
Pharmacy [][][][] DR.DAISY +27737758557 //Safe A'bo'rtion clinics// pi'l'ls i...
PDF
Abortion Pills In Bahrain?+2773-7758557#??##{{(+Abortion Pills For Sale In Jo...
PDF
ABORTION PILLS +27737758557- FOR SALE IN DUBAI MALL ABU DHABI DEIRA AJMAN SHA...
PDF
On Topic +27737758557*&*We have Abortion Pills / Cytotec Tablets Available in...
PDF
Terminating Pills +2773-7758557*@ Safe Abortion Pills For Sale In Mafikeng So...
PDF
DISCUSSION ON THE WORKPROGRAMME OF ANNEX XVI
Welcome to ETSAP Workshop 2024 in Singapore
Impact of integrating societal factors on the accuracy of optimization-based ...
Impacts of scenario definitions on CO2 mitigation cost in energy system models
Welcome to ETSAP Summer Workshop 2025, Prof. Brian Ó Gallachóir
Modelling global macroeconomic impacts of a carbon constrained energy system ...
Does myopic foresight modeling better capture the real-world energy transitio...
Analysis of the relative roles of supply-side and demand-side measures in tac...
IEA-GHG project and findings to date
Addressing Uncertainty in TIMES Using Monte Carlo Analysis
"Taking on TIAM" a new user´s experience and lessons learned
Using the TIMES-Ireland Model (TIM) to understand Ireland's Carbon Budget imp...
Global decarbonisation pathways: Contribution of different options in CO2 red...
Analysis of the relative roles of demand-side measures in tackling the global...
Abortion pills in Kuwait ௵ +918133066128**) Cytotec & mifepristone use, medic...
Pharmacy [][][][] DR.DAISY +27737758557 //Safe A'bo'rtion clinics// pi'l'ls i...
Abortion Pills In Bahrain?+2773-7758557#??##{{(+Abortion Pills For Sale In Jo...
ABORTION PILLS +27737758557- FOR SALE IN DUBAI MALL ABU DHABI DEIRA AJMAN SHA...
On Topic +27737758557*&*We have Abortion Pills / Cytotec Tablets Available in...
Terminating Pills +2773-7758557*@ Safe Abortion Pills For Sale In Mafikeng So...
DISCUSSION ON THE WORKPROGRAMME OF ANNEX XVI
Ad

More from IEA-ETSAP (20)

PDF
Generalized levelized cost as a metric for explaining model behavior of linea...
PDF
TIMES2JuMP project status report; Learnings on the feasibility of Migrating T...
PDF
A platform for open, realistic, and reproducible benchmarking of solvers on e...
PDF
Integrated Long-Term Planning and Short-Term Reliability Assessment for High-...
PDF
IEA H2 TCP Task 52 Hydrogen for Iron and Steel Making
PDF
TIMES-NZ 3.0: automating upstream data processing for an open-source workflow
PDF
Towards a national integrated energy, land and food system model for long ter...
PDF
Development of an AFOLU module for TIMES
PDF
The plant-level decarbonization pathways and mitigation cost of global oil re...
PDF
Near-optimal solutions for long-term energy planning facing the possible crit...
PDF
Integrated TIMES-E3ME-PLEXOS-DASMOD Modelling Framework for Assessing The Cze...
PDF
xl2times: progress update & a proof-of-concept interactive notebook-based wor...
PDF
Liberating energy models from modelers Amit Kanudia
PDF
The potential role of alternative fuels in the decarbonization of hard-to-aba...
PDF
Future Low-Carbon Hydrogen Production Technology Penetration with Aspen-Based...
PDF
Integrating Detailed Petrochemical Processes in Global Energy System Models f...
PDF
Are Heavy-Duty Vehicles at a Crossroads? A Real Options and Innovation Diffus...
PDF
An Assessment of the Impact of Electrification for Integration of Offshore Wi...
PDF
Role of Carbon Pricing and Emissions Constraint Pathways for India’s Net-Zero...
PDF
Implications of Party Politics for Future Energy System Scenarios in Norway
Generalized levelized cost as a metric for explaining model behavior of linea...
TIMES2JuMP project status report; Learnings on the feasibility of Migrating T...
A platform for open, realistic, and reproducible benchmarking of solvers on e...
Integrated Long-Term Planning and Short-Term Reliability Assessment for High-...
IEA H2 TCP Task 52 Hydrogen for Iron and Steel Making
TIMES-NZ 3.0: automating upstream data processing for an open-source workflow
Towards a national integrated energy, land and food system model for long ter...
Development of an AFOLU module for TIMES
The plant-level decarbonization pathways and mitigation cost of global oil re...
Near-optimal solutions for long-term energy planning facing the possible crit...
Integrated TIMES-E3ME-PLEXOS-DASMOD Modelling Framework for Assessing The Cze...
xl2times: progress update & a proof-of-concept interactive notebook-based wor...
Liberating energy models from modelers Amit Kanudia
The potential role of alternative fuels in the decarbonization of hard-to-aba...
Future Low-Carbon Hydrogen Production Technology Penetration with Aspen-Based...
Integrating Detailed Petrochemical Processes in Global Energy System Models f...
Are Heavy-Duty Vehicles at a Crossroads? A Real Options and Innovation Diffus...
An Assessment of the Impact of Electrification for Integration of Offshore Wi...
Role of Carbon Pricing and Emissions Constraint Pathways for India’s Net-Zero...
Implications of Party Politics for Future Energy System Scenarios in Norway

Recently uploaded (20)

PDF
Effects of rice-husk biochar and aluminum sulfate application on rice grain q...
PDF
FMM Slides For OSH Management Requirement
PPTX
Green Modern Sustainable Living Nature Presentation_20250226_230231_0000.pptx
PDF
Effect of salinity on biochimical and anatomical characteristics of sweet pep...
PDF
Global Natural Disasters in H1 2025 by Beinsure
PDF
Cave Diggers Simplified cave survey methods and mapping
PDF
Tree Biomechanics, a concise presentation
PPTX
Office Hours on Drivers of Tree Cover Loss
PPTX
Conformity-and-Deviance module 7 ucsp grade 12
PPTX
ser tico.pptxXYDTRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRY
PDF
Bai bao Minh chứng sk2-DBTrong-003757.pdf
DOCX
Epoxy Coated Steel Bolted Tanks for Fish Farm Water Provides Reliable Water f...
PDF
Effective factors on adoption of intercropping and it’s role on development o...
PDF
Urban Hub 50: Spirits of Place - & the Souls' of Places
PPTX
FIRE SAFETY SEMINAR SAMPLE FOR EVERYONE.pptx
PPTX
structure and components of Environment.pptx
DOCX
Epoxy Coated Steel Bolted Tanks for Anaerobic Digestion (AD) Plants Core Comp...
PDF
2-Reqerwsrhfdfsfgtdrttddjdiuiversion 2.pdf
DOCX
Double Membrane Roofs for Biogas Tanks Securely store produced biogas.docx
PDF
Earthquake, learn from the past and do it now.pdf
Effects of rice-husk biochar and aluminum sulfate application on rice grain q...
FMM Slides For OSH Management Requirement
Green Modern Sustainable Living Nature Presentation_20250226_230231_0000.pptx
Effect of salinity on biochimical and anatomical characteristics of sweet pep...
Global Natural Disasters in H1 2025 by Beinsure
Cave Diggers Simplified cave survey methods and mapping
Tree Biomechanics, a concise presentation
Office Hours on Drivers of Tree Cover Loss
Conformity-and-Deviance module 7 ucsp grade 12
ser tico.pptxXYDTRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRY
Bai bao Minh chứng sk2-DBTrong-003757.pdf
Epoxy Coated Steel Bolted Tanks for Fish Farm Water Provides Reliable Water f...
Effective factors on adoption of intercropping and it’s role on development o...
Urban Hub 50: Spirits of Place - & the Souls' of Places
FIRE SAFETY SEMINAR SAMPLE FOR EVERYONE.pptx
structure and components of Environment.pptx
Epoxy Coated Steel Bolted Tanks for Anaerobic Digestion (AD) Plants Core Comp...
2-Reqerwsrhfdfsfgtdrttddjdiuiversion 2.pdf
Double Membrane Roofs for Biogas Tanks Securely store produced biogas.docx
Earthquake, learn from the past and do it now.pdf

Addressing climate change uncertainty with a monte carloversion of TIMES

  • 1. Addressing climate change uncertainty with a monte carlo version of TIMES Socrates Kypreos (PSI - Switzerland), Antti Lehtilä (VTT - Finland) & Evangelos Panos (PSI - Switzerland) James Glynn (MaREI-UCC – Ireland) RESEARCH FELLOW PhD, M.Sc. M.A.Econ, B.Eng @james_glynn | james.glynn@ucc.ie 4th June 2018 | 38th International Energy Workshop | Paris FRANCE
  • 2. Outline & Research Question • Output from an collaborative IEA-ETSAP project lead by Paul Scherrer Institute to incorporate monte carlo analysis (MCA) functionality in the TIMES source code. • Methodological proof of concept • Technical report will be available this week after the IEA-ETSAP ExCo meeting on Thursday & MCA in next source code update. • https://iea-etsap.org/index.php/documentation • Methodological paper & application with ETSAP-TIAM.v2.ucc.19.5 for; BASE, 2°C & 1°C OS scenarios in progress • Today we show an example with ETSAP-TIAM for 1000 states of world for BASE, 2C and 1p5C targets. • Typically energy system model (ESM) & IAM scenario ensembles are presented as statistically distributed results; however IAM ensemble sets do not typically aim to fully explore and statistically represent the solution space within their models reach. • We have a lack of data availability for probability distribution of uncertain input data parameters • How does the lack of presentation of uncertainty impact policy makers interpretations of single scenario results? • How does this provide robust and resilient insights for policy design. • Are we exploring the dominant uncertain variables in our typical scenario runs? • Discount Rates? Climate Sensitivity? • We develop a monte carlo analysis framework into the TIMES source-code & workflow: applying latin hypercube sampling to uncertain probability distributions of inputs and develop a probabilistic method to compute in parallel and post-process 1000’s of TIMES cases per scenario.
  • 3. Outline & Research Question • Output from an collaborative IEA-ETSAP project lead by Paul Scherrer Institute to incorporate monte carlo analysis (MCA) functionality in the TIMES source code. • Methodological proof of concept • Technical report will be available this week after the IEA-ETSAP ExCo meeting on Thursday & MCA in next source code update. • https://iea-etsap.org/index.php/documentation • Methodological paper & application with ETSAP-TIAM.v2.ucc.19.5 for; BASE, 2°C & 1°C OS scenarios in progress • Today we show an example with ETSAP-TIAM for 1000 states of world for BASE, 2C and 1p5C targets. • Typically energy system model (ESM) & IAM scenario ensembles are presented as statistically distributed results; however IAM ensemble sets do not typically aim to fully explore and statistically represent the solution space within their models reach. • We have a lack of data availability for probability distribution of uncertain input data parameters • How does the lack of presentation of uncertainty impact policy makers interpretations of single scenario results? • How does this provide robust and resilient insights for policy design. • Are we exploring the dominant uncertain variables in our typical scenario runs? • Discount Rates? Climate Sensitivity? • We develop a monte carlo analysis framework into the TIMES source-code & workflow: applying latin hypercube sampling to uncertain probability distributions of inputs and develop a probabilistic method to compute in parallel and post-process 1000’s of TIMES cases per scenario. Nice that Resources for the Future (RFF) will generate statistical distribtuions of uncertain model divers
  • 4. Context: Uncertainty across social, economic & technical inputs & the interdependant impact on the Paris goals solution space
  • 5. Method: ETSAP-TIAM overview • 15 Region linear programming bottom-up energy system model of IEA-ETSAP • Integrated model of the entire energy system • Prospective analysis on medium to long term horizon (2100) • Demand driven by exogenous energy service demands • SSP2 from OECD Env-LINKS CGE model • Regional Structural detail of the economy • Partial and dynamic equilibrium • Price-elastic demands • General Equilibrium with MACRO • Minimizes the total system cost • Or Maximises Consumption/Utility • Hybrid General Equilibrium MSA • Optimal energy technology selection • Environmental constraints • GHG, • Local Air Pollution & Health Damages • Integrated Simple Climate Model (CMIP5 or FAIR calib) • Nordhaus type climate damage function (NEW) • Myopic and Stochastic, Robust Optimisation run options • + Monte Carlo Assessment (new) Climate Module Atm. Conc. ΔForcing ΔTemp Used for reporting & setting targets Biomass Potential Renewable Potential Nuclear Fossil Fuel Reserves (oil, coal, gas) Extraction Upstream Fuels Trade Secondary Transformation OPEC/ NON-OPEC regrouping Electricity Fuels Electricity Cogeneration Heat Hydrogen production and distribution End Use Fuels Industrial Service Composition Auto Production Cogeneration Carbon capture CH4 options Carbon sequestration Terrestrial sequestration Landfills Manure Bio burning, rice, enteric ferm Wastewater CH4 options N2O options CH4 options OI**** GA**** CO**** Trade ELC*** WIN SOL GEO TDL BIO*** NUC HYD BIO*** HETHET ELC ELC SYNH2 BIO*** CO2 ELC GAS*** COA*** Industrial Tech. Commercial Tech. Transport Tech. Residential Tech. Agriculture Tech. I*** I** (6) T** (16)R** (11)C** (8)A** (1) INDELC INDELC IS** Demands IND*** COM***AGR*** TRA***RES*** Non-energy sectors (CH4) OIL***
  • 6. Moving forward from typical modelling 3-20 scenarios to incorporating uncertainty into scenario modelling workflow • Typically most TIMES scenario analysis only vary 1 or a few constraints and a few parameter variants. • The probability of an input our outcome is not known • Interpreting deterministic scenarios makes decision making difficult • How do we know if our scenario results and policy insights are robust & resilient? • The cost of uncertainty is generally unknown and unquantified. • We apply 3 scenarios in 1000 States of World (SOW) * 18 LHS perturbation distributions for dominant variables in ETSAP-TIAM • 3*1000*18 >> 54,000 uncertain variables 0.00 1.00 2.00 3.00 4.00 5.00 6.00 7.00 8.00 1950 2000 2050 2100 2150 2200 2250 2300 2350 2400 2450 Temperature(°C) Typical 3 Scenario w Climate model temperature constraints BASE_SSP2 2C_SSP2 1p5C_OS_SSP2
  • 7. Moving forward from typical modelling 3-20 scenarios to incorporating uncertainty into scenario modelling workflow • Typically most TIMES scenario analysis only vary 1 or a few constraints and a few parameter variants. • The probability of an input our outcome is not known • Interpreting deterministic scenarios makes decision making difficult • How do we know if our scenario results and policy insights are robust & resilient? • The cost of uncertainty is generally unknown and unquantified. • We apply 3 scenarios in 1000 States of World (SOW) * 18 LHS perturbation distributions for dominant variables in ETSAP-TIAM • 3*1000*18 >> 54,000 uncertain variables
  • 8. Correlated Input distributions drive statistical results outputs
  • 9. Correlated Input distributions drive statistical results outputs INPUTS OUTPUTS
  • 10. Correlated Input distributions drive statistical results outputs Note that these results show the full distribution including the tails before the 5% interval and beyond 95% INPUTS OUTPUTS
  • 11. Scenarios + Latin Hyper Cube Sample of Uncertain Variables • Typical BASE_SSP2 – “middle of the road” trends & rates of change and without climate or emission constraints • 2C upper bound on DELTA_ATM from 2020-2400 with the coupled climate model (not carbon budget) • 1.5C_OS_SSP2 upper bound from 2100-2400; (again not a cumulative carbon budget constraint)
  • 12. Scenarios + Latin Hyper Cube Sample of Uncertain Variables • Typical BASE_SSP2 – “middle of the road” trends & rates of change and without climate or emission constraints • 2C upper bound on DELTA_ATM from 2020-2400 with the coupled climate model (not carbon budget) • 1.5C_OS_SSP2 upper bound from 2100-2400; (again not a cumulative carbon budget constraint)
  • 13. Results: Root summary scenario results (t°C, PES, CO2) 0 1 2 3 4 5 6 7 8 2005 2040 2080 2140 2220 2300 2380 DELTA-ATM(t°C) BASE_SSP2 2C_SSP2 1p5C_OS_SSP2 0 200 400 600 800 1,000 1,200 1,400 2005 2030 2060 2090 2010 2040 2070 2100 2020 2050 2080 BASE_SSP2 2C_SSP2 1p5C_OS_SSP2 PrimaryEnergySupply(ExaJoules[EJ]) ALLBIOENERGY Renewable except hydro and biomass Hydro Nuclear Gas Oil -20 -10 0 10 20 30 40 50 60 70 2005 2030 2060 2090 2010 2040 2070 2100 2020 2050 2080 BASE_SSP2 2C_SSP2 1p5C_OS_SSP2 Fossil&IndustrialEmissions(GtCO2) GBL WEU USA SKO ODA MEX MEA JPN IND FSU EEU CSA CHI CAN AUS AFR
  • 14. Results: Resilient Mitigation Technology options
  • 15. Results: CO2 emissions as a f(temp (tC), CS, Price,DiscountRate)
  • 16. Key take away messages KM1 •The de-carbonization levels under SBA scenarios for the 2DS policy case have statistical mean carbon costs starting at $2005 196/tCO2 in 2030 to increase around $2005 380/tCO2 by the end of the century. This cost are high already by 2030, as the model prefers to start with immediate policy actions to restrain temperature change. The shadow price of carbon control increases with the stringency of temperature change and lowers with rapid technological change. KM2 • Another interesting conclusion is the fact that countries with high economic development and stabilized energy consumption are selected in the optimization process to reduce carbon emissions going to negative net emissions and give space to the former LDCs of high growth to continue fuelling their economic development using fossil fuels in the end-use markets.… KM3 • Most interesting is also the fact that the cost benefit index (CBI) is below one at almost zero discount rates which means the benefits of climate change control that takes place in the last decades, compensates the energy systems cost increase of the 2DS scenario relative to BASE that starts very early.… IMMEDIATE RAPID CLIMATE ACTION HAS A NET POSITIVE BENEFIT at ZERO DISCOUNT RATES KM4 •The “optimal” tax trajectory needed to avoid a high risk of overshooting critical temperature limits strongly depends on the applied PDF describing the ECS parameter, the portfolio of technological options to be made available in the future, and the discount rate applied in the utility objective function KM5 •Feasibility of the policy target of 2 °C is achieved when immediate policy actions, successful penetration of advanced technologies to decarbonize the energy and the transportation sectors, and net negative global emissions in the next decades are achieved too.
  • 17. 17
  • 18. “Unlocking the potential of our marine and renewable energy resources through the power of research and innovation” 18