Interannual Variability and Trends in Sea Surface Temperature, Lower and Middle Atmosphere Temperature at Different Latitudes for 1980–2019

: The inﬂuence of sea-surface temperature (SST) on the lower troposphere and lower stratosphere temperature in the tropical, middle, and polar latitudes is studied for 1980–2019 based on the MERRA2, ERA5, and Met Ofﬁce reanalysis data, and numerical modeling with a chemistry-climate model (CCM) of the lower and middle atmosphere. The variability of SST is analyzed according to Met Ofﬁce and ERA5 data, while the variability of atmospheric temperature is investigated according to MERRA2 and ERA5 data. Analysis of sea surface temperature trends based on reanalysis data revealed that a signiﬁcant positive SST trend of about 0.1 degrees per decade is observed over the globe. In the middle latitudes of the Northern Hemisphere, the trend (about 0.2 degrees per decade) is 2 times higher than the global average, and 5 times higher than in the Southern Hemisphere (about 0.04 degrees per decade). At polar latitudes, opposite SST trends are observed in the Arctic (posi-tive) and Antarctic (negative). The impact of the El Niño Southern Oscillation phenomenon on the temperature of the lower and middle atmosphere in the middle and polar latitudes of the Northern and Southern Hemispheres is discussed. To assess the relative inﬂuence of SST, CO 2 , and other greenhouse gases’ variability on the temperature of the lower troposphere and lower stratosphere, numerical calculations with a CCM were performed for several scenarios of accounting for the SST and carbon dioxide variability. The results of numerical experiments with a CCM demonstrated that the inﬂuence of SST prevails in the troposphere, while for the stratosphere, an increase in the CO 2 content plays the most important role.


Introduction
Weather and climatic processes are largely determined by the interaction of the atmosphere and the ocean. In particular, simultaneous long-term changes in the temperature of the sea surface and the atmosphere can influence each other, changing the radiation regime and the dynamics of the atmosphere and the sea. In addition to long-term variability on a scale of several decades, short-term changes in sea surface temperature on a scale of several years can also affect atmospheric parameters. One of the most well-known short-term changes affecting both marine and atmospheric characteristics is the El Niño-Southern Oscillation (ENSO), which has El Niño and La Niño phases. Warming of the sea surface (El Niño phase) as well as cooling of the sea surface (La Niño phase) in the equatorial Pacific Ocean not only affects regional weather and climate changes, but also controls the general circulation of the atmosphere, which determines changes in air temperature and trace atmospheric gases at different altitudes in extratropical latitudes [1,2].
The atmospheric component of El Niño is the Southern Oscillation [1][2][3], in which there are fluctuations in atmospheric pressure between the eastern and western Pacific.
database [26]. The acquisition of SST data has been going on since the 19th century [27,28]. Due to the limitations associated with the ENSO observational data, it was difficult to determine the conditions for their occurrence [29,30]. By means of modeling with the Geophysical Fluid Dynamics Laboratory Climate Model version 2.1 (GFDL CM2.1), the amplitude corresponding to ENSO was determined [31]. However, many models did not betray the significance of the effect of CO 2 on the ENSO amplitude. Studies using Coupled Model Intercomparison Project phase 5 (CMIP5) analysis also did not show a significant relationship [32]. Therefore, this interaction is based on the physics of equilibrium changes in ENSO due to an increase in CO 2 , in the absence of transient effects [19].
Studies of temperature trends at different heights of the atmosphere are a topical subject of many recent scientific studies [33][34][35][36]. Studying the features of these trends is an important component of climate monitoring [37,38]. Trends in tropospheric warming and stratospheric cooling during the 20-21 centuries were studied on the basis of measurements of radiosondes and microwave sounding of the atmosphere [37,39,40]. The results of sounding and satellite observations were analyzed using climate models [41][42][43]. Simultaneous stratospheric temperature trends and changes in chemical processes were studied by comparing satellite observations and chemistry-climate modeling [34,44,45]. In a number of works, special attention was paid to the study of the features of warming in the tropics [46,47].
Studies of the vertical features of temperature trends revealed warming of the troposphere and cooling of the stratosphere [48][49][50]. In [34], the cooling of the stratosphere and its increase with height was associated with an increase in the concentration of greenhouse gases. It was found that the main factor of cooling in the lower stratosphere is a decrease in the ozone content, while the main role in the middle stratosphere is played by an increase in the content of greenhouse gases [33,34]. Quantitative estimates reveal a decrease in stratospheric temperature by 1-3 K over the past 40 years [33]. In this case, the lower stratosphere is cooled by −0.25 K per decade, while the average stratosphere is up to −0.5 K per decade [33]. At the same time, the cooling of the stratosphere is more significant in the Northern Hemisphere than in the Southern Hemisphere [33].
The results of chemistry-climate modeling of tropospheric and stratospheric temperature variability based on the given scenarios of changes in the atmospheric content of greenhouse gases and ozone-depleting substances (ODS), as well as taking into account other factors, such as volcanic aerosol and solar activity, are generally in good quantitative agreement with satellite measurements [43]. These results help to clarify many of the problems in understanding recent trends in tropospheric and stratospheric temperature. Nevertheless, some issues of tropospheric and stratospheric temperature variability, in particular latitudinal features of trends, are still not fully understood [34].
In this work, simultaneous trends and short-term variability of the sea surface and atmospheric temperature are estimated for the different regions and altitudes based on the reanalysis data and numerical experiments with a global chemistry-climate model for the last forty years . Special attention is paid to comparing the effect of sea surface temperature variability and greenhouse gases on tropospheric and stratospheric temperature changes based on calculations with the chemistry-climate model (CCM) for several scenarios of accounting for these factors. This work is a continuation of previous works for a shorter period of time [51][52][53][54]. Section 2 describes the methodology, data, and chemistry-climate model used. The results of assessing trends and short-term variability in sea surface and atmospheric temperature at different heights are presented in Section 3. Section 4 presents the results of numerical experiments with a global chemistry-climate model, in which the relative influence of sea-surface temperature variability and greenhouse gas content on tropospheric and stratospheric temperatures at different latitudes is estimated. Section 5 contains discussion and Section 6 is conclusion.

