Urban resilience, the combinational characteristic of nature and society, that reflects the dynamic accumulation process that is multi-level and multi-dimensional. Particularly, the rational spatial distribution structure of land mixture and compactness is an effective way to improve urban resilience because the evolution of morphology and density of the urban land blocks in the process of land spatial conversion reflect the performance characteristics of complexity, diversity, stability, compactness, and connectivity. Therefore, we evaluated the relationship between urban resilience and land use and land cover (LULC) change, to find the keys to resilient urban development for urban land and space planning. In this study, taking the Chinese hilly city of Mianyang as an example, the results show: (1) the complexity of homogeneous patch shape and heterogeneous patch combination leads to the decrease of urban morphology resilience. (2) the development trend of LULC spatial layout and structure ratio were more rational with the increased of land mixing degree. (3) the speed and intensity of urban expansion were basically coordinated with the development of urban resilience. The research provides the new ideas, approaches, and toolkits for solving the intractable problems of urban spatial planning based on coordinating conflicts between urban resilience and urban land evolution.

Cities worldwide are facing the dual pressures of growing population and land expansion, leading to the intensification of conflicts in urban productive-living-ecological spaces (PLES). Therefore, the question of “how to dynamically judge the different thresholds of different indicators of PLES” plays an indispensable role in the studies of the multi-scenario simulation of land space changes and needs to be tackled in an appropriate way, given that the process simulation of key elements that affect the evolution of urban systems is yet to achieve complete coupling with PLES utilization configuration schemes. In this paper, we developed a scenario simulation framework combining the dynamic coupling model of Bagging-Cellular Automata (Bagging-CA) to generate various environmental element configuration patterns for urban PLES development. The key merit of our analytical approach is that the weights of different key driving factors under different scenarios are obtained through the automatic parameterized adjustment process, and we enrich the study cases for the vast southwest region in China, which is beneficial for balanced development between eastern and western regions in the country. Finally, we simulate the PLES with the data of finer land use classification, combining a machine learning and multi-objective scenario. Automatic parameterization of environmental elements can help planners and stakeholders understand more comprehensively the complex land space changes caused by the uncertainty of space resources and environment changes, so as to formulate appropriate policies and effectively guide the implementation of land space planning. The multi-scenario simulation method developed in this study has offered new insights and high applicability to other regions for modeling PLES.