Study on Spatiotemporal Evolution and Influencing Factors in Cultivated Land and Construction Land in Yunnan Province in the Past 20 Years Based on Remote Sensing Interpretation
Abstract
:1. Introduction
2. Literature Review
3. Materials and Methods
3.1. Overview of the Research Area
3.2. Remote Sensing Data Acquisition and Interpretation
3.3. Spatial Econometric Model
3.4. Indicator System of Influencing Factors
3.5. Research Steps
4. Results
4.1. Analysis of the Change Characteristics of Cultivated Land and Construction Land in Yunnan Province over the Last 20 Years
4.2. Analysis of Spatiotemporal Evolution Characteristics of CULA and COLA
4.2.1. Spatiotemporal Evolution of CULA and COLA
4.2.2. Temporal and Spatiotemporal Evolution of per Capita CULA and per Capita COLA
4.3. Spatial Autocorrelation Analysis and Model Selection
4.4. Analysis of the Influencing Factors in per Capita CULA
- (1)
- Development status of the primary industry (X1). The coefficient estimation of the development status of the primary industry with SEM as the optimal model is 0.0855, passing the test for 1% significance. The model indicates that when considering the spatial autocorrelation, with other factors remaining unchanged, the proportion of primary industry production value to the total number of rural residents who are registered grows by 1%, and the average per capita cultivated land area increases by 0.0855%. The estimated results of other spatial econometric models such as SARAR, SAR, and SDM are similar, with estimated values of 0.0960, 0.1101, and 0.0573, respectively. The estimated results are all positive and successfully pass the 1%, 1%, and 10% significance level tests, respectively, further indicating that the model is robust and that the development of the primary industry can significantly promote the increase in per capita CULA. Agriculture is a significant part of the development of the main industry, even though other sectors such as animal husbandry and fishing are also involved. The rapid growth of the primary industry means that the degree of agricultural development in a region has significantly improved, and the development of agriculture needs to be supported to a certain extent by reserve land resources such as farmland on a larger scale. With the further development of rural industries such as agriculture, the scale and output of grain cultivation will also be further increased, and the dependence on cultivated land will significantly increase. Therefore, the development of the primary industry can, in a certain sense, affect the planning and layout of reserve land resources such as cultivated land and can promote the growth of per capita CULA.
- (2)
- Per capita income level of rural residents (X8). The calculated coefficient of the per capita income level of rural people using SEM as the best model is −0.1333, passing the 1% significance level test. This means that, taking into account spatial autocorrelation and with other factors unchanged, when the per capita disposable income of rural residents increases by 1%, the per capita CULA decreases by 0.1333% on average. The estimated results of other spatial econometric models such as SARAR, SAR, and SDM are similar, with estimated values of −0.1193, −0.1650, and −0.0999, respectively. The estimated results are all negative, and they have passed the 1%, 1%, and 5% significance level tests, respectively, emphasizing further that the model is robust and that the increase in the per capita income level of rural residents has led to a significant decrease in the per capita CULA. Rural residents usually receive a large proportion of their income from migrant work and production and operating income. Due to the generally low food prices in China and the significantly higher prices of some crops with a high value added, including fruits and vegetables, as well as livestock breeding, the huge price gap has led to cultivated land being converted to cash crops or even ponds for fish cultivation, leading to the “non-grain” problem of cultivated land and further reducing CULA. In addition, as the number of rural residents working has increased and their income level has increased significantly, more and more cultivated land is being abandoned or contracted for other purposes, further reducing the area of per capita cultivated land. It can be seen that the increase in rural residents’ income is also a double-edged sword, which may pose a threat to national food security and the protection of cultivated land.
- (3)
- Per capita grain output (X10). Using SEM as the optimal model, the estimated coefficient of per capita grain yield is 0.0006, passing the 1% significance level test. This means that when considering spatial autocorrelation and with other factors remaining unchanged, each increase in per capita grain yield by 1 kg/person will result in an average increase of 0.06% in per capita CULA. The estimation results of other spatial econometric models, such as SARAR, SAR, and SDM, are similar, with the estimated values of 0.0005, 0.0005, and 0.0004, respectively. The predicted outcomes are all positive and all pass the 1% significance level test, further indicating that the model is robust and that the increase in per capita food production can significantly promote the increase in per capita cultivated land area. Grain production is inextricably linked to CULA. Grain production requires a large amount of cultivated land as a basic resource. Only high-quality and large-scale cultivated land can produce more food crops. With the increase in grain production, the dependence on cultivated land becomes more apparent.
- (4)
- Population urbanization level (X11). The coefficient estimation result for the population urbanization level using SEM as the optimal model is −0.0116, passing the test for 1% significance. This means that when considering spatial autocorrelation and with other factors remaining unchanged, the rate of population urbanization increases by 1%, and the average per capita CULA decreases by 0.0116%. The estimation results of other spatial econometric models, such as SARAR, SAR, and SDM, are similar, with estimated values of −0.0119, −0.0118, and −0.0106, respectively. The estimated results are all negative and pass the test for 1% significance, further indicating that the model is robust and that the increase in population urbanization rate has led to a significant decrease in the per capita CULA. The level of urbanization has multiple dimensions, including economy and population. This study uses the degree of urbanization of the population as determined by the ratio of the urban population to the overall population. As industrialization and urbanization continue to advance, the proportion of people living in cities will continue to rise, while that of people living in rural areas will continue to decline. Growing food on cultivated land requires a robust rural population. The decrease in the percentage of rural residents will potentially lead to a reduction in the area under cultivation, and cultivated land may be increasingly abandoned.
- (5)
- Population density (X12). The coefficient estimation result for population density using SEM as the optimal model is −0.0009, passing the 1% significance level test. This means that when considering spatial autocorrelation and with other factors remaining unchanged, an increase in the population urbanization rate by 1 person/km2 will lead to a decrease in the average per capita CULA of 0.09%. The estimation results of other spatial econometric models, such as SARAR, SAR, and SDM, are similar, with estimated values of −0.0009, −0.0008, and −0.0008, respectively. The predicted outcomes are all negative and pass the 1% significance level test, further indicating that the model is robust and that the increase in population density has led to a significant decrease in the per capita CULA. Based on the premise that the total amount of cultivated land remains unchanged, there is an inverse proportional relationship between population density and per capita CULA. The greater the population, the fewer resources per capita there will be. More importantly, spatial evolution analysis shows that the total area of cultivated land in Yunnan Province has decreased in the past 20 years. Although the decline is not significant, with the increase in population, the per capita CULA has decreased significantly; thus, the per capita CULA will decline as population density rises.
- (6)
- Terrain degree of excellence (X15) and climate degree of excellence (X16). The coefficient estimation results of terrain degree of excellence and climate degree of excellence using SEM as the optimal model are 0.5927 and 1.0071, respectively, and have passed the 5% and 1% significance level tests, respectively. This means that when considering spatial autocorrelation and with other factors remaining unchanged, each 0.01 unit increase in terrain degree of excellence or climate degree of excellence will result in an increase in the average per capita CULA by 0.5927% and 1.0071%, respectively. The above results indicate that terrain and climate are the key factors affecting the per capita CULA. The better the terrain and climate, the larger the per capita CULA. This conclusion is also consistent with previous analysis (Figure 5 and Figure 10). We chose the random effect model in this study because the Hausman data failed the 10% significance level test. However, if these two variables are eliminated, the Hausman statistics are significant. As the terrain and climate degrees of excellence are constant values that do not evolve with time in each county, the above results further illustrate that adding these two variables can largely explain the differential effects from individual differences. Therefore, the terrain and climate degrees of excellence are key factors affecting the per capita CULA. Yunnan has a large mountainous area. However, cultivated land is mostly distributed in areas with gentle terrain and good climatic conditions. Therefore, areas with good terrain and climatic conditions also have relatively large per capita CULA. Although the two indicators of terrain and climate do not change over time, they can largely explain the individual differences in the per capita CULA from different regions.
4.5. Analysis of the Influencing Factors in per Capita COLA
- (1)
- Light brightness in the nighttime (X5). The coefficient estimation result for light brightness in the nighttime using SEM as the optimal model is 0.0601, passing the test for 1% significance. This means that when considering spatial autocorrelation and with other factors remaining the same, every 1% increase in light brightness in the nighttime will result in an average increase of 0.0601% in per capita COLA. The estimated results of other spatial econometric models, such as SARAR, SAR, and SDM, are similar, with estimated values of 0.0608, 0.0622, and 0.0338, respectively. The estimated results are all positive, and they pass the 1%, 1%, and 5% significance level tests, respectively, further indicating that the model is robust and that the growth in light brightness in the nighttime can significantly promote an increase in the per capita COLA. Light brightness in the nighttime has a close relationship with the economic growth of a region. According to a study by Henderson et al. (2012), when the light brightness in the nighttime of a region increases by 1%, its GDP increases by approximately 0.3% [50]. Xu Kangning et al. (2015) found that simple GDP statistical variables may have a certain “virtual height”, while nighttime lighting is more authentic [51]. Due to the concrete nature of the variable of nighttime lighting, it can be used as a good alternative variable to GDP. With the continuous growth of a region’s economy, its urbanization construction will also be accelerated; the land currently designated for construction will be unable to meet the development needs of a region or a city. The increase in construction land is an unavoidable trend in economic growth. According to model estimation results, enhanced nighttime lighting is significantly associated with an increase in the per capita COLA.
- (2)
- Per capita public financial expenditure (X7). Using SEM as the optimal model, the estimated coefficient of per capita public finance expenditure is 0.0874, passing the test for 1% significance. This means that when considering spatiotemporal autocorrelation and with other factors remaining unchanged, when per capita public finance expenditure increases by 1%, the average per capita construction land area increases by 0.0874%. The estimation results of other spatial econometric models, such as SARAR, SAR, and SDM, are similar, with estimated values of 0.0746, 0.0418, and 0.0260, respectively. The estimation results are all positive, and the SARAR and SAR models pass the 5% and 10% significance level tests, respectively, further indicating that the model is robust and that the rise in per capita public finance spending is significantly associated with the rise in per capita COLA. In recent years, although an increasing proportion of public funds are used for poverty alleviation and rural construction [52], it is undeniable that a significant amount of public expenditure is also used for various aspects of construction investment, including that devoted to urban infrastructure. As per capita public financial expenditure increases, its investment in infrastructure construction, urbanization construction, and economic development also gradually increases, which means that more public financial investment will be used for the continuous enlargement of construction land, resulting in a continuous increase in per capita COLA.
- (3)
- Per capita grain output (X10). Using SEM as the optimal model, the estimated coefficient of per capita grain yield is 0.0006, passing the test for 1% significance. This means that when considering spatial autocorrelation and with other factors remaining unchanged, each increase in per capita grain yield of 1 kg/person will result in an average increase of 0.06% in per capita COLA. The estimation results of other spatial econometric models, such as SARAR, SAR, and SDM, are similar, with estimated values of 0.0006, 0.0005, and 0.0006, respectively. The estimated results are all positive and pass the test for 1% significance, further indicating that the model is robust and that the increased food output per person has a substantial positive impact on increased per capita COLA. Previous results indicate that there is a clear correlation between the increase in per capita grain production and the increase in per capita CULA, because the production of grain requires a large amount of cultivated land as a basic resource. However, the terrain and climatic conditions in areas with high per capita grain production are generally better, but this high-quality farmland can also be considered as flat land suitable for urban construction. Therefore, high-quality terrain and climatic conditions may indirectly cause the enlargement of construction land.
- (4)
- Population urbanization level (X11). The coefficient estimation result for the population urbanization level using SEM as the optimal model is −0.0102, passing the test for 1% significance. This means that when considering spatial autocorrelation and with other factors remaining unchanged, the average per capita construction land area decreases by 0.0102% for each 1% increase in the population urbanization rate. The estimated outcomes of other spatial econometric models, such as SARAR, SAR, and SDM, are similar, with estimated values of −0.0106, −0.0119, and −0.0107, respectively. The estimated results are all negative and pass the test for 1% significance, further indicating that the model is robust and that the increase in population urbanization rate has led to a significant decrease in the per capita COLA. The urbanization rate has various manifestations. This study uses the ratio of urban population to the overall population to characterize the urbanization level of the population. The reason why the population urbanization rate is inversely correlated with the per capita COLA is that the existing rural residential land accounts for a significant percentage of construction land. The higher the rural population, the more rural residential land there is, whose utilization efficiency is far lower than that of urban construction land. Although the advancement of urbanization will to some extent lead to the expansion of urban “pie spreading”, such per capita growth is smaller compared to the expansion of rural residential land. Compared to rural residential land expansion, a higher urban population results in the effective utilization of construction land and reduces the per capita COLA, Therefore, there is a significant correlation between an increased population urbanization rate and the reduction in per capita CULA.
- (5)
- Population density (X12). The coefficient estimation result for population density using SEM as the optimal model is −0.0004, passing the test for 1% significance. This means that when considering spatial autocorrelation and with other factors remaining unchanged, each increase in population density by 1 person/km2 will lead to an average decrease of 0.04% in per capita construction land. The estimated results of other spatial econometric models, such as SARAR, SAR, and SDM, are similar, with estimated values of −0.0004, −0.0003, and −0.0005, respectively. The estimated results are all negative and pass the test for 1% significance, further indicating that the model is robust and that the rise in habitational density has led to a significant decrease in the per capita COLA. Previous research has shown that in the last 20 years, the total population of Yunnan has experienced growth; however, the rate of increase in construction land area is far higher than the rate of population growth, resulting in a significant increase in the per capita COLA. Although, on the surface, the growth of population and the increase in the per capita COLA have a synchronous trend, there are many factors that affect the per capita COLA. An important advantage of econometric models is that they scientifically “decompose” these impacts and changes, so as to accurately determine true effects of each factor. Population and per capita COLA do not necessarily have a positive correlation, and model estimates indicate that increased density of the population will reduce per capita COLA. The reason for such a result can be explained from the perspective of population migration. In the last 20 years, more and more rural Chinese people have migrated to cities and central areas with rapid economic development, resulting in a significant increase in the area of urban construction land in these areas. Compared to the inefficient nature of rural residential land, the utilization efficiency of urban construction land is higher. For example, the main urban area of Kunming has experienced rapid population growth, becoming a megacity in the past 20 years, with a significant portion of the population migrating from other prefectures. Still, the per capita expansion rate of urban construction land remains limited. However, with the migration of more rural residents to cities, this inefficiently used rural residential land remains, with some areas even becoming “hollow villages”, resulting in extremely low levels of land use efficiency. As the population moves out, the population density gradually decreases, but the total area of construction land does not shrink considerably or even increase, resulting in a gradual increase in per capita construction land.
5. Discussion
6. Conclusions, Policy Recommendations, and Research Prospects
6.1. Main Conclusions of the Study
- (1)
- Using remote sensing interpretation methods, it was found that over the past 20 years, the total and per capita CULA in Yunnan Province have significantly decreased. The total CULA has decreased from 5.5108 million hectares in 2000 to 5.3956 million hectares in 2020, with a net reduction rate of 2.09%. The per capita CULA decreased from 1351.95 m2/person in 2000 to 1142.69 m2/person in 2020, with a net decrease rate of 15.48%. At the same time, the total and per capita COLA have both significantly increased, with the total amount changing from 0.6672 million hectares in 2000 to 1.2969 million hectares in 2020, with a net increase rate of 94.38%. The per capita COLA increased from 163.68 m2/person in 2000 to 274.66 m2/person in 2020, with a net increase rate of 67.80%.
- (2)
- The results of spatiotemporal evolution analysis indicate that the total and per capita CULA in various counties of Yunnan have had relatively small changes over the past 20 years, but the per capita CULA in the southeast, northeast, northwest, and central urban areas of Yunnan is relatively low. During the same period, the total expansion rate of construction land in each county was relatively high, resulting in an “explosive growth” trend in the per capita COLA. By 2020, the per capita COLA in most counties had exceeded 250 m2/person, with a generally significant increase.
- (3)
- The analysis of the influencing factors in per capita CULA shows that improvement of the development status of the primary industry and an increase in per capita grain production can have a beneficial effect on the increase in per capita CULA. However, improvement of per capita income level of rural residents, an increase in population urbanization rate, and an increase in population density will lead to a significant decrease in per capita CULA. The analysis of the influencing factors in per capita COLA shows that an increase in population urbanization rate and population density is conducive to a decrease in per capita COLA. However, an increase in nighttime lighting brightness, per capita public financial expenditure, and per capita grain production will lead to a significant increase in per capita COLA.
6.2. Suggestions for Policy Measures
- (1)
- Improve land use efficiency, optimize the layout of construction land, transform the development of cities at all levels from “extension” to “connotation”, further control the expansion speed of urban construction land, promote intensive land use, and strictly prohibit blind expansion and rough utilization.
- (2)
- Strictly protect existing cultivated land resources and implement a national strategy to prevent the “non-agriculturalization” and “non-grain growing” conversion of cultivated land.
- (3)
- Develop agricultural incentive measures, moderately develop suitable reserve land resources for cultivation, and strive to improve the quality of cultivated land.
- (4)
- Comprehensively consolidate the existing rural residential land (especially hollow villages), convert excess rural construction land into cultivated land and forest land, further stabilize the scale of cultivated land, moderately increase forest coverage, and improve the ecological environment.
6.3. Limitations of the Study and Prospects for Future Research
Author Contributions
Funding
Institutional Review Board Statement
Informed Consent Statement
Data Availability Statement
Conflicts of Interest
References
- Harris, J.R.; Todaro, M.P. Migration, Unemployment and Development: A Two-sector Analysis. Am. Econ. Rev. 1970, 60, 126–142. [Google Scholar]
- Huang, J.; Rozelle, S. Environmental Stress and Grain Yields in China. Am. J. Agric. Econ. 1995, 76, 853–864. [Google Scholar] [CrossRef]
- Wang, C.j.; Fei, X.M. Study on Industrialization, Urbanization and Mechanism of Arable Land Conversion Based on 42 Countries Samples. Areal Res. Dev. 2013, 32, 117–122. [Google Scholar]
- Yang, Z.S.; Li, Y.Y.; Zou, Z.; He, Y.M. Study on Land Resources Development and Utilization Planning of Yunnan Province in the Great Development of Western China, 1st ed.; Yunnan Science and Technology Press: Kunming, China, 2003. [Google Scholar]
- Yang, Z.S.; Zhao, Q.G.; Xin, L. Yunnan Land Resources, 1st ed.; China Science and Technology Press: Beijing, China, 2014. [Google Scholar]
- Office of the Leading Group of the Third National Land Survey of Yunnan Province; Department of Natural Resources of Yunnan Province; Statistics Bureau of Yunnan Province. Main Data Bulletin of the Third National Land Survey of Yunnan Province. Yunnan Daily. 22 December 2021. Available online: https://www.yn.gov.cn/sjfb/tjgb/202112/t20211221_231929.html (accessed on 5 January 2023).
- Li, X.B.; Sun, L.X. Driving Forces of Arable Land Conversion in China, 1st ed.; Cornell University: New York, NY, USA, 1997. [Google Scholar]
- Liu, J.Y.; Zhang, Z.X.; Xu, X.L.; Kuang, W.H.; Zhou, M.C.; Zhang, S.W.; Li, R.D.; Yan, C.Z.; Yu, D.S.; Wu, S.X.; et al. Spatial Patterns and Driving Forces of Land Use Change in China in the Early 21st Century. Acta Geogr. Sin. 2009, 64, 1411–1420. [Google Scholar] [CrossRef]
- Wang, C.J.; He, X.R.; Fei, X.M. Study on the Effect of Industrialization and Urbanization on Cultivated Land Change: An Empirical Analysis Based on International Perspective. J. Agrotech. Econ. 2012, 31, 49–56. [Google Scholar]
- Tan, M.H.; Li, X.B.; Lu, C.H. The Expansion of Construction Land and Its Occupation of Cultivated Land in Large and Medium-sized Cities of China in the 1990s. Sci. China Ser. D Earth Sci. 2004, 34, 1157–1165. [Google Scholar]
- Zhu, L.F.; Huang, J.K. Urbanization and Cultivated Land Changes in China. Econ. Res. J. 2007, 53, 137–145. [Google Scholar]
- Zhang, X.B.; Mount, T.D.; Boisvert, R.N. Industrialization, Urbanization and Land Use in China, 1st ed.; Cornell University: New York, NY, USA, 2008. [Google Scholar]
- Steven, D. Urban Growth, Rural Land Conversion and the Fiscal Well-Being of Local Municipalities, 1st ed.; Cornell University: New York, NY, USA, 2003. [Google Scholar]
- Wang, J.T.; Zhang, H.Z. The Interactive Relation between Industrialization, Urbanization and Land System: The American Experience. Asia-Pac. Econ. Rev. 2009, 25, 52–56. [Google Scholar]
- Jia, S.F.; Zhang, J.Y. Change and Experience of Cultivated Land Use in Urbanization Process in Japan. China Popul. Resour. Environ. 2003, 13, 31–34. [Google Scholar]
- Hao, S.Y.; Wang, J.T.; Zhang, C.Z. International Study of Industrialization, Urbanization and the Evolution of Rural Land System: A Case in Japan. Shanghai J. Econ. 2007, 27, 40–50. [Google Scholar]
- Sun, Q.; Cai, Y.L. Historical Experiences of Cropland Conservation and Land Management in Japan. Acta Sci. Nat. Univ. Pekin. 2008, 44, 249–256. [Google Scholar]
- Huang, D.Q.; Zheng, W.Y. Experience of Urbanization & the Protection of Cultivated Land from Other Countries and Regions. China Land Sci. 2005, 19, 38–43. [Google Scholar]
- Liu, X.W. Analysis of Land Use in Golden Development Stage in Japan, South Korea and China Taiwan. Land Resour. Inf. 2006, 3, 50–56. [Google Scholar]
- Zhao, X.L.; Zhang, Z.X.; Wang, X.; Zuo, L.J.; Liu, B.; Yi, L.; Xu, J.Y.; Wen, Q.K. Analysis of Chinese Cultivated Land’s Spatial-temporal Changes and Causes in Recent 30 Years. Trans. Chin. Soc. Agric. Eng. 2014, 30, 1–11. [Google Scholar]
- Shi, J.; Zhang, F.R.; Zhao, T.T. Spatiotemporal Changes in Cultivated Land Area in China from 1998 to 2006. Resour. Sci. 2008, 30, 1191–1198. [Google Scholar]
- Li, J.H.; Zhao, J.Z.; Chen, W.B.; Chen, F.G.; Zhao, X.M. Cultivated Land Change and Social Economic Driving Forces of Jiangxi Province over the Past 50 Years. Acta Agric. Univ. Jiangxiensis 2006, 28, 467–471. [Google Scholar]
- Du, X.B.; Zhou, W.; Si, H.J.; Cao, Y.G.; Cai, C.Y. Driving Forces of Cultivated Land Change in Qinghai Province during 2000–2008. Res. Soil Water Conserv. 2013, 20, 180–186. [Google Scholar]
- Zhao, Y.S. Principles and Methods for Analysis of Remote Sensing Applications, 2nd ed.; Science Press: Beijing, China, 2013. [Google Scholar]
- Liu, J.Y.; Zhang, Z.X.; Zhuang, D.F.; Wang, Y.M.; Zhou, W.C.; Zhang, S.W.; Li, R.D.; Jiang, N.; Wu, S.X. A Study on the Spatial-temporal Dynamic Changes of Land-use and Driving Forces Analyses of China in the 1990s. Geogr. Res. 2003, 22, 1–12. [Google Scholar]
- Liu, J.Y.; Kuang, W.H.; Zhang, Z.X.; Xu, X.L.; Qin, Y.W.; Ning, J.; Zhou, W.C.; Zhang, S.W.; Li, R.D.; Yan, C.Z.; et al. Spatiotemporal Characteristics, Patterns, and Causes of Land-use Changes in China since the Late 1980s. J. Geogr. Sci. 2014, 24, 195–210. [Google Scholar] [CrossRef]
- Liu, J.Y.; Ning, J.; Kuang, W.H.; Xu, X.L.; Zhang, S.W.; Yan, C.Z.; Li, R.D.; Wu, S.X.; Hu, Y.F.; Du, G.M.; et al. Spatio-temporal Patterns and Characteristics of Land-use Change in China during 2010–2015. Acta Geogr. Sin. 2018, 73, 789–802. [Google Scholar]
- Fu, B.J.; Chen, L.D.; Cai, Y.L.; Li, X.B. Study on Land Use Change and Sustainable Use in the Areas around Bo Hai, 1st ed.; Science Press: Beijing, China, 2004. [Google Scholar]
- Xiao, H.Y.; Liu, G.D.; Yang, Z.; Fang, N.; Dai, H.M.; Liu, K.; Han, X.M.; Zhu, H. Spatio-temporal Change in Land Uses of the Black Soil Area in Northeast China over the Past Half Century. Geophys. Geochem. Explor. 2022, 46, 1037–1049. [Google Scholar]
- Sayedain, S.A.; Maghsoudi, Y.; Eini-Zinab, S. Assessing the Use of Cross-Orbit Sentinel-1 Images in Land Cover Classification. Int. J. Remote Sens. 2020, 41, 7801–7819. [Google Scholar] [CrossRef]
- Ramanamurthy, B.V.; Victorbabu, N. Land Use Land Cover (LULC) Classification with Wasteland Demarcation Using Remote Sensing and GIS Techniques. Mater. Sci. Eng. 2021, 1025, 12–35. [Google Scholar] [CrossRef]
- Liang, J.; Xu, J.C.; Shen, H.F.; Fang, L. Land-use Classification via Constrained Extreme Learning Classifier Based on Cascaded Deep Convolutional Neural Networks. Eur. J. Remote Sens. 2020, 53, 219–232. [Google Scholar] [CrossRef]
- McCarty, D.A.; Kim, H.W.; Lee, H.K. Evaluation of Light Gradient Boosted Machine Learning Technique in Large Scale Land Use and Land Cover Classification. Environments 2020, 7, 84. [Google Scholar] [CrossRef]
- Rajendran, G.B.; Kumarasamy, U.M.; Zarro, C.; Divakarachari, P.B.; Ullo, S.L. Land-Use and Land-Cover Classification Using a Human Group-Based Particle Swarm Optimization Algorithm with an LSTM Classifier on Hybrid Pre-Processing Remote-Sensing Images. Remote Sens. 2020, 12, 4135. [Google Scholar] [CrossRef]
- Li, X.; Hu, T.; Gong, P.; Du, S.; Chen, B.; Li, X.; Dai, Q. Mapping Essential Urban Land Use Categories in Beijing with a Fast Area of Interest (AOI)-Based Method. Remote Sens. 2021, 13, 477. [Google Scholar] [CrossRef]
- Statistics Bureau of Yunnan Province. Yunnan Statistical Yearbook-2021, 1st ed.; China Statistics Press: Beijing, China, 2021. [Google Scholar]
- Office of Yunnan Provincial Agricultural Zoning Commission. Land Area of Different Climatic Zones and Slopes in Yunnan Province, 1st ed.; Yunnan Science and Technology Press: Kunming, China, 1987. [Google Scholar]
- National Bureau of Statistics of the People’s Republic of China. China Statistical Yearbook-2021, 1st ed.; China Statistics Press: Beijing, China, 2021. [Google Scholar]
- Xu, X.L.; Liu, J.Y.; Zhang, S.W.; Li, R.D.; Yan, C.Z.; Wu, S.X.; Multi-Period Land Use and Land Cover Monitoring Data Set in China (CNLUCC). Beijing: Resource and Environment Science Registration and Publication System, Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences. Available online: http://www.resdc.cn/DOI (accessed on 26 December 2018).
- Liu, J.Y.; Liu, M.L.; Deng, X.Z.; Zhuang, D.F.; Zhang, Z.X.; Luo, D. The land use and land cover change database and its relative studies in China. J. Geograph. Sci. 2002, 12, 275–282. [Google Scholar]
- Liu, J.Y.; Liu, M.L.; Zhuang, D.F.; Zhang, Z.X.; Deng, X.Z. Study on Spatial Pattern of Land-use Change in China during 1995–2000. Sci. China Ser. D 2003, 46, 373–384. [Google Scholar]
- Kuang, W.H.; Zhang, S.W.; Du, G.M.; Yan, C.Z.; Wu, S.X.; Li, R.D.; Lu, D.S.; Pan, T.; Ning, J.; Guo, C.Q.; et al. Remotely Sensed Mapping and Analysis of Spatio-temporal Patterns of Land Use Change Across China in 2015–2020. Acta Geogr. Sin. 2022, 77, 1056–1071. [Google Scholar]
- Yang, R.Y.; Zhong, C.B. Analysis on Spatio-Temporal Evolution and Influencing Factors of Air Quality Index (AQI) in China. Toxics 2022, 10, 712. [Google Scholar] [CrossRef] [PubMed]
- Tobler, W. A Computer Movie Simulating Urban Growth in the Detroit Region. Econ. Geogr. 1970, 46, 234–240. [Google Scholar] [CrossRef]
- Moran, P. Notes on Continuous Stochastic Phenomena. Biometrika 1950, 37, 17–23. [Google Scholar] [CrossRef] [PubMed]
- Yang, Z.S.; Yang, R.Y.; Liu, F.L. Spatio-temporal Evolution and Influencing Factors of Urban-rural Income Gap in Yunnan Province Based on Poverty Classification. Geogr. Res. 2021, 40, 2252–2271. [Google Scholar]
- Elhorst, J.P. Spatial Econometrics: From Cross-Sectional Data to Spatial Panels, 1st ed.; Springer: Berlin/Heidelberg, Germany, 2014. [Google Scholar]
- Yang, R.Y.; Zhong, C.B.; Yang, Z.S.; Liu, F.L.; Peng, H.Y. Analysis on Poverty Influencing Factors in Deep Poverty County of Karst Rocky-desertified Area in Southwest China. World Reg. Stud. 2022, 31, 1298–1309. [Google Scholar]
- Yang, Z.S. Study on Land Use Changes and Its Ecological Effects in Different Landform Areas in Yunnan Province Driven by China’s Project of Converting Farmland to Forest, 1st ed.; China Science and Technology Press: Beijing, China, 2011. [Google Scholar]
- Henderson, V.; Storeygard, A.; Weil, D.N. Measuring Economic Growth from Outer Space. Am. Econ. Rev. 2012, 102, 994–1028. [Google Scholar] [CrossRef]
- Xu, K.N.; Chen, F.L.; Liu, X.Y. The Truth of China Economic Growth: Evidence from Global Night-time Light Data. Econ. Res. J. 2015, 9, 17–29. [Google Scholar]
- Yang, R.Y.; Zhong, C.B.; Yang, Z.S.; Wu, Q.J. Analysis on the Effect of the Targeted Poverty Alleviation Policy on Narrowing the Urban-Rural Income Gap: An Empirical Test Based on 124 Counties in Yunnan Province. Sustainability 2022, 14, 12560. [Google Scholar] [CrossRef]
Year | Remote Sensing Image Data | Spatial Resolution |
---|---|---|
2000 | Landsat-TM/ETM remote sensing images from December 1999 to February 2000 | 30 m × 30 m |
2010 | Landsat-TM remote sensing images from December 2009 to February 2010 | 30 m × 30 m |
2020 | Landsat-8 remote sensing images from January 2020 to February 2020 | 30 m × 30 m |
Dimensions | Variables | Method of Calculation | Name | Unit |
---|---|---|---|---|
Dependent Variables | Per Capita Cultivated Land Area (CULA) | ln (total area of local cultivated land/total population) | Y1 | m2/person |
Per Capita Construction Land Area (COLA) | ln (total area of local construction land/total population) | Y2 | m2/person | |
Industrial Economy | Development Status of Primary Industry | ln (primary industry output value/rural registered residence population) | X1 | CNY/person |
Development Status of Secondary Industry | ln (secondary industry output value/urban registered residence population) | X2 | CNY/person | |
Development Status of Tertiary Industry | ln (tertiary industry output value/total population) | X3 | CNY/person | |
Industrial Structure | output value of secondary and tertiary industries/GDP × 100% | X4 | % | |
Light Brightness in the Nighttime | ln (average brightness of night lights + 0.01) | X5 | None | |
Investment Expenditure | Fixed Assets Investment | ln (total fixed assets investment/total population) | X6 | CNY/person |
Per Capita Public Financial Expenditure | ln (public finance expenditure/total population) | X7 | CNY/person | |
Rural Development | Per Capita Income Level of Rural Residents | ln (per capita disposable income of rural residents) | X8 | CNY/person |
Proportion of Rural Employees | rural employees/total rural population | X9 | % | |
Per Capita Grain Output | total grain output/total population | X10 | kg/person | |
Population Structure | Population Urbanization Level | (total population − agricultural population)/total population × 100% | X11 | % |
Population Density | total population/land area | X12 | person/km2 | |
Ecological Environment | Biological Abundance Index | Abio × (0.35 × woodland area + 0.21 × grassland area + 0.28 × water area + 0.11 × cultivated land area + 0.04 × construction land area + 0.01 ×unused land area)/total land area; where Abio = 511.2642 | X13 | None |
Normalized Difference Vegetation Index (NDVI) | average local NDVI in the current year | X14 | None | |
Geographical Conditions | Terrain Degree of Excellence | 1 × proportion of land area with slope ≤ 8° + 0.75 × proportion of land area with a slope of 8~15° + 0.5 × land area ratio with a slope of 15~25° + 0.25 × land area ratio with a slope of 25~35° + 0.9 × proportion of water area and land area | X15 | None |
Climate Degree of Excellence | 1 × lower thermosphere + 0.8 × medium warm layer + 0.5 × high cold layer + 0.9 × Other | X16 | None |
First-Level Land Use Type | Second-Level Land Use Type | Land Use Classified Areas (Unit: 1 × 104 Hectares; m2/Person) | ||
---|---|---|---|---|
2000 | 2010 | 2020 | ||
Cultivated Land | Paddy Field | 135.91 | 134.53 | 131.39 |
Dry Farmland | 415.17 | 411.43 | 408.17 | |
Total | 551.08 | 545.96 | 539.56 | |
Cultivated Land (Per Capita) | Paddy Field (Per Capita) | 333.42 | 292.35 | 278.26 |
Dry Farmland (Per Capita) | 1018.52 | 894.09 | 864.43 | |
Total (Per Capita) | 1351.95 | 1186.44 | 1142.69 | |
Construction Land | Urban Construction Land, Rural Settlement Area, and Land for Mining and Industry | 54.84 | 74.86 | 109.17 |
Other Building Land | 11.88 | 11.87 | 20.52 | |
Total | 66.72 | 86.73 | 129.69 | |
Construction Land (Per Capita) | Urban Construction Land, Rural Settlement Area, and Land for Mining and Industry (Per Capita) | 134.54 | 162.68 | 231.20 |
Other Building Land (Per Capita) | 29.14 | 25.80 | 43.46 | |
Total (Per Capita) | 163.68 | 188.48 | 274.66 |
Dimension | Using Natural Logarithm | Items | 2000 | 2010 | 2020 | Annual Average |
---|---|---|---|---|---|---|
Per Capita Cultivated Area | No | Moran’s I | 0.447 *** | 0.424 *** | 0.408 *** | 0.446 *** |
Z statistics | 8.587 | 8.147 | 7.845 | 8.575 | ||
Yes | Moran’s I | 0.474 *** | 0.506 *** | 0.551 *** | 0.524 *** | |
Z statistics | 9.226 | 9.941 | 10.819 | 10.262 | ||
Per Capita Construction Land Area | No | Moran’s I | 0.156 *** | 0.269 *** | 0.296 *** | 0.275 *** |
Z statistics | 3.689 | 5.246 | 5.740 | 5.356 | ||
Yes | Moran’s I | 0.329 *** | 0.308 *** | 0.320 *** | 0.331 *** | |
Z statistics | 6.462 | 5.990 | 6.219 | 6.414 |
Dependent Variables | Types | Statistics | All Independent Variables Are Considered | Independent Variables Except Collinear Variables Are Considered | ||
---|---|---|---|---|---|---|
Statistical Values | p Values | Statistical Values | p Values | |||
Per Capita CULA | Spatial Error | Moran’s I | 1.869 * | 0.062 | 1.800 * | 0.072 |
LM | 41.365 *** | 0.000 | 40.092 *** | 0.000 | ||
Robust LM | 39.213 *** | 0.000 | 37.842 *** | 0.000 | ||
Spatial Lag | LM | 3.061 * | 0.080 | 3.346 * | 0.067 | |
Robust LM | 0.909 | 0.341 | 1.096 | 0.295 | ||
Per Capita COLA | Spatial Error | Moran’s I | 1.643 * | 0.100 | 1.734 * | 0.083 |
LM | 27.990 *** | 0.000 | 35.981 *** | 0.000 | ||
Robust LM | 26.763 *** | 0.000 | 34.358 *** | 0.000 | ||
Spatial Lag | LM | 2.164 | 0.141 | 2.866 * | 0.090 | |
Robust LM | 0.937 | 0.333 | 1.242 | 0.265 |
Items | RE | SARAR(FE) | SAR(RE) | SEM(RE) | SDM(RE) |
---|---|---|---|---|---|
Development Status of the Primary Industry (X1) | 0.1236 ** (0.0573) | 0.0960 *** (0.0276) | 0.1101 *** (0.0341) | 0.0855 *** (0.0330) | 0.0573 * (0.0322) |
Development Status of the Secondary Industry (X2) | −0.0169 (0.0167) | 0.0035 (0.0133) | −0.0069 (0.0158) | 0.0035 (0.0157) | 0.0297 * (0.0158) |
Industrial Structure (X4) | 0.0024 (0.0015) | 0.0012 (0.0011) | 0.0019 (0.0013) | 0.0009 (0.0013) | −0.0010 (0.0013) |
Light Brightness in the Nighttime (X5) | −0.0064 (0.0153) | −0.0158 (0.0101) | −0.0065 (0.0122) | −0.0219 * (0.0122) | −0.0222 * (0.0119) |
Per Capita Public Financial Expenditure (X7) | −0.0281 (0.0216) | −0.0333 ** (0.0158) | −0.0177 (0.0182) | −0.0318 * (0.0189) | 0.0594 ** (0.0293) |
Per Capita Income Level of Rural Residents (X8) | −0.1561 *** (0.0434) | −0.1193 *** (0.0318) | −0.1650 *** (0.0343) | −0.1333 *** (0.0376) | −0.0999 ** (0.0423) |
Proportion of Rural Employees (X9) | −0.0011 (0.0012) | −0.0003 (0.0009) | −0.0018 (0.0011) | −0.0004 (0.0011) | −0.0010 (0.0010) |
Per Capita Grain Output (X10) | 0.0006 *** (0.0001) | 0.0005 *** (0.0001) | 0.0005 *** (0.0001) | 0.0006 *** (0.0001) | 0.0004 *** (0.0001) |
Population Urbanization Level (X11) | −0.0117 *** (0.0021) | −0.0119 *** (0.0010) | −0.0118 *** (0.0013) | −0.0116 *** (0.0012) | −0.0106 *** (0.0012) |
Population Density (X12) | −0.0009 *** (0.0001) | −0.0009 *** (0.0000) | −0.0008 *** (0.0001) | −0.0009 *** (0.0001) | −0.0008 *** (0.0001) |
Biological Abundance Index (X13) | 0.0048 ** (0.0022) | 0.0008 (0.0018) | 0.0036 ** (0.0018) | 0.0033 * (0.0018) | 0.0014 (0.0018) |
NDVI (X14) | 0.0447 * (0.0261) | 0.0046 (0.0351) | 0.0588 * (0.0339) | 0.0314 (0.0423) | −0.0076 (0.0560) |
Terrain Degree of Excellence (X15) | 0.5762 (0.3589) | — | 0.6307 ** (0.2708) | 0.5927 ** (0.2902) | 0.6182 (0.4504) |
Climate Degree of Excellence (X16) | 1.0085 ** (0.4662) | — | 0.7914 ** (0.3169) | 1.0071 *** (0.3388) | −0.3780 (0.7493) |
_cons | 6.0737 *** (0.6293) | — | 4.6455 *** (0.5736) | 6.3563 *** (0.5488) | 2.4399 *** (0.7895) |
Parameter ρ | — | −0.0234 (0.0536) | 0.2388 *** (0.0484) | — | 0.5168 *** (0.0588) |
Parameter λ | — | 0.5258 *** (0.0636) | — | 0.5257 *** (0.0709) | — |
LR Test: Individual Effect | — | 33.10 *** (0.0072) | — | — | — |
LR Test: Time Effect | — | −1390.18 (1.0000) | — | — | — |
Hausman Test | 16.35 (0.1759) | — | 7.70 (0.8079) | 16.35 (0.1759) | 9.09 (0.6956) |
Individual Effect | — | Yes | — | — | — |
Time Effect | — | No | — | — | — |
Within R2 | 0.8410 | 0.8379 | 0.8233 | 0.8356 | 0.8630 |
Sample Size | 387 | 387 | 387 | 387 | 387 |
Items | FE | SARAR(FE) | SAR(FE) | SEM(FE) | SDM(FE) |
---|---|---|---|---|---|
Development Status of the Primary Industry (X1) | 0.0838 (0.0572) | 0.0985 ** (0.0474) | 0.1439 *** (0.0447) | 0.0899 * (0.0464) | 0.0624 (0.0436) |
Development Status of the Secondary Industry (X2) | −0.0085 (0.0277) | −0.0285 (0.0213) | −0.0283 (0.0204) | −0.0264 (0.0213) | −0.0166 (0.0211) |
Industrial Structure (X4) | −0.0001 (0.0023) | 0.0000 (0.0018) | 0.0017 (0.0016) | −0.0003 (0.0017) | −0.0001 (0.0017) |
Light Brightness in the Nighttime (X5) | 0.0217 (0.0202) | 0.0608 *** (0.0163) | 0.0622 *** (0.0156) | 0.0601 *** (0.0163) | 0.0338 ** (0.0159) |
Per Capita Public Financial Expenditure (X7) | −0.0081 (0.0469) | 0.0746 ** (0.0294) | 0.0418 * (0.0231) | 0.0874 *** (0.0264) | 0.0260 (0.0397) |
Per Capita Income Level of Rural Residents (X8) | −0.1602 ** (0.0630) | 0.0583 (0.0527) | −0.0059 (0.0460) | 0.0653 (0.0531) | −0.1157 ** (0.0567) |
Proportion of Rural Employees (X9) | −0.0031 * (0.0017) | −0.0024 (0.0015) | −0.0042 *** (0.0014) | −0.0020 (0.0014) | −0.0026 * (0.0014) |
Per Capita Grain Output (X10) | 0.0005 *** (0.0001) | 0.0006 *** (0.0001) | 0.0005 *** (0.0001) | 0.0006 *** (0.0001) | 0.0006 *** (0.0001) |
Population Urbanization Level (X11) | −0.0107 *** (0.0020) | −0.0106 *** (0.0017) | −0.0119 *** (0.0016) | −0.0102 *** (0.0016) | −0.0107 *** (0.0016) |
Population Density (X12) | −0.0006 *** (0.0001) | −0.0004 *** (0.0001) | −0.0003 *** (0.0001) | −0.0004 *** (0.0001) | −0.0005 *** (0.0001) |
Biological Abundance Index (X13) | −0.0055 (0.0037) | 0.0000 (0.0030) | −0.0012 (0.0030) | 0.0003 (0.0029) | −0.0016 (0.0028) |
NDVI (X14) | 0.0301 (0.0521) | 0.0622 (0.0566) | 0.0767 * (0.0421) | 0.0529 (0.0587) | −0.1156 (0.0745) |
Parameter ρ | — | 0.0974 (0.1134) | 0.4117 *** (0.0499) | — | 0.4242 *** (0.0588) |
Parameter λ | — | 0.4980 *** (0.1185) | — | 0.5793 *** (0.0573) | — |
LR Test: Individual Effect | 100.48 *** (0.0000) | 56.89 *** (0.0000) | 53.15 *** (0.0000) | 56.25 *** (0.0000) | 25.58 ** (0.0292) |
LR Test: Time Effect | 681.62 *** (0.0000) | −649.43 (1.0000) | −649.48 (1.0000) | −648.82 (1.0000) | −658.33 (1.0000) |
Hausman Test | 44.32 *** (0.0000) | — | 35.84 *** (0.0003) | 24.98 ** (0.0149) | 60.88 *** (0.0000) |
Individual Effect | Yes | Yes | Yes | Yes | Yes |
Time Effect | Yes | No | No | No | No |
Within R2 | 0.8838 | 0.8383 | 0.8448 | 0.8345 | 0.8823 |
Sample Size | 387 | 387 | 387 | 387 | 387 |
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
Cai, X.; Yang, R. Study on Spatiotemporal Evolution and Influencing Factors in Cultivated Land and Construction Land in Yunnan Province in the Past 20 Years Based on Remote Sensing Interpretation. Agriculture 2023, 13, 981. https://doi.org/10.3390/agriculture13050981
Cai X, Yang R. Study on Spatiotemporal Evolution and Influencing Factors in Cultivated Land and Construction Land in Yunnan Province in the Past 20 Years Based on Remote Sensing Interpretation. Agriculture. 2023; 13(5):981. https://doi.org/10.3390/agriculture13050981
Chicago/Turabian StyleCai, Xueli, and Renyi Yang. 2023. "Study on Spatiotemporal Evolution and Influencing Factors in Cultivated Land and Construction Land in Yunnan Province in the Past 20 Years Based on Remote Sensing Interpretation" Agriculture 13, no. 5: 981. https://doi.org/10.3390/agriculture13050981
APA StyleCai, X., & Yang, R. (2023). Study on Spatiotemporal Evolution and Influencing Factors in Cultivated Land and Construction Land in Yunnan Province in the Past 20 Years Based on Remote Sensing Interpretation. Agriculture, 13(5), 981. https://doi.org/10.3390/agriculture13050981