Spatiotemporal Evaluation of Regional Land Use Dynamics and Its Potential Ecosystem Impact under Carbon Neutral Pathways in the Guangdong–Hong Kong–Macao Greater Bay Area

The spatiotemporal distribution of ecosystem service values (ESVs) and ecological risk are critical indicators to represent the regional ecological protection level and potential of sustainable development, which largely depend on land-use patterns. Aiming to contribute to global climate mitigation, China has proposed dual-carbon goals that would remarkably influence the land-use/cover change (LUCC) distribution. Based on the Landsat land cover data of 2000, 2010 and 2020 and multisource satellite products, several driving factors are integrated into the patch-generating land use simulation (PLUS) model to simulate future LUCC patterns for the Guangdong–Hong Kong–Macao Greater Bay Area (GBA) under rapid urbanization, cropland protection and carbon neutral (CN) scenarios from 2020 to 2050. Spatial–temporal ecosystem service and ESVs are allocated using INVEST and the equivalent factor method and thus ecological risks are evaluated using the entropy method. Results indicate that forest growth is the largest under the CN scenario, especially in the northwestern and northeastern GBA, exceeding 25,800 km2 in 2050, which results in both the highest habitat quality and carbon storage. The largest ESVs, reaching higher than 5210 yuan/pixel, are found in the CN scenario, particularly expanding toward the suburban area, leading to the lowest ecological risks. From 2020 to 2050, habitat quality, carbon storage and ESVs improve, while ecological risks decline in the CN scenario. This research provides implications for economic and ecological balanced development and gives references to the carbon-neutral pathway for the GBA.