Analysis of Temperature Change in Uzbekistan and the Regional Atmospheric Circulation of Middle Asia during 1961–2016

Climate change and shrinking of the Aral Sea have significantly affected the region’s temperature variations. Observed interannual changes in Uzbekistan’s air temperature compared to the duration of synoptic weather types (SWT) in Middle Asia were analyzed. Nonparametric Mann–Kendall statistical test and climate trends coefficients were used to identify trend characteristics of observed temperature from 1961–2016 to the baseline period of 1961–1990. The results showed increasing temperature trends average to 1 °C in warm and cold half years over Uzbekistan. The 1991–2016 decadal temperature trend ranged from 0.25 °C/decade in the northwest to 0.52 °C/decade in the center, especially pronounced in the oasis and Aral Sea zones. There were also significant changes in the structure of regional SWT. The main difference in the structure of SWT in Middle Asia relative to the baseline period was expressed in a decrease of cold mass invasion duration from 113.4 to 76.1 days and an increase in low-gradient baric field duration from 65.8 to 134.6 days. The process of anthropogenic warming, which began in Uzbekistan in the 1960s of the twentieth century, has accelerated from the mid-1970s with a higher mean annual air temperature than the baseline period’s climate normals (1961–1990) and is associated with changes in the regional SWT over Middle Asia.


Introduction
Analysis of the global meteorological observations has shown that during 1951-2003, air temperatures increased all year round, and precipitation indices tended to change to wetter conditions over most landscapes [1]. The Intergovernmental Panel on Climate Change's Fifth Assessment Report AR5 [2] has pointed to the rising trend in the observed temperatures and marked Central Asia as one of the significant impact zones, with one of the highest average increase in temperatures. Extreme values of air temperature and precipitation, coupled with increasing anthropogenic pressure on natural resources, have dramatically impacted the climate of the whole area of Middle Asia in the second half of the twentieth century. Giese et al. [3] concluded that the warming trend, on average of almost 2-3 • C, has been especially pronounced in the plains and valleys, while less so in mountainous regions for the long-term temperature trends in Middle Asia from 1950 to 2000. Several researchers have pointed out an upward temperature trend in Middle Asia after 1990 [3][4][5]; however, due to a shortage of available observed data, such analysis has been limited in Uzbekistan after the 2000s. This period also a.
to analyze air temperature trends of the current period 1991-2016 compared to the climate normals of the baseline period 1961-1990 and define spatial and temporal changes to the air temperature in Uzbekistan. b.
to assess changes in the regional SWT of Middle Asia over Uzbekistan and analyze climatic variations and trends connected to the atmospheric circulations as one of the main factors of climate change in Uzbekistan.
The rest of this paper is structured as follows. Section 2 describes the study area, data processing, and methodology for analysis; the SWT classification used by the UzHydromet and their description are also included. Results are presented in Section 3 for the air temperature and changes in regional synoptic processes that occurred from 1961 to 2016. Current observed climate change analysis in Uzbekistan is limited in the literature; here, we provide this information to the international community. The discussion and summary are given in Sections 4 and 5, respectively.

Study Area and Data Collection
The study area covered the whole territory of Uzbekistan, an area of 448,978 km 2 (37 • 13 -45 • 36 N, 56 • 00 -73 • 10 E) located within the Irano-Turannian lowlands of the Central Asian region. Uzbekistan is a double landlocked country that is extremely diverse in physical and geographical conditions. Most of the country is arid and semi-arid plains (80%) that border mountains (20%) from the south, southeast, and east [16].
Uzbekistan is divided into 14 climatic regions by meteorological fluctuations (Figure 1) [17].  Uzbekistan is limited in the literature; here, we provide this information to the international community. The discussion and summary are given in Sections 4 and 5, respectively.

