Afforestation of Degraded Croplands as a Water-Saving Option in Irrigated Region of the Aral Sea Basin

: Climate change is likely to decrease surface water availability in Central Asia, thereby necessitating land use adaptations in irrigated regions. The introduction of trees to marginally productive croplands with shallow groundwater was suggested for irrigation water-saving and improving the land’s productivity. Considering the possible trade-offs with water availability in large-scale afforestation, our study predicted the impacts on water balance components in the lower reaches of the Amudarya River to facilitate afforestation planning using the Soil and Water Assessment Tool (SWAT). The land-use scenarios used for modeling analysis considered the afforestation of 62% and 100% of marginally productive croplands under average and low irrigation water supply identiﬁed from historical land-use maps. The results indicate a dramatic decrease in the examined water balance components in all afforestation scenarios based largely on the reduced irrigation demand of trees compared to the main crops. Speciﬁcally, replacing current crops (mostly cotton) with trees on all marginal land (approximately 663 km 2 ) in the study region with an average water availability would save 1037 mln m 3 of gross irrigation input within the study region and lower the annual drainage discharge by 504 mln m 3 . These effects have a considerable potential to support irrigation water management and enhance drainage functions in adapting to future water supply limitations.


Introduction
Climate change predictions for Central Asia point to a probable increase in atmospheric temperature coupled with a decrease in precipitation [1], which stipulate an increase in annual evapotranspiration [2]. These changes are likely to reduce crop yields by as much as 30% and cause substantial economic losses in this region, which is highly dependent on irrigated agriculture [3]. Irrigated agriculture in Central Asia annually utilizes approximately 90% of the 120 km 3 of water in the catchments of the Amudarya and Syrdarya rivers [4], while water management remains notoriously inefficient, characterized by excessive irrigation, associated water losses, and rising groundwater tables [5,6]. The other major concern that is likely to worsen owing to climate change is cropland degradation and yield declines resulting from soil salinization. The latter is rampant in the downstream countries, affecting approximately 50% of irrigated land in Uzbekistan and 95% in Turkmenistan [7][8][9].
Under the decreasing availability of freshwater for land reclamation via salt leaching, cropping operations in highly salinized areas in the Amudarya downstream are being abandoned [10,11]. Converting degraded cropland parcels to carbon-dense systems of salttolerant tree species, as examined at field and farm scales, can improve the soil fertility and The study area's average annual rainfall of 100 mm is greatly exceeded by evapotranspiration, having required the construction of an extensive canal and drain network for the irrigation of agricultural crops (Figure 1). Khorezm's location in the lower reaches in the regional topography is flat, with an elevation difference of 55 m a.s.l. toward the northwest ( Figure 2). The hydraulic gradient is low, leading to a slow lateral groundwater flow and shallow groundwater conditions, which are exacerbated by irrigation and marginal functional drainage [5,39,40]. Given that surface runoff is prevented from irrigated fields, which are leveled and surrounded by bunds, the discharge in the drainage system is generated primarily from the groundwater recharged by irrigation and leaching water losses.
Rising groundwater levels have instigated ubiquitous secondary soil salinization that is treated with pre-seasonal salt leaching [5,41,42]. Despite these efforts, approximately 20% of the irrigated cropland area is characterized by marginal productivity [11] that affects the yields of major crops and the revenues of the regional economy [13].

SWAT Modeling
SWAT, developed by the Agricultural Research Service (ARS), United States Department of Agriculture (USDA), is a physically based, daily time step, river basin or watershed-scale, and a semi-distributed hydrological model that can be integrated with GIS. SWAT is freely available and has strong network support [32]. These characteristics have allowed SWAT modelers to effectively address the main challenges faced by hydrological The study area's average annual rainfall of 100 mm is greatly exceeded by evapotranspiration, having required the construction of an extensive canal and drain network for the irrigation of agricultural crops (Figure 1). Khorezm's location in the lower reaches in the regional topography is flat, with an elevation difference of 55 m a.s.l. toward the northwest ( Figure 2). The hydraulic gradient is low, leading to a slow lateral groundwater flow and shallow groundwater conditions, which are exacerbated by irrigation and marginal functional drainage [5,39,40]. Given that surface runoff is prevented from irrigated fields, which are leveled and surrounded by bunds, the discharge in the drainage system is generated primarily from the groundwater recharged by irrigation and leaching water losses.
Rising groundwater levels have instigated ubiquitous secondary soil salinization that is treated with pre-seasonal salt leaching [5,41,42]. Despite these efforts, approximately 20% of the irrigated cropland area is characterized by marginal productivity [11] that affects the yields of major crops and the revenues of the regional economy [13].

