The paper discusses air (T_a) and sea surface temperature (SST) year-to-year variability due to warming of the Kara Sea, using the data from regular observations at the meteorological stations Roshydromet (GMS) in 1978–2017, NOAA optimum interpolation and reanalysis data. We use the methods of cluster, correlation analysis and Empirical Orthogonal Functions (EOF). We investigate possible cause and effect relationships of these changes with the variations of the wind field components, climatic indices and the sea ice concentration field. The cluster analysis of the three main EOF components has allowed us to identify four areas on the basis of the nature of changes of the water temperature anomalies field. The climatic changes in these areas, in the coastal and island zones of the Kara Sea have manifested themselves in the steady increase of the annual air temperature at GMS from 0,47–0,77 °C/10 years on the southwest coast to 1,33–1,49 °C/10 years in the north of the sea. This is equivalent to warming from 1,9 to 6,0 °C in the last 40 years. For the open sea the value of the T_a trend is about 1,22 °C/10 years, which corresponds to an increase in the average T_a by 4,9 °C in the last 40 years. This value is approximately 3 times greater than that for all the Northern hemisphere for the same period.

Annualy, the maximal trend was observed in November and April mainly and exceeded 2–3 °C/10 years at some of the stations. We identify anomalously warm (2016 and 2012) and anomalously cold (1978, 1979, 1992 and 1998) years: the warmest year was 2012, the coldest — 1979. Positive SST trends were observed over all the sea area during the warm period of year (to 1 °C/10 years). SST increased to 2,4 °C, which is approximately 1,5 times greater than the corresponding SST values for the Northern hemisphere. The maximum SST trend (0,4 °C/10 years) was observed in the northwest and southwest parts of the sea. From June to August the trends of SST exceed the annual ones 1,5–2 times. Interannual SST and T_a variations are characterized by close correlation links. Until approximately 1998–2004 the warming was rather insignificant, and after that the growth rate of T_a and SST increased many fold. Apparently it indicates changes in the mode and the large-scale atmospheric circulation in the early 2000s. We also observed a trend of strengthening of the southern wind during the cold period of the year and the northern one — in the warm period (0,5–0,6 m/s in 40 years). It is shown that there is a close correlation between the T_a increase and the changes in the meridional component of the wind speed during the cold period of the year for all the sea areas. For the warm period it is statistically insignificant both for T_a and SST. For the cold season we observed a contribution of the large-scale mode of atmospheric circulation into the variability of V component of the wind speed. The conribution was expressed through the indeces NAO, SCAND, Pol/EUR, AZOR, ISL and the differences of ISLSIB. For the warm season this contribution is expressed through the NAO, SCAND and AO only. For the warm period we showed statistically significant correlation between the increase in SST, T_a and the processes parametrized by the AMO, EA/WR and AZOR indeces. For the cold period the indeces are AMO, Pol/Eur, SIB and ISL SIB. The interannual variations of the sea ice concentration field are characterized by close correlation with T_a changes both in the annual cycle and during the periods of ice cover formation and evolution (R = –0,7... –0,9). For these periods we showed statistically significant relationships between the first EOF mode fluctuations and two climatic indeces — AMO (R = 0,5) and Pol/Eur (R = 0,4). The relationships between the temporary variability of the sea ice concentration and the wind field characteristics are weaker and statistically significant only for the meridional component of the wind speed (R = –0,4).

The structure of the long-period variability of the ice cover of the Barents and Greenland Seas over a long series of observations from 1930 to 2017 is analyzed. In both seas, there is a significant negative linear trend of ice cover for both the winter and summer seasons. Average for the period of 1950–2016 intra-annual changes in ice coverings demonstrate the conjugacy of the seasonal cycles of the Greenland and Barents Seas, but with certain differences. Three homogeneous groups with a similar character of intra-annual changes in the ice area are identified for each sea. Identified succession in a state of ice cover for 2 years.

The conjugacy of changes in the average decadal values of sea ice cover in April and August with the average decadal indices of atmospheric circulation AO, AD, PNA, NAO and the index of the thermal state of the North Atlantic AMO is shown. Spectral analysis of the winter and summer ice cover of the Greenland and Barents Seas for the period 1930–2016 confirmed earlier received cyclical fluctuations of 22, 9–11 and 6–7 years.

Cross-correlation analysis established a close relationship between the longitudinal changes in the ice cover and the average annual values of the following astrogeophysical parameters, the longitude coordinate of the Earth pole position Y, the Earth axis nutation indices dEps and dPsi, the Earth rotation speed index lod (length of day), Sun solar activity index (annual Wolf number) , the average for six months, the distance from the Sun to Earth in the summer SX-III and the winter SX-III periods. Significant correlation coefficients are quite large (R = |0,30| – |0,56|) for both seas, comparable to the correlation coefficients between the ice cover and average annual air temperature T, show the reality of the ice cover mediated reaction to changes in astrophysical factors. Statistical equations relating the sea ice cover to hydrometeorological and astrogeophysical factors were obtained by multiple correlation. The overall correlation coefficient varies from R = 0,80 to R = 0,87 AT. The Greenland Sea, the share of astrogeophysical factors in the long-term changes in the ice cover of both the winter and summer seasons exceeded the contribution of hydrometeorological factors by 3–4 times. In the Barents Sea, the contribution to the total dispersion of astrogeophysical factors in the winter period is somewhat less than that of hydrometeorological factors, and in the summer period they exceed only 1.4 times. The authors’ approach opens up the possibility of using it to obtain statistical equations for the diagnosis and forecast of long-term and climatic changes in sea-ice cover.