Study Area and Data Collection
The study area covered the whole territory of Uzbekistan, an area of 448,978 km 2 (37°13′-45°36′ N, 56°00′-73°10′ E) located within the Irano-Turannian lowlands of the Central Asian region. Uzbekistan is a double landlocked country that is extremely diverse in physical and geographical conditions. Most of the country is arid and semi-arid plains (80%) that border mountains (20%) from the south, southeast, and east [16].
Uzbekistan is divided into 14 climatic regions by meteorological fluctuations (Figure 1) [17]. The Ustyurt, Northern Kyzylkum, Central Kyzylkum, and Southern Kyzylkum regions are characterized as desert zones. Eight regions-Aralian, Lower Amu Darya, Zeravshan, Southeastern Kyzylkum, Southern Upper Amu Darya, Jizzakh, Tashkent, and Ferghana-comprise oasis zones. The Northern Upper Amu Darya and Western Tien Shan regions are characterized as foothills and mountainous zones. There are 86 meteorological observation stations currently operating in Uzbekistan, and air temperature data from 69 stations were used for this analysis. Station data were collected from 1961 to 2016, except for Oygaing and Termez stations, where observation started in 1963 and 1964, respectively. The most representative meteorological stations located evenly across the territory of the republic for each of the 14 regions were determined and chosen for the analysis. The choice of stations was based on the analysis of the duration, continuity, and homogeneity of the observation series. These were major stations with the longest observation record and showed relatively common climatic conditions to the representative zone.  The daily data from the UzHydromet published handbook "Calendar of Synoptic Weather Types of Middle Asia" [18,19] were used to assess changes in the statistical characteristics of SWT duration in Middle Asia. Statistical characteristics of SWT duration for the baseline climate period were taken from Inogamova et al. [17]. All datasets were manually processed and prepared from the historical catalogs and checked for observation data continuity and homogeneity for the analysis under The World Meteorological Organization (WMO) guidelines.

Synoptic Weather Types Classification
The regional atmospheric circulation and the SWT classification was first suggested by Giorgio and Bugaev [20]. This classification was further broadened by Bugaev et al. [21] and Sarymsakov et al. [22]. Bugaev, founder of the Middle Asian Tashkent Institute of Weather Forecasters and World Weather Watch under WMO, defined the SWT in Middle Asia into 11 major types in the monograph "Synoptic processes of Middle Asia" [15]. In the early 1960s, Inagamova et al. [23] updated the classification by introducing additional SWT over Middle Asia, increasing them to 15. Nowadays, the following types are considered [24]: 1-South Caspian cyclone, 2-Murgab cyclone, 3-Upper Amu Darya cyclone, 4-Wide outflow of warm air, 5-North-west cold wave invasion, 6-Northern cold wave invasion, 7-Wave activity on the cold front, 8-Slow moving cyclone over Middle Asia, 9-Southwest periphery of anticyclone, 9a-Southeast periphery of anticyclone, 9b-South periphery of anticyclone, 10-Western invasion, 11-Summer thermal depression, 12-Low gradient field of high pressure, 13-Low gradient field of low pressure, 14-Western cyclone, 15-Diving cyclone. Although there are 17 main SWTs, some of them are rare and occur once in several years. The main SWT of Middle Asia is shown in Figure 2 for warm and cold half years; a short description of all processes is given in Appendix A , Table A1. The daily data from the UzHydromet published handbook "Calendar of Synoptic Weather Types of Middle Asia" [18,19] were used to assess changes in the statistical characteristics of SWT duration in Middle Asia. Statistical characteristics of SWT duration for the baseline climate period were taken from Inogamova et al. [17]. All datasets were manually processed and prepared from the historical catalogs and checked for observation data continuity and homogeneity for the analysis under The World Meteorological Organization (WMO) guidelines.

