Spatio-Temporal Evolution of Vegetation Coverage and Eco-Environmental Quality and Their Coupling Relationship: A Case Study of Southwestern Shandong Province, China

: Propelled by rapid economic growth, the southwestern Shandong urban agglomeration (SSUA) in China has become a crucial industrial hub


Introduction
In recent decades, urban green spaces have increasingly been replaced by impermeable surfaces.This shift has exacerbated the conflict between the increasing demand for green spaces from citizens and the diminishing availability of urban green areas [1], and eco-environmental conditions have been deteriorating.Taking the example of southwestern Shandong Province, the southwestern Shandong urban agglomeration (SSUA) sits in the southwestern region of Shandong Province [2].In 2020, the Shandong Provincial Government enacted "the Guiding Opinions on Accelerating the Integrated Development of the southwest Shandong province Economic Circle", which stated that the collaborative management of the ecological environment and the ecological preservation and restoration of the SSUA should be promoted.As a key breakthrough area in the Shandong Comprehensive Pilot Zone for Replacing Old Growth Drivers with New Ones, the SSUA is undergoing a transition in traditional industrial services that must align with ecological conservation and coordinated development.Vegetation, collectively referring to plant communities covering the Earth's surface [3], is a vital metric for evaluating vegetation health and alterations in ecological systems [4], offering insights into the alterations of the eco-environment to a specific degree [5].Investigating the vegetation coverage (VC) and eco-environmental quality of urban agglomerations, along with their coupling coordination relationships (CCRs), contributes to a rational resolution of conflicts between economic growth and ecological integrity.The extraction of vegetation [6,7], assessment of ecological environment quality [8], spatio-temporal patterns of vegetation and eco-environmental quality [9][10][11], and coupling analysis of vegetation coverage or eco-environmental quality with other factors [12,13] are the four major research themes in the current research on VC and ecological quality (EQ).In terms of vegetation extraction, although there are vegetation extraction indices (e.g., Enhanced Mangrove Vegetation Index (EMVI) [14]; Vegetation Index (VI) [15]) and machine learning-based algorithms such as random forests applied to different scenarios [16], the Normalized Difference Vegetation Index (NDVI) has been tested to be more universally applicable in different types of data sources [17][18][19] and at different spatial scales [20].Secondly, EQ assessment methods are mainly divided into two categories: evaluation through specific ecological indicators [18] and indirect evaluation relying on approximate indicators such as vegetation cover or urbanization level [21].Specific ecological indicators are mainly divided into the Remote Sensing Ecological Index (RSEI) [22,23] and Improved RSEI [24,25].Chen et al. [26] introduced seasonal change factors when constructing an IRSEI, and although the IRSEI is more regionally representative, the RESI avoids the interference of subjective human factors, and the results are more objective.In the examination of spatial and temporal fluctuations in the VC and eco-environmental quality, geographic spatial analysis techniques such as Moran's index, coefficient of variation, linear regression, and the Hurst exponent have been widely applied [9,23,27].Although Moran's index and the standard deviational ellipse often analyze data at county or district levels, methods like trend analysis and the Hurst exponent can extend their analysis to the resolution of source data, allowing for a more detailed observation of regional changes.Finally, most of the current analyses on the coupling of fractional vegetation cover (FVC) and ecological conditions concentrate on investigating the coupling relationship between the two and their influencing factors [28], and few studies have concentrated on the coupling relationship between vegetation cover and ecological environment, so it is groundbreaking to investigate the interplay between vegetation and ecology in the context of industry.
In summary, there are several shortcomings in the research of regional VC and ecological environment quality: (1) The data processing based on abundant sources such as Landsat and MODIS faces challenges related to the large volume of data [10].(2) The study areas of the existing studies are mostly focused on watershed wetland ecological protection areas [8], with limited studies focusing on the long-term dynamics of industrial urban agglomeration.(3) Many studies focus solely on the factors influencing VC or changes in ecological conditions [21,28].There is insufficient research on the coordination relationship between VC and ecological conditions.Over the past two decades, the Chinese government has implemented a sequence of ecological restoration and protection policies [29], aimed at restoring ecological balance and improving environmental quality.However, the early models of urban development excessively consumed ecological resources [30], causing significant negative impacts on the environment.
In this research, based on the characteristics of the SSUA, our study focuses on two main aspects: First, we investigated whether there has been an improvement in the FVC and RSEI in the SSUA from 2000 to 2020.Second, we examined whether there are significant differences in the coupling relationship between the FVC and RSEI across different cities and between urban and rural areas within the SSUA.Using the Google Earth Engine (GEE) platform, we applied Theil-Sen median slope estimation, Mann-Kendall trend analysis, and Hurst exponent persistence analysis to examine the spatial and temporal dynamics and prospective developmental trajectories of VC in the study area.Additionally, the linkage coordination correlation amid VC and environmental quality in industrial urban agglomerations was investigated using the improved coupled coordination degree model (ICCDM).Our research offers a logical foundation for urban ecological environment preservation and planning, guiding the sustainable and robust advancement of the ecoenvironmental condition in the SSUA.

Study Area
The SSUA (114  57 ′ N) covers the southwest region of Shandong Province, China, encompassing the cities of Heze, Jining, and Zaozhuang, with a total area of about 27,803 km 2 (as shown in Figure 1).The region experiences a warm temperate, semi-humid monsoon continental climate, characterized by ample sunshine.The mean yearly temperature is 13.6 • C [31].The terrain primarily comprises plains and depressions, followed by lakes, mountains, and hills [32].This diverse landscape supports a rich array of vegetation, including broad-leaved forests, coniferous forests, shrublands, and wetland vegetation.In addition to its ecological diversity, the SSUA is endowed with significant underground resources, such as coal, iron, salt minerals, and gas, which have established it as a vital industrial hub within Shandong Province.As of 2023, the regional gross domestic product (GDP) totaled CNY 1213.77billion.Given its economic and ecological importance, examining the spatio-temporal fluctuations in the VC and ecological quality in the SSUA, and their coupling coordination relationship, is vital for fostering enduring and robust urban development.policies [29], aimed at restoring ecological balance and improving environmental quality.However, the early models of urban development excessively consumed ecological resources [30], causing significant negative impacts on the environment.In this research, based on the characteristics of the SSUA, our study focuses on two main aspects: First, we investigated whether there has been an improvement in the FVC and RSEI in the SSUA from 2000 to 2020.Second, we examined whether there are significant differences in the coupling relationship between the FVC and RSEI across different cities and between urban and rural areas within the SSUA.Using the Google Earth Engine (GEE) platform, we applied Theil-Sen median slope estimation, Mann-Kendall trend analysis, and Hurst exponent persistence analysis to examine the spatial and temporal dynamics and prospective developmental trajectories of VC in the study area.Additionally, the linkage coordination correlation amid VC and environmental quality in industrial urban agglomerations was investigated using the improved coupled coordination degree model (ICCDM).Our research offers a logical foundation for urban ecological environment preservation and planning, guiding the sustainable and robust advancement of the eco-environmental condition in the SSUA.