SWAT Modeling
SWAT, developed by the Agricultural Research Service (ARS), United States Department of Agriculture (USDA), is a physically based, daily time step, river basin or watershed-scale, and a semi-distributed hydrological model that can be integrated with GIS. SWAT is freely available and has strong network support [32]. These characteristics have allowed SWAT modelers to effectively address the main challenges faced by hydrological models, such as deficient data, spatial heterogeneity, and complex processes in study catchments [37]. Nevertheless, SWAT has limited capabilities for groundwater modeling.
SWAT integrates the following main components: weather, surface runoff, return flow, percolation, evapotranspiration, transmission losses, pond and reservoir storage, crop growth and irrigation, groundwater flow, reach routing, nutrient and pesticide loading, and water transfers. The quantification of the hydrological cycle in the SWAT is based on the water balance equation: where SW t = final soil water content (mm), SW0 = initial soil water content on day i (mm), Rday = amount of precipitation on day i (mm), Qsurf is the amount of surface runoff on day i (mm), Ea is the amount of evapotranspiration on day i (mm), Wseep is the amount of water entering the vadose zone from the soil profile on day i (mm), and Qgw is the amount of return flow on day i (mm). The basic principle behind SWAT modeling is a two-step partitioning of the study watershed into sub-units. First, topographic discretization is performed by dividing the watershed into several sub-basins that are used to determine the routing structure of water and pollutants throughout a spatially heterogeneous watershed. Second, the input information for each sub-basin is grouped into the following categories: climate, hydrological response units (HRUs), ponds/wetlands, groundwater, and the main channel or reach draining the sub-basin. The HRUs are obtained by overlaying the spatial data and are lumped areas within the sub-basin that comprise a unique combination of land cover, soil type, slope, and management practices.
For this study, the SWAT model was set up between 2000 and 2010 (11 years). The simulations were run with daily time steps, which were further aggregated to annual time steps. Spatial input data, including a digital elevation model (DEM), as well as climate, land use, and soil maps available in various spatial resolutions in a grid format or as point data, were processed and prepared in the SWAT-required format (Table S1). The irrigation amount, crop growth parameters, and crop rotation data were fed into the model ( Figure 2). models, such as deficient data, spatial heterogeneity, and complex processes in study catchments [37]. Nevertheless, SWAT has limited capabilities for groundwater modeling.
SWAT integrates the following main components: weather, surface runoff, return flow, percolation, evapotranspiration, transmission losses, pond and reservoir storage, crop growth and irrigation, groundwater flow, reach routing, nutrient and pesticide loading, and water transfers. The quantification of the hydrological cycle in the SWAT is based on the water balance equation: (1) where SWt = final soil water content (mm), SW0 = initial soil water content on day i (mm), Rday = amount of precipitation on day i (mm), Qsurf is the amount of surface runoff on day i (mm), Ea is the amount of evapotranspiration on day i (mm), Wseep is the amount of water entering the vadose zone from the soil profile on day i (mm), and Qgw is the amount of return flow on day i (mm). The basic principle behind SWAT modeling is a two-step partitioning of the study watershed into sub-units. First, topographic discretization is performed by dividing the watershed into several sub-basins that are used to determine the routing structure of water and pollutants throughout a spatially heterogeneous watershed. Second, the input information for each sub-basin is grouped into the following categories: climate, hydrological response units (HRUs), ponds/wetlands, groundwater, and the main channel or reach draining the sub-basin. The HRUs are obtained by overlaying the spatial data and are lumped areas within the sub-basin that comprise a unique combination of land cover, soil type, slope, and management practices.
For this study, the SWAT model was set up between 2000 and 2010 (11 years). The simulations were run with daily time steps, which were further aggregated to annual time steps. Spatial input data, including a digital elevation model (DEM), as well as climate, land use, and soil maps available in various spatial resolutions in a grid format or as point data, were processed and prepared in the SWAT-required format (Table S1). The irrigation amount, crop growth parameters, and crop rotation data were fed into the model ( Figure 2). Computation of water balance components (water yield, AET, GWQ, percolation, soil moisture) Impact assessment of afforestation scenarios on water balance components; estimation of water saving potential Water balance components, including percolation, surface runoff, groundwater contribution to streamflow (GWQ), water yield, and actual evapotranspiration (AET) were computed using Equation (1). Using the DEM, drainage network, and river boundaries as inputs, SWAT delineated 107 sub-catchments and 6617 HRUs within the Khorezm catchment ( Figure S1). No threshold was applied to generate HRUs. The lack of ground-observed data on discharge in the drainage system precluded the calibration and validation of the SWAT model. However, the model was parameterized using the Khorezm region database and published research specific to the study region (Table S1). Plausibility checks for SWAT simulations were conducted by comparing the computed drainage ratio (the relation of simulated drainage output to the gross irrigation water supply plus precipitation) and the modeled values of water balance components with the published estimates relevant to Khorezm's irrigation and drainage system. Water balance components, including percolation, surface runoff, groundwater contribution to streamflow (GWQ), water yield, and actual evapotranspiration (AET) were computed using Equation (1). Using the DEM, drainage network, and river boundaries as inputs, SWAT delineated 107 sub-catchments and 6617 HRUs within the Khorezm catchment ( Figure S1). No threshold was applied to generate HRUs.
The lack of ground-observed data on discharge in the drainage system precluded the calibration and validation of the SWAT model. However, the model was parameterized using the Khorezm region database and published research specific to the study region (Table S1). Plausibility checks for SWAT simulations were conducted by comparing the computed drainage ratio (the relation of simulated drainage output to the gross irrigation water supply plus precipitation) and the modeled values of water balance components with the published estimates relevant to Khorezm's irrigation and drainage system.   [43].