Synoptic Weather Types Classification
The regional atmospheric circulation and the SWT classification was first suggested by Giorgio and Bugaev [20]. This classification was further broadened by Bugaev et al. [21] and Sarymsakov et al. [22]. Bugaev, founder of the Middle Asian Tashkent Institute of Weather Forecasters and World Weather Watch under WMO, defined the SWT in Middle Asia into 11 major types in the monograph "Synoptic processes of Middle Asia" [15]. In the early 1960s, Inagamova et al. [23] updated the classification by introducing additional SWT over Middle Asia, increasing them to 15. Nowadays, the following types are considered [24]: 1-South Caspian cyclone, 2-Murgab cyclone, 3-Upper Amu Darya cyclone, 4-Wide outflow of warm air, 5-North-west cold wave invasion, 6-Northern cold wave invasion, 7-Wave activity on the cold front, 8-Slow moving cyclone over Middle Asia, 9-Southwest periphery of anticyclone, 9a-Southeast periphery of anticyclone, 9b-South periphery of anticyclone, 10-Western invasion, 11-Summer thermal depression, 12-Low gradient field of high pressure, 13-Low gradient field of low pressure, 14-Western cyclone, 15-Diving cyclone. Although there are 17 main SWTs, some of them are rare and occur once in several years. The main SWT of Middle Asia is shown in Figure 2 for warm and cold half years; a short description of all processes is given in Appendix A, Table A1. The basic principle of the classification is the position of the upper-level frontal zones over the Northern Hemisphere, which determines the development of synoptic processes; it is a qualitative conclusion from the hydrodynamic theory of pressure change in the analysis of the synoptic position transformations. It falls into a statistical-stochastic approach to the study of the dynamic climatology of Middle Asia. The inclusion of the orography helps distinguish the representation of the main types of circulation from all their variety as one of the most critical factors that create regional features of synoptic processes in Middle Asia.
Such generalization in the classification of the synoptic processes does not fully represent all complexity of atmospheric flows and details of the synoptic processes development; to a certain extent, it is a simplified scheme. Of course, this highlights the subjectivity of classification, especially in controversial cases with deviations from one type to another that is smothered by simplified classification. On the other hand, Uzbekistan's orography variation from the east to west, with The width of the circulation type is based on occurrence frequency during the year. Adapted from Inogamova et al. [23,24].
The basic principle of the classification is the position of the upper-level frontal zones over the Northern Hemisphere, which determines the development of synoptic processes; it is a qualitative conclusion from the hydrodynamic theory of pressure change in the analysis of the synoptic position transformations. It falls into a statistical-stochastic approach to the study of the dynamic climatology of Middle Asia. The inclusion of the orography helps distinguish the representation of the main types of circulation from all their variety as one of the most critical factors that create regional features of synoptic processes in Middle Asia.
Such generalization in the classification of the synoptic processes does not fully represent all complexity of atmospheric flows and details of the synoptic processes development; to a certain extent, it is a simplified scheme. Of course, this highlights the subjectivity of classification, especially in controversial cases with deviations from one type to another that is smothered by simplified classification. On the other hand, Uzbekistan's orography variation from the east to west, with different climatic zones, makes this process extremely difficult with constant changes in weather types within even climatic zones. Nevertheless, synoptic processes over the territory of Uzbekistan should be considered as a part of much broader atmospheric circulations of Middle Asia, as these bigger scale atmospheric circulations cause the main impact on the temperature and precipitation in the region. We do not address the subjectivity of the synoptic weather type classification in this paper and used official classification for the analysis.
UzHydromet publishes actual SWT in Middle Asia in the daily and monthly catalogs. The catalog includes daily data for 00, 06, 12, and 18 UTC based on analyzing a large amount of aerosynoptic conditions data and input from the meteorological stations and satellite maps. First, analysis of the 500 hPa absolute baric topography map, the macro-synoptic situation over the Eurasian continent is estimated, and the upper-level frontal zone location is determined. At the next stage, a surface synoptic map is analyzed, and the actual location of an atmospheric object (a cyclone and its frontal zone, an anticyclone and its periphery, an atmospheric cold invasion front, etc.) relative to the territory of Central Asia is determined. The identified SWT is then entered into the daily catalog. For each subsequent synoptic period, the procedure is repeated, and the identified SWT is entered into the calendar of synoptic processes. The duration of a specific type of SWT in Central Asia in the operational practice of Uzhydromet is determined as follows. If the process was registered in the calendar in one synoptic period, its duration was 0.25 days, in two periods-0.5 days, in three periods-0.75 days, in four periods-1 day, etc. This method is widely used to describe the actual weather and assess the climate variability in Uzbekistan in relation to the conditions of regional atmospheric circulation [15,[23][24][25]. The catalog of SWT in Middle Asia can be used for synoptic, synoptic-statistical, climatological, agrometeorological, hydrological, and other scientific research and development.

Study Methods
Meteorological data were processed using standard statistical methods recommended by WMO [26]. The nonparametric Mann-Kendall statistical (MK) test is used to identify and evaluate trends in the temperature time series in Uzbekistan as well as the duration of SWT [27,28]. The MK test is widely used in scientific research to analyze trends in long-term meteorological series [9,12,14,[29][30][31][32][33]. For this reason, the description of the MK test was not included in the manuscript. In this research, the trend is assumed to be statistically significant if the level of significance is more than 95%.
For the trend analysis, meteorological datasets were divided into two parts, the present period from 1991 to 2016 and the baseline period of 1961-1990. The MK test and climate trend coefficients were applied for both the present and baseline periods, separately, and the entire observed period of 1961-2016. Results were compared to the climate normals of the long-term mean air temperatures for the baseline period ) published by UzHydromet [34]. Additionally, data were divided into cold half-year (CH), which included months between November and April, and warm-half year (WH) between May and October, to better understand the annual impact on the synoptic patterns and temperature variations. The month division to warm and cold half-year was based on the thermal criteria given by the Uzhydromet and WMO. This classification is officially approved and used for annual weather catalogs. Decadal, monthly, and annual mean value analysis of all representative stations for temperature and the SWT duration trends were prepared.

Air Temperature Trends
The results showed that the air temperature change in Uzbekistan is uneven, both spatially and temporally. Perennial monthly, semi-annual, and annual mean temperature values from 1961-2016 showed warming compared to the baseline period 1961-1990 at the representative stations in all climatic regions of Uzbekistan ( Figure 3). Overall, the most significant increase in the long-term mean temperature values (to 1.0 • C and more) was observed in the cold half of the year, mainly in January, February, and March over most of Uzbekistan's territory. Results show that winter and spring are becoming comparatively warmer on average to 1.3 • C to the baseline period across the entire Uzbekistan territory (Figure 3a). Except for October, with an average increase of 1.12 • C, temperature rise in other months of the autumn season was relatively smaller, on average 0.58 • C. The most significant increase in temperatures for the summer season was observed in June and August at 0.7 • C and 0.94 • C, respectively ( Figure 3c).  Table A2). Monthly long-term mean air temperatures for (a) the cold half-year; (b) the warm half-year; (c) annual, cold and warm averages.
The extent of temperature variations in all 69 meteorological stations over the territory of Uzbekistan on long-term mean temperatures for the entire observation period (1991-2016) compared to the baseline period (1961-1990) is shown in Figure 4.  Table A2). Monthly long-term mean air temperatures for (a) the cold half-year; (b) the warm half-year; (c) annual, cold and warm averages.
The extent of temperature variations in all 69 meteorological stations over the territory of Uzbekistan on long-term mean temperatures for the entire observation period (1991-2016) compared to the baseline period ) is shown in Figure 4.  Table A2). Monthly long-term mean air temperatures for (a) the cold half-year; (b) the warm half-year; (c) annual, cold and warm averages.
The extent of temperature variations in all 69 meteorological stations over the territory of Uzbekistan on long-term mean temperatures for the entire observation period (1991-2016) compared to the baseline period ) is shown in Figure 4.  The overall warming trend has been observed all over Uzbekistan, especially pronounced in the northwest region (Karakalpakiya, Chimbay), the climatic region closest to the former dry bottom of the Aral Sea. Spatially, the largest annual increase (more than 1.0 °C warming) in the long-term mean air temperature was observed in the northwest zone of Uzbekistan (Karakalpakiya, Ustyurt, Chimbay), and the southern, central, and western oasis zones (Kyzylkum, Termez, Samarkand, Bukhara, Ferghana) (Figure 4c). A moderate increase (0.6 °C) in the long-term mean air temperature The overall warming trend has been observed all over Uzbekistan, especially pronounced in the northwest region (Karakalpakiya, Chimbay), the climatic region closest to the former dry bottom of the Aral Sea. Spatially, the largest annual increase (more than 1.0 • C warming) in the long-term mean air temperature was observed in the northwest zone of Uzbekistan (Karakalpakiya, Ustyurt, Chimbay), and the southern, central, and western oasis zones (Kyzylkum, Termez, Samarkand, Bukhara, Ferghana) ( Figure 4c). A moderate increase (0.6 • C) in the long-term mean air temperature was observed in the desert and northeast zone (southeastern Kyzylkum, Djizakh, and the western Tien Shan) (Figure 4a). The cold half-year, as stated above, had the highest added impact to temperature increase, especially in the first three months (January, February, March) of the year (northwest 1.75 • C; desert 1.15 • C; oasis, center 1.36 • C; northeast 1.35 • C; oasis, west 1.56 • C) (Figure 3a). For the cold half-year in the long-term mean air temperature, the largest increase (more than 1.  Table A2), compared to the less populated zones in the desert and mountainous areas.
The interannual temperature change and linear trends for the baseline (1961-1990) and present periods (1991-2016) for 14 representative stations were prepared. Figure 5 shows graphs of the interannual temperature changes at the representative desert, oasis, and mountain meteorological stations. It is evident that the warming process that began in the baseline climate period continues today, and in some cases have a much steeper increasing trend. It is necessary to point out that the linear trends of all stations intersect climate normals of the baseline period  in the mid-1970s of the twentieth century and continue to increase.
Trend temperature increments for the baseline, present, and entire observed periods, including decadal mean, are calculated based on the obtained linear regression equations ( Figure 5). The nonparametric statistical MK test showed that the temperature change within the baseline period did not have a high level of significance throughout the territory of Uzbekistan. In contrast, the air temperature trends in the present 1991-2016 period had a level of significance of more than 95% at the following stations: in oasis zones 1.  Figure 6). However, Bukhara station showed a slightly decreasing rate in the decadal temperature trend, while Termez station, the most southern part of Uzbekistan, remained unchanged. interannual temperature changes at the representative desert, oasis, and mountain meteorological stations. It is evident that the warming process that began in the baseline climate period continues today, and in some cases have a much steeper increasing trend. It is necessary to point out that the linear trends of all stations intersect climate normals of the baseline period  in the mid-1970s of the twentieth century and continue to increase.     Table A3).