Study Area
The SSUA (114°45′ E-117°49′ E, 34°26′ N-35°57′ N) covers the southwest region of Shandong Province, China, encompassing the cities of Heze, Jining, and Zaozhuang, with a total area of about 27,803 km 2 (as shown in Figure 1).The region experiences a warm temperate, semi-humid monsoon continental climate, characterized by ample sunshine.The mean yearly temperature is 13.6 °C [31].The terrain primarily comprises plains and depressions, followed by lakes, mountains, and hills [32].This diverse landscape supports a rich array of vegetation, including broad-leaved forests, coniferous forests, shrublands, and wetland vegetation.In addition to its ecological diversity, the SSUA is endowed with significant underground resources, such as coal, iron, salt minerals, and gas, which have established it as a vital industrial hub within Shandong Province.As of 2023, the regional gross domestic product (GDP) totaled CNY 1213.77billion.Given its economic and ecological importance, examining the spatio-temporal fluctuations in the VC and ecological quality in the SSUA, and their coupling coordination relationship, is vital for fostering enduring and robust urban development.

Data Preparation
Yu et al. [10] utilized Landsat data to analyze changes in the vegetation cover in mining areas, while Xiong et al. [22] conducted a comprehensive assessment of the ecological quality of rubber plantations using Landsat data.Numerous studies have demonstrated that the Landsat dataset is broadly applicable for detailed observations over large spatial extents [33].GEE is a cloud computing platform that utilizes global-scale data for cloud analysis, which can conduct online analyses of long time-series geographic big data [32].It is particularly appropriate for the spatio-temporal analysis of vegetation cover and ecological environment quality over a 20-year time series, large spatial extent, and large volume of data in this study.Therefore, we used the GEE platform to select the MODIS dataset from 2000 to 2020, Landsat 5 Collection 2, Tier 1 data from 2000 to 2011, and Landsat 8 Collection 2, Tier 1 data from 2013 to 2020 (see Table 1) for the study area (less than 5% cloud cover).The cloud removal (mask function) [23], mosaicking, and cropping operations on the dataset were conducted on the GEE platform.Furthermore, as vegetation has a better growth pattern during the growing season, its environmental quality also tends to be better, so images between May and September annually were chosen to compile the remote sensing image collection VC for the respective year using the pixel dichotomy technique [8].Meanwhile, images between April and October annually were chosen to calculate the RSEI by calculating the Greenness index, Dryness index, and Wetness index to synthesize the ecological environment quality images for the target year.

Experimental Methodology
In this study, we use multi-temporal Landsat data and MODIS data from GEE to extract the VC and RSEI in southwestern Shandong Province.Trend analysis and persistence analysis were conducted using Theil-Sen median slope estimation [35], the Mann-Kendall significance test [35], and Hurst index [36].Additionally, the ICCDM [37] was employed to assess their coupling coordination.The workflow is illustrated in Figure 2. Firstly, Landsat and MODIS data were preprocessed on the GEE platform, involving cloud removal, mosaicking, and cropping operations, performed separately.Then, relevant characteristic parameters were computed, yielding the VC dataset and RSEI dataset for the research region for 20 years.Subsequently, utilizing Theil-Sen median slope estimation and the Mann-Kendall statistical test, we performed trend analysis, investigating the variations in the VC and environmental quality in the SSUA over the two decades.Additionally, employing the Hurst index method, we conducted persistence analysis to explore the future variations in the FVC and RSEI in the same region across the 20 years.Building on this, a coupling coordination analysis of the interrelation for FVC and the RSEI was conducted using the ICCDM, investigating the interdependent coupling relationship in the study area.

Pixel Dichotomy Technique
NDVI is the Normalized Vegetation Index, which assesses vegetation health by measuring the disparity between the near-infrared light reflected by vegetation and red light absorbed by it [38].The NDVI method, based on the pixel dichotomy model, is currently one of the most common and accurate techniques for extracting VC data [9].For instance, Zhang et al. [6] utilized this method to extract VC information for the Inner Mongolia segment of the Yellow River, and similarly, Jiang et al. [20] applied it to assess the VC in

Pixel Dichotomy Technique
NDVI is the Normalized Vegetation Index, which assesses vegetation health by measuring the disparity between the near-infrared light reflected by vegetation and red light absorbed by it [38].The NDVI method, based on the pixel dichotomy model, is currently one of the most common and accurate techniques for extracting VC data [9].For instance, Zhang et al. [6] utilized this method to extract VC information for the Inner Mongolia segment of the Yellow River, and similarly, Jiang et al. [20] applied it to assess the VC in the Nansi Lake region.These examples demonstrate the broad applicability of the NDVI method based on the pixel dichotomy model for vegetation extraction across different scales.In this method, the pixel's information includes fractions of pure vegetation and soil [1], and the mixed pixel's value is the weighted average of these two components' Vegetation Index values [39].The formula for the calculation is as follows: where NDV I soil represents the minimum value for pure soil pixels, theoretically close to 0; NDV I veg represents the maximum value for pure vegetation pixels, theoretically close to 1.

RSEI Method
The RSEI is a comprehensive ecological assessment index that effectively reflects the ecological quality status of a region [40].It is broadly utilized for ecological environment quality assessment in the field of remote sensing data [18,22,23].This comprehensive metric incorporates four key components: Greenness (NDVI), Wetness (LST), Dryness (NDBSI), and Heat (WET) [35].Given the variation in dimensions and value ranges among these indicators, normalization is necessary [41].Here, the RSEI is used to assess the ecological environmental quality in southwestern Shandong Province.The RSEI ranges between 0 and 1, where higher values indicate better ecological and environmental quality in the region.The Heat index is derived from the MODIS dataset in the GEE platform, and the other formulas are shown in Table 2.

Formulas Bibliographic Reference
Remote Sensing Ecological Index (RSEI) Greenness Index Wetness Index [43] [44] Dryness Index Here, NDV I is the Vegetation Index, WET is the Wetness index, NDBSI is the building and bare soil index; , and ρ SW IR2 are reflectance data in the blue, green, red, near-infrared (NIR), short-wave infrared 1 (SWIR1), and short-wave infrared 2 (SWIR2) bands for the TM and OLI, respectively.

Trend Analysis
The Theil-Sen Median is a robust, non-parametric statistical technique commonly used for trend analysis [47].This method is particularly effective in describing subtle trends in the improvement or degradation of parameters at the pixel level, which is essential for capturing fine-scale regional changes [9].For example, Zhang et al. [48] and Liu et al. [9] have conducted dynamic analyses of regional vegetation based on this method.In this research, the Theil-Sen median slope estimation and the Mann-Kendall significance test are utilized to assess the changing trends of the VC and eco-environmental condition in the SSUA over the past two decades.The equations for these analyses are as follows [35]: Forests 2024, 15, 1200 7 of 21 where x i and x j represent the average values for the i-th and j-th years; median denotes the median function.If T slope > 0, it signifies a statistically significant increasing trend in the dataset over the specified time period, whereas T slope < 0 indicates a decreasing trend.A larger absolute value of T slope indicates a more pronounced trend in data variability.The Mann-Kendall statistical test is employed alongside Theil-Sen median slope estimation to investigate the trend characteristics of raster data over extended time series [35].The equations for these analyses are as follows: , S < 0 (10) where n represents the length of the data series; S stands for the Mann-Kendall test statistic; x k and x j denote the corresponding data across the temporal data; sgn(x j − x k ) is the sign function.
The trend classification in this manuscript is presented in Table 3.

Persistence Analysis
In this research, the rescaled range (R/S) technique is utilized to compute the Hurst exponent [9], which is utilized to evaluate the persistence of VC and eco-environmental quality in Shandong Province [36].The Hurst exponent can characterize the persistence of lengthy time-series raster data and is extensively employed for forecasting future data [27,49].The equations for these analyses are as follows: For the given time series {Y(t)}, t = 1, 2, . . . . . ., n, the sequence of average values is explicated as follows: Accumulated dispersion: Range: Forests 2024, 15, 1200 8 of 21 Standard deviation: where H is the Hurst index; c is a constant; τ denotes the length of the time series.Where 0.5 < H < 1, the data exhibit persistence, meaning future trends align with the past; a higher H indicates stronger persistence.Conversely, when 0 < H < 0.5, the data show antipersistence, indicating that future trends oppose the past, and a lower H value indicates stronger anti-persistence.If H = 0.5 or is a null, the data behave as random fluctuations, with no connection to the past.

Improved Coupled Coordination Degree Model
The coupling relationship can effectively reveal the extent of benign interaction between two systems [50].However, traditional coupling coordination models exhibit a certain degree of subjectivity, with equal weights assigned to the contribution coefficients of each subsystem, lacking scientific rigor [50].To objectively evaluate the coupling relationship for the VC and eco-environmental condition, it is essential to consider the specific circumstances of the region.The relatively underdeveloped subsystem should be given a higher weight [51].The ICCDM can effectively narrow the gap between the two systems and provide a comprehensive evaluation of their respective contribution coefficients [37].Therefore, in this research, we employ an ICCDM to facilitate a more objective and accurate assessment of the coupling coordination relationship between the FVC and RSEI in southwestern Shandong Province.The equations for the ICCDM analyses are as follows [37]: where C represents the level of coupling for the FVC and RSEI; U and A, respectively, indicate the integrated scores of the FVC and RSEI; α ′ and β ′ , respectively, represent the improved contribution coefficients of the FVC and RSEI; T ′ and D ′ , respectively, represent the improved comprehensive levels of the FVC and RSEI and the improved coupling coordination.
Referring to relevant studies [37,52] and utilizing the specific data from the ICCM in the SSUA, the coupling coordination relationship between the VC and eco-environmental condition can be classified into three categories and eight subtypes (Table 4).As observed from Figure 3, during the period 2000 to 2020, the FVC in the western sector of the study region exhibited elevated levels compared to the eastern sector.Regions characterized by high vegetation cover were predominantly concentrated in the cities Heze and Jining, while areas with relatively higher and moderate vegetation cover were mainly concentrated in Zaozhuang City.The share of the primary sector in Heze City and Jining City was much greater than that in Zaozhuang City.Typically, the greater share of the agricultural industry, the higher the VC tends to be.This suggests that the economic focus on agriculture is a significant factor driving the higher vegetation cover in these regions.
According to Figure 4, it can be observed that over the past two decades, most regions of the SSUA have experienced improvements in the VC, especially in Zaozhuang, where the extent of high VC has notably increased, rising from 14% to 30%.On the other hand, Heze and Jining exhibited even higher levels of VC compared to Zaozhuang, both surpassing 60%.These changes indicate that the adoption of ecological policies, especially with As observed from Figure 3, during the period 2000 to 2020, the FVC in the western sector of the study region exhibited elevated levels compared to the eastern sector.Regions characterized by high vegetation cover were predominantly concentrated in the cities Heze and Jining, while areas with relatively higher and moderate vegetation cover were mainly concentrated in Zaozhuang City.The share of the primary sector in Heze City and Jining City was much greater than that in Zaozhuang City.Typically, the greater share of the agricultural industry, the higher the VC tends to be.This suggests that the economic focus on agriculture is a significant factor driving the higher vegetation cover in these regions.Simultaneously, during the 2000 to 2020 period, as economic development and infrastructure gradually improved, economic interactions among regions became increasingly intimate.The regions with low VC in Jining and Heze have gradually expanded outward towards urban economic centers and transportation routes.But VC in the SSUA exhibited a trend of gradual transition from lower to higher levels over the years.This trend reflects the impact of urbanization and improved connectivity on the spatial distribution of VC.