Methodology, Data and Model Description
The reanalysis MERRA2 [55], ERA5 [56] and Met Office [57] data were used to study the interannual variability of the sea surface temperature (SST) and the temperature of the lower troposphere and lower stratosphere. Sea surface temperature variability is analyzed using Met Office and ERA5 data, while atmospheric temperature variability is analyzed using MERRA2 and ERA5 data. The reanalysis data were interpolated into the Institute of Numerical Mathematics of the Russian Academy of Sciences (INM RAS)-Russian State Hydrometeorological University (RSHU) chemistry-climate model [58] grid (72 grids in longitude from −185 to 175 degrees with a step of 5 degrees, 45 grids in latitude from −88 to 88 degrees with a step of 4 degrees, 31 pressure levels from the surface to the mesopause with a variable step). The reanalysis data were averaged for each month, as well as for each year from 1980 to 2019.
Produced by the European Centre for Medium-Range Weather Forecasts (ECMWF) ERA Interim is a global atmospheric reanalysis that is available from 1 January 1980 to 31 December 2019 [56]. It provides hourly estimates of a large number of atmospheric, land, and oceanic climate variables. The data assimilation system used to produce ERA5 includes a 4-dimensional variational analysis (4D-Var) with a 12-h analysis window. The spatial resolution of the dataset is approximately 80 km (T255 spectral) on 60 levels in the vertical from the surface up to 0.1 hPa. ERA Interim has been superseded by the ERA5 reanalysis, which covers the Earth on a 30 km grid and resolve the atmosphere using 137 levels from the surface up to a height of 80 km [56]. Produced by the National Aeronautics and Space Administration (NASA) the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA2) provides atmospheric data from 1980 up to today [55]. It uses GEOS model with an assimilation system to produce long-term global reanalysis based on ground and space-based observations. Produced by the Met Office Hadley Centre, sea ice and sea surface temperature datasets are based on NOAA and METOP observations. Met Office sea surface temperatures are reconstructed using a two stage reduced-space optimal interpolation procedure, followed by superposition of quality-improved gridded observations onto the reconstructions to restore local detail [57].
To calculate trends in SST and atmospheric parameters, a linear regression equation and multiple regression analysis following [33], with an estimate of the trend significance, were used. The multivariate regression model includes a linear trend term and terms accounting solar cycle, ENSO, and stratospheric aerosol loading. Solar variability was presented by year mean sunspot number [59], ENSO by Met Office SST at tropical part of Pacific Ocean, and stratospheric aerosol by particle surface area [60]. For the analysis, the calculated linear trend coefficient (degrees per decade for sea surface and atmospheric temperatures) was used. For all calculated trends, their statistical significance was estimated by Fisher's criterion [61] by comparing the coefficient of the linear trend with the spread of data relative to the regression line. The probability density function (PDF) and coefficient of determination (R 2 ) were evaluated. In order to estimate linear regression significance, Fisher's criterion analysis of variance (ANOVA) was used for the sum of squares (SS). For estimation of the probability density function (PDF) that regression is significant, the cutoff value F = (N − 2)SS E /SS R , where N is number of points, SS R -sum of squares for regression, and SS E -sum of squares for residuals, with N-2 degrees of freedom, is used to calculate probability for F distribution [61]. The significance level for calculated F-value is assessed as difference between 1 and PDF. The coefficient of determination, which provides the part of total variance explained by the regression, is estimated as R 2 = 1 − SS R /SS Tot . Trends and short-term variability were analyzed for the lower troposphere (pressure levels from 1000 to 600 mb) and the lower stratosphere (pressure levels from 200 to 30 mb) of tropical, middle, polar latitudes, and on average over the globe.
To assess the role of SST and CO 2 affecting the interannual variability of the lower troposphere and low stratosphere temperature, numerical experiments with the INM RAS-RSHU CCM [58,62,63] were performed. The model has a spatial resolution of 4 × 5 degrees latitude/longitude and covers the altitude range from the Earth's surface up to the mesopause (around 90 km) with 39 sigma-pressure levels with variable pressure spacing [58,62,63]. In order to take into account an interaction between chemical and dynamical processes in the low and middle atmosphere under the frame of a CCM, a INM RAS General Circulation Model (GCM) is interactively coupled with a RSHU global chemistrytransport model (CTM), which was designed based on the three-dimensional extension of the SUNY-SPB two-dimensional model [58,62,63]. In addition to temperature and dynamical variables, the CCM calculates variability of 74 chemically active gases in the lower and middle atmospheres, which are interacting in 174 gas-phase and heterogeneous chemical reactions and 51 photodissociation processes, including oxygen, nitrogen, hydrogen, chlorine, bromine, carbon, and sulfur cycles [58,64], and takes into account gas concentrations impact on the temperature. The model considers the processes of the formation and evolution of polar stratospheric clouds (PSC) based on the stratospheric sulfate aerosol distribution [58]. Polar vortex dynamics are resolved by the CCM based on dynamical core of INM RAS GCM [58,64].
The influence of sea surface temperature on atmospheric characteristics in INM RAS-RSHU CCM is taken into account by setting a lower boundary condition in the air temperature evolution equation. For those CCM grids that have sea surface at low boundary, the temperature at the lower boundary is set by temporal interpolation of the SST monthly average reanalysis data. In addition, the area of ice coverage of the sea is also taken into account, which is also determined from the reanalysis data. To assess the influence of SST variability on the interannual variability of atmospheric characteristics in the background of an increase in the content of greenhouse gases in the atmosphere, calculations using the CCM were performed for the period 1979-2019 according to several scenarios. In the first (BASE) scenario, the variability of both SST and carbon dioxide content was taken into account, in the second (CO2FIX), atmospheric carbon dioxide content was fixed at the 1979 level, in the third (SSTFIX), SST was fixed at the 1979 level, taking into account seasonal variability, and finally, in the fourth scenario (FIX), the values of carbon dioxide and SST were fixed at the 1979 level. In all four scenarios, interannual variability of greenhouse gases (Climate System Scenario Tables [65]) and ozone-depleting substances (data of WMO-2018 [66]) was set. Results of model estimation were interpolated to the 31 vertical pressure levels (14 levels in the troposphere and 17 levels in the stratosphere and mesosphere). Figure 1 presents the variability of the sea surface temperature (SST) averaged over the entire globe (upper plot) and for the tropics (lower plot) according to the reanalyzes of Met Office and ERA5. According to all data, during the period from 1980 to 2019, both in the tropics and across the globe, there is a significant (1% significance level) positive SST trend. In the tropics, both datasets show similar trends of about 0.1 degrees per decade. At the same time, a slightly larger trend is obtained according to ERA5 data. A significant short-term increase in SST can be noted in 1982-1983, 1986-1987, 1997-1998, 2009-2010, and 2014-2016 for variability with a period of several years in the tropical zone, corresponding to the El Niño phase of the Southern Oscillation and decreased SST values 1984SST values -1985SST values , 1989SST values -1990SST values , 1996SST values , 1999SST values -2000SST values , 2008SST values , and 2010SST values -2011 to the La Niña phase of the Southern Oscillation. In accordance with the coefficient of determination, about half of the SST variance in the tropics is explained by a linear trend, and the second half is described by short-term variability. Globally, the estimates of trends differ almost twofold: The ERA5 data reveal a trend of about 0.1 degrees per decade, and the Met Office data reveal that of about 0.05 degrees per decade (Figure 1, top). All trends are significant with a 1% significance level. The short-term variability of the SST is also quite different between different reanalysis datasets. At the same time, for all reanalysis data on a global scale, periods of increase and decrease in SST are repeated corresponding to the most powerful phases of El Niño (1982-1983, 1986-1987, 1997-1998, 2009-2010, and 2014-2016) and La Niño (1984-1985, 1989-1990, 1996, 1999-2000, 2008, and 2010-2011). The largest amplitudes of SST fluctuations on a global scale are shown by the Met Office data, while in some years, the phases of the minima and maxima of Met Office and ERA5 do not coincide. Accordingly, about a half of the global mean SST variance in the Met Office data is described by short-period variability (about 51%), while for the ERA5 data, on the contrary, about 80% of the variance explained by a linear trend. In general, a comparison of the interannual variability of SST on a global scale and in the tropics indicates that the phases of the Southern Oscillation, which have a tropical nature, have a significant impact on SST changes on a global scale. Figure 2 presents the interannual variability of the SST in the middle latitudes of the Southern (bottom) and Northern (top) Hemispheres. In the Northern Hemisphere, SST trends for the period from 1980 to 2019, according to Met Office and ERA5 data (about 0.2 degrees per decade), are twice as large as in tropical latitudes and globally. In the middle latitudes of the Southern Hemisphere, on the contrary, the SST trend is significantly less than on a global scale and in the tropics. In the Southern Hemisphere, SST trends are 5 times less than in the Northern Hemisphere (0.025-0.035 degrees per decade), which, on the one hand, indicates a significant difference in ocean processes in different hemispheres, and on the other hand, shows that in addition to the anthropogenic factor of global warming, which potentially should have the same effect in different hemispheres, the natural variability, which is different in different hemispheres, probably plays an important role, and, thirdly, creates the potential for different effects of SST variability on atmospheric parameters in different hemispheres. The short-term variability of SST at midlatitudes has a rather pronounced seasonal variability, which veils interannual fluctuations. Nevertheless, against the background of seasonal variability, the influence of the Southern Oscillation is noticeable in some years. An analysis of the coefficients of determination detects that, despite the fact that SST trends are significant in both hemispheres with a 1% significance level, in the Northern Hemisphere, the trend describes more than 80% of the variance, and in the Southern Hemisphere, it explains only 17% for the ERA5 data and 55% for the Met Office data.   At the same time, after 2016, ERA5 data demonstrate a significant increase in SST in the polar zone of the Southern Hemisphere. The short-term variability in the polar regions demonstrates a decrease in the minimum SST in the Southern Hemisphere and an increase in the maximum SST in the Northern Hemisphere. This can affect dynamic conditions affecting atmospheric parameters in the lower and middle atmosphere.

