Spatial–Temporal Change in Paddy Field and Dryland in Different Topographic Gradients: A Case Study of China during 1990–2020

: As a country with a vast area and complex terrain, the differentiation between paddy ﬁeld and dryland under different topographic gradients in China is difﬁcult. Based on a land-use grid data set with an accuracy of 1 km, this study applied the Topographic Potential Index and used land-use transition matrices and landscape analysis to compare the change in dryland and paddy ﬁeld in China from 1990 to 2020 at different elevations, slopes, and slope aspects. The results indicate that paddy ﬁeld and dryland were mostly distributed in areas with better photothermal conditions. However, in recent years, the paddy ﬁeld and dryland on the “sunny” slope decreased. Speciﬁcally, the area of paddy ﬁeld and dryland on the southeast, south, and southwest slopes decreased, while they increased on the northwest, north, and northeast slopes. From 1990 to 2020, land conversion among paddy ﬁeld, dryland, and other land use was mostly concentrated in the third ladder (<500 m elevation) of China. However, the changes in paddy ﬁeld and dryland have now become active on the second ladder of China. Moreover, the change from other land to dryland on the second ladder accounted for nearly 50% of the country’s change from other land to dryland. Paddy ﬁelds and drylands in areas with low elevation and low slopes were reduced, whereas those with higher elevation and higher slopes increased, indicating that the arable land in mountainous areas increased. This indicates that the topographic conditions of arable land that become worse may aggravate soil erosion in the planting process. The landscape fragmentation of paddy ﬁeld and dryland increased. Compared with paddy ﬁeld, the dryland was more aggregated, the shape was more complex, and the land plots were more fragmented. As a result, paddy ﬁeld and dryland show signiﬁcant differences in their spatial–temporal pattern, landscape characteristics, and land-use changes, and these results can provide an important reference for the sustainable utilization of arable land resources.


Introduction
China is a huge nation with significant regional variations in land and water resources. Divided by the Qinling Mountains-Huaihe River line, a unique pattern of arable land use of "paddy fields in the south and drylands in the north" has formed. In the south, paddy fields account for the majority of the total area of paddy fields in China. The soil types in the south are mainly yellow-brown soil and red soil (on the Middle-Lower Yangtze Plain), lateritic red soil and latosol soil (on the Pearl River Delta Plain), purple soil and paddy soil (in the Sichuan Basin), and yellow soil (in the hilly southeastern region). The climate in the south is dominated by a tropical and subtropical monsoon climate, with an accumulated temperature of 4500-8000 • C, abundant precipitation, and sufficient heat, which forms a pattern of aquatic-based paddy farming. In the north, drylands account for the majority of the total drylands in China. The soil types are mainly Huangmian soil and Heilu soil (in the Loess Plateau), gray-brown desert soil (in Xinjiang), brown soil and cinnamon soil (in the North China Plain), and black soil and dark-brown soil (in the Northeast Plain). The warm-temperate continental monsoon climate in the north, with an accumulated precipitation of approximately between 400 and 800 mm, forms scarce water resources, which contribute to a pattern of wheat-and maize-based dryland farming.
However, the change in paddy field and dryland is continuously dynamic and implicit, and its long-term accumulation may lead to implicit differences in arable land use. Taking the quantity and structure of paddy field and dryland in each province of China as an example, the increase in drylands and the decrease in paddy fields occurred simultaneously over the past 30 years. The province with the largest increase in paddy fields is Heilongjiang, with a cold-temperate climate in the northeast (+27,030 km 2 ). The province with the largest increase in drylands is Xinjiang, with a dry climate in the northwest (+23,306 km 2 ). The province with the largest decrease in both paddy fields and drylands is Guangxi Province, with fertile soil in south China (paddy fields and drylands decreased by 23,306 and 14,719 km 2 , respectively). In addition, relevant studies also pointed out that arable land has decreased in south China and increased in North China [1,2]. Specifically, the paddy fields increased rapidly in the northeast, whereas the drylands increased rapidly in the northwest. The spatial barycenter of paddy field and dryland in China moved northward, which indicates that the arable land moved from better agricultural areas to worse ones [3,4]. Thus, paddy fields and drylands are mostly reduced in areas with a better natural environment, while the newly added drylands and paddy fields are mostly located in areas with worse natural endowments. For example, benefited by the warm and humid climate, the Middle-Lower Yangtze Plain contributed one of highest yields of rice, wheat, and cotton in China. On the other hand, in inland Northwest China, surface water only accounts for approximately 8% of the total runoff in China because of the dry climate. The decrease in arable land in the Yangtze Plain and the increase in arable land in Northwest China shows a spatial mismatch of natural endowment and land resources. As a result, a serious spatial mismatch exists between the distribution of paddy field and dryland and the distribution of high-quality natural resources in China [5][6][7][8][9].
Research on the distribution patterns of paddy field and dryland has increased rapidly over the past 30 years. Previous studies mostly focused on local areas with unique geomorphological characteristics [10,11] and were devoted to analyzing the spatiotemporal change, landscape characteristics, or ecological service efficiencies of paddy field and dryland from the perspective of a single topographic element [10,12,13]. In Heilongjiang Province, Li [14] and Chen [15] found that the spatial expansion of arable land showed a strong directional trend, and the topography and geomorphology conditions were the key factors influencing the changes in paddy field and dryland. Gao [16] found that new arable and lost arable land were mainly concentrated in the plains, followed by the tablelands and hills. It was found that the drylands and paddy fields in the Loess Plateau and Chongqing city were mainly distributed at lower elevations [17,18]. By comparing land-use images from 1933, 1955, 1990, and 2005, Liu [19] found that new paddy fields in the Jinjing River of Hunan Province mainly came from woodland, which was mainly affected by topography conditions (especially elevation). By using the logit model, Zhong [20] analyzed the low mountain hilly area in southeastern China from 1999 to 2006 and found that the loss of agricultural land was the highest at low altitudes, followed by medium and high altitudes. The above studies indicate the importance of topographic factors when analyzing arable land change. Topographic differences impact the configuration of surface water, fertilizer, air, and heat, which forms the spatiotemporal differences between paddy field and dryland. Therefore, the change in paddy field and dryland on different topographic gradients should become a new tool for studying the mechanism of arable land utilization. However, in a vast area with great geographical differences, such as in China, the differentiation of long-term spatial-temporal changes in paddy field and dryland based on topographic gradients is missing.
It is crucial to study the spatial-temporal evolution and internal changes in paddy field and dryland under different topographic gradients in China. As the two major subtypes of arable land, paddy field and dryland may show different changes, further reflecting the implicit transformation of arable land use. This study refines the research perspective Land 2022, 11, 1851 3 of 20 of arable land use by separately observing paddy field and dryland and discusses the following questions: I. When observing topographic conditions in terms of elevation, slope, and slope aspect, are there any differences in the spatial distribution between paddy field and dryland?
II. What is the difference in the land type changes among paddy field, dryland, and other land under different topographic conditions?
III. What do the landscape characteristics of paddy field and dryland look like from the perspective of different topographic conditions?
To address these problems, this paper took the topographic gradient as the research perspective and analyzed the spatial-temporal differentiation of two types of arable land use, paddy field and dryland, at different elevations, slopes, and slope aspects. The scientific basis for optimizing and reconstructing sustainable arable land use was provided based on the research results.

