Next Article in Journal
Extreme Months: Multidimensional Studies in the Carpathian Basin
Previous Article in Journal
Lichen Biodiversity Index (LBI) for the Assessment of Air Quality in an Industrial City in Pahang, Malaysia
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Nature and Origin of Atmospheric Annual and Semi-Annual Oscillations

by
Vincent Courtillot
1,
Jean-Louis Le Mouël
1,
Fernando Lopes
1,* and
Dominique Gibert
2
1
Institut de Physique du globe de Paris, Université de Paris, 75005 Paris, France
2
LGL-TPE-Laboratoire de Géologie de Lyon-Terre, Planètes, Environnement, 69622 Lyon, France
*
Author to whom correspondence should be addressed.
Atmosphere 2022, 13(11), 1907; https://doi.org/10.3390/atmos13111907
Submission received: 30 August 2022 / Revised: 19 October 2022 / Accepted: 14 November 2022 / Published: 15 November 2022
(This article belongs to the Section Meteorology)

Abstract

:
This paper proposes a joint analysis of variations of global sea-level pressure (SLP) and of Earth’s rotation (RP), expressed as the coordinates of the rotation pole ( m 1 , m 2 ) and length of day (lod). We retain iterative singular spectrum analysis (iSSA) as the main tool to extract the trend, periods, and quasi periods in the data time series. SLP components are a weak trend, seven quasi-periodic or periodic components (∼130, 90, 50, 22, 15, 4, 1.8 years), an annual cycle, and its first three harmonics. These periods are characteristic of the space-time evolution of the Earth’s rotation axis and are present in many characteristic features of solar and terrestrial physics. The amplitudes of the annual SLP component and its three first harmonics decrease from 93 hPa for the annual to 21 hPa for the third harmonic. In contrast, the components with pseudo-periods longer than a year range between 0.2 and 0.5 hPa. We focus mainly on the annual and, to a lesser extent, the semi-annual components. The annual RP and SLP components have a phase lag of 152 days (half the Euler period). Maps of the first three components of SLP (that together comprise 85% of the data variance) reveal interesting symmetries. The trend is very stable and forms a triskeles structure that can be modeled as Taylor–Couette flow of mode 3. The annual component is characterized by a large negative anomaly extending over Eurasia in the NH summer (and the opposite in the NH winter) and three large positive anomalies over Australia and the southern tips of South America and South Africa in the SH spring (and the opposite in the SH autumn), forming a triskeles. The semi-annual component is characterized by three positive anomalies (an irregular triskeles) in the NH spring and autumn (and the opposite in the NH summer and winter), and in the SH spring and autumn by a strong stable pattern consisting of three large negative anomalies forming a clear triskeles within the 40–60 annulus formed by the southern oceans. A large positive anomaly centered over Antarctica, with its maximum displaced toward Australia, and a smaller one centered over Southern Africa, complement the pattern. Analysis of iSSA components of global sea level pressure shows a rather simple spatial distribution with the principal forcing factor being changes in parameters of the Earth’s rotation pole and velocity. The flow can probably best be modeled as a set of coaxial cylinders arranged in groups of three (triskeles) or four and controlled by Earth topography and continent/ocean boundaries. Flow patterns suggested by maps of the three main iSSA components of SLP (trend, annual, and semi-annual) are suggestive of Taylor–Couette flow. The envelopes of the annual components of SLP and RP are offset by four decades, and there are indications that causality is present in that changes in Earth rotation axis lead force pressure variations.

1. Introduction

This paper is an attempt to obtain better constraints on the forcings of the trend and annual components of both global sea-level pressure and variations in the Earth’s rotation and to test the hypothesis that there might be a causal link between them.
Lopes et al. [1] studied the singular spectral analysis (SSA) of the evolution of mean sea-level atmospheric pressure (SLP) since 1850. In addition to dominant trends, eleven quasi-periodic components were identified, with (pseudo-) periods of ∼130, 90, 50, 22, 15, 4, 1.8, 1, 0.5, 0.33, and 0.25 years, corresponding to the Jose, Gleissberg, Hale, and Schwabe cycles, to the annual cycle, and its first three harmonics. These periods are already known to be characteristic of the space-time evolution of the Earth’s rotation axis: the rotation pole (RP) also undergoes periodic motions longer than 1 year, at least up to the Gleissberg ∼90 years cycle [2,3,4,5,6,7,8,9,10,11,12,13,14,15]. They are encountered in solar physics [16,17,18,19,20,21,22,23,24,25] and terrestrial climate [26,27,28,29,30,31,32,33,34,35].
The rotation velocity is usually expressed as the length of day (lod); it contains periods of 1 year and shorter, but also (although they are weaker) longer periods from 11 years to 18.6 years [36,37,38,39,40]. Lopes et al. [15] showed that the secular variation of lod since 1846 [41] and RP have been carried by an oscillation whose period fits the Gleissberg cycle.
RP consists to first order of three SSA components that together amount to 85% of the total signal variance [12]: the drift [4], the free oscillation [2,3] (actually a double component with periods ∼433 and ∼434 days), and the forced annual oscillation [7]. Polar drift does not have a universally accepted explanation [5,9,42,43]. The two other oscillatory components (Chandler and annual) are obtained using the Liouville–Euler system, which describes the motion of a spherical rotating solid body (e.g., [7], chapter 3, system 3.2.9).
Polar drift was discussed in [13,15], and the free Chandler wobble was addressed in [14]. In the present paper, we focus on the third major component, the forced annual oscillation. This oscillation is often assumed to be forced by the Earth’s fluid envelopes, hence also by climate variations (e.g., [7,44], chapter 7, [45,46]). The fact that the observed Chandler period (433–434 days) is not equal to the theoretical value of the Euler period (306 days) has been interpreted in terms of Earth elasticity. The annual period is forced by interactions with fluids. Variations in the Sun–Earth distance, hence of the corresponding gravitational forces, displace the fluid atmosphere and ocean; exchanges in angular moment affect RP and lod. The annual oscillation of the rotation pole RP (coordinates m 1 and m 2 ) would therefore be due to the presence of the fluid envelopes. If Earth were devoid of fluid envelopes, the annual oscillation should not exist, which would contradict both the theory and observations of motion of a Lagrange top [47]. Wilson et Haubrich [48] recall that Spitaler [49,50] demonstrated that a part of the forcing of the annual wobble was due to the migration of air masses on and off the Asian continent. The authors in [51] showed that the annual fluctuation in water storage on the continents was also important. Study [52] re-examined the sources of annual wobble excitation and concluded that the air mass effect accounted for much but not all of the annual wobble, and water storage did not explain the remainder. Reservations appear in chapter 7 of Lambeck ’s book [7], Seasonal Variations, page 146, who writes that since Jeffreys, in 1916, the seasonal oscillation has been attributed to a geographical redistribution of mass associated with meteorological causes (e.g., atmospheric and oceanic motion, precipitation, vegetation, polar ice, etc.) but, although this conclusion is still valid today, it presents some important discrepancies when we compared observed and computed annual components. In summary, the interpretation of the forced annual component of polar motion is still hypothetical and fails to be validated by a numerical model.
In a previous analysis of global sea-level pressure (SLP), we found that trends since 1950 were very stable in time and space [1] and were organized in a dominant 3-fold symmetry about the rotation axis in the northern hemisphere (NH) and a 3 or 4-fold symmetry in the southern hemisphere (SH) [1]. These features could be interpreted as resulting from Taylor–Couette flow. In this paper, we return to the pressure data and focus on the annual oscillation.

