Atlantic Oscillation indices in meridional distribution

The goal of this study was to define new oscillation indices of the Atlantic Ocean covering latitudes of both hemispheres, for the period from 1900 until 2014. Indices described sea level atmospheric pressure (SLP) oscillations between neighbouring cells of global atmospheric circulation (northern and southern Hadley, Ferrel and polar cells) and were calculated using two methods: conventional (based on SLP values in selected stations) and extreme values method (counted using maximum or minimum SLP from adequate areas). Twelve time series of new Atlantic Oscillation indices have been obtained—six conventional and six extreme, which have been presented in continuous form and through selection data by seasons: December–February (DJF), and June–August (JJA). The indices’ time series showed similarities in shared maxima and minima moments and long-lasting increasing tendency of five extreme and three conventional indices. There were strong correlations between oscillation indices describing variability in Hadley cells on both hemispheres (conventional 0.593–0.907, extreme 0.088–0.908). A strong correlation occurred between Southern Atlantic and Antarctic South Atlantic Extreme indices (0.668 to 0.979). Overall, stronger connections were found between extreme indices than the conventional ones, and most regularities were found in the Southern Hemisphere oscillations and wintertime indices (DJF for Northern and JJA for Southern Hemisphere). Power spectrum analysis showed major impact of 6- and 12-month periods, as well as less distinct significance of approximately 5.5-year-long interval and its multiples.