Data Description and Analysis
Data from five locations in Khorezm were processed to capture the spatial variability in rainfall in daily time steps for 2000-2010 as previously reported [15]. SWAT considers one meteorological value for each sub-basin based on the nearest location. The climate parameters were rainfall (mm), maximum and minimum air temperature (°C), relative humidity (%), wind speed (km h −1 ), and sunshine duration (h). Sunshine hours were converted to solar radiation (MJ m −2 ) using Cropwat 8.0. Source data were collected from the Khorezm region database [44,45].
The rainfall during the study period ranged from 54 to 167 mm, with an average of  Data from five locations in Khorezm were processed to capture the spatial variability in rainfall in daily time steps for 2000-2010 as previously reported [15]. SWAT considers one meteorological value for each sub-basin based on the nearest location. The climate parameters were rainfall (mm), maximum and minimum air temperature ( • C), relative humidity (%), wind speed (km h −1 ), and sunshine duration (h). Sunshine hours were converted to solar radiation (MJ m −2 ) using Cropwat 8.0. Source data were collected from the Khorezm region database [44,45].
The rainfall during the study period ranged from 54 to 167 mm, with an average of 99 mm (    Source: [46]. The SWAT model setup of soil physical properties required 17 soil parameters for each polygon and soil layer (168 × 3 × 17 = 8568 attributes): soil name (optional), soil hydrologic group (A,B,C,D) based on the USDA classification [47], maximum rooting depth of soil profile (mm), fraction of porosity (void space) from which anions are excluded (optional), potential or maximum crack volume of the soil profile expressed as a fraction of the total soil volume (optional), texture of the soil layer (optional), depth from the soil surface to the bottom of a layer (mm), moist bulk density (g cm −3 ), available water capacity of the soil layer (Field Capacity (FC)-Wiliting Point (WP), mm H2O mm −1 soil), saturated hydraulic conductivity (Ks, mm h −1 ), organic carbon content (% soil weight), clay content weight), silt content (% soil weight), sand content (% soil weight), rock fragment content (% total weight), moist soil albedo, and the USLE equation-based soil erodibility factor (K USLE ) (0.013 (metric ton m 2 h)/(m 3 metric ton cm)). Six of the required soil parameters were not available: soil bulk density, FC-WP, Ks, moist soil albedo, soil hydrologic group, and K USLE (soil erodibility K factor for the top layer). These missing parameters were estimated using the Soil Par version 2.0 beta software. The Campbell method (retention function pedotransfer) was used to estimate Ks and FC-WP, and the Baumer method (point and function pedotransfer) was used to estimate the bulk density. Moist soil albedo values were extracted from the Moderate Resolution Imaging Spectroradiometer (MODIS) image (500 m spatial resolution) of 1 April 2010. The extracted values ranged from 0.064-0.28. K USLE was estimated based on [48]: where f csand is a factor that gives low soil erodibility factors for coarse (sandy) soils and vice versa; f ci-si is a factor that gives low soil erodibility factors for soils with high clay to silt ratios, forgc is a factor that reduces soil erodibility for soils with high organic carbon content, and f hisand is a factor that reduces soil erodibility for soils with extremely high sand content.
where ms is the sand fraction percentage (0.05-2.00 mm diameter particles), msilt is the percentage of silt (0.002-0.05 mm diameter particles), mc is the percentage of clay (<0.002 mm diameter particles), and orgc is the organic carbon concentration of the layer (%).

