Study on the Classiﬁcation and Change Detection Methods of Drylands in Arid and Semi-Arid Regions

: The aim of this study was to clarify the distribution of irrigated drylands in arid and semi-arid areas, where complex terrain, diverse crops and staggered cultivated lands exist. This paper studied the classiﬁcation methods of irrigated drylands based on temperature, precipitation, Normalized Difference Vegetation Index (NDVI) and Normalized Difference Water Index (NDWI) from Landsat data in the one-harvest area of the northern Loess Plateau of China by using the Google Earth Engine (GEE) platform. An extraction method was proposed for irrigated drylands in arid and semi-arid regions of northwest China. In addition, the change types of irrigated and rainfed drylands in the two periods were classiﬁed, and a method was also put forward to directly classify the change types by using the image differences between the two periods combined with the classiﬁcation results of each period. It was found that combining the ratio of NDVI and NDWI with the accumulated values of temperature and precipitation of the 30 days before imaging could effectively improve the classiﬁcation accuracy. Moreover, directly classifying the gaps of remote sensing factors in the time dimension before combining spatial clustering information could yield a more accurate type of change, because the accumulation of errors in the change maps obtained from the overlay analysis of distribution maps of the two periods could be avoided. The accuracy of classiﬁcation could be improved by introducing the dynamic information of time dimension into the classiﬁcation of historical periods. This study complements the extraction method for this type of irrigated dryland, and the classiﬁcation results can improve the accuracy of existing products in terms of spatial resolution, which can ﬁll the shortage of detailed distribution data for irrigated and rainfed drylands in this region. four pieces of data, including the cumulant of 15 days before imaging, the cumulant of the 30 days before imaging, the cumulant of 60 days before imaging, and the cumulant from April 15th (sowing time) to the imaging day, for classification. The analysis showed that the NDVI and NDWI of irrigated drylands and rainfed drylands have a different relationship with precipitation and temperature. In addition, the scatter strongly confirmed the positive correlation between precipitation and NDVI and NDWI. The slope of each linear equation represented in a scatterplot regarding pre- cipitation was positive. In contrast, temperature was negatively correlated with NDVI and NDWI, and the slope of each linear equation represented in the scatterplot regarding pre- cipitation was negative. Critical Phases In July and August, the heat conditions are the best and the crops grow fastest. Most crops reach maturity and peak NDVI in August. In July and August, prolonged drought stress has a concentrated effect on the vegetation growth and water content. Some crops trends at that