Changes in the Regional Atmospheric Circulation of Middle Asia
The trends in regional atmospheric circulation across the Middle Asia climatic zone were analyzed to define interannual changes in the SWT duration. Durations of the SWT identify a course  Table A3).

Changes in the Regional Atmospheric Circulation of Middle Asia
The trends in regional atmospheric circulation across the Middle Asia climatic zone were analyzed to define interannual changes in the SWT duration. Durations of the SWT identify a course of change to a long-term climate. Analysis of these processes could, therefore, represent the circulating causes of climate change in Uzbekistan. The mean long-term duration of SWT for the cold and warm six months and the whole year represents the climate formation dynamics in Middle Asia.
In the cold half-year of the baseline period , anticyclonic group processes (types 9, 9a, 9b) dominated for 64.8 days (35.8%) on a long-term mean. Cold invasions (types 5, 6, 10) were observed over 43.7 days (24.1%) on average. They were followed by southern cyclones (types 1, 2, 3) 26.9 days (14.8%) on average. Low-gradient fields of high and low pressure (types 12, 13) had a total duration of 16.5 days during this period, which amounted to 9.1% of the cold half-year. Wave activity on the cold front (type 7) and slow-moving cyclone (type 8) had a total duration of 12.7 days (7.0%) and 9.8 days (5.4%), respectively. The remaining SWT (4,11,14,15) lasted less than four days. A comparison of changes in the SWT duration between climate normals of the baseline period ) and the present period 1991-2016 is shown in Figure 7. In the warm half-year of the baseline period , cyclonic activity duration decreased to 3.9 days or 2.1%. The cold invasions became more protracted and mostly occurred in the summer; they were observed in 37.9% of the warm six months. The second most occurring types were lowgradient baric fields (28.5%), followed by anticyclones (22.4%) ( Table 1).  In the warm half-year of the baseline period (1961-1991), cyclonic activity duration decreased to 3.9 days or 2.1%. The cold invasions became more protracted and mostly occurred in the summer; they were observed in 37.9% of the warm six months. The second most occurring types were low-gradient baric fields (28.5%), followed by anticyclones (22.4%) ( Table 1). Slight changes to the cold invasion's duration were observed in the cold half-year of the present period (1991-2016) (Figure 7). However, the total anticyclonic activity duration did not change much (61.0 days and 33.6%). Variations to the total duration of southern cyclones were also small. On the other hand, the total duration of cold invasions decreased significantly, on average, up to the eight days compared to the baseline period or 35.5 days (19.6%). A decisive role was played by the shortening of northwestern and northern invasion; their total duration was halved. Another significant change was observed in the increased duration (almost three-fold) of low-gradient baric fields; the total duration has increased to 42.9 days (23.7%) from the baseline period of 16.5 days.
In the warm half-year of the present climate period (1991-2016), the main changes occurred in the structure of cold invasions and low-gradient baric fields. The total duration of cold invasions in the warm half-year decreased to 40.6 days (22.1%), and the low-gradient fields of high and low pressure increased almost twice to 91.7 days or 49.8%. Thus, the latter became the predominant type of SWT in the warm half-year of the present climatic period (1991-2016). As a result of these changes, the major shift in the total duration of the processes to the whole year occurred from cold invasions (113.4 days in the baseline period, 76.1 days in the present) toward low-gradient fields of high and low pressure (68.8 and 134.6 days, respectively).
Note that the above duration characteristics of the SWT in Middle Asia were averaged over a long time and do not represent their time dynamics. One of the methods to analyze the total duration of SWT is checking their interannual changes. The interannual changes of the total duration of SWT in Middle Asia are shown in Figure 8, where there was a deviation of the present and entire observed period from the climate normals of the baseline climate period . This technique was applied to assess changes in the total duration of SWT. Since the durations of SWT 3, 4, 8, 11, 14, and 15 types were insignificant, no trend changes were evaluated for them.
Note that the above duration characteristics of the SWT in Middle Asia were averaged over a long time and do not represent their time dynamics. One of the methods to analyze the total duration of SWT is checking their interannual changes. The interannual changes of the total duration of SWT in Middle Asia are shown in Figure 8, where there was a deviation of the present and entire observed period from the climate normals of the baseline climate period . This technique was applied to assess changes in the total duration of SWT. Since the durations of SWT 3,4,8,11,14, and 15 types were insignificant, no trend changes were evaluated for them. In order to assess changes to the total duration of the group of southern cyclones types 1, 2, and 3, trends of their cumulative total duration of the baseline (1961-1990) and present climatic periods (1991-2016) were calculated (Figure 8a). There was a small variation in trends in the baseline period, but in the present period, all cumulative total duration trends tended to decrease to negative values with a low level of significance (Table 2). Although the cumulative total duration of the group of In order to assess changes to the total duration of the group of southern cyclones types 1, 2, and 3, trends of their cumulative total duration of the baseline  and present climatic periods (1991-2016) were calculated (Figure 8a). There was a small variation in trends in the baseline period, but in the present period, all cumulative total duration trends tended to decrease to negative values with a low level of significance (Table 2). Although the cumulative total duration of the group of southern cyclones in the present period remained below the average of the baseline period, a steady increase was seen over the past few years. It was associated with changes in the duration of the South Caspian cyclone. The northwestern and northern cold invasions (types 5 and 6) and the closely related wave activity on the cold front (type 7) in the baseline ) and present climatic periods  show that the cumulative total duration of these types of SWT tended to steadily decrease from the very beginning of the baseline period (Figure 8b). The mid-1970s of the last century was a tipping point (i.e., the annual total duration of these processes became less than the mean of the baseline climate period). The decrease in the total duration of this group's processes in the warm and cold half-year was approximately the same. These changes showed a very high level of significance (95% for the cold half-year, 99% for the warm half-year and the whole year) ( Table 2). There was a decreasing tendency in the cumulative total duration of these SWTs in the present period. However, the rate of decline is slowing. In general, for the entire study period of 1961-2016, change trends had a very high level of significance (99.9%) due to the rather sharp decrease in the duration of wave activity on the cold front (type 7).
The interannual change in the group's total duration consisting of anticyclones (types 9, 9a, and 9b) showed increasing trends with a low level of significance for the baseline period (Figure 8c, Table 2). In the present period, this has changed to the opposite direction. The interannual course of the cumulative total duration of the group of anticyclones for the whole period 1961-2016 tended to be smaller with a low level of significance. The western invasion (type 10) decreased for both half-years in the baseline period   (Figure 8d) with a low level of significance (Table 2). In the present period , it decreased in the warm half-year and in general over the year. However, in the cold half-year, it increased and has been above the climate normals of the baseline period  over the past two decades, with a low level of significance.
The most significant changes were in the total duration of low-gradient fields of high (type 12) and low (type 13) pressure (Figure 8e, Table 2). The cumulative duration of this group's processes shows a steady increase in the total duration of type 12 in the warm and type 13 in the cold half-year. The tipping point again happened in the mid-1970s of the twentieth century, when it passed the baseline climate period's mean value. The level of significance of these changes was within 99% in the baseline 1961-1990 and present 1991-2016 periods. The cumulative total duration of this group of processes for the entire period of observation , both in half-year and in the whole year, showed significant growth trends (the level of significance was 99.9%).
In the present period (1991-2016), the overall duration of groups of synoptic processes such as southern cyclones and anticyclones did not change relative to the climate normals of the baseline period of 1961-1990. Groups associated with cold air invasions tended to decrease slightly but were stable within the limits of the climate normals of the baseline period 1961-1990. In general, the duration of the western invasion tended to decrease significantly over the year. However, in the cold half-year, it steadily increased. Finally, there was also a significant increase in the duration of low-gradient fields of high and low pressure in both the warm and cold half-year.

