Spatial Identification and Evaluation of Land Use Multifunctions and Their Interrelationships Improve Territorial Space Zoning Management in Harbin, China

Quantitative assessment and trade-off/synergy analysis of land use multifunctions can effectively identify regional conflicts and dominant functions, providing decision support for promoting sustainable socio-economic and land use development. However, current research in this field still faces challenges due to coarse scale of studies and limited availability of accurate data. Taking Harbin City as a research case, this research employed an improved mutation level method, Pearson correlation analysis, and a multi-scale geographically weighted regression model to comprehensively investigate the profiling of land use multifunctions and their trade-off /synergy relationships. The comparative advantage theory was adopted to identify dominant functional zones using the NRCA index at a grid scale, in order to achieve a territorial spatial functional zoning delineation. The results showed that there were intricate trade-off/synergy relationships among production–living–ecology functions. Moreover, the types and intensity of trade-off/synergy evolved continuously with socio-economic development and regional resource endowment disparities. Due to its exceptional resource endowment, the agricultural dominated, urban dominated, and ecological dominated functional areas accounted for a significant proportion of 29%, 7%, and 26%, respectively. However, owing to the intricate trade-offs/synergies inherent in land use multifunctions, only a mere 2% (agricultural), 1% (urban), and 1% (ecological) of the area were identified as Optimization Guidance Zones. Conversely, Remediation Improvement Zones constituted the largest share at 63% of the total area, with agricultural, urban, and ecological Remediation Improvement Zones accounting for approximately 33%, 12%, and 18%, respectively. Based on the results of the type and intensity of trade-off/synergy among production–living–ecological functions, as well as the dominant zones and the integrated territorial spatial functional zoning delineation, this article provided targeted recommendations for the sustainable development of the region. These recommendations took into account both resource endowment and socio-economic development characteristics specific to the study area. The study aims to contribute to related research gaps, while providing valuable insights for other regional studies.