Quantifying the Land Use and Land Cover Changes in the Yellow River Basin while Accounting for Data Errors Based on GlobeLand30 Maps

: Land use and land cover (LULC) change inﬂuences many issues such as the climate, ecological environment, and economy. In this study, the LULC transitions in the Yellow River Basin (YRB) were analyzed based on the GlobeLand30 land use data in 2000, 2010, and 2020. The intensity analysis method with hypothetical errors calculation was used, which could explain the deviations from uniform land changes. The strength of the evidence for the deviation was revealed even though the confusion matrixes of the LULC data at each time point for the YRB were unavailable. The results showed that at the interval scale, the land transition rate increased from the ﬁrst to the second period for all of the upper, middle, and lower reaches. The exchange component was larger than the quantity and shift component, and the gross change was 4.1 times larger than the net change. The size of cultivated land decreased during both intervals. The artiﬁcial surfaces gains were active for all three reaches and had strong evidence. A hypothetical error in 93% of the 2000 data and 58% of the 2010 data can explain deviations from uniform transition given woodland gain during 2000–2010 and 2010–2020. Ecological restoration projects such as Grain for Green implemented in 2000 in the upper reaches resulted in the woodland increase.


Introduction
Insight into land use and land cover (LULC) change at a watershed scale is meaningful because of the multi-functions provided by the land. This computation is particularly important due to the influence of human activities on the watershed. Changes to watershed land use by urbanization, ecological restoration projects, and agriculture are increasingly of concern [1]. These developments lead to LULC changes that affect the ecological environment, crop yield, and water resources [2]. Therefore, there is a need to quantify the land use processes because of increased human activities in the watershed of developing countries. Assessing data errors could complement systematic analysis of LULC change [3].
The assessment of the LULC change characteristics relies on the methods applied and the data used. Some studies compared LULC data at two time points by producing a transition matrix, which is the basis for calculating metrics that are used to analyze temporal changes of each LULC type. The rows of the transition matrix show the categories at the start of the year, whereas the columns show the categories at the end of the year. Comparison of the elements in the transition matrix can summarize the net quantity changes in each time interval. However, it is unable to perform a systematic analysis of land change processes, such as identifying the processes that are more intensive than random and uniform processes. The intensity analysis approach was developed based on Land 2021, 10, 31 3 of 18 protection and high-quality development in the YRB as one of the national strategies in 2019 [36]. Systematic insight into LULC change process is important to supply scientific information for sustainable management of land. However, studies measuring land change are still lacking. This study was designed to address the following questions: (1) What are the size, intensity, and the change components of the LULC transitions for the upper, middle, and lower reaches of the YRB? (2) How could errors in the data affect the LULC change intensities of the YRB from 2000 to 2020? The assessment of LULC change in the YRB would be helpful in optimizing the ecological management strategies in the future.

Study Area
The Yellow River is the fifth largest river worldwide [37]. The Yellow River Basin is the second largest basin in China, supporting 107 million people directly [38]. The YRB is located in northern China (95 • 53 -119 • 5 E, 32 • 10 -41 • 50 N) (Figure 1a), which is also an important ecological barrier in China [39]. It spans nine provinces, and traverses the Tibetan Plateau, the Loess Plateau, Inner Mongolia Plateau, and the North China Plain [40]. The annual mean precipitation ranges from 140 mm in the north to 1100 mm in the east. The annual mean temperature varies from 1 • C in the west to 16 • C in the east [41]. The YRB is located in the typical East Asian monsoon climate zone, and the climate changes from arid and semiarid in the northwest to sub-humid and humid in the southeast. The GDP of the YRB was 23.9 trillion yuan in 2018, which was approximately 26% of the national GDP. The population urbanization rate changed from 38% in 2006 to 54% in 2018 [42]. Over the past several decades, the YRB has experienced the influence of intensive human activities and climate changes that have resulted in serious ecological and environmental issues. results showed that GlobeLand30 and FROM-GLC [28] have the finest spatial resolution at the global extent [29]. Several previous studies on the quality of GlobeLand30 have shown a high accuracy [30][31][32][33][34][35].
In recent years, the YRB has been exposed to intensive human activities and extreme climate events that have led to drastic LULC changes. Several ecological restoration programs have been carried out since 2000. The government of China set ecological protection and high-quality development in the YRB as one of the national strategies in 2019 [36]. Systematic insight into LULC change process is important to supply scientific information for sustainable management of land. However, studies measuring land change are still lacking. This study was designed to address the following questions: (1) What are the size, intensity, and the change components of the LULC transitions for the upper, middle, and lower reaches of the YRB? (2) How could errors in the data affect the LULC change intensities of the YRB from 2000 to 2020? The assessment of LULC change in the YRB would be helpful in optimizing the ecological management strategies in the future.

