Analyze the Centennial sustainable development strategy of Japan’s “Kawasaki model” and the evolution context of coastline layout, Explore the China’s shoreline sustainable development direction of heavy chemical industry: (1）The government should improves the regional technological level to solve the problem of shoreline pollution through diversified strategies, and supports pollution-intensive enterprises to become the “root vein” of innovation.（2)Straighten out the dependence between industry development and technology path, do not blindly introduce popular and gorgeous industries, adopt customized strategies to build a high value-added “The Product Space” of heavy chemical and related industries, and promote the original root enterprises to promote the shoreline industrial base advanced.（3)Establish “ The third sector” organizations led by the government, vigorously develop the S&T service industry, gather the team of professional technology transfer managers: To break the bottleneck of knowledge flow inside and outside the region, obtain the technical elements urgently needed by enterprises; To build a whole process support system for industrial chain and innovation chain to enhance the dynamic stability of strategic alliance; To build a regional horizontal digital platform, accelerate the pace of cluster innovation, and deep integration and optimization of production resources, form an open innovation pattern, promote the shoreline industrial foundation reengineering. Drawing on the Japanese “Kawasaki Model”, we will promote the path creation and shoreline sustainable development of heavy chemical industry with relatively low energy consumption and low pollution, environmental protection equipment manufacturing, new energy, life science and other emerging industries related to heavy chemical industry in China.

In this paper, the SBM model considering undesirable output is used to calculate the green innovation efficiency of high-tech industries in 30 provinces and autonomous regions (except Tibet ) in mainland China. Social network analysis and a spatial econometric model discuss the spatial correlation structure and influencing factors. The results show that: （1)The density of green innovation-related networks in the high-tech industry is significantly improved, but the evolution type is mainly low-level network, and the overall network structure is still in the initial stage. （2)Under different thresholds, the centrality of the green innovation correlation network of high-tech industries is significantly improved, but the spatial disequilibrium characteristics still exist and are more prominent under the fourth-level network threshold.（3)Governmental science and technology finance and informatization have a significant role in promoting green innovation efficiency. Industrial accumulation will lead to negative spatial spillover in promoting green innovation efficiency. The impact of environmental regulation is not apparent, and the level of opening up hinders the improvement of green innovation efficiency. We should continue to stabilize the leading position of Beijing-Tianjin-Hebei and Yangtze River Delta in green innovation of high-tech industries, improve regional integration, increase government investment in science and technology finance, improve the level of informatization, and promote the positive spatial spillover of green innovation.

Realizing the transition from factor-driven to innovation-driven is the key to my country’s sustainable and healthy development.China is already in a critical period of economic and industrial restructuring.Speeding up green technological innovation and industrial upgrading is an objective requirement for my country to transform from a manufacturing country to a manufacturing powerhouse and to achieve a rise in the industrial chain.In the current industrial upgrading and green development, it is necessary and practical to discover the role of green technological innovation and the degree of marketization in the process of industrial upgrading and evolution from agglomeration of innovative elements.Based on this, this article uses panel data from my country’s 31 provinces, autonomous regions, and municipalities from 2000 to 2019 to construct an evaluation index system for innovation factors, green technology innovation capabilities, and marketization levels, and comprehensive evaluation based on the entropy TOPSIS method.Use the mediation effect and moderation effect model test to explore the mediation effect of green technology innovation and the moderation effect of the degree of marketization in the process of the impact of the agglomeration of innovative elements on industrial upgrading.Concluded as follow: (1) Agglomeration of innovative elements can positively promote industrial upgrading. (2) The agglomeration of innovative elements can promote industrial upgrading by improving green technological innovation, and the ability of green technological innovation has a complete intermediary effect.(3) In the influence of the concentration of innovative elements on green technological innovation, the degree of marketization has a positive regulatory effect.Therefore, the government should promote the construction of factor markets, give full play to the decisive role of the market in resource allocation, and promote the formation of green technological innovation capabilities.