Trend Patterns of Vegetation Coverage
The Theil-Sen median slope estimation and Mann-Kendall significance test were employed to conduct trend analysis on VC in the southwest of Shandong Province.Figure 5 illustrates these changes, showing that approximately 34% of the region demonstrated a non-significant increase in VC over the two decades.This indicates a general trend of vegetation recovery in these areas.Conversely, areas expanding outward from urban centers showed a non-significant decreasing trend in the VC, comprising about 27% of the region.This highlights the impact of urban expansion on local vegetation.Examining specific cities, Heze and Zaozhuang largely showed notable vegetation improvements, except for some degradation near their city centers, where economic development is concentrated.In contrast, Jining exhibited substantial vegetation degradation, indicating a pronounced trend of vegetation deterioration.Urban development in Jining appears to be significantly affecting vegetation health.According to Figure 4, it can be observed that over the past two decades, most regions of the SSUA have experienced improvements in the VC, especially in Zaozhuang, where the extent of high VC has notably increased, rising from 14% to 30%.On the other hand, Heze and Jining exhibited even higher levels of VC compared to Zaozhuang, both surpassing 60%.These changes indicate that the adoption of ecological policies, especially with the gradual enforcement of initiatives like of the policy of converting cultivated land back to grove and ecological compensation, has greatly increased the VC.
Simultaneously, during the 2000 to 2020 period, as economic development and infrastructure gradually improved, economic interactions among regions became increasingly intimate.The regions with low VC in Jining and Heze have gradually expanded outward towards urban economic centers and transportation routes.But VC in the SSUA exhibited a trend of gradual transition from lower to higher levels over the years.This trend reflects the impact of urbanization and improved connectivity on the spatial distribution of VC.

Trend Patterns of Vegetation Coverage
The Theil-Sen median slope estimation and Mann-Kendall significance test were employed to conduct trend analysis on VC in the southwest of Shandong Province.Figure 5 illustrates these changes, showing that approximately 34% of the region demonstrated a non-significant increase in VC over the two decades.This indicates a general trend of vegetation recovery in these areas.Conversely, areas expanding outward from urban centers showed a non-significant decreasing trend in the VC, comprising about 27% of the region.This highlights the impact of urban expansion on local vegetation.Examining specific cities, Heze and Zaozhuang largely showed notable vegetation improvements, except for some degradation near their city centers, where economic development is concentrated.In contrast, Jining exhibited substantial vegetation degradation, indicating a pronounced trend of vegetation deterioration.Urban development in Jining appears to be significantly affecting vegetation health.

Persistence Analysis of Vegetation Coverage
Persistence analysis provides insights into future vegetation trends in southwestern Shandong Province, aiding regional coordination and sustainable development.The Hurst index, ranging from 0.06 to 0.99, was classified into different categories (0.05-0.3, strong anti-persistence; 0.3-0.49,weak anti-persistence; 0.49-0.51,random; 0.51-0.7,weak persistence; 0.7-0.99,strong persistence).Figure 6 shows VC persistence patterns and future development from 2000 to 2020 and statistical insights into these persistence categories and future VC.Understanding these persistence patterns helps predict future vegetation changes and supports informed regional planning.

Persistence Analysis of Vegetation Coverage
Persistence analysis provides insights into future vegetation trends in southwestern Shandong Province, aiding regional coordination and sustainable development.The Hurst index, ranging from 0.06 to 0.99, was classified into different categories (0.05-0.3, strong anti-persistence; 0.3-0.49,weak anti-persistence; 0.49-0.51,random; 0.51-0.7,weak persistence; 0.7-0.99,strong persistence).Figure 6 shows VC persistence patterns and future development from 2000 to 2020 and statistical insights into these persistence categories and future VC.Understanding these persistence patterns helps predict future vegetation changes and supports informed regional planning.
The analysis revealed that about 62% of southwestern Shandong exhibits a strong persistence in VC, totaling approximately 16,510 km 2 .Areas with weak persistence covered around 7597 km 2 , representing 29% of the region.Strong persistence is mainly found in areas radiating outward from city centers, indicating resilient vegetation.This suggests that regions with strong persistence are better equipped to maintain or improve vegetation over time.
In Heze, Jining, and Zaozhuang, the persistence distribution is relatively uniform, with strong persistence concentrated around urban centers.Overall, future predictions indicate that 53.1% of the area (about 14,060 km 2 ) will see vegetation improvements, while 43.5% (approximately 11,504 km 2 ) might experience degradation, particularly in expanding economic centers.The challenge lies in balancing urban expansion with ecological preservation to ensure sustainable development.
Controlling soil desertification and restoring vegetation are key future trends.Southwestern Shandong's urban agglomeration must align with these trends, promoting industrial structural reform to harmonize economic and ecological development.Successful integration of these priorities will determine the region's sustainable future.
Hurst index, ranging from 0.06 to 0.99, was classified into different categories (0.05-0.3, strong anti-persistence; 0.3-0.49,weak anti-persistence; 0.49-0.51,random; 0.51-0.7,weak persistence; 0.7-0.99,strong persistence).Figure 6 shows VC persistence patterns and future development from 2000 to 2020 and statistical insights into these persistence categories and future VC.Understanding these persistence patterns helps predict future vegetation changes and supports informed regional planning.