Study Area
The Yellow River is the fifth largest river worldwide [37]. The Yellow River Basin is the second largest basin in China, supporting 107 million people directly [38]. The YRB is located in northern China (95°53′-119°5′ E, 32°10′-41°50′ N) (Figure 1a), which is also an important ecological barrier in China [39]. It spans nine provinces, and traverses the Tibetan Plateau, the Loess Plateau, Inner Mongolia Plateau, and the North China Plain [40]. The annual mean precipitation ranges from 140 mm in the north to 1100 mm in the east. The annual mean temperature varies from 1 °C in the west to 16 °C in the east [41]. The YRB is located in the typical East Asian monsoon climate zone, and the climate changes from arid and semiarid in the northwest to sub-humid and humid in the southeast. The GDP of the YRB was 23.9 trillion yuan in 2018, which was approximately 26% of the national GDP. The population urbanization rate changed from 38% in 2006 to 54% in 2018 [42]. Over the past several decades, the YRB has experienced the influence of intensive human activities and climate changes that have resulted in serious ecological and environmental issues.

Land Use and Land Cover Data
The LULC data sets for 2000, 2010, and 2020 at a 30 m resolution from GlobeLand30 were used in this study. The spatial resolution was relatively high for such a large scale. The GlobeLand30 data sets were extracted by the hybrid-pixel and object-based method based on Landsat and HJ-1 satellite imageries [17,43]. In this study, each time point was composed of 11 tiles that were downloaded from the website http://www.globallandcover.com/. The tiles were mosaicked to cover the whole study area for each time point. The GlobeLand30 comprises 10 land cover types: cultivated land, forest, grassland, shrub land, wetland, Land 2021, 10, 31 4 of 18 water bodies, tundra, artificial surfaces, bare land, and permanent snow and ice [44]. In this study, these 10 land cover types were reclassified into six aggregate types, including cultivated land, grassland, woodland, waters, artificial surfaces, and unused land ( Figure 1). The reclassification scheme meant that the forest and shrub land were aggregated into woodland, the water bodies and wetlands were aggregated into waters, and the bare land and permanent snow and ice were aggregated into unused land.

Intensity Analysis
The intensity analysis framework was used in this study. The language R was applied to calculate and plot the data. Detailed information can be found in studies by Aldwaik and Pontius [4,20], Pontius and Santacruz [18], and Pontius [45]. The flow chart ( Figure 2) and the main equations are presented as follows.

Land Use and Land Cover Data
The LULC data sets for 2000, 2010, and 2020 at a 30 m resolution from GlobeLand30 were used in this study. The spatial resolution was relatively high for such a large scale. The GlobeLand30 data sets were extracted by the hybrid-pixel and object-based method based on Landsat and HJ-1 satellite imageries [17,43]. In this study, each time point was composed of 11 tiles that were downloaded from the website http://www.globallandcover.com/. The tiles were mosaicked to cover the whole study area for each time point. The GlobeLand30 comprises 10 land cover types: cultivated land, forest, grassland, shrub land, wetland, water bodies, tundra, artificial surfaces, bare land, and permanent snow and ice [44]. In this study, these 10 land cover types were reclassified into six aggregate types, including cultivated land, grassland, woodland, waters, artificial surfaces, and unused land ( Figure 1). The reclassification scheme meant that the forest and shrub land were aggregated into woodland, the water bodies and wetlands were aggregated into waters, and the bare land and permanent snow and ice were aggregated into unused land.