Based on the panel data of 104 prefecture-level cities in the Yangtze River Economic Belt from 2010 to 2019, this paper used the principal component analysis to calculate the comprehensive index of informatization development and the new-type urbanization development. And then, a spatial econometric model was constructed to empirically examine the direct effect and spatial spillover effect of informatization development on the new-type urbanization. The results show that: (1) From the global regression results, the development of informatization has a significant positive impact on the new-type urbanization, while the spatial spillover effect is not obvious. Among them, information infrastructure has the greatest impact on new-type urbanization, followed by information service consumption, while the information industry has no significant impact on the development of new-type urbanization, but there is a negative spatial spillover effect. (2) From the results of subdomain regression, there are obvious spatial differences in the impact of informatization development on new-type urbanization. The upstream area has the greatest impact, followed by the middle and lower reaches. (3) Other influencing factors: Both physical capital investment and regional consumption level have a positive effect on the development of new urbanization, and there is no spatial spillover effect. The degree of opening to the outside world has no significant impact on the development of new urbanization, and there is no obvious spatial spillover effect.

In this paper, taking 11 provinces and cities in the Yangtze River Economic Belt as the research object, the panel data from 2009 to 2018 is selected, and the green innovation efficiency of high-tech industries in various provinces and cities is calculated by Using Matlab Software based on the Super-SBM model considering undesirable output, constructing a comprehensive index of Internet development through principal component analysis. Finally, the two-way fixed effect model is used as a benchmark model to explore the role of Internet development and environmental regulation on the efficiency of green innovation in high-tech industries. Studies show that: (1) Under the condition of only considering regional effects and not adding environmental regulatory variables, the development of the Internet plays a positive role in the green innovation efficiency of high-tech industries.  (2) When considering regional effects and time effects, environmental regulations inhibit the improvement of green innovation efficiency of high-tech industries, while the development of the Internet also plays a role in inhibition, and the synergy between Internet development and environmental regulation can promote the improvement of green innovation efficiency of high-tech industries.  (3) Finally, it proposes to make full use of Internet technology, adopt differentiated environmental regulatory policies and integrate the development of the Internet and environmental regulations, with a view to improving the green innovation efficiency of high-tech industries in the Yangtze River Economic Belt.

The convergence of urban productivity in the Yangtze River Delta plays an important role in promoting the convergence of economic growth and is related to the overall development of high-quality integration in the Yangtze River Delta. This paper constructs a research framework of “measurement, decomposition, spatial differentiation and convergence identification” to reveal the phased characteristics, dynamic evolution law, regional differences, and convergence characteristics of urban productivity growth in the Yangtze River Delta. The results show that: from 2000 to 2018, the urban intensive economic growth driven by total factor productivity in the Yangtze River Delta is insufficient. From the decomposition results, the change of urban productivity in the Yangtze River Delta is more dependent on the “level effect” of technological efficiency improvement, and lacks the “growth effect” of technological progress. In the Yangtze River Delta, the productivity of Jiangsu, Zhejiang and Anhui provinces fluctuated in a mixed manner, and the over-variable density was the main source of the overall difference, but its contribution showed a steady downward trend. At present, the productivity of cities in Yangtze River Delta does not have the characteristic of convergence. This paper has a unique benchmark value and demonstration significance for promoting the high-quality integrated development of the Yangtze River Delta and realizing regional coordinated and sustainable development.

