1. Introduction
The 1930s, which were part of the early-twentieth-century Arctic warming (ETCAW, usually defined in the literature as the period 1921–50; see, e.g., Yamanouchi 2011), were the warmest decade in the twentieth century in the Arctic (Przybylak 2007). On the other hand, 10-yr running means of surface air temperature (SAT) calculated for the period since the mid-1990s, when the contemporary Arctic warming (CAW) began, exceeded the values observed in the 1930s (see Przybylak 2007, 2016; Przybylak and Wyszyński 2020). The ETCAW was a prevalent subject of study and analysis among scientists of the time. Wood and Overland (2010) reported that in the first half of the twentieth century, over 1000 articles were published. The literature of this period is reviewed in, for example, Jensen (1939), Ahlmann (1948), Veryard (1963), Wállen (1984), Ellsaesser et al. (1986), Przybylak (1996, 2000a, 2002a, 2003, 2016), Wood and Overland (2010), Brönnimann (2015), and Przybylak et al. (2021). Since the 1980s there has been increased interest in the ETCAW period, and particularly in the reasons for the warming (e.g., Grant et al. 2009; Wood and Overland 2010; Yamanouchi 2011; Hegerl et al. 2018; Svendsen et al. 2018; Xiao et al. 2020, and for more details see Table S1 and supplement 1 in the online supplemental material).
Table S1 shows that there is no general agreement on the causes of the ETCAW, despite nearly a century of effort (see also Bokuchava and Semenov 2021). This is particularly true when articles published in the last 30 years are analyzed (Fig. 1). Conversely, the main driving mechanisms had previously generally been agreed upon among scientists. They believed that the ETCAW was an effect of natural forcing—mainly of changes in atmospheric circulation (71.6%) and, to a lesser degree, oceanic circulation (20.6%). Only Callendar (1938) connected this warming to increased CO2 concentrations (Table S1, Fig. 1). One possible reason for the lack of consensus over the driving mechanisms still ascribed as having caused the ETCAW is that we have limited knowledge about (i) the weather and climate (including extremes) in this period, (ii) other components of the Arctic climate system (e.g., sea ice, glaciers, SST, aerosols, vegetation cover), and (iii) some forcings (in particular, solar and volcanic) (see, e.g., Bengtsson et al. 2004; Wood and Overland 2010; Suo et al. 2013; Klaus et al. 2018; Przybylak et al. 2021).

Development of views on the reasons for occurrence of ETCAW in papers published in (a) 1929–2020, (b) 1929–90, and (c) 1991–2020. Key: Numbers denote percentage of publications in which a given view was presented. Natural forcings: AC, atmospheric circulation; OC, oceanic circulation; SIE, sea ice extent; SST, sea surface temperature; SA, solar irradiation or solar activity; VA, volcanic activity; CV, internal climate variability; OT, other. Anthropogenic forcings: GG, greenhouse gases; AA, atmospheric aerosols.
Citation: Journal of Climate 35, 7; 10.1175/JCLI-D-21-0162.1

Development of views on the reasons for occurrence of ETCAW in papers published in (a) 1929–2020, (b) 1929–90, and (c) 1991–2020. Key: Numbers denote percentage of publications in which a given view was presented. Natural forcings: AC, atmospheric circulation; OC, oceanic circulation; SIE, sea ice extent; SST, sea surface temperature; SA, solar irradiation or solar activity; VA, volcanic activity; CV, internal climate variability; OT, other. Anthropogenic forcings: GG, greenhouse gases; AA, atmospheric aerosols.
Citation: Journal of Climate 35, 7; 10.1175/JCLI-D-21-0162.1
Development of views on the reasons for occurrence of ETCAW in papers published in (a) 1929–2020, (b) 1929–90, and (c) 1991–2020. Key: Numbers denote percentage of publications in which a given view was presented. Natural forcings: AC, atmospheric circulation; OC, oceanic circulation; SIE, sea ice extent; SST, sea surface temperature; SA, solar irradiation or solar activity; VA, volcanic activity; CV, internal climate variability; OT, other. Anthropogenic forcings: GG, greenhouse gases; AA, atmospheric aerosols.
Citation: Journal of Climate 35, 7; 10.1175/JCLI-D-21-0162.1
Most up-to-date published studies also try to (or their results allow readers to) compare the scale of warming during the ETCAW against the CAW in the Arctic using annual mean values and, more rarely, monthly or seasonal ones (see, e.g., Hanssen-Bauer and Forland 1998; Jones et al. 1999; Delworth and Knutson 2000; Jones and Moberg 2003; Polyakov et al. 2003; Kuzmina et al. 2008; Box et al. 2009; Wood and Overland 2010; Wood et al. 2010; Yamanouchi 2011; Fyfe et al. 2013; Johannessen et al. 2016; Łupikasza and Niedźwiedź 2019; Brennan et al. 2020; Nordli et al. 2020). Meanwhile, Wood and Overland (2010), Yamanouchi (2011), and recently also Wegmann et al. (2017) suggest that a reliable explanation of the ETCAW requires data with greater resolution. To date, however, we have few studies that use measurement data with greater resolution (daily or subdaily) for climate comparison between the studied periods. The direct reason is the limited availability of such data because the instrumental records of Arctic temperature are brief, incomplete, and geographically sparse (Przybylak 2000a). Indeed, only six records extend back to the second half of the nineteenth century (Upernavik: date of start 1874; Jakobshavn: 1874; Godthảb: 1876; Ivigtut: 1880; Malye Karmakuly: 1882, and Angmagssalik: 1895). In the twentieth century, the first station was established in Spitsbergen in 1911 (Green Harbour). In the 1920s, seven stations came into operation, mainly in the Atlantic region of the Arctic. Following the Second Polar Year (1932/33), most Russian stations were established, while most Canadian stations were founded after World War II. As a result, the existing papers using daily and subdaily data analyze temperature in the ETCAW for only a small area, or even a single site, of the Arctic, and for a short time (usually one year) (see, e.g., Klaus et al. 2018; Przybylak et al. 2018; Araźny et al. 2019; Schweiger et al. 2019; Nordli et al. 2020). For the first time, the present paper uses these kinds of data for sites representing a significantly greater area of the Arctic and a more extended period of observations.
The methods used are similar to those that we used earlier for analyses conducted for historical periods before the onset of the ETCAW (Przybylak and Vízi 2005; Przybylak and Wyszyński 2017) and for the ETCAW period itself (Przybylak et al. 2018). Such an approach provides a longer perspective of weather and climate changes in the Arctic. The World Climate Research Program recently recommended an increased effort to extend the historical observational records and analysis to improve, among others, event attribution capabilities (see National Academies of Sciences Engineering and Medicine 2016). Weather and climate information (including the probability, magnitude, and circumstances of events) is also urgently needed for the so-called counterfactual world (in which human influence on climate is absent). In the Arctic, anthropogenic influences are minimal before the mid-twentieth century, which covers the entire ETCAW.
The main objective of this article is to provide new (or improved) knowledge about rarely studied aspects of SAT conditions (including extremes) in the Arctic during the ETCAW. High-resolution data (daily and subdaily) for selected meteorological stations in five climatic Arctic regions (Treshnikov 1985) have been used. The improved and widened knowledge should be helpful for further investigations, including (i) for the development and validation of climatic models [Latonin et al. (2021) showed that even the multimodel ensemble means in the new generation of high-resolution CMIP6 models do not correctly reproduce the ETCAW (expressed in terms of winter and annual means)] and (ii) for better understanding of the causes and mechanisms of that warming.
This article is structured as follows. Section 2 describes the area, as well as the data and methods used. Section 3 presents the results of the climate comparison between the two analyzed periods (ETCAW and CAW). Section 4 discusses the obtained results, while section 5 presents the most important conclusions and final remarks.
2. Area, data, and methods
The present analysis uses daily and subdaily SAT from six meteorological stations [Barentsburg, Kanin Nos, Tiksi, Ostrov Vrangel, Coppermine, and Ilulissat (Jacobshavn)] located in the five climatic regions (Atlantic, Siberian, Pacific, Canadian, and Baffin Bay) that Treshnikov (1985) distinguished for the Arctic. The analysis characterizes the climate during the ETCAW and its changes in comparison to the CAW periods (Table 1, Fig. 2). The data come directly from the national institutes responsible for meteorological measurements in the Arctic (see Table 1) and therefore are most reliable. The southern boundary of the Arctic and its climatic regions were delimited by the authors of the Atlas Arktiki based on (i) analysis of the distribution of mean long-term seasonal and annual fields of all meteorological variables (in case of SAT and sea level pressure, data were taken from period 1881–1965), which covers the study period analyzed in the paper, and (ii) the character of variability of these variables and factors controlling the climate of the Arctic [for more details, see Treshnikov (1985)]. This means that the weather and climate conditions within the distinguished climatic regions must be similar. Przybylak (1997) found that the radius of the extent of statistically significant correlation coefficients of mean seasonal SAT changes at some Arctic stations [Svalbard Lufthavn (Norwegian Arctic), Ostrov Kotelny (Russian Arctic), and Resolute Airport (Canadian Arctic)] in period 1951–90 was equal to 2000–2500 km in winter and 1500–2000 km in summer. Similar results were also obtained by, among others, Eserkepova et al. (1982), Smirnova and Subbotin (1983), and Subbotin (1983). Spatial correlation maps of mean monthly and annual SATs between most of the analyzed stations and the rest of the Arctic (area > 50°N) were drawn using data from 20CRv3 for the periods 1931–50 and 1991–2010. In Figs. S1a and S1b (see supplement 2 in the online supplemental material) we present some of them for January, July, and the full year. The results confirm the significant correlation of SAT in the Arctic during the two periods up to an average radius of about 1000–1400 km around the selected stations in winter and 700–1000 km in summer (except for Barentsburg in the period 1991–2010, where the radius was about 450 km). For the annual means, these ranges oscillated between 800 and 1400 km (period 1931–50) and 1400–2100 km (1991–2010). Rigor et al. (2000) also found a significant correlation between 12-hourly SAT observations from land stations and over ice (data from buoys, the International Arctic Buoy Program/POLES) from 1979 to 1998 at ranges of up to 1000 km in winter when most of the bias occurs, and 400–450 km in summer. Donat et al. (2013) estimated the same statistics also (which they called “decorrelation length scales”) for the entire globe, but for an extreme temperature index (maximum monthly value of daily maximum temperature) for January, July, and the year. Their results revealed that the correlation of this temperature index is greatest in the polar regions, particularly in the Arctic. The decorrelation length scale values in the Arctic reached about 1100 km for January and 700–800 km for July and the year. The well-documented significant spatial correlation of SAT in the Arctic allows it to be concluded that the different SAT characteristics analyzed for the six stations selected in the paper are representative for large parts of the regions where they lie [i.e., on average within a radius of at least 400–500 km in summer and 1000 km in winter; we also obtained such results for data taken from the ERA5 (Hersbach et al. 2020) reanalysis for period 1991–2010; not shown].