Intensity Analysis
The intensity analysis framework was used in this study. The language R was applied to calculate and plot the data. Detailed information can be found in studies by Aldwaik and Pontius [4,20], Pontius and Santacruz [18], and Pontius [45]. The flow chart ( Figure 2) and the main equations are presented as follows. At the time interval scale, the observed change intensity during year t to year t + 1 was calculated by Equation (1). The uniform change intensity during all time intervals was calculated by Equation (2) where Yt and Yt+1 are years at time point t and t + 1. St is the observed change, and U is the uniform change. If St was greater than U, then the commission error was computed by Equation (3). If St was lower than U, then the omission error was computed by Equation (4). At the time interval scale, the observed change intensity during year t to year t + 1 was calculated by Equation (1). The uniform change intensity during all time intervals was calculated by Equation (2).
where Y t and Y t+1 are years at time point t and t + 1. S t is the observed change, and U is the uniform change. If S t was greater than U, then the commission error was computed by Equation (3). If S t was lower than U, then the omission error was computed by Equation (4).
where E S t is the commission error, and O S t is the omission error. At the categorical level, the annual gain and loss intensity of the categories were computed by Equations (5) and (6), respectively: where G tj and L ti are the annual gain intensity of category j and the annual loss intensity of category i, respectively. If the observed gain of j was larger than the uniform change during [Y t , Y t+1 ], then the hypothetical commission of j error at t + 1 (denoted as E G tj ) was computed using Equation (7). On the contrary, if the observed gain of j was less than the uniform change during [Y t , Y t+1 ], then the hypothetical omission of j error at t + 1 (denoted as O G tj ) was computed by Equation (8).
C tij is the area that transfers from class i to class j during [Y t , Y t+1 ]. If the observed loss of i was larger than the uniform change during [Y t , Y t+1 ], then the hypothetical commission of i error at t (denoted as E L ti ) was computed by Equation (9). On the contrary, if the observed loss of i was less than the uniform change during [Y t , Y t+1 ], then the hypothetical omission of i error at t (denoted as O L ti ) was computed by Equation (10).
At the transition level, when the transitions to category n are examined, Equation (11) gives the observed transition intensity from category j to n (denoted as R tin ). Whereas Equation (12) gives the uniform intensity for class n if n gains uniformly across the study area (denoted as W tn ).
When the transitions from category m are examined, Equation (13) gives the observed transition intensity from category m to j (denoted as Q tmj ). Whereas Equation (14) gives the uniform intensity for class m if m gains uniformly across the study area (denoted as V tm ).
If the transition from i was larger than the uniform gain of n, then the commission of category i error at t (denoted as E R tin ) was calculated by Equation (15). On the contrary, if the transition from i was less than the uniform gain of n, then the omission of category i error at t (denoted as Q R tin ) was calculated by Equation (16).
Similarly, if the transition to j was larger than the uniform loss of n, then the commission error at t + 1 (denoted as E Q tmj ) was calculated by Equation (17). On the contrary, if the Land 2021, 10, 31 6 of 18 transition to j was less than the uniform loss of n, then the omission error at t + 1 (denoted as O Q tmj ) was calculated by Equation (18).
In addition, the overall difference between land use data at two time points comprises three components: quantity, shift, and exchange. The exchange and shift components are allocation differences. Exchange difference is caused by pairwise confusion, whereas shift component is caused by non-pairwise confusion. A clear description on the calculation processes is illustrated in previous studies [18,45].