Land Use
The land use statistics for SWAT modeling were derived using multi-temporal land use classification using satellite data [49]. The MODIS-based Normalized Difference Vegetation Index (NDVI) time series (geometric resolution of 250 m) for 2000-2010 was used in combination with a random forest model classification approach. A quality assessment was applied to the MODIS time series, and data affected by clouds or data production problems were linearly interpolated over time. The comparison of different satellite-based vegetation index products showed that the NDVI based on 8-day surface reflectance (MOD09Q1) was the most suitable for vegetation monitoring in Khorezm [49]. The land use maps contained three non-agricultural classes, namely settlements, water, and fallow (barren) land, and five agricultural land use types: cotton, rice, wheat-fallow, wheat-rice, and wheat-other crops ( Figure 6).  31.8%, and that of rice also decreased from 5.8% to 2.0%, whereas fallow areas increased from 13.2% to 27.3%. The percentage of the wheat-only class increased slightly from 1.4% to 3.4% because the reduced water supply in 2008 limited the cultivation of the second crop after wheat.

Irrigation
The surface water applications to different land use classes were estimated based on local standards for irrigation requirements depending on crop type (predominantly cotton, followed by rice, winter wheat, and crop rotations with winter wheat) and the prevalent hydro-module zone (Table 1; [15]).  [50,51]; irrigation network efficiency of 63% for Khorezm according to [15].
All hydro-module zones in Khorezm are typified by groundwater levels of 1-2 m below the soil surface [52]. The crop irrigation requirements were further adjusted for conveyance losses in irrigation canals by considering the overall irrigation network efficiency of 63% in Khorezm [15]. In SWAT, the actual amounts of irrigation water applied (considering the irrigation supply and irrigation losses) were assigned to different crops based on the local standards for crop irrigation.
The irrigation demand of trees, i.e., Russian olive (Elaeagnus angustifolia L.), Euphrates poplar (Populus euphratica Oliv.), and Siberian elm (Ulmus pumila L.) comprised less than 2000 m 3 ha −1 in the first two growing seasons based on field-scale (<2.5 ha) afforestation experiments on marginal cropland in Khorezm [53]. After the irrigation was discontinued, the tree plantations fully relied on shallow groundwater use until the cutting age [54]