2. The Pressure Data and Method of Analysis

The pressure data are maintained by the Met Office Hadley Centre (https://www.metoffice.gov.uk/hadobs/hadslp2/data/download.html, accessed on 10 October 2020) and can be accessed as maps of global pressure, every month from 1850 up to the present. The sampling is 5 × 5 ; the data we use are labeled HadSLP2r. Allan and Ansell ([53]) built this series starting with ground observations whose number and location are shown in their figure 2. These are subjected to a quality check, corrected for various local effects, then homogenized with the empirical mode decomposition (EMD) filter.
As in [1], we have used the iterative singular spectrum analysis algorithm (iSSA) to extract the components, mainly the annual one, in each grid cell and as a function of time. For the iSSA method, see [54]; for properties of the Hankel and Toeplitz matrices, see [55]; and for the singular value decomposition (SVD) algorithm, see [56]. We will now summarize the four steps of the SSA method.
Let us consider a discrete (non zero) time series ( X N ) of length N (N > 2):
X N = ( x 1 , , x N ) .

2.1. Step 1 (Embedding Step)

The time series X N is divided into K segments of length L in order to build a matrix X with dimension K × N , where K = N L + 1 will condition our decomposition. This is the first “tuning knob”. Integrating X yields a Hankel matrix:
X = x 1 x 2 x 3 x K x 2 x 3 x 4 x K + 1 x 3 x 4 x 5 x K + 2 x L x L + 1 x L + 2 x N
Embedding, the first step in an SSA, consists of projecting the one-dimensional time series in a multidimensional space of series X N such that vectors X i = ( x i , , x i + L 1 ) t belong to R L , where K = N L + 1 . The parameter that controls the embedding is L, the size of the analyzing window, an integer between 2 and N 1 . The Hankel matrix has a number of symmetry properties: its transpose X t , called the trajectory matrix, has dimension K. Embedding is a compulsory step in the analysis of non-linear series. It consists formally in the empirical evaluation of all pairs of distances between two offset vectors, delayed (lagged) in order to calculate the correlation dimension of the series.

2.2. Step 2 (Decomposition in Singular Values (SVD))

The SVD [56] of non-zero trajectory matrix X (dimensions L × K ) takes the shape:
X = i = 1 d λ i U i V i t
where the eigenvalues λ i ( i = 1 , , L ) of matrix S = XX T are arranged in order of decreasing amplitudes. Eigenvectors U i and V i are given by:
V i = X T U i / λ i
The V i form an orthonormal basis and are arranged in the same order as the λ i . Let X i be a part of matrix X such that:
X i = λ i U i V i t .
Embedding matrix X can then be represented as a simple linear sum of elementary matrices X i . If all eigenvalues are equal to 1, then decomposition of X into a sum of unitary matrices is:
X = X 1 + X 2 + + X d
with d being the rank of X ( d = rank X = m a x { i | λ i > 0 } ). SVD allows one to write X as a sum of d unitary matrices, defined in a univocal way.
Let us now discuss the nature and the characteristics of the embedding matrix. Its rows and columns are sub-series of the original time series (or signal). The eigenvectors U i and V i have a time structure, and they can be considered as a representation of temporal data. Let X be a suite of L lagged parts of ( X and X 1 , , X K ), the linear basis of its eigenvectors. If we let
Z i = i = 1 d λ i V i
with i = 1 , , d , then relation (5) can be written:
X = i = 1 d U i Z i t
that is, for the jth elementary matrix:
X j = i = 1 d z j i U i
where z j i is a component of vector Z i . This means that vector Z i is composed of the ith components of vector X j . In the same way, if we let
Y i = i = 1 d λ i U i
we obtain for the transposed trajectory matrix:
X j t = i = 1 d U i Y i t
that corresponds to a representation of the K lagged vectors in the orthogonal basis ( V 1 , …, V d ). One sees why SVD is a very good choice for the analysis of the embedding matrix, as it provides two different geometrical descriptions.

2.3. Step 3 (Reconstruction)

As we have seen, X i matrices are unit matrices, and one can “re-group” these matrices into a physically homogeneous quantity. This is the second “tuning knob” of SSA. In order to regroup the unit matrices, one divides the set of indices i { 1 , , d } into m disjoint subsets of indices { I 1 , , I m } .
Let I be the grouping of p indices of I = { i 1 , i 2 , , i p } ; because (6) is linear, then the resulting matrix X I that regroups indices I can be written:
X I = X I 1 + X I 2 + + X I m
This step is called regrouping the eigen-triplets ( λ , U, and V). In the limit case m = d , (12) becomes exactly (6), and we find again the unit matrices.
Next, how can one associate pairs of eigen-triplets? This means separating the additive components of a time series. One must first consider the concept of separability.
Let X be the sum of two time series X ( 1 ) and X ( 2 ) such that x i = x i ( 1 ) + x i ( 2 ) for any i [ 1 , N ] . Let L be the analyzing window (with fixed length), X, X ( 1 ) , and X ( 2 ) the embedding matrices of series X , X ( 1 ) , and X ( 2 ) . These two sub-series are separable (even weakly) in Equation (6) if there is a collection of indices I { 1 , , d } such that X ( 1 ) = i I X i , respectively, if there is a collection of indices such that X ( 1 ) = i I X i .
In the case when separability does exist, the contribution of X ( 1 ) , for instance, corresponds to the ratio of associated eigenvalues ( i I λ i ) to total eigenvalues ( i = 1 d λ i ).
Therefore, regrouping SVD components can be summarized by the decomposition into several elementary matrices, whose structure must be as close as possible to a Hankel matrix of the initial trajectory matrix (this is true on paper only; in reality things are much more difficult).

2.4. Step 4 (The Diagonal Mean, Aka the Hankelization Step)

The next, final step consists in going back to data space, that is, to calculate time series with length N associated with sub-matrices X I . Let Y be a matrix with dimension L * K and for each element y i j we have 1 i L and 1 j K . Let L * be the minimum and K * be the maximum. One always has N = L + K 1 . Finally, let y i j * = y i j if L < K and y i j * = y j i otherwise. The diagonal average applied to kth index of time series y associated with matrix Y gives:
y k = 1 k m = 1 k y m , k m + 1 * 1 k L * 1 L * m = 1 L * y m , k m + 1 * L * k K * 1 N K + 1 m = k K * + 1 N K * + 1 y m , k m + 1 * K * k N *
The relation (13) corresponds to the mean of the element on the anti-diagonal i + j = k + 1 of the matrix. For k = 1 , y 1 = y 1 , 1 . For k = 2 , y 2 = ( y 1 , 2 + y 2 , 1 ) / 2 , etc.
As noted above, step 3 is the most difficult part. We have chosen one approach among many others: iterative SSA (iSSA). Because relation (6) is linear, we can iterate the decomposition. We start with a small value of L (we are looking for the longest period) that we increase until getting a quasi-Hankel matrix (step 1 and 2). We then extract the corresponding lowest frequency component that it subtracted from the original signal. We increase again the value of L to find the next component (shortest period). The algorithm stops when no pseudo-cycle can be detected or extracted. In this way, we scan the series from low to high frequencies.

3. The SSA Pressure Components

We obtain the following components:
The trend (from ∼1009.15 hPa in 1850 to ∼1008.35 hPa at present) is the first and largest iSSA component. It represents more than 70% of the total variance of the original series. The sequence of the quasi-periodic components is, in decreasing order of periods (amplitudes are in hectoPascal, hPa):
  • ∼130 years (∼0.7 hPa). Compatible with the Jose [17] cycle,
  • ∼90 years (∼21 hPa). We recognize the Gleissberg [16] cycle [21],
  • ∼50 years (∼0.2 hPa),
  • ∼22 years (∼0.50 hPa). We recognize the Hale cycle [57],
  • ∼15 years (∼0.2 hPa). We may recognize an upper bound of the Schwabe cycle [58],
  • ∼4 years (∼0.3 hPa),
  • ∼1.8 years (∼0.3 hPa).
Then:
  • 1 year ( ∼93 hPa),
  • 0.5 years (∼65 hPa),
  • 0.33 years (∼44 hPa),
  • 0.25 years (∼21 hPa).

4. Annual and Semi-Annual SSA Pressure Components and Variations in Polar Motion: Time Analysis

Many studies (e.g., [7,51,59]) attempt to relate forced (annual) oscillations to polar motion RP ( m 1 , m 2 ) or length of day (lod) (data file EOPC01IAU2000, [60], maintained by IERS, https://www.iers.org/IERS/EN/DataProducts/EarthOrientationData/eop.html, accessed on 12 October 2020). The data are shown and discussed in [14]. The annual iSSA components of m 1 and m 2 are displayed in Figure 1. Both undergo some amount of modulation. That modulation is significant for m 1 , particularly between 1850 and 1930. For m 2 , it is smaller and more regular. An m 1 vs. m 2 diagram (see Appendix A) shows the classical Lissajou elliptical shape of annual polar motion. Note that the date of 1930 is when the Chandler wobble undergoes a phase jump of π (e.g., [61,62]). The lod data start only in 1962 (not shown, see [15]).
The annual and semi-annual iSSA components of lod are modulated (Figure 2a,b). The annual component grows slightly until 1990, flattens until 2013, and then starts growing again. The semi-annual component decreases regularly over the sixty years of available data.
Figure 2. Annual and semi-annual components extracted from length of day: (a) iSSA annual component of length of day from 1850 to the present; (b) same as (a) but for the semi-annual component.
Figure 2. Annual and semi-annual components extracted from length of day: (a) iSSA annual component of length of day from 1850 to the present; (b) same as (a) but for the semi-annual component.
Atmosphere 13 01907 g002aAtmosphere 13 01907 g002b
Figure 3 shows the annual iSSA component since 1850, extracted from the SLP pressure series. Its modulation looks very much like that of m 2 , with a phase shift (Figure 1). We have checked that the phase shift is constant over the 70 years of data. Its mean value is 152.31 ± 2.68 days. We have already obtained (exactly) such a phase shift and note that it is (exactly) half of the Euler period of 306 days: a reason for this remains to be found.
Using the Hilbert transform, we have determined the envelopes of the two annual components of polar motion and atmospheric pressure (Figure 4a). Polar motion clearly leads pressure by at least 40 years. In Figure 4b, we see that the trend of polar motion is close to the envelope of pressure, in general preceding it slightly. Given the mathematical properties of the Liouville–Euler set of equations, the potential direction of causality between the two phenomena is in the sense of polar motion leading atmospheric pressure. If correct, this is an important result that confirms Laplace’s statements in his 1799 Traité de Mécanique Céleste (volume 5, chapter 1, page 347) [63].
Thus, after applying the Kepler’s law of conservation of areas, Laplace concludes that atmospheric motions do not affect pole rotation.

5. The Annual and Semi-Annual SSA Pressure Components: Spatial Analysis

We show in Figure 5 polar maps of the mean amplitude of the annual component (oscillation) for each season since 1850, and in Figure 6a–d polar maps of the mean of the semi-annual component.
In Figure 6a–d, in the spring SH, there is, in addition, a very regular pattern at the intersections of the 20 S to 50 S annulus and the southern tips of the three southern continents (Africa, Australia, and South America). There is a smaller dipole with its negative part centered almost perfectly on the South Pole and its positive part centered on 70 S, 105 W. In the fall, the pattern is the same but with a reversal in sign. In the summer, there is a positive band between 15 and 25° S latitude from Madagascar to New Caledonia (40°–170° E) and a small negative spot near the South Pole (winter is the same with a sign reversal). The positive annulus is actually also seen, much weaker, in the Atlantic and Pacific oceans. The NH is rather different in the spring, with a large negative area over Asia from the Red Sea to Japan and extending in latitude from 80 N to 15 N. This large “anomaly” (using geophysical language) extends negative arms towards Senegal and over the western USA. Over the year, the “anomaly” changes sign. Still, in the summer, there are two sizable positive anomalies over the northern Atlantic and Pacific Oceans between 40 N and 60 N. In the SH, the large mean annual pressure variations occur over the southern continents; in the NH, the Asian continental “anomaly” dominates, although two significant features occur over the northern parts of the oceans. Because the Asian “anomaly” and the two (opposite sign) anomalies in the North Atlantic and Pacific are separated by almost 180° in longitude, their maxima being always in phase, their contributions to the seasonal excitation function tend to cancel [64].
The semi-annual component (see Figure 7) has stronger patterns of symmetry. The SH features a very stable 3-fold symmetry (triskeles); three large anomalies form an equilateral triangle centered on 70 E, 30 W, and 210 W, between latitudes 40 S and 60 S, a large anomaly of opposite sign is centered on the South Pole (the maximum is actually offset from the SP by 20 E towards the 120 E meridian), and an anomaly with the same sign as that over the pole, but located over southern Africa. The NH has weaker anomalies forming an irregular triangle with apices near the North Pacific (40 N, 165 E), Baffin Sea (40 N, 50 W), and Iran (30 N, 60 W).
The bottom two maps in Figure 8b–c are the spring maps for the annual and semi-annual components.
It is interesting to compare the maps of Figure 8 with the maps of mean pressure gradients shown in Figure 2 of Lopes et al. [1]. The zonal structure appears clearly, with sharp changes in the sign of anomalies at continent/ocean edges. The zonal structure is reminiscent of the bands seen in the atmosphere of the giant planets; the continent/ocean boundaries interrupt the zonal belts, leading to a segmented structure.
We can compare the maps of Figure 5, Figure 6, Figure 7 and Figure 8 for the annual, semi-annual, and trend components with corresponding figures from the reference treatise of [64]. Recall that, taken together, these three components suffice to capture the total signal variance. Figures 4–12 and Figures 4–13 ([64], pages 156–157 and reproduced here Figure 5) show Lamb’s representation of the average mean sea level pressure over (respectively) the northern and southern hemispheres in the 1950s.
Figure 5 shows an overlay of the SSA triskeles patterns (trends) and annual oscillation of SLP from Lopes et al. ([1]) on top of maps from Lamb (1972, page 157, Figures 4–13). The left column is for the northern hemisphere, the right one for the southern hemisphere. Lamb’s maps are calculated for a time span of 40 years; our SSA results are for 170 years. There is very good agreement for the northern hemisphere: the January map from Lamb and our winter map show the dominant “dipole” with the large high pressure (HP) over Asia and the low pressure (LP) over the northern Pacific. The April and spring maps do not agree as well. The July map from Lamb and our summer map both show a large LP over Asia and two HPs over the northern Pacific and northern Atlantic. The October map from Lamb and our autumn map show a weaker HP over Asia. For the southern hemisphere, agreement between the two sets of maps is significantly lower; the winter vs. January maps both feature a HP centered on Antarctica and a LP extending from Madagascar to Australia, and the reverse in summer/July. Our spring (respectively, autumn) maps feature strong HPs (respectively, LPs) on the southern tips of the southern continents. This pattern is not as sharp in Lamb’s (1972) maps.
Lamb’s Figures 3–17 ([64] and Figure 5) represent the annual mean distribution of atmospheric pressure at sea level for the time span 1900–1940 for the northern hemisphere and 1900–1950s for the southern hemisphere. This can be compared to our maps of the trend ([1], Figure 6 top row and Figure 7). The agreement is quite good: the 3-fold symmetry of the northern hemisphere and the 4-fold symmetry of the southern hemisphere are even clearer in the [1] maps using the iSSA method to extract the trend (i.e., the main component) from the data. As concluded by [64], the effects of geography modify the circumpolar circulation much more in the northern than in the southern hemisphere; “The broad zonal character of the long-period mean circulation is as permanent as the circumpolar arrangement of the heating and cooling zones and the circumpolar vortex aloft”. The patterns we obtain in 2020 and those Lamb (1972) [64] obtained 50 years before are so similar that one can conclude the multidecadal stability of the atmospheric circulation. Our maps are based on 170 vs. 40 years of data, which lends credence to the sharper and simpler features of our maps. It appears to us that such maps from the southern hemisphere are particularly interesting, and so is the triskeles in Figure 5 (top left). Specialists may wish to analyze these maps in more detail; they can be provided on request.

6. Discussion

6.1. Symmetries and Forcings

We pointed out in a preliminary paper on symmetries of global sea level pressure (Lopes et al. (2022b) [1]) that the pattern of Figure 8 was remarkably stable over the 150 years of the data and that the observed geometry could be modeled as Taylor–Couette flow of mode 3 (NH) or 4 (SH). The remarkable regularity and order three symmetry of the NH triskeles occurs despite the lack of cylindrical symmetry of the northern continents. The stronger intensity and larger size of features in the SH is linked to the presence of the annular currents. Following Kepler and Laplace (see above), it appears that the molecules in the atmosphere cannot perturb polar axis motion. The atmosphere behaves as a rotating cylinder undergoing stationary flow. The top of the cylinder is a free surface; the bottom is the surface of the solid crust and/or fluid ocean. This topography interacts strongly with the atmosphere.

6.2. Sun–Earth Distance and Phases

The seasonal change in sign agrees with changes in the Sun–Earth distance. The present values of the ellipticity of the Earth’s orbit and the obliquity of its rotation axis lead to changes of orbital velocity with summer in the southern hemisphere at perihelion and summer in the northern hemisphere at aphelion. As discussed in [14,15], forces acting on the polar axis do so through a torque ( m 1 , m 2 ), obeying the Liouville-Euler system of equations (from Laplace (1799) [63]). The influence of these torques is in quadrature with the length of day (lod). Changes in orbital velocity of our planet lead to a variation in the solar torque, which acts against the sense of rotation of the planet. At perihelion, orbital velocity is 30.29 km/s and at aphelion 29.29 km/s. The torque associated with the Sun changes as this velocity changes. Thus, annual variations in lod should be in opposition to changes in the Sun–Earth distance, and they are (Figure 9a, bottom). In the northern latitudes, both the Earth’s rotation and the winds go from west to east, the Earth’s rotation slows in summer and accelerates in winter, hence there are positive pressure anomalies in winter, negative in summer. In other words, the pressure anomalies are due to the variation in rotation velocity (tracked through lod changes and winds). The phase difference between the annual components of lod and SLP is constant at 32 ± 1.2 days (Figure 9b). Relative variations of amplitudes should be similar, and indeed they are (Figure 9b).

7. Sketch of a Mechanism

Because more than 70% of the total signal variance is captured by the trend of the atmospheric pressure, which varies by only 1 per mil since 1850, one can resort to the theory of turbulent flow in the case of infinitesimal perturbations (cf. [65]). The non-stationary v 1 ( r , t ) perturbation of the stationary solution v 0 ( r ) verifies the Navier–Stokes equation for an incompressible fluid:
v 1 t + ( v 0 ) v 1 + ( v 1 ) v 0 = p 1 ρ + ν . grad v 1 , div v 1 = 0 .
The general solution can be written as the sum of particular solutions where the dependence of v 1 on time is of the form e i ω t . There is still no theoretical basis for a mathematical solution of the stability of flow about finite dimension bodies plunged in a rotating fluid (e.g., [65,66,67]). Only experiments show that flow stability depends on the relative value of the Reynolds number (R), each type of flow having its critical Reynolds number ( R c ), beyond which instability sets in (cf. [68]). It is still difficult to solve turbulent flows other than in the case of concentric cylinders (e.g., [69,70,71,72,73,74,75]). Fortunately, this is useful in the case of atmospheric and oceanic flows on the spherical Earth, as is suggested (among other evidence) by the annular structures shown in Figure 5, Figure 6, Figure 7 and Figure 8.
Ref. [76] tackled the study of stability of a fluid comprised between two concentric rotating cylinders at high R. It was generalized for any value of R by [77] that showed that particular independent solutions were given by:
v 1 ( r , φ , z ) = f ( r ) * e i ( n φ + k z ω t )
with n being an integer (n=0 for axial symmetry), k the real wavelength of the instability, and ω an acceptable frequency solving the equations in plane z = constant with v 1 = 0 on cylinders with radii r = R 1 and r = R 2 (see also the paper on triskeles by [15]). The Reynolds number is evaluated by Ω 1 R 1 2 / ν (or Ω 2 R 2 2 / ν ). For fixed n and k, the ω become a discrete set of ω n ( k ) values. In the case of the Earth, v 1 = 0 at the solid surface and at the atmosphere free surface, and the only observed solutions are 24 h and 1 year. We therefore expect to find harmonics of these two values. The detailed study for two concentric cylinders can be found in [66]. The v 1 flow is stationary and consists of toric vortices (or Taylor vortices) regularly spaced along the cylinder generatrices. If the series of periods ω n ( k ) linked to the annual oscillation is 1, 1/2, 1/3, 1/4 years (see Section 3), then one expects to find experimentally perturbations of order n = 3 or 4, as was already the case for trends ([1]; Figure 8a, top). We find in the present study that indeed such is the case for the annual (Figure 8b, middle) and semi-annual (Figure 8c, bottom) oscillations.

8. Conclusions

In this paper, we focus on the sources of the trend and annual components of both global sea-level pressure (SLP) and variations in the Earth’s rotation (RP–coordinates of the rotation pole, and lod—length of day) and test the hypothesis that there might be a causal link between them. In a previous study [1], it was shown, using singular spectrum analysis (SSA), that the mean SLP contains, in addition to a weak trend, a dozen quasi-periodic or periodic components. These periods are characteristic of the space-time evolution of the Earth’s rotation and are found in many characteristic features of solar and terrestrial physics. Polar drift and the free Chandler wobble have been studied in earlier papers ([13,14,15]). With the use of iterative SSA (iSSA), we have extracted the components of SLP in each grid cell and as a function of time since 1850. The trend averages 1009 hPa and varies by only 0.7 hPa over the 170 year period with available data. The amplitudes of the annual component and its harmonics decrease from 93 hPa for the annual to 21 hPa for the third harmonic. In contrast, the components with pseudo-periods longer than a year range between 0.2 and 0.5 hPa. The trend is as large as the annual component (21 hPa) and could be part of the 90 years Gleissberg cycle.
We have further analyzed the components of polar motion RP (coordinates m 1 and m 2 ) and length of day (lod). The annual components of RP and SLP have a phase difference of 152 ± 2 days, which is constant over 70 years; these components are modulated in the same way, growing in amplitude between ∼1870 and ∼1960. The phase difference happens to be half of the Euler period of 306 days, for reasons we do not yet know. There is a phase lag of ∼40 years between RP and SLP. With the Euler–Liouville system of equations [14,63] in mind, we propose that there may be a causal link between the two.
We have mapped the first three iSSA components of global SLP that account for more than 85% of the total data variance. The trend pattern is stable, and the observed geometry can be modeled as Taylor–Couette flow of mode 3 (NH) or 4 (SH). The stronger intensity and larger size of features in the SH is linked to the presence of the annular mode (SAM). The annual component is characterized by a large negative anomaly extending over all of Eurasia in the NH summer (and the opposite in the NH winter) and three large positive anomalies over Australia and the southern tips of South America and South Africa in the SH spring (and the opposite in the SH autumn). The semi-annual component is characterized by three positive anomalies extending over Iran and the surrounding central Asia, the northwest Atlantic, and the northern Pacific in the NH spring and autumn (and the opposite in the NH summer and winter) and in the SH spring and autumn by a strong stable pattern consisting of three large negative anomalies forming a triskeles within the 40°–60° annulus formed by the southern oceans. A large positive anomaly centered over Antarctica, with its maximum actually displaced toward Australia, and a smaller one centered over Southern Africa complement the pattern. The pattern is opposite in the NH summer and winter.
The main features that appear in the maps underline the importance of the geographical patterns of continent–ocean boundaries and of the occurrence of (zonal) annular features (particularly at latitudes south of 20° S).
To the first order, large scale atmospheric motions can be modeled as rotating cylinders, stable in time and space, in the shape of 3 or 4 mode triskeles. A component of flow forced by the Sun–Earth distance modifies slightly polar rotation leading to seasonality of pressure anomalies. Because lod is physically linked to the solid Earth, geographical regions with strong topography are most affected by these variations.
The present values of the ellipticity of the Earth’s orbit and the obliquity of its rotation axis lead to changes in orbital velocity. Forces acting on the polar axis do so through the torque ( m 1 , m 2 ). The torque changes as orbital velocity changes. The influence of this torque is in quadrature with the length of day (lod). Annual variations in lod are in opposition to changes in the Sun–Earth distance (Figure 9a, bottom). The phase difference between the annual components of lod and SLP is constant at 32 ± 1.2 days. Relative variations of amplitudes should be similar, and indeed they are (Figure 9b).
We propose that the pressure and rotational variations in the annual to multi-decadal range evidenced in increasing number and detail are, in significant part, driven by commensurabilities of the giant Jovian planets forcing variations on the Sun and in some cases directly on the inclination of the Earth’s rotation axis and velocity. These variations, in turn, excite variations in a number of Earth parts (atmospheric pressure, sea-level, sea ice, and probably many more). They can be considered as a form of high-frequency Milankovitch “cycles”. We intend to develop this point in a forthcoming paper.

Author Contributions

V.C., J.-L.L.M., F.L. and D.G. contributed to conceptualization, formal analysis, interpretation, and writing. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Université de Paris, IPGP and the LGL-TPE de Lyon.

Data Availability Statement

The used data are freely available at the following address: Met Office Hadley Centre: https://www.metoffice.gov.uk/hadobs/hadslp2/data/download.html (accessed on 10 October 2020). IERS: https://www.iers.org/IERS/EN/DataProducts/EarthOrientationData/eop.html (accessed on 12 October 2020).

Acknowledgments

We thank the anonymous reviewers for their helpful comments on the manuscript.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Appendix A

Figure A1. Spatio temporal evolution of the annual component of the rotation pole: (a) Lissajou pattern of the m 1 vs. m 2 coordinates of the iSSA annual component of the rotation pole; (b) same as in (a), with the added dimension of time (1850–2020).
Figure A1. Spatio temporal evolution of the annual component of the rotation pole: (a) Lissajou pattern of the m 1 vs. m 2 coordinates of the iSSA annual component of the rotation pole; (b) same as in (a), with the added dimension of time (1850–2020).
Atmosphere 13 01907 g0a1aAtmosphere 13 01907 g0a2b

References

  1. Lopes, F.; Courtillot, V.; Le Mouël, J.L.; Gibert, D. Triskeles and Symmetries of Mean Global Sea-Level Pressure. Atmosphere 2022, 13, 1354. [Google Scholar] [CrossRef]
  2. Chandler, S.C. On the variation of latitude, I. Astron. J. 1891, 11, 59–61. [Google Scholar] [CrossRef]
  3. Chandler, S.C. On the variation of latitude, II. Astron. J. 1891, 11, 65–70. [Google Scholar] [CrossRef]
  4. Markowitz, W. Concurrent astronomical observations for studying continental drift, polar motion, and the rotation of the Earth. In Symposium-International Astronomical Union; Cambridge University Press: Cambridge, UK, 1968; Volume 32, pp. 25–32. [Google Scholar]
  5. Stoyko, A. Mouvement seculaire du pole et la variation des latitudes des stations du SIL. In Symposium-International Astronomical Union; Cambridge University Press: Cambridge, UK, 1968; Volume 32, pp. 52–56. [Google Scholar]
  6. Kirov, B.; Georgieva, K.; Javaraiah, J. 22-year periodicity in solar rotation, solar wind parameters and earth rotation. In Solar variability: From core to outer frontiers, Proceedings of the 10th European Solar Physics Meeting, Prague, Czech Republic, 9–14 September 2002; Wilson, A., Ed.; ESA SP-506; ESA Publications Division: Noordwijk, The Netherlands, 2002; Volume 1, pp. 149–152. ISBN 92-9092-816-6. [Google Scholar]
  7. Lambeck, K. The Earth’s Variable Rotation: Geophysical Causes and Consequences; Cambridge University Press: Cambridge, UK, 2005; ISBN 05-2167-330-5. [Google Scholar]
  8. Zotov, L.; Bizouard, C. On modulations of the Chandler wobble excitation. J. Geodyn. 2012, 62, 30–34. [Google Scholar] [CrossRef]
  9. Markowitz, W.; Guinot, B. Continental Drift, Secular Motion of the Pole, and Rotation of the Earth; Springer Science & Business Media: Dordrecht, Holland, 2013; Volume 32. [Google Scholar]
  10. Chao, B.F.; Chung, W.; Shih, Z.; Hsieh, Y. Earth’s rotation variations: A wavelet analysis. Terra Nova 2014, 26, 260–264. [Google Scholar] [CrossRef]
  11. Zotov, L.; Bizouard, C.; Shum, C.K. A possible interrelation between Earth rotation and climatic variability at decadal time-scale. Geod. Geodyn. 2016, 7, 216–222. [Google Scholar] [CrossRef] [Green Version]
  12. Lopes, F.; Le Mouël, J.L.; Gibert, D. The mantle rotation pole position. A solar component. Comptes Rendus Geosci. 2017, 349, 159–164. [Google Scholar] [CrossRef]
  13. Le Mouël, J.L.; Lopes, F.; Courtillot, V. Sea-Level Change at the Brest (France) Tide Gauge and the Markowitz Component of Earth’s Rotation. J. Coast. Res. 2021, 37, 683–690. [Google Scholar] [CrossRef]
  14. Lopes, F.; Le Mouël, J.L.; Courtillot, V.; Gibert, D. On the shoulders of Laplace. Phys. Earth Planet. Inter. 2021, 316, 106693. [Google Scholar] [CrossRef]
  15. Lopes, F.; Courtillot, V.; Le Mouël, J.L.; Gibert, D. On two formulations of polar motion and identification of its sources. arXiv 2022, arXiv:2204.11611v1. [Google Scholar] [CrossRef]
  16. Gleissberg, W. A long-periodic fluctuation of the sunspot numbers. Observatory 1939, 62, 158. [Google Scholar]
  17. Jose, P.D. Sun’s motion and sunspots. Astron. J. 1965, 70, 193–200. [Google Scholar] [CrossRef]
  18. Coles, W.A.; Rickett, B.J.; Rumsey, V.H.; Kaufman, J.J.; Turley, D.G.; Ananthakrishnan, S.A.J.W.; Sime, D.G. Solar cycle changes in the polar solar wind. Nature 1980, 286, 239–241. [Google Scholar] [CrossRef]
  19. Charvatova, I.; Strestik, J. Long-term variations in duration of solar cycles. Bull. Astron. Inst. Czechoslov. 1991, 42, 90–97. [Google Scholar]
  20. Scafetta, N. Empirical evidence for a celestial origin of the climate oscillations and its implications. J. Atmos. Sol.-Terr. Phys. 2010, 72, 951–970. [Google Scholar] [CrossRef] [Green Version]
  21. Le Mouël, J.L.; Lopes, F.; Courtillot, V. Identification of Gleissberg cycles and a rising trend in a 315-year-long series of sunspot numbers. Sol. Phys. 2017, 292, 1–9. [Google Scholar] [CrossRef]
  22. Usoskin, I.G. A history of solar activity over millennia. Living Rev. Sol. Phys. 2017, 14, 1–97. [Google Scholar] [CrossRef] [Green Version]
  23. Scafetta, N. Solar oscillations and the orbital invariant inequalities of the solar system. Sol. Phys. 2020, 295, 1–19. [Google Scholar] [CrossRef]
  24. Courtillot, V.; Lopes, F.; Le Mouël, J.L. On the prediction of solar cycles. Sol. Phys. 2021, 296, 1–23. [Google Scholar] [CrossRef]
  25. Scafetta, N. Reconstruction of the interannual to millennial scale patterns of the global surface temperature. Atmosphere 2021, 12, 147. [Google Scholar] [CrossRef]
  26. Wood, C.A.; Lovett, R.R. Rainfall, drought and the solar cycle. Nature 1974, 251, 594–596. [Google Scholar] [CrossRef]
  27. Mörth, H.T.; Schlamminger, L. Planetary motion, sunspots and climate. In Solar-Terrestrial Influences on Weather and Climate; Springer: Dordrecht, The Netherlands, 1979; pp. 193–207. [Google Scholar]
  28. Mörner, N.A. Planetary, solar, atmospheric, hydrospheric and endogene processes as origin of climatic changes on the Earth. In Climatic Changes on a Yearly to Millennial Basis; Springer: Dordrecht, The Netherlands, 1984; pp. 483–507. [Google Scholar]
  29. Schlesinger, M.E.; Ramankutty, N. An oscillation in the global climate system of period 65–70 years. Nature 1994, 367, 723–726. [Google Scholar] [CrossRef]
  30. Lau, K.M.; Weng, H. Climate signal detection using wavelet transform: How to maketime series sing. Bull. Am. Meteorol. Soc. 1995, 76, 2391–2402. [Google Scholar] [CrossRef]
  31. Courtillot, V.; Gallet, Y.; Le Mouël, J.L.; Fluteau, F.; Genevey, A. Are there connections between the Earth’s magnetic field and climate? Earth Planet. Sci. Lett. 2007, 253, 328–339. [Google Scholar] [CrossRef]
  32. Courtillot, V.; Le Mouël, J.L.; Kossobokov, V.; Gibert, D.; Lopes, F. Multi-decadal trends of global surface temperature: A broken line with alternating ˜30 years linear segments? Atmos. Clim. Sci. 2013, 3, 34080. [Google Scholar] [CrossRef] [Green Version]
  33. Le Mouël, J.L.; Lopes, F.; Courtillot, V. A solar signature in many climate indices. J. Geophys. Res. Atmos. 2019, 124, 2600–2619. [Google Scholar] [CrossRef]
  34. Scafetta, N.; Milani, F.; Bianchini, A. A 60-year cycle in the Meteorite fall frequency suggests a possible interplanetary dust forcing of the Earth’s climate driven by planetary oscillations. Geophys. Res. Lett. 2020, 47, e2020GL089954. [Google Scholar] [CrossRef]
  35. Connolly, R.; Soon, W.; Connolly, M.; Baliunas, S.; Berglund, J.; Butler, C.J.; Zhang, W. How much has the Sun influenced Northern Hemisphere temperature trends? An ongoing debate. Res. Astron. Astrophys. 2021, 21. [Google Scholar] [CrossRef]
  36. Guinot, B. Variation du pôle et de la vitesse de rotation de la Terre. In Traité de Géophysique Interne, Tome I: Sismologie et Pesanteur; Masson & Cie, Éditeurs: Paris, France, 1973. [Google Scholar]
  37. Ray, R.D.; Erofeeva, S.Y. Long-period tidal variations in the length of day. J. Geophys. Res. Solid Earth 2014, 119, 1498–1509. [Google Scholar] [CrossRef]
  38. Le Mouël, J.L.; Lopes, F.; Courtillot, V.; Gibert, D. On forcings of length of day changes: From 9-day to 18.6-year oscillations. Phys. Earth Planet. Inter. 2019, 292, 1–11. [Google Scholar] [CrossRef]
  39. Dumont, S.; Silveira, G.; Custódio, S.; Guéhenneux, Y. Response of Fogo volcano (Cape Verde) to lunisolar gravitational forces during the 2014–2015 eruption. Phys. Earth Planet. Inter. 2021, 312, 106659. [Google Scholar] [CrossRef]
  40. Petrosino, S.; Dumont, S. Tidal modulation of hydrothermal tremor: Examples from Ischia and Campi Flegrei volcanoes, Italy. Front. Earth Sci. 2022, 9, 775269. [Google Scholar] [CrossRef]
  41. Stephenson, F.R.; Morrison, L.V. Long-term changes in the rotation of the Earth: 700 BC to AD 1980. Philos. Trans. R. Soc. A 1984, 313, 47–70. [Google Scholar] [CrossRef]
  42. Hulot, G.; Le Huy, M.; Le Mouël, J.L. Influence of core flows on the decade variations of the polar motion. Geophys. Astrophys. Fluid Dyn. 1996, 82, 35–67. [Google Scholar] [CrossRef]
  43. Deng, S.; Liu, S.; Mo, X.; Jiang, L.; Bauer-Gottwein, P. Polar drift in the 1990s explained by terrestrial water storage changes. Geophys. Res. Lett. 2021, 48, e2020GL092114. [Google Scholar] [CrossRef]
  44. Gross, R.S.; Fukumori, I.; Menemenlis, D. Atmospheric and oceanic excitation of the Earth’s wobbles during 1980–2000. J. Geophys. Res. Solid Earth 2003, 108. [Google Scholar] [CrossRef] [Green Version]
  45. Bizouard, C.; Seoane, L. Atmospheric and oceanic forcing of the rapid polar motion. J. Geod. 2010, 84, 19–30. [Google Scholar] [CrossRef]
  46. Chen, J.; Wilson, C.R.; Kuang, W.; Chao, B.F. Interannual oscillations in earth rotation. J. Geophys. Res. Solid Earth 2019, 124, 13404–13414. [Google Scholar] [CrossRef]
  47. Landau, L.D.; Lifchitz, E.M. Mechanism: Course of Theoretical Physics; Mir Edition: Moscow, Russia, 1984. [Google Scholar]
  48. Wilson, C.R.; Haubrich, R.A. Meteorological excitation of the Earth’s wobble. Geophys. J. Int. 1976, 46, 707–743. [Google Scholar] [CrossRef] [Green Version]
  49. Spitaler, R. Die ursache der breitenschwankungen; Wien, K.-K., Ed.; Hof- und staatsdruckerei: Vienna, Austria, 1897. [Google Scholar]
  50. Spitaler, R. Die periodischen Luftmassenvershiebungen und ihr Einfluss auf dieLagenanderungen der Erdachse (Breitenschwankungen). Petermanns Mitteilungen Erganxungsband 1901, 29, 137. [Google Scholar]
  51. Jeffreys, H. Causes contributory to the Annual Variation of Latitude. (Plate 8). Mon. Not. R. Astron. Soc. 1916, 76, 499–525. [Google Scholar] [CrossRef] [Green Version]
  52. Munk, W.; Mohamed, M. Atmospheric excitation of the Earth’s wobble. Geophys. J. Int. 1961, 4, 339–358. [Google Scholar] [CrossRef] [Green Version]
  53. Allan, R.; Ansell, T. A new globally complete monthly historical gridded mean sea level pressure dataset (HadSLP2): 1850–2004. J. Clim. 2006, 19, 5816–5842. [Google Scholar] [CrossRef]
  54. Golyandina, N.; Zhigljavsky, A. Singular Spectrum Analysis for Time Series; Springer: Berlin, Germany, 2013; ISBN 978-3642349126. [Google Scholar]
  55. Lemmerling, P.; Van Huffel, S. Analysis of the structured total least squares problem for hankel/toeplitz matrices. Numer. Algorithms 2001, 27, 89–114. [Google Scholar] [CrossRef]
  56. Golub, G.H.; Reinsch, C. Singular value decomposition and least squares solutions. In Linear Algebra; Springer: Berlin/Heidelberg, Germany, 1971; pp. 134–151. [Google Scholar]
  57. Usoskin, I.G.; Solanki, S.K.; Schüssler, M.; Mursula, K.; Alanko, K. Millennium-scale sunspot number reconstruction: Evidence for an unusually active Sun since the 1940s. Phys. Rev. Let. 2017, 91, 211101. [Google Scholar] [CrossRef] [Green Version]
  58. Schwabe, H. Sonnenbeobachtungen im Jahre 1843. Von Herrn Hofrath Schwabe in Dessau. Astron. Nachrichten 1844, 21, 233. [Google Scholar]
  59. Adhikari, S.; Ivins, E.R. Climate-driven polar motion: 2003–2015. Sci. Adv. 2016, 2, e1501693. [Google Scholar] [CrossRef] [Green Version]
  60. Vondrák, J.; Ron, C.; Pesek, I.; Cepek, A. New global solution of Earth orientation parameters from optical astrometry in 1900–1990. Astron. Astrophys. 1995, 297, 899. [Google Scholar]
  61. Gibert, D.; Holschneider, M.; Le Mouël, J.L. Wavelet analysis of the Chandler wobble. J. Geophys. Res. Solid Earth 1998, 103, 27069–27089. [Google Scholar] [CrossRef]
  62. Gibert, D.; Le Mouël, J.L. Inversion of polar motion data: Chandler wobble, phase jumps, and geomagnetic jerks. J. Geophys. Res. Solid Earth 2008, 113. [Google Scholar] [CrossRef] [Green Version]
  63. Laplace, P.S. Traité de Mécanique Céleste; de l’Imprimerie de Crapelet: Paris, France, 1799. [Google Scholar]
  64. Lamb, H.H. Climate: Present, Past and Future. 1. Fundamentals and Climate Now; Methuen Publishing: London, UK, 1972. [Google Scholar]
  65. Landau, L.D.; Lifshitz, E.M. Fluid Mechanics: Courses of Theoretical Physics; Pergamon Press: Oxford, UK, 1987; Volume 6. [Google Scholar]
  66. Chandrasekhar, S. Hydrodynamic and Hydromagnetic Stability; Oxford University Press: Oxford, UK, 1961; ISBN 978-0486640716. [Google Scholar]
  67. Frisch, U. Turbulence: The Legacy of AN Kolmogorov; Cambridge University Press: Cambridge, UK, 1995; ISBN 0521457130. [Google Scholar]
  68. Landau, L.D. On the problem of turbulence. Dokl. Akad. Nauk USSR 1944, 44, 311. [Google Scholar]
  69. Schrauf, G. The first instability in spherical Taylor-Couette flow. J. Fluid Mech. 1986, 166, 287–303. [Google Scholar] [CrossRef] [Green Version]
  70. Mamun, C.K.; Tuckerman, L.S. Asymmetry and Hopf bifurcation in spherical Couette flow. Phys. Fluids 1995, 7, 80–91. [Google Scholar] [CrossRef] [Green Version]
  71. Nakabayashi, K.; Tsuchida, Y. Flow-history effect on higher modes in the spherical Couette system. J. Fluid Mech. 1995, 295, 43–60. [Google Scholar] [CrossRef]
  72. Hollerbach, R.; Junk, M.; Egbers, C. Non-axisymmetric instabilities in basic state spherical Couette flow. Fluid Dyn Res. 2006, 38, 257. [Google Scholar] [CrossRef]
  73. Mahloul, M.; Mahamdia, A.; Kristiawan, M. The spherical Taylor–Couette flow. Eur. J. Mech.-B/Fluids 2016, 59, 1–6. [Google Scholar] [CrossRef]
  74. Garcia, F.; Seilmayer, M.; Giesecke, A.; Stefani, F. Modulated rotating waves in the magnetised spherical Couette system. J. Nonlinear Sci. 2019, 29, 2735–2759. [Google Scholar] [CrossRef] [Green Version]
  75. Mannix, P.M.; Mestel, A.J. Bistability and hysteresis of axisymmetric thermalconvection between differentially rotating spheres. J. Fluid Mech. 2021, 911. [Google Scholar] [CrossRef]
  76. Rayleigh, L. On convection currents in a horizontal layer of fluid, when the higher temperature is on the under side. Lond. Edinb. Dublin Philos. Mag. J. Sci. 1916, 32, 529–546. [Google Scholar] [CrossRef]
  77. Taylor, G.I., VIII. Stability of a viscous liquid contained between two rotating cylinders. Philos. Trans. R. Soc. A 1923, 223, 289–343. [Google Scholar]
Figure 1. The iSSA annual component of rotation pole coordinates m 1 and m 2 from 1850 to the present.
Figure 1. The iSSA annual component of rotation pole coordinates m 1 and m 2 from 1850 to the present.
Atmosphere 13 01907 g001
Figure 3. The iSSA annual component of sea level pressure SLP from 1850 to the present.
Figure 3. The iSSA annual component of sea level pressure SLP from 1850 to the present.
Atmosphere 13 01907 g003
Figure 4. Annual envelopes and trends extracted from pole movement and SLP. (a) Envelopes of oscillations of iSSA annual components of polar motion m (blue curve, right scale) and global sea-level pressure SLP (black curve, left scale). (b) Trend of the pole movement (blue curve) and envelope of iSSA component of atmospheric pressure SLP (black curve).
Figure 4. Annual envelopes and trends extracted from pole movement and SLP. (a) Envelopes of oscillations of iSSA annual components of polar motion m (blue curve, right scale) and global sea-level pressure SLP (black curve, left scale). (b) Trend of the pole movement (blue curve) and envelope of iSSA component of atmospheric pressure SLP (black curve).
Atmosphere 13 01907 g004
Figure 5. Overlay of the triskeles patterns (SSA trend, top) and SSA annual oscillation of SLP (bottom) from Lopes et al. (2022) with maps from Lamb (1972) page 157, Figures 4–13. Left, southern hemisphere; right, southern hemisphere. Lamb’s maps are calculated for a time span of 40 years; our SSA results are for 170 years.
Figure 5. Overlay of the triskeles patterns (SSA trend, top) and SSA annual oscillation of SLP (bottom) from Lopes et al. (2022) with maps from Lamb (1972) page 157, Figures 4–13. Left, southern hemisphere; right, southern hemisphere. Lamb’s maps are calculated for a time span of 40 years; our SSA results are for 170 years.
Atmosphere 13 01907 g005
Figure 6. Polar stereographic projection of the northern hemisphere (left column) and southern hemisphere (right column) showing the mean (1850–2020) of the annual oscillation (iSSA component 2) for all (from top to bottom) springs, summers, autumns, and winters since 1850: (a) during springs; (b) during summers; (c) during autumns; (d) during winters.
Figure 6. Polar stereographic projection of the northern hemisphere (left column) and southern hemisphere (right column) showing the mean (1850–2020) of the annual oscillation (iSSA component 2) for all (from top to bottom) springs, summers, autumns, and winters since 1850: (a) during springs; (b) during summers; (c) during autumns; (d) during winters.
Atmosphere 13 01907 g006aAtmosphere 13 01907 g006b
Figure 7. Same as Figure 5 but for the semi-annual oscillation (SSA component 3): (a) during springs; (b) during summers; (c) during autumns; (d) during winters.
Figure 7. Same as Figure 5 but for the semi-annual oscillation (SSA component 3): (a) during springs; (b) during summers; (c) during autumns; (d) during winters.
Atmosphere 13 01907 g007aAtmosphere 13 01907 g007b
Figure 8. Polar views of northern hemisphere (left column) and southern hemisphere (right column) showing the trend (iSSA component 1 from Lopes et al. [1]), and maps of iSSA components 2 (annual) and 3 (semi-annual) for all springs since 1850.
Figure 8. Polar views of northern hemisphere (left column) and southern hemisphere (right column) showing the trend (iSSA component 1 from Lopes et al. [1]), and maps of iSSA components 2 (annual) and 3 (semi-annual) for all springs since 1850.
Atmosphere 13 01907 g008aAtmosphere 13 01907 g008b
Figure 9. Astronomical and geophysical annual oscillation components: (a) (top) lod data (blue) and their annual component (red); (middle) variation in the Sun–Earth distance (red); (bottom) the two curves above superimposed, showing phase opposition; (b) in green, the annual SSA component 2 of global sea level pressure, and in blue, the annual component of lod. The phase difference is constant at 32 ± 1.2 days.
Figure 9. Astronomical and geophysical annual oscillation components: (a) (top) lod data (blue) and their annual component (red); (middle) variation in the Sun–Earth distance (red); (bottom) the two curves above superimposed, showing phase opposition; (b) in green, the annual SSA component 2 of global sea level pressure, and in blue, the annual component of lod. The phase difference is constant at 32 ± 1.2 days.
Atmosphere 13 01907 g009aAtmosphere 13 01907 g009b
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Courtillot, V.; Le Mouël, J.-L.; Lopes, F.; Gibert, D. On the Nature and Origin of Atmospheric Annual and Semi-Annual Oscillations. Atmosphere 2022, 13, 1907. https://doi.org/10.3390/atmos13111907

AMA Style

Courtillot V, Le Mouël J-L, Lopes F, Gibert D. On the Nature and Origin of Atmospheric Annual and Semi-Annual Oscillations. Atmosphere. 2022; 13(11):1907. https://doi.org/10.3390/atmos13111907

Chicago/Turabian Style

Courtillot, Vincent, Jean-Louis Le Mouël, Fernando Lopes, and Dominique Gibert. 2022. "On the Nature and Origin of Atmospheric Annual and Semi-Annual Oscillations" Atmosphere 13, no. 11: 1907. https://doi.org/10.3390/atmos13111907

APA Style

Courtillot, V., Le Mouël, J. -L., Lopes, F., & Gibert, D. (2022). On the Nature and Origin of Atmospheric Annual and Semi-Annual Oscillations. Atmosphere, 13(11), 1907. https://doi.org/10.3390/atmos13111907

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop