Assessing the Potential Impacts of Climate Change on Drought in Uzbekistan: Findings from RCP and SSP Scenarios

: Future climate change and its impact on drought is critical for Uzbekistan, located in Central Asia, the world’s largest arid zone. This study examines the evolving intensity of climate change and drought events using multi-model ensembles (MMEs) derived from the Coupled Model Intercomparison Project Phase 5 and 6 (CMIP5 and CMIP6) simulated under the Representative Concentration Pathway and Shared Socioeconomic Pathway (RCP and SSP) scenarios. The projections show different rates of increase in temperature and precipitation under the RCPs and SSPs. Projected temperature increases are expected to reach up to 2–2.5 ◦ C under SSP1-2.6, SSP2-4.5, and SSP3-7.0, by mid-century. By 2080–2099, an increase is projected of 2–3 ◦ C in monthly mean temperatures throughout the year (SSP1-2.6), and a more pronounced increase in summer up to 3–4 ◦ C (SSP2-4.5) and 4–6 ◦ C (SSP3-7.0), with a marked contrast in conditions between the mountainous and desert regions of Uzbekistan. Regional changes in precipitation over the study periods show relatively little variability, except for FD, where notable trends are found. Under SSP1-2.6 and SSP2-4.5, the increase in precipitation is relatively modest, whereas the changes in SSP3-7.0 are more substantial, with some regions experiencing variations of up to 10–20 mm per period. The Standardized Precipitation Evapotranspiration Index (SPEI), calculated based on the projected temperature and precipitation, provides an estimate of future drought trends. Our results show increasing aridity under all scenarios by mid-century, with longer-term projections indicating stabilization around different SPEI values by 2100: RCP2.6 and SSP1-1.9 stabilize around − 1.0; RCP4.5, RCP6.0, SSP2-4.5, and SSP3-7.0 stabilize around − 1.5; while RCP8.5 and SSP5-8.5 scenarios project values of − 2 or less by 2100. Notable differences in the SPEI index are found between lowland and foothill regions. In view of Uzbekistan’s heavy reliance on agriculture and irrigation, which are the sectors that are expected to be mostly affected by climate change, our study provides a scientific basis for informed policy decision-making. This includes various aspects such as planning and management water resources, as well as the broader socioeconomic development of the country.