Afforestation Scenarios
The delineation of the marginal cropland, which is characterized by sufficient surface water availability to satisfy the initial irrigation demand of tree plantations [15] defined the total cropland area of low productivity but suitable for afforestation. In total, 67% (approximately 663 km 2 ) of the marginal cropland (or 16.6% of the total area in Khorezm), mostly in the southwestern part of the region, was found to be suitable for conversion to tree plantations (Figure 7). Of these, "moderately marginal" croplands were dominant, followed by "slightly marginal," "marginal," and "highly marginal. These spatial data were integrated with land use maps of 2005 (representing average conditions) and 2008 (representing drought) (Figure 6a,d) considering the following afforestation scenarios: 1.
"100% afforestation": Afforestation of the entire marginal cropland characterized by sufficient water availability for trees (67% of the total marginal cropland area, that is, approximately 663 km 2 ). The two related scenarios (LU2005-100% and LU2008-100%) are depicted in Figure 6c,f for average and drought years, respectively. 2.
These scenarios were used in the SWAT to assess the impact of afforestation on water balance components (Figures 2 and 7).

Land Use Allocation under Simulation Scenarios
Land use statistics for the years 2005 and 2008 revealed the predominant share of total and marginal croplands growing cotton, followed by those containing wheat rota- These scenarios were used in the SWAT to assess the impact of afforestation on water balance components (Figures 2 and 7).

Land Use Allocation under Simulation Scenarios
Land use statistics for the years 2005 and 2008 revealed the predominant share of total and marginal croplands growing cotton, followed by those containing wheat rotations with other crops in either year. Considering the land use changes in afforestation scenarios based on the land use map of 2005, the "100% afforestation" scenario would result in approximately 438 km 2 (10.8%) of (marginal) cotton land converted to tree plantation under average water availability conditions whereas only 263 km 2 (6.5%) of the land would be afforested in the "62% afforestation" scenario ( Figure 8). Based on the land use map of 2008, the "100% afforestation" would result in approximately 234 km 2 (5.8%) of cotton land converted to tree plantations. In the "62% afforestation" scenario, this land use change would occur only on 119 km 2 (2.9%) of cotton land.

Impact of Afforestation on Water Balance Components
The SWAT-based simulations considering the average water availability conditions (land use map of 2005) showed declining values for all the examined water balance components, namely GWQ, percolation, soil moisture, water yields (surface runoff + groundwater lateral flow), and AET as a result of both the 62% and 100% afforestation for marginal croplands (Figure 9). The GWQ was reduced by 62 mm and 116 mm in the 62% and 100% afforestation scenarios, respectively, under average water availability conditions. Similarly, percolation was reduced by 67 mm and 126 mm, whereas soil moisture showed only a small decline. Notably, AET would also decrease by 64 mm in 100% afforestation. The water yield was reduced by 66 mm and 124 mm in the 62% and 100% afforestation scenarios, respectively.
Under drought conditions, as reflected in the 2008 land use map, a similar declining trend was observed for all computed water balance components in either afforestation scenario for marginal land ( Figure 10). However, with decreased water availability and a consequently reduced cropping area and crop water use, the reductions in water balance component values were also smaller. For example, in the 100% afforestation scenario within the land use map of 2008, the reduction in AET was approximately 19 mm less, and that the water yield was approximately 44 mm less than the simulated values in the 2005 land use map.

Impact of Afforestation on Water Balance Components
The SWAT-based simulations considering the average water availability conditions (land use map of 2005) showed declining values for all the examined water balance components, namely GWQ, percolation, soil moisture, water yields (surface runoff + groundwater lateral flow), and AET as a result of both the 62% and 100% afforestation for marginal croplands (Figure 9). The GWQ was reduced by 62 mm and 116 mm in the 62% and 100% afforestation scenarios, respectively, under average water availability conditions. Similarly, percolation was reduced by 67 mm and 126 mm, whereas soil moisture showed only a small decline. Notably, AET would also decrease by 64 mm in 100% afforestation.
The water yield was reduced by 66 mm and 124 mm in the 62% and 100% afforestation scenarios, respectively.  Under drought conditions, as reflected in the 2008 land use map, a similar declining trend was observed for all computed water balance components in either afforestation scenario for marginal land ( Figure 10). However, with decreased water availability and a consequently reduced cropping area and crop water use, the reductions in water balance component values were also smaller. For example, in the 100% afforestation scenario within the land use map of 2008, the reduction in AET was approximately 19 mm less, and that the water yield was approximately 44 mm less than the simulated values in the 2005 land use map.

