el ni ño and how to get rid of it
DESCRIPTION
El Ni ño and how to get rid of it. With thanks to Prashant SardeshmukhLudmila Matrosova Ping ChangMoritz Fl ügel Brian EwaldRoger Temam The Climate Diagnostics CenterOGP. C écile Penland. Review of Linear Inverse Modeling. Assume linear dynamics: d x /dt = B x + x - PowerPoint PPT PresentationTRANSCRIPT
![Page 1: El Ni ño and how to get rid of it](https://reader035.vdocument.in/reader035/viewer/2022081519/56812d15550346895d91fe55/html5/thumbnails/1.jpg)
El Niño and how to get rid of it
Cécile Penland
With thanks to
Prashant Sardeshmukh Ludmila Matrosova
Ping Chang Moritz Flügel
Brian Ewald Roger Temam
The Climate Diagnostics Center OGP
![Page 2: El Ni ño and how to get rid of it](https://reader035.vdocument.in/reader035/viewer/2022081519/56812d15550346895d91fe55/html5/thumbnails/2.jpg)
Review of Linear Inverse Modeling
Assume linear dynamics:
dx/dt = Bx +
Diagnose Green function from data:
G() = exp(B) = <x(t+)xT(t) >< x(t)xT(t) >-1 .
Eigenvectors of G() are the “normal” modes {ui}.
Most probable prediction: x’(t+) = G() x(t)
The neat thing: G() ={G() }
![Page 3: El Ni ño and how to get rid of it](https://reader035.vdocument.in/reader035/viewer/2022081519/56812d15550346895d91fe55/html5/thumbnails/3.jpg)
SST Data used:
• COADS (1950-2000) SSTs in 30E-70W, 30N – 30S consolidated onto a 4x10-degree grid.
• Subjected to 3-month running mean.• Projected onto 20 EOFs (eigenvectors of <xxT>)
containing 66% of the variance.• x, then, represents the vector of SST anomalies,
each component representing a location, or else it represents the vector of Principal Components.
![Page 4: El Ni ño and how to get rid of it](https://reader035.vdocument.in/reader035/viewer/2022081519/56812d15550346895d91fe55/html5/thumbnails/4.jpg)
El Niño can be described this way.
If LIM is successful, prediction error does not depend on the lag at which the covariance matrices are evaluated. This is true for El Niño; it is not true for the chaotic Lorenz system. Below, different colors correspond to different lags used to identify the parameters.
![Page 5: El Ni ño and how to get rid of it](https://reader035.vdocument.in/reader035/viewer/2022081519/56812d15550346895d91fe55/html5/thumbnails/5.jpg)
The annual cycle:
dx/dt = Bx + (t)T (t) = Q(t)
Given stationary B use (time-dependent) conservation of variance to diagnose Q(t).
Result: The annual cycle of Q(t) looks nothing like the phase locking of either El Niño or the optimal structure to the annual cycle.
But
A model generated with the stationary B and the stochastic forcing with cyclic statistics Q(t) does reproduce the correct phase-locking in both.
![Page 6: El Ni ño and how to get rid of it](https://reader035.vdocument.in/reader035/viewer/2022081519/56812d15550346895d91fe55/html5/thumbnails/6.jpg)
What do models say?
• Hybrid coupled model (Chang 1994)• Dynamical core of GFDL Ocean model,
considerably simplified• Statistical atmosphere based on EOFs of
observed wind stress• Interactive annual cycle• Strength of coupling determines dynamical
regime (Note: an artificial parameter )
![Page 7: El Ni ño and how to get rid of it](https://reader035.vdocument.in/reader035/viewer/2022081519/56812d15550346895d91fe55/html5/thumbnails/7.jpg)
![Page 8: El Ni ño and how to get rid of it](https://reader035.vdocument.in/reader035/viewer/2022081519/56812d15550346895d91fe55/html5/thumbnails/8.jpg)
![Page 9: El Ni ño and how to get rid of it](https://reader035.vdocument.in/reader035/viewer/2022081519/56812d15550346895d91fe55/html5/thumbnails/9.jpg)
![Page 10: El Ni ño and how to get rid of it](https://reader035.vdocument.in/reader035/viewer/2022081519/56812d15550346895d91fe55/html5/thumbnails/10.jpg)
Optimal initial structure for growth over lead time :
Right singular vector of G() (eigenvector of GTG())
Growth factor over lead time :
Eigenvalue of GTG().
![Page 11: El Ni ño and how to get rid of it](https://reader035.vdocument.in/reader035/viewer/2022081519/56812d15550346895d91fe55/html5/thumbnails/11.jpg)
(a)
The transient growth possible in a multidimensional linear system occurs when an El Niño develops. LIM predicts that an optimal pattern (a) precedes a mature El Niño pattern (b) by about 8 months.
(b)
![Page 12: El Ni ño and how to get rid of it](https://reader035.vdocument.in/reader035/viewer/2022081519/56812d15550346895d91fe55/html5/thumbnails/12.jpg)
Does it? Judge for yourself! (c and d)
c)d)
(c) (d)
In (c), the red line is the time series of pattern correlations between pattern (a) and the sea surface temperature pattern 8 months earlier. The blue line is a time series index of how strong pattern (b) is at the date shown; the blue line is an index of El Niño when it is positive and of La Niña when it is negative. In (d) we see a scatter plot of the El Niño index and pattern correlations shown in (c).Pattern (a) really does precede El Niño! Pattern (a) with the signs reversed precedes La Niña!
![Page 13: El Ni ño and how to get rid of it](https://reader035.vdocument.in/reader035/viewer/2022081519/56812d15550346895d91fe55/html5/thumbnails/13.jpg)
![Page 14: El Ni ño and how to get rid of it](https://reader035.vdocument.in/reader035/viewer/2022081519/56812d15550346895d91fe55/html5/thumbnails/14.jpg)
Decay mode, = 31 months
0
0.5
1
1.5
0 5 10 15 20 25Mode number
momo
momo
momo
decay timeT = Period
Projection of adjoints onto O.S. and modal timescales.
![Page 15: El Ni ño and how to get rid of it](https://reader035.vdocument.in/reader035/viewer/2022081519/56812d15550346895d91fe55/html5/thumbnails/15.jpg)
Location of indices: N3.4, IND, NTA, EA, and STA.
![Page 16: El Ni ño and how to get rid of it](https://reader035.vdocument.in/reader035/viewer/2022081519/56812d15550346895d91fe55/html5/thumbnails/16.jpg)
-3
-2
-1
0
1
2
3
1960 1970 1980 1990 2000
Nino3.4
months
-3
-2
-1
0
1
2
3
1960 1970 1980 1990 2000
Filtered Nino3.4
months
![Page 17: El Ni ño and how to get rid of it](https://reader035.vdocument.in/reader035/viewer/2022081519/56812d15550346895d91fe55/html5/thumbnails/17.jpg)
10-4
10-3
10-2
10-1
100
101
1 10 100 1000
Period (months)
![Page 18: El Ni ño and how to get rid of it](https://reader035.vdocument.in/reader035/viewer/2022081519/56812d15550346895d91fe55/html5/thumbnails/18.jpg)
-5
0
5
10
15
20
25
1 10 100 1000
Period (months)
66.4 mo39.9 mo
18.1 mo
15.3 mo
![Page 19: El Ni ño and how to get rid of it](https://reader035.vdocument.in/reader035/viewer/2022081519/56812d15550346895d91fe55/html5/thumbnails/19.jpg)
-505
101520253035
100 101 102 103 104
43.9 mo
16.5 mo
5.2 mo
Period (weeks)
![Page 20: El Ni ño and how to get rid of it](https://reader035.vdocument.in/reader035/viewer/2022081519/56812d15550346895d91fe55/html5/thumbnails/20.jpg)
-1
-0.5
0
0.5
1
1950 1960 1970 1980 1990 2000Date
R(Unfiltered, El Nino) = 0.36
-1.5
-1
-0.5
0
0.5
1
1.5
1950 1960 1970 1980 1990 2000Date
R(Unfiltered, El Nino) = 0.44
-1.5
-1
-0.5
0
0.5
1
1.5
1950 1960 1970 1980 1990 2000Date
R(Unfiltered, El Nino) = 0.45
-1.5
-1
-0.5
0
0.5
1
1.5
1950 1960 1970 1980 1990 2000
R(Unfiltered, El Nino) = 0.61
EA
ST
AIN
D
NT
A
R = 0.36 R = 0.45
R = 0.44 R = 0.61
Indices. Black: Unfiltered data. Red: El Niño signal.
![Page 21: El Ni ño and how to get rid of it](https://reader035.vdocument.in/reader035/viewer/2022081519/56812d15550346895d91fe55/html5/thumbnails/21.jpg)
-0.8
-0.6
-0.4
-0.2
0
0.2
-100 -50 0 50 100Lead (months)
8 months
STA leads PC1 leads
-0.2
-0.1
0
0.1
0.2
0.3
0.4
0.5
0.6
-100 -50 0 50 100Lead (months)
IND leads PC1 leads
-0.6
-0.4
-0.2
0
0.2
0.4
-100 -50 0 50 100Lead (months)
EA leads PC1 leads
9 months
-0.2
-0.1
0
0.1
0.2
0.3
0.4
0.5
-100 -50 0 50 100Lead (months)
NTA leads PC1 leads
Lagged correlation between El Niño indices and PC 1.
STA leads PC1 leads PC1 leads
PC1 leads PC1 leads
EA leads
IND leads NTA leads
![Page 22: El Ni ño and how to get rid of it](https://reader035.vdocument.in/reader035/viewer/2022081519/56812d15550346895d91fe55/html5/thumbnails/22.jpg)
Decay mode, = 31 months
0
0.5
1
1.5
0 5 10 15 20 25Mode number
momo
momo
momo
decay timeT = Period
Projection of adjoints onto O.S. and modal timescales.
![Page 23: El Ni ño and how to get rid of it](https://reader035.vdocument.in/reader035/viewer/2022081519/56812d15550346895d91fe55/html5/thumbnails/23.jpg)
EOF 1 of Residual
-15
-10
-5
0
5
10
1950 1960 1970 1980 1990 2000 2010Date
u1 of un-filtered data
The pattern correlation between the longest-lived mode of the unfiltered data and the leading EOF of the residual data is 0.81.
![Page 24: El Ni ño and how to get rid of it](https://reader035.vdocument.in/reader035/viewer/2022081519/56812d15550346895d91fe55/html5/thumbnails/24.jpg)
-1
-0.5
0
0.5
1
1950 1960 1970 1980 1990 2000
R(Unfiltered, El Nino +Trend) = 0.75
-1.5
-1
-0.5
0
0.5
1
1.5
1950 1960 1970 1980 1990 2000Date
R(Unfiltered, El Nino+Trend) = 0.79
-1.5
-1
-0.5
0
0.5
1
1.5
1950 1960 1970 1980 1990 2000Date
R(Unfiltered, El Nino+Trend) = 0.77
-1.5
-1
-0.5
0
0.5
1
1.5
1950 1960 1970 1980 1990 2000
R(Unfiltered, El Nino + Trend) = 0.62
EA
SS
TA
(C
)
ST
A S
ST
A (
C)
IND
SS
TA
(C
)
NT
A S
ST
A (
C)
R = 0.75 R = 0.77
R = 0.79 R = 0.62
Indices. Black: Unfiltered data. Green: El Niño signal + Trend.
![Page 25: El Ni ño and how to get rid of it](https://reader035.vdocument.in/reader035/viewer/2022081519/56812d15550346895d91fe55/html5/thumbnails/25.jpg)
All this results from SST dynamics being essentially linear.
But linear dynamics implies symmetry between El Niño and La Niña events.
SST anomalies appear to be positively skewed.
Is the skew significant?
![Page 26: El Ni ño and how to get rid of it](https://reader035.vdocument.in/reader035/viewer/2022081519/56812d15550346895d91fe55/html5/thumbnails/26.jpg)
Additive Noise Model
dx/dt = Bx +
Multiplicative Noise Model 1
dx/dt=(B*+A*)x +
Multiplicative Noise Model 2
dx/dt=B’(I+I’)x +
![Page 27: El Ni ño and how to get rid of it](https://reader035.vdocument.in/reader035/viewer/2022081519/56812d15550346895d91fe55/html5/thumbnails/27.jpg)
Conclusions• El Niño appears to be a damped system forced by
stochastic noise.• There is evidence that the phase locking of El
Niño to the annual cycle is due to the annually-varying statistics of the stochastic noise.
• An optimal initial structure for growth precedes a mature El Niño event by 6 to 9 months.
• The nonnormal dynamics are dominated by 3 nonorthogonal modes.
• The essential linearity of the system allows isolation of the El Niño signal.
![Page 28: El Ni ño and how to get rid of it](https://reader035.vdocument.in/reader035/viewer/2022081519/56812d15550346895d91fe55/html5/thumbnails/28.jpg)
• El Niño indices in the equatorial Indian Ocean and the North Tropical Atlantic Ocean are very similar.
• El Niño and the (parabolic) trend dominate the SST variability in the Indian Ocean as well as in the Equatorial and S. Tropical Atlantic Ocean.
• The observed skew in El Niño –La Niña events IS NOT significant compared with realistic null hypothses.
• The trend IS significant compared with realistic null hypotheses.
![Page 29: El Ni ño and how to get rid of it](https://reader035.vdocument.in/reader035/viewer/2022081519/56812d15550346895d91fe55/html5/thumbnails/29.jpg)
![Page 30: El Ni ño and how to get rid of it](https://reader035.vdocument.in/reader035/viewer/2022081519/56812d15550346895d91fe55/html5/thumbnails/30.jpg)