Hanjiang river ecological economy belt is an emerging economic belt across the east and west of the Yangtze River Basin，which has a unique and irreplaceable strategic position in the pattern of regional coordination development. Based on the panel data of 13 prefecture-level cities in Hanjiang river ecological economy belt from 2009 to 2019，this paper introduces the coupling coordination model and constructs a coupling coordination index system of traffic conditions-regional economy-ecological environment in Hanjiang river ecological economy belt to make a comprehensive empirical analysis on coupling coordination development of the three systems in Hanjiang river ecological economy belt. The results show that：（1）the overall coupling coordination of traffic condition-regional economy-ecological environment in Hanjiang river ecological economy belt shows a fluctuating upward trend，and the evaluation level of traffic condition and regional economy lags behind the evaluation level of ecological environment. （2）The comprehensive development level of the three systems between cities is uneven. For the evaluation level of traffic condition and regional economic，Wuhan is far ahead，followed by Luoyang，Xiangyang and Nanyang，and the rest of the cities are relatively backward. While the level of ecological environmental assessment in each city fluctuates greatly with a overall rising-falling-rising trend. (3）For each prefecture-level city，the coupling coordination degree of the three systems experienced stable or fluctuant ascent，generally developing in the direction of benign coordination and eased regional differences.

The industrial ecology is an important way and inevitable choice to realize a sustainable development of China’s economy and society. When taking a look at the concept of modern industrial ecology, we must closely follow the nature of industrial development to build an index system. It is able to scientifically measure the level of industrial ecologization with the method of factor analysis which deeply analyzes the spatial differentiation characteristics and spatial correlation of the industrial ecology with the ESDA method. Then mining the key influencing factors of industrial ecologization using geographic detector model. The results show that: (1) The industrial ecology of the Chongqing Metropolitan Area has significant advantages, the development power of Chongqing northeast urban agglomeration is insufficient, and  there is a short board effect in the industrial ecology of the “two clusters region”. (2) The city’s industrial ecological level presents a significant “core-edge” distribution characteristic, and the links of industrial energy conservation and reduction in emissions in the counties around the main urban area are relatively weak. (3) The city’s scale power and the spatial pattern of energy conservation and reduction in emissions are relatively stable. The radiation effect of the industrial ecology in the main urban area was found to be significant, and the gap of the industrial output efficiency has been gradually shortened. (4) There is a significant spatial correlation in the industrial ecology, and the spatial agglomeration effect has been increasing gradually. (5) The main influencing factors are the following. Industrial waste gas emissions per unit of the GDP, GDP per capita, total increase in secondary and tertiary industries, and the number of R&D personnel per one million people.

A relatively mature system of nature protected areas has been formed in Yangtze Basin. Revealing spatial distribution characteristics of nature protected areas and influencing factors can provide a scientific basis for the declaration, management, protection and sustainable development of nature protected areas in the watershed. Taking five types of national nature protected areas with a large number of Yangtze Basin as the research object, the spatial distribution differences of their quantity, area, geographic concentration and kernel density were analyzed, and the impact factors were discussed in depth. The results showed that: (1) There are a large number of national nature protected areas in Hunan, Sichuan, Jiangxi, Hubei, and Guizhou, with 152, 141, 125, 98, and 73 respectively. The total area of national nature protected areas of Qinghai, Sichuan, Hunan, Hubei, and Jiangxi has advantages, reaching 19 741 617, 3 893 913, 924 177, 739 134, and 643 504 km2 respectively. (2) National nature reserves (7.10), national scenic spots (7.16), and national water conservancy scenic areas (7.84) have the highest geographic concentration index in Hunan. National forest parks and national geological parks have the highest geographic concentration index in Jiangxi (7.24) and Guizhou (4.01) respectively. (3) The bandwidth value of the kernel density analysis for national nature reserves, national scenic spots, national forest parks, national water conservancy scenic spots, and national geological parks determined by the Ripley’s K method in Yangtze Basin are 20, 20, 20, 20, 25 km respectively. The spatial distribution of national nature reserves presents a pattern of one horizontal and two vertical strips. The spatial distribution of national scenic spots presents a pattern of one cluster in the lower reaches, two groups in the middle reaches, and a single hot core in the upper reaches. The spatial distribution of national forest parks presents a pattern of one cluster in the lower reaches, three clusters in the middle reaches, and a single core in the upper reaches. The spatial distribution of national water conservancy scenic spots presents a pattern of a single core in the lower reaches, two clusters in the middle reaches, and two cores in the upper reaches. The spatial distribution of national geological parks presents a pattern of one cluster in the lower reaches and one horizontal and two vertical strips in the middle and upper reaches. (4) The spatial distribution of various types of national nature protected areas is comprehensively affected by many different factors such as topography, vegetation, water resources, biodiversity, population, land, capital, transportation, and policy.

As one of the first national park trials in China, Mastering the ecosystem health status and historical developing process of Shennongjia forest，that is conducive to formulating and implementing Shennongjia forest development policies. This paper collected social economy, social economy and policy measures data from 2013 to 2019 in Shennongjia Forest region, and constructed a comprehensive evaluation index system of ecosystem health and a coordinated evaluation index system of socio-economic and natural ecological coupling in Shennongjia forest region. The ecosystem health status, development index and coupling coordination degree of Shennongjia forest were evaluated respectively. The results indicated that the ecosystem health status, development coordination degree and system development indexes of shennongjia forest region showed a positive trend from 2013 to 2019, and the social economy and natural ecology dimensions kept a high coupling degree.The evaluation results are in line with the expectation of national policies and measures, and can be used as a reference to evaluate the achievements of shennongjia National Park pilot work.

To explore the pattern of spatialvariations in phytoplankton community structure and the characteristics of floracomposition in the main stem of the Yangtze River basin，we set 74 sampling sectionsat five researchregions(i.e.，headwater region，Jinsha Riverregion，Three Gorge reservoir，middle and lower reaches of the Yangtze River) in the main stem to conduct the field survey on phytoplankton in July-September of 2019. A total of 112 species，belonging to 6 phyla in 53 genera，were identified during the survey. Bacillariophyta was the dominant taxa except the lower reach area. Diatoms were predominant at five research areas，e.g.，Synedra acus，Diatoma vulgare，Melosira granulata var. Angustissima. The mean density and biomass of phytoplankton throughout the main stem were 1.29×106 cells/L and 1.34 mg/L，respectively. Statisticallysignificant differences were foundin thespecies composition，Diatom Quotient and the diversity indices among five research areas. Specifically，the phytoplankton species in the headwater region were extremely dissimilar to those at other regions. The Diatom Quotient atthe lower reaches were significantly higher than those in other regions (P<0.05). The diversity indices in the Jinsha River area and Three Gorges Reservoir region were significantly higher than those atother regions (P<0.05). The indices of Shannon-Wiener diversity (H’)，Margalef richness (D) and Pielou evenness (J) were employed to evaluate the ecological health of water bodies in the Yangtze River basin. The results show that the overall water quality in the main stem was amoderate pollution level. Pearson correlation analysis shows that WT (water temperature)，TN (total nitrogen)and NH3-N (ammonium)were the primaryfactors affecting the densityand biomass of phytoplankton，while DO (dissolved oxygen)was the most important factor affecting the phytoplankton diversity. Our results will provide systematic basic data for the implementation of the “10-year Yangtze fishing ban” ，and supply withscientific basis for the management and protection of aquatic ecology in the main stem of the Yangtze River.

In recent years, the water resources security and ecological impact of the middle route of the South-to-North Water Diversion Project has been widely concerned, and the sedimentation problem gradually appears in the slow flow and backwater areas of the middle route, and the suspended solids is an important source in the main canal. To solve this problem, this study carried out research on the spatial distribution and dynamics of suspended solids in the main channel of the middle route of the South-to-North Water Diversion Project in November 2020 and July 2021. The concentration and composition of total suspended solids showed obvious spatial and seasonal differences. The concentration of volatile suspended solids was higher than that of non-volatile suspended solids in the north area of the water diversion channel. The generalized additive model was used to analyze the relation between phytoplankton cell density and suspended solids composition, the result showed the concentration of volatile suspended solids is significantly related to phytoplankton cell density, meaning phytoplankton account for a large proportion of volatile suspended solids. According to phytoplankton community analysis, Aphanocapsa sp and Melosira granulata were the dominant species in November 2020; Synedraacus and Aphanocapsa pulchra were the dominant species in July 2021. The cell density in July 2021 is higher than November 2020, the Bacillariophyta, Cyanophyta, Chlorophyta accounts for the largest proportion of biovolume and cell density in the community, they all positively correlated with dissolved oxygen and water temperature. Thus, it is necessary to pay attention to the growth and reproduction of Bacillariophyta species in the main canal to prevent them from increasing the concentration of volatile suspended solids and affecting the water receiving in the northern area of the main canal

In order to study the function of the fish habitat conservation area of the tributary Bachu River after the construction of Batang hydropower station in the upper reaches of Jinsha River, the monitoring of fish resources in the main stream below Batang dam and the tributary Bachu River was carried out in 2021, and the fish species composition and community structure were analyzed by using diversity index and multivariate statistical methods. A total of 11 fish species were collected in the main stream and tributary, including 4 endemic species and 3 alien species in the upper reaches of the Yangtze River. The dominant species in the main stream are Schizothorax, and the dominant species in the tributary are two small-sized fishes, Schizopygopsis malacanthus and Triplophysa stenura. Fish community of nine sampling locations were partitioned into two groups at the 43.53% Bray Curtis similarity level. We attributed the difference between these two groups to the different dominant species in the main river and the tributary Bachu river according to SIPMER analysis. The biodiversity index of the two sampling locations near the confluence of Bachu river were the highest, which may be related to the complexity of the habitat at the confluence. It is suggested to strengthen the protection of the habitat at the confluence. At the same time, restoring the connectivity of the upper reach of the confluence of Bachu river may prolong the migratory distance of Schizothorax, which is conducive to expand its population scale.

Surface ozone (O₃) pollution in the Yangtze River Economic Belt (YREB) region is becoming increasingly severe, which is not beneficial to the integrated regional O₃ prevention and prevention. Based on the trend analysis, hotspot analysis, spatial clustering and spatio-temporal geographically weighted regression model (GTWR), we investigated the complex nonlinear relationships between the spatiotemporal distribution, evolution patterns and drivers of O₃ concentrations in the YREB region from 2013 to 2020.The results show that: (1) The annual average O₃ concentration in the YREB is over 100 μg/m3, with a significant increase trend and an average increase rate of 2.88 μg/m³/year (p<0.05); Spatially, the average O₃ concentration shows a distribution pattern of high in the northeast and poor in the west. (2) The hot spot regions of O₃ concentration variations rate were distributed in Hubei, Anhui, Jiangsu, Shanghai, Zhejiang, Jiangxi and Hunan regions; the cold spot regions were distributed in eastern Sichuan, eastern Chongqing, western Hubei and Hunan and eastern Guizhou. (3) The proportion of secondary industry in GDP was the dominant socioeconomic driver of O₃ concentration in the YREB, with a regression coefficient of 0.328. Precipitation was the primary meteorological driver of O₃ concentration, followed by temperature, relative humidity, wind speed, mean visibility and pressure; O₃ concentration was significantly positively associated with temperature and negatively correlated with precipitation, relative humidity, wind speed, mean visibility and barometric pressure. This investigation could provide a benchmark for the prevention and control of air pollution in the YREB.

The average wind speed in the southeastern Wuhan is dominated by soft breeze, which decreases slowly in central urban area with heat island being concentrated, while keeping stable fluctuation in suburbs with no heat island in most areas. The study constructs a marco ventilation corridor system in the southeastern Wuhan, and then performs numerical sensitivity experiments to examine impacts of three-typed corridors on the meteorological fields around the surrounding areas by WRF. The simulation period is set on 18th July 2018 with the prevailing southeast wind. Two ventilation corridors are first set up, which are parallel and perpendicular to the prevailing wind, respectively. The results indicate that the urban ventilation corrisor can decrease temperature, increase humidity and wind speed. Compared with no corridors, the corridor parallel to the prevailing wind can make the near surface temperature dropping by 1.4℃, the relative humidity increasing by 8% and the wind speed increasing by 1.4 m·s^-1. The effect of cooling and humidifying is obvious within a vetical height of 100 m, and the wind speed increase by 1 m·s^-1 at 120 m. The effect of the corridor perpendicular to the prevailing wind on temperature, humidity and wind speed is weaker than that of the parallel one, and the effect of the composite corridors is basically the same as that of the single one. The influence of ventilation corridors on temperature and humidity extends downstream along the wind, but little on wind speed. The temperature decreases by 0.98℃ and the relative humidity increases by 5.93% for 1 km downwind of the corridor parallel to the prevailing wind, while 0.46℃ and 3.12% respectively for 2 km. The study shows that the ventilation corridors are designed to connect the suburban green space and water body with the central heat island area by the prevailing east wind, which introduce the cold natural wind and export the urban hot air to form local convection for alleviating urban heat island.

In this paper, we analyzed the contents and δ34S values of total sulfur, carbon-bonded sulfur, ester-bonded sulfate, SO₄^2-, total reduced inorganic sulfur and stream SO₄^2-, and SRB quantity at different depths in calcareous soil and yellow soil profiles in Chenqi and Houzhaihe small karst catchments of Puding County in Guizhou Province, with the aim of elucidating the response of soils in karst catchments to atmospheric acid deposition. The results showed that under the same acid deposition level, there were significant differences in the sulfur content and δ34S values of different sulfur forms and SRB numbers in the calcareous soil and yellow soil profiles. Also more than 70% of the SO₄^2- in the stream water of sampled slope came from soils at different depths in calcareous soil and yellow soil profiles. These results indicate the different response of calcareous soil and yellow soil to atmospheric acid deposition. The formation of organic sulfur and dissimilatory SO₄^2- reduction are sulfur retention processes that exist in both calcareous soil and yellow soil profiles. SO₄^2- adsorption exists in the yellow soil profile and accumulates in the bottom soil layers; therefore, SO₄^2- adsorption is the main cause of sulfur retention in yellow soil. The SO42- exported from the calcareous soil and yellow soil profiles is mainly from SO₄^2- input from atmospheric acid deposition and not from organic sulfur mineralization of the soil itself. However, organic sulfur is the dominant sulfur form in calcareous soil and yellow soil, and organic sulfur mineralization in these two types of profiles may release a large amount of SO₄^2- into the surface water after the sulfur deposition flux is significantly reduced.

The heavy metals accumulated in the sediments is a potential threat to the urban river ecology. In August 2019, sediments were collected in Qingshui River, located in the main urban area of Chongqing. The content of heavy metals in the samples and the environmental indicators of adjacent surface water were monitored. The ecological risks and sources of heavy metals were evaluated by various statistical methods. The results showed that the average value of As, Cd, Cr, Cu, Hg, Ni, Pb and Zn in the sediments of the main stream of Qingshui River were 7.55±8.25, 0.238±0.056, 86.43±5.14, 78.21±18.22, 2.49±2.00, 33.86±5.29, 42.99±9.21 and 198.71±43.55 mg·kg^-1, respectively. According to the result of Geoaccumulation Index (Igeo) evaluation system, Hg reached the Severe Pollution Level, while Cu, Zn, Cd and Pb were at Slightly Polluted Level. The heavy metal content of sediments at most sites exceeded the current agricultural land soil pollution risk control standards, which indicated the dredged silt should be carefully used in agricultural land. The mean value of sediment Potential Ecological Risk Index (RI) was 1 984.54±1 515.92 and Hg contributed 90.78% for it, showing a great ecological risk of the sediments in Qingshui River. The exogenous sources of the heavy metals could be apportioned to the headwater and higher pH of the upstream surface water., which resulting in a gradual decrease trend of heavy metals (except As and Ni) from upstream to downstream. Due to the high gradient of mountainous rivers, the exogenous inputs of heavy metals at the headwater are able to affect the content of heavy metals in the sediments covering the entire upper and middle reaches. Principal component analysis and cluster analysis showed that Cr, Cu, Hg, Pb, and Zn in the sediments mainly originated from industrial point sources at the headwater, and there were three types of subdivisions derived from different types of industries. As and Ni mainly originated from sewage treatment station/urban non-point sources. Cd accumulation in the sediments were both significantly affected by the mentioned two sources.

This paper studies the effects of subalpine coniferous forest degradation caused by dwarf bamboo expansion on soil K and Mg contents. We used community diversity indexes, soil pH, soil organic matter and the contents of N, K and Mg in litter and soil to compare among primary subalpine coniferous forest (PF), mixed broadleaf-conifer forest (SF) and dwarf bamboo bush (BB). In addition, stepwise regression was used to explore the main factors affecting soil K and Mg contents. The results showed that: (1) plant numbers in tree and shrub layers, species richness and diversity indexes decreased significantly with forest degradation; (2) the contents of TNlitter, TKlitter and TMglitter displayed significant differences among the three degradation levels, which was highest in PF, and lowest in BB; (3) with the forest degraded to bamboo bush stage, SOM, TNsoil, ANsoil, TKsoil, AKsoil and TMgsoil significantly decreased; (4) stepwise regression analysis showed that TNsoil contributed 57.0% and 68.0% to TKsoil and AKsoil, respectively, SOM contributed 35.9% to AKsoil, and ANsoil contributed 47.1% to TMgsoil, indicating that soil N and SOM content were the main factors affecting soil K and Mg contents, followed by community diversity. Our results proposed that the degradation of subalpine coniferous forest caused by dwarf bamboo expansion not only reduced soil C and N content, but also led to soil K and Mg loss.

In order to clarify the trade-off / synergy relationship between the main ecosystem services of typical karst mountain area, taking Wenshan city of Yunnan Province as the research object, based on theoretical framework of Millennium Ecosystem Assessment (MA), four ecosystem services including food supply, net primary productivity (NPP), soil conservation and water yield were selected to evaluate the temporal and spatial changes of land use and ecosystem services with the application of CASA model, InVEST model, RUSLE model and related analysis methods. The results show that from 2000 to 2017, the area of forest and unused land in Wenshan City has decreased, and other land use types have increased relatively, of which construction land has increased the most; The annual average value of ecosystem service in Wenshan City has shown an overall increasing trend in 17 years, with inter-annual changes increasing and decreasing; There is a synergistic relationship between water yield and food supply, both of them are in a trade-off relationship with NPP and soil conservation, and there is a synergistic relationship between NPP and soil conservation; In terms of spatial distribution, the ecosystem services are mainly synergistic, the spatial patterns of trade-off / synergy relationship are different, which show that Wenshan downtown area has significant synergy or significant trade-off relationship, while other areas in the study area are scattered.