Plausibility of the SWAT Simulation and Estimated Water Saving
A plausibility check of the SWAT simulations was conducted for 2005, representing an average water supply and demand. Relating the simulated annual water yield to the gross irrigation input plus precipitation (2663/(4741 + 450) in mln m 3 ) resulted in a drainage ratio of 51%. This value is somewhat lower than those based on the monitoring by [6,55]. However, reference [55] reported a drainage ratio of 58% for Khorezm, but only in the vegetation (thus only in the irrigated) period of 2005; Reference [6] reported a ratio of 55% for the Water Consumers Association Shomokhulun in the vegetation period of 2007.
In the estimation of water-saving potential from replacing crops on marginal land with trees, we considered (i) the LU2005-100% scenario, that is, the change in cropping pattern in 2005 following the afforestation of all marginal land (approximately 663 km 2 ; Figure 8), (ii) the gross irrigation water input based on the crop-specific norms stated in Table 1, (iii) the water amount for pre-seasonal salt leaching reported by [41,42], that is, 500 mm for cotton, 300 mm for rice, 250 mm for wheat, and 200 mm for wheat + other, and (iv) the network efficiency of 63% [15]. The estimated reduction in gross irrigation water input amounts to 1037 mln m 3 and is the approximation for the water-saving potential. This amount equals approximately 22% of the water withdrawal for irrigation in Khorezm from the Amudarya River (approximately 4741 mln m 3 in 2005).
Furthermore, the SWAT-simulated water yield of 657 mm (Figure 9e) approximated the total drainage flow and was equal to a discharge sum of 2663 mln m 3 in 2005. As simulated by SWAT in the LU2005-100% scenario, this drainage output can be lowered to 2159 mln m 3 (533 mm; Figure 9e) because of the reduced irrigation water input following afforestation.

Limitations of SWAT-Based Analysis
Simulations based on SWAT modeling allowed for evaluating the alternative conception of land and water use at a regional scale because an analysis of the monitored data alone could not include the effect of planned strategies. The major limitation of this study is the lack of observed catchment discharge data to ensure the adequate calibration and validation of SWAT. However, we parameterized the crop and soil input data using values specific to the study area. The SWAT-computed water balance components based on land and water use in 2005, the year representing average conditions, showed a close agreement with the observed values available for Khorezm. Specifically, the drainage ratio integrating both the water yield (mainly as drainage output) and water supply (irrigation water input and precipitation) at the regional scale was matched by the SWAT simulations. This indicates the overall plausibility of the model application. The modeled PET value (1269 mm) falls within the previously reported range of 1150-1500 mm in agricultural Khorezm [40,56]. The SWAT-computed AET value (462 mm) was lower than the remote sensing-based estimate of seasonal AET (March to October) averaging 627 mm during 2003-2011 for Khorezm [57]). The other water balance components, besides soil moisture (highly variable under irrigated conditions), have not been previously reported for Khorezm, thereby precluding comparisons with the modeled data; these are therefore the first and only available estimates to date.
Another major limitation of our analysis is the lack of elaborate groundwater modeling. The groundwater contribution is roughly accounted for in the recommended irrigation norms for the hydro-module zone; groundwater also constitutes the main source of water by tree plantations introduced on marginal land. In this view, the groundwater component is important in simulating the impact of afforestation on water balance; however, it is currently not comprehensively developed in SWAT. For each HRU, the aquifer system is only represented by two layers: shallow/unconfined and deep/confined aquifers, and a linear reservoir model is used to simulate the groundwater flow. The water from the shallow aquifer could move to the deep aquifer, whereas the reverse process was not allowed, which may lead to inaccuracies in simulations for the areas where the groundwater system is more complicated and the drainage basin is predominately recharged by groundwater [58].