Simultaneous Change of Sea Surface and Air Temperature Trend and Short-Term Variability Based on Reanalysis Data
The polar regions are characterized by large differences between ERA5 and Met Office data. Trends according to Met Office have low values in both the Northern and Southern Hemispheres, although they are significant at a 1% significance level, but with a low coefficient of determination. At the same time, short-term variability explains more than 80% of the variance in the Northern Hemisphere and more than 70% of the variance in the Southern Hemisphere. At the same time, the trend estimated from the ERA5 data describes almost 90% of the variance in the Northern Hemisphere, but about a half in the Southern Hemisphere. Wherein, similar to middle latitudes, the ERA5 data have a strong increase in SST in the Southern Hemisphere after 2016, which increases the contribution of short-term variability to the total variance. Figure 4 shows the SST variability in the Southern Oscillation region (equatorial part of the Pacific Ocean) according to the Met Office and ERA5 reanalysis data (top) and the change in air temperature in the lower troposphere in the tropics (middle) and globally (bottom) according to the ERA5 and MERRA2 reanalysis data. In the South Oscillation region, the SST trend is insignificant if one focuses on the 1 or 5% significance level, and the Southern Oscillation phases correspond to the SST maximum (El Niño phase) in 1982-1983, 1986-1987, 1991-1992, 1994-1995, 1997-1998, 2002-2003, 2004-2005, 2007, 2009-2010, and 2014-2016, and the minimum SST (La Niña phase) in 1989, 1996, 2000, 2006,2008, 2011, and 2018 [51,67-69]. Almost all variance is explained by the short-term variability mostly associated with ENSO. The change in air temperature in the tropical zone shown in Figure 4 (middle) demonstrates a positive trend, significant at a significance level of 1%, with values of 0.2 degrees per decade for the MERRA2 reanalysis and 0.13 degrees per decade for the ERA5 reanalysis. The short-term variability of air temperature in the lower troposphere of the tropics clearly follows the phases of the Southern Oscillation. It can be noted that the ERA5 data exceed the MERRA2 data by 0.5 degrees in the early 1980s of the XX century and by 0.1 degrees at the end of the 2010s of the XXI century. The difference is decreasing due to the fact that the trend according to MERRA2 data is almost twice as large as according to ERA5.
The global variability of the temperature of the lower troposphere (bottom of Figure 4) is close for the MERRA2 and ERA5 reanalysis data. The trend is about 0.15 degrees per decade, which is slightly less than ERA5 tropical trend and slightly larger than MERRA2 tropical trend, and explains about 65-70% of the total variance. The short-term variability in the temperature of the lower troposphere on a global scale, as well as in the tropics, follows the SST variability in the ENSO region and describes about 30% of the total variance. Figure 5 presents the variability of air temperature in the lower troposphere of the middle latitudes of the Northern Hemisphere (top) and the Southern Hemisphere (bottom) according to the MERRA2 and ERA5 reanalysis data. Estimation of trends in surface air temperature shows their positive values of slightly more than 0.2 degrees per decade in the Northern Hemisphere and two times less in the Southern Hemisphere (0.07-0.09 degrees per decade). Trends are significant with a 1% significance level. The values of the air temperature trends in the lower troposphere at middle latitudes are in good agreement with the SST trends in the middle latitudes of the Northern Hemisphere and are twice as large as the SST trends in the middle latitudes of the Southern Hemisphere ( Figure 2). The short-term variability of the surface air temperature in the middle latitudes of the Northern Hemisphere also, on the whole, coincides with the short-term variability of the SST in the middle latitudes, with maximums in 1981, 1983, 1990, 1998, 2005, and 2014 and minimums in 1985, 1992-1993, 1996, 2009 The interannual variability of the surface air temperature in the polar regions according to the MERRA2 and ERA5 reanalysis data is shown in Figure 6. In the Arctic region, there are maximum positive temperature trends of about 0.45 degrees per decade, which is more than two times higher than the global average for both air temperature ( Figure 4) and SST (Figure 1). The maximum trends in the Arctic correspond to the phenomenon of the Arctic amplification [52,64,65,69]. At the same time, the temperature trend in the Arctic is three times higher than the maximum local SST trend according to ERA5 data and 10 times higher than the minimum SST trend according to Met Office data ( Figure 3). The short-term variability of air temperature in the lower troposphere of the Arctic has maximums in 1981, 1984-1985, 1988, 1990, 1995, 1998, 2002-2003, 2005, 2010, 2012, and 2016. Compared to the El Niño phases of the Southern Oscillation, these maximums coincide for the El Niño phases of 1988, 1995, 1998, 2003, 2010, and 2016. As in the middle and tropical latitudes, ERA5 data systematically exceed MERRA2 values by about 0.5 degrees.
In Antarctica, the surface air trend is weakly negative and insignificant for the 5% significance level. More than 99% of the total variance is explained by short-term variability in the Southern polar region, while in the Northern polar region, more than 65% of total variance is described by the trend. The short-term variability of the temperature of the lower troposphere in the Antarctic region has maximums in 1984, 1988, 1991, 1996, 1998, 2002, 2005, 2007, 2011, and 2016. These highs correspond to El Niño events in 1988, 1988, and 2016. Unlike other areas of the globe, in Antarctica, MERRA's surface air temperature data outperform ERA5 data by more than two degrees.
The temperature of the lower stratosphere ( Figure 7) as a whole across the globe is characterized by a cooling during the period from 1980 to 2019. The magnitude of the trend is about −0.15 degrees per decade for MERRA2 data and −0.3 degrees per decade according to ERA5 data. Trends are significant ata 1% significance level. More than half of the total variance is described by short-term variability both in the tropics and on a global scale. The short-term variability of the temperature of the lower stratosphere on a global scale is characterized by the presence of maxima in 1982, 1988, 1992, 1998, 2002, 2009, and 2012-2013 and 2016. At the end of the 20th century, the highest temperature maxima of the lower stratosphere, both in the tropics and on a global scale, fall on the periods after the major volcanic eruptions El Chichon (Mexico, 1982) and Pinatubo (Philippines, 1991). A sharp increase in the content of stratospheric aerosol after these largest volcanic eruptions led to a trapping of the incoming solar radiation in the stratosphere and to its heating. In the XXI, in some years, the opposite phase change in the temperature of the lower stratosphere in the tropics and on average over the globe is noted. In this regard, 2002, 2011, 2012, and 2015 can be noted. The El Niño phenomenon manifests itself in the global average temperature maximums of the lower stratosphere in 1988, 1998, and 2016. In the lower stratosphere of the tropics, no significant connection with the El Niño phenomenon was noted. In contrast to the temperature of the lower troposphere, where the MERRA2 and ERA5 data have a systematic difference of up to 0.5 degrees, in the stratosphere, the data of the two considered reanalyzes practically do not differ.  The interannual variability of the temperature of the lower stratosphere in the middle latitudes of the northern (top) and Southern (bottom) Hemispheres according to the MERRA2 and ERA5 reanalysis data is presented in Figure 8. For both hemispheres, the temperature trends in the lower stratosphere are negative at about −0.2 degrees per decade and significant at a 1% significance level. Trends according to ERA5 data are slightly high than trends according to MERRA2 data and amount to −0.2-0.3 degrees per decade. The short-term variability is characterized by the presence of temperature maxima in the lower stratosphere of the middle latitudes of both hemispheres in 1982-1983, 1991-1992, 2005, 2007, 2009, and 2014. An increase in the temperature of the lower stratosphere as a result of large volcanic eruptions at the end of the 20th century is also noted in mid-latitudes, with El Chichon more influencing the temperature in the Northern Hemisphere, and Pinatubo in the Southern Hemisphere. At the same time, the global maximum of 2002 (Figure 7) appears only in the temperature of the lower stratosphere of the middle latitudes of the Southern Hemisphere. As well as for the whole globe and in the tropics, the data on the temperature of the lower stratosphere of middle latitudes differ little for the MERRA2 and ERA5 reanalysis data. The interannual variability of the temperature of the lower stratosphere of the polar regions according to the MERRA2 and ERA5 reanalysis data is shown in Figure 9. Temperature trends are insignificant for the 10% significance level and, at the same time, weakly positive for the Northern Hemisphere and negative for the Southern Hemisphere. More than 95% of total variance is explained by the short-term variability in polar regions. Among the global maximums in the temperature of the lower stratosphere, one can note the maximums of 1991-1992, 1998, 2009, 2012-2013, and 2016, which are manifested in both hemispheres and the maximums of 1998 and 2013 in the polar zone of the Northern Hemisphere and the maximum of 2002 in the Southern Hemisphere. Compared with the phases of the El Niño phenomenon, its most pronounced manifestation in 1998 is characterized by the presence of a temperature maximum in the lower stratosphere of the Arctic and a minimum in the lower stratosphere of the Antarctic. A similar picture is observed in 2010, and in 2016, the maximum temperature of the lower stratosphere appears in the polar lower stratosphere of both hemispheres. In general, it can be noted that, despite the fact that the temperature of the lower stratosphere of the polar regions is characterized by a well-pronounced interannual variability determined by dynamic processes affecting the stability of the polar vortex, such as wave activity and sudden stratospheric warming events, global variability, largely determined by processes in the tropical zone, also manifests itself in polar regions. The short-term variability in the polar regions is largely determined by the dynamics of the polar vortex, the stable state of which leads to a cooling of the polar stratosphere due to the isolation of the polar region from warmer middle latitudes, and the unstable polar vortex leads to a warmer stratosphere [13]. The polar vortex is stable almost every year throughout the entire winter in Antarctica and unstable in the Arctic, which is caused by the higher wave activity in the Arctic, leading to the rapid destruction of the polar vortex that forms at the beginning of winter. Sudden stratospheric warming events (SSW), which regularly occur in the Arctic and rarely in the Antarctic, can also lead to the destruction of the polar vortex. As a result, the average temperature of the lower stratosphere in the Antarctica is 6-7 degrees lower than in the Arctic (Figure 9). Against the background of a higher temperature of the lower stratosphere in the Arctic, individual years with a stable polar vortex (1997 and 2011) stand out, when the temperature of the lower stratosphere is significantly lower than the average climatic temperature (Figure 9, top) during the winter-spring period, which led to the formation of ozone anomalies [53]. In Antarctica, on the contrary, there are several years when the temperature is significantly higher than the average climatic, mainly as a result of the SSW, which was the strongest in 2002 (Figure 9, bottom). It should be noted that during the last decade (2010-2020), the destruction of the polar vortex, leading to a sharp increase in the temperature of the lower stratosphere, occurs in Antarctica more often than in previous decades (Figure 9, bottom).