Interval Level
In this study, we examined the LULC change in the upper, middle, and lower reaches of the YRB during two time intervals. Each interval had the same time duration (2000-2010, and 2010-2020). The change rate was faster from 2010 to 2020 than from 2000 to 2010 for all of the three reaches ( Figure 3a). The change rate was fastest in the lower reaches during 2010 to 2020 with a change percentage of 14.58% of the lower reaches. The LULC change in the upper reaches was the smallest among the three reaches. The commission of change errors were 1.42%, 4.57%, and 1.94% for the domain for the upper, middle, and lower reaches, which can explain the change rate from 2010 to 2020 being larger than the uniform change rate. Instead, the omission of change errors were 1.19%, 2.49%, and 1.44% of the domain for the upper, middle, and lower reaches, which could explain the change during 2000-2010 going slower than the uniform change rate. Figure 3b shows the gross changes during each interval for the three reaches, including the three components: quantity, exchange, and shift. Nearly fifteen percent of the middle and lower reaches showed an overall difference during the second period, and the temporal changes of the upper reaches were relatively less. The exchange was the largest component during both intervals for all reaches except the lower reaches during the second interval. The shift was the smallest component during the two intervals for all three reaches. The quantity components for the upper reaches were much larger than the other two reaches during both intervals.  Figure 4 shows the change intensity of the gains and losses of each category for the upper reaches of the YRB. If a category's bar ends before the uniform intensity line, then the loss or gain intensity of this category was dormant. On the contrary, if a category's intensity is larger than the uniform intensity line, then the change intensity of this category was active. Figure 5 shows the size of the gains and losses of all LULC types along with commission or omission errors at the category level for the upper reaches. The observed change for each active category's gains and losses was the union of the  Figure 4 shows the change intensity of the gains and losses of each category for the upper reaches of the YRB. If a category's bar ends before the uniform intensity line, then the loss or gain intensity of this category was dormant. On the contrary, if a category's intensity is larger than the uniform intensity line, then the change intensity of this category was active. Figure 5 shows the size of the gains and losses of all LULC types along with commission or omission errors at the category level for the upper reaches. The observed change for each active category's gains and losses was the union of the commission error and uniform change. Meanwhile, the union of observed annual gain size and hypothesized omission error size was uniform change size. The grassland experienced larger gain and loss than other categories in terms of size. However, the gain and loss intensities of grassland during both intervals were dormant. The reason was due to the large area of grassland, which accounts for 66.0%, 64.9%, and 62.2% of the total area of the upper reaches in 2000, 2010, and 2020, respectively. The large size of grassland in the denominator produced the small intensity. The gain of artificial surfaces was the most intensively active category during 2010-2020, while the gain of woodland was the most intensively active category during 2000-2010. Furthermore, the error analysis showed that artificial commission error during 2010-2020 was the largest (Figure 5b). That is to say, the evidence for the conclusion of active gain of artificial surfaces was strong. For the upper reaches of the YRB, the change in intensity of woodland was larger than cultivated land during both intervals (Figures 4 and 5).   The categorical level change intensities and sizes for the middle reaches of the YRB are shown in Figure 6 and Figure 7, respectively. The uniform change intensity was 1.46 during 2010-2020, which was 2.65 times higher than for the period 2000-2010 (0.55). During 2000-2010, the intensities of waters' losses and gains were larger than other changes, but the change in the sizes of waters was small. The reason was that the   The categorical level change intensities and sizes for the middle reaches of the YRB are shown in Figure 6 and Figure 7, respectively. The uniform change intensity was 1.46 during 2010-2020, which was 2.65 times higher than for the period 2000-2010 (0.55). During 2000-2010, the intensities of waters' losses and gains were larger than other changes, but the change in the sizes of waters was small. The reason was that the The categorical level change intensities and sizes for the middle reaches of the YRB are shown in Figures 6 and 7, respectively. The uniform change intensity was 1.46 during 2010-2020, which was 2.65 times higher than for the period 2000-2010 (0.55). During 2000-2010, the intensities of waters' losses and gains were larger than other changes, but the change in the sizes of waters was small. The reason was that the proportion of waters in the total area of middle reaches was small (0.53%, 0.51%, and 0.59% for 2000, 2010, and 2020, respectively). It was the same as the upper reaches, where the artificial surfaces were also the most intensively active category during 2010-2020 for the middle reaches. The left side of Figure 7 demonstrates where a class error of 50.4% in the 2010 data can explain the differences between the observed gain intensities and the uniform intensity   The categorical level change intensity and change area on gains and losses are presented in Figure 8 and Figure 9, respectively, for the lower reach of YRB. Four dormant categorical changes during 2000-2010: cultivated land gain, unused land loss, cultivated land loss, and unused land loss are shown in Figure 8a. All other categorical changes were active. During 2010-2020, three dormant categorical changes were indicated: cultivated land gain, cultivated land loss, and artificial surface loss. All other     The categorical level change intensity and change area on gains and losses are presented in Figures 8 and 9, respectively, for the lower reach of YRB. Four dormant categorical changes during 2000-2010: cultivated land gain, unused land loss, cultivated land loss, and unused land loss are shown in Figure 8a. All other categorical changes were active. During 2010-2020, three dormant categorical changes were indicated: cultivated land gain, cultivated land loss, and artificial surface loss. All other changes were shown to be active. Therefore, cultivated land was the only category which had a gain and loss that were dormant during both time intervals in the lower reaches. During 2010-2020, the intensity of artificial surfaces gain was larger than the gains of other categories. The cultivated land loss was the largest in terms of area size. However, the intensity of cultivated land loss was less than the other four categories. The reason was due to the large area of cultivated land, which accounted for 80.6%, 76.8%, and 71.0% of the total area of the lower reaches in 2000, 2010, and 2020, respectively.        For the middle reaches, artificial surfaces change during 2010-2020 was the only category that had a quantity component that was larger than the exchange component. For other LULC types, exchange components were larger than the quantity and shift components during both time intervals (Figure 11). During 2000-2010, cultivated land experienced the largest change due to the large cultivated land gain and even larger cultivated land loss. Whereas during 2010-2020, the change size of grassland was the largest due to the fact of its large gain and even larger loss. For the time interval 2010-2020, the quantity overall line showed that the quantity overall was 11.2% of difference overall, and the quantity + exchange overall was 95.3% of difference overall. Therefore, the exchange component and shift accounted for 84.1% and 4.7% of the difference overall. Figure 11. Size (a) and Intensity (b) of differences for the middle reaches of the YRB. The sign "+" denotes net gain and the sign "−" denotes net loss.