Correlation Analysis of the Air Temperature and SWT
The Pearson correlation coefficient (PCC) was used to analyze the air temperature and SWT duration correlation. The PCC is commonly used to measure the linear correlation between two sets of variables, in this case, the relationship of the distribution of the mean annual air temperature and SWT durations over Uzbekistan territory. PCC is a linear correlation value and shows a positive or negative linear correlation between −1 and +1.
Mean annual year temperature values over 1961-2016 were compared to the SWTs total durations during the year shown in the Figure 9. Advection of the cold airflows from the north and northwest in SWT (5, 6, 7) showsa negative correlation (PCC from −0.34 to −0.69) and positive correlation for the anticyclonic group (SWT 12,13) (PCC from 0.34 to 0.60). Interestingly, SWT 10, from the same group of advection of the cold airflows from the north and northwest, and anticyclonic group for 9, 9a, 9b SWT showed no correlation with almost the same duration during the year with SWT 5, 6, 7. The same was true for the group of cyclones from south SWT 1, 2, 3, which also showed no correlation. Mean annual year temperature values over 1961-2016 were compared to the SWTs total durations during the year shown in the figure 9. Advection of the cold airflows from the north and northwest in SWT (5, 6, 7) showsa negative correlation (PCC from −0.34 to −0.69) and positive correlation for the anticyclonic group (SWT 12,13) (PCC from 0.34 to 0.60). Interestingly, SWT 10, from the same group of advection of the cold airflows from the north and northwest, and anticyclonic group for 9, 9a, 9b SWT showed no correlation with almost the same duration during the year with SWT 5, 6, 7. The same was true for the group of cyclones from south SWT 1, 2, 3, which also showed no correlation.
The duration of the northwestern (type 5) and northern (type 6) invasions, and the wave activity on the cold front (type 7) were at their possible minimum values in the present period 1991-2016, relative to the baseline period of . From the circulation point of view, reduction in the North Atlantic (type 5), Arctic (type 6), and the cold front of the SWT 7 durations over the year caused a decrease in masses of cold air moving into Middle Asia. In turn, the SWT 12 and 13 duration increased with the formation of a low-gradient baric field with slightly cloudy weather (Figure 7). The air mass over the region rapidly transformed in response to the radiation factor, and the surface heat day-byday directly affected air temperature near the surface.