The data on snow the temperature which was monitored to a depth of 10 m in the vicinity of Vostok Station by the TAUTO autonomous system in 2010–2017 are presented. By analyzing seasonal temperature variations at different depth with the aid of a heat-transfer model we have inferred a relationship between relative thermal conductivity of snow and its porosity at this site. The same approach was also applied to analyze similar data obtained at Dome Fuji station in 1995–1997. It was found that the thermal conductivity of snow layers with identical density is noticeably lower at Dome Fuji than at Vostok, which point to a difference in structural characteristics of snow that determine its thermophysical properties. We demonstrate that the conduction is the dominant heat-transport mechanism which controls the temperature distribution in snow pack on the Antarctic plateau. The obtained parameters of the heat-transfer model can be used for reconstructing the past surface temperature variations from the long-term temperature measurements in the upper 100 m thick layer of the ice sheet.

New data about geomorphological structure and neotectonic movements of the Lena Delta is presented. Thet are based on results of the Russian-German expeditions Lena-2013, Lena-2014, Lena-2015 and Lena-2018 and include geomorphological profiling with high quality satellite instruments and a number of radiocarbon dates on the Sobo-Sise, Kurungnakh, Jangylakh-Sis and Khardang-Sise islands. These islands consist of the Late Pleistocene Ice Complex (IC) remnants eroded by river and sea, and the first terrace of the Delta, which adjoins the remnants. The first terrace started to form 8,000 years ago in the western part of the Delta and was finally formed in the last millennium in the eastern part of the Delta. From the previous works, it is known that the western part of the Delta is higher than the eastern part. In our work, we explain it by the eneven movements of the Earth’s crust in this region. The aim of the paper is to study the quantitative characteristics of the tectonic movements in the Lena River Delta. For this purpose we present the geomorphological schemes and descriptions of the islands mentioned and compare the terraces heights in the different parts of the Delta. In the Late Pleistocene, according to the heights of the IC remnants, the western part of the Delta rose 1 mm per year faster than its eastern part. In the Holocene the speed difference increased to approximately 2 mm per year, which is shown by the terrace surface’s altitude. The amplitude between the western and eastern parts of this surface is about 4 m. Finally, according to 60-years observation period of the water level in the Laptev Sea, the modern speed difference of the western and eastern parts movement in the Lena Delta is 2 mm per year. As a result, we can observe changes in the main flow direction in the Delta channels from the Olenekskaya branch at the beginning of Holocene to the Bykovskaya branch today.

The paper discusses the results of the ground-penetrating radar (GPR) survey carried out in February 2019 in the area of Thala Bay (Larsemann Hills, East Antarctica). Thala Bay is one of the strategic facilities of the Russian Antarctic Expedition (RAE) in the Progress station area as since 2019 heavy cargo has been unloaded here intended for the construction of new facilities at the Vostok station. Transportation of goods to the point of formation of logistic traverses takes place on ice tracks, whose safety must be evaluated taking into account the expanded system of crevasses. In addition, the current track is characterized by a significant slope of the terrain, which also complicates the relocation of heavy equipment.

In February 2019, a GPR survey was carried out within the Thala Bay area to assess the possibility of organizing an alternative section of the route within it. According to the visual observations, this area was characterized by an extensive system of crevasses, the width of which at the surface reached 20-30 cm, and the prevailing longitudinal direction coincided with the direction of the route. The task of the geophysical survey was to map the crevasses not identified by visual inspection and to determine their morphology. According to the GPR data, it was shown that the crevasses within the site are located to the firn layer and are characterized by an irregular shape, significantly expanding at the deeper levels and reaching a width of 6 m. The results of the survey are illustrated with the scheme of the firn thickness which shows location of the crevasses. According to the recommendations of the authors, the section of the glacier is suitable for operation provided the glaciological situation using the GPR method is monitored annualy.

Wilkes Land is a key region for Gondwana reconstruction, however it remains one of the largest regions on Earth with poorest knowledge of geology. This study comprehensively reviews the ICECAP/ IceBridge geophysical data for the Law Dome region including Vanderford and Totten adjacent glaciers over Wilkes Land and their role in obtaining new insight on the East Antarctic geology hidden under the ice cover. We analyzed more than 100,000 line kilometers of new magnetic, gravity and subglacial bedrock topography data that are available through the National Snow and Ice Data Center (USA). The newly acquired data supports our previous idea of the continuous rift structure existence at the southern boundary of Law Dome that runs between Vanderford and Totten Glaciers. The rift length exceeds 400 km and width varies from 50 to 100 km. In accordance with results of depth to Moho estimations and density modelling, for axial part of the rift it is characteristic an essential thinning of the Earth crust thickness, it is raised up to 24–26 km and continue to be elevated along entire length of this structure. The thickness of sedimentary rocks within the rift exceeds 3 km, their high density probably evidence that they were formed during Late Paleozoic – Early Mesozoic. The results of our investigations support tectonic nature of this structure as continuous rift developed since the Mesozoic extension phase (~160 Ma) of the Wilkes Land continental margin. Second distinctive structure is the strong reversely magnetized Law Dome magnetic anomaly with an area of about 9,500 km². This anomaly would map out one of the largest mafic/ultramafic intrusions of the Earth, similar in extent to Norway’s Bjerkreim-Sokndal layered intrusion, the Coompana Block gabbro in Australia, or even the granitic-gneiss complex in the Adirondack Mountains of North America.