Category Level
For the lower reaches of the YRB, the cultivated land change was largest in terms of size during the two periods, followed by artificial surfaces. The quantity component of cultivated land was more intensive than the quantity overall during the first time interval but became less intensive than the quantity overall during the second time interval (Figure 12). The net changes of all LULC types were stationary through time except grassland and unused land. The quantity changes of grassland and unused land were gains during the first time interval and turned to loss during the latter time interval. For the middle reaches, artificial surfaces change during 2010-2020 was the only category that had a quantity component that was larger than the exchange component. For other LULC types, exchange components were larger than the quantity and shift components during both time intervals (Figure 11). During 2000-2010, cultivated land experienced the largest change due to the large cultivated land gain and even larger cultivated land loss. Whereas during 2010-2020, the change size of grassland was the largest due to the fact of its large gain and even larger loss. For the time interval 2010-2020, the quantity overall line showed that the quantity overall was 11.2% of difference overall, and the quantity + exchange overall was 95.3% of difference overall. Therefore, the exchange component and shift accounted for 84.1% and 4.7% of the difference overall.  For the middle reaches, artificial surfaces change during 2010-2020 was the only category that had a quantity component that was larger than the exchange component. For other LULC types, exchange components were larger than the quantity and shift components during both time intervals (Figure 11). During 2000-2010, cultivated land experienced the largest change due to the large cultivated land gain and even larger cultivated land loss. Whereas during 2010-2020, the change size of grassland was the largest due to the fact of its large gain and even larger loss. For the time interval 2010-2020, the quantity overall line showed that the quantity overall was 11.2% of difference overall, and the quantity + exchange overall was 95.3% of difference overall. Therefore, the exchange component and shift accounted for 84.1% and 4.7% of the difference overall. Figure 11. Size (a) and Intensity (b) of differences for the middle reaches of the YRB. The sign "+" denotes net gain and the sign "−" denotes net loss.
For the lower reaches of the YRB, the cultivated land change was largest in terms of size during the two periods, followed by artificial surfaces. The quantity component of cultivated land was more intensive than the quantity overall during the first time interval but became less intensive than the quantity overall during the second time interval (Figure 12). The net changes of all LULC types were stationary through time except grassland and unused land. The quantity changes of grassland and unused land were gains during the first time interval and turned to loss during the latter time interval. For the lower reaches of the YRB, the cultivated land change was largest in terms of size during the two periods, followed by artificial surfaces. The quantity component of cultivated land was more intensive than the quantity overall during the first time interval but became less intensive than the quantity overall during the second time interval (Figure 12). The net changes of all LULC types were stationary through time except grassland and unused land. The quantity changes of grassland and unused land were gains during the first time interval and turned to loss during the latter time interval.