Introduction
The profound and multifaceted challenges posed by global warming and climate change have increasingly captured the attention of the scientific community [1][2][3].The Intergovernmental Panel on Climate Change (IPCC)'s Sixth Assessment Report paints a stark picture, documenting a remarkable increase in global surface temperature.According to the report, from 1850-1900 to 2001-2020, the global surface temperature increased by 0.99 • C, accelerating to 1.09 • C when comparing 1850-1900 to 2011-2020 [4].This increase, as described by various climate models [5], is expected to amplify the hydrological cycle and cause an uneven global distribution of precipitation.Evaporation rates are expected to increase with temperature, potentially outpacing precipitation increases and escalating drought risks, especially in water-scarce regions.Since the mid-20th century, the frequency and severity of droughts have increased in many regions of Africa, East Asia, and South Asia, affecting agriculture and wild ecosystems [6].Understanding future dry-wet climate trends is critical to addressing the challenges posed by global warming.
Scenarios describing future human-induced climate change trends are a key focus of climate research.The Coupled Model Intercomparison Project (CMIP), operational since 1995, coordinates international climate model experiments to improve the understanding of historical and future climate dynamics [7].CMIP5 (initiated in 2008) [8] and CMIP6 (initiated in 2016) [9,10] have conducted extensive simulations using updated models and scenarios to project future climate change.These efforts contribute significantly to the IPCC assessments by providing insights into temperature patterns, precipitation trends, sea level variations, and climate mechanisms that are essential for the development of mitigation and adaptation strategies [7,[9][10][11][12][13].Despite the inherent uncertainties in model simulations, the multi-model ensemble mean (MME) is quite robust.Bias correction methods adjust model outputs, resulting in more accurate climate predictions and enabling policy makers to make informed decisions [14,15].
Central Asia (CA) has become a region of significant concern in the face of global warming, primarily because of its low precipitation and consequent constraints on water availability [16].Therefore, understanding the changes in natural water resources induced by climate change is of paramount importance.Research shows a remarkable trend in global dryland temperatures, with a significant increase since the late 1940s, exceeding both global and Northern Hemisphere averages [17].CA exhibits a worrisome divergence in water dynamics, with potential evapotranspiration increasing by 5.2 mm per decade, far outpacing the moderate increase in precipitation of 4.0 mm per decade [18,19].This relentless trend of warming and drying, documented between 1901 and 2002, presents a growing threat to water security in the region, particularly for western Uzbekistan, central Kazakhstan, Turkmenistan, central Mongolia, and southern Xinjiang in China, where the impacts are most pronounced [20].Projections based on CMIP5/6 scenarios indicate a continued upward trend in temperature, precipitation, and evapotranspiration throughout the 21st century [19,[21][22][23][24][25].Nevertheless, the expected increase in evapotranspiration that exceeds the concomitant increase in precipitation by 2030 will lead to an escalation of drought conditions [25,26].This underscores the critical need for vigilant monitoring of evolving climate dynamics in CA to mitigate potential impacts on water resources.
Previous research has provided estimates of the impact of climate change on drought intensity on a global [27] or regional scale [1,23,28].However, there is an urgent need for more detailed assessments focusing on sub-regions within countries.The fragile arid and semi-arid ecosystem of Uzbekistan faces an existential threat from climate change.Recently, the region has experienced a severe ecological crisis including the drastic reduction of the Aral Sea [29], increased land desertification [30,31] and an increase in extreme events, including droughts [32], sand and dust storms [33].These ecological crises have severely disrupted local industrial and agricultural production, with profound implications for the livelihoods of the population.
This study examines the evolving characteristics of future drought events in Uzbekistan and highlights their spatiotemporal patterns under different climate scenarios.Using the data from the CMIP5 and CMIP6 MMEs, we assess projected temperature and precipitation trends in the different regions of Uzbekistan under the Representative Concentration Pathway (RCP) and Shared Socioeconomic Pathway (SSP) scenarios.To comprehensively capture the evolving drought picture, we use the Standardized Precipitation Evapotranspiration Index (SPEI), which incorporates variations in temperature and precipitation [34][35][36].Our analysis reveals spatial and temporal shifts in drought intensity, frequency, and duration for two critical time periods: 2021-2060 and 2061-2100, commonly referred to as the "near future" and the "far future", respectively.The remaining sections of this paper are organized as follows: Section 2 provides a description of the datasets, models, and methods used.Section 3 presents the comprehensive assessment of drought characteristics in Uzbekistan under different climate scenarios using CMIP5 and CMIP6 data.Section 4 discusses these findings in comparison with existing research.Finally, Section 5 synthesizes the main key research findings and insightful conclusions for future drought resilience strategies in Uzbekistan.

Study Domain
Uzbekistan, located in CA (Figure 1), has extensive areas of sandy and rocky desert landscapes.The most prominent of these deserts are the Kyzyl Kum Desert, the Turan Depression, and the Ustyurt Plateau.The country is crossed by the Amu Darya and Syr Darya rivers, which are important sources of water for irrigation and agriculture development.The Tien Shan and Pamir Mountain ranges also cross the country, with the highest point being Khazret Sultan, which rises to 4643 m above sea level [16,37].Uzbekistan has a continental climate with marked seasonal temperature variations.This climate is influenced by continental air masses driven by cold northerly winds in winter and warm southerly winds in the summer months.The average annual temperature is 14.5 • C, with January temperatures around −3 • C and July mean temperatures exceeding +29 • C. The country's location within the subtropical high-pressure belt also contributes to the formation of a dry climate, characterized by low rainfall and high evaporation rates.The average annual rainfall is about 200 mm, while the average annual evaporation is about 2500 mm.
The Republic of Uzbekistan is administratively divided into 12 provinces, which are strategically grouped into five regions based on geographical homogeneity in terms of climatic conditions and water resources.This grouping is particularly aligned with the major river basins, and the regions share similarities in vulnerability levels, especially in key sectors such as water management and agriculture [38].As shown in Figure 1, the five planning zones delineated within the major river basins of Uzbekistan are the Amu Darya basin, which includes the Southern Zone (SZ), the Middle Course of the Amu Darya (MCA), and the Lower Course of the Amu Darya (LCA); and the Syr Darya basin, which includes the Fergana Valley (FV) and the Middle Course of the Syr Darya (MCS).
Uzbekistan faces significant challenges from climate change, making it one of the most vulnerable countries in Eurasia.The country is experiencing a warming trend above the global average, leading to accelerated glacier melting, reduced water resources, and increased extreme weather events like sandstorms and droughts [39,40].These impacts harm the national economy and well-being of the local population.

The Climate Research Unit Observation Data
The Climate Research Unit (CRU) observations are a series of long-term climate datasets compiled and analyzed by the Climate Research Unit at the University of East Anglia in the United Kingdom [41,42].The CRU TS v. 4.01 dataset, released in 2017, is constructed from meteorological stations contributing interpolated data values distributed over a 0.5 • longitude/latitude grid.This dataset serves as a valuable resource for providing comprehensive meteorological information for various analyses and studies, facilitating insights into climate patterns and trends on a global scale.These datasets provide information on temperature, precipitation, atmospheric circulation, and other climate variables collected from a variety of sources, including weather stations, satellites, and reanalysis models.The more recent CRU dataset has been used extensively in Asian dryland research [43][44][45][46].However, it is important to note that the pre-1960 CRU dataset may not be suitable for applications requiring high data accuracy, primarily due to the limited inclusion of ground station data.Therefore, precipitation and temperature data prior to 1960 are considered to be relatively less robust and reliable.

The ERA5 Reanalysis
ERA5 is the fifth generation ECMWF (European Centre for Medium-Range Weather Forecasts) atmospheric reanalysis of the global climate.This dataset provides a consistent and standardized compilation of a wide range of meteorological variables, including temperature, precipitation, and wind.It has a high temporal and spatial resolution of 1 h and f 0.1 • [47].Constructed through a fusion of observational data and numerical weather forecast models, ERA5 spans from 1940 to near-real time, ensuring a comprehensive and contemporary perspective.Stringent quality control procedures are an integral part of ERA5, ensuring accuracy and consistency in its many applications in meteorology, climatology, and atmospheric science [48][49][50].This dataset is used extensively in academic research and policy analysis, serving as a key tool to examine the variability and trends of meteorological variables and to understand the multiple impacts of climate change.In regions such as Central Asia, which are marked by a lack of dense measurement networks, ERA5 reanalysis data are proving to be a credible substitute for direct observations [38].

RCP Scenarios
RCPs are prominent trajectories characterizing greenhouse gas (GHG) concentrations used by IPCC for advanced climate modeling and empirical research [51].These trajectories delineate distinctive climate projections that depend on the amount of GHGs emitted in the coming years.The nomenclature of the main RCPs, namely RCP2.6,RCP4.5, RCP6, and RCP8.5, reflects the potential range of radiative forcing values projected for the year 2100, measured in watts per square meter (W/m 2 ).Each indicated pathway is delineated by a potential range of radiative forcing values by the end of the century, specifically 2.6, 4.5, 6, and 8.5 W/m 2 .These numerical representations encapsulate the projected impact on the Earth's energy balance resulting from the concentration of GHGs in the atmosphere.As a key metric, these values play a critical role in assessing potential climate outcomes associated with different emission levels and mitigation strategies.

SSP Scenarios
SSPs delineate global pathways that project different climate change mitigation and adaptation challenges.A conceptual framework underpinning the SSPs includes five alternative narratives that provide qualitative representations of the key characteristics of these pathways in terms of national population, urbanization, and gross domestic product (GDP) per capita [9,13].Broadly speaking, SSPs articulate potential developments in future societies.SSP1 [52] and SSP5 [53] represent relatively optimistic trajectories for human development, characterized by substantial investments in education and health, rapid economic growth, and the presence of robust institutions.However, SSP5 envisions an energy-intensive, fossil-based economic structure, while SSP1 illustrates a gradual shift toward sustainable practices.In contrast, SSP3 [54] and SSP4 [55] project more pessimistic development trends, characterized by limited investment in education and health, burgeoning populations, and escalating inequalities.SSP3 emphasizes regional security, while SSP4 depicts scenarios dominated by extensive inequalities within and between countries, culminating in societies highly vulnerable to the impacts of climate change.SSP2 [56] outlines a central pathway in which prevailing trends continue along historical trajectories without significant deviation.This pathway represents a moderate and balanced narrative within the spectrum of socioeconomic development scenarios.
Using the SSP framework, the IPCC Sixth Assessment Report assessed expected temperature outcomes based on a set of five scenarios [57].These scenarios are named by combining the respective SSP (SSP1-SSP5) with the expected level of radiative forcing projected for the year 2100 (ranging from 1.9 to 8.5 W/m 2 ).Consequently, the nomenclature of the main SSPs results in the following scenario names: (1) SSP1-1.9:Significantly low GHG emissions, characterized by the achievement of carbon dioxide (CO 2 ) emissions reduced to net zero around the year 2050.(2) SSP1-2.6:Indicating low GHG emissions, with a focus on achieving net zero CO 2 emissions around the year 2075.(3) SSP2-4.5:shows intermediate GHG emissions, with a CO 2 emission projected to continue at current levels through 2050, followed by a decline, but not to net zero by 2100.(4) SSP3-7.0:shows high GHG emissions, accompanied by a scenario in which CO 2 emissions are projected to double by 2100.(5): SSP5-8.5:Reflects very high GHG emissions, with a scenario projecting a tripling of CO 2 emissions by 2075.

CMIP5/CMIP6 Model Dataset
The study uses output data from two different generations of climate models, CMIP5 and CMIP6.The output data from the CMIP5 models are available at https://aims2.llnl.gov/search/cmip5/ (accessed on 10 July 2024), with a selected set including scenario experiments for RCP2.6/4.5/8.5 across 30 climate models (see Table A1, Appendix A).Meanwhile, CMIP6 model output data are available at https://aims2.llnl.gov/search/cmip6/ (accessed on 10 July 2024).Of the 58 available CMIP6 models, 27 were carefully selected based on specific criteria related to the availability and quality of the required parameters, as shown in Table A2, Appendix A. In particular, the selection process considers only the use of concentration-based experiments where concentrations, as opposed to GHG emissions, are prescribed.In addition, a single ensemble member is selected from each model in an effort to balance the weights across models, regardless of the existence of multiple versions with different initial conditions and different representations of the physical process.This approach ensures a balanced and consistent integration of the selected models for the calculation of multi-model ensembles.The number of models used for each RCP/SSP scenario may vary.

The SPEI Drought Index
SPEI has been proposed as an improved drought index, with a specific focus on investigating the effects of global warming on drought severity [34].The SPEI is distinguished by its capability to consider the impact of reference evapotranspiration on the severity of drought and its ability to discriminate between different droughts and how they affect different systems.[35].Here, we calculate SPEI over a 12-month period (SPEI12), which provides a perspective on the accumulation and depletion of water within this specific period.The SPEI12 is often used to examine interannual variability in dry and wet characteristics and is critical for identifying patterns in climatic conditions.The classification system that describes the degree of drought and wetness is shown in Table 1 [34].This categorization assists in interpreting the severity of drought events and assessing moisture conditions based on the calculated SPEI values, providing a comprehensive framework for understanding the climatic dynamics under study.
Table 1.The SPEI drought index categories [34].In our paper, the calculation of SPEI is implemented based on the R package SPEI (http://cran.r-project.org/web/packages/SPEI,accessed on 6 April 2021).Among several options available in the package, the logarithmic probability density function distribution recommended by [35], the unbiased PWM (Probability Weighted Moment), and the Hargreaves method [58] were used to calculate the monthly potential evapotranspiration, which was calculated in a modified form presented in [59].

Definition of Specific Periods
The baseline periods refer to the historical simulations for 1986-2005 and 1995-2014 for the CMIP5 and CMIP6 cases, respectively.The projected periods under different emission scenarios are 2020-2039, 2040-2059, 2060-2079, and 2080-2099.

Climate Model Bias Correction
In the context of climate modeling, a bias correction in CMIP5/6 refers to a statistical adjustment applied to correct systematic errors or biases in the output of climate models [8,11,60,61].These biases can arise from various sources, such as limitations in model physics, numerical approximations, or inaccuracies in the representation of processes affecting the climate system.The goal of bias correction is to improve the reliability and accuracy of climate model simulations by making them more consistent with observed climate data.This process aims to reduce discrepancies and ensure greater consistency between modeled results and actual climate observations.This correction typically involves comparing the model's historical simulations with corresponding observations for the same time and location.Any discrepancies or biases between the model's output and the observed data are quantified and then adjusted through statistical methods.Bias correction can be applied to various climate variables, including temperature and precipitation [8,11,60,61].By implementing bias correction using ERA5 and CRU, we aim to produce model simulations that provide a more realistic representation of historical climate conditions (Figures 2 and 3).We found that the original model values overestimate (underestimate) precipitation in the FV (MCS, SZ, and MCA) regions, and overestimate (underestimate) temperature in the SZ and MCA (FV and MCS) regions.This phenomenon highlights the influence of orography, with the model showing minimal error in the flattest region (i.e., LCA).The divergent temperature and model errors observed in these specific regions emphasize the need to consider them individually.It is important to note that while averages for the entire territory of Uzbekistan may be significantly erroneous due to aggregation effects, a region-specific approach allows for a more accurate representation of climatic conditions.

Time Series Analysis of Historical and Climate Projection Datasets
Since the 1970s, the pronounced warming trend in the Northern Hemisphere, including the CA regions, has been striking [62].In particular, the increasing warming, with an average rate of 0.39 • C per decade from 1979 to 2011 in CA, exceeds the average rate of temperature change observed in global land areas [63] and several other regions [64].In addition, the rate of warming in CA during the first two decades of the 21st century exceeds that of previous decades (see Figure 2).This increased rate in CA during the early 21st century is consistent with the averaged rates observed for Russia and China and exceeds that of Europe [29,[63][64][65].
In response to this rapid warming, precipitation across the vast arid expanse of CA is undergoing a significant increase, characterized by greater magnitudes of variability [66].A sustained increase in precipitation has been observed since the late 1970s, possibly related to persistent shifts in the mid-latitude westerly circulation (see Figure 3).Of note is the extensive interannual variability in precipitation patterns within CA, marked by distinct cycles of 5-6 years and 2-3 years.This recurrent pattern plays an important role as the dominant interannual climate signal, influencing both atmospheric circulation and precipitation over the mid-latitude Asian region.In particular, the 2-3-year cycle emerges as a key driver influencing atmospheric dynamics, evident in various climate variables in the troposphere and within the mid-latitude westerly circulation and Asian summer monsoon boundary regions.Further extensive study of this dominant cycle is essential to fully understand its mechanisms and implications for the broader climate system in mid-latitude Asia.

Spatial Variations of Projected Temperature and Precipitation
The RCPs were completed in time to be used in the IPCC Fifth Assessment Report [51] and replaced by the SSP scenarios, which are widely applied in the scientific community and are expected to be used in the future.Therefore, we will continue our analysis of temperature and precipitation based on SSPs only.It is unclear which SSP is more likely or unlikely, as it depends on policy choices, technological advances, and other factors that influence future developments.However, SSP1-1.9 and SSP5-8.5 are the most divergent scenarios in terms of the future trajectory of human activities and GHG emissions.SSP1 is considered a scenario of low population growth and rapid economic development, with strong policies to reduce GHG emissions, and is considered the most optimistic scenario for achieving the goals of the Paris Agreement.On the other hand, SSP5 is a scenario of high population growth and low economic development, with weak policy actions to reduce GHG emissions, and it is seen as the least optimistic scenario for achieving the goals of the Paris Agreement.Therefore, RCPSSP1-1.9and SSP5-8.5 are not considered here as they are less likely.
Compared to the historical reference period (1995-2014), both precipitation and temperature in the SSP scenarios show a predominant increasing trend in the short, medium, and long term, which is particularly noticeable over time and in scenarios with higher emissions (see Figures 4 and 5).In particular, the increase in precipitation is found mainly in mountainous areas.With respect to temperature, there is a consistent pattern of more pronounced increases across all scenarios in the LCA and SZ regions.Based on the assessments of all scenarios, the temperature in LCA is projected to increase by an average of 1 • C in the period 2020-2039.In the subsequent period of 2040-2059, the SSP2.6 and SSP4.5 scenarios project an average temperature increase of 1.5 • C, while SSP7.0 projects a temperature increases of nearly 2 • C. Despite the SSP scenario, the projections indicate a worsening situation by 2100.Of particular concern is the escalating temperatures outside of Uzbekistan in mountainous regions that are critical to the formation of river basins, which may negatively affect the state of perennial snow and glaciers in the region.

Seasonal Anomaly According to Climate Projection Simulations
Projections for Uzbekistan under different SSP scenarios show a consistent trend of increasing temperatures over all future periods (2020-2099), with more severe impacts under higher emission scenarios (Figure 6).By mid-century (2040-2059), temperature increases of up to 2-2.5 • C during summer months are expected under SSP1-2.6,SSP2-4.5, and SSP3-7.0, with stronger effects in winter and spring for SSP1-2.6 and notable temperature increases in summer for SSP2-4.5 and SSP3-7.0.In addition, there are regional differences, with mountainous areas showing a more significant temperature increase in winter and spring, while LCA, MCA, and SZ show a more significant temperature increase in summer.By 2080-2099, the scenarios show different temperature changes.SSP1-2.6 shows a steady increase of 2-3 • C in monthly mean temperatures throughout the year.SSP2-4.5 shows more pronounced increases, especially in summer (June-August), reaching up to 3-4 • C. SSP3-7.0, the high emissions scenario, predicts the most drastic warming of 4-6 • C in certain months, especially in summer.
Across all future periods (2020-2099), there is a general trend of increasing precipitation, particularly notable during the winter months (Figure 7).Under SSP1-2.6 and SSP2-4.5, the increase in precipitation is relatively modest, while SSP3-7.0indicates more substantial changes, with some regions experiencing variations up to 10-20 mm.The figures reveal that mountainous regions are likely to experience more pronounced increases in precipitation compared to the lowland areas.Seasonal variations are evident, with the spring and autumn months showing smaller increases or even decreases in some scenarios.The rates of temperature increase show different patterns over the periods (Figure 8).SSP1-2.6 reveals a notable acceleration in the first half of the century, particularly in summer, followed by a more gradual increase from 2040.As the century progresses, the winter and spring seasons experience more pronounced warming.In contrast, SSP2-4.5 shows a steady and uniform temperature increase over the specified periods, indicating a sustained upward trend throughout the century.SSP3-7.0 shows a more pronounced temperature increase, particularly during the summer months.Regional changes in precipitation over the study periods identify relatively low variability, except for FD, where notable trends are projected (Figure 9).In SSP1-2.6,there is an increase in precipitation that occurs mainly in the winter and spring months, with a maximum in March.In SSP2-4.5, there is a significant increase in precipitation at the beginning of the year, from January to April.Conversely, SSP3-7.0 shows a distinct pattern with an increase in winter precipitation from December to February, followed by a decrease in March and April.

Projected Variations of the SPEI Index
The SPEI index serves as a valuable component in the assessment of impending drought trends [35].The SPEI provides a thorough assessment, considering variations in temperature, precipitation, and potential evapotranspiration to identify essential characteristics of drought events, including their onset, end, duration, and intensity, among other parameters.The identification of drought events depends on three primary criteria: (1) a consistent SPEI value persistently below 0; (2) the occurrence of the lowest SPEI value within the designated time frame falling below −1; and (3) the cumulative duration reaching or exceeding three months [35,36].
In the short term, all scenarios indicate an increase in drought by 2060.In the long term (2060-2100), however, different trends emerge.Specifically, the RCP2.6 and SSP1-1.9scenarios stabilize around −1.0; the RCP4.5, RCP6.0,SSP2-4.5, and SSP3-7.0scenarios stabilize around −1.5; while the RCP8.5 and SSP5-8.5 scenarios reach values of −2 and even lower by 2100.Although the annual mean temperature and precipitation values are largely similar under the RCP and SSP scenarios, notable differences in the SPEI index occur between the western and eastern regions, as shown in Figure 10.This divergence can be attributed to the finer modeling of the temperature contrast due to the higher resolution of the model grids and a more accurate representation of the near-surface layer microphysics.Consequently, the foothills show a stabilization of about −0.5 for all SSP scenarios except the most extreme one (SSP5-8.5).In addition, the FD region could potentially return to conditions like those of 2015-2020, mainly due to a significant increase in precipitation.

Discussion
Understanding the variability and intensity of droughts over long periods of time is critical for water resource management, agricultural planning, and mitigating the effects of climate change.SPEI is a widely used metric that combines precipitation and potential evapotranspiration to assess drought conditions.This index is particularly valuable because it captures the effects of temperature on drought severity, providing a comprehensive measure of moisture availability [36, [67][68][69].The SPEI has been successfully applied to a wide range of climatic conditions and time periods, demonstrating its robustness and reliability [69][70][71].This extensive application in different scenarios ensures that the method is well validated and suitable for the specific requirements of this study, providing confidence in its results and conclusions.
The need for bias correction in GCM outputs is critical due to their inherent systematic biases and coarse spatial resolutions [8,11,60,61], which make them unsuitable for various applications without adjustments [72][73][74][75][76]. Methodological choices in bias correction can significantly influence the representation of temperature and precipitation, especially for extreme events, leading to variability in impact assessments related to hydrology and ecosystem dynamics.Mixed conclusions in studies assessing the importance of bias correction highlight the complexity of this process.For example, some studies suggest that the choice of bias correction methods does not significantly affect projection ranges [76].In contrast, others warn that reliance on a single set of bias-corrected GCMs may lead to overconfidence in projections [74].These discrepancies are often due to differences in the geographic focus of the studies, the meteorological or hydrological variables examined, and the specific GCMs and scenarios used.In this study, a region-specific approach to bias correction improves the accuracy and fidelity of climate state representations by adapting correction methods to the unique geographic and environmental factors that influence local climate data.However, some uncertainty may remain despite recent advances in GCM bias correction methods, which support the credibility of climate projections by effectively mitigating systematic biases and enhancing the reliability of future climate predictions.Therefore, continued improvements in bias correction techniques are crucial for refining climate models and ensuring thorough and robust assessments of climate impacts.
The collective results of several large-scale studies consistently confirm the presence of a robust warming trend, which is particularly evident over CA, and has significantly intensified in recent years [17][18][19][21][22][23][24][25].Our seasonal analysis highlights winter, spring, and summer as experiencing the most pronounced shifts in temperature and precipitation.Notably, previous studies are consistent with our results, consistently identifying winter as the season experiencing the most rapid warming in the region [21][22][23][24][25].The apparent emphasis on colder seasons, which experience enhanced warming trends in semi-arid areas, is noteworthy.Discrepancies in trend magnitudes may be due to differences in regional coverage, time periods, and data sources used.These methodological differences potentially affect findings not only within CA studies, but also in broader global studies.
While complex irrigation systems and water diversion between river basins may potentially mitigate future drought effects, this research primarily examines long-term climate change patterns in Uzbekistan and analyzes regional drought characteristics.Evaluating drought in this context is complex and requires a comprehensive regional approach rather than a narrow national focus.However, the primary objective of this article is to address broader climate change trends that are unique to Uzbekistan, rather than to provide a regional drought assessment.
Considering these discernible climate trends, it is imperative to recognize the expected significant impacts on Uzbekistan's dryland water cycle.Foreseeable consequences include increased drought problems and an increase in the frequency of dust storms.These changes are expected to exacerbate the degradation of the region's forest and grassland ecosystems, putting further pressure on local agricultural production.Given the interconnectedness of these challenges, it is critical to foster international cooperation and collaboration among the countries of the region.Addressing these myriad challenges will require the implementation of comprehensive climate change adaptation and disaster risk reduction strategies.Such proactive measures are essential to ensure the resilience and sustainability of Uzbekistan's ecosystems and agriculture in the face of a dynamically changing climate.By adopting adaptation strategies and cooperating at the regional and international levels, Uzbekistan can work towards mitigating the adverse impacts of climate change and promoting a more sustainable future for its environment and agricultural sectors.

Conclusions
This study examined the impacts of climate change on Uzbekistan MMEs derived from CMIP5 and CMIP6 simulated under the RCP and SSP scenarios.The long-term climate datasets from the CRU and the ERA5 reanalysis tracking historical climate patterns and trends were used in the bias correction method applied to the MME outputs to improve their reliability and accuracy.The corrected models showed that the original values often over-or underestimated key climate variables, highlighting the importance of regionspecific adjustments.
The scenarios ranged from optimistic trajectories with strong mitigation efforts (RCP-2.6,SSP1-1.9) to pessimistic trajectories with high emissions (RCP-8.5,SSP5-8.5).These projections consistently indicated increasing temperatures and precipitation across Uzbekistan, with more severe impacts under higher emission scenarios.The study highlighted regional differences, with the LCA and SZ regions experiencing more pronounced temperature increases.
By mid-century (2040-2059), temperature increases of up to 2-2.5 • C during the summer months are expected under SSP1-2.6,SSP2-4.5, and SSP3-7.0, with stronger effects during winter and spring for SSP1-2.6 and notable temperature increases in summer for SSP2-4.5 and SSP3-7.0.In addition, there are regional differences, with mountainous areas showing stronger temperature increases in winter and spring, while LCA, MCA and SZ show increases in summer.By 2080-2099, the scenarios show different temperature changes.SSP1-2.6 shows a steady increase of 2-3 • C in monthly mean temperatures throughout the year.SSP2-4.5 shows more pronounced increases, especially in summer (June-August), reaching up to 3-4 • C. SSP3-7.0, the high emissions scenario, predicts the most drastic warming of 4-6 • C in certain months, especially in summer.
Under SSP1-2.6 and SSP2-4.5, the increase in precipitation is relatively modest, while SSP3-7.0indicates more substantial changes, with some regions experiencing variations of up to 10-20 mm.The figures show that mountainous regions are likely to experience more pronounced increases in precipitation than lowland areas.Seasonal variations are evident, with spring and fall months showing smaller increases or even decreases in some scenarios.Regional changes in precipitation over the study periods show relatively little variability, except for FD, where notable trends are found.In SSP1-2.6,there is an increase in precipitation that occurs mainly in the winter and spring months, with a maximum in March.In SSP2-4.5, there is a significant increase in precipitation at the beginning of the year, from January to April.Conversely, SSP3-7.0 shows a distinct pattern with an increase in winter precipitation from December to February, followed by a decrease in March and April.
Finally, the SPEI index proves indispensable for predicting future drought trends, providing a nuanced assessment of drought characteristics based on temperature, precipitation, and moisture factors.Our results show increasing aridity under all scenarios through 2060, with longer-term projections indicating stabilization around various SPEI values through 2100: RCP2.6 and SSP1-1.9stabilize around −1.0; RCP4.5, RCP6.0,SSP2-4.5, and SSP3-7.0stabilize around −1.5; while RCP8.5 and SSP5-8.5 scenarios project values of −2 or less by 2100.While annual mean temperature and precipitation values show little variation across RCP and SSP scenarios, notable differences in the SPEI index are between western and eastern regions.For example, the foothill regions show stabilization around −0.5 in all SSP scenarios except SSP5-8.5.In addition, the FD region could potentially return to conditions like 2015-2020, primarily due to a significant increase in precipitation. in this research.

Conflicts of Interest:
The authors declare no conflicts of interest.The following designations are used here:

Appendix A
• AORI (Atmosphere and Ocean Research Institute, Japan);

Figure 1 .
Figure 1.The study areas.The 5 regions are defined as the Southern Zone (SZ), Middle Course of Amu Darya (MCA), Lower Course of Amu Darya (LCA), Fergana Valley (FV), and Middle Course of Syr Darya (MCS).

Figure 2 .
Figure 2. Annual mean temperature ( • C) time series for Uzbekistan derived from the CRU observations for 1901-2021, the ERA5 reanalysis for 1970-2020, and the CMIP5 and CMIP6 historical (for 1986-2005 and 1995-2014, respectively) and climate projection (for 2006-2100 and 2015-2100, respectively) simulations.The left and right panels show the original and bias-corrected model results, respectively.

Figure 3 .
Figure 3.Time series of annual precipitation (mm) for Uzbekistan derived from the CRU observations for 1901-2021, the ERA5 reanalysis for 1970-2020, and the CMIP5 and CMIP6 historical (for 1986-2005 and 1995-2014, respectively) and climate projection (for 2006-2100 and 2015-2100, respectively) simulations.The left and right panels show the original and bias-corrected model results, respectively.

Table A1 .
A list of CMIP5 models selected for analysis and compilation of the ensemble.

Table A2 .
A list of CMIP6 models selected for analysis and compilation of the ensemble.