dynamic systems sea-level rise model: …mesl.ce.gatech.edu/research/forcing function...
TRANSCRIPT
![Page 1: Dynamic Systems Sea-Level Rise Model: …mesl.ce.gatech.edu/RESEARCH/Forcing Function Analysis_2012.pdfSea-Level Rise Model: Interpretation of the Forcing Function ... (SPM) released](https://reader031.vdocument.in/reader031/viewer/2022022508/5acfa9867f8b9a1d328d52e4/html5/thumbnails/1.jpg)
Dynamic Systems Sea-Level Rise Model: Interpretation of the Forcing Function
Jiabao Guan, Biao Chang and Mustafa M. Aral Multimedia Environmental Simulations Laboratory School of Civil and Environmental Engineering Georgia Institute of Technology
![Page 2: Dynamic Systems Sea-Level Rise Model: …mesl.ce.gatech.edu/RESEARCH/Forcing Function Analysis_2012.pdfSea-Level Rise Model: Interpretation of the Forcing Function ... (SPM) released](https://reader031.vdocument.in/reader031/viewer/2022022508/5acfa9867f8b9a1d328d52e4/html5/thumbnails/2.jpg)
Climate Change / Sea-level Rise DATA on: Temperature Sea-Level Rise Radiative forces (CO2 Emissions)
![Page 3: Dynamic Systems Sea-Level Rise Model: …mesl.ce.gatech.edu/RESEARCH/Forcing Function Analysis_2012.pdfSea-Level Rise Model: Interpretation of the Forcing Function ... (SPM) released](https://reader031.vdocument.in/reader031/viewer/2022022508/5acfa9867f8b9a1d328d52e4/html5/thumbnails/3.jpg)
Global Temperature Change:
![Page 4: Dynamic Systems Sea-Level Rise Model: …mesl.ce.gatech.edu/RESEARCH/Forcing Function Analysis_2012.pdfSea-Level Rise Model: Interpretation of the Forcing Function ... (SPM) released](https://reader031.vdocument.in/reader031/viewer/2022022508/5acfa9867f8b9a1d328d52e4/html5/thumbnails/4.jpg)
Global Sea-Level Rise:
![Page 5: Dynamic Systems Sea-Level Rise Model: …mesl.ce.gatech.edu/RESEARCH/Forcing Function Analysis_2012.pdfSea-Level Rise Model: Interpretation of the Forcing Function ... (SPM) released](https://reader031.vdocument.in/reader031/viewer/2022022508/5acfa9867f8b9a1d328d52e4/html5/thumbnails/5.jpg)
Data on CO2 Emissions:
![Page 6: Dynamic Systems Sea-Level Rise Model: …mesl.ce.gatech.edu/RESEARCH/Forcing Function Analysis_2012.pdfSea-Level Rise Model: Interpretation of the Forcing Function ... (SPM) released](https://reader031.vdocument.in/reader031/viewer/2022022508/5acfa9867f8b9a1d328d52e4/html5/thumbnails/6.jpg)
Sea Level Rise (m at 2090-2099 relative to 1980-1999)
Case Model-based range excluding future rapid dynamical changes in ice flow
B1 scenario 0.18 – 0.38 A1T scenario 0.20 – 0.45 B2 scenario 0.20 – 0.43 A1B scenario 0.21 – 0.48 A2 scenario 0.23 – 0.51 A1FI scenario 0.26 – 0.59
The Summary for Policy Makers (SPM) released recently provide the following table of sea level rise projections (IPCC4th Framework Report, 2007):
The IPCC Study:
![Page 7: Dynamic Systems Sea-Level Rise Model: …mesl.ce.gatech.edu/RESEARCH/Forcing Function Analysis_2012.pdfSea-Level Rise Model: Interpretation of the Forcing Function ... (SPM) released](https://reader031.vdocument.in/reader031/viewer/2022022508/5acfa9867f8b9a1d328d52e4/html5/thumbnails/7.jpg)
IPCC estimates:
![Page 8: Dynamic Systems Sea-Level Rise Model: …mesl.ce.gatech.edu/RESEARCH/Forcing Function Analysis_2012.pdfSea-Level Rise Model: Interpretation of the Forcing Function ... (SPM) released](https://reader031.vdocument.in/reader031/viewer/2022022508/5acfa9867f8b9a1d328d52e4/html5/thumbnails/8.jpg)
Semi-Empirical Models: Rahmstorf’s Study (Science, Vol. 315 pp.19, 2007) and others:
![Page 9: Dynamic Systems Sea-Level Rise Model: …mesl.ce.gatech.edu/RESEARCH/Forcing Function Analysis_2012.pdfSea-Level Rise Model: Interpretation of the Forcing Function ... (SPM) released](https://reader031.vdocument.in/reader031/viewer/2022022508/5acfa9867f8b9a1d328d52e4/html5/thumbnails/9.jpg)
{ }0.5 1.4m−
Semi-Empirical Models: Rahmstorf’s Study (Science, Vol. 315 pp.19, 2007) and others.
Sea-level Rise Interval Predicted above 1990 level:
IPCC Interval Predicted above 1990 level:
{ }0.09 0.88m−
![Page 10: Dynamic Systems Sea-Level Rise Model: …mesl.ce.gatech.edu/RESEARCH/Forcing Function Analysis_2012.pdfSea-Level Rise Model: Interpretation of the Forcing Function ... (SPM) released](https://reader031.vdocument.in/reader031/viewer/2022022508/5acfa9867f8b9a1d328d52e4/html5/thumbnails/10.jpg)
Major limitations of the previous empirical models
No feedback of SLR on temperature change
Zero-dimensional, thus does not capture spatial variations in SLR
Impact of external forcing is not considered
![Page 11: Dynamic Systems Sea-Level Rise Model: …mesl.ce.gatech.edu/RESEARCH/Forcing Function Analysis_2012.pdfSea-Level Rise Model: Interpretation of the Forcing Function ... (SPM) released](https://reader031.vdocument.in/reader031/viewer/2022022508/5acfa9867f8b9a1d328d52e4/html5/thumbnails/11.jpg)
A Dynamic Systems Model
![Page 12: Dynamic Systems Sea-Level Rise Model: …mesl.ce.gatech.edu/RESEARCH/Forcing Function Analysis_2012.pdfSea-Level Rise Model: Interpretation of the Forcing Function ... (SPM) released](https://reader031.vdocument.in/reader031/viewer/2022022508/5acfa9867f8b9a1d328d52e4/html5/thumbnails/12.jpg)
Earlier studies showed that the relationship between T & H is linear.
Our Hypothesis: (Dynamic Systems Model):
Temperature:
Sea-Level:
( )1 1, , ,T f T H U c=
( )2 2, , ,H f T H U c=
Expected Relationship:
![Page 13: Dynamic Systems Sea-Level Rise Model: …mesl.ce.gatech.edu/RESEARCH/Forcing Function Analysis_2012.pdfSea-Level Rise Model: Interpretation of the Forcing Function ... (SPM) released](https://reader031.vdocument.in/reader031/viewer/2022022508/5acfa9867f8b9a1d328d52e4/html5/thumbnails/13.jpg)
11 12 1
21 22 2
( ) ( ) ( )
( ) ( ) ( )
dT t a T t a H t cdt
dH t a T t a H t cdt
= + +
= + +
Initial Model Proposed:
(Simplified)
![Page 14: Dynamic Systems Sea-Level Rise Model: …mesl.ce.gatech.edu/RESEARCH/Forcing Function Analysis_2012.pdfSea-Level Rise Model: Interpretation of the Forcing Function ... (SPM) released](https://reader031.vdocument.in/reader031/viewer/2022022508/5acfa9867f8b9a1d328d52e4/html5/thumbnails/14.jpg)
( ) ( ) ( ),d t dT t dH tdt dt dt
τ =
X
( ) ( )d t tdt
= +X AX C
11 12
21 22
a aa a
=
A 1
2
cc
=
C
Initial Model Proposed:
(Simplified)
![Page 15: Dynamic Systems Sea-Level Rise Model: …mesl.ce.gatech.edu/RESEARCH/Forcing Function Analysis_2012.pdfSea-Level Rise Model: Interpretation of the Forcing Function ... (SPM) released](https://reader031.vdocument.in/reader031/viewer/2022022508/5acfa9867f8b9a1d328d52e4/html5/thumbnails/15.jpg)
Initial Model Proposed:
(Simplified)
![Page 16: Dynamic Systems Sea-Level Rise Model: …mesl.ce.gatech.edu/RESEARCH/Forcing Function Analysis_2012.pdfSea-Level Rise Model: Interpretation of the Forcing Function ... (SPM) released](https://reader031.vdocument.in/reader031/viewer/2022022508/5acfa9867f8b9a1d328d52e4/html5/thumbnails/16.jpg)
Dynamic Systems Model:
( )tU
( )tXddtX
( )tU
( )tX ddtX
System
![Page 17: Dynamic Systems Sea-Level Rise Model: …mesl.ce.gatech.edu/RESEARCH/Forcing Function Analysis_2012.pdfSea-Level Rise Model: Interpretation of the Forcing Function ... (SPM) released](https://reader031.vdocument.in/reader031/viewer/2022022508/5acfa9867f8b9a1d328d52e4/html5/thumbnails/17.jpg)
( )
( )
11 12 13 1
21 22 23 2
( ) ( ) ( )
( ) ( ) ( )
i ii
i ii
dT t a T t a H t a U t cdt
dH t a T t a H t a U t cdt
= + + +
= + + +
∑
∑
Proposed Model:
![Page 18: Dynamic Systems Sea-Level Rise Model: …mesl.ce.gatech.edu/RESEARCH/Forcing Function Analysis_2012.pdfSea-Level Rise Model: Interpretation of the Forcing Function ... (SPM) released](https://reader031.vdocument.in/reader031/viewer/2022022508/5acfa9867f8b9a1d328d52e4/html5/thumbnails/18.jpg)
more flexible; may answer more questions; may provide control analysis perspective; and, hopefully will be more useful. potential drawback may require more data.
The Proposed model is:
![Page 19: Dynamic Systems Sea-Level Rise Model: …mesl.ce.gatech.edu/RESEARCH/Forcing Function Analysis_2012.pdfSea-Level Rise Model: Interpretation of the Forcing Function ... (SPM) released](https://reader031.vdocument.in/reader031/viewer/2022022508/5acfa9867f8b9a1d328d52e4/html5/thumbnails/19.jpg)
( ) ( ) ( ) ( )d t t t tdt
= + + +X AX BU C w
11 12
21 22
a aa a
=
A 11 12 1
21 22 2
m
m
b b bb b b
=
B
1
2
cc
=
C
Proposed Model:
![Page 20: Dynamic Systems Sea-Level Rise Model: …mesl.ce.gatech.edu/RESEARCH/Forcing Function Analysis_2012.pdfSea-Level Rise Model: Interpretation of the Forcing Function ... (SPM) released](https://reader031.vdocument.in/reader031/viewer/2022022508/5acfa9867f8b9a1d328d52e4/html5/thumbnails/20.jpg)
( 1) ( ) ( ) ( )k k k k+ = + + + 1XΦX ΓU Ω w
; ;t t t= + ∆ = ∆ = ∆Φ I A Γ B Ω C
Discrete form of the Proposed Model:
![Page 21: Dynamic Systems Sea-Level Rise Model: …mesl.ce.gatech.edu/RESEARCH/Forcing Function Analysis_2012.pdfSea-Level Rise Model: Interpretation of the Forcing Function ... (SPM) released](https://reader031.vdocument.in/reader031/viewer/2022022508/5acfa9867f8b9a1d328d52e4/html5/thumbnails/21.jpg)
( ) ( ){ }* minimizei
i i i iF τ= − −φ
YΛφ Y Λφ
LSM Model to determine aij , bij & ci :
![Page 22: Dynamic Systems Sea-Level Rise Model: …mesl.ce.gatech.edu/RESEARCH/Forcing Function Analysis_2012.pdfSea-Level Rise Model: Interpretation of the Forcing Function ... (SPM) released](https://reader031.vdocument.in/reader031/viewer/2022022508/5acfa9867f8b9a1d328d52e4/html5/thumbnails/22.jpg)
Confidence Interval:
( )( )( )( )
12/2, 4
12/2, 4
ˆ ˆ ˆ ˆˆ( ) ( ) 1 ( ) ( )
ˆ ˆ ˆ ˆˆ( ) ( ) 1 ( ) ( )
CI n T
CI n H
T k T k t Z k Z k
H k H k t Z k Z k
τ τα
τ τα
σ
σ
−
−
−
−
= ± +
= ±
+
Λ Λ
Λ Λ
(Wadsworth, 1998)
![Page 23: Dynamic Systems Sea-Level Rise Model: …mesl.ce.gatech.edu/RESEARCH/Forcing Function Analysis_2012.pdfSea-Level Rise Model: Interpretation of the Forcing Function ... (SPM) released](https://reader031.vdocument.in/reader031/viewer/2022022508/5acfa9867f8b9a1d328d52e4/html5/thumbnails/23.jpg)
Application:
2-year moving average outcome is used for
both state variables.
![Page 24: Dynamic Systems Sea-Level Rise Model: …mesl.ce.gatech.edu/RESEARCH/Forcing Function Analysis_2012.pdfSea-Level Rise Model: Interpretation of the Forcing Function ... (SPM) released](https://reader031.vdocument.in/reader031/viewer/2022022508/5acfa9867f8b9a1d328d52e4/html5/thumbnails/24.jpg)
Temperature Data:
![Page 25: Dynamic Systems Sea-Level Rise Model: …mesl.ce.gatech.edu/RESEARCH/Forcing Function Analysis_2012.pdfSea-Level Rise Model: Interpretation of the Forcing Function ... (SPM) released](https://reader031.vdocument.in/reader031/viewer/2022022508/5acfa9867f8b9a1d328d52e4/html5/thumbnails/25.jpg)
Sea-Level Data:
![Page 26: Dynamic Systems Sea-Level Rise Model: …mesl.ce.gatech.edu/RESEARCH/Forcing Function Analysis_2012.pdfSea-Level Rise Model: Interpretation of the Forcing Function ... (SPM) released](https://reader031.vdocument.in/reader031/viewer/2022022508/5acfa9867f8b9a1d328d52e4/html5/thumbnails/26.jpg)
Resulting Matrix Coefficients:
System System matrix Control matrix
Constant vector
Discrete
Continuous
0.75042 -0.000530.41015 0.99568
0.002450
0.857440.25863−
-0.24958 -0.000530.41015 -0.00432
0.002450
0.857440.25863−
![Page 27: Dynamic Systems Sea-Level Rise Model: …mesl.ce.gatech.edu/RESEARCH/Forcing Function Analysis_2012.pdfSea-Level Rise Model: Interpretation of the Forcing Function ... (SPM) released](https://reader031.vdocument.in/reader031/viewer/2022022508/5acfa9867f8b9a1d328d52e4/html5/thumbnails/27.jpg)
-1
-0.8
-0.6
-0.4
-0.2
0
0.2
0.4
Tem
pera
ture
(o C)
ObservationReconstructionConfidence interval
1880 1900 1920 1940 1960 1980 2000Year
-20
-15
-10
-5
0
5
10
Sea
-leve
l ris
e (c
m)
ObservationReconstructionConfidence interval
(a)
(b)
Global CO2 Impact: Calibration
R2 = 0.80
R2 = 0.97
![Page 28: Dynamic Systems Sea-Level Rise Model: …mesl.ce.gatech.edu/RESEARCH/Forcing Function Analysis_2012.pdfSea-Level Rise Model: Interpretation of the Forcing Function ... (SPM) released](https://reader031.vdocument.in/reader031/viewer/2022022508/5acfa9867f8b9a1d328d52e4/html5/thumbnails/28.jpg)
-1
-0.5
0
0.5
1
1.5
2
2.5Te
mpe
ratu
re (o C
)
ObservationProjectionConfidence interval
1880 1900 1920 1940 1960 1980 2000 2020 2040 2060 2080 2100Year
-20
-10
0
10
20
30
40
50
60
70
Sea
-leve
l ris
e (c
m)
ObservationProjectionConfidence intervalRahmstorf's results
(a)
(b)
Global CO2 Impact: Prediction for “2 oC” Scenario (Scenario is designed to limit global warming in 2100 2oC above the temperature in 2000 )
(Hansen et al., 2000)
![Page 29: Dynamic Systems Sea-Level Rise Model: …mesl.ce.gatech.edu/RESEARCH/Forcing Function Analysis_2012.pdfSea-Level Rise Model: Interpretation of the Forcing Function ... (SPM) released](https://reader031.vdocument.in/reader031/viewer/2022022508/5acfa9867f8b9a1d328d52e4/html5/thumbnails/29.jpg)
1880 1900 1920 1940 1960 1980 2000 2020 2040 2060 2080 2100Year
0
200
400
600
800
1000
1200
CO
2 e
mis
sio
n (
pp
m)
SCENARIOSA1FIA1BA1TA2B1B2
IPCC Global CO2 Emission Scenarios:
![Page 30: Dynamic Systems Sea-Level Rise Model: …mesl.ce.gatech.edu/RESEARCH/Forcing Function Analysis_2012.pdfSea-Level Rise Model: Interpretation of the Forcing Function ... (SPM) released](https://reader031.vdocument.in/reader031/viewer/2022022508/5acfa9867f8b9a1d328d52e4/html5/thumbnails/30.jpg)
-1
0
1
2
3
4
5
6
Tem
pera
ture
(o C
)
ObservationA1FIA1BA1TA2B1B2
1880 1900 1920 1940 1960 1980 2000 2020 2040 2060 2080 2100Year
-20
0
20
40
60
80
100
120
Sea
-leve
l ris
e (c
m)
ObservationA1FIA1BA1TA2B1B2
(a)
(b)
Global CO2 Impact: Comparison of results
Solid lines are Dyn. Sys. Model results. Dashed lines are IPCC results.
![Page 31: Dynamic Systems Sea-Level Rise Model: …mesl.ce.gatech.edu/RESEARCH/Forcing Function Analysis_2012.pdfSea-Level Rise Model: Interpretation of the Forcing Function ... (SPM) released](https://reader031.vdocument.in/reader031/viewer/2022022508/5acfa9867f8b9a1d328d52e4/html5/thumbnails/31.jpg)
For example: If we want to restrict the global temperature and sea-level rise to zero-growth, the global CO2 emission should be controlled with the relationship given by
Decision Making on Global CO2 Emissions:
2( ) 101.87 ( ) 0.22 ( ) 349.98COu t T t H t= + +
2( )COu t represents the amount of yearly CO2 emission (ppm).
where