Transition Level
Considering the rapid urbanization and increased human activities, this study focused on the transition to artificial surfaces and cultivated land from other LULC types. The transition from cultivated land, woodland, grassland, waters, unused land to artificial surfaces in the upper reaches of the YRB is shown in Figure 13. Figure 13a indicates the transition size along with hypothetical errors in the LULC data at the initial year that could account for the difference between transition intensity and the uniform intensity line for the upper reaches. Figure 13b shows the transition intensity. If the bar ends before the uniform intensity line, then the gain of artificial surfaces avoids the corresponding category. On the other hand, if the bar extends beyond the uniform intensity line, then the gain of artificial surfaces targets that category. Gains in artificial surfaces targeted cultivated land and avoided the other four land use types during the two time intervals in the upper and middle reaches (Figures 13 and 14). For the lower reaches, the gains in artificial surfaces targeted both the cultivated land and unused land ( Figure 15). The uniform transition and the errors in the initial years of the time intervals that can account for the departures from uniform transition intensity are shown in Figures 13-18a.

Transition Level
Considering the rapid urbanization and increased human activities, this study focused on the transition to artificial surfaces and cultivated land from other LULC types. The transition from cultivated land, woodland, grassland, waters, unused land to artificial surfaces in the upper reaches of the YRB is shown in Figure 13. Figure 13a indicates the transition size along with hypothetical errors in the LULC data at the initial year that could account for the difference between transition intensity and the uniform intensity line for the upper reaches. Figure 13b shows the transition intensity. If the bar ends before the uniform intensity line, then the gain of artificial surfaces avoids the corresponding category. On the other hand, if the bar extends beyond the uniform intensity line, then the gain of artificial surfaces targets that category. Gains in artificial surfaces targeted cultivated land and avoided the other four land use types during the two time intervals in the upper and middle reaches (Figures 13 and 14). For the lower reaches, the gains in artificial surfaces targeted both the cultivated land and unused land ( Figure 15). The uniform transition and the errors in the initial years of the time intervals that can account for the departures from uniform transition intensity are shown in Figures 13a, 14a, 15a, 16a, 17a and 18a. Land 2021, 10, x FOR PEER REVIEW 12 of 19 Figure 12. Size (a) and Intensity (b) of differences for the lower reaches of the YRB. The sign "+" denotes net gain and the sign "-" denotes net loss.

Transition Level
Considering the rapid urbanization and increased human activities, this study focused on the transition to artificial surfaces and cultivated land from other LULC types. The transition from cultivated land, woodland, grassland, waters, unused land to artificial surfaces in the upper reaches of the YRB is shown in Figure 13. Figure 13a indicates the transition size along with hypothetical errors in the LULC data at the initial year that could account for the difference between transition intensity and the uniform intensity line for the upper reaches. Figure 13b shows the transition intensity. If the bar ends before the uniform intensity line, then the gain of artificial surfaces avoids the corresponding category. On the other hand, if the bar extends beyond the uniform intensity line, then the gain of artificial surfaces targets that category. Gains in artificial surfaces targeted cultivated land and avoided the other four land use types during the two time intervals in the upper and middle reaches (Figures 13 and 14). For the lower reaches, the gains in artificial surfaces targeted both the cultivated land and unused land ( Figure 15). The uniform transition and the errors in the initial years of the time intervals that can account for the departures from uniform transition intensity are shown in Figures 13-18a.                 The transitions from cultivated land to woodland, grassland, and other land use types in the upper and middle reaches of the YRB are shown in Figure 16 and Figure 17, respectively. The loss of cultivated land targeted grassland and artificial surfaces more intensively than other categories in the upper reaches. Whereas in the middle reaches, the annual transition sizes from cultivated land to grassland and woodland were larger than the other categories (Figure 17a). The hypothetical error in the land use data for 2000, which accounts for the divergences from the uniform intensity for the period of 2000-2010, and the hypothetical error in the LULC data for 2010, which accounts for the divergences from the uniform intensity for the period 2010-2020 are shown in Figures  16a and 17a