Data Sources
This paper used land-use raster data on China with a resolution of 1 km in 1990, 1995, 2000, 2005, 2010, 2015, and 2020. These data came from the Chinese Academy of Sciences, Resource, and Environmental Science Data Center (https://www.resdc.cn/, accessed on 10 February 2021). The land-use raster data were established based on remote sensing satellite imagery data (Landsat 8 OLI and GF-2). Using a high-resolution remote sensing drone ground survey observation system, the spatial distribution of the land-use cover was interpreted artificially by manual visual interpretation. The testing involved a large number of samples and achieved > 94.3% accuracy [1], and such testing has played an important role in national land research. The land-use classification included six land-use types: arable land, woodland, grassland, water area, construction land, and unutilized land. These data extracted the arable land in the data set and included the two major subtypes of arable land used in the analysis: dryland and paddy field. Specifically, dryland included rain-fed land and land with irrigation facilities that grows dryland crops; paddy field contained land with irrigation facilities that grows aquatic crops.

The Identification of the Topographic Conditions of Paddy field and Dryland
This study used 250 m DEM raster data to calculate the elevations and slopes of paddy field and dryland using the ArcGIS Toolbox. We extracted by mask the 1 km land-use data by the reclassified DEM information to obtain the different topographic conditions of the paddy field and dryland distribution. Based on the suitability of growing staple crops, the elevations were divided into 7 gradients including 0-200, 200-500, 500-1000, 1000-1500, 1500-2500, 2500-3500, and >3500 m. Similarly, the slopes were divided into 5 grades: 0-2 • , 2-6 • , 6-15 • , 15-25 • , and >25 • . The slope directions were divided into north, northeast, east, southeast, south, southwest, west, and northwest categories.

The Transformation between Paddy Field and Dryland
This paper calculated the land-use transformation matrix among paddy field, dryland, and other land-use types based on the toolbox in ArcGIS 10.2. The calculations were as follows: Firstly, paddy field, dryland, and other land-use types in the A th year were reclassified as 1, 2, and 3. Paddy field, dryland, and other land-use types in the B th year were reclassified as 10, 20, and 30.
Secondly, the land-use layers of the A th year and the B th year were summed in the raster calculator, and the attribute data of the newly exported layer indicated the number of rasters that changed from one land-use type in the A th year to another one in the B th year. For instance, 12 indicates the land-use change from a paddy field in the A th year to dryland in the B th year. Detailed attributed data and the implication of the land-use change are listed in Table 1. Thirdly, the total area of land changed from one land-use type to another was calculated by multiplying the number of the changed rasters by the area of the rasters (i.e., 1 km 2 in this study). In land use, different land units combine into different spatial patterns such as land patch densities, shape characteristics, and aggregation degrees. These features constitute the landscape differences between paddy field and dryland. This paper evaluated the landscape characteristic of paddy field and dryland in terms of three aspects: intensity, shape, and vergence. Eight indexes were further analyzed: patch number, landscape aggregation degree, etc. (Table 2). By importing the raster data of paddy field and dryland into Fragstats software, this study calculated the values of the landscape indexes of paddy field and dryland from 1990 to 2020 [21,22]. Table 2. Introduction of the landscape indexes.

Indexes Title Introduction
Intensity NP Number of patches NP is the total number of all patches in the landscape, which can reflect the landscape spatial pattern. PD Patch density PD is the number of patches per unit area.

LPI Largest patch index (%)
LPI is the proportion of the largest patch of a land-use type in the whole landscape, which is used to measure the characteristics of the dominant landscape patches.

LSI
Landscape shape index LSI is the shape complexity of arable land patches, which was measured by the deviation of the shape of the land patch from a circle or a square of the same area.
AWMSI Area-weighted mean shape index AWMSI is the sum of the average shape factor of each land patch multiplied by the weight (the land patch area in the total landscape area). Bigger landscape patches have a higher weight than smaller patches.

AI Aggregation index (%)
AI refers to the aggregation degree of the landscape based on the common boundary length of patches of the same type of landscape.

DIVISION
Landscape isolation DIVISION refers to the individual isolated distribution of different patches in the landscape types.

SPLIT Splitting Index
SPLIT refers to the ratio of landscape fragmentation to the total landscape area index, which is used to describe the dispersion of the landscape pattern.

Topographic Potential Index
Because the distribution of land-use types is highly related to elevation and slope, this paper adopts the Topographic Position Index composed of slope and elevation to reflect the comprehensive topographic features of the paddy field and dryland.
where TPI refers to the Topographic Potential Index; E and E 0 are the elevation and the average elevation of an area, respectively; S and S 0 are the slope and the average slope of an area, respectively. The T value increases along with the increase in slope and elevation. When the slope increases while the elevation decreases, or when the elevation increases while the slope decreases, the value of TPI is in the middle of the two. Referring to the relevant literature [23], this paper classified the topographic position index (abbreviated as TPI) into three levels: high (>0.8), medium (0.4-0.8), and low (<0.4).

Changes in the Topographic Characteristics of Paddy Field and Dryland
The results showed some similarities in the distribution of the TPI of paddy field and dryland ( Figure 1). Firstly, paddy fields and drylands with a low TPI were mainly distributed in the low-elevation areas of eastern China. Specifically, the paddy fields with low a TPI were mainly distributed in the Yangtze Plain areas (including Jiangsu, Shanghai, Anhui, Zhejiang, Hunan, Hubei, and Sichuan), the southern areas of China (Guangdong and Taiwan), the northeastern area (Liaoning and Heilongjiang), and the northeast area (Ningxia Province), while the drylands with a low TPI were mainly distributed in the northern area (Shandong, Henan, Hebei, and Shaanxi provinces), the northeastern area (Liaoning, Jilin, Heilongjiang, and Inner Mongolia provinces), the northwestern area (Xinjiang), and the southern area (Guangdong Province). Secondly, the paddy fields and drylands with a medium TPI were mainly scattered in the second ladder (1000-2000 m elevation) area and small parts of the first ladder (>4000 m elevation) and third ladder (<500 m elevation) in China. Specifically, the paddy fields with a medium TPI were distributed in the south of Shaanxi Province, and the drylands were distributed in Shaanxi, Shanxi, Gansu, Ningxia, Chongqing, Sichuan, Guizhou, Guangxi, Yunnan, and other provinces. Thirdly, the paddy fields and drylands with a high TPI were scattered in the southwest of China, and the drylands area far surpassed that of paddy fields.
An obvious change occurred in the topographic condition of arable land. The paddy fields with a low TPI decreased (Table S1 in the Supplementary Materials), while the paddy fields and drylands with a medium and high TPI increased. Specifically, the proportion of paddy fields with a low TPI decreased from 86.9% to 84.9%, and the proportion of drylands with a low TPI decreased from 76.5% to 75.4%. Because the low TPI areas are mainly distributed in the third ladder of China, which is flat and low in elevation, the hydrothermal conditions and the development degree of its paddy fields and drylands are better, the loss of paddy fields and drylands with a low TPI due to the fact of urban expansion is regrettable, which may lead to a further decline in the suitability of agricultural production [24].
From the perspective of elevation ( Figure 2), 40% of the drylands and 64% of paddy fields were distributed at 0-200 m elevation, which decreased over time. The dryland area was greater than that of paddy field at the same elevation. From 1990 to 2020, the most obvious change in dryland and paddy field was concentrated at an elevation of 0-200 m. When looking at the change in dryland, this study found that the dryland increased in all elevations except for 0-200 m, with the highest annual growth rate (1.01%) in the >3500 m elevation area and the lowest annual growth rate (0.03%) in the 1500-2500 m elevation area. When looking at the change in paddy field, the results showed that the paddy field area increased at 200-500 m, 2500-3500 m, and >3500 m elevations but decreased at the other elevations, with the highest average annual growth rate (2.25%) at >3500 m elevations and the lowest average annual growth rate (0.01%) at 200-500 m elevations.
In terms of the slope condition ( Figure 3), most of the paddy fields and drylands were concentrated at the 0-2 • slope and decreased along with the increase in the slopes. From 1990 to 2020, paddy fields and drylands gradually moved from lower to higher slopes; the area of paddy fields and drylands in flatter areas (0-6 • slope) decreased, whereas the area of paddy fields and drylands in areas with a slope > 6 • increased. When looking at the change trends, dryland showed a fluctuating change (first increased and then decreased at 0-2 • and 2-6 • slopes), while paddy field showed a steadily decreasing trend at 0-2 • and 2-6 • slopes. When looking at the obvious change, dryland changed the most at a 6-15 • slope (increased by 15,026 km 2 ), whereas the paddy field area changed the most at a 0-2 • Land 2022, 11, 1851 6 of 20 slope (decreased by 18,846 km 2 ). When looking at the fastest change, the dryland area and paddy field area increased the fastest (0.6% and 1.09% a year) at a slope of >25 • . In summary, the increase in sloping arable land changed significantly, and the emergence of steeply sloping paddy fields and dryland is especially alarming. Because sloping paddy field and dryland have limits regarding soil and water conservation [17], to guarantee the crop yield sustainably, more slope treatment practices, such as drainage ditches, protection forest, silt arresters, and contour ploughing (terrace field), are needed in this area. From the perspective of elevation (Figure 2), 40% of the drylands and 64% of paddy fields were distributed at 0-200 m elevation, which decreased over time. The dryland area was greater than that of paddy field at the same elevation. From 1990 to 2020, the most obvious change in dryland and paddy field was concentrated at an elevation of 0-200 m. When looking at the change in dryland, this study found that the dryland increased in all  area increased at 200-500 m, 2500-3500 m, and >3500 m elevations but decreased at the other elevations, with the highest average annual growth rate (2.25%) at >3500 m elevations and the lowest average annual growth rate (0.01%) at 200-500 m elevations.  Table S2 in the Supplementary Materials).
In terms of the slope condition ( Figure 3), most of the paddy fields and drylands were concentrated at the 0-2° slope and decreased along with the increase in the slopes. From 1990 to 2020, paddy fields and drylands gradually moved from lower to higher slopes; the area of paddy fields and drylands in flatter areas (0-6° slope) decreased, whereas the area of paddy fields and drylands in areas with a slope >6° increased. When looking at the change trends, dryland showed a fluctuating change (first increased and then decreased at 0-2° and 2-6° slopes), while paddy field showed a steadily decreasing trend at 0-2° and 2-6° slopes. When looking at the obvious change, dryland changed the most at a 6-15° slope (increased by 15,026 km 2 ), whereas the paddy field area changed the most at a 0-2° slope (decreased by 18,846 km 2 ). When looking at the fastest change, the dryland area and paddy field area increased the fastest (0.6% and 1.09% a year) at a slope of >25°. In summary, the increase in sloping arable land changed significantly, and the emergence of steeply sloping paddy fields and dryland is especially alarming. Because sloping paddy field and dryland have limits regarding soil and water conservation [17], to guarantee the crop yield sustainably, more slope treatment practices, such as drainage ditches, protection forest, silt arresters, and contour ploughing (terrace field), are needed in this area.    Table S2 in the Supplementary Materials).
In terms of the slope condition ( Figure 3), most of the paddy fields and drylands were concentrated at the 0-2° slope and decreased along with the increase in the slopes. From 1990 to 2020, paddy fields and drylands gradually moved from lower to higher slopes; the area of paddy fields and drylands in flatter areas (0-6° slope) decreased, whereas the area of paddy fields and drylands in areas with a slope >6° increased. When looking at the change trends, dryland showed a fluctuating change (first increased and then decreased at 0-2° and 2-6° slopes), while paddy field showed a steadily decreasing trend at 0-2° and 2-6° slopes. When looking at the obvious change, dryland changed the most at a 6-15° slope (increased by 15,026 km 2 ), whereas the paddy field area changed the most at a 0-2° slope (decreased by 18,846 km 2 ). When looking at the fastest change, the dryland area and paddy field area increased the fastest (0.6% and 1.09% a year) at a slope of >25°. In summary, the increase in sloping arable land changed significantly, and the emergence of steeply sloping paddy fields and dryland is especially alarming. Because sloping paddy field and dryland have limits regarding soil and water conservation [17], to guarantee the crop yield sustainably, more slope treatment practices, such as drainage ditches, protection forest, silt arresters, and contour ploughing (terrace field), are needed in this area.  When looking at the slope aspects (Figure 4), the results showed that the paddy fields and drylands were mainly distributed on the southern, southwestern, southeastern, eastern, and western slopes of hills or mountains. Moreover, important to note is that the area of paddy fields and drylands on southeastern, southern, and southwestern slopes generally decreased, while the area of paddy fields and drylands on northwestern, northern, and northeastern slopes generally increased. Moreover, the area of paddy fields and drylands increased the most in the north, whereas the area of paddy fields decreased most in the south, and drylands decreased most in the southeast. The results indicate that the area of paddy fields and drylands with better light conditions decreased. Under the influence of man-made disturbance in the urbanization process, the advantages of ideal light and heat resources on arable land were ignored, and the area of paddy fields and drylands with superior natural conditions were reduced. and northeastern slopes generally increased. Moreover, the area of paddy fields and drylands increased the most in the north, whereas the area of paddy fields decreased most in the south, and drylands decreased most in the southeast. The results indicate that the area of paddy fields and drylands with better light conditions decreased. Under the influence of man-made disturbance in the urbanization process, the advantages of ideal light and heat resources on arable land were ignored, and the area of paddy fields and drylands with superior natural conditions were reduced.  Table  S4 in the Supplementary Materials).