Characteristics of Spatio-Temporal Distribution of Environmental Condition
The RSEI provides a clear view of the fluctuations in eco-environmental quality in the SSUA.We categorized the RSEI into five degrees: poor (0-0.2),fair (0.2-0.4), moderate (0.4-0.6), good (0.6-0.8), and excellent (0.8-1).These classifications helped visualize the eco-environmental quality in southwestern Shandong Province (Figure 7) and the changes over time in specific cities like Heze, Jining, and Zaozhuang (Figure 8).This classification system is crucial for understanding and monitoring ecological health.
Over the past two decades, the western part of the SSUA has consistently shown a better eco-environmental quality than the eastern part.This trend is closely linked to the industrial structure in the region.Zaozhuang, with its high industrial activity, exhibited a lower eco-environmental quality due to significant air pollution and reduced ecological land.In contrast, Heze and Jining generally maintained a higher eco-environmental quality.Industrial activity directly impacts regional eco-environmental quality, with more pollutionintensive industries leading to poorer outcomes.
Throughout the past two decades, most regions in southwestern Shandong have seen improvements in their eco-environmental quality.Heze and Jining have shifted from a moderate to good eco-environmental quality, while Zaozhuang has moved from fair to moderate.The prevalence of good and moderate eco-environmental qualities around urban economic centers emphasizes the uneven distribution of environmental benefits.The overall trend suggests a positive shift towards better eco-environmental quality, especially in less industrialized areas.The RSEI provides a clear view of the fluctuations in eco-environmental quality in the SSUA.We categorized the RSEI into five degrees: poor (0-0.2),fair (0.2-0.4), moderate (0.4-0.6), good (0.6-0.8), and excellent (0.8-1).These classifications helped visualize the eco-environmental quality in southwestern Shandong Province (Figure 7) and the changes over time in specific cities like Heze, Jining, and Zaozhuang (Figure 8).This classification system is crucial for understanding and monitoring ecological health.Over the past two decades, the western part of the SSUA has consistently shown a better eco-environmental quality than the eastern part.This trend is closely linked to the industrial structure in the region.Zaozhuang, with its high industrial activity, exhibited a lower eco-environmental quality due to significant air pollution and reduced ecological land.In contrast, Heze and Jining generally maintained a higher eco-environmental quality.Industrial activity directly impacts regional eco-environmental quality, with more pollution-intensive industries leading to poorer outcomes.With ongoing ecological projects, the ecological condition in the SSUA has generally enhanced, transitioning from lower to higher levels.Areas with good and moderate qualities are expanding, while regions with poor quality are shrinking.However, urban expansion and industrial development continue to pose challenges.Balancing ecological restoration with urban and industrial growth is essential for sustainable development in the region.

Trend Characteristics of Eco-Environmental Condition
In this study, Theil-Sen median slope estimation and the Mann-Kendall significance test are used to conduct the trend analysis and statistics on the RSEI of the SSUA to obtain the changing trends in environmental conditions during the twenty years of study.As shown in Figure 9, most regions exhibited an improving trend in eco-environmental quality, covering approximately 11,460.563km 2 or 41.30% of the total area.Moving outward from the urban centers, there was a decreasing trend in the VC, though most declines were nonsignificant, covering about 5389.625 km 2 or 19.42% of the area.In Heze, eco-environmental quality degradation was more pronounced but still largely non-significant, aligning with urban development trends.This indicates that urban expansion has had a more noticeable influence on surrounding areas.In this study, Theil-Sen median slope estimation and the Mann-Kendall significance test are used to conduct the trend analysis and statistics on the RSEI of the SSUA to obtain the changing trends in environmental conditions during the twenty years of study.As shown in Figure 9, most regions exhibited an improving trend in eco-environmental quality, covering approximately 11,460.563km 2 or 41.30% of the total area.Moving outward from the urban centers, there was a decreasing trend in the VC, though most declines were non-significant, covering about 5389.625 km 2 or 19.42% of the area.In Heze, eco-environmental quality degradation was more pronounced but still largely non-significant, aligning with urban development trends.This indicates that urban expansion has had a more noticeable influence on surrounding areas.Heze experienced the most severe degradation in eco-environmental quality among the cities studied, which is largely due to its higher proportion of agricultural land and less industrialization compared to Jining and Zaozhuang.Jining's degradation ranked second, while Zaozhuang mostly saw improvements.Heze's predominantly agricultural base has made it more vulnerable to eco-environmental decline amidst urbanization.
Overall, while the majority of southwestern Shandong Province shows an increasing trend in eco-environmental quality, urban zones have experienced a decline.This trend is closely tied to the contrasting impacts of urbanization and industrialization.Balancing urban development with ecological preservation remains a key challenge for sustainable growth in the region.

Persistence Analysis of Eco-Environmental Condition
With the aim of forecasting the future trajectory of ecological conditions in southwestern Shandong Province, we performed a persistence examination of the RSEI using the Hurst index, which varied from 0.06 to 0.99.This index was categorized into five levels: strong anti-persistence (0.06-0.3), weak anti-persistence (0.3-0.49), random (0.49-0.51), weak persistence (0.51-0.7), and strong persistence (0.7-0.99).The results, shown in Figure 10, illustrate the persistence of ecological conditions from 2000 to 2020 and the projected future trends and statistical examination of the temporal persistence and projected future alterations in eco-environmental quality.Heze experienced the most severe degradation in eco-environmental quality among the cities studied, which is largely due to its higher proportion of agricultural land and less industrialization compared to Jining and Zaozhuang.Jining's degradation ranked second, while Zaozhuang mostly saw improvements.Heze's predominantly agricultural base has made it more vulnerable to eco-environmental decline amidst urbanization.
Overall, while the majority of southwestern Shandong Province shows an increasing trend in eco-environmental quality, urban zones have experienced a decline.This trend is closely tied to the contrasting impacts of urbanization and industrialization.Balancing urban development with ecological preservation remains a key challenge for sustainable growth in the region.

