the key to the future lies in the past.ntthung/wp-content/uploads/...Β Β· 2018. 1. 17.Β Β· linear...
TRANSCRIPT
18 January 2018
Stefano Galelli
people.sutd.edu.sg/~stefano_galelli/
Resilient Water Systems Group
REVEALS HISTORY OF REGIME SHIFTS
STREAMFLOW RECONSTRUCTION
IN NORTHERN THAILAND
Nguyen Tan Thai Hung
people.sutd.edu.sg/~ntthung/
A LINEAR DYNAMICAL SYSTEMS APPROACH
TO
The key to the future lies in the past.
3
Paleohydrology
4
ΠαλαιΟΟ = old, ancient
Paleohydrology
Proxy data
β’ Tree rings
β’ Ice core
β’ Corals
β’ β¦
Instrumental data
β’ Streamflow
β’ Precipitation
β’ Drought index
β’ β¦
Model
Paleoreconstructed
data
ΠαλαιΟΟ = old, ancient
Study site: Ping River
Monsoon Asia Drought Atlas (MADA)
Cook, E. R., Anchukaitis, K. J., Buckley, B. M., DβArrigo, R. D., Jacoby, G. C., & Wright, W. E. (2010). Asian Monsoon Failure and Megadrought
During the Last Millennium. Science, 328(5977), 486β489. http://doi.org/10.1126/science.1185188
Figure 1B in Cook et al (2010)
Temporal resolution Annual
Spatial resolution 2.5o x 2.5o
Temporal range 1300 β 2005
Gridded time series of the Palmerβs Drought Severity Index
The conventional method
β’ How do we model catchment dynamics?
β’ Will a dynamic model be more accurate?
β’ What more insights can we gain with a
dynamic model?
8
π¦π‘ = πΌ + π½π’π‘ + ππ‘
Linear dynamical systems
π₯π‘+1 = π΄π₯π‘ + π΅π’π‘ +π€π‘
π¦π‘ = πΆπ₯π‘ + π·π’π‘ + π£π‘
π€π‘ βΌ π© 0,ππ£π‘ βΌ π©(0,π )π₯1 βΌ π©(π1, π1)
Systemπ₯
Inputπ’
Outputπ¦
π₯ β βπ system state
π¦ β βπ system output
π’ β βπ system input
π΄ β βπΓπ state transition matrix
π΅ β βπΓπ input-state matrix
πΆ β βπΓπ observation matrix
π· β βπΓπ input-observation matrix
π β βπΓπ covariance matrix of the state noise
π β βπΓπ covariance matrix of the observation noise
Learning: Expectation-Maximization
Shumway, R. H., & Stoffer, D. S. (1982). An Approach to The Time Series Smoothing and Forecasting Using the EM Algorithm. Journal of Time Series Analysis, 3(4), 253β264. https://doi.org/10.1111/j.1467-9892.1982.tb00349.x
Ghahramani, Z., & Hinton, G. E. (1996). Parameter Estimation for Linear Dynamical Systems. Technical Report CRG-TR-96-2. https://doi.org/10.1080/00207177208932224
Cheng, S., & Sabes, P. N. (2006). Modeling Sensorimotor Learning with Linear Dynamical Systems. Neural Computation, 18(4), 760β793. https://doi.org/10.1162/089976606775774651
E-Step
αππ+1 = arg max β π| π, αππ
M-Step
π αππ = πΌ π|π, αππ
Forward pass:
Kalman filter
Backward pass:
RTS recursion
Maximum
likelihood
ΰ·π₯π‘|π‘ = πΌ π₯π‘|π¦1, β¦ , π¦π‘, αππ
ΰ·π₯π‘|π = πΌ π₯π‘|π¦1, β¦ , π¦π , αππ
Kalman filter
Forward pass:
Kalman filter
Backward pass:
RTS recursion
Maximum
likelihood
ΰ·π₯π‘|π‘ = πΌ π₯π‘|π¦1, β¦ , π¦π‘, αππ
Kalman, R. E. (1960). A New Approach to Linear Filtering and Prediction Problems. Journal of Basic Engineering, 82(1), 35. https://doi.org/10.1115/1.3662552
Faragher, R. (2012). Understanding the basis of the Kalman filter via a simple and intuitive derivation [lecture notes]. IEEE Signal Processing Magazine, 29(5), 128β132. https://doi.org/10.1109/MSP.2012.2203621
Figure 5 in Faragher (2012)
ΰ·π₯π‘|π‘β1 = π΄ΰ·π₯π‘β1|π‘β1 + π΅π’π‘ΰ·π¦π‘|π‘β1 = πΆ ΰ·π₯π‘|π‘β1 +π·π’π‘ππ‘|π‘β1 = π΄ ππ‘β1|π‘β1π΄β² + π
πΎπ‘ = ππ‘|π‘β1πΆβ² πΆ ππ‘|π‘β1πΆ
β² + π β1
ΰ·π₯π‘|π‘ = ΰ·π₯π‘|π‘β1 + πΎπ‘ π¦π‘ β ΰ·π¦π‘|π‘β1ππ‘|π‘ = πΌ β πΎπ‘πΆ ππ‘|π‘β1
For π‘ = 2,β¦ , π
RTS recursion
Forward pass:
Kalman filter
Backward pass:
RTS recursion
Maximum
likelihood
ΰ·π₯π‘|π = πΌ π₯π‘|π¦1, β¦ , π¦π , αππ
π½π‘ = ππ‘|π‘π΄ ππ‘+1|π‘β1
ΰ·π₯π‘|π = ΰ·π₯π‘|π‘ + π½π‘ ΰ·π₯π‘+1|π β ΰ·π₯π‘+1|π‘ππ‘|π = ππ‘|π‘ + π½π‘ ππ‘+1|π β ππ‘+1|π‘ π½π‘
β²
ΰ·π¦π‘|π = πΆ ΰ·π₯π‘|π + π·π’π‘
Rauch, H. E., Tung, F., & Striebel, C. T. (1965). Maximum likelihood estimates of linear dynamic systems. AIAA Journal, 3(8), 1445β1450. https://doi.org/10.2514/3.3166
For π‘ = π β 1,β¦ , 1
Maximum likelihood estimation
Forward pass:
Kalman filter
Backward pass:
RTS recursion
Maximum
likelihood
Quadratic terms only Analytical solutions
Algorithm 1: LDS-EM
14
π‘ = π,β¦ , 1
Simultaneous learningβreconstruction
15
π¦π‘ β ΰ·π¦π‘|π
π¦π‘ β ΰ·π¦π‘|π‘β1
Replace missing π¦π‘ with its best available estimate
Forward pass
M-step
π₯1
Rationale for SLR
16
The substitution turns all terms related to missing π¦π‘ into zero
ΰ·π₯π‘|π‘ = ΰ·π₯π‘|π‘β1 + πΎπ‘ π¦π‘ β ΰ·π¦π‘|π‘β1
E-Step
M-Step
Algorithm 2: SLR
17
π‘ = π, β¦ , 1
Model performance
18
π πΈ = 1 β
π‘βπ±π¦π‘ β ΰ·π¦π‘
2
Οπ‘βπ± π¦π‘ β π¦π
2
πΆπΈ = 1 βΟπ‘βπ± π¦π‘ β ΰ·π¦π‘
2
Οπ‘βπ± π¦π‘ β π¦π£
2
Model performance
19
Residual analysis
20
A reconstructed history of the Ping
21
Figure 2 in Cook et al (2010)
Stochastic replicates
22
Conclusions
β’ Replacement for conventional method
Better model performance and desirable features
β’ A more conservative policy for the Bhumibol
There seems to be less water in the system
β’ Regional hydrological understanding (complementing the MADA)
History of regime shifts
β’ Direct application: regime-informed reservoir operation
Stochastic replicates of both streamflow and regime
23
APPENDICES
DendrochronologyΞ΄ΞΞ½Ξ΄ΟΞΏΞ½ (tree limb) + ΟΟΟΞ½ΞΏΟ (time) = tree dating
Other reconstructions
Woodhouse et al,
2006
Gangopadhyay et al,
2009Devineni et al, 2013 Patskoski et al, 2015 Ho et al, 2016
Lo
cati
on
&D
ata
Colorado River
β’ 4 stations
β’ 62 chronologies
Colorado River at Lees
Ferry, Arizona
β’ 62 chronologies
Upper Delaware River
Basin
β’ 5 stations
β’ 8 chronologies
South-eastern US (NC,
SC, GA, FL)
β’ 8 stations
β’ 7 chronologies
Missouri River Basin
β’ 55 stations
β’ LBDA
Pe
rfo
rman
ce
β’ RE ~ 0.65 - 0.8
β’ nRMSE ~ 0.14
β’ adjusted R2
~ 0.7 -
0.8
β’ R2
= 0.76β’ RE ~ 0.2 - 0.5
β’ CE ~ 0.1 - 0.5
β’ Adjusted R2 ~ 0.1 -
0.4
β’ Normalized RMSE
~0.25 - 0.5
β’ NSE (positive /
negative, average
positive)
β’ Reduction of error
(mostly positive,
average around ~0.2
β’ Adjusted R2
~0.5 -
0.9
Search radius
27
Spatial correlation
Site Distance r p-value
LS001 406.71438 -0.2225 0.0407
LS002 438.74650 -0.1447 0.1863
TH001 55.37224 0.2024 0.0632
TH002 354.28653 0.1293 0.2757
TH003 369.80428 -0.0365 0.7589
TH004 423.49371 0.1829 0.0919
TH006 85.10499 -0.0358 0.7464
MADA Tree rings
M-step solution
29
Wavelet analysis
30
Applications
β’ Drought adaptation planning
β Agriculture & Agri-Food Canada
β Prairie Provinces Water Board
β Denver Water Board
β’ Informing the public (Colorado River)
β’ Reliability of urban water supply
β Cities of Calgary and Edmonton
31