Location of selected meteorological stations operating in the Arctic during the ETCAW and CAW study periods for which good-quality daily and subdaily air temperature data are available. The Arctic and its climatic regions are defined according to Treshnikov (1985). The southern Arctic boundary is delimited based on analysis using long-term averages of all meteorological variables, their seasonal cycles, and variability characteristics. The thick line is the boundary of the Arctic; thin lines are boundaries between climatic regions.
Citation: Journal of Climate 35, 7; 10.1175/JCLI-D-21-0162.1

Location of selected meteorological stations operating in the Arctic during the ETCAW and CAW study periods for which good-quality daily and subdaily air temperature data are available. The Arctic and its climatic regions are defined according to Treshnikov (1985). The southern Arctic boundary is delimited based on analysis using long-term averages of all meteorological variables, their seasonal cycles, and variability characteristics. The thick line is the boundary of the Arctic; thin lines are boundaries between climatic regions.
Citation: Journal of Climate 35, 7; 10.1175/JCLI-D-21-0162.1
Location of selected meteorological stations operating in the Arctic during the ETCAW and CAW study periods for which good-quality daily and subdaily air temperature data are available. The Arctic and its climatic regions are defined according to Treshnikov (1985). The southern Arctic boundary is delimited based on analysis using long-term averages of all meteorological variables, their seasonal cycles, and variability characteristics. The thick line is the boundary of the Arctic; thin lines are boundaries between climatic regions.
Citation: Journal of Climate 35, 7; 10.1175/JCLI-D-21-0162.1
Location of Arctic meteorological stations whose daily data series are used in the present article, and their sources.


The source data were mean daily air temperature (MDAT) and maximum and minimum daily temperature (TMAX and TMIN, respectively) for all stations except Ilulissat, for which only TMAX and TMIN were available. MDAT was calculated using the following simple formula for this station: MDAT = (TMAX + TMIN)/2. Only the Ilulissat and Kanin Nos stations cover the entire ETCAW, but there are many gaps in the analyzed series, particularly for the 1920s (both stations) and 1930s (only Ilulissat). Other stations were established at the end of the 1920s (Coppermine and Ostrov Vrangel) or shortly after the Second International Polar Year 1932/33 (Barentsburg and Tiksi). For greater comparability of SAT data between the mentioned stations, data were extracted from 1934 to 1950. This period represents the second part of the ETCAW, which is its warmest and climatically most stable subperiod. However, there are also some crucial gaps in the data during World War II, which are particularly large for Barentsburg station (August 1941–November 1947). Because of these gaps, there was no period of at least 10 years for which a continuous sequence existed for all stations. We were able to find only data from 10 individual years. Unfortunately, the gaps also made it impossible to find data for a single set of years common to all stations (see Table 1). These differences in groups of utilized years may have introduced some biases.
The almost complete data series for 1934–50 for Kanin Nos station at the edge of the European Arctic was used to estimate the scale of possible biases. Using these data, we estimated that biases in mean annual and seasonal SAT values at all stations (except Ilulissat) do not exceed ±0.1° and ±0.4°C, respectively. On the other hand, the set of years used for Ilulissat data should introduce a cold bias of about 0.5°C in annual values, on the condition that year-to-year SAT changes in both regions (the European Arctic and the western coast of Greenland) were similar. Such an assumption for these two regions is probably wrong (differences in the influence of the NAO; see Przybylak 2000a). Therefore, we also calculated the differences in mean seasonal and annual SAT between the period 1941–45 and other years from the ETCAW period that were available for Ilulissat station. For both periods, the mean yearly SATs at Ilulissat were the same. Small differences (up to ±0.2°C) were recorded in spring, summer, and fall; winter was 0.4°C colder during World War II. We also checked the influence of these different sets of years on diurnal temperature range (DTR) values in Kanin Nos and Ilulissat. The biases were significantly smaller than those for mean values and did not exceed ±0.2°C for seasonal and annual means. Negligible biases were also found for continentality and oceanicity indices (about 1.5%–2.0% and 0.1, respectively).
The most significant yearly number of gaps in daily data during the ETCAW was recorded at Ilulissat but did not exceed 20 days. Gaps of longer than 10 days occurred only in two years: 1940 and 1948. In the entire 10-yr set of daily data at that station, gaps constitute only about 2%. The only other very small gaps (0.1%–0.2%) were found for Barentsburg and Coppermine, while the data from the other stations are complete. We estimate thus that the influence of those gaps on the presented results is negligible. In accordance with expectation data availability is better for the CAW than for the ETCAW. For the CAW, a continuous 10-yr series of data was used for comparison purposes. Because there is a large gap in the years 2005 and 2006 for Ilulissat station, we took the first available period for analysis (i.e., 2007–16). Detailed dates of gaps in the four thermal parameters are shown in Table S2 (supplement 1).
The following statistics have been calculated using the abovementioned thermal parameters:
-
Monthly, seasonal, and annual means.
-
Annual temperature range (ATR), by subtracting the coldest month’s mean temperature from the warmest month’s mean temperature.
-
Continentality/oceanicity of climate, using indices proposed by Ewert (1972) and Marsz (1995). The continentality index (K) proposed by Ewert was calculated according to the formula: K = [ATR − (3.81 sinφ + 0.1)/(38.39 sinφ + 7.47)] × 100%, while the oceanicity index of climate (Oc) was calculated according to the formula Oc = (0.732φ + 1.767)/ATR. In both formulas, ATR is the annual temperature range, and φ geographical latitude. Based on the distribution of values of Oc around the world, Marsz (1995) distinguished the following types of climate: ultracontinental (Oc < 1), continental (Oc 1.00–1.99), suboceanic (Oc 2.00–2.99), oceanic (Oc 3.00–3.99), and ultraoceanic (Oc > 3.99). For more details, see Przybylak et al. (2014).
-
Day-to-day MDAT variability (DDTV); two methods were used to describe this characteristic: 1) modulus of MDAT change from one day to the next, and 2) standard deviation (SD) calculated for series of MDAT data taken from each month and season (MAM, JJA, etc.). To investigate the reasons for the highest observed change in DDTV during the ETCAW, we utilized the Twentieth Century Reanalysis version 3 (20CRv3; Slivinski et al. 2019) as synoptic background.
-
Number of days with DTR of ≥5°C, ≥10°C, and ≥15°C; DTR was calculated by subtracting TMIN from TMAX.
-
Number of so-called characteristic days, distinguished according to criteria proposed originally by Przybylak and Vízi (2005) and slightly modified for the present paper:
-
TMAX > 15°C = exceptionally warm day
-
TMAX > 10°C = very warm day
-
TMAX > 5°C = warm day
-
TMIN > 0°C = no-frost day
-
TMAX > 0°C and TMIN ≤ 0°C = slight frost day
-
TMAX < 0°C = frost day
-
TMAX < −10°C = cold day
-
TMAX < −20°C = very cold day
-
TMAX < −30°C = exceptionally cold day
-
-
Thermal seasons (winter, spring, summer, and fall), using the following criteria proposed by Baranowski (1968):
-
winter: MDAT ≤ −2.5°C
-
spring and fall: −2.5°C < MDAT < 2.5°C
-
summer: MDAT ≥ 2.5°C
The threshold values were distinguished based on analysis of the yearly cycle of MDATs in Hornsund (1957–60) and their monthly distributions stratified into 2°C intervals.
The onset and end of a particular thermal season were calculated using the following formulas constructed by Gumiński (1948): -
where
(1) is the formula used for rising air temperature in the annual cycle,
(2) is the formula used for falling air temperature in the annual cycle,
t1 = mean air temperature in the month preceding occurrence of threshold temperature,
t2 = mean air temperature in the month following the occurrence of threshold temperature,
tp = value of threshold air temperature, and
x = number of days between the day of threshold air temperature and the 15th day of the preceding month.
Method assumptions:
Monthly mean air temperature falls on the 15th day of the month
Every month has 30 days
Air temperature changes steadily (rise, fall) from month to month
The number of days (x) calculated from both formulas is added to the 15th day of the preceding month. If the value is greater than 15, the real number of days in this month (28, 30, or 31) is used in the process of adding. In this way, the sought date is obtained: it will be the date of the end of a given season or the date of the onset of the next season, depending on which seasonal threshold the air temperature is being counted for.
3. Results
a. Monthly resolution
Looking at spatially averaged SAT changes in the Arctic (65°–90°N) in 1900–2014, it is apparent that the two study periods (1934–50 and 2007–16) represent the warmest phases of the ETCAW and CAW, respectively [see Fig. 1 in Bekryaev et al. (2010) or Fig. 1 in Johannessen et al. (2016)]. This means that the period of greatest warming rate (in historical time or in the present) can be more correctly estimated.
Annual cycles of mean monthly SATs in the analyzed Arctic stations during the ETCAW and CAW periods are presented, alongside their differences, in Table S3 (supplement 1) and Fig. 3. Almost all monthly and seasonal means were greater in the CAW than in the ETCAW, except at the Ilulissat station (Table 2, Fig. 3). The evidently most significant SAT increases between the ETCAW and the CAW occurred in the Pacific region (2.7°C on average) and the Canadian region (1.9°C) (Table 2). Warming in these regions was recorded mainly in the cold half-year, at more than 3° and 2°C, respectively. By contrast, the smallest warming was recorded in the Baffin Bay region (on average only 0.2°C) and occurred only in winter (being of 1.2°C) and summer (0.5°C) (Table 2).

(a),(c),(e),(g),(i),(k) Annual cycles of 10-yr mean monthly air temperatures (TMEAN) in the analyzed Arctic stations during ETCAW (solid red line) and CAW (solid black line) periods and (b),(d),(f),(h),(j),(l) their differences (bars) between ETCAW and CAW periods. Data from CAW were subtracted from ETCAW. Variability of mean monthly air temperature shown using SD calculated from a 30-yr period (1990–2019), with distances of ±1 and ±2 SD from mean presented as dotted and dashed lines, respectively.
Citation: Journal of Climate 35, 7; 10.1175/JCLI-D-21-0162.1

(a),(c),(e),(g),(i),(k) Annual cycles of 10-yr mean monthly air temperatures (TMEAN) in the analyzed Arctic stations during ETCAW (solid red line) and CAW (solid black line) periods and (b),(d),(f),(h),(j),(l) their differences (bars) between ETCAW and CAW periods. Data from CAW were subtracted from ETCAW. Variability of mean monthly air temperature shown using SD calculated from a 30-yr period (1990–2019), with distances of ±1 and ±2 SD from mean presented as dotted and dashed lines, respectively.
Citation: Journal of Climate 35, 7; 10.1175/JCLI-D-21-0162.1
(a),(c),(e),(g),(i),(k) Annual cycles of 10-yr mean monthly air temperatures (TMEAN) in the analyzed Arctic stations during ETCAW (solid red line) and CAW (solid black line) periods and (b),(d),(f),(h),(j),(l) their differences (bars) between ETCAW and CAW periods. Data from CAW were subtracted from ETCAW. Variability of mean monthly air temperature shown using SD calculated from a 30-yr period (1990–2019), with distances of ±1 and ±2 SD from mean presented as dotted and dashed lines, respectively.
Citation: Journal of Climate 35, 7; 10.1175/JCLI-D-21-0162.1
10-yr seasonal and annual (YEAR) means of air temperature (°C) in analyzed Arctic stations in ETCAW and CAW periods, and differences (ETCAW − CAW).