Persistence Analysis of Eco-Environmental Condition
With the aim of forecasting the future trajectory of ecological conditions in southwestern Shandong Province, we performed a persistence examination of the RSEI using the Hurst index, which varied from 0.06 to 0.99.This index was categorized into five levels: strong anti-persistence (0.06-0.3), weak anti-persistence (0.3-0.49), random (0.49-0.51), weak persistence (0.51-0.7), and strong persistence (0.7-0.99).The results, shown in Figure 10, illustrate the persistence of ecological conditions from 2000 to 2020 and the projected future trends and statistical examination of the temporal persistence and projected future alterations in eco-environmental quality.Over the past two decades, most areas in the southwest of Shandong Province exhibited strong persistence in eco-environmental quality, covering approximately 114,862 km 2 or 55% of the entire region.This suggests that more than half of the region has maintained a stable and strong eco-environmental quality.
Weak persistence was observed in a significant portion of the region, covering about 8718 km 2 or 32% of the area.Areas with anti-persistence were smaller and mainly situated in the northern sector of Heze City and the western zone of Jining City.Zaozhuang showed minimal anti-persistence.This pattern suggests varying degrees of eco-environmental stability across different regions.
Looking ahead, most areas in the southwest of Shandong Province are expected to improve in eco-environmental quality, covering around 17,070.5 km 2 or 84.8% of the entire area.Improvements will primarily occur in the eastern sector of Heze and the northern part of Jining.This optimistic forecast highlights the effectiveness of ongoing environmental restoration efforts.
However, some areas are predicted to experience eco-environmental quality degradation, affecting approximately 1841.8 km 2 or 9.1% of the total area.These areas are primarily around economic centers and transportation hubs, with the most severe degradation near Jining City.Targeted measures are needed to address environmental challenges in these vulnerable areas.

Assessment of the Coupling Coordination between Vegetation Coverage and Eco-Environmental Condition
Based on the ICCDM, the coordination analysis of VC and environmental conditions in southwestern Shandong Province was conducted, and this effectively explores the interactive relationship between VC and environmental quality in the study region.The results of this analysis in southwestern Shandong are depicted in Figure 11.The results indicate that over the past twenty years, the overall coupling coordination between the VC Over the past two decades, most areas in the southwest of Shandong Province exhibited strong persistence in eco-environmental quality, covering approximately 114,862 km 2 or 55% of the entire region.This suggests that more than half of the region has maintained a stable and strong eco-environmental quality.
Weak persistence was observed in a significant portion of the region, covering about 8718 km 2 or 32% of the area.Areas with anti-persistence were smaller and mainly situated in the northern sector of Heze City and the western zone of Jining City.Zaozhuang showed minimal anti-persistence.This pattern suggests varying degrees of eco-environmental stability across different regions.
Looking ahead, most areas in the southwest of Shandong Province are expected to improve in eco-environmental quality, covering around 17,070.5 km 2 or 84.8% of the entire area.Improvements will primarily occur in the eastern sector of Heze and the northern part of Jining.This optimistic forecast highlights the effectiveness of ongoing environmental restoration efforts.
However, some areas are predicted to experience eco-environmental quality degradation, affecting approximately 1841.8 km 2 or 9.1% of the total area.These areas are primarily around economic centers and transportation hubs, with the most severe degradation near Jining City.Targeted measures are needed to address environmental challenges in these vulnerable areas.

Assessment of the Coupling Coordination between Vegetation Coverage and Eco-Environmental Condition
Based on the ICCDM, the coordination analysis of VC and environmental conditions in southwestern Shandong Province was conducted, and this effectively explores the interactive relationship between VC and environmental quality in the study region.The results of this analysis in southwestern Shandong are depicted in Figure 11.The results indicate that over the past twenty years, the overall coupling coordination between the VC and RESI in southwestern Shandong has remained stable in an orderly and coordinated development state, with mean values transformed from 0.72 (2000) to 0.74 (2020).Specifically, in terms of time trends (Figure 11a,b), during the 2000 to 2020 period, various regions, such as Zones I and II, exhibited a state of joint expansion in their coupling coordination status, with their development direction aligning closely with urban planning direction.With the expansion of areas exhibiting coupling coordination imbalance in the economic core region of southwestern Shandong, the coupling coordination status in the Zaozhuang area gradually improved.The dominant status shifted from Zones I, II, III, and IV to Zones V, VI, and VII, indicating that with the implementation of the environmental regulations and the ecological civilization construction concept, further improvements had been made in the modification of the regional landscape and industrial composition, promoting a noticeable improvement trend in both the VC and eco-environmental quality.This progress indicated a further development in the beneficial interaction between the two.
Geographically, the coupling coordination relationship between the FVC and RSEI in the SSUA exhibited a network-like distribution pattern characterized by low coupling coordination degrees among various economic development centers (Figure 11a,b).Regions with imbalanced development in coupling coordination were mainly concentrated around the economic centers of each county.On one hand, with the progress of urbanization, urbanization and industrialization activities in the economic core areas, accompanied by pollution and land degradation, may have negative impacts on both VC and environmental quality to some degree.On the other hand, unreasonable land-use practices such as excessive cultivation and soil erosion may destroy habitats for vegetation and organisms, exacerbating the expansion of areas with imbalanced coupling coordination between economic centers and VC and ecological environment quality.Specifically, in terms of time trends (Figure 11a,b), during the 2000 to 2020 period, various regions, such as Zones I and II, exhibited a state of joint expansion in their coupling coordination status, with their development direction aligning closely with urban planning direction.With the expansion of areas exhibiting coupling coordination imbalance in the economic core region of southwestern Shandong, the coupling coordination status in the Zaozhuang area gradually improved.The dominant status shifted from Zones I, II, III, and IV to Zones V, VI, and VII, indicating that with the implementation of the environmental regulations and the ecological civilization construction concept, further improvements had been made in the modification of the regional landscape and industrial composition, promoting a noticeable improvement trend in both the VC and eco-environmental quality.This progress indicated a further development in the beneficial interaction between the two.
Geographically, the coupling coordination relationship between the FVC and RSEI in the SSUA exhibited a network-like distribution pattern characterized by low coupling coordination degrees among various economic development centers (Figure 11a,b).Regions with imbalanced development in coupling coordination were mainly concentrated around the economic centers of each county.On one hand, with the progress of urbanization, urbanization and industrialization activities in the economic core areas, accompanied by pollution and land degradation, may have negative impacts on both VC and environmental quality to some degree.On the other hand, unreasonable land-use practices such as excessive cultivation and soil erosion may destroy habitats for vegetation and organisms, exacerbating the expansion of areas with imbalanced coupling coordination between economic centers and VC and ecological environment quality.
When analyzed at the county level (Figure 11c,d), from 2000 to 2020, the dynamic coupling interaction between the FVC and RSEI in southwestern Shandong Province was mainly in an orderly and coordinated development state, primarily distributed in Heze City and Jining City, while most areas of Zaozhuang City exhibited a slightly coordinated development state.There was no significant change in the coupling relationship between the FVC and RSEI in each county of Heze City.The coupling relationship between the FVC and RSEI in Liangshan County and Wenshang County of Jining City had shifted from an orderly coordinated development to a dominant complementary and coordinated development, while Zoucheng City had transitioned from a scarcely coordinated development to a marginally coordinated development.However, there were also instances of degradation in the coupling coordination relationship between the VC and eco-environmental condition in some parts of Jining City, where Yutai County had shifted from a dominant complementary coordinated development to an orderly coordinated development.In Zaozhuang City, the coupling coordination relationship between the VC and eco-environmental quality in Shanting County had changed from barely coordinated to slightly coordinated, and in Tengzhou City, Yicheng County, and Taierzhuang County, the coupling coordination relationship had shifted from slightly coordinated development to orderly coordinated development.Zaozhuang City is an industrial city mainly focusing on chemical, machinery, and building materials industries, and it is also one of the important energy bases in East China.Although Heze City and Jining City had relatively diverse industrial structures, their agricultural industries still account for higher proportions compared to Zaozhuang City, leading to a lower level of economic development.Agricultural activities result in relatively lower damage to the ecological environment, and cities with higher proportions of agriculture often have a higher VC and ecological quality compared to cities with higher levels of industrialization.Therefore, the level of coordination between the FVC and RSEI in Zaozhuang City is slightly lower than that in the other two cities.With the improvement of regional vegetation conditions and environmental quality, there is also further improvement in the coupling coordination relationship among various counties.

