Spatio-Temporal Reorganization of Cropland Development in Central Asia during the Post-Soviet Era: A Sustainable Implication in Kazakhstan

The abandonment and cultivation of croplands in the Eurasian Steppe has become the focus of global and regional food security and agricultural policy-making. A large area of cropland in some post-Soviet countries has proven to be abandoned with the disintegration of the Soviet Union; however, it is unclear as to whether Kazakhstan also experienced a similar change as one of the main food providers for the former Soviet Union. In this study, we used the annual land cover dataset (1992–2015) from the European Space Agency, Climate Change Initiative (ESA-CCI) to detect spatio-temporal characteristics of rainfed and irrigated cropland changes in Kazakhstan. The Mann–Kendall test and regime shift analysis showed that rainfed and irrigated cropland at national level had a significant increasing trend with a significant rising up to 1999 and stagnation during 2000–2015, which was further confirmed with analysis at 14 regions. The greatest contributor to rainfed and irrigated croplands was steppes, followed by shrinkage of water bodies since 2005 to a great extent, rapid urbanization process resulted in losses of a part of irrigated oases. The trend surface analysis indicated that reorganized stable pattern characterized by rainfed cropland in north and irrigated cropland in south was driven by the strategy of the gradual agricultural development of oases. The nonexistence of cycle between the abandonment and recultivation proved that newly-gained cropland from steppe may be less degraded and more productive for sustainable land use in Kazakhstan. In conclusion, this study can provide strong evidence for sustainable land use and a basis for food security policy-making in Kazakhstan, and even all of the Central Asian countries in the future.


Introduction
The Eurasian Steppe has historically served as the home for pastoral nomads [1][2][3]. Since the last Ice Age, this large inland area had been disturbed by the encroachment of sedentary civilization, which was accompanied by intensified grazing and cropland cultivation [4,5]; however, this change was gradual and fluctuated. Drastic changes started at the beginning of the 20th century when increasingly larger areas of the Eurasian Steppe became parts of communist states, especially the former Soviet Union [6][7][8], and peaked between the 1950s and 1960s with the Virgin Land Program in the former Soviet Union [9] and the Great Leap Campaign in China [10]. These newly rising country, with temperatures as high as 43°C in the southwestern deserts, while the winter temperatures often fall to −18°C and can be as cold as −46°C. The average annual number of days of sunshine for most of Kazakhstan is 200 days [35,36]. There are few scattered xerophytic (drought-resistant) plants that exist in the arid desert area.
Steppe is the main ecosystem type in semiarid plain regions, and woodlands occur mainly in the cooler and wetter mountainous areas. More than 500 native plant species and dozens of wild animals grow only in Kazakhstan. An estimated 150 species of fish inhabit Kazakhstan's several thousand lakes. At present, Kazakhstan also has nine nature reserves [18].
Approximately 12% of Kazakhstan has soils that are suitable for agriculture. In much of the western, southern, and central portions of the country, soils are alkaline, sandy, and often saline-very poorly suited for agriculture. In the eastern and southeastern upland areas, the soil dominated by loams is of good quality, which presents an opportunity for high-yield farming. Spring wheat and barley are the principal crops grown using a dryland cultivation system of rotation with fallow every 3-4 years [22,37].
Kazakhstan is one of the least densely populated countries with a 16.7 million population, land is the largest natural resource but receives little investment from the Republic and local government; for example, among an overall budget for seven years (2004)(2005)(2006)(2007)(2008)(2009)(2010) amounted to KZT 976 million (USD 7 million) for rural territory, only KZT 2 million (USD 14 thousand) was allocated to land issues [22].Obviously, agriculture is a huge and neglected industry in Kazakhstan.