Discussion
This research showed that the most significant temperature changes, exceeding 1.0 °C relative to the baseline climatic period, occurred over Uzbekistan in the cold half-year, mainly in January, February, March, and October. Warming trends were especially pronounced in the northwest part The duration of the northwestern (type 5) and northern (type 6) invasions, and the wave activity on the cold front (type 7) were at their possible minimum values in the present period 1991-2016, relative to the baseline period of . From the circulation point of view, reduction in the North Atlantic (type 5), Arctic (type 6), and the cold front of the SWT 7 durations over the year caused a decrease in masses of cold air moving into Middle Asia. In turn, the SWT 12 and 13 duration increased with the formation of a low-gradient baric field with slightly cloudy weather (Figure 7). The air mass over the region rapidly transformed in response to the radiation factor, and the surface heat day-by-day directly affected air temperature near the surface.

Discussion
This research showed that the most significant temperature changes, exceeding 1.0 • C relative to the baseline climatic period, occurred over Uzbekistan in the cold half-year, mainly in January, February, March, and October. Warming trends were especially pronounced in the northwest part and almost all oasis zones over Uzbekistan, showing added anthropogenic impact compared to trends of the less-populated desert and mountainous zones. Trends in the mean annual temperature in Tashkent were 1.  [14,29,30,35]. Ayakagitma, a station in the desert zone, showed a negative decadal trend for −0.20 • C/decade, which is due to the station's location in the Ayakagitma depression zone, a relatively low area surrounded by high land. The values of the Oygaing station located in the glacial zone tended to have high trends compared to neighboring mountain zones. In typical glacial zones of China-Altai, Tien Shan, Kilian, western and eastern Kunlun, Tanggul, and the Himalayas-the mean annual air temperature over the past 58 years (1960-2017) has tended to increase at a rate of 0.38 • C/decade, 0.31 • C/decade, 0.33 • C/decade, 0.26 • C/decade, 0.37 • C/decade, 0.63 • C/decade, and 0.47 • C/decade, respectively [36].
Large-scale changes in atmospheric circulation analysis by Wen et al. [14] showed that climate change in Gansu and Northwest China is associated with increased anticyclonic circulation and rapid warming on the Eurasian continent. The regional circulation of the atmosphere over Middle Asia is an integral part of the general circulation of the Northern Hemisphere and is changing under the influence of global climate warming in the Northern Hemisphere. The major shift in the total annual duration of the SWT occurred from the interchange of cold invasions to low-gradient fields of high and low pressure, while durations of the remaining SWT groups were unchanged. Trends in the group of cold invasions (types 5, 6, and 7) and low-gradient baric fields (types 12 and 13) had a very high significance level (99.9%) over the entire period  in both warm and cold half-years. A combined analysis of the structure of SWT of Middle Asia and the temperature regime showed that the leading cause of climate warming in Uzbekistan is connected to anthropogenic and circulation factors. The air temperature increase in the baseline period was associated with a reduction in the duration of the northwestern and northern invasion and a significant increase in the low-gradient fields of high-and low-pressure durations. The low-gradient baric fields continued the warming trend in the present climatic period, and the northwestern and northern invasions are at the minimum duration value of the baseline period 1961-1990. These circumstances point to the possible contribution of the warm half-year to the warming trend in subsequent periods in Uzbekistan.
Drying up of the Aral Sea and consequently created saline desert could have probably influenced the circulation process. Uzbekistan's northern region (Ustyurt, Chimbay) has observed the largest increase in temperature both for the warm and cold half-year. Decrease in the northern (northern, northwestern) cold invasions have reduced the number of cold days, and possibly precipitation in the region. These processes accelerate the increase of days with extreme temperatures and reduce the natural vegetation cover that will negatively affect the desert ecosystems.
Although significant changes have been found in the northern part of Uzbekistan, it should be noted that this region has an insufficient density of meteorological observations to identify spatial features of climate change, especially the Ustyurt and the Northern Kyzylkum. These areas are remote and less populated, and according to Holmes et al. [37], this could alter outcomes of analysis as the same condition is seen in western China (deserts and plains). A denser observation network could provide additional data for analysis of the temperature regime changes in relation to the circulation conditions in the Central Asian region, which is undoubtedly of great importance in preventing and mitigating the adverse effects of climate change on the sustainable use and management of natural resources, especially water availability in the region.