Numerical Modeling of the Simultaneous Influence of Sea Surface Temperature Variability and an Increase in the Content of Greenhouse Gases on the Atmospheric Temperature during 1979-2019
To numerically simulate the simultaneous effect of sea surface temperature variability and an increase in the content of greenhouse gases on the atmospheric temperature, the chemistry-climate model of the INM RAS-RSHU was used. The calculations were performed for the period from 1979 to 2019 [52,54]. To take into account the influence of greenhouse gases for this period, the variability of their surface concentrations was set in accordance with the IPCC Climate System Scenario Tables [50,54,65,66], and to take into account the effect of sea surface temperature (SST) and sea ice coverage (SIC), the reanalysis data were used adopted to the CCM grid. The variability of the chemical composition of the atmosphere was taken into account by setting the surface concentrations of chemically active components in accordance with the data of WMO-UNEP 2018. In the baseline scenario, two calculations were performed with SST and SIC from two different re-analysis data: Met Office and ERA5. To assess the impact of climatic changes over the past 40 years on the temperatures of the lower troposphere and lower stratosphere, the model simulated values of the temperatures of the lower troposphere (925 hPa) and lower stratosphere (20 hPa), averaged over 11 years (first and last solar cycle) at the beginning of the 11-year climatic period (solar cycle) (1980)(1981)(1982)(1983)(1984)(1985)(1986)(1987)(1988)(1989)(1990) and at the end of the climatic period (2009-2019) were compared. The first period (1980)(1981)(1982)(1983)(1984)(1985)(1986)(1987)(1988)(1989)(1990) is characterized by lower values of SST, atmospheric loading of greenhouse gases (CO 2 , CH 4 , N 2 O), and emissions of ozone-depleting substances (ODS) than in the second period (2009-2019). The results of the model calculations were compared with the re-analysis data from MERRA2 and ERA5. The difference between the investigated characteristics is shown in Figure 10. The results of the calculations performed with the CCM showed that in the lower troposphere in the Arctic there is a much stronger warming than in the Antarctic, which corresponds to the well-known hypothesis of "Arctic amplification". For both SST and SIC reanalysis sets, calculations with a chemistry-climate model show a maximum increase in the temperature of the lower troposphere in the Greenland sector, as well as in the eastern part of Siberia, and a decrease in SST and air temperature west of South America. The reanalysis data (Figure 11) also support the hypothesis of Arctic amplification and local cooling in the Antarctic lower troposphere.
In the Antarctic, in contrast to the Arctic, according to the results of calculations with a chemistry-climate model and according to reanalysis data, the warming in the lower troposphere is the lowest compared to other latitudes, and even a cooling is noted in some areas of the southern polar zone. The smallest changes in the surface temperature are obtained in the calculations with the CCM using the SST and SIC from the data of the ERA5 reanalysis. The significant difference between the two polar regions suggests that, on the one hand, the dynamical factors that determine the difference between the polar zones of the northern and Southern Hemispheres play a critical role in the variability of the temperature of the lower troposphere in the Arctic and Antarctic, and, on the other hand, the ability of chemical-climatic model with a given variability of influencing parameters to describe the observed features of temperature variations in the surface layer of the polar atmosphere. In the lower stratosphere, the results of model calculations, as well as the reanalysis data, show, in contrast to the lower troposphere, a cooling in most of the globe. At the same time, as well as in the surface layer, there is a significant difference between the northern and southern polar zones. While in Antarctica the least cooling of the stratosphere, and in some places even warming, is recorded in comparison with other latitudes, in the Arctic, there is a significant difference between the western and eastern parts. In the eastern part of the polar zone of the Northern Hemisphere, the results of model calculations show the maximum cooling of the stratosphere for the entire globe, and in the western zone there is minimum cooling, while in the north of the American continent and Greenland, there is even warming of the stratosphere. At the same time, while in the lower troposphere the use of different reanalysis data for SST and SIC did not lead to a significant difference in the results of model calculations, in the stratosphere, the results of modeling temperature variability for different SST and SIC differ significantly, especially for the polar regions. This suggests that due to the important role of dynamic processes, which are significantly influenced by the exchange of heat, mass, and angular momentum between the ocean and the atmosphere, the variability of the SST and SIC is important not only for the troposphere, which directly borders the ocean, but also for the stratosphere, which is more than 10 km away from the ocean. It should be noted that the ERA5 and MERRA2 reanalysis data also differ more in the stratosphere than in the troposphere, both among themselves and with the results of model calculations with different SST.
Numerical modeling of the effect of the sea surface temperature and sea ice coverage on the variability of the composition and structure of the stratosphere at different latitudes with different data from the reanalysis of SST and SIC showed that these parameters play an important role for the temperature not only in the troposphere, but in the stratosphere as well.
For an in-depth study of the influence of SST and SIC on the content of atmospheric gases in the Arctic and Antarctic and highlighting their influence on the variability of temperature against the background of other influencing factors, additional calculations were carried out for three additional scenarios with the assignment of SST and SIC according to the Met Office reanalysis data. In the first scenario, the values of carbon dioxide loading and SST/SIC were fixed at 1979 level. This scenario model simulation was performed to assess the impact of other but CO 2 and SST/SIC changes during 1979-2019. Among these factors are emissions of other greenhouse (CH 4 and N 2 O) and ozone-depleting gases, aerosols, and solar activity. In the second additional model experiment, the CO 2 year to year variability in accordance with the IPCC Climate System Scenario Tables was taken into account, while the SST and SIC remained at the 1979 level. The purpose of this model experiment was to assess the separate effect of an increase in carbon dioxide content on tropospheric and stratospheric temperatures. In the third additional scenario, carbon dioxide loading was fixed at the 1979 level, while SST and SIC were changed according to Met Office data. This experiment made it possible to evaluate the separate SST/SIC impact on the atmospheric circulation and, through it, on the tropospheric and stratospheric temperature.
In the experiment with fixed CO 2 and SST/SIC ( Figure 12A,B) in the lower troposphere, the temperature increased fairly uniformly across the globe with the centers of maximum increase in the Arctic and Subarctic and cooling in the Antarctic ( Figure 12A). In the stratosphere, at the same time, there is a slight uniform cooling within −0.5 degrees in tropical and middle latitudes, as well as centers of temperature decrease of more than −1.5 degrees in Eastern Siberia and the Far East, insignificant temperature changes in Northern Europe and Antarctica ( Figure 12B). In addition, in some regions of Antarctica, under the influence of factors not related to CO 2 and SST/SIC, an increase in the temperature of the lower stratosphere by more than 0.5 degrees is noted. The results of this model experiment showed that many local features of changes in temperature in polar regions, in particular in the north of the American continent, Siberia, and Antarctica, are not a consequence of changes in the content of CO 2 and SST/SIC but are associated with other global factors and features of local variability, in particular, circulation in the regions.
In the experiment with varying CO 2 emissions and fixed SST/SIC, the change in the temperature of the lower troposphere ( Figure 12C) is similar to that noted in the first experiment ( Figure 12A). At the same time, warming due to an increase in CO 2 loading is slightly higher than due to other greenhouse gases and other factors, and the regions of maximum warming and cooling are also located in the Arctic and Antarctic. Since the content of carbon dioxide increases uniformly over the entire globe, such localization may be a consequence of local dynamic conditions superimposed on the change in the radiation balance with an increase in the content of CO 2 . In the lower stratosphere, an increase in carbon dioxide content leads to a significant cooling with maximum values of up to −4 degrees in the Antarctic and more than −2 degrees in certain regions of the Arctic ( Figure 12D).
In the scenario with fixed CO 2 and varying SST/SIC in the lower troposphere, the air temperature increases more than in the previous scenarios ( Figure 12E). The greatest warming is noted in the Arctic (Arctic amplification), but in Antarctica, for this scenario, the temperature of the lower troposphere increases more than in middle and tropical latitudes. In the lower stratosphere ( Figure 12G), a slight cooling persists, similar to that noted in the first scenario, but in addition, areas of temperature decrease in the tropics and subtropics appear. At the same time, warming in some regions of Antarctica increases, while in the Arctic and Subarctic, the area of weak temperature changes increases in comparison with the first scenario. In addition, areas of increased temperature appear in eastern Eurasia and Alaska.
In the baseline scenario, when all factors are taken into account ( Figure 12G,H), the Arctic amplification is well expressed in the lower troposphere, which manifests itself in maximum warming in the Arctic and minimum warming, and even cooling in some regions, in Antarctica. Thus, the results of numerical experiments with the CCM revealed that only taking into account the simultaneous influence of the CO 2 and SST/SIC variability makes it possible to sufficiently reliably reproduce the observed changes in the temperature of the lower troposphere and the differences between the northern and Southern Hemispheres ( Figure 12G). In the lower stratosphere, for the baseline scenario, SST/SIC variability partially compensates the cooling due to increased CO 2 ( Figure 12D,H).