The Land Conversion of Paddy Field, Dryland, and Other Land
The spatial distribution of the land-use changes among paddy field, dryland, and other land from 1990 to 2020 is shown in Figure 5.
In terms of the loss of paddy fields, the transition from paddy field to dryland was concentrated in Sichuan, Hebei, and Liaoning provinces. The seasonal water shortage and the decline in groundwater levels may be the major reasons for this change. On the North China Plain, long-term planting and industrial production induces sever groundwater shortages, which form a huge groundwater "funnel" in China. To reduce massive water usage, the government is actively encouraging farmers to change paddy field to dryland by providing subsidies [25]. In Sichuan, seasonal drought occurs frequently due to the spatial-temporal variability of rainfall. Eighty percent of the rainfall in Sichuan is concentrated in May to September. In spring and canicular days, the possibility of droughts is higher than 60% [26]. The transition from paddy field to other land was mostly concentrated in areas with higher urbanization levels, such as the Yangtze River Delta, the Pearl River Delta, and Chengdu-Chongqing. Taking the Yangtze River Delta as an example, fertile land and a temperate climate are beneficial for agricultural production. As one of the richest regions, the Yangtze River Delta contributes more than 20% GDP yearly to the whole country. However, most construction land is transformed from paddy field in this area [27]. Thus, a spatial overlap can be seen between the high-quality paddy fields and economically developed areas [5,9]. As a result, massive losses of paddy fields with better hydrothermal conditions should be protected against with stronger practices, especially in rapidly urbanized areas.  Table S4 in the Supplementary Materials).