Conclusions
Characteristics of interannual changes in air temperature in Uzbekistan, as well as the duration of SWT in Middle Asia, was analyzed using the nonparametric statistical MK test and the coefficients of climatic trends to identify changing rate. All regions showed an increase in the mean long-term monthly air temperature relative to the baseline period. Changes in air temperature in Uzbekistan had high spatio-temporal variability; a portion of the cold half-year in annual mean air temperature increase was more significant than that of the warm half-year. Climate warming is characterized by increasing air temperature on average to 1 • C in both warm and cold half-years over Uzbekistan for the period of 1991-2016. Decadal temperature changes in the present climatic period have increased relative to the baseline period and range from 0.22 • C/decade in the desert zone to 0.52 • C/decade in oasis zones.
The major shift in the duration of SWT in Middle Asia relative to the baseline period occurred from cold invasions, which have decreased to a total of 37.3 days toward low-gradient fields of high and low pressure, which increased to a total of 65.8 days. The trend in the total duration of cold invasions (types 5, 6, and 7) in the present period 1961-1990 had a very high level of significance (95% for a cold half-year and 99% for the warm half-year and whole year). However, the decreasing trend of cold invasions duration slowed down in the present . The level of significance in changes of low-gradient baric fields was within 99% in the baseline  and present period (1991-2016) and reached its maximum value of 99.9% in the whole observed period .
The geographical distribution of air temperature deviations was mosaic and mesoscale and is due to the peculiarities of the physical and geographical conditions of the climatic regions of Uzbekistan (relief, hydrography, soil type, etc.). Even within a relatively uniform landscape, climatic changes are quite complicated, both in time and meteorological values. The process of modern warming, which began in Uzbekistan in the 1960s of the twentieth century, intensified from the mid-1970s with higher values than the climate normals of the baseline period  and is associated with changes in regional SWT over Middle Asia. The anthropogenic impact on the temperature increase in Uzbekistan needs to be studied, as it has been shown that the oasis areas showed the highest increment in air temperature increase compared to less populated areas.

