Climate Teleconnections with a Time Lag

Original: "Lagged teleconnections of climate variables identified via complex rotated Maximum Covariance Analysis"
arXiv:2105.04618 · 2021-05-10 · CC BY 4.0 · 3 min
A new method hears the pauses in the ocean–land conversation, even when the reply comes months later.

In 1960, a meteorologist at a station in Peru records water temperature, while a farmer in India watches the sky waiting for the monsoon. There is no direct link between these two events, but months later one will echo in the other. Earth's climate is like a conversation in which the speakers do not reply immediately but after a pause: the ocean warms, and only months later does the rainfall pattern change on another continent. Old methods only heard simultaneous cues, so half the conversation remained unintelligible.

Classical maximum covariance analysis looks for common patterns in two data fields but assumes changes happen at the same time. If the signal lags, the method misses it.

The authors of the new study decided to teach statistics to hear pauses. They turned each time series of temperature and precipitation into an analytic signal using the Hilbert transform, which shifts oscillations by a quarter period. This allowed them to account for the fact that the peak of warmth in one region may occur earlier than the peak of rain in another. To avoid distortions at the edges of the sixty-year record, they extended the data into the past and future using the Theta model, then trimmed the excess. After that, they computed the covariance matrix between the two complex fields and applied singular value decomposition. The first 150 modes, explaining 99.82% of the covariance, were subjected to Varimax rotation to make the patterns more local and interpretable.

The method is implemented in the Python package xmca and is ready for use in finding delayed teleconnections between any climate fields.

The result resembled a transcription of a polyphonic choir. The method identified six main modes of coupled ocean–precipitation variability. The first mode (62.7% of covariance) describes the seasonal cycle: the phase function allowed calculating that the seasonal temperature maximum in the North Pacific follows the South with a delay of 196 days, while in the Indian Ocean precipitation lags temperature by 81 days. The second mode (7.0%) corresponds to El Niño: the correlation with the ONI index was 0.941. The third mode (6.9%) reflects global warming with a correlation of 0.781 with the ocean warming index. The fourth (1.6%) reproduces the Pacific Decadal Oscillation (correlation 0.565), the fifth (1.6%) El Niño Modoki (0.506), and the sixth (1.1%) the Atlantic Meridional Mode (0.741).

Previously, these teleconnections were found individually in scattered regional studies, and each time one had to assume in advance which delay mattered. The new method is the first to combine a complex signal representation with Varimax rotation within maximum covariance analysis, extracting all major modes together with their time shifts in a single procedure. It does not require prior trend removal, which is critical in a changing climate.

Now the ocean–land conversation can be heard in full, with all its pauses. The farmer in India and the meteorologist in Peru are connected by an invisible thread of delayed signal, and the new method teaches us to discern it amid the noise. Ahead lies application to the Madden–Julian oscillation and other variable pairs, from ocean salinity to soil moisture. Perhaps we will hear new lines in this ancient dialogue.

The seasonal temperature maximum in the northern Pacific occurs 196 days later than in the southern Pacific, which is almost exactly half a year. This means the two ocean hemispheres work in antiphase, like a seesaw with a very slow period.

The idea of delayed teleconnections recalls the 'butterfly effect' from Ray Bradbury's story 'A Sound of Thunder': a small event in one place, through a chain of connections, leads to major consequences elsewhere. Only here the connections are not chaotic but statistically stable.

\hat{X}_s = X_s + i H(X_s)
Each time series is turned into a complex one, where the imaginary part is phase-shifted by a quarter period
A_s = (V_s \odot \bar{V}_s)^{\odot 1/2}
Shows in which regions a given mode is most pronounced
\Theta_s = \arctan2(\mathrm{Im}(V_s), \mathrm{Re}(V_s))
Gives the relative temporal shift at each point; for narrowband signals it translates into a delay
Concepts
Teleconnection reanalysis Atlantic multidecadal oscillation multifractal detrended fluctuation analysis ocean-atmosphere coupling decadal variability ERA5 climate network North Atlantic Oscillation climate marine heatwave covariance matrix coupled general circulation model seasonal forcing seasonal cycle empirical orthogonal function mixed mode oscillations GODAS natural variability principal component analysis