Discussion
For all of the three reaches of the YRB, the total land change rate of the second period was faster than the first period. The values of the uniform change intensity for 2010-2020 were 1.40, 2.65, and 1.36 times higher than for the period 2000-2010 for the upper, middle, and lower reaches, respectively. For the middle reaches, the uniform change intensity during the second period was almost three times higher than during the first period. There are several reasons for the differences in the uniform change intensity in region and time. First, there were obvious territorial disparities between the upper, middle, and lower YRB. The economic development of the middle reaches lagged behind that of the lower reaches, which led to the slower land transitions in the earlier period and the greater land transitions in the latter period. Second, the land resources of the middle reaches were much richer than the lower reaches, and the contradictions between artificial surfaces, cultivated land, and ecological land were relatively moderate. Accordingly, the restrictions on land use transition were looser than the lower reaches. Third, some ecological restoration and land consolidation projects have been carried out over the last decades. For example, several provinces in the middle reaches such as Shanxi, Inner Mongolia, and Shaanxi have launched the reclamation of abandoned industrial and mining areas from 2012, which contributed to the land transitions in the middle reaches.
In all three reaches, the largest quantity component of change was for artificial surfaces gain in terms of intensity. Whereas in terms of size, the largest quantity components of change were for cultivated land loss in the middle and lower reaches, and for grassland loss in the upper reaches. The temporal land change rate of the upper reaches was slower than the lower and middle reaches. The results are in line with the background of economic development and urbanization in China [46].
The gross changes were 3.6, 6.3, and 2.4 times larger than the net changes in the upper reaches, middle reaches, and lower reaches, respectively. The exchange component was the largest, which is in line with Sbafizadeh-Moghadam et al. [17], who found that the exchange component accounted for three quarters of the overall temporal change during 2000-2010 in Asia based on the GlobeLand30 data in 2000 and 2010. Sohl et al. [47] assessed the modeled LULC and the existing data sources for the conterminous United States, and they found that the exchange component was the largest component of difference between the modeled and the existing LULC data. Cissell et al. [48] mapped and quantified mangrove forest cover changes along Campeche's coast from 1999 to 2016, and the results showed that the exchange component was the largest overall for the land use change in their study. In addition, some studies analyzed the quantity and allocation components, without further dividing the allocation component into exchange and shift components [49].
The programs, such as the large-scale development strategy for the western region, put forward the expansion of artificial surfaces in the YRB. As a result, artificial surfaces were the most intensively active category during 2000-2010 for the three reaches. The increase in the artificial surface led to the occupation of cultivated land and ecological land. Cultivated land continued to be lost for all three reaches during the two intervals, which was caused by rapid industrialization and urbanization over the last two decades. In addition, in the upper and middle reaches, another important reason for the cultivated land loss was the Grain to Green Program and the Slope Cropland Control Projects. In the upper and middle reaches, the ecological environment is fragile and soil erosion is severe, and the cultivation conditions are relatively poor. Therefore, ecological restoration projects such as the Grain to Green Program were implemented from 2000, which aim to return cultivated land to forest and grassland in the ecologically vulnerable areas. As Figure 17 shows, the transition size from cultivated land to grassland and woodland was larger than other categories in the middle reaches, which provides evidence for the effects of the Grain to Green Program. Figure 4 also shows that the woodland area increased most intensively among all of the categories during 2000-2010, which was an achievement of the ecological project.
Although ecological restoration programs have been carried out for almost 20 years, our results still find that the net area of grassland in the upper and middle reaches kept reducing during 2000-2020. This was different from the results from Zhang et al. [19], who found that the area of grassland increased from 2000-2018 markedly. The reason may be that the land use data set used in Zhang et al. and in this study are different. The GlobeLand30 data set with a 30 × 30 m 2 spatial resolution was used in this study, while Zhang et al. used another land use data set with a 1 × 1 km 2 spatial resolution [50]. In addition, some grassland was transferred to woodland with the implementation of the projects of Construction of Windblown Sand Protection Forest [51]. Furthermore, the effectiveness of ecological restoration programs in protecting vegetation land is uncertain in some regions. Jiang et al. [52] found that the soil erosion for the northern dryland of China was substantially reduced; however, the vegetation cover presented decreasing trends significantly in certain sandstorm regions. In these areas, several drought events restrained plant growth and, finally, offset afforestation measures to some extent since 2005. Therefore, identifying suitable time and environmental conditions for implementing ecological restoration programs is critical to achieving targets.
A large amount of high-quality cultivated land is distributed in the lower reaches of the YRB such as in Shandong Province and Hunan Province. The occupation of cultivated land by artificial surfaces would threaten food security. Therefore, the cultivated land Red Line and Requisition Compensation Balance policy was put forward. According to the policies, the losses in cultivated land to artificial surfaces should be replenished with new cultivated land elsewhere, aiming for a constant area of cultivated land in the long term [53]. These cultivated land protection policies were implemented strictly in the provinces of the lower reaches, since these provinces have high-quality cultivated land and are important agricultural supply areas. As Figure 18 shows, the gain in cultivated land targeted unused land mostly in terms of size and intensity. In the lower reaches of the YRB, the unused land was mainly bare land with vegetation cover lower than 10%, including desert, sandy fields, Gobi, bare rocks, saline, and alkaline land, etc. [44]. Therefore, it is important to pay attention to the quality of the new cultivated land that is used to supplement the quantity of cultivated land. The policy aim should be quality equilibrium instead of just quantity equilibrium.
In the upper reaches, the woodland commission hypothetical error was larger than the other errors during 2000-2010, and the artificial surfaces commission hypothesis error was largest during 2010-2020. Therefore, there was powerful evidence for the conclusion of active gain in woodland during 2000-2010 and active gain in artificial surfaces during 2010-2020. During 2000-2010, grassland commission error in the middle reaches was lower; therefore, the evidence for active grassland loss was weak. Similarly, in the middle reaches, the cultivated land omission error was largest, so there was powerful evidence for the conclusion that cultivated land was dormant regarding gains rather than the other conclusion with respect to gains. In the lower reaches, the artificial surfaces commission and cultivated land omission hypothetical errors were larger; therefore, there was powerful evidence that artificial surfaces gain was active and cultivated land gain was dormant.

Conclusions
The application of an intensity analysis framework to GlobeLand30 data of the YRB revealed the size and intensity of LULC changes in this study. The results showed that the overall LULC change in the three reaches of the YRB has been accelerating, reflecting the socioeconomic development of this region during 2000-2020. The increase in the uniform land change intensity was especially large for the middle reaches, owing to the factors of territorial characteristics, land resource enrichment, and environmental projects. The exchange component was the largest overall for each of the reaches during both time intervals, and this result was consistent with several previous studies worldwide. The upper reaches are the environmentally fragile zone, and ecological projects have been carried out since 2000. As a result, the woodland experienced the largest gain intensity during 2000-2010. For the middle reaches, grassland and cultivated land had the largest gains and losses, but cultivated land was dormant during the two intervals owing to the large area it covers. The land resource in the lower reaches is mostly limited, and the contradiction between the artificial surfaces and cultivated land is outstanding. The new cultivated land was mainly transferred from artificial and unused lands under the Requisition-Compensation Balance policy. The size of the minimum hypothetical errors indicated the strength of the evidence for the deviation from the uniform intensity. Based on this research, the causal relationship between LULC change and the potential driving factors for the YRB will be examined quantitatively using statistical methods and machine learning in the future.