The Land Conversion of Paddy Field, Dryland, and Other Land
The spatial distribution of the land-use changes among paddy field, dryland, and other land from 1990 to 2020 is shown in Figure 5.

Land-Type Change in Paddy Field and Dryland at Different Elevations
Land conversion between paddy field, dryland, and other land from 1990 to 2020 was mainly concentrated at an elevation of 0-200 m ( Figure 6). Among them, the greatest area of land conversion occurred between dryland and other land, followed by the conversion between paddy field and other land. The area of the conversion between paddy field and dryland was the lowest.
It is worth noting that when looking at the interconversion between paddy field and dryland, the results show that at a 0-200 m elevation, the area changed from dryland to paddy field was much greater than that of paddy field to dryland. On the contrary, at a >200 elevation, the change from paddy field to dryland was more than that of dryland to paddy field. In terms of the loss of paddy fields, the transition from paddy field to dryland was concentrated in Sichuan, Hebei, and Liaoning provinces. The seasonal water shortage and the decline in groundwater levels may be the major reasons for this change. On the North China Plain, long-term planting and industrial production induces sever groundwater shortages, which form a huge groundwater "funnel" in China. To reduce massive water usage, the government is actively encouraging farmers to change paddy field to dryland by providing subsidies [25]. In Sichuan, seasonal drought occurs frequently due to the spatialtemporal variability of rainfall. Eighty percent of the rainfall in Sichuan is concentrated in May to September. In spring and canicular days, the possibility of droughts is higher than 60% [26]. The transition from paddy field to other land was mostly concentrated in areas with higher urbanization levels, such as the Yangtze River Delta, the Pearl River Delta, and Chengdu-Chongqing. Taking the Yangtze River Delta as an example, fertile land and a temperate climate are beneficial for agricultural production. As one of the richest regions, the Yangtze River Delta contributes more than 20% GDP yearly to the whole country. However, most construction land is transformed from paddy field in this area [27]. Thus, a spatial overlap can be seen between the high-quality paddy fields and economically developed areas [5,9]. As a result, massive losses of paddy fields with better hydrothermal conditions should be protected against with stronger practices, especially in rapidly urbanized areas.
In terms of newly added paddy fields, most of the new paddy fields in the north came from the old drylands, especially in central Henan Province, central Jiangsu Province, and eastern Heilongjiang Province (Table S5 in the Supplementary Materials). In Henan Province, as the hub of the ecological protection and high-quality development strategy of the Yellow River Basin and the Mid-Line of the South-North Water Diversion Project [28], the local paddy field had a good geographical advantage, which helped supply water resources and plain farming and provided the basis for the transition from dryland to paddy fields. In Jiangsu Province, local land consolidation planning policies encouraged the change from dryland plantations to a rice-wheat rotation to increase the economic value of arable land output. This kind of policy support provided a strong incentive for transitioning from local dryland to paddy field. In eastern Heilongjiang Province, the long-term conservation of massive amounts of fertile black soil on the plain area and the flood control projects for rice fields provided sufficient advantages for the transition from dryland to paddy field. Moreover, the consolidation projects of low-yield fields helped to improve the centralization of paddy field patches. In general, paddy field has higher economic benefits (for instance, rice has a higher price compared to wheat and corn), which stimulates farmers to change their dryland to paddy field to obtain a higher income. Because paddy field is a non-negligible source of carbon emissions [29], the transition from dryland to paddy field may lead to an increase in greenhouse gas emissions, which needs further monitoring in future research.
In terms of drylands loss, the transition from dryland to paddy field was concentrated in the central part of Henan Province, the central part of Jiangsu Province, and the eastern part of Heilongjiang Province. The transition from dryland to other land was mainly concentrated in the periphery of urban agglomerations such as the Beijing-Tianjin-Hebei Region. It can be seen that the larger-scale reduced dryland was concentrated in zones with a better natural resource endowment and a stronger economic location.
In terms of the newly added drylands, the transition from paddy field to dryland was concentrated in Sichuan, Hebei, and Liaoning, and a larger scale of converting other land to dryland was more concentrated in Xinjiang and Guangxi. Among them, converting other land to dryland in Guangxi was mainly concentrated in the west of Guangxi, and the primary source of new, local dryland was the reclamation of abandoned industrial and mining land. The new dryland came from areas with poor suitability for development and utilization. Land conversion between paddy field, dryland, and other land from 1990 to 2020 was mainly concentrated at an elevation of 0-200 m ( Figure 6). Among them, the greatest area of land conversion occurred between dryland and other land, followed by the conversion between paddy field and other land. The area of the conversion between paddy field and dryland was the lowest.