Discussion
The results of numerical experiments have shown that the air temperature in the lower troposphere is basically determined by the variability of SST. Correlation analysis confirms this conclusion. The correlation coefficient between SST anomaly and air temperature anomaly at the level of 925 hPa is 0.85, both for the reanalysis data and for the simulation results. If we estimate the correlation coefficients for some regions, then for the Pacific Ocean it is equal to 0.86 and for the tropical part of the Pacific Ocean (ENSO region, 5 • S-5 • N, 170 • W-120 • W) it is 0.90. For the Atlantic Ocean it is 0.82 and for the Indian Ocean it is 0.8. For the Arctic and Antarctic, it is 0.81 and 0.7, respectively. For areas over land, the correlation is less-the correlation coefficients are 0.65 and 0.7 for the troposphere over Eurasia and America. The correlation between the SST anomaly in the El Niño region and the air temperature anomaly in the troposphere of the Arctic and Antarctic are −0.4 and −0.03, respectively, and the air temperature anomaly in the Northern and Southern Hemispheres-−0.4 and 0.17, respectively.
In the stratosphere, the SST anomaly variability does not significantly affect the air temperature anomaly-the correlation coefficients range from 0.2 to 0.3. For the Northern Hemisphere, the correlation coefficient between SST anomaly and air temperature anomaly in the stratosphere is −0.37. At the same time, in El Niño years, the correlation coefficients between SST anomaly and air temperature anomaly in the Arctic stratosphere are 0.4-0.5, in La Niña years-0.4, in years of the neutral phase-from −0.45 to 0.4. For the Southern Hemisphere, the correlation coefficient is 0.2.
As for the influence of the level of carbon dioxide on air temperature, in the troposphere it has a significant effect on changes in air temperature. So, the correlation coefficient between the concentration of carbon dioxide and air temperature anomaly is 0.65 for the whole planet. The largest value of this coefficient was obtained for the Northern Hemisphere-the coefficient is 0.75, which is explained by the presence of a large number of carbon dioxide sources, mainly of anthropogenic origin [16][17][18]. For the Southern Hemisphere, it is 0.65. For the stratosphere, the correlation coefficient is −0.5-the correlation is negative, which corresponds to the cooling of the stratosphere, while it is more significant than between the air temperature in the stratosphere and the SST. At the same time, significant differences in the distribution of air temperature in the stratosphere for experiments with fixed and non-fixed CO 2 levels indicates a greater significance of changes in the level of CO 2 for the stratosphere than changes in SST [33,34,52].
Comparison of the results of additional scenarios of model calculations with the baseline scenario, in which the variability of all influencing parameters is set, allows us to conclude that, due to the peculiarities of circulation in the polar regions, an increase in the carbon dioxide content leads to less cooling of the stratosphere in the Arctic compared to other latitudes and more cooling in the Antarctic, if we do not take into account the influence of variability in sea surface temperature and ice surface area on circulation in the north polar region. The cumulative effect of climatic changes associated with the influence of the variability of carbon dioxide and the temperature of the sea surface and the area of ice coverage leads to the localization of the circumpolar vortex, as a result of which in the eastern hemisphere, especially in northern Siberia, the strongest cooling of the stratosphere occurs [33,34,70].