Discussion
In terms of research methodology, we employed Theil-Sen trend analysis [35] and Hurst index analysis [36] to explore the spatio-temporal variations and future trends in vegetation and environmental conditions in the SSUA.Compared to spatial analysis methods such as elliptical variance and spatial autocorrelation analysis, which typically focus on city, county, or small-scale regions [23,27], these methods offer greater precision, allowing for a more detailed observation of subtle changes within the region [36].However, the limited resolution of Landsat imagery ignores, to some extent, the fragmented regional vegetation cover and ecological quality [30].Future research could utilize multisource high-resolution data to conduct more refined monitoring at a smaller scale.Moreover, most studies have been confined to the singular coupling relationship between VC, ecological environment quality, and their driving factors [28].This research enhanced our understanding by employing an improved coupled coordination model to explore the interaction mechanisms between the FVC and RSEI, thereby supplementing the understanding of vegetation ecological relationships.Nevertheless, the current findings on the coupled effects in the industrial city cluster of southwestern Shandong lack sufficient representativeness.Future studies could compare multiple typical industrial city clusters to explore the coupling mechanisms and internal dynamics of vegetation and ecology in greater depth.
In terms of the research findings, the spatio-temporal evolution of FVC and the RSEI across different regions of the SSUA in the past 20 years exhibits a general trend of improvement in VC and ecological quality.This observation aligns with the results reported by He [53] and Dong [54].The improvements can be largely attributed to policies such as the change of cultivated land to forest or grassland and the afforestation of barren hills [29].However, it is notable that areas around urban economic centers tend to have lower VC and ecological quality, with some areas even experiencing significant declines.This phenomenon mirrors findings in similar studies, where high ecological quality is typically found in urban green spaces and remnants of natural landscapes, while lower values are observed in areas undergoing extensive urban development [55].The increase in construction land has resulted in a corresponding decrease in ecological land use [56].Intense urbanization processes hinder vegetation and ecological development and tend to accelerate these issues with the expansion of urban areas [57].Urban construction also leads to landscape fragmentation, disrupting existing ecological connectivity.Overall, there is a spatially joint pattern of high vegetation cover and good ecological quality in the southwestern Shandong urban agglomeration (SSUA).Regions with high vegetation cover generally exhibit better ecological quality.This positive correlation between vegetation and ecological quality is supported by studies by Aizizi [56] and Yuan [58].However, significant spatial heterogeneity exists between FVC and the RSEI, with Jining and Heze showing markedly higher ecological quality and vegetation cover compared to Zaozhuang.The primary reason for this disparity is the geographical positioning of Jining and Heze at the base of the Taihang Mountains.The orographic lift effect of the mountains causes moist air masses to rise and cool, leading to increased precipitation on the windward slopes.Studies indicate that vegetation cover responds sensitively and timely to precipitation variations [59], and there is a positive correlation between rainfall and ecological quality [58].Additionally, Heze and Jining's predominant agricultural land use further enhances the regions' responsiveness to climatic changes, amplifying this effect [59].Regarding the coupling and coordination relationship between vegetation cover and ecological quality, most regions exhibit orderly and coordinated development.However, economic centers often show signs of discoordination.In megacities and supercities, higher levels of coordination are observed [60].This suggests a bipolar impact of economic development on regional coordination.In economically underdeveloped areas, economic activities may adversely affect regional coordination.This pattern indicates potential conflicts between economic activities and ecological systems in many industrially transitioning cities, necessitating further intervention by relevant authorities to promote harmonious development.
In the context of rapid urbanization, this research provides a fresh perspective on ecological research in the SSUA.It offers robust data support for the orderly and coordinated development and planning of heavy industrial agglomeration and urban areas, promoting the integration of ecological protection with economic growth.Future strategies should fully account for regional differences and be adaptable to the evolving environmental conditions.