Impact of Afforestation on Irrigation Water Savings and Drainage
Forests and trees are important modulators of global, regional, and local hydrological cycles and patterns [22,59,60]. Negative hydrological impacts, such as the depletion of soil water and reduced annual runoff and water yields, are often observed for dryland afforestation [18,20,24], especially in large-scale planting schemes [19,61]. However, with the excessive cropland irrigation in Khorezm [5]), the reduction in all examined water balance components following the afforestation of marginal land reflects the reduced irrigation inputs and lowered drainage output. These impacts have positive implications for irrigation water saving [62] and mitigation of elevated groundwater tables [16,17]. Water-saving strategies are essential for irrigated areas in the downstream countries of Central Asia where the water consumption rates are among the highest worldwide and both irrigation water-use efficiency and water productivity remain low [9]. Converting marginal cropland, which is otherwise cultivated with low-yielding and water-consumptive crops, to tree plantations [13], offers advantages for sustaining long-term irrigated crop production. This outcome is particularly desirable as climate change increases precarity to water supplies [1,2] and sharpens competition for water in the transboundary Amudarya Basin [9].
The magnitude of reduction in water balance component values (although favorable for reducing the irrigation water losses in the region) is not likely to cause dramatic declines in soil moisture and the groundwater level, which would have jeopardized the sustainability of crop and tree water use within the study region. Even in the larger-scale afforestation scenario (LU2005-100%), only approximately 11% of cotton-cropped land can be converted to tree plantations. This land use conversion would result in financial gains considering the low yields of common crops versus the high production value of trees on marginal land, as well as the possibility of utilizing water savings from the afforestation for the cultivation of alternative, commercially attractive crops on productive cropland [13]. In Central Asia, afforestation of degraded croplands could increase low forest cover, especially in downstream areas dominated by annual crops [63], while not competing for crop production on prime farming land.
In other irrigated drylands experiencing soil salinization, such as in India and Australia, reducing annual runoff and enhancing AET via tree planting (bio-drainage) has been a deliberate strategy for controlling elevated groundwater tables and soil salinization [16,17]. To ensure effective salinity control, this bio-drainage option requires the strategic placement of trees within the landscape (possibly in combination with conventional drainage facilities), which is not consistent with the criteria used for spatial targeting of afforestation on marginal land [11,15]. In addition, according to the SWAT simulations, the AET decreased in all afforestation scenarios, in contrast to the more common expectation from trees to enhance evapotranspiration (e.g., [23]. A reduced AET can be viewed as a desirable outcome of afforestation in adaptation to climate change-induced drought events (LU2008 scenario), but it might limit the bio-drainage capacity. Overall, even if reductions in soil salinity through afforestation might not be expected in the study region, the increased productivity of marginal cropland planted with salt-tolerant tree species capable of utilizing relatively untapped groundwater resources is still a viable option for managing salinized croplands.

Conclusions
The SWAT modeling results indicated a dramatic decrease in the examined water balance components (i.e., GWQ, percolation, AET, soil moisture, and water yield) in afforestation scenarios considering both low and average water availability for irrigation. The decrease would largely be caused by the reduced irrigation demand of trees as opposed to the main crop of cotton. These findings support the hydrological feasibility of afforestation at a regional scale, that is, on 100% of marginal cropland (approximately 663 km 2 ) in the study region. Moreover, conversion of degraded cropland characterized by over-irrigation and shallow groundwater tables to a tree plantation may contribute to a portfolio of watersaving options to reduce Amudarya water withdrawals and improve the efficiency of water and land use in the lower reaches. Further modeling analysis must integrate the dynamics of groundwater components in the assessment of long-term hydrological impacts of afforestation.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/w13101433/s1, Figure S1: SWAT-delineated 107 sub-basins in the study region of Khorezm, Uzbekistan; Table S1: Input data used in SWAT model for estimation of water balance components.

Data Availability Statement:
The data presented in this study are available on request from the authors.