Group Code Type Description
Cyclones from the south and southwest 1 South Caspian cyclone Develops in the south of the Caspian Sea, a well-formed baric system with heavy precipitation in mountainous areas, winds in a desert area, and rising temperatures in winter. This cyclone usually covers the entire troposphere with its circulation, the center is in the north of Central Asia, but its formation in other regions is not excluded. Cool in summer with rainfall in mountainous areas, cold temperatures in winter.

Western invasion
The middle-latitudes air, less arctic air, from the west, coming through the Caucasus ridges and the Caspian Sea, with cool temperatures, strong winds, and rainfalls in summer and winter. 15 Diving cyclone Develop over the northern seas (Norwegian Sea, Barents, Kara) and move from north to south. Cold weather with precipitation and strong windsAnticyclonic weather.
Anticyclonic weather 9 Southwest periphery of anticyclone The periphery of the Siberian high anticyclone. Usually clear and dry weather and low winds. 9a Southeast periphery of anticyclone The stationary anticyclonic process above the Usturt plateau, the Lower Volga region, or Western Kazakhstan. Cold and clear weather in winter, cool in summer.
9b South periphery of anticyclone Part of Siberian high anticyclone covering the eastern regions of Russia's European territory, the Volga region, Western Kazakhstan.
Mostly cold and dry conditions with cool and clear conditions in summer.

Summer thermal depression
Emerge from the southeast of Middle Asia, a frontless low-pressure area with a well-pronounced cyclonic circulation in the lower troposphere. Clear, dry, hot, and hazy weather on plains.
12 Low gradient field of high pressure An area of high pressure uniting the Siberian anticyclone with the anticyclone over the European territory of Russia. Cold, dry, and clear weather with light winds.

13
Low gradient field of low pressure The low-pressure area with meridional orientation over Central Asia. Dry and warm weather conditions with low rainfall. Local heavy rainfall in summer due to convective clouds.
Mid-latitude cyclone 14 Western cyclone Strong winds, precipitations in mountainous areas.