Land-Type Change in Paddy Field and Dryland at Different Slopes
From 1990 to 2020, land conversion among paddy field, dryland, and other land was mainly located on 0-2° slopes, and the area of land-use change decreased along with the increase in slope (Figure 7). Among them, the land conversion area between dryland and other land was the largest, followed by the conversion between paddy field and other land, whereas the conversion between paddy field and dryland was the lowest. In terms of the interconversion between paddy field and dryland, the area changed from dryland to paddy field at a 0-2° slope was more than that of paddy field to dryland, but on other slopes, the area changed from paddy field to dryland was larger. In terms of the conversion between dryland and other land, the dryland area changed to other land on a >25° slope was more than that of other land to dryland, but on other slopes, the area changed from other land to dryland was higher. In terms of the conversion between paddy field and other land, the area changed from paddy field to other land was more than that of It is worth noting that when looking at the interconversion between paddy field and dryland, the results show that at a 0-200 m elevation, the area changed from dryland to paddy field was much greater than that of paddy field to dryland. On the contrary, at a >200 elevation, the change from paddy field to dryland was more than that of dryland to paddy field.
In terms of the conversion between paddy field and other land, the area of paddy field to other land at each elevation was more than that of other land to paddy field. Unlike paddy field, the conversion between dryland and other land was not only active at the 0-200 m elevation but also significant at the 1000-1500 m elevation (the "second ladder" of China). This discovery confirms the importance of further analyzing the change in dryland and other lands at higher elevations in western China.

Land-Type Change in Paddy Field and Dryland at Different Slopes
From 1990 to 2020, land conversion among paddy field, dryland, and other land was mainly located on 0-2 • slopes, and the area of land-use change decreased along with the increase in slope (Figure 7). Among them, the land conversion area between dryland and other land was the largest, followed by the conversion between paddy field and other land, whereas the conversion between paddy field and dryland was the lowest. In terms of the interconversion between paddy field and dryland, the area changed from dryland to paddy field at a 0-2 • slope was more than that of paddy field to dryland, but on other slopes, the area changed from paddy field to dryland was larger. In terms of the conversion between dryland and other land, the dryland area changed to other land on a >25 • slope was more than that of other land to dryland, but on other slopes, the area changed from other land to dryland was higher. In terms of the conversion between paddy field and other land, the area changed from paddy field to other land was more than that of other land to paddy field on slopes, which indicates that the decrease in the paddy field area was significant whether the terrain conditions were steep or flat. other land to paddy field on slopes, which indicates that the decrease in the paddy field area was significant whether the terrain conditions were steep or flat.

Land-Type Change in Paddy Field and Dryland on Different Slope Aspects
On the eastern slope, the conversion from paddy field to dryland, the conversion from dryland to other land, and the conversion from paddy field to other land were the largest. On the western slope, the change from dryland to paddy field, the change from other land to dryland, and the change from other land to paddy field were the greatest.
The conversion between dryland and other land was the most prominent, the conversion between paddy field and other land was the second-most prominent, and the conversion between paddy field and dryland was the least prominent.
It is worth noting that, when comparing the southern slope (with better photothermal conditions) and the northern slope (with worse light conditions), the results show that the loss of paddy field and dryland was greater on the southern slope than on the northern slope, while the increase in paddy field and dryland was greater on the northern slope than on the southern slope. In other words, the paddy fields and drylands with better light conditions in mountainous and hilly areas decreased (Figure 8).

Land-Type Change in Paddy Field and Dryland on Different Slope Aspects
On the eastern slope, the conversion from paddy field to dryland, the conversion from dryland to other land, and the conversion from paddy field to other land were the largest. On the western slope, the change from dryland to paddy field, the change from other land to dryland, and the change from other land to paddy field were the greatest.
The conversion between dryland and other land was the most prominent, the conversion between paddy field and other land was the second-most prominent, and the conversion between paddy field and dryland was the least prominent.
It is worth noting that, when comparing the southern slope (with better photothermal conditions) and the northern slope (with worse light conditions), the results show that the loss of paddy field and dryland was greater on the southern slope than on the northern slope, while the increase in paddy field and dryland was greater on the northern slope than on the southern slope. In other words, the paddy fields and drylands with better light conditions in mountainous and hilly areas decreased (Figure 8).