Introduction
Irrigated cropland, which accounts for 18% of the world's arable land, produces about 40% of the world's total grain output, and the grain yield per unit area of irrigated cropland is more than twice that of dry croplands [1][2][3]. This gap is particularly dramatic in arid and semi-arid regions. Irrigation allows more crop selectivity in this region, so that more profitable crops are usually grown. According to surveys, in normal harvest years, the yield of naked oats can reach 750−1500 kg per hectare, while that of corn can reach 15,000 kg after irrigation. In the vast arid and semi-arid regions of northwest China, irrigated farmland is the main source of grain production [4,5]. Therefore, developing irrigated cropland not only has a huge impact on agricultural production and product types, but also has a special impact on the economic benefits of farmers and the rural economy. Due to different crop types and crop management, the effects of irrigated and rainfed croplands on the ecological environment are quite different. Agricultural irrigation in arid and semi-arid regions consumes more than 50% of the total water supply [2,6] which exacerbates local water shortages, leads to a reduction in groundwater, and even affects the climate and ecological systems. In addition, the average amount of fertilizer applied to irrigated farmlands is about six times that of rainfed farmlands. The large amount of fertilizer applied with irrigation water increases the risk of agricultural nonpoint source pollution. Irrigation, especially flood and overirrigation, further aggravates the soil erosion of cultivated land [7,8]. Due to the difference in yields between irrigated and rainfed croplands, irrigation facilities have been widely built in this region to increase crop yields and thus improve the living standards of farmers. However, the region is water scarce and its Standardized Precipitation Evapotranspiration Index (SPEI) ranges between −2 and 0 [9][10][11]. Desertification and land degradation have been exacerbated by the fragile ecology because of the decline in groundwater level caused by excessive water use. The contradictions between increasing farmers' income and lack of water resources, and between economic development and ecological environment, have caused great concern worldwide. Clarifying the quantity and spatial pattern of irrigated and rainfed croplands can provide data to balance the contradictions.
In addition to the sophisticated hydrothermal and soil conditions, the northern part of the Loess Plateau of China has a unique spatial pattern of agriculture. Irrigated drylands are found both on contiguous plains and on scattered terraces and slopes. A wide variety of crops are grown, including some intercropping types. Due to the small-scale farmer economy with the family as the basic unit and flexible crop cultivation practices, the spatial pattern of agriculture in the research area is variable rather than fixed. Therefore, the rapid and timely identification of changes in the irrigated drylands in the study area is considered as a great challenge and of great practical value for improving crop yields and protecting the environment. Previous studies include irrigated cropland distribution purely derived from remote sensing data, or by fusing remote sensing data and statistical data. The earliest large-scale distribution map of irrigated cropland was generated by spatialized statistical data. Such data always lacked real-time information and accurate spatial information. Some studies aimed to synthesize various data to obtain more accurate datasets [12][13][14], which showed higher consistency in large irrigated areas and lower consistency in small irrigated fields.
In recent years, supervised classification methods using machine learning techniques, such as random forest and decision tree, have been used to extract irrigated and rainfed croplands. Many of these studies have been conducted at a large scale and low spatial resolution. For example, GRIPC (global rainfed, irrigated, and paddy croplands) [15] used MODIS (Moderate Resolution Imaging Spectroradiometer) data with a spatial resolution of 500 m, the global database of Nagaraj et al. with a resolution of 9 km [16], and GMIA (Global Map of Irrigation Areas) with a resolution of 5 arcmin [17,18]. In addition, most of these researchers focused on irrigated croplands, including paddy fields, rather than irrigated drylands (not including paddy fields). Paddy fields are usually distributed in areas with sufficient water resources (southeast China). Irrigation poses relatively less threat to water resources, while irrigated drylands, especially in arid areas, pose a greater threat to the environment. Therefore, research on irrigated drylands has received more and more attention. Samasse et al. extracted irrigated croplands and rainfed croplands in parts of Africa at a resolution of 30 m [19]. Liu Liangyun developed a global landuse database with a resolution of 30 m, including irrigated croplands and non-irrigated croplands [20]. Although these studies have a higher spatial resolution than before, the samples were mainly selected from representative agricultural areas around the world. With a small number of samples of irrigated drylands, there were not enough studies on the irrigated and rainfed drylands in arid and semi-arid regions. It was also noted that Facundo Rojas [21] studied the spatial-temporal changes of irrigated drylands, but most of his study area was irrigated drylands. Moreover, there was no staggered distribution with rainfed drylands. The above characteristics also appear in studies in China. They focus on the extraction of large-scale irrigated areas, rather than small irrigated fields. For example, Xiufang Zhu obtained the distribution accuracy of irrigated croplands in China in 2000, but the resolution was only 10 km and 1 km [22]. Afterwards, Dong Tingting, Liu Yizhu and Kunlun Xiang did this work using MODIS, and increased the resolution to 250 m, which was still insufficient [1,23,24]. In Kunlun Xiang's research on irrigated farmlands in northeast China, irrigation was determined by comparing the surface humidity between farmlands and the surrounding forest. However, this study also focused on paddy fields. Irrigated drylands do not have such distinct characteristics, and there are not enough forests in the northwest for comparison [25]. There are some data regarding the distribution of irrigated croplands in existing studies, but they cannot achieve high agreement, especially in west China [13]. It was difficult to obtain data with high accuracy and reliability in west China, because current research generally focuses on typical irrigated areas, such as the North China Plain and Sichuan Basin. Although they are nationwide studies, less attention is paid to the arid and semi-arid regions. As a result, these regions lack reliable and detailed data. Irrigated drylands have received attention in some parts of China, but they always have limitations. Some scholars have studied extraction methods using relatively high value durations of NDVI and NDWI and spatial-temporal changes of irrigated and rainfed drylands in Zhangjiakou, but the method had limitations and cannot be directly used elsewhere because it only considered the local environment [26].
For the studies of land cover and changes, the overlay analysis of multi-period status quo is mostly adopted. This method is simple and easy to use, but it may accumulate errors in the status quo extracted from two periods. Due to the infeasibility of sampling in the field and inaccurate memory in questionnaires, the lack of samples can lead to more errors in historical years. For arid and semi-arid regions in northwest China, the classification accuracies of irrigated and rainfed drylands are not as high as those in large irrigated areas. The fragmentation and staggered distribution of irrigated and rainfed drylands increased the difficulty of classification. Moreover, the characteristics of remote sensing factors of irrigated drylands are not as obvious as those of paddy fields. Therefore, the overlay analysis of two periods would bring more errors. So, the change types of irrigated and rainfed drylands need to be studied further, including lands with no change, noncultivated lands converted to irrigated drylands, non-cultivated lands converted to rainfed drylands, irrigated drylands converted to rainfed drylands, rainfed drylands converted to non-cultivated lands, irrigated drylands converted to non-cultivated lands and rainfed drylands converted to irrigated drylands.
There is no specific research on the distribution of irrigated and rainfed drylands in the arid and semi-arid regions of northwest China. At present, existing databases including this region always fail to extract small patches of irrigated dryland effectively. However, there are lots of irrigated drylands of a small size because of the terrain and small-scale family farming, resulting in irrigated drylands being reflected insufficiently on maps [15]. Moreover, there are few targeted studies on the spatial-temporal changes of irrigated drylands in this area. However, irrigated drylands and their changes cannot be ignored for the local food production and ecological environment. This paper aims to study a classification method and a change detection method for irrigated and rainfed drylands suitable for arid and semi-arid regions. The one-harvest area in the northern part of the Loess Plateau was selected as the study area, and the GEE platform was used to implement the research. Irrigated drylands were separated from rainfed drylands in 2015 by using the multi-Temporal data of Landsat with a resolution of 30 meters, combined with precipitation and temperature data during the growing season. Based on the differences between the data of 2005 and 2015, and combined with the classifications in 2005 and 2015, change types can be classified regarding temporal and spatial differences. This paper can provide inspiration for the classification and change detection of irrigated and rainfed drylands in similar areas. The distribution pattern and changes of irrigated drylands in the one-harvest Remote Sens. 2022, 14, 1256 4 of 21 area of the northern Loess Plateau of China are provided, which can support irrigation policies, water conservancy projects and agricultural policies in this region.

Study Area
In this paper, the northern part of the Loess Plateau of China was selected as the study area (Figure 1), including the south of Bayannur and Baotou, Hohhot and Ordos in Inner Mongolia, Shuozhou, Xinzhou, Luliang and Taiyuan in Shanxi Province, and Yulin in Shaanxi Province. Located in the transition zone from farming to grazing and from semi-arid to arid, this area presents a rich variety of land cover types (e.g., grassland, desert, farmland, forest and bare soil) [27]. The topography of the northern Loess Plateau is also quite complex, with the elevations ranging from 800 to 1200 m. In particular, the mountains (e.g., Taihang Mountains, Yin Mountains and Lvliang Mountains), rivers (e.g., Fen River, Hutuo River, Sanggan River and Heihe River) and basins (e.g., Taiyuan Basin, Xinding Basin and Datong Basin) synergize to form a unique spatial pattern of agriculture. The growing season starts in mid-April and ends in mid-October, with only one harvest [28]. A variety of crops are grown, including corn, wheat, potatoes, naked oats, flax, millet, sorghum, buckwheat, rice, various beans and forage grass. Moreover, there are some intercropping types such as apricot trees, other fruit trees, walnuts or jujubes mixed with food or vegetable crops. There is a relatively low annual precipitation of 150-500 mm on average. In response to the lack of water resources, irrigation has become one of the vital auxiliary means for farming. The irrigated drylands in this region are distributed over continuous plains and scattered terraces and slopes. Spatially, the northwest of the study area generally has lower precipitation compared to the southeast, and is therefore accompanied by more irrigated drylands due to the greater need for water replenishment [29][30][31]. Along with the population growth and economic development in the study area in recent years, the issue of increasing food production and security here has attracted widespread attention. In addition to natural rivers and lakes, the groundwater has also been inevitably pumped away for irrigating crops, which leads to the over-exploitation of groundwater, and seriously threatens the local ecological balance. Subsequently, the conflict between supply and demand of water resources has further intensified. In addition, the soil structure of the Loess Plateau is too loose to resist frequent irrigation. Therefore, soil erosion is becoming a potential crisis in the study area [32].