A significant feature of climate classification procedures is a description of the degree of climate continentality/oceanicity. To roughly estimate climate continentality/oceanicity, the ATR is usually used. Still, for scientific purposes, climate continentality/oceanicity needs to be evaluated precisely by calculating indices (see section 2 for more details) that eliminate the influence of latitude on the ATR value.
The highest continentalities and lowest oceanicities of climate in both study periods occurred in the Siberian region (Tiksi) and the Canadian region (Coppermine), at about 80%–85% and 1.3–1.4, respectively. By contrast, the smallest continentality and greatest oceanicity of climate were recorded in the Atlantic region; at Barentsburg, calculated indices are about 30% and 3.2–3.6, respectively (Table 3). The climate continentality index in the Arctic was usually 1%–4% greater during the ETCAW than during the CAW, except in the Russian Arctic.
10-yr averages of annual air temperature range (ATR), thermal climate continentality index (K) according to Ewert’s (1972) formula, and thermal oceanicity index (Oc) according to Marsz’s (1995) formula during ETCAW and CAW periods, and differences (ETCAW − CAW).


b. Daily and subdaily resolution
It is well known that in the process of averaging “raw” meteorological data for increasingly longer periods, important climatic information (including climate change detection) may very often be lost (Robinson et al. 1995). Therefore, as in our previous papers analyzing temperature changes in the Arctic from historical periods to the present (Przybylak and Vízi 2005; Przybylak and Wyszyński 2017; Przybylak et al. 2018), we decided to analyze more precisely the SAT regime for the Arctic using different parameters of daily data [TMEAN (mean monthly air temperatures), TMAX, TMIN, and DTR].
1) Annual cycles
A comparison of 10-yr mean annual cycles of MDAT in both analyzed periods and their differences is shown in Fig. 4. On average, most days in the ETCAW were colder than in the CAW, except at Ilulissat station, where a more or less equal number of positive and negative daily differences is recorded. It is very apparent that areas of greatest climate continentality (the Siberian and Canadian regions) have a less stable range of differences than the places with the smallest continentality (the Atlantic and Pacific regions). There is a clear domination of negative differences in these latter areas, with only minor exceptions (Figs. 4b,d,h). On the other hand, in areas with large climate continentality, many positive differences are recorded, particularly in February and March (see Figs. 4f,j).

(a),(c),(e),(g),(i),(k) 10-yr mean annual cycles of MDAT in ETCAW (solid red line) and CAW (solid black line) periods, and their (b),(d),(f),(h),(j),(l) mean daily differences (bars) between ETCAW and CAW periods. Data from CAW were subtracted from ETCAW. Dotted and dashed lines indicate ±1 SD and ±2 SD, respectively. SDs are calculated on basis of data for 2007–16.
Citation: Journal of Climate 35, 7; 10.1175/JCLI-D-21-0162.1

(a),(c),(e),(g),(i),(k) 10-yr mean annual cycles of MDAT in ETCAW (solid red line) and CAW (solid black line) periods, and their (b),(d),(f),(h),(j),(l) mean daily differences (bars) between ETCAW and CAW periods. Data from CAW were subtracted from ETCAW. Dotted and dashed lines indicate ±1 SD and ±2 SD, respectively. SDs are calculated on basis of data for 2007–16.
Citation: Journal of Climate 35, 7; 10.1175/JCLI-D-21-0162.1
(a),(c),(e),(g),(i),(k) 10-yr mean annual cycles of MDAT in ETCAW (solid red line) and CAW (solid black line) periods, and their (b),(d),(f),(h),(j),(l) mean daily differences (bars) between ETCAW and CAW periods. Data from CAW were subtracted from ETCAW. Dotted and dashed lines indicate ±1 SD and ±2 SD, respectively. SDs are calculated on basis of data for 2007–16.
Citation: Journal of Climate 35, 7; 10.1175/JCLI-D-21-0162.1
2) Day-to-day MDAT variability
Several papers based on climate model simulations argue that high-frequency temperature variability should decrease in a warmer climate (e.g., Houghton et al. 1990, 1992, 1996; Karl et al. 1995, and references cited therein; Mearns et al. 1995; Zwiers and Kharin 1998; Moberg et al. 2000; Screen 2014). In the present paper, intraseasonal variability of MDAT is analyzed using DDTV and SD. Earlier in the article, we proved that the ETCAW was colder than the CAW, and therefore greater MDAT variability is expected in the ETCAW than in the CAW. Analysis of Table 4 and Fig. 5 generally confirms the existence of such a tendency, except in the Baffin Bay region.

(a) Average 10-yr monthly means of DDTV in the Arctic stations during the ETCAW period, and (b) their differences in relation to the CAW period. Bars from left to right denote Barentsburg (dark red), Kanin Nos (light red), Tiksi (orange), Ostrov Vrangel (yellow), Coppermine (dark blue), and Ilulissat light blue), (c),(d) As in (a) and (b), but for SD.
Citation: Journal of Climate 35, 7; 10.1175/JCLI-D-21-0162.1

(a) Average 10-yr monthly means of DDTV in the Arctic stations during the ETCAW period, and (b) their differences in relation to the CAW period. Bars from left to right denote Barentsburg (dark red), Kanin Nos (light red), Tiksi (orange), Ostrov Vrangel (yellow), Coppermine (dark blue), and Ilulissat light blue), (c),(d) As in (a) and (b), but for SD.
Citation: Journal of Climate 35, 7; 10.1175/JCLI-D-21-0162.1
(a) Average 10-yr monthly means of DDTV in the Arctic stations during the ETCAW period, and (b) their differences in relation to the CAW period. Bars from left to right denote Barentsburg (dark red), Kanin Nos (light red), Tiksi (orange), Ostrov Vrangel (yellow), Coppermine (dark blue), and Ilulissat light blue), (c),(d) As in (a) and (b), but for SD.
Citation: Journal of Climate 35, 7; 10.1175/JCLI-D-21-0162.1
Mean seasonal and annual values of DDTV and SD in analyzed Arctic stations in the ETCAW and CAW periods.


In both study periods, DDTV and SD are greatest in winter and smallest in summer (Table 4, Figs. 5a,c). Such a pattern is typical for the Arctic (e.g., Przybylak 2002b; Schweiger et al. 2019). Differences in MDAT variability between the ETCAW and CAW periods, both positive and negative, are greatest from October to April and smallest in July and August (see Fig. 5b,d). It is worth noting that the decrease in MDAT variability in a warmer climate is markedly registered. Its rate usually increases in areas with smaller (greater) climate continentality (oceanicity). Therefore, the most significant changes between ETCAW and CAW periods were observed in stations lying on islands (Barentsburg and Ostrov Vrangel) in the Atlantic and Pacific regions.
The annual cycles of extreme DDTV (Fig. S2, supplement 2) show similar features, as do the mean values of DDTV (Fig. 5). In both analyzed periods, markedly greater extremes of DDTV were recorded in the cold half-year than in the warm half-year. The contrast between them is particularly huge in the areas with a degree of climate continentality (Barentsburg and Ostrov Vrangel) (see Fig. S2). An unclear pattern of changes in extreme DDTV (both rises and drops) between the ETCAW and CAW is generally recorded (see red and blue values in Fig. S2).
In the cold half-year (November to April), extreme DDTV (both drops and rises) usually exceeds 10°C, and very often also 15°C. The highest observed jump in DDTV (of 25.8°C) occurred at Barentsburg from 14 to 15 March 1948 (Table S4, supplement 1). Araźny (2008) also documented similar large jumps in DDTV in the Norwegian Arctic in contemporary times; for example, at Svalbard Airport (approximately 35 km northeast of Barentsburg), MDAT increased from −28.2°C on 3 March 1976 to −0.9°C the next day. Przybylak (1992) found that the greatest DDTV changes in the Arctic are associated with nonperiodic changes in air temperature caused by changes in atmospheric circulation.
In the analyzed case study of extreme DDTV (14–15 March 1948), within 24 h atmospheric pressure fell by about 30 hPa, wind speed increased from 1 to 10 m s−1, and wind direction changed from north to southeast at the grid point near Barentsburg (see Fig. 6). The synoptic situation that mainly explains the described dramatic, highly atypical weather changes in Svalbard is shown as an animation in online supplement 3. A deep cyclone with pressure oscillating around 970 hPa in the center of the Greenland Sea pushed a high pressure system situated over Svalbard eastward, causing the fast advection of very warm (5°–10°C) air masses from the south to reach Svalbard. Therefore, it can be stated that 20CRv3 correctly reproduced this extreme weather pattern, similarly as was earlier documented for other extreme weather cases from lower latitudes (see Fig. 16 in Slivinski et al. 2019, and Fig. 3 in Slivinski et al. 2021).

Courses of (a) sea level pressure (SLP; hPa) and air temperature at 2 m AGL (T2m; °C) and (b) wind speed at 10 m AGL (wspd; m s−1) and wind direction (dir; 0°–360°) in the nearest 20CRv3 (Slivinski et al. 2019) grid point (78.00°N, 14.00°E) to Barentsburg (78.07°N, 14.25°E) from 12 Mar 1948 to 17 Mar 1948.
Citation: Journal of Climate 35, 7; 10.1175/JCLI-D-21-0162.1