The Landscape Characteristics of Paddy Field and Dryland under Different Topographic Conditions
The landscape indexes of paddy field from 1990 to 2020 are shown in Table 3. In terms of the intensity indices, the NP and PD increased, which indicates that the landscape of the paddy field became more fragmented. The LPI increased, which indicates that the dominance of the main paddy field patches increased. Specifically (Tables S6-S8 in the Supplementary Materials), the largest NPs appeared at 200-500 m, on a 2-6° slope, and on the western slope aspect; the largest increase in the NP appeared at 200-500 m, on a 15-25° slope, and on the northern slope aspect; the largest decrease in the NP appeared at 0-200 m, on a 2-6° slope, and on the southern slope aspect. The largest PDs were distributed at 2500-3500 m, on a 15-25° slope, and on the northeastern slope aspect. The largest increase in the PD appeared at >3500 m and on the southern slope aspect, and the largest decrease in the PD appeared on a 6-15° slope and the northwestern slope aspect. The largest LPIs were distributed at >3500 m DEM, on a 0-2° slope, and on the southeastern slope aspect; the largest increase in the LPI was distributed at 0-200 m and on the southern slope aspect; the largest decrease in the LPI was distributed at >3500 m, on a 6-15° slope, and on the southeastern slope aspect.
In terms of the shape indices, the AWMSI and LSI experienced fluctuated growth, which indicates that the shape of the paddy patches tended to be complex and irregular. Specifically (Tables S6-S8 in the Supplementary Materials), the largest LSIs were distributed at 0-200 m, on a 2-6° slope, and on the western slope aspect; the largest increase in the LSI appeared at 200-500 m, on a 15-25° slope, and on the northern slope aspect; the largest decrease in the LSI appeared at 0-200 m, on a 2-6° slope, and on the southern slope aspect. The largest AWMSIs were distributed at 0-200 m, on a 0-2° slope, and on the southeastern slope aspect; the largest increase in the AWMSI was distributed at 0-200 m; the largest decrease in the AWMSI was distributed at 200-500 m, on a 0-2° slope, and on the southeastern slope aspect.
In terms of the vergence indices, the AI decreased while the DIVISION and SPLIT increased, which meant that the landscape of paddy field tended to be more dispersed. Specifically (Tables S6-S8 in

The Landscape Characteristics of Paddy Field and Dryland under Different Topographic Conditions
The landscape indexes of paddy field from 1990 to 2020 are shown in Table 3. In terms of the intensity indices, the NP and PD increased, which indicates that the landscape of the paddy field became more fragmented. The LPI increased, which indicates that the dominance of the main paddy field patches increased. Specifically (Tables S6-S8 in the Supplementary Materials), the largest NPs appeared at 200-500 m, on a 2-6 • slope, and on the western slope aspect; the largest increase in the NP appeared at 200-500 m, on a 15-25 • slope, and on the northern slope aspect; the largest decrease in the NP appeared at 0-200 m, on a 2-6 • slope, and on the southern slope aspect. The largest PDs were distributed at 2500-3500 m, on a 15-25 • slope, and on the northeastern slope aspect. The largest increase in the PD appeared at >3500 m and on the southern slope aspect, and the largest decrease in the PD appeared on a 6-15 • slope and the northwestern slope aspect. The largest LPIs were distributed at >3500 m DEM, on a 0-2 • slope, and on the southeastern slope aspect; the largest increase in the LPI was distributed at 0-200 m and on the southern slope aspect; the largest decrease in the LPI was distributed at >3500 m, on a 6-15 • slope, and on the southeastern slope aspect.
In terms of the shape indices, the AWMSI and LSI experienced fluctuated growth, which indicates that the shape of the paddy patches tended to be complex and irregular. Specifically (Tables S6-S8 in  the largest decrease in the AWMSI was distributed at 200-500 m, on a 0-2 • slope, and on the southeastern slope aspect. In terms of the vergence indices, the AI decreased while the DIVISION and SPLIT increased, which meant that the landscape of paddy field tended to be more dispersed. Specifically (Tables S6-S8 in the Supplementary Materials), DIVISION decreased at 0-200 m and 1000-1500 m but increased at other elevations, slopes, and slope aspects. The largest SPLITs were distributed at 200-500 m, on a 2-6 • slope, and on the western slope aspect, whereas the largest increase in the number of SPLIT was distributed at 200-500 m, on a 6-15 • slope, and on the northern slope aspect, and the largest decrease in the number of SPLIT was distributed at 1000-1500 m, on a 2-6 • slope, and on the southern slope aspect. The largest AIs were distributed at 0-200 m, on a 0-2 • slope, and on the western slope aspect, the largest increase in the AI was distributed on a >25 • slope and on the northwestern slope aspect, and the largest decrease in the number of SPLIT was distributed at >3500 m, on a 0-2 • slope, and on the southern slope aspect. After 2015, the AI increased while the DIVISION and SPLIT decreased slightly, which indicates that along with the land consolidation projects, the landscape of paddy field appeared to be more concentrated. The intensity indices of dryland are shown in Table 4. The NP and PD increased while the LPI decreased, which indicates that the landscape fragmentation of dryland became more significant. Specifically (Tables S9-S11 in the Supplementary Materials), the largest NPs were distributed at 500-1000 m, on a 2-6 • slope, and on the eastern slope aspect; the largest increase in the NP appeared at 500-1000 m, on a 6-15 • slope, and on the northern slope aspect; the largest decrease in the NP appeared at 0-200 m, on a 0-2 • slope, and on the southeastern slope aspect. The largest PDs were distributed at >3500 m, on a >25 • slope, and on the western slope aspect; the largest increase in the PD was distributed at 1500-2500 m and on the eastern slope aspect; the largest increase in the PD was distributed at >3500 m, on a 15-25 • slope, and on the northern slope aspect. The largest LPIs were distributed at 0-200 m, on a 0-2 • slope, and on the northern slope aspect; the largest increase in the number of LPIs was distributed at 0-200 m, on a 6-15 • slope, and on the northern slope aspect; the largest decrease in the number of LPIs was distributed at >3500 m, on a 0-2 • slope, and on the southern slope aspect.
In terms of the shape indices, the AWMSI and LSI increased, which indicates that the shape of the dryland tended to be more complex and irregular. Specifically (Tables S9-S11 in the Supplementary Materials), the largest LSIs were distributed at 1000-1500 m, on a 2-6 • slope, and on the eastern slope aspect; the largest increase in the LSI was distributed at 500-1000 m, on a 15-25 • slope, and on the northern slope aspect; the largest decrease in the LSI was distributed at 1000-1500 m, on a 2-6 • slope, and on the southeastern slope aspect. The largest AWMSIs were distributed at 0-200 m, on a 0-2 • slope, and on the northern slope aspect; the largest increase in AWMSI was distributed at 0-200 m, on a 0-2 • slope, and on the northern slope aspect; the largest decrease in the AWMSI was distributed at 1000-1500 m, on a 2-6 • slope, and on the southern slope aspect.
In terms of the vergence indices, the AI decreased while the DIVISION and SPLIT increased, which proves the fragmentation and diversification characteristics of dryland landscape. Specifically (Tables S9-S11 in the Supplementary Materials), DIVISION was larger at other elevations, slopes, and slope aspects; the largest increase in the DIVISION was distributed at >3500 m, on a 0-2 • slope. The largest SPLITs were distributed at 2500-3500 m, on a 15-25 • slope, and on the western slope aspect; the most significant increase in SPLIT was distributed at >3500 m, on a >25 • slope, and on the western slope aspect; the largest decrease in the SPLIT was distributed at 200-500 m, on a 6-15 • slope, and on the northern slope aspect. The largest AIs was distributed at 0-200 m, on a 0-2 • slope, and on the northern slope aspect. In general, the landscape fragmentation of paddy field and dryland increased, but dryland showed a higher aggregate degree and more obvious change towards complex shapes and fragmentation of land plots. Furthermore, most of the fragmentation and complexity in the shape of paddy field and dryland were concentrated in low-elevation areas with flat terrain, which emphasizes the loss of high-quality arable land with superior topographic conditions. In addition, the dryland in high-elevation areas tended to be more dispersed and more complicated, which needs further observation for long-term agricultural production.

Changes in Microterrain Factors Led to a Decrease in Paddy Field and Dryland with Good Photothermal Conditions
As the primary resources of arable land use, water, heat, light, and other natural factors impact the growth of crops. Scholars have found that as a microterrain factor, the slope aspect has a significant effect on the growth and spatial distribution of plants [30]. The slope aspect determines the photothermal conditions of vegetation through the reception of solar radiation and hydrological processes which, in turn, affects crop growth. In countries in the Northern Hemisphere, the northern slope aspect has lower temperatures, lower light intensity, higher relative humidity, and more abundant soil nutrients [31][32][33], while the southern slope aspect has sufficient light, less moisture, and large diurnal temperature differences. Therefore, most studies indicate that plants on the southern slope aspect have a higher photosynthetic level and are more productive [34][35][36]. Similarly, this paper showed that 40.56% of paddy fields and 38.48% of drylands in China were distributed on the southern, southwestern, and southeastern slope aspects, which further proves the higher suitability of paddy field and dryland use on southern slopes.
However, the results of this paper found that the area of paddy field and dryland with better light conditions in China was decreasing. The area of paddy fields and drylands on the southeastern, southern, and southwestern slopes generally decreased from 1990 to 2020 but increased on the northwestern, northern, and northeastern slopes. Furthermore, conversion from paddy field and dryland to other land was more active on the southern, southeastern, eastern, and northeastern slopes. These results indicate the worsened solarthermal conditions of paddy field and dryland, which may cause an explicit potential decline in the productivity of arable land. Though some studies have paid attention to the change in arable land area and spatial distribution, the change in productivity caused by the microterrain changes in arable land is rarely discussed, which needs further attention in future research.

In Addition to the "Third Ladder", the Changes in Paddy Field and Dryland Have Become
Active on the "Second Ladder" of China The terrain of China is high in the west and low in the east. According to the elevation from west to east, the whole country can be roughly divided by three ladders including the first ladder (>4000 m) in the west, the second ladder (1000-2000 m) in the middle, and the third ladder (<500 m) in the east. This study found that the land conversion of paddy field and dryland in China mainly occurred in the eastern plain on the third ladder. Similarly, relative studies have supported this conclusion such as the large-scale conversion from paddy field to dryland and construction land in the Beijing-Tianjin-Hebei Region [37], the conversion from other land to dryland and paddy field in the Northeast China Plain [4,15], and the conversion from paddy field and dryland to construction land in Yangtze Plain areas [27].
It is worth emphasizing that this article has a new finding. In addition to the third ladder, the transition between paddy field, dryland, and other land on the second ladder was also active (such as at an elevation of 1000-1500 m). Coincidentally, Chi [38] and Dong [39] also found that a large area of the land conversion of paddy field, dryland, and other land occurred on the second ladder, such as Inner Mongolia, in recent years. Therefore, the land-use change in paddy field and dryland became active on the second ladder. Although climatic conditions, population density, and economic development on the second ladder are poorer than that on the third ladder, the use of arable land at higher elevations is proven to have increased.
The reasons for this can be found by tracing the regional development of the second ladder. Since the implementation of the China Western Development Strategy in 2000, the Chinese government has increased its support for land consolidation and basic arable land construction in the west and has carried out a series of ecological restoration projects to prevent desertification and improve the ecological environment of arable land. Such policy supports resulted in a significant increase in new agricultural modernization, agricultural capital investment, and per capita arable land area [40]. After the agricultural tax was abolished by the No. 1 Central Document in 2006, "The Guidance on Promoting Sustainable Development of Agriculture and Animal Husbandry in Northwest Arid Regions" was issued, which increased agricultural subsidies in poor areas in West China. These have further motivated the enthusiasm for farm production on the second ladder.
It is worth mentioning that in 2017, the government began to permit trans-provincial "land-ticket transactions" to keep an arable land requisition-compensation balance (i.e., arable land occupied by urban construction in provinces with a land shortage can be replenished by provinces with abundant arable land reserve resources). Benefited by abundant land resources, the provinces on the second ladder developed and replenished a large amount of new arable land in "land-ticket transactions", which has also contributed to the dramatic land conversion between paddy field, dryland, and other land at the second ladder.
However, because of the special arid climate environment of the second ladder, the increase in arable land in this area may further intensify the contradiction of water use [41], which will cause the degradation of the ecosystem and the quality of arable land and lead to an increase in desertification over the long term [42][43][44]. Therefore, faced with the active conversion of paddy field and dryland on the second ladder, the local government needs to focus on pushing forward the delineation of the "three lines" (i.e., urban development boundary, permanent basic arable land, and ecological protection red line) and the "three areas" (i.e., urban space, agricultural space, and ecological space) in territorial space planning. Unsuitable arable land with high costs and ecological risks should be retired and the reclamation of arable land in ecologically fragile areas should be limited to rationalize the land-use structure at a higher elevation. Moreover, large-scale land use in Europe has shown some adverse effects such as the persistence of pesticides and other agricultural inputs. Taking a cue from this, the government in China has new opportunities and challenges [45]. The government should further strengthen the arable land-use intensity and high-standard cropland construction and actively promote the protection of arable land such as by letting land lie fallow and through reasonable crop rotation [46]. To improve the acceptance of environmentally friendly techniques, the government should provide more subsidies for eco-friendly agriculture to farmers.

Increased Paddy Field and Dryland on Slopes Exacerbated the Erosion Risk
From the perspective of land suitability, arable land in low-elevation plains is rich in water and heat resources, which is more suitable for agricultural production. This study found that paddy fields at low elevations (0-200 m) accounted for 64.6% of the total paddy field area, and paddy fields on low slope (0-2 • ) areas accounted for 67.62% of the total paddy field area. Drylands at low elevations (0-200 m) and low slopes (0-2 • ) accounted for 39.78% and 61.46% of the total dryland area, respectively. Compared with dryland, paddy field is more sensitive to topography and have higher requirements for water retention. However, this study found that the paddy field and dryland went "up the hill", i.e., an increasing number of paddy fields and drylands appeared at higher elevations (>200 m) and steeper slopes (>6 • ). Because of the large quick water flow on slopes, the phenomenon of soil and water erosion caused by the increase in paddy field and dryland on slopes needs more attention [47]. Rice cultivation in China and Southeast Asia is an important source of farmers' income. The increasingly unideal topographic conditions of paddy fields may bring about more socioeconomic uncertainties, which need further study on the consequences by researchers.
In addition, in terraced fields on slopes, paddy fields tend to exhibit a higher water and fertilizer retention capacity and lower soil erosion intensity than drylands [48]. For example, Gao [49] and Xiao [50] found that paddy fields were less affected by soil erosion than drylands in karst areas, and Chen [51] pointed out that the paddy fields in Taiwan's mountains were less affected by rainfall erosion than drylands. However, this paper found that the transition from paddy field to dryland within a >2 • slope and >200 m elevation was greater than the transition from dryland to paddy field. Taking the mountainous southwest China (including Chongqing, Sichuan, Guizhou, and Yunnan) as an example, the area of paddy field to dryland within a >2 • slope and >200 m elevation was 24,080 and 25,981 km 2 , respectively, which is much higher than the area from dryland to paddy field (15,022 and 16,561 km 2 ). As a result, considering the massive change from paddy field to dryland, stronger water-soil conservation is highly needed to improve the productivity of mountainous arable land.

The Landscape Fragmentation of Paddy Field and Dryland Emphasize the Importance of Land Consolidation
This study states that paddy field and dryland both tended to be dispersed, fragmented, and complex in shape, but compared with paddy field, dryland was more aggregated and showed a stronger change towards complex and fragmented shapes.
Similarly, previous studies showed that Heilongjiang had increased fragmentation in terms of the paddy field landscape [14,52]; landscape fragmentation and the heterogeneity of dryland and paddy field were intensified in Guangdong Province [21]. Moreover, the paddy field in the Jinjing River of Hunan [19] and the paddy field and dryland in the Dongting Lake area [22] showed increased fragmentation and dispersed distribution. These findings are consistent with the conclusions of this paper. Inversely, the decrease in landscape fragmentation of dryland in the Horqin region [53], the fluctuating fragmentation of paddy field and dryland in the Loess Plateau [54], and the decline in the fragmentation of dryland and paddy field decreased in the Karst area of Chishui City [55] and Zhoushan Island [56]. These conclusions are partly different in this study. As a result, the landscape characteristics of paddy field and dryland have zonal differences, but both show fragmentation, irregular shape, and complex distribution. Therefore, increasing the concentration level and contiguity of land plots and strengthening the infrastructure of water, roads, and forests via targeted land consolidation projects should be the major measures used to alleviate the inefficiency of arable land production caused by the unpredictable landscape characteristics of paddy field and dryland.

Conclusions
Based on 1 km land-use grid data from 1990 to 2020, this study applied a landscape analysis, the land-use change matrix for different topographic grades of China, to determine the divergent patterns and the spatial-temporal changes in paddy field and dryland resources under different terrain gradients. The results show that: First, although paddy field and dryland were mostly distributed in areas with better photothermal conditions, the area of paddy fields and drylands on sunny slopes significantly decreased in recent years. The area of paddy fields and drylands on the southeastern, southern, and southwestern slopes decreased, while the area of paddy fields and drylands on the northwestern, northern, and northeast slopes increased.
Second, from 1990 to 2020, land conversion between paddy field, dryland, and other land was mainly concentrated in the third ladder (0-200 m DEM), but the changes in the second ladder (500-2500 m DEM) became prominent.
Third, paddy field and dryland went "up the hill". Paddy field and dryland at lowelevation, low-sloping areas were reduced, and those at higher-elevation, higher-sloping areas increased, which may trigger potential soil erosion. Meanwhile, at high slopes and high elevations, the transition from paddy field to dryland was more prominent.
Finally, there were differences in the landscape characteristics of paddy field and dryland. Although the fragmentation of patches, dispersion of distribution, and complexity of shapes of paddy field landscapes increased, the new paddy field showed the characteristics of aggregation in some areas (such as the northeast and Sichuan). In contrast, although the fragmentation of patches, dispersion of distribution, and complexity of shapes of dryland landscapes increased, the landscape agglomeration characteristics were higher compared to paddy field, decreasing with time.
Based on the research findings, this paper puts forward some suggestions on the utilization of paddy field and dryland. Firstly, the government should strengthen the protection of arable land on sunny slopes and control the conversion from arable land to construction land on sunny slopes in traditional agricultural production areas. Secondly, the protection of water resources in the utilization of arable land in the second ladder regions (such as Xinjiang and Gansu) should be strengthened, and unsuitable arable land should be retired from agricultural production. Thirdly, the government should protect the arable land on the plain area and limit the extensive use of paddy field and dryland on sloping areas to reduce the aggravation of soil erosion. Lastly, land consolidation should emphasize reducing the landscape fragmentation and dispersed distribution of arable land patches.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/land11101851/s1. Table S1: The area of paddy field and dryland under different TPI in China in 1990 and 2020; Table S2: Change in the elevations of paddy field and dryland; Table S3: Change in the slopes of paddy field and dryland; Table S4: Change in slope aspects of paddy field and dryland; Table S5: Land use transformation matrices for paddy field, dryland, and other land among the 32 provinces; Table S6: Landscape indices of paddy field in different elevations; Table S7: Landscape indices of paddy field in different slopes; Table S8: Landscape indices of paddy field in different slope aspects; Table S9: Landscape indices of dryland in different elevations; Table S10: Landscape indices of dryland in different slopes; Table S11: Landscape indices of dryland in different slope aspects.
Author Contributions: Conceptualization, S.X.; methodology, S.X.; software, S.X.; formal analysis, S.X.; investigation, W.W., Q.S. and Z.Z.; data curation, S.X. and W.W.; writing-original draft preparation, S.X.; writing-review and editing, G.Y.; visualization, S.X.; supervision, G.Y.; project administration, S.X.; funding acquisition, G.Y. All authors have read and agreed to the published version of the manuscript. Data Availability Statement: Data are available in a publicly accessible repository that does not issue DOIs. Publicly available data sets were analyzed in this study. These data can be found here: https://www.resdc.cn/, accessed on 10 February 2021.