Data Introduction
The Google Earth Engine (GEE) platform was used in this study because of its rich data, diverse algorithms and strong calculation capacity. With long time series and midto high-resolution data available on GEE, its fast calculation capacity can help us obtain

Data Introduction
The Google Earth Engine (GEE) platform was used in this study because of its rich data, diverse algorithms and strong calculation capacity. With long time series and midto high-resolution data available on GEE, its fast calculation capacity can help us obtain results in a short time frame.
Landsat imagery is the primary dataset used in this study. In the arid and semi-arid region of northwest China, irrigated and rainfed drylands are interlaced. Considering the growth rate of local crops and small plots of cultivated land in irrigated drylands, Landsat images with a temporal resolution of half a month and a geographical resolution of 30 m were selected as the major datasets. Landsat 8-Day NDVI Composite (LC8_L1T_8DAY_NDVI of GEE) and Landsat 8-Day NDWI Composite (LC8_L1T_8DAY_NDWI of GEE) were used [33]. In this paper, these data were synthesized into a semimonthly period. The first half of a month was marked as -1 and the second half of a month was marked as −2, e.g., 5−2 means the second half of May. In order to reduce the effect of data quality (cloud occlusion, These databases have been calibrated radiometrically, and they can provide comparable NDVI and NDWI [33]. Finally, the semi-monthly NDVI and NDWI data of the growing season (May-September) were formed.
The formulas of NDVI and NDWI are as follows [34,35]: where N IR is the surface reflectance of near infrared band; RED is the surface reflectance of red band; and SW IR is the surface reflectance of short-wave infrared band. The MYD11A1 V6 product was used. It provides the daily land surface temperature (LST) within the kilometer grid. The temperature value was derived from the MOD11_L2 swath product [36]. The PERSIANN-CDR was used too [37,38]. It is a daily quasi-global precipitation product that spans the period from 1 January 1983, to the present. The data have good temporal resolution and high space resolution sufficient to represent the spatial differences in temperature [36] and precipitation [37,38].
Land cover data were the basis of this study, and irrigated and rainfed drylands were further extracted from the farmland category. The land cover data were extracted from the national land-use/cover database of China (hereinafter referred to as NLUD-C) at 1:100,000 scale that contains Chinese land-use/cover data of six periods (the 1980s, 1995, 2000, 2005, 2010 and 2015). It was visually interpreted from Landsat data by Chinese Academy of Sciences, and the accuracy was above 95% [39].
The proportions of irrigated and rainfed drylands in each region (county and city) were calculated by using the area of cultivated land, irrigated dryland, rainfed dryland, paddy field and the area of effective irrigated cropland (the area of effective irrigated cropland is the sum of the area of paddy field and irrigated dryland) in the statistical yearbook of each region [40][41][42][43][44][45][46][47][48][49]. These data were used to help select samples and verify the classification results.
Samples ( Figure 1) in 2015 were obtained through field surveys, questionnaires, and visual interpretations using Google Earth images. The questionnaire survey included locations, proportion of irrigated drylands and frequency of irrigation. The texture features and location distribution of samples from field and questionnaire surveys were analyzed to help obtain more samples through visual interpretation from Google Images. Other samples were selected according to the proportion of irrigated or rainfed drylands in the statistical data. For example, if irrigated or rainfed drylands accounted for more than 95% in an administrative unit, characteristics of the main cropland in the region were determined as irrigated or rainfed drylands combined with Google Images.
In order to improve the evaluation of methods and avoid imbalances in the number of category validations, the number of validation samples was increased, especially in irrigated drylands. In 2015, a total of 77 irrigated drylands and 212 rainfed drylands samples were obtained for training, and 180 irrigated drylands and 278 rainfed drylands samples were obtained for verification. However, the sample of 2005 could not be obtained from the field survey and questionnaire survey, and only a small number of historical images are available on Google Earth. Therefore, it was largely be obtained through statistical data from areas with relatively homogeneous types of cultivated land. In 2005, 91 irrigated drylands and 128 rainfed drylands were selected as training samples. The validation samples included 109 irrigated drylands and 103 rainfed drylands. Among them, 44 and 67 samples of irrigated drylands and rainfed drylands did not change, respectively, which were mainly used to determine the floating range of NDVI and NDWI from 2005 to 2015, and 113 samples with changes were used to verify the accuracy of the classification of change types.

Methodology
In this study, irrigated and rainfed drylands were further extracted from the farmlands. Dryland areas were clipped from land cover data in 2015 and 2005, respectively.

Feature Selection
The diverse environments within the study area had a great impact on the effectiveness of irrigation. For example, water available for vegetation growth of irrigated drylands in arid areas may not be better than rainfed drylands in relatively humid areas. In warmer areas, vegetation may not grow well because of intensified drought. In various studies, remote sensing factors and environmental factors are used to differentiate between irrigated and rainfed farmlands. Salmon et al. used remote sensing, climate and agricultural data to generate global maps of irrigation, rainfed and paddy fields through supervised classification. [15]. Nagaraj et al. [16] applied artificial intelligence methods to produce global irrigated areas using remote sensing and monthly climate data from 2001 to 2015. Although the study area in this paper was not as large as the above study, this study also combined the remote sensing and environmental factors as classification features because of the complexity of terrain, diversity of environment, and variety of irrigated dryland types.

Features by Remote Sensing
Because water is the greatest restraint for the growth of crops in the arid and semi-arid regions of Northwest China, the difference in soil water content is largely reflected by the growth of vegetation [50]. Therefore, vegetation growth rather than soil moisture (which is blocked by vegetation and changes too quickly to be captured by remote sensing) was chosen to differentiate irrigation. NDVI is regarded as the most common and stable index of vegetation growth, and NDWI is the water index that best reflects vegetation water content [34]. When the vegetation canopy is under water stress, NDWI can respond in time, which is of great significance in drought monitoring [35,51]. Therefore, these two factors were used as remote sensing factors.

Environmental Factors
NDVI and NDWI can reflect the moisture status of vegetation, but cannot identify the contribution from precipitation or irrigation, especially in areas with large differences in precipitation. Therefore, precipitation is used as a classification factor to reduce the error caused by different precipitation. NDVI and NDWI of rainfed drylands may be positively correlated with precipitation because they are fully dependent on precipitation. However, this relationship is not obvious for irrigated drylands. Case studies have shown that EVI anomalies, which can be obtained from two steps-first, subtracting the seasonal cycle and then the long-term trend from the raw time series of EVI [52]-are significantly correlated with drought anomalies in arid and semi-arid rainfed drylands. Irrigation can well alleviate the adverse impact of drought on irrigated drylands [52,53]. The contribution of unit precipitation to NDVI is reflected by the ratio of NDVI and NDWI to precipitation. On the one hand, the NDVI of rainfed drylands is relatively low due to the stress of water scarcity, so the ratio of NDVI to precipitation is lower. On the other hand, according to the existing studies, irrigated drylands have higher utilization efficiency of precipitation than rainfed drylands, because irrigated drylands are largely more intensively cultivated and more fertilized. Therefore, the soil comprehensive situation of water, fertilizer, and dissolved oxygen and heat in irrigated drylands is better than that of rainfed drylands [54,55]. Dengke et al (2015) [56] conducted a detailed study on the water-use efficiency (WUE) of crops in rainfed drylands in north China. The results showed that intensive tillage could significantly improve crop growth and WUE, and the crop yield and WUE were significantly improved after irrigation compared with no irrigation [56]. Therefore, the ratio of NDVI and NDWI to precipitation can distinguish irrigated drylands from rainfed drylands.
Similarly, the higher the temperature, the stronger the evaporation, and the water content of soil and vegetation will be affected [57]. Under drought stress in rainfed drylands, the higher the temperature, the lower the water content, and the worse the vegetation growth [58]. In contrast, in irrigated drylands without drought stress, this relationship is not obvious, and generally in this area the higher the temperature the better the growth, due to better heat. Therefore, NDVI and NDWI may be negatively correlated with temperature, and the ratio of NDVI and NDWI to temperature is used to differentiate irrigated drylands with rainfed drylands. In total, the scatter plots of NDVI, NDWI, temperature and precipitation were analyzed. For precipitation and temperature, we used four pieces of data, including the cumulant of 15 days before imaging, the cumulant of the 30 days before imaging, the cumulant of 60 days before imaging, and the cumulant from 15 April (sowing time) to the imaging day, for classification.
The analysis showed that the NDVI and NDWI of irrigated drylands and rainfed drylands have a different relationship with precipitation and temperature. In addition, the scatter strongly confirmed the positive correlation between precipitation and NDVI and NDWI. The slope of each linear equation represented in a scatterplot regarding precipitation was positive. In contrast, temperature was negatively correlated with NDVI and NDWI, and the slope of each linear equation represented in the scatterplot regarding precipitation was negative.

Critical Phases
In July and August, the heat conditions are the best and the crops grow fastest. Most crops reach maturity and peak NDVI in August. In July and August, prolonged drought stress has a concentrated effect on the vegetation growth and water content. Some crops enter the harvest stage in September, and the NDVI and NDWI of crops exhibit downward trends at that point. Moreover, according to the scatterplots (Figure 2), irrigated drylands and rainfed drylands have a good distinction from the second half of June to the second half of August, while the relationship is not obvious in other periods. In order to select the critical classification phases, the data combinations of different phases were classified and compared. For comparative analysis, they were classified and validated using the same samples and methods.
enter the harvest stage in September, and the NDVI and NDWI of crops exhibit downward trends at that point. Moreover, according to the scatterplots (Figure 2), irrigated drylands and rainfed drylands have a good distinction from the second half of June to the second half of August, while the relationship is not obvious in other periods. In order to select the critical classification phases, the data combinations of different phases were classified and compared. For comparative analysis, they were classified and validated using the same samples and methods.

Classifier
We used random forest as classifier, which is a type of machine learning. Machine learning is a very common and effective classification method, and random forest is widely used in land cover classification and irrigated cropland identification [59]. Random forest is a combination of tree-like predictors. Each tree depends on the value of a random vector sampled independently, and all trees in the forest have the same distribution [60]. Random forest can improve the prediction accuracy without significantly increasing the computation time. The results are robust to missing and unbalanced data,

Classifier
We used random forest as classifier, which is a type of machine learning. Machine learning is a very common and effective classification method, and random forest is widely used in land cover classification and irrigated cropland identification [59]. Random forest is a combination of tree-like predictors. Each tree depends on the value of a random vector sampled independently, and all trees in the forest have the same distribution [60]. Random forest can improve the prediction accuracy without significantly increasing the computation time. The results are robust to missing and unbalanced data, and can well predict the effects of many explanatory variables [60]. Random forest can explain the contribution of each classification factor, which also helps us to refine the feature factors and optimize classification methods. We chose the ee.Classifier.SmileRandomForest algorithm of the GEE platform and set the parameter to 10 trees according to the experiments in this study, because it shows that 10, 20 and 30 trees have a similar effect on the results, respectively.
In the Importance of Features Table (Table 1) from the explanation of the classifier, each feature had periods of relatively high importance. Moreover, each period had features of relatively high importance. There was no obvious characteristic or rule. Therefore, it is considered that the features and phases selected in classification were all important.

Change Types Detection Method
The normal approach to obtain the spatiotemporal dynamics of land-use types is to perform overlay analysis. However, errors of the classification results of different years would be accumulated inevitably, resulting in large errors in the dynamics. In order to reduce such errors, the changes of NDVI and NDWI were classified firstly, and then the classification results of 2015 and 2005 were combined to obtain different types of change, containing no change, non-cultivated lands to irrigated drylands, non-cultivated lands to rainfed drylands, irrigated drylands to rainfed drylands, rainfed drylands to non-cultivated lands, irrigated drylands to non-cultivated lands and rainfed drylands to irrigated drylands. Specific steps are as follows.
Step 1. The classification of 2015 was used as the benchmark. The period of 2015 was considered the best, because the data were the most complete and the sample number was the largest; Step 2. The mean values of NDVI and NDWI for each half month in 2015 were subtracted from those in 2005 to obtain the gaps between the two periods; Step 3. The floating range of NDVI and NDWI in the unchanged area was obtained by the statistical analysis of the unchanged samples; Step 4. The types of change were classified according to the gaps between the two periods.

Verification Method
In this study, three methods were used to verify the precision of the extracted irrigated drylands. Firstly, the results were verified by using a confusion matrix of the verification samples. Secondly, they were verified by statistical data. We used the proportion of the area of the county-level irrigated drylands to compare with the results. They were also evaluated by STDEV (standard deviation), R 2 (coefficient of determination) and RMSE (root mean squared error). Thirdly, they were compared with existing research results.

Results
The random forest classification experiments of 10 trees were carried out by combining the following remote sensing data, environmental data and phases. Confusion matrix verification, statistical data comparison and comparison with existing results were used for validation. Table 2 shows that a good result could be obtained by using only NDVI and NDWI of the whole growing season (10 periods) in the above classification. The overall accuracy and Kappa coefficient were 0.824 and 0.632, respectively. When temperature and precipitation were added, the accuracy reached 0.857 and 0.697, respectively. The best results appeared in the cumulative values of temperature and LST of the 30 days before imaging. Therefore, this parameter will be used in next experiments. In the classification of each critical period, experiments were conducted mainly for the characteristic period from June to August selected in Figure 2. The overall accuracy was 0.847 and the Kappa coefficient was 0.684, which was close to the classification accuracy of the whole growing season. Then, data of June, July and August were used for classification separately. It can be seen that July and August had higher accuracy and were the most discriminative characteristic periods. Therefore, the overall accuracy and Kappa coefficient were higher using data from four periods in July and August and data of 7-2 and 8-1 (two periods), reaching 0.851, 0.687 and 0.860, 0.710, respectively.  The rest were classified as follows:

Confusion Matrix Verification
If the gaps of the two indicators were both positive, it was considered that the moisture situation had become better. Therefore, if the pixels were rainfed drylands in 2015, the year 2005 must be rainfed drylands. If the pixels were irrigated drylands in 2015, they would be judged by the classification of 2005.
If the gaps of the two indicators were both negative, it was considered that the moisture situation had become worse. Therefore, if the pixels were irrigated drylands in 2015, then they must be irrigated drylands in the year 2005. If the pixels were rainfed drylands in 2015, they would be judged by the classification of 2005.
If the gaps of the NDVI and NDWI were one positive and one negative, this could not be identified from the time axis. Then, they would be judged by the classification of 2005.

Verification Method
In this study, three methods were used to verify the precision of the extracted irrigated drylands. Firstly, the results were verified by using a confusion matrix of the verification samples. Secondly, they were verified by statistical data. We used the proportion of the area of the county-level irrigated drylands to compare with the results. They were also evaluated by STDEV (standard deviation), R 2 (coefficient of determination) and RMSE (root mean squared error). Thirdly, they were compared with existing research results.

Results
The random forest classification experiments of 10 trees were carried out by combining the following remote sensing data, environmental data and phases. Confusion matrix verification, statistical data comparison and comparison with existing results were used for validation. Table 2 shows that a good result could be obtained by using only NDVI and NDWI of the whole growing season (10 periods) in the above classification. The overall accuracy and Kappa coefficient were 0.824 and 0.632, respectively. When temperature and precipitation were added, the accuracy reached 0.857 and 0.697, respectively. The best results appeared in the cumulative values of temperature and LST of the 30 days before imaging. Therefore, this parameter will be used in next experiments. In the classification of each critical period, experiments were conducted mainly for the characteristic period from June to August selected in Figure 2. The overall accuracy was 0.847 and the Kappa coefficient was 0.684, which was close to the classification accuracy of the whole growing season. Then, data of June, July and August were used for classification separately. It can be seen that July and August had higher accuracy and were the most discriminative characteristic periods. Therefore, the overall accuracy and Kappa coefficient were higher using data from four periods in July and August and data of 7-2 and 8-1 (two periods), reaching 0.851, 0.687 and 0.860, 0.710, respectively. In comparison, all the classification methods with temperature and precipitation were more accurate than those without temperature and precipitation. The overall accuracy was increased by 4.03% on average, while the Kappa coefficient was increased by 8.54% on average, indicating that the classification accuracy of irrigated drylands was improved after adding environmental factors. Moreover, the accuracies of irrigated drylands and rainfed drylands also led to a significant increase in Kappa coefficient.

Confusion Matrix Verification
In general, the three classification methods using the ratio of NDVI and NDWI to temperature and precipitation from 6-2 to 8-1, 7-1 to 8-2 (4 phases) and 7-2 to 8-1 (2 phases) were better, and the classification accuracies were close to that of the whole growing season. The classification accuracies of using only NDVI and NDWI in critical phases were worse than that of the whole growing season. It can be seen that adding precipitation and temperature can be used for several phases and still maintain a good classification result.

Comparison of Classification Results and Statistical Data
The results of 6, 8 and 16 classifications with high accuracy in the confusion matrix were compared with the statistical results in Figure 4. The proportions of irrigated drylands in non-paddy fields in the counties of Bayannaoer in Inner Mongolia, Xinzhou in Shanxi Province, and Yulin in Shaanxi Province were selected for comparison. Bayannur, located in the northwest of the study area, is a typical Yellow River irrigation area, of which more than 97% is irrigated drylands. Xinzhou, located in the eastern part of the study area, has both small inter-mountain basins with irrigated drylands and rainfed drylands on the surrounding upland. Yulin, located in the southwest, has typical irrigated drylands developed on wind-sand beaches. The sources of statistical data were the Inner Mongolia Statistical Yearbook 2016, Yulin Statistical Yearbook 2017, and the third agricultural census (Shanxi Province) [40,41,[45][46][47]. The year of data survey was 2015 or 2016. result.

Comparison of Classification Results and Statistical Data
The results of 6, 8 and 16 classifications with high accuracy in the confusion matrix were compared with the statistical results in Figure 4. The proportions of irrigated drylands in non-paddy fields in the counties of Bayannaoer in Inner Mongolia, Xinzhou in Shanxi Province, and Yulin in Shaanxi Province were selected for comparison. Bayannur, located in the northwest of the study area, is a typical Yellow River irrigation area, of which more than 97% is irrigated drylands. Xinzhou, located in the eastern part of the study area, has both small inter-mountain basins with irrigated drylands and rainfed drylands on the surrounding upland. Yulin, located in the southwest, has typical irrigated drylands developed on wind-sand beaches. The sources of statistical data were the Inner Mongolia Statistical Yearbook 2016, Yulin Statistical Yearbook 2017, and the third agricultural census (Shanxi Province) [40,41,[45][46][47]. The year of data survey was 2015 or 2016.
From the perspective of districts and counties (Figure 4), the classifications of all three methods maintained high consistency with the statistical data, with the fitted R 2 reaching 0.9 and the RMSE around 0.1. The proportion of misclassification compared with statistical data were mostly kept within 10%. However, it was noticed that a large gap of 30% occurred in Wutai and Yuyang. For this reason, these two areas were specially analyzed. According to the classification results, the irrigated drylands in Wutai County are mainly distributed in Dongye Town, Taicheng Town and Doucun Town. There are low and flat river valleys areas,  From the perspective of districts and counties (Figure 4), the classifications of all three methods maintained high consistency with the statistical data, with the fitted R 2 reaching 0.9 and the RMSE around 0.1. The proportion of misclassification compared with statistical data were mostly kept within 10%. However, it was noticed that a large gap of 30% occurred in Wutai and Yuyang.
For this reason, these two areas were specially analyzed. According to the classification results, the irrigated drylands in Wutai County are mainly distributed in Dongye Town, Taicheng Town and Doucun Town. There are low and flat river valleys areas, which are rare in Wutai County. Furthermore, it is confirmed by several documents that there are several irrigation projects in these small river valleys, including Xiaoyin River and Lvsi River irrigation areas [61][62][63][64]. In Yuyang District Statistical Yearbook (2017) [40], the area of arable land was 69, 540 hectares at the beginning of the year, and 85, 314 hectares at the end of the year. The gap reached 22.68%. Therefore, the results of these two counties were not considered for evaluation.
From the scale of prefecture-level cities, the misclassified area was 14.93% at most, which was Ordos City. On the whole, mountainous and desert areas are prone to misclassification, for example, Wutai County, located in the Wutai Mountains, and Ordos and Yulin located in Mu Us Desert. The complex topography of the mountainous areas and the more fragmented plots of cropland result in more mixed pixels. These cropland pixels mixed with forest or grassland are more difficult to identify. Here, the probability of misclassifying rainfed drylands as irrigated drylands is slightly higher. It is also partly due to the orchards. At the edge of the desert, on the one hand, irrigated drylands may be classified as rainfed drylands due to insufficient and unstable irrigation water sources, while some statistics classify the unguaranteed fields as irrigated drylands; on the other hand, changes in croplands are very common in the region, and some statistical data and literature data also vary greatly. Due to the inconsistency of the statistical standard, it is normal to have some discrepancies with statistical data. In general, the distribution trend of irrigated and rainfed drylands was consistent with the statistical data.

4.3.Distribution and Comparison with Existing Data
Compared with GMIA [18], the overall distribution trend of irrigated drylands was consistent. As can be seen from Figure 8 and Figure 5, irrigated drylands in the study area were mainly distributed in small plains with lower elevations, including Hetao Plain in the northwest, Tumochuan Plain in the north, as well as Datong Basin, Xinding Basin and

Distribution and Comparison with Existing Data
Compared with GMIA [18], the overall distribution trend of irrigated drylands was consistent. As can be seen from Figures 5 and 8, irrigated drylands in the study area were mainly distributed in small plains with lower elevations, including Hetao Plain in the northwest, Tumochuan Plain in the north, as well as Datong Basin, Xinding Basin and Taiyuan Basin in the east. In the southern edge of Mu Us Desert in the west and south, there are less cultivated lands, but they are mainly irrigated drylands. The Loess hilly area in the middle is dominated by rainfed drylands, with only a small number of irrigated drylands found in valleys near rivers.

Changes Distribution of Irrigated and Rainfed Drylands
As Table 3 and Figure 7 show, the transformation of large areas of cultivated land into irrigated drylands occurs mainly in the west and northwest, such asOrdos (11,425.43 km 2 ), and Bayannaoer (6127.62 km 2 ), where some arable lands are reclaimed following the development of irrigation facilities. The transformation of irrigated drylands into nonarable lands is very rare, and is mainly concentrated in the northwest of the sandy beach, where arable lands have been abandoned because water sources cannot be guaranteed. The mutual transformation between rainfed drylands and non-arable lands occurs mainly in the mountainous areas in the central and southern parts of the study area, such as Yulin (29,332.19 km 2 ) and Lvliang (15,607.43 km 2 ), where farming is difficult and less profitable. Farmland is often abandoned, but not perpetually. When farmers are able to cultivate the lands or have no other source of income (e.g., go out to work in cities), they will cultivate them again. The transformation from irrigated into rainfed drylands occurs mainly in the east and southwest, such as in Yulin (11,748.95 km 2 ) and Xinzhou (5032.12 km 2 ). Some irrigated drylands can no longer be effectively irrigated due to the disrepair of the engineering facilities in the old irrigated areas, coupled with groundwater recession and increasing irrigation costs.  GMIA has few irrigated drylands in the southwest, few cultivated lands in some areas in the south, and a small number of irrigated drylands in the central and southern Lvliang Mountains. This is inconsistent with the results of this study. According to the results of this study, there is a small number of cultivated lands in the southern edge of the Mu Us Desert in the west and south, but irrigated drylands are dominant, with a small portion of rainfed drylands. The central Loess hilly region is dominated by rainfed drylands, with only a few irrigated drylands in river valleys, and a large number of cultivated lands in the mountainous areas of the southern Loess Plateau. However, because the cultivated lands in these areas are relatively fragmented, they cannot be reflected in the classification map with low spatial ratio.

Changes Distribution of Irrigated and Rainfed Drylands
As Table 3 and Figure 7 show, the transformation of large areas of cultivated land into irrigated drylands occurs mainly in the west and northwest, such asOrdos (11,425.43 km 2 ), and Bayannaoer (6127.62 km 2 ), where some arable lands are reclaimed following the development of irrigation facilities. The transformation of irrigated drylands into nonarable lands is very rare, and is mainly concentrated in the northwest of the sandy beach, where arable lands have been abandoned because water sources cannot be guaranteed. The mutual transformation between rainfed drylands and non-arable lands occurs mainly in the mountainous areas in the central and southern parts of the study area, such as Yulin (29,332.19 km 2 ) and Lvliang (15,607.43 km 2 ), where farming is difficult and less profitable. Farmland is often abandoned, but not perpetually. When farmers are able to cultivate the lands or have no other source of income (e.g., go out to work in cities), they will cultivate them again. The transformation from irrigated into rainfed drylands occurs mainly in the east and southwest, such as in Yulin (11,748.95 km 2 ) and Xinzhou (5032.12 km 2 ). Some irrigated drylands can no longer be effectively irrigated due to the disrepair of the engineering facilities in the old irrigated areas, coupled with groundwater recession and increasing irrigation costs.

Classification Method of Irrigated and Rainfed Drylands
This paper holds that NDWI and NDVI can reflect vegetation growth and vegetation moisture content intensively, and it is unnecessary to use other visible bands of Landsat images. Many of other studies have selected dozens of classification factors, including but not limited to various vegetation indices, such as EVI SVI, etc., [15] and visible bands. Samasse et al. [19] used six bands of surface reflectance and four vegetation indices. The inclusion of environmental factors can reduce, to a certain extent, the index differences caused by regional differences, so that regions with complex natural conditions can be classified using the same method. We only included the sum of precipitation and temperature of the 30 days before imaging. In other studies, precipitation, temperature and evapotranspiration were generally selected as environmental factors [15], as well as wind speed, atmospheric pressure, albedo, etc. For example, Nagaraj et al. selected 11 characteristic bands [16]. In our study, the most classic vegetation indices and the most basic environmental factors were chosen. The data were available on GEE. Different combinations of factors should be studied in the future. In terms of the relationship between remote sensing factors and environmental factors, only positive and negative correlation was reflected. In the future, more accurate numerical relationships should be studied.
In order to obtain the difference of time series curves between irrigated and rainfed drylands, series images with high temporal resolution (every half month during the growing season) were used. Reducing the number of features by selecting the data of critical phases not only reduced the difficulty of classification implementation, but also reduced information interference. For example, from May to June, there was no significant difference between these factors for irrigated drylands and rainfed drylands. From August to September, the characteristics of irrigated and rainfed drylands tended to be the same again. If they were all involved in the operation, the classification accuracy may be affected. Temporal resolution and critical phases should be singled out according to the research object. Therefore, it was confirmed that the characteristic period of growth in arid and semi-arid regions in northwest China is from June to August. In other studies, some used annual averages, and some used averages of dry and wet seasons. Compared to them, this paper has a higher temporal resolution and is more advantageous for capturing irrigated drylands with less distinct characteristics.

Classification Results of Irrigated and Rainfed Drylands
In this paper, the irrigated drylands in drylands were distinguished. In contrast to the extraction of irrigated croplands in other studies, there was no difference in crop types and no significant water characteristics between irrigated and rainfed drylands. In addition, this area was fragmented with small fields and more heterogeneous crops due to the topography and household management patterns, making it more difficult to be extracted than paddy fields or irrigated drylands in the plain region. The validation accuracy of the confusion matrix exceeded 85%, a higher level in this region compared with existing data.
Irrigated drylands in this region mainly included the following categories: cultivated lands located on the edge of the desert, those located on small flood plains, and around rivers. Classification effects of these typical irrigated drylands are shown in Figures 9-11. Irrigated drylands located at the southeastern edge of Mu Us Sandy Land in the northwest of Yuyang District in Figure 9 show that the drylands equipped with circulating sprinkler irrigation can be obviously extracted, based on the data with 30 m spatial resolution. Its shapes were well extracted with reference to the high-resolution image from Google in 2015. Figure 10 shows the twisted Yellow River beach in Hequ County, Shanxi Province. Compared with Google's 2015 HD image, the croplands on the river beach were completely classified as irrigated drylands, while those on the ridges, tablelands, and hills farther away from the Yellow River were classified as rainfed drylands. Figure 11 shows the Yangwu River alluvial fan belonging to the Hutuo River basin, which is the Yangwu river irrigation area in the Yuanping County of Xinzhou City. It is a large irrigated dryland area. The shape of the alluvial fan of the irrigated dryland could be completely extracted. Irrigated drylands in this region mainly included the following categories: cultivated lands located on the edge of the desert, those located on small flood plains, and around rivers. Classification effects of these typical irrigated drylands are shown in Figures 9-11. Irrigated drylands located at the southeastern edge of Mu Us Sandy Land in the northwest of Yuyang District in Figure 9 show that the drylands equipped with circulating sprinkler irrigation can be obviously extracted, based on the data with 30 m spatial resolution. Its shapes were well extracted with reference to the high-resolution image from Google in 2015. Figure 10 shows the twisted Yellow River beach in Hequ County, Shanxi Province. Compared with Google's 2015 HD image, the croplands on the river beach were completely classified as irrigated drylands, while those on the ridges, tablelands, and hills farther away from the Yellow River were classified as rainfed drylands. Figure 11 shows the Yangwu River alluvial fan belonging to the Hutuo River basin, which is the Yangwu river irrigation area in the Yuanping County of Xinzhou City. It is a large irrigated dryland area. The shape of the alluvial fan of the irrigated dryland could be completely extracted.   Irrigated drylands in this region mainly included the following categories: cultivated lands located on the edge of the desert, those located on small flood plains, and around rivers. Classification effects of these typical irrigated drylands are shown in Figures 9-11. Irrigated drylands located at the southeastern edge of Mu Us Sandy Land in the northwest of Yuyang District in Figure 9 show that the drylands equipped with circulating sprinkler irrigation can be obviously extracted, based on the data with 30 m spatial resolution. Its shapes were well extracted with reference to the high-resolution image from Google in 2015. Figure 10 shows the twisted Yellow River beach in Hequ County, Shanxi Province. Compared with Google's 2015 HD image, the croplands on the river beach were completely classified as irrigated drylands, while those on the ridges, tablelands, and hills farther away from the Yellow River were classified as rainfed drylands. Figure 11 shows the Yangwu River alluvial fan belonging to the Hutuo River basin, which is the Yangwu river irrigation area in the Yuanping County of Xinzhou City. It is a large irrigated dryland area. The shape of the alluvial fan of the irrigated dryland could be completely extracted.

Classification of Change Types
In order to avoid the accumulation of classification errors, this study adopted the direct classification of the factor gaps between the two periods and considered the respective spatial clustering characteristics at the same time.
In this study, the difference between the mean values of NDVI and NDWI data of two periods were used for the classification of changes. The dynamic classification index was relatively single. The peak value or variance may also have had great classification significance, which could make the classification more accurate. Currently, sample-based supervised classification is still the main method for the classification of land cover and land use, but it is difficult to study long time series due to the difficulty of collecting historical samples. The gap-based detection is mainly used for quantitative analysis, while the transformation classification based on the time dimension information can assist in the status classification and improve the accuracy of a single period, especially in historical classification. In this study, the classification of 2005 was calibrated against the dynamic data after the direct classification. The distribution of 2005, the dynamic data and the distribution of 2015 could be combined to help with the classification. In other words, the classification in 2015 should also consider the dynamic information of the time dimension, not only with the help of the changes from 2005 to 2015, but also with the help of the changes from 2015 to 2020.

Conclusion
In this study, a method for extracting irrigated drylands based on the GEE platform and Landsat images was proposed. By selecting a one-harvest region in the northern Loess Plateau of China as the study area, the method was proven feasible for extracting irrigated drylands in arid and semi-arid regions. The method was easy to implement, because only four pieces of data, the NDVI, NDWI, temperature, and precipitation, were involved. In addition, a method was proposed to directly classify the types of changes in the two periods using the data gaps and then combining the results of each classification. The conclusions are as follows.
1. The classification accuracy can be effectively improved by using the ratio of NDVI and NDWI to the accumulative values of temperature and precipitation of the 30 days before imaging. The selection of data from the critical crop growing phases can ensure accuracy and greatly reduce the workload. In this study, data of the environmental factors were easily available, and the discriminative calculation relationship was obtained by studying the correlation between remote sensing factors and environmental factors. This provides inspiration for the classification of irrigated drylands in arid and semi-arid regions in northwest China and other regions. At the same time, research on the extraction method of irrigated drylands (no paddy fields), rather than irrigated croplands has been strengthened. While the former is uncommon, they are crucial in water-scarce areas where

Classification of Change Types
In order to avoid the accumulation of classification errors, this study adopted the direct classification of the factor gaps between the two periods and considered the respective spatial clustering characteristics at the same time.
In this study, the difference between the mean values of NDVI and NDWI data of two periods were used for the classification of changes. The dynamic classification index was relatively single. The peak value or variance may also have had great classification significance, which could make the classification more accurate. Currently, sample-based supervised classification is still the main method for the classification of land cover and land use, but it is difficult to study long time series due to the difficulty of collecting historical samples. The gap-based detection is mainly used for quantitative analysis, while the transformation classification based on the time dimension information can assist in the status classification and improve the accuracy of a single period, especially in historical classification. In this study, the classification of 2005 was calibrated against the dynamic data after the direct classification. The distribution of 2005, the dynamic data and the distribution of 2015 could be combined to help with the classification. In other words, the classification in 2015 should also consider the dynamic information of the time dimension, not only with the help of the changes from 2005 to 2015, but also with the help of the changes from 2015 to 2020.

Conclusions
In this study, a method for extracting irrigated drylands based on the GEE platform and Landsat images was proposed. By selecting a one-harvest region in the northern Loess Plateau of China as the study area, the method was proven feasible for extracting irrigated drylands in arid and semi-arid regions. The method was easy to implement, because only four pieces of data, the NDVI, NDWI, temperature, and precipitation, were involved. In addition, a method was proposed to directly classify the types of changes in the two periods using the data gaps and then combining the results of each classification. The conclusions are as follows.
1. The classification accuracy can be effectively improved by using the ratio of NDVI and NDWI to the accumulative values of temperature and precipitation of the 30 days before imaging. The selection of data from the critical crop growing phases can ensure accuracy and greatly reduce the workload. In this study, data of the environmental factors were easily available, and the discriminative calculation relationship was obtained by studying the correlation between remote sensing factors and environmental factors. This provides inspiration for the classification of irrigated drylands in arid and semi-arid regions in northwest China and other regions. At the same time, research on the extraction method of irrigated drylands (no paddy fields), rather than irrigated croplands has been strengthened. While the former is uncommon, they are crucial in water-scarce areas where there are few paddy fields. In addition, this also provides references for land cover and land-use classifications in areas with complex natural environments.
2. In this paper, the classifications of the one-harvest region of the northern Loess Plateau of China complemented the current data, improved the accuracy of existing products in terms of spatial resolution, and provided a more detailed distribution of irrigated and rainfed drylands. This is very helpful for water resource management, ecological environment protection and local agricultural policy making. Compared with humid regions, the local ecological environment is more fragile, and the proportion of agriculture is larger in semi-arid areas, so the harmonious development of agriculture and ecological environment is particularly important. The distribution pattern of irrigated and rainfed drylands can be adjusted according to the development and changes in laws without reducing economic benefits.
3. Directly classifying the gaps of remote sensing factors in specific time dimensions and then combining these with spatial clustering information can lead to more accurate types of changes. Meanwhile, it can also avoid error accumulation in the change maps obtained by the overlaying analysis of two-phase distribution. The introduction of the dynamic information of time dimensions into the classification of historical periods can improve classification accuracy, and compensate for the problem of insufficient samples. This study provides a feasible method for detecting the changes in rainfed and irrigated drylands in arid and semi-arid regions. This method can also be used to study temporal and spatial changes in land use and land cover.