Conclusions
This research sought to investigate the changes over time and across different areas of vegetation cover and ecological quality in the SSUA in China, as well as the coupling coordination relationship between these two factors.The principal results are summarized as follows: (1) Over the past two decades, the vegetation coverage and eco-environmental quality in the study area has a certain degree of improvement.These aspects have demonstrated a consistent shift from lower to higher levels.This discovery backs the hypothesis that urban and regional development policies have positively influenced ecological quality in the area.(2) Following the trend analysis results, most areas in the SSUA have exhibited an improving trend in both VC and ecological conditions over the past twenty years.This confirmation suggests that while overall trends are positive, urban expansion can negatively impact specific regions.(3) There will be a tendency towards improved VC and ecological conditions in the future in the majority of regions within the SSUA.However, there are still some regions where degradation is anticipated.Future changes in vegetation cover and ecological quality require researchers to continuously monitor and develop adaptive management strategies.(4) The coupling relationship between vegetation cover and ecological conditions in the southwestern Shandong region predominantly manifests as an orderly and coordinated state of development.However, there is a growing trend of discoordination near economic centers, indicating that urbanization and economic activities increasingly challenge the ecological balance.Areas with high economic activity are more susceptible to ecological imbalances.
This study still has several limitations that may influence the comprehensiveness and precision of the findings.Primarily, the reliance on a single dataset with a uniform resolution of 250 m after resampling may not capture the fine-scale variations in vegetation cover and eco-environmental quality.Additionally, while summer data highlight peak vegetation growth, they do not fully capture annual ecological trends.In regions like Shandong Province, with significant seasonal variations, this approach may miss crucial seasonal dynamics that are essential for a holistic understanding of the ecosystem.
To address these limitations, future research should integrate higher-resolution and more diverse remote sensing datasets to enhance the granularity of spatial analysis.Additionally, expanding the temporal scope to include data from all seasons will provide a more comprehensive view of annual ecological dynamics.

Figure 1 .
Figure 1.Study area in southwest Shandong Province.(a) Geographic location, (b) location of study area in Shandong Province, and (c) surrounding provinces.

Figure 1 .
Figure 1.Study area in southwest Shandong Province.(a) Geographic location, (b) location of study area in Shandong Province, and (c) surrounding provinces.

3 . 1 .
Characteristics of Vegetation Coverage in SSUA 3.1.1.Patterns of Spatio-Temporal Fluctuations in Vegetation Coverage To visually depict the spatial and temporal fluctuations in vegetation across SSUA, this study conducted a classification assessment of the VC.The classification was based on five categories: 0-0.2 representing low FVC, 0.2-0.4representing relatively low FVC, 0.4-0.6 representing moderate FVC, 0.6-0.8representing relatively high FVC, and 0.8-1 representing high FVC. Figure 3 presents the classified distribution of the vegetation density in the SSUA. Figure 4 depicts the distribution of vegetation cover in the cities of Heze, Jining, and Zaozhuang within the southwest of Shandong Province.

Forests 2024 , 22 Figure 5 .
Figure 5. Spatial distribution and statistical analysis of FVC trends in southwestern Shandong Province from 2000 to 2020.

Figure 6 .
Figure 6.Persistence characteristics of FVC in southwestern Shandong Province from 2000 to 2020.(a) Spatial distribution of FVC Hurst index, (b) future trends of FVC, (c) statistics on the grading of FVC Hurst index, and (d) statistics on the future trend of FVC.

Figure 5 .
Figure 5. Spatial distribution and statistical analysis of FVC trends in southwestern Shandong Province from 2000 to 2020.

Figure 6 .
Figure 6.Persistence characteristics of FVC in southwestern Shandong Province from 2000 to 2020.(a) Spatial distribution of FVC Hurst index, (b) future trends of FVC, (c) statistics on the grading of FVC Hurst index, and (d) statistics on the future trend of FVC.

Figure 6 .
Figure 6.Persistence characteristics of FVC in southwestern Shandong Province from 2000 to 2020.(a) Spatial distribution of FVC Hurst index, (b) future trends of FVC, (c) statistics on the grading of FVC Hurst index, and (d) statistics on the future trend of FVC.

Figure 9 .
Figure 9. Spatial distribution and statistical analysis of RSEI trends in southwestern Shandong Province from 2000 to 2020.

Figure 9 .
Figure 9. Spatial distribution and statistical analysis of RSEI trends in southwestern Shandong Province from 2000 to 2020.

Figure 10 .
Figure 10.Persistence characteristics of RSEI in southwestern Shandong Province from 2000 to 2020.(a) Spatial distribution of RSEI Hurst index, (b) future trends of RSEI, (c) statistics of the grading of RSEI Hurst index, and (d) statistics of the future trend of RSEI.

Figure 10 .
Figure 10.Persistence characteristics of RSEI in southwestern Shandong Province from 2000 to 2020.(a) Spatial distribution of RSEI Hurst index, (b) future trends of RSEI, (c) statistics of the grading of RSEI Hurst index, and (d) statistics of the future trend of RSEI.

Forests 2024 ,
15, x FOR PEER REVIEW 17 of 22 and RESI in southwestern Shandong has remained stable in an orderly and coordinated development state, with mean values transformed from 0.72 (2000) to 0.74 (2020).

Figure 11 .
Figure 11.Coupling coordination relationship results of FVC and RSEI in southwestern Shandong Province from 2000 to 2020.(a,b) Spatial distribution of coupling and coordination relationship between FVC and RSEI in 2000 and 2020, respectively.(c,d) Spatial distribution of coupling and coordination relationship between FVC and RSEI at county level in 2000 and 2020, respectively.

Figure 11 .
Figure 11.Coupling coordination relationship results of FVC and RSEI in southwestern Shandong Province from 2000 to 2020.(a,b) Spatial distribution of coupling and coordination relationship between FVC and RSEI in 2000 and 2020, respectively.(c,d) Spatial distribution of coupling and coordination relationship between FVC and RSEI at county level in 2000 and 2020, respectively.
Author Contributions: Conceptualization, D.M.; methodology, D.M. and Q.W.; software, Q.H. and Y.Y.; data curation, Z.L.; writing-original draft preparation, D.M.; writing-review and editing, Q.W.; visualization, Q.W. and Q.H.; supervision, Z.L.; project administration, D.M.; funding acquisition, D.M.All authors have read and agreed to the published version of the manuscript.Funding: This research was supported by the National Natural Science Foundation of China (grant number [42171435]); the Natural Science Foundation of Shandong Province (grant number [ZR2020MD025]); and the Doctoral Fund Projects of Shandong Jianzhu University (grant number [X21079Z]).

Table 1 .
Remote sensing data information and sources.

Table 3 .
[35]ciples for classifying the degree of change in the FVC and RSEI[35].