Courses of (a) sea level pressure (SLP; hPa) and air temperature at 2 m AGL (T2m; °C) and (b) wind speed at 10 m AGL (wspd; m s−1) and wind direction (dir; 0°–360°) in the nearest 20CRv3 (Slivinski et al. 2019) grid point (78.00°N, 14.00°E) to Barentsburg (78.07°N, 14.25°E) from 12 Mar 1948 to 17 Mar 1948.
Citation: Journal of Climate 35, 7; 10.1175/JCLI-D-21-0162.1
Courses of (a) sea level pressure (SLP; hPa) and air temperature at 2 m AGL (T2m; °C) and (b) wind speed at 10 m AGL (wspd; m s−1) and wind direction (dir; 0°–360°) in the nearest 20CRv3 (Slivinski et al. 2019) grid point (78.00°N, 14.00°E) to Barentsburg (78.07°N, 14.25°E) from 12 Mar 1948 to 17 Mar 1948.
Citation: Journal of Climate 35, 7; 10.1175/JCLI-D-21-0162.1
3) Diurnal temperature range
Diurnal temperature range is a very useful and valuable weather characteristic. For the Arctic its high values in the warm half-year, and particularly in the summer months, are a proxy for more or less stable, sunny weather with calm or weak wind resulting from the occurrence of a high pressure system (anticyclone). Only in this period does a statistically significant negative correlation exist between cloudiness and DTR, although not in all Arctic regions (Przybylak 1997). On the other hand, a slight positive correlation prevailed in winter that was not statistically significant, except for one station located in the Baffin Bay region. This means that, in this season, other factors controlled the magnitude of the DTR. Baranowski (1968), Przybylak (1992), Araźny (2008, 2019), and Araźny et al. (2018) used temperature data from Svalbard to show that the DTR in winter, early spring, and late fall (when the solar radiation is low or not present) was shaped mainly by nonperiodic day-to-day changes in air temperature. Conversely, the latter changes were controlled largely by the advection of warm and humid air masses associated with synoptic-scale cyclones and anticyclones.
The highest values of the DTR in the Arctic during the ETCAW were observed in March and April (except at stations in the Atlantic region), while the lowest occurred almost exclusively in September (Fig. 7). Continental parts of the Arctic have the highest mean annual DTR values (Tiksi and Coppermine: about 7°–8°C), while the most oceanic had the lowest (Barentsburg and Kanin Nos: 4°–5°C) (Table 5). In the mentioned continental parts of the Arctic, mean seasonal values were highest in both analyzed periods in spring (8.5°–9.0°C) and lowest in fall (6.1°–6.6°C). In turn, in the Atlantic region, DTR is highest in winter and lowest in summer (Barentsburg) or shows two maxima, one in summer and one in winter (Kanin Nos) (Table 5).

(a) Average 10-yr monthly means of diurnal temperature range (DTR) in the Arctic stations during the ETCAW period, and (b) their differences in relation to the CAW period. Bars from left to right denote Barentsburg (dark red), Kanin Nos (light red), Tiksi (orange), Ostrov Vrangel (yellow), Coppermine (dark blue), and Ilulissat (light blue).
Citation: Journal of Climate 35, 7; 10.1175/JCLI-D-21-0162.1

(a) Average 10-yr monthly means of diurnal temperature range (DTR) in the Arctic stations during the ETCAW period, and (b) their differences in relation to the CAW period. Bars from left to right denote Barentsburg (dark red), Kanin Nos (light red), Tiksi (orange), Ostrov Vrangel (yellow), Coppermine (dark blue), and Ilulissat (light blue).
Citation: Journal of Climate 35, 7; 10.1175/JCLI-D-21-0162.1
(a) Average 10-yr monthly means of diurnal temperature range (DTR) in the Arctic stations during the ETCAW period, and (b) their differences in relation to the CAW period. Bars from left to right denote Barentsburg (dark red), Kanin Nos (light red), Tiksi (orange), Ostrov Vrangel (yellow), Coppermine (dark blue), and Ilulissat (light blue).
Citation: Journal of Climate 35, 7; 10.1175/JCLI-D-21-0162.1
10-yr seasonal and annual means of DTR (°C) in analyzed Arctic stations in the ETCAW and CAW periods, and differences (ETCAW − CAW). Key: highest (lowest) seasonal means are shown in bold (italic).


A comparison of mean monthly, seasonal, and annual values of DTR between the ETCAW and CAW is shown in Table 5 and Fig. 7. Their analysis reveals that, during the ETCAW, the DTR was usually greater than in recent times, excluding in the most continental parts of the Arctic where, in annual terms, there was no change (Tiksi) or DTR was lower (Coppermine) (Table 5). Such a pattern was also usually observed in all seasons, except summer, when rises in DTRs dominated between the ETCAW and CAW periods (except for the Atlantic region) (Table 5, Fig. 7).
The occurrence of large and extremely large DTR values reaching at least 5°, 10°, or 15°C during the ETCAW period reveal a clear annual cycle similar to that shown for mean values (cf. Figs. 7 and 8). The probability of occurrence of such defined DTRs, particularly those reaching at least 10°C, is definitely greater during the cold half-year (October–April) than the warm half-year. This pattern is evident in stations representing the most oceanic climate, where in summer extremely large DTRs (at least ≥10°C) were observed either not at all (Barentsburg) or very rarely (Kanin Nos, Ilulissat, and Ostrov Vrangel; <10%) (see Table S5, supplement 1).

(a),(c),(e),(g),(i),(k) Annual cycles of relative frequency of occurrence (in %) of characteristic days with DTR of ≥15°C, ≥10°C, and ≥5°C) in the analyzed Arctic stations during the ETCAW period, and (b),(d),(f),(h),(j),(l) their differences. Data from CAW period were subtracted from ETCAW period.
Citation: Journal of Climate 35, 7; 10.1175/JCLI-D-21-0162.1

(a),(c),(e),(g),(i),(k) Annual cycles of relative frequency of occurrence (in %) of characteristic days with DTR of ≥15°C, ≥10°C, and ≥5°C) in the analyzed Arctic stations during the ETCAW period, and (b),(d),(f),(h),(j),(l) their differences. Data from CAW period were subtracted from ETCAW period.
Citation: Journal of Climate 35, 7; 10.1175/JCLI-D-21-0162.1
(a),(c),(e),(g),(i),(k) Annual cycles of relative frequency of occurrence (in %) of characteristic days with DTR of ≥15°C, ≥10°C, and ≥5°C) in the analyzed Arctic stations during the ETCAW period, and (b),(d),(f),(h),(j),(l) their differences. Data from CAW period were subtracted from ETCAW period.
Citation: Journal of Climate 35, 7; 10.1175/JCLI-D-21-0162.1
Both large (≥5°C) and extremely large (≥10°C and ≥15°C) DTRs in areas of oceanic climate mainly showed the greatest average frequency of occurrence in the ETCAW period rather than in present times in all seasons except summer (see Table S5 and Fig. 8). At Ilulissat, large DTRs occurred during the ETCAW period far less frequently than at present (on average 11.2% less), while extremely large DTRs showed opposite relations, except in summer (see Table S5 and Fig. 8). For stations representing continental climate (Tiksi and Coppermine), the picture is more complicated, particularly for extremely large DTRs. Generally, changes are smaller and more stable here than in the oceanic stations; their monthly values usually do not exceed 10% (Fig. 8).
4) Characteristic days
To estimate in more depth the thermal conditions during the ETCAW, the frequencies of occurrence of days that exceed defined thresholds have been calculated (for details, see section 2) and presented in Table S6 (supplement 1), Fig. S4 (supplement 2), and Fig. 9.

Annual courses of differences between number of characteristic days (in %) in the analyzed Arctic stations. Data from CAW period were subtracted from ETCAW data. Order of characteristic days from coldest to warmest is shown from left to right.
Citation: Journal of Climate 35, 7; 10.1175/JCLI-D-21-0162.1

Annual courses of differences between number of characteristic days (in %) in the analyzed Arctic stations. Data from CAW period were subtracted from ETCAW data. Order of characteristic days from coldest to warmest is shown from left to right.
Citation: Journal of Climate 35, 7; 10.1175/JCLI-D-21-0162.1
Annual courses of differences between number of characteristic days (in %) in the analyzed Arctic stations. Data from CAW period were subtracted from ETCAW data. Order of characteristic days from coldest to warmest is shown from left to right.
Citation: Journal of Climate 35, 7; 10.1175/JCLI-D-21-0162.1
Exceptionally cold days (TMAX < −30°C) were not present or were recorded very rarely during the ETCAW period at stations representing oceanic climate (Barentsburg, Kanin Nos, Ostrov Vrangel, and Ilulissat), whereas in the continental parts of the Arctic (Tiksi and Coppermine) they occurred from November to March with a maximum frequency of about 20%–35% in January and February (Fig. S3, supplement 2). Very cold days (TMAX < −20°C) were significantly more common in the latter group of stations (50%–80%), while in the former they were not frequent (less than 10%), except at Ostrov Vrangel (40%–50%) (Table S6, Fig. S3). Such a frequency pattern of highest/lowest values in continental/oceanic parts of the Arctic is also seen for cold days (TMAX < −10°C). On the other hand, frost days were not registered only in July and August (except at Ostrov Vrangel), whereas in the cold half-year (November–April) they were the norm.
Slight frost days (TMAX > 0°C and TMIN ≤ 0°C) and no-frost days (TMIN > 0°C) were recorded in almost all months in stations lying in oceanic parts of the Arctic (except Ostrov Vrangel) (Fig. S3). At the rest of the stations, these days were not usually present from November to March.
TMAX above 5°C was recorded mostly from April/May (except Ilulissat, where they occurred in all months) until October. Their greatest frequency at all stations was recorded in July and August (80%–100%), except at Ostrov Vrangel (40%–50%). Temperatures above 10°C were not very frequent at Svalbard and Ostrov Vrangel (usually much below 10% of days in a month). At the other stations, they occurred with a frequency oscillating between 40% and 70% in July and August (Fig. S3). During the ETCAW exceptionally warm days (TMAX >15°C) were recorded only at stations representing a continental climate and at Ilulissat from May to September, with a maximum usually in July (Fig. S3).
Table S6 and Fig. 9 summarize the results of changes that occurred in the described frequency of characteristic days between the ETCAW and CAW periods. At present, we usually have fewer days of all categories of cold days and more warm days than in the ETCAW, except at Ilulissat, where no significant changes occurred (Table S6, Fig. 9). The direction of changes (decrease/increase) in the number of slight frost days depends on the month, but this category of days was generally more common during the ETCAW than today.
5) Thermal seasons
Baranowski (1968) proposed to use two thresholds of MDAT values for Arctic climate conditions (−2.5° and 2.5°C) to delimit four seasons (for more details, see section 2). His proposition was used in our previous papers studying historical climate in Novaya Zemlya (Przybylak and Wyszyński 2017) and Svalbard (Przybylak et al. 2018). Such an approach allows us to study other valuable aspects of thermal seasons (besides their mean values) such as onset, end, and duration. These additional aspects of thermal seasons are also crucial, as are their mean values, and therefore more research is needed.
Winter is the longest season, exceeding 200 days everywhere except at Kanin Nos (Table 6). However, from the ETCAW to the CAW, the winter became shorter by 2%–8% (8–28 days), except at Ilulissat, where it is presently two days longer (Fig. 10). Summer in all analyzed stations is usually the second-longest season. The summer duration between ETCAW and CAW lengthened by between 12 days (Tiksi) and 63 days (Ostrov Vrangel), except at Ilulissat station, where at present it is two days shorter (Table 6). Changes in duration were smallest for the transitional seasons, except at Ostrov Vrangel in fall (Table 6, Fig. 10). In conclusion, one must say that the ETCAW was characterized by longer winters, springs, and falls, and shorter summers (Fig. 10). The greatest changes in almost all areas were recorded in summer. In contrast, in terms of mean values the most significant changes occurred in winter at most stations. The direction of these changes agrees well with the changes expected due to the increasing warming of the Arctic, which has been particularly great in the twenty-first century.

