The relationship and obstacles between sustainable development (SD) and ecosystem services (ESs) are crucial factors for SD decision-making and ecological conservation strategies. The production-living-ecological (PLE) functional perspective provides a new research entry point to analyze the interrelationship between the SD of human society and ESs. In view of this, based on the Sustainable Development Goals (SDGs), this study established an SD localization evaluation framework from the perspective of the PLE functions and subsequently evaluated the SD levels of 16 cities in Yunnan Province from 2005 to 2020. Furthermore, changes in the ecosystem service value (ESV) were explored by combining the PLE land classification and dynamic equivalence methods. The relationship between SD and ESV was then analyzed using the modified coupling coordination degree (CCD) model. Finally, the obstacle degree model was introduced to identify the main factors affecting the CCD level of the two systems. The results revealed the SD level in Yunnan Province to be high in the central region and low in other areas. From 2005 to 2020, the province’s SD level increased from 0.146 to 0.341. Furthermore, the main contribution of ESV originated from ecological land (over 90%). The ESV was higher in the southwest of Yunnan Province compared to that in the northeast. From 2005 to 2020, the ESV decreased by 2.231 billion CNY and was dominated by the regulation services. The CCD level increased significantly from 2005 to 2020. The CCD subcategory changed from SD lag to ESV lag and ESV-SD balanced type. Moreover, the dominant obstacle factors did not change significantly during the study period. The ESV and SD systems were dominated by regulation services and sustainable production subsystems, respectively. The research results have important guiding significance for further optimization of regional development strategies and ecological protection measures.

In the Qinghai-Tibet Plateau, fractures and fault zones have developed in this area, where earthquakes and extreme rainfall frequently induce debris flow disasters, which seriously threaten the safety of the people and properties. In this study, Zhujiagou, Minxian County, Gansu Province in China, has a typical debris flow channel in a small watershed on the northeastern edge of the Qinghai-Tibet Plateau, which has been used as the case study for risk assessment. A debris flow-risk assessment method has been developed in this study. The fusion algorithm of entropy weight and coefficient of variation has been used to evaluate the susceptibility of debris flow in each tributary channel in the watershed. Further, numerical simulation of the debris flow events at the main channel and at the high-prone debris flow tributary channel has been carried out. The “7.18” Zhujiagou debris flow event has been used to verify that the accuracy of the numerical simulation is higher than 74.12%. This method has been applied to design the Zhujiagou debris flow-risk zoning under extreme rainfall conditions, aimed at forming a set of debris flow-risk evaluation system suitable for this type of small watershed. The results show that the peak single-width clear water flow at the mouth of the channel, the length of the channel, and the shallow surface rock formation are the main factors influencing the susceptibility of debris flow. The main tributaries of the watershed are Hagu channel, Zhuling channel, Songshu channel, and Langjia channel which are all high-prone debris flow channels. The Zhujiagou debris flow accumulation fan under the designed 1% rainfall frequency will block the Taohe, which will threaten the safety of the residents and properties near the channel, at the mouth of the channel and in the urban area on the opposite bank of the Taohe. Compared with Hagugou, Zhulinggou, and Songshugou, the low-risk area is Langjiagou. The research results of this study can be used for evaluation methods and as a basis for preventing debris flow disasters on the northeastern edge of the Qinghai-Tibet Plateau.

This paper takes the surrounding rock of deep tunnel as the research object and considers the action mechanism under the influence of seepage. Based on the Mohr-Coulomb criterion, the stress mechanism of surrounding rock of deep buried tunnel is analyzed by a convergence constraint method. Based on the elastic-plastic solution, the nonlinear elastic-plastic solution of the interaction between surrounding rock and lining structure considering the effect of seepage force is proposed, and the radius of surrounding rock plastic zone is obtained. The relationship between surrounding rock stress and displacement, radial deformation of lining, and support reaction force was observed. At the same time, considering the effects of seepage, strain softening, and intermediate principal stress, the surrounding rock is divided into a plastic residual zone, plastic softening zone, and elastic zone, and the stress distribution expressions of the plastic zone and each zone of surrounding rock of circular tunnel are derived. The results show that with the change of nonuniform permeability coefficient, the seepage shows anisotropy in different directions, and the closer to the horizontal or vertical direction, the more obvious the influence of nonuniform permeability coefficient on pore water pressure distribution. Seepage and material softening have different effects on the distribution of surrounding rock stress field and the size of plastic zone. Material softening is more unfavorable to the stability of surrounding rock than seepage. The intermediate principal stress coefficient has a significant impact on the tangential stress and plastic zone of surrounding rock. When the intermediate principal stress effect is not considered, the calculation results are relatively conservative and cannot give full play to the strength of surrounding rock effectively. The research conclusion can provide a theoretical reference for studying the stability of surrounding rock in tunnel excavation under water-bearing rock.