Dataset Preprocessing
In this study, 24 consistent global annual land cover maps , with a 300 m spatial resolution, were chosen to analyze cropland dynamics. The coordinate reference system (CRS) used for the global land cover databases is a geographic coordinate system (GCS) based on the World Geodetic System 84 (WGS84) reference ellipsoid. This dataset was provided by the CCI-LC (Climate Change Initiative, Land Cover) project, ESA (European Space Agency) (http://maps.elie.ucl.ac.be/ CCI/viewer/download.php), and has been applied in cropland, forest, and plant functional type research at the regional and global scales [38][39][40]. We masked this global land cover dataset with There are few scattered xerophytic (drought-resistant) plants that exist in the arid desert area.
Steppe is the main ecosystem type in semiarid plain regions, and woodlands occur mainly in the cooler and wetter mountainous areas. More than 500 native plant species and dozens of wild animals grow only in Kazakhstan. An estimated 150 species of fish inhabit Kazakhstan's several thousand lakes. At present, Kazakhstan also has nine nature reserves [18].
Approximately 12% of Kazakhstan has soils that are suitable for agriculture. In much of the western, southern, and central portions of the country, soils are alkaline, sandy, and often saline-very poorly suited for agriculture. In the eastern and southeastern upland areas, the soil dominated by loams is of good quality, which presents an opportunity for high-yield farming. Spring wheat and barley are the principal crops grown using a dryland cultivation system of rotation with fallow every 3-4 years [22,37].
Kazakhstan is one of the least densely populated countries with a 16.7 million population, land is the largest natural resource but receives little investment from the Republic and local government; for example, among an overall budget for seven years (2004-2010) amounted to KZT 976 million (USD 7 million) for rural territory, only KZT 2 million (USD 14 thousand) was allocated to land issues [22].Obviously, agriculture is a huge and neglected industry in Kazakhstan.

Dataset Preprocessing
In this study, 24 consistent global annual land cover maps , with a 300 m spatial resolution, were chosen to analyze cropland dynamics. The coordinate reference system (CRS) used for the global land cover databases is a geographic coordinate system (GCS) based on the World Geodetic System 84 (WGS84) reference ellipsoid. This dataset was provided by the CCI-LC (Climate Change Initiative, Land Cover) project, ESA (European Space Agency) (http://maps.elie.ucl.ac.be/ CCI/viewer/download.php), and has been applied in cropland, forest, and plant functional type research at the regional and global scales [38][39][40]. We masked this global land cover dataset with the administrative boundary of Kazakhstan and its 14 regions, which was acquired from the Global Administrative Areas (GADM) dataset (V3.6, released on 6 May 2018, https://gadm.org/index.html) with longitude/latitude the coordinate reference system and the WGS84 as the datum. Then, we aggregated the 37 subdatasets of the ESA land cover classes for Kazakhstan into 8 classes according to the UNCCD (United Nations Convention to Combat Desertification and Drought) land cover category with the following two important differences: (1) Since the majority of the Kazakhstan territory is desertification-prone with high-frequent drought [41,42], water was the most significant limiting factor for dryland agriculture [43,44]; hence, croplands were subdivided into rainfed cropland and irrigated cropland in this study. (2) The word steppe was used instead of grasslands in this study. As one Russian word for "grassy plain", steppe is mostly dominated by short grasses without trees [2].
The new classification for the ESA-CCI land cover dataset is listed in Table 1.

Analysis Methods
First, we divided the whole of Kazakhstan into 14 administrative regions at spatial scale and the whole study period (24 years) into five subperiods (SPs) with a five year interval except the first subperiod with four years (1992)(1993)(1994)(1995)  , which is a suitable temporal scale for time series related analysis [45,46]. Then we conducted a time series analysis to detect significant changes in cropland (rainfed and irrigated) at the national and regional levels during the whole study period. In addition, FAO (Food and Agriculture Organization of the United Nations) arable land statistics for Kazakhstan (available from http://www.fao.org/faostat/en/#country/108) were also analyzed for comparison. Here, the Mann-Kendall test method was employed for time series analysis, which is an often-used nonparametric test method to analyze trends in meteorological data [47,48] because it is less sensitive to outliers [49][50][51], and has been broadly applied in hydrology and water resource research [52,53]. The test is based on the Kendall's S statistic as follows: And var(S) = n(n − 1)(2n where S is Kendall's statistic, which is expected to have N(0,var(S)) distribution with the null hypothesis H0 that there is no trend displayed by the time series. x k and x i are sequential data values representing the annual values in years k and i, respectively; n is the length of the dataset; sgn(θ) = −1 if θ < 0, sgn(θ) = 0 if θ = 0, sgn(θ) = 1 if θ > 0; e i is the number of ties of extent i. The Mann-Kendall test statistic Z c is computed by: where, Z c is a standard normal variable, The null hypotheses that Z c is not statistically significant or has no significant trend will be accepted if where ±Z 1−α/2 ±Z1 − α /2 are the standard normal deviates and α α is the significance level for the test [54]. A positive Z c indicates an upward trend while a negative Z c indicates a downward trend [55]. The calculations were performed using the 'rkt' package [56] in the R software [57]. Sequential t tests for an analysis of regime shifts (STARS) [58] were employed simultaneously to test mean and variance values for statistically significant shifts between alternate states at national and regional levels during the whole study period. Attribute to its easy-to-use and statistically powerful analysis for shifts near the ends of the time series [59], this method has been applied mainly in oceanography [60][61][62] and expanded to terrestrial ecosystem applications recently [63][64][65]. In this method, when a new observation is added into a time series data, it is tested if there is a statistically significant deviation from the mean value of the current regime, if it is true, that time for this observation is marked as a potential change point c, and subsequent observations are used to accept or reject this hypothesis, which is tested by the regime shift index (RSI) as follows: where m = 0, 1, . . . , l − 1 (i.e., number of years since the start of a new regime), l is the cut-off length of the regimes tested, and σ l is the average standard deviation for allone-year intervals in the time-series. RSI is a cumulative sum of normalized deviations x * i from thehypothetical mean level for the new regime (x new ), for which the D, difference from the mean level for the current regime (x cur ), is statistically significant according to a Student's t-test: where t is the value of the t-distribution with 2l − 2 degrees of freedom at the given probability level p. If, at any time from the start of the new regime, RSI becomes negative, the test fails and a zero value is assigned. If RSI remains positive throughout l − 1, then c is declared to be the time of a regime shift at the level ≤ p. The search for the next regime shift starts with c + 1 to ensure that its timing is detected correctly even if the actual duration of the new regime is less than 1 year. In this analysis, STARS software was used to detect a regime shift of rainfed and irrigated cropland development periods from 1992 to 2015, and the software was parameterized as a probability equal to 0.05 to coincide with the Mann-Kendall test, cutoff length equal to five to coincide with the time steps of subperiods, and a Huber parameter equal to one for the default value [58,59]. Second, the Land Change Modeler (LCM) provided by the IDRISI TerrSet software, which had been applied to study LULCC in many studies [66,67], was used to evaluate gains and losses in land cover classes, land cover persistence, and specific transitions between selected categories based on cross-tabulation principles [68]. Using the aggregated land cover maps from 1992 to 2015 as input parameters, we identified the locations and magnitude of the rainfed and irrigated cropland changes and persistence. Then, the other land cover contributions to the net areal change in rainfed and irrigated cropland were analyzed and mapped with a cross-tabulation method for the whole study period and five subperiods.
Finally, we defined the transition from cropland to steppe as abandonment, the transition from steppe to cropland as cultivation, and the transition from new-gained steppe from abandoned cropland back to cropland as recultivation. We performed the grid-based calculation for each of the five subperiods to detect if a cycle of abandonment and recultivation exists using IDRIS TerrSet software analysis tools. In order to illuminate the possible transition between steppe and cropland directly, we used trend surface analysis (TSA) in LCM of the IDRISI TerrSet software to estimate the spatial trends of the transitions from steppe to rainfed or irrigated cropland. TSA, as an interpolation technique disaggregating the broad regional patterns from the nonsystematic, fine-scale variation within the data [68,69], assumes the inherent spatial trend of the observed data can be reasonably represented by a polynomial surface of closest fit to the observed ones by minimizing the difference between the interpolated value at a data location and its original value [70]. This is calculated mathematically as: where Z stands for the transition from steppe to cropland (rainfed or irrigated). α represents the polynomial coefficients, and U and V represent the location coordinates. The TSA surfaces are calculated by assigning the pixels of a specific transition a value of 1 and pixels of no change a value of 0; they are then treated as if they were continuous values [68,69]. We performed the third-order polynomial TSA for transitions from steppe to cropland (rainfed or irrigated) to visualize the general locations of transitions in the whole study period and five subperiods. LCM in the IDRISI TerrSet software was used again to complete this calculation.

Spatio-Temporal
Change of Cropland Area at National and Regional Scales

National Level
Between 1992 and 2015, areas of rainfed cropland increased 18.59% from 333,169 km 2 to 395,093 km 2 , or an increasing rate of 18.59%. The Mann-Kendall test showed that the increasing trend was significant (p < 0.05) ( Table 2), and significant change occurred in 1999 (RSI 1999 = 2.456) (Figure 2a), whereas areas of irrigated cropland increased from 69,750 km 2 to 76,223 km 2 , or an increasing rate of 9.28%, with significant increasing trend (p < 0.05) ( Table 2) and two significant change points in 1999 and 2003 respectively (RSI 1999 = 2.319, RSI 2003 = 0.408) (Figure 2b). However, data from FAO statistics showed that the area of arable land (defined as land under temporary crops, temporary meadows for mowing or pasture, land under market and kitchen gardens and land temporarily fallow for less than five years) in Kazakhstan decreased from 350,550 km 2 to 293,950 km 2 , a decreasing rate of 16.15% at the same period, with a significant decreasing trend (p < 0.05) ( Table 2) and two significant change points in 1996 and 2000 respectively (RSI 1996 = −1.283, RSI 2000 = −2.001) (Figure 2c). The reason for this contrast trend may be the effects from the coarse resolution and inaccurate classification for the former (300 m spatial resolution) and from data sources of FAO Sustainability 2018, 10, 4042 7 of 20 estimates with official publication for the latter. Here, we considered that the former trend seems to be more reliable and could yield a more objective and reasonable trend of cropland change in Kazakhstan. former trend seems to be more reliable and could yield a more objective and reasonable trend of cropland change in Kazakhstan.

Regional Level
In the western, northern, eastern and central parts of the country, all of the rainfed cropland areas showed a significant increase (p < 0.05) ( Table 2) (Figure 3b-f,h-j) with one or several significant change points (most in 1999-2000), except in Aqtobe with insignificant increasing trend (p < 0.05) ( Table 2) as well as a significant increasing time point in 1999 and a significant decreasing point in 2009 (Figure 3a). Irrigated cropland in western and northern parts had smaller areas with less than 100 km 2 , even several km 2 , their significant changes were not tested in this study, while those in eastern and central parts showed a relatively complex condition, there was a relatively large area of irrigated cropland in East Kazakhstan with a significant increasing trend (p < 0.05) ( Table 2) (Figure 3h). An exception to Mangghystau near to the Caspian Sea in the western part is that almost no rainfed or irrigated cropland existed (Figure 3g), which could be attribute to arid climate, complex geomorphology, higher desertification risk, and the lowest population density [71][72][73].
In the southern part, due to rich water provision from glaciers and snowmelt on the north side of the Tien Shan Range [74], large irrigated oases developed. Here, the areas of irrigated cropland was more than that of rainfed cropland and showed significant increasing trends (p < 0.05) in Qyzylorda, South Kazakhstan, Zhambyl and Almaty (Table 2, Figure 3k

Net Areal Changes of Cropland and Their Contributors
Rainfed cropland was mapped, and its distribution was mainly in the northern part of Kazakhstan (with the majority in Qostanay, North Kazakhstan, Aqmola and Pavlodar), along the boundary regions with Russia (except Aqmola) and a small area in southeastern Kazakhstan, which was mosaiced with irrigated cropland (Almaty, Zhambyl and South Kazakhstan) (Figure 4). Over the whole study period , the persistent areas of rainfed and irrigated cropland were maintained as 327,180 and 68,402 km 2 , respectively, which accounted for 12.58% and 2.63% of the total territory (2.6 million km 2 ), respectively, in which the gain areas were 67,913 and 7820 km 2 , which accounted for 20.76% and 11.43% of the persistence areas of rainfed and irrigated cropland, respectively. Loss areas amount to 5989 and 1347 km 2 , which accounted for 1.83% and 1.97% of the persistence areas of rainfed and irrigated cropland, respectively (Figure 4a). This result reflects the continuous increasing trend and reduced abandonment of rainfed and irrigated cropland during this period.  In general, the area of rainfed cropland between 1992-2015 was increasing in the whole study period and increasing in the five subperiods, with most of the contribution from steppe (Figure 4).  In general, the area of rainfed cropland between 1992-2015 was increasing in the whole study period and increasing in the five subperiods, with most of the contribution from steppe (Figure 4). In addition, water bodies had a great contribution to the increase of rainfed cropland in 2005-2010 and 2010-2015 (Figure 4e,F). The situation is obviously different for irrigated cropland. During the whole study period and five subperiods, steppe was still the greatest contributor to irrigated cropland, and bareland was the second largest. Similarly to rainfed cropland, water bodies had a nonnegligible contribution to irrigated cropland in 2005-2010 and 2010-2015 (Figure 5e,f). In addition, irrigated cropland became an important provider for artificial areas since 2000 (Figure 5c-f)

Abandonment and Recultivation Cycle Analysis
For rainfed cropland, it was found that the abandoned areas were 120, 435 and 838 km 2 during the first three subperiods (SP-I, SP-II and SP-III) respectively, among which 4, 15 and 8 km 2 in SP-IV were recultivated into rainfed cropland, which accounted for 3.11, 3.45 and 0.93% of the abandoned area in the first three subperiods. 1, 12 and 4 km 2 in SP-V were recultivated into rainfed cropland, which accounted for 0.96, 2.86 and 0.43% of abandoned area in the first three subperiods. The total recultivated area in the whole study period was 27 and 17 km 2 , which accounted for 0.385% and 1.063% of the cultivated areas in next two subperiods (SP-IV and SP-V), respectively (Table 3). Obviously, the transition ratios of recultivated area to abandoned area in the first three subperiods and the cultivated areas in the next two subperiods were insignificant. We also inspected the abandonment and recultivation cycle between irrigated cropland and steppe. We determined that the abandoned areas were 11, 36, 235, 665, and 148 km 2 , and the cultivated areas were 1392, 2773, 2183, 682 and 538 km 2 in the five subsequent subperiods; however, no recultivated area was detected.

Spatial Pattern of Rainfed and Irrigated Cropland
Viewed from Figure 6, the general trend of transition from steppe to rainfed cropland in the whole study period was located in the northern and eastern parts of the study area, including North Kazakhstan, Aqumolar, Pavlodar, Qostanay, Qaraghandy and East Kazakhstan (Figure 6a). The major transitions from steppe to rainfed cropland were detected only in the northern part (North Kazakhstan, Aqumolar, Pavlodar) in SP-I (Figure 6b) and expanded to the northern and central parts (Qostanay, Qaraghandy and East Kazakhstan included) in SP-II (similar to those for the whole period) (Figure 6c). The transition expanded further to the entire eastern and central parts

Abandonment and Recultivation Cycle Analysis
For rainfed cropland, it was found that the abandoned areas were 120, 435 and 838 km 2 during the first three subperiods (SP-I, SP-II and SP-III) respectively, among which 4, 15 and 8 km 2 in SP-IV were recultivated into rainfed cropland, which accounted for 3.11, 3.45 and 0.93% of the abandoned area in the first three subperiods. 1, 12 and 4 km 2 in SP-V were recultivated into rainfed cropland, which accounted for 0.96, 2.86 and 0.43% of abandoned area in the first three subperiods. The total recultivated area in the whole study period was 27 and 17 km 2 , which accounted for 0.385% and 1.063% of the cultivated areas in next two subperiods (SP-IV and SP-V), respectively (Table 3). Obviously, the transition ratios of recultivated area to abandoned area in the first three subperiods and the cultivated areas in the next two subperiods were insignificant. We also inspected the abandonment and recultivation cycle between irrigated cropland and steppe. We determined that the abandoned areas were 11, 36, 235, 665, and 148 km 2 , and the cultivated areas were 1392, 2773, 2183, 682 and 538 km 2 in the five subsequent subperiods; however, no recultivated area was detected. Table 3. Analysis results of cycle between abandonment and recultivation during 5 subperiods (km 2 ).

Spatial Pattern of Rainfed and Irrigated Cropland
Viewed from Figure 6, the general trend of transition from steppe to rainfed cropland in the whole study period was located in the northern and eastern parts of the study area, including North Kazakhstan, Aqumolar, Pavlodar, Qostanay, Qaraghandy and East Kazakhstan (Figure 6a

Cropland Area Change in the Post-Soviet Era
The breakdown of the Soviet Union resulted in the abandonment of a large area of cropland in some post-Soviet and ex-Soviet countries in the Eurasian steppe [14,16,24,85,86]. However, as one of the primary food providers for the former Soviet Union [87,88], it was unclear whether Kazakhstan had also experienced a similar change. According to the incomplete statistics from FAO, there has been a significant decreasing trend in arable land in Kazakhstan since 1995 ( Figure 2); however, our results showed rainfed and irrigated cropland at the national level had a significant increasing trend with a significant increasing stage before 1999-2000 and a relatively stable stage during 2000-2015 (Figure 2), and the Mann-Kendall test and regime shift analysis at 14 regions also confirmed this result (Table 2, Figure 3). As an important traditional industry in Kazakhstan, agriculture had formed a relatively stable structure during the Soviet Union Era [89], which was affected very little by political surge. Due to its huge virgin land resources, especially the large area of the steppe [87,88,90], traditionally agricultural development is far easier than the development of other industries such as oil and gas [14,91]. Hence, a significant increase of cropland area in an earlier stage of newly sovereign Kazakhstan is a necessary consequence, with the coming of gradual foreign investment for resource exploitation, the agriculture in Kazakhstan lost its primary importance and cropland development showed a relatively stable trend [22,92,93].

Cropland Gains and Spatial Reorganization
While the cropland area exhibited a significant increase in Kazakhstan, it was unknown what the country-scale transition types and locations of change were. It was revealed that almost all cropland was gained from steppe in the whole study period and five subperiods; this can be viewed as a virgin land development strategy that will be continued [9,19]. The results from the trend surface analysis showed the cropland distribution was changing from a disordered expansion to ordered stability, while a small area of rainfed cropland was disappearing from the mosaic of irrigated cropland in the southern region. This may mainly result from the historical legacy of the former Soviet Union's agricultural policies [90,94], and a reorganized pattern was gradually formed as rainfed cropland in the North and irrigated cropland in the South, which reflected the implementation of an oasis agricultural development strategy in the whole country [95][96][97]. In addition, a greater contribution from water bodies to rainfed and irrigated cropland since 2005 (SP-IV and SP-V) reflected lake shrinkage to a great extent with climate change and less water from the upper reaches of many transboundary rivers in recent years [96,98,99]. Meanwhile, irrigated cropland has been an important provider for artificial areas since 2000, which reflects the rapid urbanization process at the expense of irrigated oases with good agricultural infrastructure [100][101][102][103].

Cycle of Abandonment and Recultivation between Cropland and Steppe
We inspected the transition cycle between cropland and steppe for the greatest contributor to cropland, and the results showed that there is a negligible transition from rainfed cropland to steppe then back to rainfed cropland during the five subperiods. Obviously, the soil fertility of the original steppe was higher than that of the abandoned rainfed cropland [104,105], and local farmers should prefer to explore the original steppe as rainfed cropland as opposed to recultivating abandoned or fallow cropland.
For irrigated cropland, no transition from irrigated cropland to steppe and then back to irrigated cropland was observed during the five subperiods. The main reason for this result was that the abandonment of irrigated cropland in Kazakhstan may result from the severe secondary salinization [106,107]. Hence, this abandoned land was not able to be recultivated without cost-effective desalinization measures. Obviously, the nonexistence of cycle between the abandonment and recultivation proved that newly-gained cropland from steppe may be less degraded and more productive for sustainable land use in Kazakhstan [11].

Conclusions and Implications for Sustainability
Kazakhstan's cropland development experienced a period of rapid increase and then a long-term stable period, forming a reorganization pattern with rainfed cropland in western, northern, eastern, and central parts and irrigated cropland in the southern part. The majority of new-developed cropland, rainfed or irrigated, comes from the original steppe rather than the recultivation of abandoned cropland. Rainfed cropland showed a significant increasing trend and maintained a recent long-term stable stage, dryland vulnerability to climate change can be reduced by diversified cropping systems, conservation tillage, water saving technologies, and less water-demanding crop species [108,109]. A great amount of rainfed cropland abandonment did not occur, hence we can deduce that the rainfed cropland system under the current situation has high sustainability. For irrigated cropland, the main impact on its areas and productivity came from water provision capacity of inland rivers; the amount of water was determined by precipitation, glacier-and snow-melt and the upstream water utilization amount. The former two were affected by climate change, the latter one affected by transboundary cooperation [96]. Similar to rainfed cropland, irrigated cropland also showed a significant increasing trend and maintained a recent long-term stable stage, a great amount of irrigated cropland abandonment did not occur, hence we can deduce that the cropland system under the current situation has medium sustainability limited mainly by transboundary water allocation. Undoubtedly, the results of the analysis can provide strong evidence for sustainable land use and a basis for food security policy-making in Kazakhstan and even all of Central Asia in the future [110]. In addition. It is undeniable that some significant changes at the regional scale may be neglected or omitted using ESA-CCI land cover dataset with an annual temporal resolution and 300 m spatial resolution; hence, some satellite images with high temporal resolution, such as MODIS (Moderate Resolution Imaging Spectroradiometer), and high spatial resolution, such as Landsate TM (Thematic Mapper), can be used to discover the regional changes in croplands in Kazakhstan when combined with several popular software packages, such as Landtrendr [26,111] and BFAST [112,113], which can characterize land use changes and their causes more directly and concretely.