Changes in duration of thermal seasons between ETCAW and CAW periods (ETCAW − CAW) expressed in percentage of their duration in an entire year.
Citation: Journal of Climate 35, 7; 10.1175/JCLI-D-21-0162.1

Changes in duration of thermal seasons between ETCAW and CAW periods (ETCAW − CAW) expressed in percentage of their duration in an entire year.
Citation: Journal of Climate 35, 7; 10.1175/JCLI-D-21-0162.1
Changes in duration of thermal seasons between ETCAW and CAW periods (ETCAW − CAW) expressed in percentage of their duration in an entire year.
Citation: Journal of Climate 35, 7; 10.1175/JCLI-D-21-0162.1
Dates of onset, end, and duration of each thermal season in analyzed Arctic stations during ETCAW and CAW periods.


4. Discussion
The current warming of the Arctic, which is progressing at a rate comparable to that recorded in the first half of the twentieth century, has undoubtedly inspired many researchers to analyze in detail the climatic conditions (especially SAT) of this second period. They have been interested not only in studying the magnitude of the warming but also its causes. The literature review presented in the introduction shows that although many publications have appeared, many of the thermal characteristics analyzed in this paper (other than the standard monthly, seasonal, or annual averages) have been little studied for the Arctic.
This paper analyzes several important thermal characteristics based on daily and subdaily data (including extremes). All climate models simulating future climate indicate that extremes will play a more significant role than they do now. They are expected to increase in both frequency and magnitude. The models reveal that, in the future, average values of extreme temperature will increase by double the warming in mean temperature. Hegerl et al. (2004) and Portmann et al. (2009) found that changes between means and extremes, minimum and maximum, and upper and lower tail can be significantly different in rate, and even opposite. Zhang et al. (2011) also underline that information about the far tails of the SAT distribution is more relevant to society and the environment than information about temperatures characterizing more frequently occurring aspects of the distribution. For these reasons, the climate models need to simulate both the mean states of the future atmosphere and its extreme conditions (Sillmann et al. 2013). However, the models should first be validated using both present and historical climate data. For validation, however, good-quality data are needed. Such data are particularly crucial for the Arctic, where, even today, the network of stations is very scarce (see section 2).
Latonin et al. (2021) showed that even the multimodel ensemble means in the new generation of high-resolution CMIP6 models do not correctly reproduce the ETCAW (expressed in terms of winter and annual means). The HIRHAM5 regional climate model, too, was not able correctly to reconstruct the annual cycle of daily temperature (overwintering 1930–31) in Franz Joseph Land, particularly in the cold half of the year (Klaus et al. 2018). We can thus expect that temperature extremes for the ETCAW period will also not be captured well, and those discrepancies in observations will be greater than those for mean values. Such a conclusion is strongly supported by results obtained by Sillmann et al. (2013). They showed a large discrepancy between globally averaged extreme temperature indices (hottest and coldest day of a year) selected from four reanalyses (ERA-40, ERA-Interim, NCEP–NCAR, and NCEP–DOE) for the period 1948–2010. Substantial differences between the models (CMIP3 and CMIP5) and the mentioned reanalyses are also evident in their study. The models and reanalyses also disagree with HadEX2 (gridded dataset of indices, updated version; Donat et al. 2013) for the abovementioned temperature extreme indices and the DTR. This discrepancy is particularly evident in Alaska and Greenland (two Arctic regions they analyzed). Sillmann et al. (2013) also suggest that HadEX2 values may be biased by poor observational network coverage in these regions. This weakness of the network is confirmed by the results of climate extremes (data taken from five regional and global reanalyses from the period 2000–16) across the North American Arctic presented recently by Avila-Diaz et al. (2021). Relative to observations, the reanalyses reveal the weakest performance over far northern basins (e.g., the Arctic Ocean and Hudson Bay basins), where observation networks are less dense. However, whatever the shortcomings of the present network, it is nevertheless significantly denser now than it was during the ETCAW.
Investigations conducted by, among others, Brönnimann et al. (2018, 2019) and Slivinski et al. (2019, 2021) have documented the importance of historical data for improving the quality of reanalyses and simulating past climates using climate models, and for identifying factors controlling climate changes. Therefore, there is an urgent need to study climate change in the Arctic in a longer perspective—as Avila-Diaz et al. (2021) described it: “not just in an average sense, but in extreme events that have significant impacts on people and places” (p. 2385). The climate analysis presented in the paper using the many different temperature characteristics (including extremes) for the Arctic for the ETCAW is unique in terms of the characteristics’ comprehensiveness and the number of daily and subdaily data used from meteorological stations. Based on the review of recently published results done by Przybylak and Wyszyński (2020) and the present analysis, we can conclude that good-quality data series from stations should still be used to analyze long-term SAT trends and other temperature characteristics (including extremes) in the Arctic.
As was mentioned in the introduction, there exist very few studies that analyze the climate of the ETCAW using daily and subdaily data. As a result, the comparison of the results presented here and in other studies is limited. Similar to other authors (e.g., Johannessen et al. 2016), we confirmed that the rate and magnitude of warming are greater in the present than in the historical period. From the ETCAW to the CAW, the scale of warming in the Arctic was not spatially homogeneous. The greatest warming occurred in the Pacific and Canadian Arctic regions, as Przybylak (2007) noted. During the ETCAW, the warming was greatest in the Atlantic region (Przybylak 2000a, 2002a). This ETCAW–CAW change in the spatial pattern of warming is connected with the decrease in sea ice extent being greatest in the Pacific sector of the Arctic (Shalina et al. 2020). The smallest changes between the ETCAW and the CAW occurred in the western coast of Greenland, which is in line with results presented for Greenland as a whole by Box et al. (2009). They also stated that the rate of warming in 1919–32 and 1994–2007, directly preceding the warmest parts of the ETCAW and the CAW analyzed in the present paper, was 0.6°C greater in the earlier period than in the later one.
Climate continentality is one of the most important indicators of climate change. Hirschi et al. (2007) found that the most significant changes (decreases) in continentality due to current global warming (1948–2005) occurred in the polar regions, particularly in the Arctic. Furthermore, Stonevicius et al. (2018), analyzing climate continentality (Conrad’s continentality index) in the Northern Hemisphere above 30°N in 1950–2015, confirm this finding. Our results using another precise metric of climate continentality (the Ewert index) confirm this tendency for most of the analyzed Arctic regions. Such changes in the Arctic are in line with tendencies observed in recent decades and centuries in many parts of the world and for the world as a whole (Sadowski 1991; Braganza et al. 2003; Przybylak et al. 2005; Hirschi et al. 2007; Wypych 2010; Stonevicius et al. 2018), while opposite changes recorded in the Russian Arctic (rise in continentality from the ETCAW to the CAW) are atypical. Some local influences and/or circulation changes probably increased the ATR in the CAW relative to the ETCAW. It seems likely that the decline in sea ice extent observed in the Arctic in recent decades led to the decrease in climate continentality. At present, air masses developing over open seas in the Arctic are more frequent than in more cold periods. Thus, the degree of continentality can indirectly tell us about sea ice conditions. It is widely known that sea ice is a significant forcing. Knowledge about its occurrence during the ETCAW is incomplete and not very reliable [for more details, see, e.g., the recently published book Sea Ice in the Arctic (Johannessen et al. 2020)].
Another important temperature characteristic analyzed in the paper is intraseasonal DDTV. This aspect of the climate of the Arctic has not previously been investigated for the ETCAW. We found that most Arctic areas (except the Baffin Bay region) show a decrease in DDTV between the ETCAW and the CAW. The most pronounced reductions occurred in the Atlantic and Pacific regions, where the greatest increases in oceanicity of climate were also recorded. Trends in DDTV in 1951–90 were positive in the Norwegian Arctic and eastern Greenland and negative in the Canadian and Russian Arctic, but everywhere were not statistically significant (Przybylak 2002b). Moberg et al. (2000) identified decreasing trends in winter, spring and fall DDTV for northern Europe in a series of 275 years. Screen (2014) also confirmed that intraseasonal DDTV in the cold season significantly decreased over the mid- to high-latitude Northern Hemisphere in 1979–2013. He explained this by arguing that cold days, being the result of advection of winds from the north, warmed more than warm days (advection from the south), and therefore Arctic amplification reduced subseasonal temperature variance. Furthermore, Collow et al. (2019) found significant decreases in cold-season (November–April) intraseasonal temperature anomaly variability for two areas centered over northern Eurasia and northern North America from 1981–90 to 2005–14. They generally confirm Screen’s (2014) conception. Collow et al. (2019) also conclude that the weakened latitudinal temperature gradient results from decreased sea ice in the Arctic.
Out of all temperature extremes analyzed in the present paper, DTR shows the poorest performance over the Arctic in the reanalyses, including the newest ones used by Avila-Diaz et al. (2021), even for recent data (2000–16). The reanalyses usually tend to underestimate the value of DTR. Good knowledge about DTR changes in stations near the coast is crucial because they are sensitive to changes in their surroundings (e.g., sea ice changes). As a result, they can potentially serve as a proxy to reconstruct sea ice presence/absence. Isaksen et al. (2016) found a close relationship between the land-based SAT at Spitsbergen and local and regional sea ice concentration, particularly during the occurrence of anticyclonic circulation. It seems likely that a highly significant relationship also exists between sea ice and DTR.
DTR data for the ETCAW period (the expedition year 1944/45) in the Arctic are available only for Haudegen station in the north of Nordastlandet Island (Przybylak et al. 2018). An apparent decrease in DTR is recorded between ETCAW and the present time (2011–16 and 2014–17, data from Verlegenhuken and Rijpfjorden stations, respectively). Greater DTR in the ETCAW in comparison to today’s values was also found by us for Barentsburg (the closest station to Haudegen) based on 10-yr periods of observations. Sea ice reductions (more open water) observed recently in the vicinity of the Svalbard archipelago probably led to the decrease in DTR directly or through atmospheric circulation (more vigorous and more frequent advection of warm and humid air from southern latitudes). Based on Chernokulsky and Esau’s (2019) studies, the role of cloudiness in this DTR reduction should be excluded. Observed changes in cloudiness (no change in total cloudiness, but increase in convective and decrease in stratiform clouds) in the western part of the Eurasian Arctic should instead enlarge the DTR. Most continental parts of the Arctic show no change (Tiksi) or even an increase (Coppermine) in DTR between the ETCAW and the CAW. Such behavior agrees with the positive DTR trends observed in some parts of the Canadian Arctic in the periods 1951–90 and 1961–90 (particularly in summer) (see Przybylak 2000b). In the Siberian region, the decrease in DTR in winter and summer is correlated with rising cloudiness over 1936–2012. At the same time, an increase in DTR in spring and autumn is probably related to the rising frequency of convective and falling frequency of stratiform clouds, as Chenokulsky and Esau (2019) document. The spatial pattern of the DTR values in the Arctic during the ETCAW was also similar to the DTR patterns presented by Przybylak (2000b) for 1951–90, based on data from 39 stations.
In this paper, we analyze significantly more temperature extreme indices than are proposed by the Expert Team on Climate Change Detection and Indices (Donat et al. 2013). Some are similar, like DTR and selected characteristic days (e.g., frost days). Only one paper analyzing the occurrence of characteristic days is available for the Arctic for the ETCAW (Przybylak et al. 2018). Data come from the Haudegen station (northern part of Nordaustlandet Island, Svalbard) from the 1944/45 expedition year. Analysis reveals no apparent changes in the occurrence of characteristic days between this year and the present period (the increase in cold days was balanced out by the increase in warm days). We obtained more reliable results for Barentsburg based on 10 years of observations during the ETCAW. The frequencies of all cold (warm) days categories were greater (smaller) during the ETCAW than now. However, Donat et al. (2013) found that the hottest days (defined as a monthly maximum value of daily Tmax) in the 1930s had similarly high values as today. The trend of this index in the Arctic in 1951–2010 was negative, excluding the Canadian Arctic. On the other hand, the coldest nights (monthly minimum value of daily TMIN) revealed significant warming in the entire globe, including the Arctic (Donat et al. 2013).
In climatological studies, the year is usually divided into four seasons (DJF, MAM, etc.). Such a division has a physical meaning for moderate latitudes (it more or less correctly resembles the annual temperature cycle), but not for the Arctic, where the winter is much longer than the other seasons. Moreover, such an arbitrary division of the year omits such useful information as the season’s onset, end, and duration (Baranowski 1968). Such knowledge is essential for, among other things, vegetation, economic activity, and the lives of indigenous people. That is why we utilized the criteria proposed by Baranowski (1968) to delimit four standard seasons. Our findings (i.e., a shortening of winter and lengthening of summer) between the ETCAW and the CAW agree well with similar analyses available for Novaya Zemlya (Przybylak and Wyszyński 2017) and northern Svalbard (Przybylak et al. 2018). However, a rather unexpected result, given that the warming is greater in winter than in summer (Bintanja and van der Linden 2013; Przybylak and Wyszyński 2020), is that the changes in the duration of seasons between ETCAW and CAW were more significant in summer than in winter.
This decade will see us pass 100 years since the beginning of continuous research by thousands of scientists trying to solve some major challenges associated with the ETCAW. One such challenge is to unambiguously establish a reliable explanation for the main mechanisms driving this great warming. In Table S1 and Fig. 1 we show how frequently opinions have different and conflicted among scientists throughout the 90-yr history of research, but particularly in articles published in the last 30 years. By contrast, before this time, the main driving mechanisms were generally agreed upon among scientists. Most of them (97.1%) presented the opinion that the ETCAW was an effect of natural forcings—mainly changes in atmospheric circulation (71.6%) and to a lesser degree in oceanic circulation (20.6%). Only Callendar (1938) connected this warming to an increase in CO2 concentrations (2.9%). This cause is now more commonly proposed (13.9%), but still the ascribing of the early twentieth-century warming to natural forcings predominates (80%). In recent papers, however, scientists have proposed a greater set of forcings. Besides atmospheric (24.5%) and oceanic (12.0%) forcings, there are also usually mentions of internal variability (15.9%), solar radiation or solar activity (12.4%), and volcanic activity (7.0%) (Fig. 1). It seems that this evolution of views between those presented in papers published in 1929–90 and those from 1991 to 2020 is due to the improvement in knowledge of the ETCAW with regard to (i) the Arctic climate system (including the development of new global databases, in particular different types of reanalyses and gridded products), and (ii) forcings (in particular solar and volcanic). In recent decades, the dynamic development of global and regional climate models has also played an important role in this process. These are the most efficient tools commonly now in use for detection and attribution analyses. However, improvements in all the mentioned categories of our knowledge are still necessary.
5. Conclusions and final remarks
The analysis of new aspects of Arctic climate during the ETCAW done based on daily and subdaily data gives new or improved knowledge about the character of climate in that period. From the literature review, we know that the scale of the ETCAW was comparable to the scale of the CAW until the beginning of the twenty-first century. The results we present here also show that, in the last 10–15 years, the warming has become greater than during the ETCAW, particularly in the Pacific and Canadian regions (Table 2).
Our investigations also support the conclusion that, in most study areas, the climate was a little more continental (Table 3) and less stable during the ETCAW than at present (Tables 4 and 5, Figs. 5 and 7). It is essential to underline that the greatest decreases in climate continentality between the ETCAW and the CAW in the Arctic occurred in areas with the smallest degree of climate continentality (Atlantic and Pacific regions), but also where the greatest decreases in sea ice extent have recently occurred. The greater amount of open seawater now than in the ETCAW period is also responsible for the more stable character of weather and climate recently observed in the Arctic, as we documented using the DDTV and DTR indices. Again, changes were the most remarkable in those parts of the Arctic most sensitive to sea ice changes (i.e., the Atlantic and Pacific regions). Our results showed that the reactions to global warming differ between continental and oceanic parts of the Arctic and are sometimes even opposite to expectations. Examples include (i) the increase in climate continentality in the Siberian region and (ii) changes in weather variability between the ETCAW and the CAW being smaller in the Canadian Arctic region than elsewhere, or even absent. In conclusion, we must underline that any Arctic-wide averaging of results for these aspects of climate must be treated with caution; at best, this should be done separately for continental and oceanic parts of the Arctic.
Another important finding is the lengthening of summer from the ETCAW to the CAW being greater than the shortening of winter, mainly in the Pacific and Atlantic regions (Table 6, Fig. 10). This seasonal pattern contrasts with the intensity of warming being greater in winter than in summer in oceanic parts of the Arctic between the two periods and the opposite relation in continental parts (see Table 2).
As with other studies of historical periods in the Arctic, this study has certain limitations due to the scarcity of data. First, it is limited by the discontinuity of the data series, such that no continuous, 10-yr period of the ETCAW is represented by full data for any set of stations. Moreover, no set of 10 non-continuous years for the ETCAW could be found for which there was data for all stations. In light of the impossibility of generating historical datasets that simply do not exist and the pressing need to perform the best analysis that the actual data permit, the selected stations are represented by 10 non-continuous years that do not entirely coincide, but that were selected to maximize overlap. Other measures taken to reduce the negative impact of this problem were also detailed in the study. A second limitation is that the scarcity of data restricts the capacity for any set of stations to be fully representative of their sector. Thus an assessment was made of the geographical range for which each selected station can be deemed to be representative, which varies by season.
Almost 100 papers were reviewed to document the evolution of views about the reasons for the ETCAW (Fig. 1). The markedly dominant opinions/findings (for details see Table S1) presented since 1929 (the first paper) have been that the ETCAW was caused by natural forcings (97.1% and 80.0% in articles published in 1929–90 and 1991–2020, respectively), particularly by atmospheric circulation (71.3% and 24.5%, respectively). It is also clear that there has recently been significantly greater diversity of opinions—probably due to, for example, better knowledge about climate and forcings, newly available datasets, and the use of statistical and modeling tools. As a result, no explanation of the main driving mechanisms responsible for the ETCAW is yet definitive. A workshop devoted to this issue could be one of the best solutions.
We hope that the expanded and improved knowledge about climate during the ETCAW presented in the paper will be useful in research on, for example, 1) state and changes observed in other components of the Arctic climate system from ETCAW to present times, and 2) causes of the early-twentieth-century warming using both statistical and modeling tools. However, there is still an urgent need for more investigations similar to that presented in this paper [i.e., based on data of finer-than-monthly resolution (daily and subdaily)]. To date, however, as our review shows, only a very limited number of such studies are available for the ETCAW period. Data-rescue activity aiming to collect, copy, and digitize “old” data from the entire world (including the Arctic) should be intensified [for details, see Brönnimann et al. (2019)]. Such activity is necessary to study the historical climate and environmental changes in every area of the world. We need more data available in libraries, archives, and private collections that have not yet seen the light of day. More efforts and more scientists are therefore urgently needed to obtain these data for analysis.
Acknowledgments.
The study was carried out as a part of the project entitled “Causes of the early 20th century Arctic warming” funded by the National Science Centre, Poland (Grant 2015/19/B/ST10/02933), and the Research University–Initiative of Excellence: the Emerging Field “Global Environmental Changes”, “Climate Change Research Unit” at Nicolaus Copernicus University in Toruń. The Twentieth Century Reanalysis V3 data were provided by the NOAA/OAR/ESRL PSL, Boulder, Colorado, USA, from their website https://psl.noaa.gov/data/gridded/data.20thC_ReanV3.html. Support for the Twentieth Century Reanalysis Project version 3 dataset is provided by the U.S. Department of Energy, Office of Science Biological and Environmental Research (BER), by the National Oceanic and Atmospheric Administration Climate Program Office, and by the NOAA Physical Sciences Laboratory. The authors declare that they have no conflict of interest.
Author contributions.
Study design by RP. Data collection and selection by PW and RP. Literature review by RP. Statistical analysis and visualization by PW and AA. Interpretation of results by RP, PW, and AA. Preparation of manuscript by RP with contributions from all co-authors.
Data availability statement.
Datasets for this research were derived from the following public domain resources: 1) All-Russia Research Institute of Hydrometeorological Information–World Data Centre (RIHMI-WDC), http://meteo.ru/; 2) The Government of Canada (Environment and Climate Change Canada), https://climate.weather.gc.ca/; 3) Danish Meteorological Institute (DMI), https://www.dmi.dk/publikationer/ [as cited in Cappelen (2020)].
REFERENCES
Ahlmann, H. W., 1948: The present climatic fluctuation. Geogr. J., 112, 165–193, https://doi.org/10.2307/1789696.
Araźny, A., 2008: Bioclimatic Conditions and Their Variability in the Norwegian Arctic for the Period 1971–2000 (in Polish). Wydawnictwo Uniwersytetu Mikołaja Kopernika, 215 pp.
Araźny, A., 2019: Temporal and spatial variability of thermal and humidity stimuli in the Hornsund area (Svalbard). Pol. Polar Res., 40, 29–53, http://dx.doi.org/10.24425/ppr.2019.126346.
Araźny, A., R. Przybylak, P. Wyszyński, T. Wawrzyniak, A. Nawrot, and T. Budzik, 2018: Spatial variations in air temperature and humidity over Hornsund fjord (Spitsbergen) from 1 July 2014 to 30 June 2015. Geogr. Ann., 100A, 27–43, https://doi.org/10.1080/04353676.2017.1368832.
Araźny, A., P. Wyszyński, and R. Przybylak, 2019: A comparison of bioclimatic conditions on Franz Josef Land (the Arctic) between the turn of the nineteenth to twentieth century and present day. Theor. Appl. Climatol., 137, 2623–2638, https://doi.org/10.1007/s00704-018-02763-y.
Avila-Diaz, A., D. H. Bromwich, A. B. Wilson, F. Justino, and S.-H. Wang, 2021: Climate extremes across the North American Arctic in modern reanalyses. J. Climate, 34, 2385–2410, https://doi.org/10.1175/JCLI-D-20-0093.1.
Baranowski, S., 1968: Thermals of the Periglacial Tundra SW Spitsbergen (in Polish). Acta Universitatis Wratislaviensis Series, Vol. 68, Państwowe Wydawnictwo Naukowe, 77 pp.
Bekryaev, R. V., I. Polyakov, and V. A. Alexeev, 2010: Role of polar amplification in long-term surface air temperature variations and modern Arctic warming. J. Climate, 23, 3888–3906, https://doi.org/10.1175/2010JCLI3297.1.
Bengtsson, L., V. Semenov, and O. Johannessen, 2004: The early twentieth-century warming in the Arctic—A possible mechanism. J. Climate, 17, 4045–4057, https://doi.org/10.1175/1520-0442(2004)017<4045:TETWIT>2.0.CO;2.
Bintanja, R., and E. C. van der Linden, 2013: The changing seasonal climate in the Arctic. Sci. Rep., 3, 1556, https://doi.org/10.1038/srep01556.
Bokuchava, D. D., and V. A. Semenov, 2021: Mechanisms of the early 20th century warming in the Arctic. Earth-Sci. Rev., 222, 103820, https://doi.org/10.1016/j.earscirev.2021.103820.
Box, J. E., L. Yang, D. H. Bromwich, and L.-S. Bai, 2009: Greenland ice sheet surface air temperature variability: 1840–2007. J. Climate, 22, 4029–4049, https://doi.org/10.1175/2009JCLI2816.1.
Braganza, K., D. J. Karoly, A. C. Hirst, M. E. Mann, P. Stott, R. J. Stouffer, and S. F. B. Tett, 2003: Simple indices of global climate variability and change: Part I—Variability and correlation structure. Climate Dyn., 20, 491–502, https://doi.org/10.1007/s00382-002-0286-0.
Brennan, M. K., G. J. Hakim, and E. Blanchard-Wrigglesworth, 2020: Arctic sea-ice variability during the instrumental era. Geophys. Res. Lett., 47, e2019GL086843, https://doi.org/10.1029/2019GL086843.
Brönnimann, S., 2015: Climate Changes since 1700. Springer, 360 pp., https://doi.org/10.1007/978-3-319-19042-6.
Brönnimann, S., and Coauthors, 2018: A roadmap to climate data rescue services. Geosci. Data J., 5, 28–39, https://doi.org/10.1002/gdj3.56.
Brönnimann, S., and Coauthors, 2019: Unlocking pre-1850 instrumental meteorological records: A global inventory. Bull. Amer. Meteor. Soc., 100, ES389–ES413, https://doi.org/10.1175/BAMS-D-19-0040.1.
Callendar, G. S., 1938: The artificial production of carbon dioxide and its influence on temperatures. Quart. J. Roy. Meteor. Soc., 64, 223–227, https://doi.org/10.1002/qj.49706427503.
Cappelen, J., Ed., 2020: Greenland–DMI historical climate data collection 1784–2019. DMI Rep. 20-04, Danish Meteorological Institute, 119 pp.
Chernokulsky, A., and I. Esau, 2019: Cloud cover and cloud types in the Eurasian Arctic in 1936–2012. Int. J. Climatol., 39, 5771–5790, https://doi.org/10.1002/joc.6187.
Collow, T. W., W. Wang, and A. Kumar, 2019: Reduction in northern midlatitude 2-m temperature variability due to Arctic sea ice loss. J. Climate, 32, 5021–5035, https://doi.org/10.1175/JCLI-D-18-0692.1.
Delworth, T. L., and T. R. Knutson, 2000: Simulation of early 20th century global warming. Science, 287, 2246–2250, https://doi.org/10.1126/science.287.5461.2246.
Donat, M. G., and Coauthors, 2013: Updated analyses of temperature and precipitation extreme indices since the beginning of the twentieth century: The HadEX2 dataset. J. Geophys. Res. Atmos., 118, 2098–2118, https://doi.org/10.1002/jgrd.50150.
Ellsaesser, H. W., M. C. MacCracken, J. J. Walton, and S. L. Grotch, 1986: Global climatic trends as revealed by the recorded data. Rev. Geophys., 24, 745–792, https://doi.org/10.1029/RG024i004p00745.
Eserkepova, I. B., K. M. Lugina, I. A. Speranskaya, P. D. Kagan, L. P. Smirnova, and V. V. Subbotin, 1982: Ispolzovanije novych metodov v monitoringe termićeskogo reźima Arktiki. Issledovanija Arktiki, Antarktiki i Mirovogo okeana, Gidrometeoizdat, 56–64.
Ewert, A., 1972: About calculating thermic continentality of climate (in Polish). Prz. Geogr., 44, 273–328.
Fyfe, J. C., K. von Salzen, N. P. Gillett, V. K. Arora, G. M. Flato, and J. R. McConnell, 2013: One hundred years of Arctic surface temperature variation due to anthropogenic influence. Sci. Rep., 3, 2645, https://doi.org/10.1038/srep02645.
Grant, A. N., S. Brönnimann, T. Ewen, T. Griesser, and A. Stickler, 2009: The early twentieth century warm period in the European Arctic. Meteor. Z., 18, 425–432, https://doi.org/10.1127/0941-2948/2009/0391.
Gumiński, R., 1948: Concerning the establishment of agricultural climatic regions in Poland (in Polish). Przegląd Meteor. Hydrol., 1 (1), 7–20.
Hanssen-Bauer, I., and E. Forland, 1998: Long-term trends in precipitation and temperature in the Norwegian Arctic: Can they be explained by changes in atmospheric circulation patterns? Climate Res., 10, 143–153, https://doi.org/10.3354/cr010143.
Hegerl, G. C., F. W. Zwiers, P. A. Stott, and V. V. Kharin, 2004: Detectability of anthropogenic changes in annual temperature and precipitation extremes. J. Climate, 17, 3683–3700, https://doi.org/10.1175/1520-0442(2004)017<3683:DOACIA>2.0.CO;2.
Hegerl, G. C., S. Brönnimann, A. Schurer, and T. Cowan, 2018: The early 20th century warming: Anomalies, causes, and consequences. Wiley Interdiscip. Rev.: Climate Change, 9, e522, https://doi.org/10.1002/wcc.522.
Hersbach, H., and Coauthors, 2020: The ERA5 global reanalysis. Quart. J. Roy. Meteor. Soc., 146, 1999–2049, https://doi.org/10.1002/qj.3803.
Hirschi, J., B. Sinha, and S. A. Josey, 2007: Global warming and changes of continentality since 1948. Weather, 62, 215–221, https://doi.org/10.1002/wea.88.
Houghton, J. T., G. J. Jenkins, and J. J. Ephraums, Eds., 1990: Climate Change: The IPCC Scientific Assessment. Cambridge University Press, 365 pp.
Houghton, J. T., B. A. Callander, and S. K. Varney, Eds., 1992: Climate Change 1992: The Supplementary Report to the IPCC Scientific Assessment. Cambridge University Press, 200 pp.
Houghton, J. T., L. G. Meila Filho, B. A. Callander, N. Harris, A. Kattenberg, and K. Maskell, Eds., 1996: Climate Change 1995: The Science of Climate Change. Cambridge University Press, 572 pp.
Isaksen, K., Ø. Nordli, E. J. Førland, E. Łupikasza, S. Eastwood, and T. Niedźwiedź, 2016: Recent warming on Spitsbergen—Influence of atmospheric circulation and sea ice cover. J. Geophys. Res. Atmos., 121, 11 913–11 931, https://doi.org/10.1002/2016JD025606.
Jensen, A. S., 1939: Concerning a change of climate during recent decades in the Arctic and subarctic regions, from Greenland in the West to Eurasia in the East, and contemporary biological and geophysical changes. Det Kgl. Danske Videnskabernes Selskab. Biologiske Medd, Vol. XIV, Ejnar Munksgaard, 76 pp.
Johannessen, O. M., S. L. Kuzmina, L. P. Bobylev, and M. W. Miles, 2016: Surface air temperature variability and trends in the Arctic: New amplification assessment and regionalisation. Tellus, 68A, 28234, https://doi.org/10.3402/tellusa.v68.28234.
Johannessen, O. M., L. P. Bobylev, E. V. Shalina, and S. Sandven, Eds., 2020: Sea Ice in the Arctic: Past, Present and Future. Springer, 579 pp.
Jones, P. D., and A. Moberg, 2003: Hemispheric and large-scale surface air temperature variations: An extensive revision and an update to 2001. J. Climate, 16, 206–223, https://doi.org/10.1175/1520-0442(2003)016<0206:HALSSA>2.0.CO;2.
Jones, P. D., M. New, D. E. Parker, S. Martin, and I. G. Rigor, 1999: Surface air temperature and its changes over the past 150 years. Rev. Geophys., 37, 173–199, https://doi.org/10.1029/1999RG900002.
Karl, T. R., R. W. Knight, and N. Plummer, 1995: Trends in high-frequency climate variability in the twentieth century. Nature, 377, 217–220, https://doi.org/10.1038/377217a0.
Klaus, D., P. Wyszyński, K. Dethloff, R. Przybylak, and A. Rinke, 2018: Evaluation of 20CR reanalysis data and model results based on historical (1930–1940) observations from Franz Josef Land. Pol. Polar Res., 36, 225–254, https://doi.org/10.24425/118747.
Kuzmina, S. I., O. M. Johannessen, L. Negtsson, O. G. Aniskina, and L. P. Bobylev, 2008: High northern latitude surface air temperature: Comparison of existing data and creation of a new gridded data set 1900–2000. Tellus, 60A, 289–304, https://doi.org/10.1111/j.1600-0870.2008.00303.x.
Latonin, M.M., I.L. Bashmachnikov, L.P. Bobylev, and R. Davy, 2021: Multimodel ensemble mean of global climate models fails to reproduce early twentieth century Arctic warming. Polar Sci., 30, 100677, https://doi.org/10.1016/j.polar.2021.100677.
Łupikasza, E. B., and T. Niedźwiedź, 2019: The influence of mesoscale atmospheric circulation on Spitsbergen air temperature in periods of Arctic warming and cooling. J. Geophys. Res. Atmos., 124, 5233–5250, https://doi.org/10.1029/2018JD029443.
Marsz, A., 1995: Index of Thermic Oceanity as a Measure of Climatic Works in Ocean–Atmosphere–Continents System (in Polish). Wyższa Szkoła Morska, Wydział Nawigacyjny, 110 pp.
Mearns, L. O., F. Giorgi, L. McDaniel, and C. Shields, 1995: Analysis of variability and diurnal range of daily temperature in a nested regional climate model: Comparison with observations and doubled CO2 results. Climate Dyn., 11, 193–209, https://doi.org/10.1007/BF00215007.
Moberg, A., and Coauthors, 2000: Day-to-day variability trends in 160- to 275-year-long European instrumental records. J. Geophys. Res., 105, 22 849–22 868, https://doi.org/10.1029/2000JD900300.
National Academies of Sciences, Engineering, and Medicine, 2016: Attribution of Extreme Weather Events in the Context of Climate Change. The National Academies Press., 93 pp., https://doi.org/10.17226/21852.
Nordli, Ø., P. Wyszyński, H. M. Gjelten, K. Isaksen, E. Łupikasza, T. Niedźwiedź, and R. Przybylak, 2020: Revisiting the extended Svalbard Airport monthly temperature series, and the compiled corresponding daily series 1898–2018. Polar Res., 39, 3614, https://doi.org/10.33265/polar.v39.3614.
Polyakov, I. V., R. V. Bekryaev, G. V. Alekseev, U. S. Bhatt, R. L. Colony, M. A. Johnson, A. P. Makshtas, and D. Walsh, 2003: Variability and trends of air temperature and pressure in the maritime Arctic, 1875–2000. J. Climate, 16, 2067–2077, https://doi.org/10.1175/1520-0442(2003)016<2067:VATOAT>2.0.CO;2.
Portmann, R. W., S. Solomon, and G. C. Hegerl, 2009: Spatial and seasonal patterns in climate change, temperatures, and precipitation across the United States. Proc. Natl. Acad. Sci. USA, 106, 7324–7329, https://doi.org/10.1073/pnas.0808533106.
Przybylak, R., 1992: Thermic and humidity relations against the background of circulations conditions in Hornsund (Spitsbergen) in the period 1978–1983 (in Polish). Dokumentacja Geograficzna, 2, 1–105.
Przybylak, R., 1996: Variability of Air Temperature and Precipitation over the Period of Instrumental Observations in the Arctic (in Polish). Uniwersytet Mikołaja Kopernika, Rozprawy, 280 pp.
Przybylak, R., 1997: Spatial and temporal changes in extreme air temperatures in the Arctic over the period 1951–1990. Int. J. Climatol., 17, 615–634, https://doi.org/10.1002/(SICI)1097-0088(199705)17:6<615::AID-JOC136>3.0.CO;2-0.
Przybylak, R., 2000a: Temporal and spatial variation of air temperature over the period of instrumental observations in the Arctic. Int. J. Climatol., 20, 587–614, https://doi.org/10.1002/(SICI)1097-0088(200005)20:6<587::AID-JOC480>3.0.CO;2-H.
Przybylak, R., 2000b: Diurnal temperature range in the Arctic and its relation to hemispheric and Arctic circulation patterns. Int. J. Climatol., 20, 231–253, https://doi.org/10.1002/(SICI)1097-0088(20000315)20:3<231::AID-JOC468>3.0.CO;2-U.
Przybylak, R., 2002a: Variability of Air Temperature and Atmospheric Precipitation during a Period of Instrumental Observation in the Arctic. Kluwer Academic, 330 pp.
Przybylak, R., 2002b: Changes in seasonal and annual high-frequency air temperature variability in the Arctic from 1951 to 1990. Int. J. Climatol., 22, 1017–1032, https://doi.org/10.1002/joc.793.
Przybylak, R., 2003: The Climate of the Arctic. Kluwer Academic, 270 pp.
Przybylak, R., 2007: Recent air-temperature changes in the Arctic. Ann. Glaciol., 46, 316–324, https://doi.org/10.3189/172756407782871666.
Przybylak, R., 2016: The Climate of the Arctic. 2nd ed. Springer, 287 pp.
Przybylak, R., and Z. Vízi, 2005: Air temperature changes in the Canadian Arctic from the early instrumental period to modern times. Int. J. Climatol., 25, 1507–1522, https://doi.org/10.1002/joc.1213.
Przybylak, R., and P. Wyszyński, 2017: Air temperature in Novaya Zemlya Archipelago and Vaygach Island from 1832 to 1920 in the light of early instrumental data. Int. J. Climatol., 37, 3491–3508, https://doi.org/10.1002/joc.4934.
Przybylak, R., and P. Wyszyński, 2020: Air temperature changes in the Arctic in the period 1951–2015 in the light of observational and reanalysis data. Theor. Appl. Climatol., 139, 75–94, https://doi.org/10.1007/s00704-019-02952-3.
Przybylak, R., and Coauthors, 2005: Temperature changes in Poland from the 16th to the 20th centuries. Int. J. Climatol., 25, 773–791, https://doi.org/10.1002/joc.1149.
Przybylak, R., and Coauthors, 2014: Spatial distribution of air temperature on Svalbard during 1 year with campaign measurements. Int. J. Climatol., 34, 3702–3719, https://doi.org/10.1002/joc.3937.
Przybylak, R., P. Wyszyński, and M. Woźniak, 2018: Air temperature conditions in northern Nordaustlandet (NE Svalbard) at the end of World War II. Int. J. Climatol., 38, 2775–2791, https://doi.org/10.1002/joc.5459.
Przybylak, R., P. N. Svyashchennikov, J. Uscka-Kowalkowska, and P. Wyszyński, 2021: Solar radiation in the Arctic during the early twentieth century warming (1921–1950), presenting a compilation of newly available data. J. Climate, 33, 21–37, https://doi.org/10.1175/JCLI-D-20-0257.1.
Rigor, I. G., R. L. Colony, and S. Martin, 2000: Variations in surface air temperature observations in the Arctic, 1979–1997. J. Climate, 13, 896–914, https://doi.org/10.1175/1520-0442(2000)013<0896:VISATO>2.0.CO;2.
Robinson, D. A., D. J. Leathers, M. A. Palecki, and K. F. Dewey, 1995: Some observations on climate variability as seen in daily temperature structure. Atmos. Res., 37, 119–131, https://doi.org/10.1016/0169-8095(94)00073-M.
Sadowski, M., 1991: Variability of extreme climatic events in central Europe since the 13th century. Z. Meteor., 41, 350–356.
Schweiger, A. J., K. R. Wood, and J. Zhang, 2019: Arctic sea ice volume variability over 1901–2010: A model-based reconstruction. J. Climate, 32, 4731–4752, https://doi.org/10.1175/JCLI-D-19-0008.1.
Screen, J. A., 2014: Arctic amplification decreases temperature variance in northern mid- to high-latitudes. Nat. Climate Change, 4, 577–582, https://doi.org/10.1038/nclimate2268.
Shalina, E. V., O. M. Johannessen, and S. Sandven, 2020: Changes in Arctic sea ice cover in the twentieth and twenty-first centuries. Sea Ice in the Arctic: Past, Present and Future, Johannessen et al., Eds., Springer, 93–166.
Sillmann, J., V. V. Kharin, X. Zhang, F. W. Zwiers, and D. Bronaugh, 2013: Climate extremes indices in the CMIP5 multimodel ensemble: Part 1. Model evaluation in the present climate. J. Geophys. Res. Atmos., 118, 1716–1733, https://doi.org/10.1002/jgrd.50203.
Slivinski, L. C., and Coauthors, 2019: Towards a more reliable historical reanalysis: Improvements for version 3 of the Twentieth Century Reanalysis System. Quart. J. Roy. Meteor. Soc., 145, 2876–2908, https://doi.org/10.1002/qj.3598.
Slivinski, L. C., and Coauthors, 2021: An evaluation of the performance of the Twentieth Century Reanalysis version 3. J. Climate, 34, 1417–1438, https://doi.org/10.1175/JCLI-D-20-0505.1.
Smirnova, L. P., and V. V. Subbotin, 1983. Prostranstvenno-vremiennaja struktura termićeskich polej Arktiki. Obektivnaja otsenka meteorologićeskich informacji pri obespećenii poletov vozdusnych sudov, OŁAGA, 106–117.
Stonevicius, E., G. Stankunavicius, and E. Rimkus, 2018: Continentality and oceanity in the mid and high latitudes of the Northern Hemisphere and their links to atmospheric circulation. Adv. Meteor., 2018, 5746191, https://doi.org/10.1155/2018/5746191.
Subbotin, V. V., 1983: Ispol’zovanije regresionnogo analiza dla polućenija nepreryvnych rjadov meteorologićeskich elementov po nabljudenijam drejfujuśćich stancij. Trudy AANII, 381, 111–121.
Suo, L., O. H. Ottera, M. Bentsen, Y. Gao, and O. M. Johannessen, 2013: External forcing of the early 20th century Arctic warming. Tellus, 65A, 20578, https://doi.org/10.3402/tellusa.v65i0.20578.
Svendsen, L., N. Keenlyside, I. Bethke, Y. Gao, and N. E. Omrani, 2018: Pacific contribution to the early twentieth-century warming in the Arctic. Nat. Climate Change, 8, 793–797, https://doi.org/10.1038/s41558-018-0247-1.
Treshnikov, A. F., Ed., 1985: Atlas Arktiki. Glavnoye Upravlenye Geodeziy i Kartografiy, 204 pp.
Veryard, H. G., 1963: A review of studies on climatic fluctuations during the period of the meteorological record. Changes of Climate: Proceedings of the Rome Symposium Organised by UNESCO and WMO, UNESCO, 3–15.
Wállen, C. C., 1984: Present Century Climate Fluctuations in the Northern Hemisphere and their Impact. World Climate Impact Studies Programme WMO/TD-9, 85 pp.
Wegmann, M., S. Brönnimann, and G. P. Compo, 2017: Tropospheric circulation during the early twentieth century Arctic warming. Climate Dyn., 48, 2405–2418, https://doi.org/10.1007/s00382-016-3212-6.
Wood, K. R., and J. E. Overland, 2010: Early 20th century Arctic warming in retrospect. Int. J. Climatol., 30, 1269–1279, https://doi.org/10.1002/joc.1973.
Wood, K. R., J. E. Overland, T. Jónsson, and B. V. Smoliak, 2010: Air temperature variations on the Atlantic–Arctic boundary since 1802. Geophys. Res. Lett., 37, L17708, https://doi.org/10.1029/2010GL044176.
Wypych, A., 2010: Variability of the European climate on the basis of differentiation of indicators of continentalism. The Polish Climate in the European Context: An Historical Overview, R. Przybylak et al., Eds., Springer, 473–484.
Xiao, H., F. Zhang, L. Miao, X. S. Liang, K. Wu, and R. Liu, 2020: Long-term trends in Arctic surface temperature and potential causality over the last 100 years. Climate Dyn., 55, 1443–1456, https://doi.org/10.1007/s00382-020-05330-2.
Yamanouchi, T., 2011: Early 20th century warming in the Arctic: A review. Polar Sci., 5, 53–71, https://doi.org/10.1016/j.polar.2010.10.002.
Zhang, X., L. Alexander, G. C. Hegerl, P. Jones, A. K. Tank, T. C. Peterson, B. Trewin, and F. W. Zwiers, 2011: Indices for monitoring changes in extremes based on daily temperature and precipitation data. Wiley Interdiscip. Rev.: Climate Change, 2, 851–870, https://doi.org/10.1002/wcc.147.