Conclusions
Reanalysis data, as well as numerical simulations with a global chemistry-climate model of the lower and middle atmosphere were used to study the effect of sea surface temperature and carbon dioxide variability on the temperature of the lower troposphere and low stratosphere during the period from 1980 to 2019.
Analysis of sea surface temperature trends based on reanalysis data revealed that, on average, a significant positive SST trend of about 0.01 degrees per year is observed over the globe. At the same time, in the middle latitudes of the Northern Hemisphere, the trend (about 0.02 degrees per year) is 2 times higher than the global average, and 5 times higher than in the Southern Hemisphere (about 0.004 degrees per year). At polar latitudes, opposite SST trends are observed in the Arctic (positive) and Antarctic (negative). The short-term variability of the SST on average over the globe coincides with a good degree of accuracy with the phases of the Southern Oscillation in the tropics, while in the middle and polar latitudes, the phases of the Southern Oscillation are noticeable for the most powerful El Niño and La Niña phenomena in the short-period variability in the Northern Hemisphere and are practically invisible in the Southern Hemisphere.
The air temperature in the lower troposphere is characterized by a positive significant trend on average over the globe (about 0.015 degrees per year), slightly exceeding the global SST trend. For some latitudinal areas, the maximum trend in the temperature of the lower troposphere is noted in the Arctic (about 0.45 degrees per year), and the minimum-in the Antarctic (insignificant negative trend). In the middle latitudes of the Northern Hemisphere, the trend in the temperature of the lower troposphere (about 0.02 degrees per year) is almost three times higher than the trend in the temperature of the lower troposphere of the Southern Hemisphere (about 0.007 degrees per year). The short-term variability of the mean global temperature of the lower troposphere is in good agreement with the phases of the southern oscillation, which has a tropical nature. At the same time, in middle and, especially, polar latitudes, the dependence of short-period fluctuations in the temperature of the lower troposphere on the phases of the southern oscillation is much less than for tropical latitudes and on average over the globe. The reanalysis data for the temperature of the lower troposphere have systematic differences: Everywhere, except for the Antarctic region, the values according to the ERA5 data exceed the values according to the MERRA2 data.
In the lower stratosphere, a significant cooling is observed during 1980-2019 period, i.e., the direction of the temperature trend is opposite to its trend in the lower troposphere [32,33,70]. The cooling rate in the lower stratosphere on a global scale is about −0.02 degrees per year, and these values are close both in the tropical zone and in the middle latitudes of both hemispheres and in Antarctica. Only in the lower stratosphere of the Arctic is there an insignificant, weakly positive trend. The short-term variability of the temperature of the lower stratosphere is less dependent on the phases of the southern oscillation, and to a greater extent is determined by the circulation features of different hemispheres and latitudinal zones. The reanalysis data for the temperature of the lower stratosphere differ much less from each other than was observed in the lower troposphere. Only a slightly lower linear trend coefficient can be noted for ERA5 data compared to MERRA2 data.
The results of model calculations with a global chemistry-climate model demonstrated that in the lower troposphere, an increase in air temperature is determined to a greater extent by a change in the temperature of the sea surface and its ice coverage than by an increase in the content of carbon dioxide [33,34,70]. At the same time, the Arctic amplification is determined by the cumulative effect of both the sea surface temperature and sea ice coverage, and an increase in the content of greenhouse gases. On the other hand, if CO 2 changes are not taken into account, then the air temperature in the stratosphere will change little. Thus, the study showed that changes in the CO 2 content are a significant factor for the observed cooling of the stratosphere, the fixation of which makes the temperature of the stratosphere almost constant during 1980-2019, and the change in SST has only an insignificant effect on it.  Data Availability Statement: Reanalysis data are available from RSHU data service (ra.rshu.ru/ra (accessed on 23 March 2021)). Results of model calculations are available upon request (smyshl@rshu.ru).