| Literature DB >> 36231918 |
Zihao Li1,2, Xihang Xie3, Xinyue Yan4, Tingting Bai5, Dong Xu6.
Abstract
The market entry of rural collective operating construction land (MERCOCL) is an important way for the Chinese government to promote the marketization of rural land. However, the impact of China's Rural Land Marketization on the ecological environment quality (EEQ) remains to be understood. Understanding these mechanisms is necessary for regional sustainable development and rational resource allocation. Therefore, a universal assessment model of China's regional EEQ was built based on the Landsat 5/8 and the national ecological index (EI) provided by the Ministry of Ecology and Environment at the national district and county scale. A total of 229 counties (32 pilot counties and other counties in the pilot cities) in China from 2011 to 2018 were taken as the research object. This paper empirically studied the evolution process, driving mechanism and spatial heterogeneity of EEQ from the perspective of MERCOCL. The study shows that China's EEQ presented a spatial distribution pattern of "high in the south, low in the north, high in the east and low in the west". When a county implemented the MERCOCL policy, its EEQ index increased by 0.342, with the improvement effect occurring in the second year after the MERCOCL implementation. Regarding the mechanism, MERCOCL mainly improved the EEQ by promoting industrial structure optimization and increasing urban population aggregation. From the perspective of spatial heterogeneity, the improvement effect of MERCOCL on EEQ was more significant in regions with lower economic development levels and latitudes (southern China). This study will facilitate an understanding of the impact of China's rural land marketization on the EEQ and provide scientific data support for government departments to formulate sustainable urban development policies that meet local conditions.Entities:
Keywords: ecological environment quality; market entry of rural collective operating construction land; rural land marketization
Mesh:
Year: 2022 PMID: 36231918 PMCID: PMC9566321 DOI: 10.3390/ijerph191912619
Source DB: PubMed Journal: Int J Environ Res Public Health ISSN: 1660-4601 Impact factor: 4.614
Detailed description of remote sensing data.
| Data | Format | Spatial Resolution | Temporal Resolution | Source |
|---|---|---|---|---|
| Landsat 5 | TIFF | 30 m | 16-day | USGS |
| Landsat 7 | TIFF | 30 m | 16-day | USGS |
| Landsat 8 | TIFF | 30 m | 16-day | USGS |
| MCD12Q1 | HDF | 500 m | Annual | NASA |
NASA: National Aeronautics and Space Administration (https://www.nasa.gov/, accessed on 2 June 2022). USGS: United States Geological Survey (https://www.usgs.gov/, accessed on 2 June 2022).
Figure 1Geographical distribution of pilot areas.
Descriptive statistics of variables.
| Variables | Indicator Meaning | Observation | Mean | S.D. | Minimum | Maximum |
|---|---|---|---|---|---|---|
| EEQ | Ecological environment quality (-) | 1832 | 54.324 | 12.196 | 24.896 | 89.521 |
| ey | Economic growth (10,000 $) | 1832 | 1.105 | 1.061 | 0 | 6.3 |
| s | Industrial structure (%) | 1832 | 44.932 | 15.426 | 4.67 | 92.85 |
| den | Population density (-) | 1832 | 5.740 | 0.923 | 2.792 | 8.657 |
| pop | Population agglomeration (%) | 1832 | 30.587 | 19.024 | 8.18 | 99.67 |
| str | Industrial structure optimization (%) | 1832 | 36.819 | 10.816 | 10.3 | 75.26 |
| gov | Government intervention (%) | 1832 | 24.132 | 16.699 | 2.26 | 105.58 |
| temp | Temperature (°C) | 1832 | 14.39 | 4.54 | 3.01 | 24.76 |
| rain | Rainfall (mm) | 1832 | 990.34 | 460.44 | 197.17 | 2516.57 |
Figure 2Accuracy verification of CHEQ and RSEI based on ecological index (EI) provided by the Ministry of Ecology and Environment of the People’s Republic of China (MEEPRC): (a) Accuracy verification of CHEQ and RSEI in northeast China; (b) Accuracy verification of CHEQ and RSEI in north China; (c) Accuracy verification of CHEQ and RSEI in east China; (d) Accuracy verification of CHEQ and RSEI in northwest China; (e) Accuracy verification of CHEQ and RSEI in southeast China; (f) Accuracy verification of CHEQ and RSEI in central south China.
Figure 3Spatial distribution maps of EEQ (CHEQ) in China in 2018. Spatial distribution map of EEQ in China in 2018 (a) at the 1000 m grid scale and (b) at the country scale.
Figure 4Spatial and temporal variation maps of EEQ in 229 counties: (a–h) Spatial distribution maps of EEQ in 229 counties from 2011 to 2018; (b) EEQ heat map of 32 pilot counties from 2011 to 2018 (i).
Results analysis of MERCOCL on the effect of EEQ.
| Variables | (1) | (2) | (3) | (4) | (5) | (6) | (7) |
|---|---|---|---|---|---|---|---|
| EEQ | EEQ | EEQ | EEQ | EEQ | EEQ | EEQ | |
| dt | 0.342 ** | 0.392 ** | 0.389 ** | 0.454 *** | 0.459 *** | 0.400 ** | 0.444 ** |
| (0.171) | (0.167) | (0.168) | (0.169) | (0.169) | (0.173) | (0.180) | |
| ey | 1.169 *** | 2.518 *** | 2.047 *** | 2.042 *** | 1.786 *** | 1.862 *** | |
| (0.245) | (0.577) | (0.603) | (0.610) | (0.603) | (0.615) | ||
| ey2 | −0.284 *** | −0.224 ** | −0.229 ** | −0.202 * | −0.210 * | ||
| (0.109) | (0.109) | (0.111) | (0.108) | (0.111) | |||
| s | 0.022 * | 0.022 * | 0.026 ** | 0.022 * | |||
| (0.013) | (0.013) | (0.013) | (0.013) | ||||
| den | −0.696 | −0.603 | −0.536 | ||||
| (0.851) | (0.900) | (0.885) | |||||
| rain | 0.0008 *** | 0.0008 *** | |||||
| (0.0001) | (0.0001) | ||||||
| temp | −0.114 | ||||||
| (0.085) | |||||||
| Intercept term | 54.297 *** | 53.001 *** | 52.177 *** | 51.575 *** | 55.591 *** | 54.295 *** | 55.655 *** |
| (0.014) | (0.272) | (0.420) | (0.582) | (5.030) | (5.342) | (5.434) | |
| Year fixed | YES | YES | YES | YES | YES | YES | YES |
| County fixed | YES | YES | YES | YES | YES | YES | YES |
| Observation | 1832 | 1832 | 1832 | 1832 | 1832 | 1832 | 1832 |
|
| 0.003 | 0.017 | 0.023 | 0.027 | 0.027 | 0.034 | 0.035 |
Note: ***, ** and * represent the significance levels of 1%, 5%, and 10%, respectively.
Applicability test of the PSM-DID method (common support hypothesis).
| Variables | Mean of the Treatment Group | Mean of the Control Group | Difference | ||
|---|---|---|---|---|---|
| ey | 1.325 | 1.178 | 0.147 | 1.22 | 0.222 |
| ey2 | 2.883 | 2.322 | 0.561 | 0.99 | 0.324 |
| s | 46.954 | 47.445 | −0.491 | −0.30 | 0.765 |
| den | 5.953 | 5.903 | 0.050 | 0.50 | 0.621 |
| temp | 13.836 | 13.903 | −0.067 | −0.12 | 0.902 |
Results analysis of the impact of MERCOCL on EEQ using PSM-DID.
| Variables | (1) | (2) | (3) | (4) | (5) | (6) | (7) |
|---|---|---|---|---|---|---|---|
| EEQ | EEQ | EEQ | EEQ | EEQ | EEQ | EEQ | |
| dt | 0.397 ** | 0.419 ** | 0.403 ** | 0.498 *** | 0.510 *** | 0.461 *** | 0.490 *** |
| (0.169) | (0.167) | (0.169) | (0.171) | (0.171) | (0.175) | (0.182) | |
| ey | 0.905 *** | 2.534 *** | 1.800 *** | 1.787 *** | 1.651 *** | 1.742 *** | |
| (0.227) | (0.541) | (0.575) | (0.578) | (0.590) | (0.595) | ||
| ey2 | −0.365 *** | −0.263 *** | −0.269 *** | −0.253 *** | −0.266 *** | ||
| (0.096) | (0.092) | (0.090) | (0.091) | (0.092) | |||
| s | 0.030 ** | 0.030 ** | 0.032 ** | 0.029 ** | |||
| (0.014) | (0.014) | (0.014) | (0.015) | ||||
| den | −1.324 | −1.353 | −1.248 | ||||
| (1.083) | (1.110) | (1.114) | |||||
| rain | 0.0005 *** | 0.0006 *** | |||||
| (0.0002) | (0.0002) | ||||||
| temp | −0.090 | ||||||
| (0.091) | |||||||
| Intercept term | 53.555 *** | 52.580 *** | 51.585 *** | 50.768 *** | 58.502 *** | 58.109 *** | 58.846 *** |
| (0.015) | (0.246) | (0.406) | (0.601) | (6.465) | (6.614) | (6.707) | |
| Year fixed | YES | YES | YES | YES | YES | YES | YES |
| County fixed | YES | YES | YES | YES | YES | YES | YES |
| Observation | 1646 | 1646 | 1646 | 1646 | 1646 | 1646 | 1646 |
|
| 0.004 | 0.012 | 0.021 | 0.027 | 0.029 | 0.032 | 0.033 |
Note: Control variables are gradually added to Columns (1)–(7) in the same order as Table 3. ***, ** represent the significance levels of 1% and 5%, respectively.
Driving mechanism test of MERCOCL on EEQ.
| Variables | Benchmark Regression | Structure Optimization | Government Intervention | Population Agglomeration | |||
|---|---|---|---|---|---|---|---|
| (1) | (2) | (3) | (4) | (5) | (6) | (7) | |
| EEQ | str | EEQ | gov | EEQ | pop | EEQ | |
| dt | 0.444 ** | 1.279 * | 0.439 ** | −0.829 | 0.443 ** | 4.794 *** | 0.307 * |
| (0.180) | (0.721) | (0.180) | (0.697) | (0.180) | (1.674) | (0.182) | |
| str | 0.003 ** | ||||||
| (0.0015) | |||||||
| gov | −0.001 | ||||||
| (0.008) | |||||||
| pop | 0.029 *** | ||||||
| (0.007) | |||||||
| Intercept term | 55.655 *** | 25.544 | 55.568 *** | 40.425 *** | 55.609 *** | −132.54 ** | 59.436 *** |
| (5.435) | (20.997) | (10.26) | (14.343) | (5.387) | (65.431) | (4.347) | |
| Control variable | YES | YES | YES | YES | YES | YES | YES |
| Year fixed | YES | YES | YES | YES | YES | YES | YES |
| County fixed | YES | YES | YES | YES | YES | YES | YES |
| Observations | 1832 | 1832 | 1832 | 1832 | 1832 | 1832 | 1832 |
|
| 0.035 | 0.093 | 0.036 | 0.095 | 0.035 | 0.163 | 0.067 |
Note: ***, ** and * represent the significance levels of 1%, 5%, and 10%, respectively.
Figure 5Parallel trend test.
Time sensitivity test.
| Variables | 2014–2016 | 2013–2017 | 2012–2018 |
|---|---|---|---|
| EEQ | EEQ | EEQ | |
| dt | 0.273 * | 0.298 ** | 0.493 *** |
| (0.145) | (0.150) | (0.182) | |
| Intercept term | 60.239 *** | 54.299 *** | 56.339 *** |
| (5.810) | (6.306) | (6.244) | |
| Control Variable | YES | YES | YES |
| Year fixed | YES | YES | YES |
| County fixed | YES | YES | YES |
| Observations | 687 | 1145 | 1603 |
|
| 0.043 | 0.078 | 0.039 |
Note: ***, ** and * represent the significance levels of 1%, 5%, and 10%, respectively.
Policy interference test.
| Variables | (1) | (2) | (3) | (4) | (5) |
|---|---|---|---|---|---|
| EEQ | EEQ | EEQ | EEQ | EEQ | |
| dt | 0.444 ** | 0.444 ** | 0.437 ** | 0.427 ** | 0.424 ** |
| (0.180) | (0.181) | (0.180) | (0.184) | (0.185) | |
| talk | −0.295 | −0.287 | |||
| (0.196) | (0.194) | ||||
| low | 0.113 | 0.041 | |||
| (0.133) | (0.135) | ||||
| cet | 0.555 *** | 0.548 *** | |||
| (0.145) | (0.147) | ||||
| Intercept term | 55.655 *** | 54.957 *** | 55.859 *** | 57.421 *** | 56.793 *** |
| (5.435) | (5.573) | (5.470) | (5.362) | (5.523) | |
| Control variable | YES | YES | YES | YES | YES |
| Year fixed | YES | YES | YES | YES | YES |
| County fixed | YES | YES | YES | YES | YES |
| Observations | 1832 | 1832 | 1832 | 1832 | 1832 |
|
| 0.035 | 0.038 | 0.036 | 0.041 | 0.043 |
Notice: Column (1) shows the model that only includes the MERCOCL; Column (2) shows the model that includes talk (Environmental Protection Interview); Column (3) adds low (low-carbon city pilots); Column (4) adds cet (carbon emission trading); Column (1) shows the model that includes all policies. ***, ** represent the significance levels of 1% and 5%, respectively.
Heterogeneity analysis of regional characteristics.
| Variables | Higher Economic Development Level | Lower Economic Development Level | Southern Region | Northern Region |
|---|---|---|---|---|
| EEQ | EEQ | EEQ | EEQ | |
|
| 0.337 ** | 0.390 | 0.385 ** | 0.436 |
| (0.168) | (0.321) | (0.191) | (0.304) | |
| Intercept term | 51.623 *** | 55.921 *** | 49.052 *** | 62.749 *** |
| (3.603) | (11.608) | (4.141) | (9.722) | |
| Control variable | YES | YES | YES | YES |
| Year fixed | YES | YES | YES | YES |
| County fixed | YES | YES | YES | YES |
| 7 | 616 | 1216 | 1016 | 816 |
|
| 0.070 | 0.071 | 0.077 | 0.031 |
Note: ***, ** represent the significance levels of 1% and 5%, respectively.