Introduction
In idealised Earth thermal circulation model to maintain stable conditions, it is necessary to retain radiative equilibrium in which energy absorbed by atmosphere, hydrosphere and lithosphere is counterbalanced by the energy radiated by them. Since solar energy influx on the globe differs significantly depending on latitudes, hemisphere, earth's surface type and time, the imbalance of heat distribution is moderated through air mass movements in the atmosphere and water currents in the world's ocean (Vallis 2006).
Global atmosphere circulation from the equator to the poles can be described through a simplified system of northern and southern circulation cells: Hadley, Ferrel and polar (Fig. 1). In the tropics, air heated up near the earth surface rises to the upper layers of the troposphere and is moved in the direction of poles and the subtropical region cools down. A portion of this air mass drops to the sea level and trade winds drive circulation back to the equator, closing up the Hadley cell. A share of descending air from the upper troposphere is transported further in the polar directions, where it is deflected to the right in the Northern Hemisphere and to the left in the Southern Hemisphere, creating westerly winds. Mainly due to impact of low pressure areas, in the temperate regions, air is lifted to the upper troposphere and moved in the direction of the equator, creating the Ferrel cell. Part of the air from the upper stream of Ferrel cell is distorted and becomes part of the polar cell in which it is transported further to the poles and cooled down. Eventually, it drops to the surface and is carried back to lower latitudes (Fig. 1).
Variability of air mass circulation in troposphere and adjoined upper level of hydrosphere may be observed through oscillation indices that indicate energy exchange interactions through air pressure systems and oceanic vortices. Phases of oscillation indices (positive, neutral or negative) correspond to specific weather occurrences.
Currently, fluctuations in oscillation indices' time series are analysed in different scales and localisations. These research were pioneered by Sir Gilbert Walker who in the 1930s defined North Atlantic Oscillation (NAO), a bipolar variability of sea level pressure (SLP) between two constant pressure centres: Islandic Low and Azores High (Walker and Bliss 1932). High NAO index values are a sign of relatively large pressure differences between two centres that induce westerly winds carrying wet air masses over Europe and consequently pushing cold arctic air to North America. Low NAO index means weakening of this mechanism causing continental weather in northern Europe.
Since the first half of the twentieth century, studies on the North Atlantic Oscillation (NAO) have been based on twopoint time series of atmospheric pressure fluctuations in the area of the subtropical Azores High and the Icelandic Low. The research method focused on pairing different measuring stations and extending the observation series (Rogers 1984, Hurrell 1995, Jones et al. 1997, Osborn 2006, Cornes et al. 2013). Thanks to increase in available data and more sophisticated methodology, NAO index time series have been extended, firstly to 1823 (Jones et al. 1997), then to 1692 (Cornes et al. 2013), up to 1500 , and with a larger degree of uncertainty for the whole duration of Holocene (Rimbu et al. 2003).
Two-point indicators have the advantage of a wide range of historical data availability, but they do not reflect dynamic changes in pressure action centres, so the currently accepted notion is that spatially defined indices (most often based on Empirical Orthogonal Functions (EOF)) provide a more precise description of NAO fluctuation (Jones et al. 1997;Li and Wang 2003). Other approaches to overcome this obstacle include the NAOI index introduced by Li and Wang (2003) using normalised averaged zonally SLP or Borström methodology (2012) determining NAO index with the highest from the Azores High and lowest from the Icelandic Low region SLP values. The advantage of this last solution was used in the submitted work. Nevertheless, the station-based indices were not abandoned, since the correlation coefficient between them and EOF indices was proven to be 0.92 (Hurrell et al. 2003).
In the last four decades of the twentieth century, a significant warming in extratropical areas in the Northern Hemisphere, including Europe, has been observed, coinciding with an unprecedented upward trend in the course of the NAO index (Hurrell 1995;Osborn 2006). Naturally, this increase in the NAO indicator has been linked to the global warming of the area, so it has strongly increased interest in its studies since the 1990s.
Since then, NAO has been considered a factor in the number of occurring anomalies, among others, of surface air temperature (SAT) (Hurrell and Van Loon 1997;Folland et al. 2008;Osborn 2011;Deser et al. 2017), SST (Hurrell 1995;Hurrell and Van Loon 1997;Czaja and Frankignoul 2001), precipitation (Hurrell 1995;Hurrell and Van Loon 1997;Jones et al. 1997;Hurrell et al. 2003;Folland et al. 2008;Deser et al. 2017), wind tracks and magnitude (Hurrell 1995, Hurrell and Van Loon 1997, Osborn 2006, Delworth et al. 2016, cloudiness (Folland et al. 2008), sea ice melting (Hurrell and Deser 2010;Delworth et al. 2016), ocean mixing (Hurrell andDeser 2010, Delworth et al. 2016) and even ecology and economy (Hurrell et al. 2003) of Europe and North America. Despite these efforts, the mechanisms controlling NAO fluctuations are not yet fully understood (Hurrell et al. 2003;Hurrell and Deser 2010) with a multitude of factors being considered, including increasingly anthropogenic impact (Hurrell 1995;Osborn 2006;Folland et al. 2008;Dong et al. 2011;Deser et al. 2017). Oscillation of NAO showed either quasi-variability in periods of 2 to 3 years (Hurrell et al. 2003) and more or less a decade (Hurrell et al. 2003, Li and Wang 2003. There is a consensus that North Atlantic Oscillation is a most important pattern influencing weather variability in this region (Osborn 2006;Wanner et al. 2001), but it has been suspected to be a regional emanation for a more complex oscillation in the Atmosphere-Ocean system, so it was analysed in comparison with other known oscillations (definitions in Table 1). With ENSO (Rogers 1984), NAO shared impact over western Atlantic, but had no noticeable communal variability. In other studies, North Atlantic Oscillation has been considered to be a subpart of Arctic Oscillation (Hurrell et al. 2003;Osborn 2006;Dong et al. 2011) especially since 1970s the relation between their indices has risen (Dong et al. 2011).
In the southern part of the Atlantic, Morioka et al. (2011) identified a new mechanism for the growth and decay of SST anomalies in the Southern Atlantic Subtropical Dipol (SASD) between the northeast and the southwest of the South Atlantic, showing growing dipole when the subtropical High South Atlantic migrates south and becomes stronger in late spring. It then causes a positive anomaly of the latent heat stream and reduces the thickness of the mixing layer in early summer. Also, in the southern part of the Atlantic, an oscillation in dominant Southern Atlantic Subtropical High (SASH) was used for the climatic description (Sun et al. 2017). The main study results revealed that SASH achieves the largest available surface ranges and pressure extremes during the solstice months (in June and December). In June, during the winter in the Southern Hemisphere, the SASH is the weakest, closest to the equator and located on the west side of the South Atlantic basin. In December, during the austral summer, it moves farthest towards the South Pole and reaches maximum pressure. This mode occurs in sync with La Niña phase in the Southern Oscillation (ENSO) as well as the positive phase of Southern Annular Mode (SAM) (Sun et al. 2017).
Through analysis in atmosphere fluctuation over other oceans, a multitude of different oscillations and teleconnections have been described, from which the most considered are summed up in Table 1. Continued advancement of the field flourished in using different then SLP variables to describe its lability such as sea surface temperature (SST), mass air movement tendencies and shifts in geopotential height.
It is apparent from Table 1 and referred literature that there is a disproportion between the number of defined oscillation indices on both hemispheres. NAO AMO, PNA and PDO describe only the Northern Hemisphere's dynamics; ENSO, MJO and IOD refer to oscillations in the vicinity of the equator, and only the Antarctic Oscillation represents variations on the southern half of the globe. South Atlantic Subtropical Dipole that describes SST between north-eastern and western parts of the South Atlantic is a fluctuation driven by ocean mixed layer thickness, not directly explained by atmospheric circulation over the ocean. South Subtropical High location and magnitude analysed by Sun et al. (2017) needs further investigation for it to be represented as a numerical descriptor. The lack of precisely defined oscillation indices of southern Atlantic is particularly contrasting with the well-examined northern part of this ocean. This is important since in spite of the simplified scheme from Fig. 1 atmospheric circulation in the Southern Hemisphere is not identical to that of the Northern, mainly due to an asymmetric distribution of solar energy, uneven thermohaline circulation and disproportion in quantity of land and water between the two.
The overview from Table 1 also shows diversity of the parameters used for calculation of particular indices as well as different methodologies used, which disallows direct comparison of their values, phases and coherence coefficients.
Therefore, there is a demand for oceanic oscillation indices on varied latitudes, calculated using the same variables. Identical methodology and data source promise comparability of homogenous results.
In response to above stipulations, the goal of this study was to calculate oscillation indices on latitudes covering both hemispheres of the Atlantic, consistent in both methodology and data source. The indices were calculated using two earlier formulas: first using standard methods introduced by Walker and Bliss (1932) and simplified by Rogers (1984) and formulated by Hurrell (1995), and the second using the highest SLP values from high pressure centres and lowest SLP values from low pressure centres from selected areas (Bostrӧm 2012).
The data source and its validation was also described. In the results section, new indices' time series are presented: their basic statistic, analysis of their long-lasting trends, extreme moments and power spectra. Then, the relations between particular oscillation indices were shown through correlation and coherence investigation. In the final part of the article, the discussion with previous works and conclusions from the whole analysis were summed up and further goals outlined.

Materials and methods
In this study, the fluctuations in oscillation indices on Atlantic Ocean were examined throughout 115 years, calculated with a formula (Hurrell 1995): Where, SLP high Sea level pressure on the high pressure station SLP high115 Mean sea level pressure on the high pressure station in the 115-year period SLP low Sea level pressure on the low pressure station  Barnston and Livezey 1987, Lin et al. 2015, Madden and Julian 1971, Mantua et al. 1997, Mantua and Hare 2002, Marshall 2003, Rogers and Van Loon 1982, Sae-Rim and Kwang-Yul 2015, Saji et al. 1999, Schlesinger and Ramankutty 1994, Thompson and Wallace 1998, Walker and Bliss 1932, Wallace and Gutzler 1981 SLP low115 Mean sea level pressure on the low pressure station in the 115-year period σStandard deviation These indices were named conventional and were calculated using SLP values from chosen stations located in the areas of constant pressure centres (alternating lows and highs) in the Atlantic (Fig. 2, Table 2). Therefore, this index shows differences between the standardised air pressure values between stations representing anticyclone and cyclone regions. The localisation of pressure centres overlaps the climatic zones along the meridian. The nomenclature of the indices was rooted in naming introduced by Walker and Bliss (1932).
Taking into account the dynamic variability of pressure centres, Bostrӧm (2012) proposed an extended formula (2) using the highest values of SLP in the high pressure centres and minimal SLP values from low pressure centres from general areas of constant highs and lows: Where, max(SLP high _ centre )Maximum SLP in the strongest high centres of high-pressure zone max(SLP high _ centre _ 115 )Mean maximum SLP in the strongest high centre of high-pressure zone in a 115-year period min(SLP low _ centre )Minimum SLP in the deepest low centre of low-pressure zone min(SLP low _ centre _ 115 )Mean minimum SLP in the deepest low centres of low-pressure zone in a 115-year period σStandard deviation These extreme indices, calculated with formula (2), were named identically as their conventional counterparts but with a word "Extreme" in the end of them, e.g. NAOE-North Atlantic Oscillation Extreme. Extreme SLP values were gathered from Atlantic areas divided in accordance with the climatic zones, in the meridional belt between 70°W and 20°E (Fig. 2, Table 2).
Since in previous works wintertime NAO indices were proved to be superior fluctuation describers (Hurrell 1995, Jones et al. 1997, Li and Wang 2003, Hurrell et al. 2003, Osborn 2006, three types of each time series have been presented: continuous oscillations, boreal wintertime/austral summertime indices (mean values from December to February (DJF)) and boreal summertime/austral wintertime indices (mean values from June to August (JJA)). Summertime oscillations are also valuable according to results presented by Hurrell et al. (2003) and Folland et al. (2008).
For index calculation, the SLP values were downloaded from the Twentieth Century Reanalysis (V2c) provided by National Oceanic and Atmospheric Administration (NOAA) and Cooperative Institute for Research in Environmental  1981). The longest times series of observed data were acquired from the measurement stations in Reykjavik, Ponta Delgada, Porto Alegre and Stanley, Falklands. Values from the reanalysis-modelled data-were investigated through quantitative analysis by comparison with observational datasets. Data validation was done by calculating dependency ratio, root mean square error (RMSE) and correlation coefficients (all of the values were statistically significant, tested by T test, with α = 0.05) between observational and modelled data (Table 3).
Additionally, new North Atlantic Oscillation Indices were compared with officially available NAO time series from the NOAA's Climate Prediction Centre website (NOAA 2005) for further confirmation of the selected methodology (Fig. 3). All of NAO indices show parallelism, convergent amplitude and similar fluctuation trends. The correlation coefficients between mean yearly NAO index from NOAA (calculated using EOF method) and new conventional and extreme NAO series were respectively 0.814 and 0.846, in boreal wintertime 0.930 and 0.835, and in boreal summertime 0.510 and 0.681. These results are in agreement with outcomes presented by Hurrell et al. (2003), confirming validity of Northern Hemisphere wintertime stationary indices by their 0.92 correlation with EOF results. Disparities between time series of NOAA NAO and new NAO indices seen on Fig. 3 are most noticeable in the years before 1970, so when available data for analysis was reduced. Also, NOAA time series has been normalised with standard deviation of each month which is observable since the orange plot is less wavering than others. Other discrepancies may be attributed to less compatible results depending on season between stationary and EOF indices, since NAO's impact on weather variability in this region during summer was shown to be less significant (Hurrell et al. 2003, Li andWang 2003).
Conducted validation has proven satisfactory and confirmed using chosen methodology for calculating new Atlantic Oscillation Indices.
Basic statistics were calculated for mean monthly values of new oscillation indices: mean and extreme values, standard deviation and 95% confidence interval.  Analysis of series characteristics was performed through comparison of their overall fluctuations, extreme moments and amplitudes as well as their individual and collective periodicity using fast Fourier transform (FFT) (Mathworks 2017). In reference to results from Li and Wang (2003) showing that decadal variations of Atlantic index are stronger using data from all months, power spectrum analysis was performed on continuous datasets, not seasonal ones. Due to intense variability of new indices' series, they were monthly averaged, and later normalised with moving means of different lengths: 1, 3, 5, 11 and 25 years long.
In FFT analysis, the longest period investigated was 19.2 years which accounts to one-sixth of the total time of 115 years, although it is recommended that the longest cycle should be repeated ten times within the analysed time. This means that periods longer than 11.5 years may not be evaluated to a satisfactory degree (MathWorks 2017).
To investigate relations between new indices in particular climatic zones, Pearson correlation (Łomnicki 2003) and coherence (Kowalik 1968) coefficients were calculated for monthly averaged sequences values as well as for all moving-mean normalised series.

Basic statistics
In the indices' time series, long-lasting trends were noticeable with differing amplitudes indicating range of fluctuations. Annual mean values of all indices were close to 0; the largest deviations from it were in extreme index series, maximally 0.0056 in NAOE index. Also, standard deviation values were larger in the extreme index series with maximum value of 1.082 for NAOE. Predictably, with two exceptions, both minimal and maximal values were more extreme for the extreme indices. Confidence intervals were narrower for conventional indices with exception of North Atlantic Tropical indices, and ranged from 0.025 for ANAO to 0.057 for both NAOE and SAOE (Table 4).

New Atlantic Oscillation index fluctuations
On Fig. 4, new Atlantic Oscillation indices' variability from 1900 until 2014 was presented. Figure 4 was divided into six rows of plots each representing, in meridional order, the time series of overlapping conventional and extreme indices in three columns showing separately their annual, December-February and June-August fluctuations. Shown plots are mean yearly index values normalised with 3-year moving mean. Additionally, in Table 5, general characteristics of overall fluctuations are presented: their approximate amplitude, range and most noticeably extreme positive and negative moments. Table 5 also shows which global atmospheric circulation cells are represented by particular indices, although it is important to remember this is only a vastly simplified concept.
Oscillation amplitudes were generally smaller in the conventional indices than in extreme ones apart from the North Atlantic Tropic Oscillation. Wintertime values for both hemispheres (DJF for the north and JJA for the south) showed more similarities to annual oscillations and higher values. A regularity is shown that annual oscillations in general fluctuating around 0 are a composite of positive wintertime and negative summertime values. In boreal summertime (JJA), ANAO/ ANAOE and NAO/NAOE variability plots were flattened in comparison with their wintertime and annual counterpart and lacked common maxima. In the south, ASAO/ASAOE and SAO/SAOE variability is also reduced and has fewer shared characteristics of the annual fluctuations; however, this tendency is fainter than that in the north.   'Further from the equator indices' plots seem to be less parallel but repeated extremes still occur: negative in 1920 for ANAO, ANAOE, NAO, NAOE, NATO, NATOE, SATO, SATOE, SAO, SAOE and ASAOE;in 1940 for NAO, NAOE, NATO, NATOE, SAO, SAOE and ASAOE;in 1955-1965 for NAO, NAOE, NATO NATOE, SATO andSATOE;positive in 1910 for SAO, SAOE and ASAOE;in 1995 for NAO, NAOE, SAO, SAOE and ASAOE, often more vividly announced in the winter season (DJF for the Northern Hemisphere, JJA for the Southern Hemisphere).
A less obvious close association exists between fluctuations of SAO, SAOE and ASAOE, in which annual plots all have similar extreme points and are overall parallel despite SAOE values ranging in a wider spectrum. This is mimicked in both JJA and DJF plots, with JJA showing more similarities to annual changes.
In general, in many indices' time series exists a longrunning increasing tendency, apparent in plots of the following: ANAOE, NAOE, SATOE, SAOE and ASAOE, so all but NATOE extreme indices. This could signify a meaningful discrepancy between two used methods-either extreme formula due to a wider range of analysed data is more responsive to an important, secular trend or it is somehow defective.
It is also worth noting that more parallel variabilities of indices are represented in the Southern Hemisphere which may suggest that Atmosphere-Ocean circulation in the southern Atlantic is more systematic.

Periodicity in time series of new Atlantic Oscillation indices
Investigating periodicity in annual new Atlantic Oscillation indices' time series was performed using the FFT method. In power spectrum plots of all sequences, by far, the largest extremes were assigned to 6-and 12-month periods (Table 6). This indicates that most meaningful variability of SLP values on the Atlantic took place in annual and semi-annual cycles. In 70% of cases, yearly maxima overweighed 6-month periods by at least one order of magnitude, and its impact was most significant in extreme oscillation indices (Table 6).

Relationships between new Atlantic Oscillation indices
Connections between newly calculated oscillation indices were examined through correlation and coherence analyses on its monthly average values, as well as its moving-mean normalised sequences.
In Tables 8 and 9, Pearson's correlation coefficients for yearly mean values are presented for conventional and extreme indices, respectively, annually and separately for DJF and JJA months. Most significant relations were further investigated using power spectrum comparison and coherence analysis.
Tropical indices NATO-SATO and NATOE-SATOE describe interconnection between northern and southern Hadley cells and their relations were further examined on different levels of data normalisation. Calculated correlation values of annual fluctuations ranged from 0.593 to 0.907 for conventional indices and from 0.088 to 0.908 for extreme indices, with high correlations also in DJF and JJA period analyses with    (Table 9). Their power spectra both had maxima in 0.5-, 1-and 5.5-year periods (Figs. 5 and 6). Coherence above significance level of 0.4 (red line on right side of Fig. 5) between NATO and SATO occurred over ten times, occasionally drawing near 1 which indicates correlation of the signals and similarities in its determination. Between NATOE and SATOE, significant coherence happened more rarely and was weaker than per conventional indices (right side of Fig. 6).
For indices describing corresponding Ferrel and polar cells on both hemispheres (NAO-SAO, NAOE-SAOE, ANAO-ASAO, ANAOE-ASAOE), correlation coefficients for differently smoothened data series, annually and separately for DJF and JJA periods, did not show significant relations.
Very strong correlations, from 0.668 to 0.978, were noted between extreme indices of South Atlantic Oscillation (SAOE) and Antarctic South Atlantic Oscillation ASAOE, for both analysed seasons and continuous data. This is a proof of strong connection between circulation in austral Ferrel and polar cells. Their power spectra had shared maxima in 0.5-, 1-, 5.5-and 10.5-year periods. Coherence above 0.4 level occurred over ten times, occasionally getting close to a value of 1 suggesting deterministic connection of signals (Fig. 7).

Possible data errors and uncertainties
Firstly, possible uncertainties may occur through data acquisition process. Reanalysed values are always only modelled approximations and especially in secluded areas and in distant past, where there was no sufficient amount of observational data for validating the model, its values may be untrue. Therefore, results from the beginning of the twentieth century and from polar regions should be treated with special caution. This particularly concerns ASAO index series, since, due to lack of stations further south, the southernmost measurement point used for calculating conventional index, Vernadsky Antarctic station, is located in the general area of Subantarctic Low not Antarctic High (Fig. 2).
Chosen index calculation methods may also cause data errors; partly because of imperfection of available formulas, two solutions were used in this study. Conventional method is limited in its undynamic data source localisation, whereas extreme formula has not been well tested and validated in the scientific community. Visible in the results section are discrepancies between extreme and conventional indices, such as general increasing tendency in extreme sequences not noticeable in conventional plots, which may indicate inferiority of one of the methods. Data smoothing applied for sake of plot readability could have hidden valuable information.
Additionally, sequence analysis tools are also imperfect. Since only 115 years of validated data was available, FFT power spectrum examination can appropriately check only periods shorter than 11.5 years; however, it is possible that a longer cycle could be especially significant.
Correlation and coherence coefficients cannot detect all forms of relations between data series.

What the newly obtained results reveal: discussion
The correlation coefficients of NAO, NAOE, NAO DJF and NAOE DJF indices with NAO NOAA index were 0.814 and 0.846 for annual indices and 0.930 and 0.835 for winter ones; therefore, we concluded that the obtained time series represent real processes in the Atlantic Ocean climate system (Figs. 2 and 4). This positive validation confirms the correctness of chosen methodology of determining extreme and conventional indices.
The variability of new indices' oscillation showed some resemblance to previously reported data. Rogers (1984), Hurrell (1995), Hurrell and Van Loon (1997) and Li and Wang (2003) reported NAO index values decreasing from around the 1920s until the 1960s to the 1970s. This trend is distinguishable in new Atlantic index time series (Fig. 4), although a more general downward tendency from the beginning to more or less the middle of twentieth century is shown on NATO, NATOE, SATO, SATOE, SAOE and ASAOE plots. A long upward trend from the 1960s to the 1980s in NAO was reported (Hurrell 1995, and Hurrell and Van Loon 1997, Jones et al. 1997, Li and Wang 2003, which is also true for variability of new: NAO, NAOE, NATO, NATOE, SATO and SATOE, SAOE and ASAOE (Fig. 4). Increasing NAO trend in the second half of twentieth century was interrupted in the winter of 1995 (Jones et al. 1997) and particularly low values of NAO have since been recorded in the winter of 2009/2010 (Osborn 2011), but a more general increase is alleged in the near future (Folland et al. 2008) due to its possible relation to rising greenhouse gases present in the Earth's atmosphere (Dong et al.  . This presumption may be confirmed by escalating values of a majority of new Atlantic Oscillation indices (Fig. 4).
A power spectrum analysis has proved strong significance of 6-month and yearly periods for all new indices' fluctuations, as well as showed repeated importance of 5.5-year-long period multiples: The interannual and decadal variabilities in atmospheric oscillations have been identified for NAO: Hurrell and Van Loon (1997) and Hurrell et al. (2003) (2-4-, and 6-10-year periods), also recognised by Luterbacher et al. (2001), and is believed to describe South Atlantic Subtropical High position and magnitude (Sun et al. 2017). Some of the new Atlantic Oscillation indices share close to decadal variability with period of 10.5 years, possibly connected to previously shown decadal regularities in the results of Hurrell (1995) and Hurrell and Van Loon (1997). Rogers (1984) found that covariance between winter NAO and ENSO indices has been strongest in periods between 5 and 6 years in accordance with the 5.5-year cycle in a majority of new indices but the nature of this parallel is unknown. The 25-30-year cyclicity of SST thermal anomalies was identified in the western subtropical Southern Ocean (Wainer and Venegas 2002). The obtained 5.5-year cycles and their multiples of 10.5 and 16.4 years of conventional and extreme indices of the Southern Hemisphere in the Atlantic sector refer to this two decades. Therefore, the undertaken works have the character of "kinematic descriptions of features" of one of the main components of the South Atlantic oscillations through time series of SATO, SATOE, SAO, SAOE, ASAO and ASAOE.
The determined correlations between the new oscillation indices showed a significant and high correlation in tropical and the combination of tropical and temperate zones in winter seasons. The relationship between NATO and SATO conventional indices is clearly shaped at 0.774 correlation. The connections between the NAO DJF and NATO DJF showed a correlation of 0.736 and between SAO JJA and SATO JJA of 0.522 (Tab. 8). This means that tropical zones, located almost symmetrically on both sides of the equator, appear to be strongly associated with Hadley cells; however, indices of moderate adjacent zones (Ferrel cells) correlate only during winter periods of atmospheric activity.
The high correlations between the ASAOE and SAOE indices suggest regularities in fluctuations of the atmosphere in the southern polar cell. During the winter in the Southern Hemisphere, the Subtropical South Atlantic High is displaced towards the pole in the La Niña years with coinciding positive phase of the Southern Annual Mode (SAM/AAO) (Sun et al. 2017). This shows that during the austral winter the atmosphere variability in the Southern Hemisphere is interconnected not only in directly neighbouring zones but possibly in the entire Southern Hemisphere. Further analysis into new Atlantic Oscillation indices could bring more information about these consistencies.
High correlations were confirmed by significant coherence in adjacent zones.
New conventional and extreme indices have been calculated for the area of the whole ocean and cautiously may be treated as proxies for oscillations in particular cells of global atmospheric circulation (northern and southern Hadley, Ferrel and polar cells).

Conclusions
Twelve new Atlantic Oscillation indices were established: (conventional) ANAO, NAO, NATO, SATO SAO, ASAO; (extreme) ANAOE, NAOE, NATOE, SATOE, SAOE and ASAOE. New indices' fluctuations had coherent tendencies, which in some respect are in accordance with NAO fluctuations recognised in previous studies, with the most apparent and common upward trend from the middle of the twentieth century.
An overall close proximity between NATO-SATO, NATOE-SATOE and SAOE-ASAOE has been detected, giving evidence for close relations between fluctuations in two Hadley cells, and southern polar and Hadley cells. This is linked to a general conclusion that more apparent regularities exist in southern indices suggesting uniform fluctuations in the austral part of the Atlantic.
More positive index values, and accentuated general tendencies shown in annual variations, are attributed to wintertime index variations, DJF for northern and JJA for Southern Hemisphere. This is a confirmation of previous works suggesting superiority of wintertime-based NAO index (Hurrell 1995, Jones et al. 1997, Li and Wang 2003, Hurrell et al. 2003, Osborn 2006; however, it is the first time this regularity has been shown for oscillation across all latitudes. FFT power spectrum analysis showed predominant impact of 0.5-and 1-year periods on all indices' variability, confirming results of Li and Wang (2003). Moreover, maxima around 5.5-year period and its multiplicities occurred, suggesting significance of this cycle for Atlantic Oscillation.
Correlation analysis showed stronger relations between extreme indices than conventional ones possibly indicating the superiority of this method. Strong correlations adjoining all Southern Hemisphere and extreme JJA indices imply localisation and approach that provide data with most recognised regularities.
Strong correlation and promising coherence between tropical indices NATO-SATO and NATOE-SATOE suggested symmetry in dynamics of Hadley cells on both hemispheres and a very strong correlation was visible between extreme indices of SAOE and ASAOE, describing southern Ferrel and polar cells, proving their cohesion.
In further studies, Empirical Orthogonal Function methodology should be applied to examine Atlantic Oscillation variability from even more dynamic perspective. Weather phenomena connected to specific phases of NAO have been thoroughly examined, and the same scrutiny is necessary to understand significance of new Atlantic Oscillation indices, especially in the southern part of the ocean. Moreover, past search for oscillations is mainly focused on the troposphere and its interaction with oceans. Increased interest in this approach and focus on the North Atlantic resulted in some forgetfulness of the direction initiated by Rossby and Willett (1948), who proposed a deterministic approach taking into account covariation in the troposphere and the stratosphere. With these in mind, further investigation may bring desirable information about Earth's climate system and possibilities of its real forecasting.